fly/flightgear
1
0
Fork 0
Code Activity

Latest JSBSim changes -- EGT fixes.

Browse source
This commit is contained in:
david 2001-12-15 13:53:53 +00:00
parent 0b5594d50d
commit 0505201362
1 changed files with 15 additions and 12 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

27
src/FDM/JSBSim/FGPiston.cpp
View file

@ -396,21 +396,24 @@ void FGPiston::doEnginePower(void)
void FGPiston::doEGT(void) void FGPiston::doEGT(void)
{ {
double delta_T_exhaust = 0.0; double delta_T_exhaust;
double heat_capacity_exhaust;
double enthalpy_exhaust; double enthalpy_exhaust;
double heat_capacity_exhaust;
double dEGTdt;
combustion_efficiency = Lookup_Combustion_Efficiency->GetValue(equivalence_ratio); if ((Running) && (m_dot_air > 0.0)) { // do the energy balance
enthalpy_exhaust = m_dot_fuel * calorific_value_fuel * combustion_efficiency * 0.33; combustion_efficiency = Lookup_Combustion_Efficiency->GetValue(equivalence_ratio);
heat_capacity_exhaust = (Cp_air * m_dot_air) + (Cp_fuel * m_dot_fuel); enthalpy_exhaust = m_dot_fuel * calorific_value_fuel *
combustion_efficiency * 0.33;
if (heat_capacity_exhaust >= 0.0000001) heat_capacity_exhaust = (Cp_air * m_dot_air) + (Cp_fuel * m_dot_fuel);
delta_T_exhaust = enthalpy_exhaust / heat_capacity_exhaust; delta_T_exhaust = enthalpy_exhaust / heat_capacity_exhaust;
else ExhaustGasTemp_degK = T_amb + delta_T_exhaust;
delta_T_exhaust = 0.0; ExhaustGasTemp_degK *= 0.444 + ((0.544 - 0.444) * Percentage_Power / 100.0);
} else { // Drop towards ambient - guess an appropriate time constant for now
ExhaustGasTemp_degK = T_amb + delta_T_exhaust; dEGTdt = (298.0 - ExhaustGasTemp_degK) / 100.0;
ExhaustGasTemp_degK *= 0.444 + ((0.544 - 0.444) * Percentage_Power / 100.0); delta_T_exhaust = dEGTdt * dt;
ExhaustGasTemp_degK += delta_T_exhaust;
}
} }
//%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% //%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%