modelGasolineEngineChp_EngineModel

Internal combustion engine model for CHP-applications.

Information

The model of the internal combustion engine is the centrepiece of the BHKW model developed. It is divided into a thermal and a mechanical-chemical part.

The energy balance of the combustion engine is used to determine the maximum heat released. The mechanical power is calculated using a mean value model, which is characterized by short calculation times and simple parameterization. An empirical approach is used to describe the specific enthalpy of the exhaust gas. The required mass flows of the combustion products are determined by assuming complete combustion with a known fuel composition. Various liquid and gaseous fuels have been implemented on the basis of the existing material models and are available to the user.

Due to the complexity of a combustion engine, assumptions have to be made. These are listed below to understand the function and applicability of the model.


Assumptions


Assumptions made and resulting limitations of the internal combustion engine model:

- The nominal power point of the power unit is known and the modulationof the operation point is achieved by a reduction of the introduced fuel.

- The indicated mean pressure is assumed to be constant as a necessary measure for the calculation of engine power. This corresponds to a constant thermodynamic combustion process within the cylinders.

- The engine must be started with an electric machine. It is controlled by the release of the fuel quantity from a minimum speed (800rpm). The speed then increases to equilibrium with the counteracting generator torque.

- Complete and superstoichiometric combustion is assumed to solve the gross reaction equation

- Entry of air and fuel at ambient conditions and constant amount of fuel and air per combustion cycle

-> Only conditionally with turbocharging of the engines, since then the cylinder filling can vary depending on the boost pressure (slight consideration due to stored rated performance data)


Air ratio or residual oxygen in the exhaust gas is known to estimate the combustion process. So this is a necessary assumption for the calculation of material flows (mass flows, composition of the exhaust gas).

The mean specific heat capacity of the exhaust gas for a temperature range from 0 ° C to the maximum adiabatic combustion temperature is used to calculate the exhaust gas temperature.

The mean specific heat capacity is determinated with a potency approach according to Müller (1968).

Frictional losses that can be calculated based on a known friction mean pressure at a speed of 3000rpm (if not known, default average values ​​from VK1 by S.Pischinger) are converted into usable heat.

The calculation of the exhaust gas enthalpy according to an empirical approach is based on investigations by R.Pischinger which uses a reference point temperature of 25 ° C (initial state of the reaction educts from combustion air and fuel).

-> Consideration of the chemical and thermal proportions of the enthalpy

-> Limited accuracy for diesel engine (non-premixed) processes

Parameters

