WOLFRAM SYSTEM MODELER

TurbojetEngine

Turbojet engine: unparameterized

Diagram

Wolfram Language

In[1]:=
SystemModel["Aircraft.Physical.FixedWing.Parts.Propulsions.Engines.TurbojetEngines.TurbojetEngine"]
Out[1]:=

Information

This turbojet engine model extends the EngineBase model and models the thrust, drag and mass properties (if weight estimation method is used) of a turbojet engine.

Parameters

The required parameters for the turbojet engine model are listed in the Propulsion tab, as shown in Figure 1.


Figure 1: Required input parameters for turbojet engine model.

Should the value of any of the estimated properties listed in the Estimated Properties parameter tab be known (such as engine mass, nacelle length or diameter), the equations can be bypassed by entering the known value directly to the respective parameter input field. Similarly, if the surface roughness height of the nacelles is known, its value can be entered inside the Aerodynamic Coefficients tab. Otherwise its value is propagated from the global surface roughness height declared in AircraftBase.

The remaining parameters listed in the Propagated Properties tab are propagated from the AircraftBase as the propulsion system is compiled in TurbojetPropulsion, and thus they can and should be left unchanged.

Thrust Available

The thrust available (Tavail) for turbojet engines is solved through the method presented by Gudmundsson [1], which defines the value of Tavail  as a function of the engine static thrust (Tstatic), Mach number (M), temperature ratio (θ), pressure ratio (δ) and throttle ratio (TR) as

.

The derivation of θ, δ and TR can be found in section 3.2 in Reference [2]. In short, it can be concluded that the thrust available decreases as the Mach number and altitude increase.

Fuel Consumption

The parameter for thrust-specific fuel consumption (TSFC) is propagated to the energySystem component, where fuel consumption calculations are performed.

Drag Due to Nacelle

For solving the parasite drag coefficient (CD,0,nac) of different components in the aircraft, including the engine nacelles, the component buildup method presented by Raymer [3] is used, defined as

,

where the nacelle form factor (FFnac) and the nacelle wetted area (Swet,nac) are functions of the estimated (or entered, if known) values of nacelle diameter and length. The nacelle dimensions are estimated through empirical relationships as a function of the engine static thrust. Thus, the more powerful the engine, the larger the parasite drag coefficient.

The skin friction coefficient (Cf,nac) is a function of the nacelle surface roughness height, Mach number and the Reynolds number for the flow over the nacelle, and thus the air compressibility effects are included in the drag calculations. The complete derivation of the CD,0,nac can be found in section 3.3.1 in Reference [2].

The magnitude of the drag force itself (Dnac) is the parasite drag coefficient times the dynamic pressure (q) times the main wing reference area (Sref,w)

.

In this model, only the parasite drag of the nacelle is considered, and thus the drag force generated by the nacelle is invariable with the angle of attack. Consequently, the lift generated by the nacelle due to angle of attack is also omitted.

Mass Properties

If the weight estimation method is used and no mass properties are entered by the user, the engine mass is solved by using an empirical relationship based on the given Tstatic parameter value, as presented by Isikveren [4]. Additionally, the existence of a thrust reverser, here triggered by entering a nonzero negThrust parameter value, will slightly increase the engine mass.

The engine moments of inertia are estimated as if the engine were a solid cylinder with its estimated length and diameter. Hence, the center of mass of the engine is also estimated to be at the geometric center of this cylinder.

The full derivations of the estimated engine mass and moments of inertia are also shown in Reference [2] in sections 3.4.4 and 3.6.1, respectively. 

References

[1]  Gudmundsson, S. (2016). General Aviation Aircraft Design. Elsevier.

[2]  Erä-Esko, N. (2022). "Development and Use of System Modeler 6DOF Flight Mechanics Model in Aircraft Conceptual Design."
      Available atmodelica://Aircraft/Resources/Documents/EraeEskoThesis.pdf.

[3]  Raymer, D. P. (1992). Aircraft Design: A Conceptual Approach, 2nd Ed. American Institute of Aeronautics and Astronautics.

[4]  Isikveren, A. T. (2002). "Quasi-analytical Modelling and Optimisation Techniques for Transport Aircraft Design," Thesis, Institutionen för
      flygteknik, Department of Aeronautics.

Parameters (26)

weightEst

Value:

Type: Boolean

Description: true, if weight estimation method is used for masses, center of mass location and inertia tensor

bWing

Value:

Type: Length (m)

Description: Main wing span

SrefWing

Value:

Type: Area (m²)

Description: Main wing reference area

CADshapes

Value:

Type: Boolean

Description: true, if external CAD files are used for animation

negThrust

Value:

Type: Real

Description: Maximum negative thrust (0 to 1 of thrust available)

T0

Value:

Type: Temperature (K)

Description: Temperature at sea-level

p0

Value:

Type: Pressure (Pa)

Description: Static pressure at sea-level

gammaAir

Value:

Type: Real

Description: Adiabatic index for air

Tstatic

Value:

Type: Force (N)

Description: Static thrust of one engine at sea level

TSFC

Value:

