Уважаемые форумчане, скажите пожалуйста, возможно ли применение следующей модели для расчёта максимальной скорости пули при вылете её из гладкоствольного ружья?
Имеется гладкий ствол заданной длины

и цилиндрическая пуля массой

и радиусом

.
Положим, что в результате взрыва пороха, между ударником и пулей образуется газ, давление которого задано и равно, скажем,

. считая процесс расширения
пороховых газов адиабатным примем:

тогда в итоге, по второму закону Ньютона:

Здесь

- некий первоначальный "зазор" между ударником и пулей.
Далее, умножим обе части уравнения на

и проинтегрируем их:

Здесь

. Так как

И окончательно:

Подставляя найденные мною в Интернете значения некоторых постоянных получаются более менее реальные цифры.
Остаётся вопрос - можно ли вообще использовать столь грубый способ для решения подобных задач?
Всем заранее спасибо за комментарии и высказанные мнения.