TypeNameDefaultDescription
Modelica.Units.SI.VolumeVCylCHPEngData.VEng/CHPEngData.zCylinder displacement
RotationSpeednEngNominal25.583Nominal engine speed at operating point
Modelica.Units.SI.PowerP_mecNominalCHPEngData.P_mecNominalMecanical power output at nominal operating point
Modelica.Units.SI.TemperatureT_Amb298.15Ambient temperature (matches to fuel and combustion air temperature)
GasConstantR8.31446Gasconstant for calculation purposes
RealQuoDCylCHPEngData.QuoDCyl
BooleanFuelTypeMedium1.isGasTrue = Gasoline fuel, False = Liquid fuel
Modelica.Units.SI.MassFlowRatem_MaxExhCHPEngData.P_FueNominal/H_U*(1 + Lambda*L_St)Maximal exhaust gas flow based on the fuel and combustion properties
Modelica.Units.SI.Massm_FueEngRotCHPEngData.P_FueNominal*60/(H_U*CHPEngData.nEngMax*CHPEngData.i)Injected fuel mass per engine rotation(presumed as constant)
Modelica.Units.SI.Pressurep_Amb101325Ambient pressure
Modelica.Units.SI.Pressurep_mip_mfNominal + p_meNominalConstant indicated mean effective cylinder pressure
Modelica.Units.SI.Pressurep_meNominalCHPEngData.p_meNominalNominal mean effective cylinder pressure
Modelica.Units.SI.Pressureref_p_mfNominalCHPEngData.ref_p_mfNominalFriction mean pressure of reference engine for calculation(dCyl=91mm & nEng=3000rpm & TEng=90°C)
Modelica.Units.SI.Pressurep_mfNominalref_p_mfNominal*QuoDCyl^(-0.3)Nominal friction mean pressure
Modelica.Units.SI.TemperatureT_ExhOutCHPEngData.T_ExhPowUniOutAssumed exhaust gas outlet temperature of the CHP unit for heat calculations
Modelica.Units.SI.SpecificEnergyH_UMedium1.H_USpecific calorific value of the fuel
RealLambdaCHPEngData.LambdaCombustion air ratio
RealL_StMedium1.L_stStoichiometric air consumption per mass fuel
Reall_MinL_St*MM_Fuel/MM_AirMinimum molar air consumption per mole fuel
Modelica.Units.SI.MolarMassMM_FuelMedium1.MMMolar mass of the fuel
Modelica.Units.SI.MolarMassMM_AirMedium2.MMMolar mass of the combustion air
Modelica.Units.SI.MolarMass[:]MM_ComExhMedium3.data[:].MMMolar masses of the combustion products: N2, O2, H2O, CO2
Real[:]expFacCpComExh{0.11, 0.15, 0.20, 0.30}Exponential factor for calculating the specific heat capacity of N2, O2, H2O, CO2
Modelica.Units.SI.SpecificHeatCapacity[:]cpRefComExh{1000, 900, 1750, 840}Specific heat capacities of the combustion products at reference state at 0°C
Modelica.Units.SI.TemperatureRefT_Com1473.15Reference combustion temperature for calculation purposes
Realn_N2Exhif FuelType then Medium1.moleFractions_Gas[1] + Lambda*l_Min*Medium2.moleFractions_Air[1] else Lambda*l_Min*Medium2.moleFractions_Air[1]Exhaust: Number of molecules Nitrogen per mole of fuel
Realn_O2Exh(Lambda - 1)*l_Min*Medium2.moleFractions_Air[2]Exhaust: Number of molecules Oxygen per mole of fuel
Realn_H2OExhif FuelType then 0.5*sum(Medium1.moleFractions_Gas[i]*Medium1.Fuel.nue_H[i] for i in 1:size(Medium1.Fuel.nue_H, 1)) else 0.5*(Medium1.Fuel.Xi_liq[2]*Medium1.MM/Medium1.Fuel.MMi_liq[2])Exhaust: Number of molecules H20 per mole of fuel
Realn_CO2Exhif FuelType then sum(Medium1.moleFractions_Gas[i]*Medium1.Fuel.nue_C[i] for i in 1:size(Medium1.Fuel.nue_C, 1)) else Medium1.Fuel.Xi_liq[1]*Medium1.MM/Medium1.Fuel.MMi_liq[1]Exhaust: Number of molecules CO2 per mole of fuel
Real[:]n_ComExh{n_N2Exh, n_O2Exh, n_H2OExh, n_CO2Exh}
Realn_Exhsum(n_ComExh[j] for j in 1:size(n_ComExh, 1))Number of exhaust gas molecules per mole of fuel
Modelica.Units.SI.MolarMassMM_Exhsum(n_ComExh[i]*MM_ComExh[i] for i in 1:size(n_ComExh, 1))/sum(n_ComExh[i] for i in 1:size(n_ComExh, 1))Molar mass of the exhaust gas
Modelica.Units.SI.MassFractionX_N2ExhMM_ComExh[1]*n_ComExh[1]/(MM_Exh*n_Exh)Mass fraction of N2 in the exhaust gas
Modelica.Units.SI.MassFractionX_O2ExhMM_ComExh[2]*n_ComExh[2]/(MM_Exh*n_Exh)Mass fraction of O2 in the exhaust gas
Modelica.Units.SI.MassFractionX_H2OExhMM_ComExh[3]*n_ComExh[3]/(MM_Exh*n_Exh)Mass fraction of H2O in the exhaust gas
Modelica.Units.SI.MassFractionX_CO2ExhMM_ComExh[4]*n_ComExh[4]/(MM_Exh*n_Exh)Mass fraction of CO2 in the exhaust gas
Modelica.Units.SI.MassFraction[size(n_ComExh, 1)]Xi_Exh{X_N2Exh, X_O2Exh, X_H2OExh, X_CO2Exh}
Unit properties
AixLib.DataBase.CHP.ModularCHPEngineData.CHPEngDataBaseRecordCHPEngDataDataBase.CHP.ModularCHPEngineData.CHP_SenerTecDachsG5_5()Needed engine data for calculations

Connectors

TypeNameDefaultDescription
Modelica.Fluid.Interfaces.FluidPort_bport_exh
Modelica.Mechanics.Rotational.Interfaces.Flange_aflange_aMechanical port of the engines output drive
Modelica.Blocks.Interfaces.RealInputexhaustGasTemperatureExhaust gas temperature calculated in the thermal engine model
AixLib.Controls.Interfaces.CHPControlBuscHPEngineBusSignal bus of the mechanical combustion engine model

Components

