public void SimulateAeroProperties(out Vector3 aeroForce, out Vector3 aeroTorque, Vector3 velocityWorldVector, double altitude)
{
// Rodhern: It seems that this method, 'SimulateAeroProperties', is only used in FARAPI, which in turn can be used by say KSPTrajectories.
// The parameter 'FARCenterQuery dummy' is from a code fix by Benjamin Chung (commit 18fbb9d29431679a4de9dfc22a443f400d2d4f8b).
FARCenterQuery center = new FARCenterQuery();
FARCenterQuery dummy = new FARCenterQuery();
float pressure;
float density;
float temperature;
float speedOfSound;
CelestialBody body = vessel.mainBody; //Calculate main gas properties
pressure = (float)body.GetPressure(altitude);
temperature = (float)body.GetTemperature(altitude);
density = (float)body.GetDensity(pressure, temperature);
speedOfSound = (float)body.GetSpeedOfSound(pressure, density);
if (pressure <= 0 || temperature <= 0 || density <= 0 || speedOfSound <= 0)
{
aeroForce = Vector3.zero;
aeroTorque = Vector3.zero;
return;
}
float velocityMag = velocityWorldVector.magnitude;
float machNumber = velocityMag / speedOfSound;
float reynoldsNumber = (float)FARAeroUtil.CalculateReynoldsNumber(density, Length, velocityMag, machNumber, temperature, body.atmosphereAdiabaticIndex);
float reynoldsPerLength = reynoldsNumber / (float)Length;
float pseudoKnudsenNumber = machNumber / (reynoldsNumber + machNumber);
float skinFriction = (float)FARAeroUtil.SkinFrictionDrag(reynoldsNumber, machNumber);
FlightEnv fenv = FlightEnv.NewPredicted(vessel.mainBody, altitude, machNumber);
if (_currentAeroSections != null)
{
for (int i = 0; i < _currentAeroSections.Count; i++)
{
FARAeroSection curSection = _currentAeroSections[i];
if (curSection != null)
{
curSection.PredictionCalculateAeroForces(density, machNumber, reynoldsPerLength, pseudoKnudsenNumber, skinFriction, velocityWorldVector, center);
}
}
for (int i = 0; i < _legacyWingModels.Count; i++)
{
FARWingAerodynamicModel curWing = _legacyWingModels[i];
if ((object)curWing != null)
{
Vector3d force = curWing.PrecomputeCenterOfLift(velocityWorldVector, fenv, dummy);
center.AddForce(curWing.transform.position, force);
}
}
}
aeroForce = center.force;
aeroTorque = center.TorqueAt(vessel.CoM);
}
public void SimulateAeroProperties(out Vector3 aeroForce, out Vector3 aeroTorque, Vector3 velocityWorldVector, double altitude)
{
FARCenterQuery center = new FARCenterQuery();
float pressure;
float density;
float temperature;
float speedOfSound;
CelestialBody body = vessel.mainBody; //Calculate main gas properties
pressure = (float)body.GetPressure(altitude);
temperature = (float)body.GetTemperature(altitude);
density = (float)body.GetDensity(pressure, temperature);
speedOfSound = (float)body.GetSpeedOfSound(pressure, density);
if (pressure <= 0 || temperature <= 0 || density <= 0 || speedOfSound <= 0)
{
aeroForce = Vector3.zero;
aeroTorque = Vector3.zero;
return;
}
float velocityMag = velocityWorldVector.magnitude;
float machNumber = velocityMag / speedOfSound;
float reynoldsNumber = (float)FARAeroUtil.CalculateReynoldsNumber(density, Length, velocityMag, machNumber, temperature, body.atmosphereAdiabaticIndex);
float reynoldsPerLength = reynoldsNumber / (float)Length;
float skinFriction = (float)FARAeroUtil.SkinFrictionDrag(reynoldsNumber, machNumber);
float pseudoKnudsenNumber = machNumber / (reynoldsNumber + machNumber);
if (_currentAeroSections != null)
{
for (int i = 0; i < _currentAeroSections.Count; i++)
{
FARAeroSection curSection = _currentAeroSections[i];
if (curSection != null)
{
curSection.PredictionCalculateAeroForces(density, machNumber, reynoldsPerLength, pseudoKnudsenNumber, skinFriction, velocityWorldVector, center);
}
}
for (int i = 0; i < _legacyWingModels.Count; i++)
{
FARWingAerodynamicModel curWing = _legacyWingModels[i];
if ((object)curWing != null)
{
curWing.PrecomputeCenterOfLift(velocityWorldVector, machNumber, density, center);
}
}
}
aeroForce = center.force;
aeroTorque = center.TorqueAt(vessel.CoM);
}