From b9bd32cca226658a1405b6a435ad45b1fea26f31 Mon Sep 17 00:00:00 2001 From: Bill Kuker Date: Sun, 22 Jan 2012 19:41:30 +0000 Subject: [PATCH] Track fuel burnt at each step --- src/com/billkuker/rocketry/motorsim/Burn.java | 6 +++++- 1 file changed, 5 insertions(+), 1 deletion(-) diff --git a/src/com/billkuker/rocketry/motorsim/Burn.java b/src/com/billkuker/rocketry/motorsim/Burn.java index cea51d5..a87ba62 100644 --- a/src/com/billkuker/rocketry/motorsim/Burn.java +++ b/src/com/billkuker/rocketry/motorsim/Burn.java @@ -113,6 +113,7 @@ public class Burn { public Amount chamberPressure; Amount chamberProduct; public Amount thrust; + public Amount fuelBurnt; public String toString(){ return time + " " + dt + " " + regression + " " + chamberPressure + " " + chamberProduct; @@ -173,6 +174,7 @@ public class Burn { initial.chamberPressure = settings.getAtmosphereicPressure(); initial.chamberProduct = Amount.valueOf(0, SI.KILOGRAM); initial.thrust = Amount.valueOf(0, SI.NEWTON); + initial.fuelBurnt = Amount.valueOf(0, SI.KILOGRAM); data.put(Amount.valueOf(0, SI.SECOND), initial); @@ -226,7 +228,9 @@ public class Burn { assert(positive(volumeBurnt)); //log.info("Volume Burnt: " + volumeBurnt.to(SI.MILLIMETER.pow(3))); - Amount mGenRate = volumeBurnt.times(motor.getFuel().getIdealDensity().times(motor.getFuel().getDensityRatio())).divide(dt).to(MassFlowRate.UNIT); + Amount massBurnt = volumeBurnt.times(motor.getFuel().getIdealDensity().times(motor.getFuel().getDensityRatio())).to(Mass.UNIT); + next.fuelBurnt = massBurnt; + Amount mGenRate = massBurnt.divide(dt).to(MassFlowRate.UNIT); assert(positive(mGenRate)); //log.debug("Mass Gen Rate: " + mGenRate); -- 2.47.2