TypeNameDefaultDescription
BooleanSwitchOnOfftrueOperation switch of the CHP unit (true=On, false=Off)
RotationSpeednEngCurrent engine speed
Modelica.Units.SI.MassFlowRatem_flow_ExhMass flow rate of exhaust gas
Modelica.Units.SI.MassFlowRatem_flow_CO2ExhMass flow rate of CO2 in the exhaust gas
Modelica.Units.SI.MassFlowRatem_flow_FueMass flow rate of fuel
Modelica.Units.SI.MassFlowRatem_flow_AirMass flow rate of combustion air
Modelica.Units.SI.SpecificHeatCapacity[size(n_ComExh, 1)]meanCpComExhCalculated specific heat capacities of the exhaust gas components for the calculated combustion temperature
Modelica.Units.SI.SpecificHeatCapacitymeanCpExhCalculated specific heat capacity of the exhaust gas for the calculated combustion temperature
Modelica.Units.SI.SpecificEnergyh_Exh1000*(-286 + 1.011*T_ExhCHPOut - 27.29*Lambda + 0.000136*T_ExhCHPOut^2 - 0.0255*T_ExhCHPOut*Lambda + 6.425*Lambda^2)Specific enthalpy of the exhaust gas
Modelica.Units.SI.PowerP_effEffective(mechanical) engine power
Modelica.Units.SI.PowerP_Fuem_flow_Fue*H_UFuel expenses at operating point
Modelica.Units.SI.PowerH_ExhEnthalpy stream of the exhaust gas
Modelica.Units.SI.PowerCalQ_thermCalculated heat from engine combustion
Modelica.Units.SI.PowerQ_thermTotal heat from engine combustion
Modelica.Units.SI.TorqueMmotCalculated engine torque
Modelica.Units.SI.TemperatureT_logEngCool356.15Logarithmic mean temperature of coolant inside the engine
Modelica.Units.SI.TemperatureT_ComTemperature of the combustion gases
Modelica.Units.SI.TemperatureT_ExhCHPOut383.15Exhaust gas outlet temperature of CHP unit
RealmodFac1Modulation factor for energy outuput control of the Chp unit
RealA01.0895 - 1.079*10^(-2)*(T_logEngCool - 273.15) + 5.525*10^(-5)*(T_logEngCool - 273.15)^2
RealA14.68*10^(-4) - 5.904*10^(-6)*(T_logEngCool - 273.15) + 1.88*10^(-8)*(T_logEngCool - 273.15)^2
RealA2-4.35*10^(-8) + 1.12*10^(-9)*(T_logEngCool - 273.15) - 4.79*10^(-12)*(T_logEngCool - 273.15)^2
RealB0-2.625*10^(-3) + 3.75*10^(-7)*(nEng*60) + 1.75*10^(-5)*(T_logEngCool - 273.15) + 2.5*10^(-9)*(T_logEngCool - 273.15)*(nEng*60)
RealB18.95*10^(-3) + 1.5*10^(-7)*(nEng*60) + 7*10^(-6)*(T_logEngCool - 273.15) - 10^(-9)*(T_logEngCool - 273.15)*(nEng*60)
Modelica.Units.SI.Pressurep_mfp_mfNominal*((A0 + A1*(nEng*60) + A2*(nEng*60)^2) + (B0 + B1*(p_meNominal/100000)))Current friction mean pressure at operating point
Modelica.Units.SI.Pressurep_me(modFac*p_mi) - p_mfCurrent mean effective pressure at operating point
RealetaMecp_me/p_miCurrent percentage of usable mechanical power compared to inner cylinder power from combustion
Modelica.Fluid.Sources.MassFlowSource_TexhaustFlowFlow source of the exhaust gas from the engine combustion
Modelica.Blocks.Sources.RealExpressionmassFlowExhaust
Modelica.Blocks.Sources.RealExpressioneffectiveMechanicalTorque
Modelica.Mechanics.Rotational.Sources.TorqueengineTorqueCalculated engine torque
Modelica.Mechanics.Rotational.Components.InertiainertiaInertia model of the engine
Modelica.Blocks.Sources.RealExpressionengineSpeed
Modelica.Blocks.Sources.RealExpressionfuelPower
Modelica.Blocks.Sources.RealExpressionthermalPower
Modelica.Blocks.Sources.RealExpressionoutputTorque
Modelica.Blocks.Sources.RealExpressionairFlow
Modelica.Blocks.Sources.RealExpressioncarbondioxideFlow
Modelica.Blocks.Sources.RealExpressionSpecificHeatCapExhGas
Modelica.Blocks.Sources.RealExpressionfuelFlow

Contents

NameDescription
Medium1
Medium2
Medium3
RotationSpeed
GasConstant

Revisions

  • April, 2019  by Julian Matthes:
    First implementation (see issue #667)