Main Content

Air Intake

Combustion cylinder air supply intake

  • Library:
  • Simscape / Driveline / Engines & Motors / Engine Subcomponents

Description

The Air Intake block represents an abstract air intake system that computes lumped-parameter air flow dynamics. This block is suitable for internal combustion engines. The block computes the intake manifold pressure from the throttle and engine speed.

Equations

The Air Intake block computes the intake manifold pressure as

ddtPim=RairTambVim(m˙thrm˙cyl),

where

  • Rair is the Dry air specific gas constant parameter.

  • Tamb is the Ambient temperature parameter.

  • Vim is the Intake manifold volume parameter.

The block calculates the combustion cylinder mass flow rate, cyl, as

m˙cyl=ωVdPim2RairTamb,

which is a time derivative formulation of the ideal gas law. ω is the engine speed the block receives from port EngSpd, and Vd is the Engine displacement volume parameter.

The block calculates the throttle mass flow rate, thr as

m˙thr={SthrTambc1(PambPim)                            LaminarSthrTambPambc2[(PambPim)cb(PambPim)ca]  SubsonicSthrTambc3Pamb                                        Choked

depending on the flow conditions. In this equation, Sthr is the area of the throttle opening, and Pamb is the ambient pressure. The block defines the equation constants as

ca=κ+1κcb=2κc1=c2(ΠlamcbΠlamca)1Πlamc2=2κRair(κ1)c3=κRairΠcrca

where

  • κ is the Dry air specific heat ratio parameter.

  • Πlam is the laminar flow threshold constant.

  • Πcr is the critical pressure ratio.

The block finds Πcr as

Πcr=(2κ+1)κκ1.

Throttle Open Area

The block treats the throttle open area, Sthr, as a third-order polynomial that functions with respect to the throttle angle, θ. This method simplifies the block computation to improve performance. The block computes the throttle open area as

Sthr(θ)=k0+k1θ+k2θ2+k3θ3,

where the block derives the coefficients, ki, using these quantities:

  • DplateThrottle plate diameter parameter

  • ArestThrottle rest angle parameter

  • Apeak — 85°

  • SleakThrottle leakage area parameter

The throttle open area uses these constraints:

Sthr(0)=Seng(0)Sthr(π2)=Seng(Apeak)ddθSthr(Arest)=ddθSeng(Arest)ddθthrSthr(Apeak)=0

where Seng is the engineering formula for the throttle open area, such that

Seng(θ)=π(Dplate2)2(1cos(θ)cos(Arest))+Sleak.

This figure shows a comparison of the polynomial result and the engineering result.

Throttle open area with respect to throttle angle. The polynomial approximation has a different shape than the engineering equation, but they trend close to one another.

Ports

Input

expand all

Physical signal input associated with the normalized throttle command, where 0 represents a fully closed throttle, and 1 represents a fully open throttle. The block maps this input to a throttle-open angle between 0 and 90 degrees.

Physical signal input associated with the rotational speed, in rad/s.

Output

expand all

Physical signal input associated with the intake manifold, in kPa.

Parameters

expand all

Specific gas constant of dry air.

Specific heat ratio of dry air.

Pressure of the air entering the intake.

Temperature of the air entering the intake.

Radial distance from the resting throttle position to 0 degrees.

Diameter of the butterfly valve disc in the throttle body.

Cumulative area of throttle leakage.

Volumetric capacity of the intake manifold.

Combined displacement of all engine pistons.

Time constant associated with the response rate of the engine to the throttle command.

Version History

Introduced in R2022a