hongshaorou/openrocket
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Avoid divide by 0 in TubeCalc when velocity is 0 by setting CD to 0 in this case.

Browse Source
This commit is contained in:
JoePfeiffer 2023-10-26 14:06:42 -06:00
parent 1690f96c3f
commit 01f1dd1dc6
1 changed files with 5 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
core/src/net/sf/openrocket/aerodynamics/barrowman/TubeCalc.java
View File

@ -40,6 +40,11 @@ public abstract class TubeCalc extends RocketComponentCalc {
public double calculatePressureCD(FlightConditions conditions, public double calculatePressureCD(FlightConditions conditions,
double stagnationCD, double baseCD, WarningSet warnings) { double stagnationCD, double baseCD, WarningSet warnings) {
// If we aren't moving, treat CD as 0
final double v = conditions.getVelocity();
if (v < MathUtil.EPSILON)
return 0;
// Need to check for tube inner area 0 in case of rockets using launch lugs with // Need to check for tube inner area 0 in case of rockets using launch lugs with
// an inner radius of 0 to emulate rail guides (or just weird rockets, of course) // an inner radius of 0 to emulate rail guides (or just weird rockets, of course)
double tubeCD = 0.0; double tubeCD = 0.0;
@ -49,7 +54,6 @@ public abstract class TubeCalc extends RocketComponentCalc {
final double p = conditions.getAtmosphericConditions().getPressure(); final double p = conditions.getAtmosphericConditions().getPressure();
final double t = conditions.getAtmosphericConditions().getTemperature(); final double t = conditions.getAtmosphericConditions().getTemperature();
final double rho = conditions.getAtmosphericConditions().getDensity(); final double rho = conditions.getAtmosphericConditions().getDensity();
final double v = conditions.getVelocity();
// Reynolds number (note Reynolds number for the interior of a pipe is based on diameter, // Reynolds number (note Reynolds number for the interior of a pipe is based on diameter,
// not length (t)) // not length (t))