Type: ThrustSpecificFuelConsumption (kg⋅N⁻¹⋅s⁻¹)

Description: Thrust-specific fuel consumption in cruise

altEngDes

Value:

Type: Height (m)

Description: Engine design altitude

machEngDes

Value:

Type: Real

Description: Engine design Mach number

kSkinNac

Value:

Type: Length (m)

Description: Nacelle surface roughness height

mNac

Value: if negThrust > 0 then 0.345 * 1.18 * mEngDry else 0.345 * mEngDry

Type: Mass (kg)

Description: Nacelle mass

mPyl

Value: 0.574 * mEngDry ^ 0.736

Type: Mass (kg)

Description: Pylon mass

mEngDry

Value: 0.0117 * Tstatic ^ 1.0572

Type: Mass (kg)

Description: Engine dry mass

mEngInstalled

Value: mEngDry + mNac + mPyl

Type: Mass (kg)

Description: Installed engine mass

xEngCM

Value: lNac / 2

Type: Length (m)

Description: Engine center of mass x-coordinate w.r.t. engine rear end

IxxEng

Value: mEngInstalled * (dNac / 2) ^ 2 / 2

Type: MomentOfInertia (kg⋅m²)

Description: Engine roll moment of inertia

IyyEng

Value: mEngInstalled * (dNac / 2) ^ 2 / 4 + mEngInstalled * lNac ^ 2 / 12

Type: MomentOfInertia (kg⋅m²)

Description: Engine pitch moment of inertia

IzzEng

Value: mEngInstalled * (dNac / 2) ^ 2 / 4 + mEngInstalled * lNac ^ 2 / 12

Type: MomentOfInertia (kg⋅m²)

Description: Engine yaw moment of inertia

dNac

Value: 4 * (0.0625 + 1 / (4 * sqrt(2)) * max(1.730 * log(max(Tstatic / 1000, 1)) - Modelica.Constants.pi, 0) ^ 0.5)

Type: Length (m)

Description: Nacelle diameter

lNac

Value: 5 * (Tstatic / 1000) ^ 0.9839 / (6 * Modelica.Constants.pi * (1.730 * log(max(Tstatic / 1000, 1)) - Modelica.Constants.pi))

Type: Length (m)

Description: Nacelle length

SwetNac

Value: 2 * Modelica.Constants.pi ^ 2 * 0.2028 * dNac * ((0.20571 * lNac ^ 2 + 0.04661 * dNac ^ 2) ^ 0.5 + (0.1853 * lNac ^ 2 + 0.07557 * dNac ^ 2) ^ 0.5 - (0.005077 * lNac ^ 2 + 0.01611 * dNac ^ 2) ^ 0.5 - (0.01651 * lNac ^ 2 + 0.03666 * dNac ^ 2) ^ 0.5)

Type: Area (m²)

Description: Nacelle wetted area

FFnac

Value: 1.17 * (1 + 0.35 * (dNac / max(lNac, 0.01)))

Type: Real

Description: Nacelle form factor

throtR

Value: TengDes / T0 * (1 + (gammaAir - 1) / 2 * machEngDes ^ 2)

Type: Real

Description: Throttle ratio

Inputs (1)

flightData

Type: FlightData

Description: Global flight data variables

Connectors (3)

deltaThrotCmd

Type: RealInput

Description: Throttle position command

engineRP

Type: Frame_b

Description: Connector to engine reference point

yConsumption

Type: RealOutput

Description: Consumed electric energy or fuel in an engine

Components (15)

flightData

Type: FlightData

Description: Global flight data variables

energySystem

Type: FuelSystemThrustSpecific

Description: Model to calculate the consumed electric energy or fuel in an engine

thrustEngine

Type: WorldForce

Description: Engine net thrust

engShape

Type: FixedShape

Description: Engine visualization

thrustAvailable

Type: RealExpression

Description: Thrust available of one engine

engineThrustDmd

Type: Product

Description: Thrust demand for one engine

engineDynamics

Type: CriticalDamping

Description: Simplified model of the engine dynamics

throttleOnGround

Type: Limiter

Description: Limits the throttle from maximum negative thrust to 1 once on the ground

dragNacelle3D

Type: RealExpression[3]

Description: Drag due to one nacelle

dragNacelle

Type: WorldForce

Description: Drag of nacelle

translEng

Type: FixedTranslation

Description: Translation from engine rear end to the point where thrust and drag are acting on (currently no translation)

bodyEng

Type: Body

Description: Mass and inertia of engine

switch

Type: Switch

Description: Switch between two Real signals

throttleInAir

Type: Limiter

Description: Limits the throttle from 0 to 1 once airborne

onGround

Type: BooleanExpression

Description: Boolean expression for indicating aircraft being on ground

Extended by (2)

RollsRoyceRB16286

Aircraft.Physical.FixedWing.Parts.Propulsions.Engines.TurbojetEngines

Turbojet engine: Rolls-Royce RB.162-86

PrattWhitneyJT4A11

Aircraft.Physical.FixedWing.Parts.Propulsions.Engines.TurbojetEngines

Turbojet engine: Pratt & Whitney JT4A-11