import com.billkuker.rocketry.motorsim.Validating.ValidationException;\r
\r
public class Burn {\r
- //Some constants to tune adaptive regression step\r
- private static final double regStepIncreaseFactor = 1.01;\r
- private static final double regStepDecreaseFactor = .5;\r
- private static final Amount<Pressure> chamberPressureMaxDelta = Amount.valueOf(.5, SI.MEGA(SI.PASCAL));\r
+ private static Logger log = Logger.getLogger(Burn.class);\r
\r
- private static final Amount<Pressure> endPressure = Amount.valueOf(.1, RocketScience.PSI);\r
+ private static final BurnSettings settings = new BurnSettings();\r
+ public static final BurnSettings getBurnSettings(){\r
+ return settings;\r
+ }\r
+ \r
+ /**\r
+ * A class representing all the settigns one can change on a burn\r
+ * @author bkuker\r
+ */\r
+ public static class BurnSettings {\r
+ private BurnSettings(){};\r
+\r
+ public enum BurnVolumeMethod {\r
+ DeltaVolume,\r
+ SurfaceTimesRegression;\r
+ }\r
+ \r
+ private BurnVolumeMethod volumeMethod = BurnVolumeMethod.SurfaceTimesRegression;\r
+ private double regStepIncreaseFactor = 1.01;\r
+ private double regStepDecreaseFactor = .5;\r
+ private Amount<Pressure> chamberPressureMaxDelta = Amount.valueOf(.5, SI.MEGA(SI.PASCAL));\r
+ private Amount<Pressure> endPressure = Amount.valueOf(.1, RocketScience.PSI);\r
+ private Amount<Pressure> atmosphereicPressure = Amount.valueOf(101000, SI.PASCAL);\r
+ \r
+ public void setVolumeMethod(BurnVolumeMethod volumeMethod) {\r
+ this.volumeMethod = volumeMethod;\r
+ }\r
+ public BurnVolumeMethod getVolumeMethod() {\r
+ return volumeMethod;\r
+ }\r
+ public double getRegStepIncreaseFactor() {\r
+ return regStepIncreaseFactor;\r
+ }\r
+ public void setRegStepIncreaseFactor(double regStepIncreaseFactor) {\r
+ this.regStepIncreaseFactor = regStepIncreaseFactor;\r
+ }\r
+ public double getRegStepDecreaseFactor() {\r
+ return regStepDecreaseFactor;\r
+ }\r
+ public void setRegStepDecreaseFactor(double regStepDecreaseFactor) {\r
+ this.regStepDecreaseFactor = regStepDecreaseFactor;\r
+ }\r
+ public Amount<Pressure> getChamberPressureMaxDelta() {\r
+ return chamberPressureMaxDelta;\r
+ }\r
+ public void setChamberPressureMaxDelta(Amount<Pressure> chamberPressureMaxDelta) {\r
+ this.chamberPressureMaxDelta = chamberPressureMaxDelta;\r
+ }\r
+ public Amount<Pressure> getEndPressure() {\r
+ return endPressure;\r
+ }\r
+ public void setEndPressure(Amount<Pressure> endPressure) {\r
+ this.endPressure = endPressure;\r
+ }\r
+ public Amount<Pressure> getAtmosphereicPressure() {\r
+ return atmosphereicPressure;\r
+ }\r
+ public void setAtmosphereicPressure(Amount<Pressure> atmosphereicPressure) {\r
+ this.atmosphereicPressure = atmosphereicPressure;\r
+ }\r
+ }\r
\r
- private static Logger log = Logger.getLogger(Burn.class);\r
protected final Motor motor;\r
\r
private boolean burning = false;\r
public void setProgress(float p);\r
public void burnComplete();\r
}\r
- \r
+ \r
private Set<BurnProgressListener> bpls = new HashSet<Burn.BurnProgressListener>();\r
- \r
- private static final Amount<Pressure> atmosphereicPressure = Amount.valueOf(101000, SI.PASCAL);\r
- \r
+\r
public class Interval{\r
public Amount<Duration> time;\r
public Amount<Duration> dt;\r
\r
public void addBurnProgressListener( BurnProgressListener bpl ){\r
bpls.add(bpl);\r
+ if ( done )\r
+ bpl.burnComplete();\r
}\r
\r
public void burn(){\r
burning = true;\r
}\r
log.info("Starting burn...");\r
+ for (BurnProgressListener bpl : bpls ){\r
+ bpl.setProgress(0);\r
+ }\r
+ \r
int endPressureSteps = 0;\r
long start = new Date().getTime();\r
\r
initial.time = Amount.valueOf(0, SI.SECOND);\r
initial.dt = Amount.valueOf(0, SI.SECOND);\r
initial.regression = Amount.valueOf(0, SI.MILLIMETER);\r
- initial.chamberPressure = atmosphereicPressure;\r
+ initial.chamberPressure = settings.getAtmosphereicPressure();\r
initial.chamberProduct = Amount.valueOf(0, SI.KILOGRAM);\r
initial.thrust = Amount.valueOf(0, SI.NEWTON);\r
\r
step:\r
for ( int i = 0; i < 5000; i++ ) {\r
assert(positive(regStep));\r
- regStep = regStep.times(regStepIncreaseFactor);\r
+ regStep = regStep.times(settings.getRegStepIncreaseFactor());\r
\r
Interval prev = data.get(data.lastKey());\r
log.debug(prev);\r
\r
//log.debug("Vnew: " + motor.getGrain().volume(next.regression).to(SI.MILLIMETER.pow(3)));\r
\r
- //TODO Amount<Volume> volumeBurnt = motor.getGrain().volume(prev.regression).minus(motor.getGrain().volume(next.regression));\r
- Amount<Volume> volumeBurnt = motor.getGrain().surfaceArea(prev.regression).times(regStep).to(Volume.UNIT);\r
+ Amount<Volume> volumeBurnt;\r
+ if ( settings.getVolumeMethod() == BurnSettings.BurnVolumeMethod.DeltaVolume ){\r
+ volumeBurnt = motor.getGrain().volume(prev.regression).minus(motor.getGrain().volume(next.regression));\r
+ } else {\r
+ volumeBurnt = motor.getGrain().surfaceArea(prev.regression).times(regStep).to(Volume.UNIT);\r
+ }\r
assert(positive(volumeBurnt));\r
//log.info("Volume Burnt: " + volumeBurnt.to(SI.MILLIMETER.pow(3)));\r
\r
\r
Amount<MassFlowRate> mNozzle;\r
{\r
- Amount<Pressure> pDiff = prev.chamberPressure.minus(atmosphereicPressure);\r
+ Amount<Pressure> pDiff = prev.chamberPressure.minus(settings.getAtmosphereicPressure());\r
//log.debug("Pdiff: " + pDiff);\r
Amount<Area> aStar = motor.getNozzle().throatArea();\r
double k = motor.getFuel().getCombustionProduct().getRatioOfSpecificHeats();\r
//Product can not go negative!\r
if ( !positive(next.chamberProduct) ){\r
log.warn("ChamberProduct Negative on step " + i + "!, Adjusting regstep down and repeating step!");\r
- regStep = regStep.times(regStepDecreaseFactor);\r
+ regStep = regStep.times(settings.getRegStepDecreaseFactor());\r
continue step;\r
}\r
assert(positive(next.chamberProduct));\r
\r
log.debug("Product Density: " + combustionProductDensity);\r
\r
- next.chamberPressure = combustionProductDensity.times(specificGasConstant).times(chamberTemp).plus(atmosphereicPressure).to(Pressure.UNIT);\r
+ next.chamberPressure = combustionProductDensity.times(specificGasConstant).times(chamberTemp).plus(settings.getAtmosphereicPressure()).to(Pressure.UNIT);\r
assert(positive(next.chamberPressure));\r
\r
next.chamberPressure = Amount.valueOf(\r
SI.PASCAL);\r
\r
Amount<Pressure> dp = next.chamberPressure.minus(prev.chamberPressure);\r
- if ( dp.abs().isGreaterThan(chamberPressureMaxDelta)){\r
+ if ( dp.abs().isGreaterThan(settings.getChamberPressureMaxDelta())){\r
log.warn("DP " + dp + " too big!, Adjusting regstep down and repeating step!");\r
- regStep = regStep.times(regStepDecreaseFactor);\r
+ regStep = regStep.times(settings.getRegStepDecreaseFactor());\r
continue step;\r
}\r
\r
- next.thrust = motor.getNozzle().thrust(next.chamberPressure, atmosphereicPressure, atmosphereicPressure, motor.getFuel().getCombustionProduct().getRatioOfSpecificHeats2Phase());\r
+ next.thrust = motor.getNozzle().thrust(next.chamberPressure, settings.getAtmosphereicPressure(), settings.getAtmosphereicPressure(), motor.getFuel().getCombustionProduct().getRatioOfSpecificHeats2Phase());\r
assert(positive(next.thrust));\r
\r
- if ( i > 100 && next.chamberPressure.minus(atmosphereicPressure).abs().isLessThan(endPressure)){\r
+ if ( i > 100 && next.chamberPressure.minus(settings.getAtmosphereicPressure()).abs().isLessThan(settings.getEndPressure())){\r
log.info("Pressure at ~Patm on step " + i);\r
endPressureSteps++;\r
if ( endPressureSteps > 5 )\r