I'm trying to model atmospheric ballistic trajectories. I think I have a good approximation using excel (but not so much on the calculus). I modeled acceleration, thrust, weight, drag and have an approximation to drag through the sound barrier. I have assigned a drag coefficient that seems to work. My model works pretty good launch and through sound barrier. After SRB burnout I model a ballistic trajectory to descent back through the sound barrier as far as chutes out and my sim is over.
Problem : Ascent Velocity from immediately past sound barrier > 360m/s to engine burnout and Ap isnt matching my model and is out buy about 5%. My model is faster than flight test comming out of sound barrier
I've got TWR correct ie model matches KE readings. It took a bit but eventually I got a fudge for Thrust efficiency against atmospheric pressure with a decent stdev. I used a logistic curve to approximate the efficiency curve that KE reports in the VAB (yeah whatever but it works)
I've modeled Kerbin Atmosphere and drag F_Drag=1/2.P.v^2.K where K is a drag coefficient I tweek till matches "reality" (well flight recording), I did flight recording by screenshotting every few seconds and entering the data back to excel. (maybe there is a quicker way but its ok).
I know its not that I got the sound barrier approximation wrong from 300-360m/s because my model matches recording all the way to 360m/s then starts diverging.
The brunt is : is there something happening with drag as you go supersonic (>360m/s) that I need to consider?
Thanks ahead to any help, I would like to show the graphs but I'm not sure how to post pics, Im hoping there is something silly I forgot.
PS I'm new to the game its my first and only career , never did sandbox have a ton of mods, killed three Kerballs and have started a science program that so I dont loose more greenies. I play very slow but I'm really loving getting back into basic physics its really wonderful.
.