Troels Dyhr Pedersen.indd - Solid Mechanics
Troels Dyhr Pedersen.indd - Solid Mechanics
Troels Dyhr Pedersen.indd - Solid Mechanics
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
dQ<br />
dt<br />
= h ⋅ A ⋅<br />
- 32 - -<br />
( T − T )<br />
Where h is the convective heat transfer coefficient with units W/m 2 K. Woschni<br />
developed an expression for this coefficient, based on the empirical correlations for flow.<br />
With HCCI combustion this expression was in need of a revision to account for the<br />
difference with lean combustion. Junseok et al [20] made an extensive experimental study<br />
on the instantaneous heat transfer in HCCI combustion and proposed a modified<br />
expression for the heat transfer coefficient:<br />
gas<br />
sur<br />
−0.<br />
73<br />
−0.<br />
2 0.<br />
8<br />
h = α ⋅ L ⋅ p ⋅Tgas<br />
⋅ v<br />
Here, is a scaling factor which is adjusted to satisfy the total energy balance. The rest of<br />
the left hand terms are instantaneous values. L is the characteristic length, which is<br />
changed from being the cylinder bore in the original expression, to the instantaneous<br />
chamber height. p and T are the pressure and temperature, respectively. v is the gas<br />
velocity, which is divided into two separate terms, with the first being a mean gas<br />
velocity and the second a combustion induced velocity:<br />
C2<br />
VdT<br />
r<br />
v = C S p + ( p − p<br />
6 p V<br />
1 mot<br />
r r<br />
Sp is the average piston speed and Vd is the displacement volume. Tr, pr and Vr are values<br />
of temperature, pressure and volume at a reference location, such as intake valve closing.<br />
p is the pressure and pmot is the motored pressure at the same reference location. The<br />
modification to this expression is that C2 is divided by 6, since the combustion induced<br />
velocity is much lower with HCCI combustion than with SI combustion.<br />
The constants C1 and C2 are engine specific, but may be approximated by the<br />
expressions:<br />
πBw<br />
p<br />
C1<br />
= 2. 28 + 3.<br />
08 ; C2<br />
= 0.<br />
00324<br />
S<br />
Where wp is the swirl factor, which must be estimated if not provided by the engine<br />
manufacturer.<br />
p<br />
)