SpectraSim Documentation
Parameter · 05 Boundary Conditions

Outer environment

wall.thermal.environment Step 05 · Boundary Conditions

What surrounds the outside of a heat-loss wall, which decides the outer convective heat-transfer coefficient hh: natural convection in still air, forced cross flow from the vehicle's motion, or a typed hh.

Last updated

Theory

The outer surface is treated as a horizontal cylinder of diameter DoutD_{out}: the gas-side characteristic diameter 4V/Awalls4V/A_{walls} (exact for a pipe, an equivalent length for a can-shaped muffler) plus twice the layer stack. The air properties are those of dry air at the film temperature Tf=(Ts+Ta)/2T_f = (T_s + T_a)/2 and the ambient pressure: ideal-gas density, Sutherland viscosity, cp(T)c_p(T) of air, k=μcp/Prk = \mu c_p/Pr with PrPr of the Material step, β=1/Tf\beta = 1/T_f.

The solver receives this convective coefficient; radiation is added by the boundary condition itself through the emissivity. Because hh depends on the outer surface temperature TsT_s, it is evaluated first at a guessed TsT_s and updated once from the solved heat flux after the start-up phase (Thermal condition of a wall).

Equations

Nu=[0.60+0.387 Ra1/6(1+(0.559/Pr)9/16)8/27]2,Ra=g β ∣Ts−Ta∣ Dout3ν αNu = \left[0.60 + \frac{0.387\,Ra^{1/6}}{\left(1 + (0.559/Pr)^{9/16}\right)^{8/27}}\right]^2, \qquad Ra = \frac{g\,\beta\,\lvert T_s - T_a\rvert\,D_{out}^3}{\nu\,\alpha}
(1)
Nu=0.3+0.62 Re1/2Pr1/3(1+(0.4/Pr)2/3)1/4[1+(Re282000)5/8]4/5,Re=VDoutνNu = 0.3 + \frac{0.62\,Re^{1/2}Pr^{1/3}}{\left(1 + (0.4/Pr)^{2/3}\right)^{1/4}}\left[1 + \left(\frac{Re}{282000}\right)^{5/8}\right]^{4/5}, \qquad Re = \frac{V D_{out}}{\nu}
(2)
h=Nu kDout,Dout=4VfluidAwalls+2∑itih = \frac{Nu\,k}{D_{out}}, \qquad D_{out} = \frac{4V_{fluid}}{A_{walls}} + 2\sum_i t_i
(3)
wheremeaningunit
NuNuNusselt number-
RaRaRayleigh number-
ReReReynolds number of the cross flow-
PrPrPrandtl number (Material step)-
ggGravitational acceleration, 9.80665m/s²
β\betaExpansion coefficient, 1/Tf1/T_f1/K
ν,α\nu, \alphaKinematic viscosity and thermal diffusivity of air at TfT_fm²/s
kkConductivity of air at TfT_fW/(m K)
VVVehicle speedm/s
DoutD_{out}Outer characteristic diameterm
VfluidV_{fluid}Fluid volume of the modelm³
AwallsA_{walls}Area of the outer wall patchesm²
tit_iThicknesses of shell and insulationm

Default and rule

Default
Still airComputed by the program
Rule

Still air. The automatic hh needs the exported geometry (volume and wall area); until then it is pending.

Source

Correlations as cited. Representing the whole body by a single cylinder of equal volume-to-wall-area ratio is the program's simplification.

Valid range Still air | Vehicle moving | Manual h

One cylinder

The whole muffler is one horizontal cylinder of diameter 4V/A4V/A; a flat-sided or vertical body has a different natural convection.

Cross flow

Churchill-Bernstein is for flow across a cylinder (Re Pr>0.2Re\,Pr > 0.2); flow along the muffler's axis, or under a vehicle floor, is not modelled.

References 2

  1. Churchill, S. W. and Chu, H. H. S. (1975). Correlating equations for laminar and turbulent free convection from a horizontal cylinder. International Journal of Heat and Mass Transfer, 18(9), 1049-1053. DOI 10.1016/0017-9310(75)90222-7
  2. Churchill, S. W. and Bernstein, M. (1977). A correlating equation for forced convection from gases and liquids to a circular cylinder in crossflow. Journal of Heat Transfer, 99(2), 300-306. DOI 10.1115/1.3450685