I'm writing a simulation to determine the most efficient thrust profile for a rocket to achieve the highest altitude with the least amount of fuel. Given the following:
mass (as a function of time) diameter Cd air density
The results were surprising. But before I state the results I got (in order to have a "blind" confirmation), has anyone else attempted this kind of calculation before? And if so, what results did you get?
Thanks in advance! Dave
Didn't find your answer? Ask the community — no account required.
J
Jerry Irvine
Yes.
The results are in books called "Flight sheet guides" and "Malewicki Charts."
Jerry
formatting link
A
Alan Jones
Hint: this is a classic problem problem known as the Goddard problem.
What did you find surprising?
Alan
D
dave.harper
I'm not asking to find a ideal weight for a given engine... What I'm talking about is doing it in reverse. Say you're given a rocket with parameters Cd and D. Now you're given "x" amount of fuel (with an isp of 'you pick')... if you could control the thrust curve, what's the highest you can make it go?
For example, say if I use a K1100 with a rocket with a diameter of
5.5", a Cd of .75, and an initial mass of 5.8kg (which is about optimum for that engine). That'll get you an altitude of about 3872 ft. But using an optimized thrust curve for the same rocket with same fuel (and isp), I'm predicting an altitude of 5013 ft.
It's not as applicable to solids as it is to hybrids (or liquids). It's not a really big difference until you start looking at optimizing curves for spaceshots, where fuel is the majority of your weight and atmospheric density changes significantly.
Dave
Jerry Irv> >
D
David Weinshenker
For typical model/HPR conditions, optimum thrust/time for altitude seems to be a bit on the long duration side, compared to what you really want for a clean hard liftoff - but a "kicker" motor generates high velocity at low altitude, which coasts off rapidly... some motors seem to hit a "sweet spot" (in "typical" consumer rocketry airframes) - Cesaroni Pro38's, Econojet G38FJ, and Apogee 18mm D10 mini composites come to mind: they come off the pad clean and hard, and then keep going and going like the energizer bunny...
-dave w
"dave.harper" wrote:
D
dave.harper
That the thrust curve that maximizes altitude seeks to satisfy: V(t)=Terminal Velocity (C,A,rho(t))
i.e. the most efficient vertical velocity is the same as terminal velocity at a given altitude for a given rocket.
Dave
J
Jerry Irvine
The results are in books called "Flight sheet guides"
What is your test case?
B
Brett Buck
Assuming you don't care about practical considerations, and are only doing
1-dimensional analysis - slow, so you don't waste any more energy on drag than necessary.
In practice it doesn't quite work the same way, since you have to handle perturbations.
Brett
A
Alan Jones
Close enough. The optimal thrust for most MR and HPR is bang-sigular arc-coast. In practical terms it is max trust (ideally an impulse) to get you up to speed, followed by approximately T=2W until burn out, followed by coast to apogee.
Now then, assuming a case where your optimal trust time exceeds 15 seconds, how does the optimal thrust change when a 15 second burn time limit is imposed?
Alan
D
dave.harper
Yeah, T=2W gives you terminal velocity, only in the positive direction. As you start getting into the higher atmosphere, that ratio starts to increase, and in a vaccuum it'd be infinite.
Given the same amount of fuel as the thrust curve that exceeded 15 seconds? I imagine it would be boost, sustain at terminal velocity (i.e. T=2W), then at 14.9 seconds, an impulse of all remaining fuel.
Dave
D
dave.harper
My test case? You mean this?
optimum
Dave
J
Jerry Irvine
5.5 x 80" rocket CDr 0.55 8-20lbs 80 degrees 29 inhg at ground level. Neutral Thrust Motor CDA (IN^2)=13.066677 ORBIT98.BAS COPYRIGHT JERRY IRVINE WEIGHT ALTITUDE BURNOUT BURNOUT MACH ALTITUDE COAST (LBS) (FT) ALT (FT) VEL(FPS) MAX (MILES) TIME (S)
5.5 x 80" rocket CDr 0.55 8-20lbs 80 degrees 29 inhg at ground level. Progressive Thrust Motor CDA (IN^2)=13.066677 ORBIT98.BAS COPYRIGHT JERRY IRVINE WEIGHT ALTITUDE BURNOUT BURNOUT MACH ALTITUDE COAST (LBS) (FT) ALT (FT) VEL(FPS) MAX (MILES) TIME (S)
5.5 x 80" rocket CDr 0.55 8-20lbs 80 degrees 29 inhg at ground level. Stepped Thrust Motor (realistic) CDA (IN^2)=13.066677 ORBIT98.BAS COPYRIGHT JERRY IRVINE WEIGHT ALTITUDE BURNOUT BURNOUT MACH ALTITUDE COAST (LBS) (FT) ALT (FT) VEL(FPS) MAX (MILES) TIME (S)
First, are you accounting for reduction in mass due to propellant burn? Just checking to make sure we're using the same methods.
Anyway, using this one as an example, I get 7966 ft max altitude, 5816 ft at burnout, and a velocity of 608 ft/s at burnout, so we're relatively close (within about 5%).
Using an optimized thrust profile with the same amount of propellant, same mass, Cd, dia, isp, etc... I get a max alt of 9487 ft, burnout alt of 8565 ft, and a burnout velocity of 292 ft/s. That's almost a 20% increase in altitude using the same amount of fuel (and using mine as the reference). That's what I mean by optimizing the thrust.
Dave
A
Alan Jones
Well, yes T=-2W would not get you much positive altitude. "Terminal velocity" is a poor description in this case. I usually refer to the optimal cruise (or singular arc) speed to distinguish it from the initial optimal acceleration phase, and the final coasting phase.
That is true, but even large HPR flights seldom exceed 20K Ft. T=2W is quite practical for typical MR, HPR flights. There are two interesting consequences of thrusting to stay on the optimal singular arc. On the one hand, as mass burns off, optimal thrust tends to decrease, while on the other hand as density decreases with altitude, trust tends to increase to accelerate the rocket and maintain "optimal drag' (or even "terminal velocity"). For some things like MR altitude with an AT E6 sustainer, constant thrust can be just about right.
Yes, and assume constant ISP for all thrusting, just to keep it simple. The 15 seconds is cumulative for all thrust phases, but impulse thrusting is OK. Actually I'd like to suggest an example of a
62.5 gram, 200 sec ISP, propellant MR motor with less than 80N average thrust, but of course you also need a rocket mass and Cd such at that the optimal thrust is well over 15 seconds, making the 15 second limit active in the thrust time lime limited optimization problem.
OK. I had a more impractical solution in mind. I should probably just solve this problem, but it is more fun to pose it as an open problem, especially when you said that you had a simulation program set up to optimize thrust.
Alan
D
dave.harper
I'd have to add some coding to make it a simulation program to "optimize thrust" given specific limitations... the simulation I use right now only limits the max thrust of the engine.
Dave
J
Jerry Irvine
Yes using the thrust curve as a proxy.
But you do not say what that thrust profile is.
I do.
Also my run is a practical buildable motor.
Your run if ever disclosed could probably be approximated by a buttkicker booster motor cluster (0.4s) and a central endburner with a very low probabibility of vertical flight.
A
Alan Jones
I find your maximum thrust limitation odd, and interesting.
A related problem is to assume a high thrust motor that is not throttleable, but that can be turned on and off at will. E.g. an 80N motor that would accelerate the model to speed and then be turned on and off in some duty cycle to maintain the desired speed.
Alan
BTW, which Dave Harper are you? Is you dad's name Roy? Where do you live now?
D
dave.harper
COAST
TIME
203.796914
propellant,
Well, I'll post the full profile when I get home tonight, if that'll make it more believable.
Considering I haven't stated what my thrust profile is (as you mention above), this is a big assumption, isn't it?
Actually, it is pretty similar (in profile), not nessesarily in a that configuration... Much like your "Stepped Thrust Motor" profile, only with a lower, longer sustained flight (which is what the optimization simulation calculates). And it hits 200mph in under a second, so I don't know why you say it'd probably fail to go vertical...?
Dave
J
Jerry Irvine
We'll look forward to it.
Nope. A unlateral statement.
The low thrust combined with gravity turn and greatly accelerated by wind perturbation.
Tech Jerry
D
dave.harper
program
If I didn't limit thrust, the rocket would experience over 1000g's (10,000+ lbs of thrust) during the first loop iteration. The optimizing routine would try to get the rocket to it's most efficient vertical speed in 1 step. Secondly, if it did accelerate to it's most efficient velocity in one step, it would leave drag out of the picture during that first step, since drag is determined by the initial velocity at the begining of the loop. That's why I put the limitation in.
That's certainly possible... similar to pulse detonation engines...
Sorry, but no, my dad's name is not Roy. And I'm in Georgia, but Birmingham originally.
Dave
Join the Discussion
Have something to add? Share your thoughts — no account required.
Didn't find your answer?
Ask the community — no account required
Report Content
You are reporting this content to the moderators. They will look at it
ASAP.