A Novel Way to Control Wall Temperature Distribution by Grading Blowing Rate

Zheng, Zihao; Bai, Xiaohui; Nakayama, Akira

doi:10.1115/1.4055889

Abstract

The Graetz problem in a transpiration-cooled channel was analytically attacked so as to explore the developing temperature field due to a sudden change in wall temperature of the channel subject to an arbitrary distribution of the local mass flux over the porous wall. Analytical expressions for the developments of the thermal boundary layer thickness, wall temperature, and Nusselt number were obtained for the thermal entrance region, assuming hydrodynamically forced convective flow in a channel with a locally variable blowing mass flux. When the blowing mass flux is kept constant over the wall surface, the cooling by the coolant is less effective near the entrance, thus, exposing to danger of thermal damage. This study reveals that the blowing mass flux graded inversely proportional to one-third power of the axial distance is quite effective to keep the wall temperature uniform. Numerical calculations based on finite volume method were also carried out to verify the analysis. The findings from this study can be applied to possible thermal managements of heat generating stacks such as in EV batteries and PEMFC, in which temperature uniformity is essential for product longevity.

1 Introduction

The Graetz problem for the high Prandtl number case refers to the problem of developing temperature field due to a sudden change in wall temperature, in which the flow is already hydrodynamically fully developed, while a thin thermal boundary layer just starts to form there [1]. In fact, the Graetz model serves as a reasonable convective entrance model for all forced convective flows of high Prandtl number fluids in conduits, since the velocity field of such high Prandtl number fluid develops much faster than the temperature field. In these cases, one may assume that the flow in a channel spreads its fully developed velocity profile almost immediately after the entrance, whereas the temperature profile starts to develop right from the entrance. The Graetz problems associated with laminar and turbulent forced convection channel flows are found in many engineering processes and thus have been widely investigated theoretically and experimentally. Furthermore, recent investigations of microchannels have driven an increasing interest in thermal entrance characteristics in laminar flow regimes. A considerable number of reviews on this topic may be found in the literature, e.g., [2–5].

Due to the advent of advanced permeable wall materials, transpiration cooling is now considered to be one of the most effective ways to thermally protect hot wall surfaces. There exist many cases in which inner wall surfaces are exposed to a high temperature environment, such as in jet engines, rocket engines, and turbines. In such cases, a hot gas enters in a passage whose wall surface can be thermally damaged due to high thermal loads unless some positive cooling is used. In a transpiration cooling system, a coolant fluid from a cold reservoir passes through a permeable wall, effectively absorbing the conductive heat within the permeable wall structure, whose surface is exposed to the hot gas environment. As the coolant exits the permeable wall and then covers the hot inner wall surface, it insulates the wall from the hot gas, reducing the heat flux entering into the wall. The transpiration cooling is considered to be more thermally effective to the film cooling in which the coolant is injected through a limited number of rows of discrete holes on the wall surface.

In this study, the Graetz version of forced convective channel flows subject to locally variable blowing ratio will be treated in consideration of possible engineering applications to controlled transpiration cooling systems. It is well-known that, for the case of uniform blowing over the uniformly porous wall, the fully developed flow field is established in which the velocity profile becomes independent of axial distance, while its bulk velocity increases linearly as the coolant fluid is blown at a constant blowing ratio into the channel [6]. The momentum transfer in channel and duct flows with uniformly porous walls with fluid injection (or suction) has been theoretically attacked by a considerable number of researchers, including Berman [7], Sellars [8], Yuan [9], Donoughe [10], Eckert et al. [11], and White [12]. Zhou and Catton [13] and Zhou et al. [14] introduced a volume averaging theory based modeling and its closure for cooling problems encountered in engineering applications. A volume averaging procedure similar to theirs have been successfully applied for transpiration cooling in a channel to attack the corresponding energy balance equation by Zhang et al. [15] and Bai et al. [16].

Recently, a two-domain analytical approach has been proposed by Nakayama et al. [17] to reveal the fully developed temperature fields in forced convective flows in a channel with uniform porous walls subject to transpiration cooling. Both temperature fields in the porous wall domain and the channel domain with uniform porous walls are treated simultaneously and matched on the interstitial porous wall. However, little has been reported in the literature with thermal entry problem of forced convective flow in a transpiration cooled channel. In many cases of transpiration cooling, the wall is partially cooled by transpiration, where the thermal damage is expected unless some active cooling is taken. Thus, the Graetz version of forced convective channel flows subject to locally graded mass flux must be fully explored if one is to apply such a cooling system to control the wall surface temperature as desired.

In this study, an analytical solution for the temperature field in the thermal entrance region is obtained assuming hydrodynamically fully developed forced convective flow in a channel in which the coolant fluid is supplied through the porous wall at a locally variable blowing mass flux. When the blowing mass flux is kept constant over the wall surface, the cooling by the coolant is less effective near the entrance since the coolant film covering the wall is thinner, and thus the wall temperature is higher there. The task in this study is to find the cooling effectiveness as a function of the level of the blowing mass flux and also a local distribution of the blowing mass flux which makes the wall temperature uniform over the entire transpiration-cooled wall section. Such a local distribution of the blowing mass flux can be experimentally implemented by grading the blowing rate locally by means of variable permeability or local pressure distribution within the coolant reservoir.

Numerical calculations based on finite volume method will also be carried out to verify the analysis. The results in this study can be quite useful for effective cooling systems for heat generating stacks such in EV batteries and PEMFC. In these cases, it is required to control its average temperature to the desired working temperature as close as possible, and at the same time, to keep the temperature distribution as uniform as possible, within a certain allowable deviation from its average [18,19].

2 Physical Model for Forced Convective Flow in a Channel Subject to Blowing

Thermal entrance problem known as Graetz problem is considered for the case of transpiration-cooled channel flow. Hydrodynamically fully developed hot fluid of uniform temperature ${T_{F}}_{B} (0)$ enters at the mass flowrate per unit area $G_{H} (kg / m^{2} s)$ into a channel of height 2H, as shown in Fig. 1. Both upper and lower porous walls are subject to blowing at locally graded mass flux $G (x) ($ kg/s $m^{2})$ (though the case of constant G is illustrated in the figure for brevity). The coolant fluid of the thermal properties as the same as those of the hot fluid is supplied from the cold reservoir of constant temperature $T_{c}$ outside, and enters into the porous wall at the variable mass flux $G (x)$ ⁠. The temperature of the coolant fluid rises from $T_{c}$ to $T_{w}$ as the coolant passes vertically through the porous wall and exits from the inner porous wall surface to the main stream of hot fluid. Thus, the heat transferred from the flowing hot fluid to the porous wall surface, conducts through the porous material, heating up the coolant fluid, as it travels through the porous wall. The unknown temperature distribution over the porous surface $T_{w} (x)$ will be determined in consideration of the heat balance in both porous wall and channel domains.

Fig. 1

View large Download slide

Energy balance for forced convective flow in a channel subject to blowing

2.1 Momentum Balance.

The continuity and momentum equations for the hot fluid in a channel can be given in terms of the velocity components u and v as

ρ_{f} \frac{\partial u}{\partial x} {+ ρ}_{f} \frac{\partial v}{\partial y} = 0

(1)

ρ_{f} u \frac{\partial u}{\partial x} {+ ρ}_{f} v \frac{\partial u}{\partial y} = - \frac{\partial p}{\partial x} + μ_{f} \frac{\partial^{2} u}{\partial y^{2}}

(2)

White [6] and many others pointed out that similarity velocity profiles exist for hydrodynamically fully developed flows with constant blowing mass flux. Their arguments can be extended to the cases of variable blowing mass flux $G (x)$ ⁠.

Thus, with the help of the mass conservation equation , one may assume that the streamwise velocity distribution is given by

u (x, y^{*}) = u_{B} (x) f (η^{*}) = u_{B} (0) (1 + \int_{0}^{x} \frac{F (x)}{H} d x) f (η^{*})

(3)

where

u_{B} (x)

is the bulk mean velocity, and

η^{*} = \frac{y^{*}}{H}

(4)

is the dimensionless vertical coordinate measured from the inner porous wall surface, while

F (x) \equiv \frac{G (x)}{G_{H}} = \frac{G (x)}{ρ_{f} u_{B} (0)}

(5)

is the blowing ratio (i.e., injection ratio). The mass conservation equation may be integrated as

ρ_{f} v - G (x) = - \int_{0}^{y^{*}} ρ_{f} \frac{\partial u}{\partial x} d y^{*} = - ρ_{f} \frac{{d u}_{B} (x)}{d x} H \int_{0}^{η^{*}} f (η^{*}) d η^{*} = - G (x) \int_{0}^{η^{*}} f (η^{*}) d η^{*}

(6)

which satisfies

{v|}_{η^{*} = 1} = 0

at the channel center. Equation can be substituted into Eq. to eliminate

v

as

Re (f^{2} + (1 - \int_{0}^{η^{*}} f (η^{*}) d η^{*}) \frac{d f}{d η^{*}}) = Λ + \frac{d^{2} f}{d {η^{*}}^{2}}

(7)

where

Re (x) = \frac{G (x) H}{μ_{f}}

(8)

is the Reynolds number based on the blowing mass flux. The dimensionless pressure gradient

Λ

is defined by

Λ = - \frac{H^{2}}{μ_{f} u_{B} (x)} \frac{d p}{d x}

(9)

Equation (7) is exact only for the case of constant $Re$ ⁠, since the velocity profile function $f$ cannot be a function of $η^{*}$ alone, when $Re$ varies with $x$ ⁠. However, a fairly accurate integral solution with locally variable $Re$ is possible as long as the blowing ratio is moderate. In fact, it will be shown shortly that the streamwise velocity profile function $f (η^{*})$ is quite insensitive to the variation of $Re$ ⁠, and almost solely depends on $η^{*}$ ⁠.

Thus, an integral solution is sought for the velocity profile function

f (η^{*})

⁠, in which the pressure drops according to the distribution of the blowing ratio

F (x)

p (x) - p (0) = - Λ \frac{μ_{f} u_{B} (0)}{ε H^{2}} (x + \int_{0}^{x} \int_{0}^{x} \frac{F (x)}{H} d x d x)

(10)

where

Λ

can be determined for given

Re

⁠. The differential equation can be solved by a standard numerical integration scheme such as Runge–Kutta–Gill method (e.g., Ref. [20]), using the boundary and compatibility conditions as follows:

{\frac{d f}{d η}|}_{η^{*} = 1} = 0

(11a)

{f|}_{η^{*} = 0} = 0

(11b)

and

\int_{0}^{1} f d η^{*} = 1

(11c)

Alternatively, Yi [21] demonstrated that sufficiently accurate velocity profiles valid for any finite blowing ratio are obtainable, using the integral version of the momentum equation with the foregoing appropriate conditions. To this end, the following velocity profile function with the unknown velocity gradient at the wall,

f^{'} (0) \equiv ({d f / d η^{*})|}_{η^{*} = 0},

is introduced:

f (η^{*}) = f^{'} (0) η^{*} + \frac{12 - 5 f^{'} (0)}{2} {η^{*}}^{2} - \frac{12 - 4 f^{'} (0)}{3} {η^{*}}^{3}

(12)

which automatically satisfies the foregoing appropriate conditions given by Eqs. –. The dimensionless pressure gradient

Λ

may be expressed in terms of

f^{'} (0)

by writing the momentum equation at the channel wall

Λ = - {\frac{d^{2} f}{d η^{2}}|}_{η^{*} = 0} + Re {\frac{d f}{d η}|}_{η^{*} = 0} = (5 + Re) f^{'} (0) - 12

(13)

where the conditions – are used. Thus, the unknown parameter

f^{'} (0)

can be determined using an integral form of the momentum equation with

Λ given by

the forgoing equation

Re \int_{0}^{1} (f^{2} + (1 - \int_{0}^{η^{*}} f (η^{*}) d η^{*}) \frac{d f}{d η^{*}}) d η^{*} = - {\frac{d f}{d η}|}_{η^{*} = 0} + (5 + Re) f^{'} (0) - 12

(14a)

which under the conditions given by Eqs. – reduces to

2 R e (\int_{0}^{1} f^{2} d η) = (4 + Re) f^{'} (0) - 12

(14b)

This integral momentum equation, when the velocity profile function is substituted, yields the following quadratic equation in terms of the unknown parameter

f^{'} (0)

⁠:

2 R e {(f^{'} (0))}^{2} - (131 R e + 420) f^{'} (0) + 12 (26 R e + 105) = 0

(15)

which can readily be solved for given

Re

as

f^{'} (0; Re) = \frac{420 + 131 R e - 420 \sqrt{1 + \frac{17}{30} R e + \frac{419}{5040} {Re}^{2}}}{4 R e}

(16)

It is well-known (Ref. [6]) that the velocity profiles

f (η^{*})

for the cases of infinitely small and large F are given by

f (η^{*}) = \frac{3}{2} (1 - {(1 - η^{*})}^{2}) : F \to 0

(17a)

f (η^{*}) = \frac{π}{2} cos (\frac{π}{2} (1 - η^{*})) : F \to \infty

(17b)

As presented in Fig. 2, Eq. (16) yields $f^{'} (0; 0) =$ 3.0, reducing Eq. (12) to the exact profile, namely, Eq. (17a), for $Re \to 0$ (i.e., $F \to 0$ ⁠), whereas it gives $f' (0; \infty) =$ 2.475 for $Re \to \infty$ (i.e., $F \to \infty$ ⁠), which is very close to the exact value ${(π / 2)}^{2} =$ 2.467 based on Eq. (17b).

Fig. 2

View large Download slide

Velocity profiles for various blowing Reynolds number

As can be observed in Fig. 2, the velocity gradients $f' (0)$ at the wall for these limiting cases, namely, $f^{'} (0; 0) =$ 3.0 for $F \to 0$ and $f' (0; \infty) =$ 2.475 for $F \to \infty$ ⁠, respectively, are fairly close to each other, despite of being the two limiting cases. In fact, the streamwise velocity profiles are quite alike for all possible cases of moderate blowing ratio. Since the blowing flow ratio F (⁠ $\propto Re$ ⁠) in the case of transpiration is typically less than a few percent, the streamwise velocity profiles may be approximated by Eq. (17a) (i.e., $f^{'} (0; 0) =$ 3.0) for the all cases of laminar flow in the transpiration-cooled channels.

2.2 Energy Balance.

Focusing on the thermal boundary layer developing over the porous wall subject to blowing, the energy conservation equation for the hot fluid flowing in the channel is given by

ρ_{f} c_{p_{f}} \frac{\partial {u T}_{F}}{\partial x} {+ ρ}_{f} c_{p_{f}} \frac{\partial v T_{F}}{\partial y^{*}} = \frac{\partial}{\partial y^{*}} k_{f} \frac{\partial T_{F}}{\partial y^{*}}

(18)

where the properties of the hot fluid, the density

ρ_{f}

⁠, thermal conductivity

k_{f}

⁠, and specific heat capacity at constant pressure

{c_{p}}_{f}

of the hot fluid are assumed to be the same as those of the coolant fluid.

u

and

v

are the hot fluid velocity components in the x and y directions, while

T_{F}

is the temperature of the hot fluid. The thermal boundary conditions for the energy equation are given by

{T_{F}|}_{x = 0} = {T_{F}}_{B} (0)

(19a)

and

{T_{F}|}_{y^{*} = 0} = T_{w} (x)

(19b)

where the wall temperature

T_{w} (x)

must be determined such that it satisfies the heat balance between the heat flux on the wall and the enthalpy of the coolant fluid supplied from the reservoir of constant temperature

T_{c}

⁠. von Wolfersdorf [22] pointed out that conductive heat absorbed into the porous wall, heats up the coolant gas, increasing its temperature from

T_{c}

to

T_{w}

as it travels through the porous wall. Thus, the following heat balance holds:

G c_{p_{f}} (T_{w} - T_{c}) = k_{f} {\frac{\partial T_{F}}{\partial y^{*}}|}_{y^{*} = 0}

(20)

As done in the momentum equation, one may apply an integral procedure also to this thermal entry problem, and integrate the energy equation across the developing thermal boundary layer thickness

δ

⁠. Thus, integrating Eq. from

y^{*} = 0

to

δ

⁠, with reference to Fig. 1, along with the continuity equation and boundary conditions and , one obtains the integral energy equation as follows:

\frac{d}{d x} \int_{0}^{δ} ρ_{f} c_{p_{f}} u (T_{F} - {T_{F}}_{B} (0)) d y^{*} + G c_{p_{f}} ({{T_{F}}_{B} (0) - T}_{w}) = - k_{f} {\frac{\partial T_{F}}{\partial y^{*}}|}_{y^{*} = 0}

(21a)

which can be combined and integrated with the heat balance equation to give

\int_{0}^{δ} ρ_{f} u (T_{F} - {T_{F}}_{B} (0)) d y^{*} = (T_{c} - {T_{F}}_{B} (0)) \int_{0}^{x} Gdx

(21b)

where

u

in this Graetz problem can be approximated by the linear velocity profile since the thermal boundary layer is much thinner than the channel height

u (x, y^{*}) = u_{B} (x) f (η^{*}) = u_{B} (0) (1 + \int_{0}^{x} \frac{F (x)}{H} d x) f (η^{*}) ≅ u_{B} (0) (1 + \int_{0}^{x} \frac{F (x)}{H} d x) f^{'} (0; 0) \frac{y^{*}}{H}

(22)

Hence, by substituting Eq. into Eq. ), one obtains

3 (1 + \int_{0}^{x} \frac{F (x)}{H} d x) {(\frac{δ}{H})}^{2} \int_{0}^{1} ζ θ (ζ) d ζ = \int_{0}^{x} \frac{F (x)}{H} d x

(23)

where the approximation $f^{'} (0; Re) ≅ f^{'} (0; 0) =$ 3.0 is used since $f^{'} (0; Re)$ is insensitive to $Re$ (⁠ $\propto F$ ⁠) in a moderate range of the blowing ratio, as confirmed in Fig. 2.

The temperature profile may be assumed as

θ (ζ) \equiv \frac{T_{F} - {T_{F}}_{B} (0)}{T_{c} - {T_{F}}_{B} (0)} = θ_{w} ({(1 - ζ)}^{2} + c {(1 - ζ)}^{3} - c {(1 - ζ)}^{4})

(24)

where

ζ = \frac{y^{*}}{δ}

(25)

The temperature profile

θ (ς)

automatically satisfies

θ (0) = θ_{w}

(26a)

θ (1) = 0

(26b)

and

{\frac{d θ}{d ζ}|}_{ζ = 1} = 0

(26c)

Writing the original energy equation at the wall as

c_{p_{f}} G {\frac{\partial T_{F}}{\partial y^{*}}|}_{y^{*} = 0} = k_{f} {\frac{\partial^{2} T_{F}}{{\partial y^{*}}^{2}}|}_{y^{*} = 0}

(27)

which, when the temperature profile is applied, reduces to

c = \frac{2 ({Pe}_{δ} + 1)}{{Pe}_{δ} + 6}

(28)

where

{Pe}_{δ} = \frac{c_{p_{f}} G δ}{k_{f}}

(29)

is the Peclet numbers based on $δ$ and the blowing mass flux $G$ ⁠.

Likewise, the heat balance equation gives

θ_{w} = \frac{T_{w} - {T_{F}}_{B} (0)}{T_{c} - {T_{F}}_{B} (0)} = \frac{{Pe}_{δ}}{{Pe}_{δ} + 2 - c} = \frac{{Pe}_{δ} ({Pe}_{δ} + 6)}{{Pe}_{δ}^{2} + 6 {Pe}_{δ} + 10} : Overall cooling effectiveness

(30)

which is one of the most important indices for the performance evaluation of transpiration cooling, namely, the overall cooling effectiveness. The foregoing equation in terms of ${Pe}_{δ} \propto F$ clearly indicates that $θ_{w}$ is zero (i.e., $T_{w} = {T_{F}}_{B} (0)$ ⁠) for $F \to 0$ ⁠, and unity (i.e., $T_{w} = T_{c}$ ⁠) for $F \to \infty$ ⁠.

Thus, the integral

\int_{0}^{1} ζ θ (ζ) d ζ

can be evaluated as

\int_{0}^{1} ζ θ (ζ) d ζ = \int_{0}^{1} ζ θ_{w} ({(1 - ζ)}^{2} + c {(1 - ζ)}^{3} - c {(1 - ζ)}^{4}) d ζ = \frac{θ_{w}}{12} (1 + \frac{c}{5}) = \frac{{Pe}_{δ} (7 {Pe}_{δ} + 32)}{60 ({Pe}_{δ}^{2} + 6 {Pe}_{δ} + 10)}

(31)

which is substituted into Eq. to find

\frac{{Pe}_{δ}^{3} (7 {Pe}_{δ} + 32)}{{Pe}_{δ}^{2} + 6 {Pe}_{δ} + 10} = \frac{20 {Pe}^{2} \int_{0}^{ξ} Pe d ξ}{1 + \int_{0}^{ξ} Pe d ξ}

(32)

where

Pe (x) = RePr = \frac{c_{p_{f}} G (x) H}{k_{f}} = \frac{{c_{p_{f}} G}_{H} H}{k_{f}} F (x)

(33)

is the blowing Peclet number based on the blowing mass flux, so that

{Pe}_{δ} = Pe (δ / H) = {(c_{p_{f}} G}_{H} H / k_{f}) (δ / H) F

⁠, and

ξ = \frac{(\frac{x}{H})}{(\frac{{c_{p_{f}} G}_{H} H}{k_{f}})} \propto x

(34)

is the Graetz number (note that some define its inverse as Graetz number). According to Bai and Nakayama [23], the linear (Leveque) approximation of the velocity profile provides reasonably accurate results for the entrance region given by $0 < ξ < 0.1,$ namely, $(x / H) < 0.1 ({c_{p_{f}} G}_{H} H / k_{f})$ ⁠. This covers a sufficiently large entrance region especially when the channel Peclet number $({c_{p_{f}} G}_{H} H / k_{f})$ is high.

For given

ξ and Pe (x) (

i.e., graded distribution of the blowing mass flux

G (x)

or

F (x)

⁠), the value of the right-hand side in Eq. can readily be calculated. Then, the local value

{Pe}_{δ} (x)

on the left-hand side in Eq. can be determined iteratively. Once

{Pe}_{δ} (x)

is found, the overall cooling effectiveness

θ_{w} (x)

is given by Eq. , while Nusselt number

{Nu}_{H}

of interest can be evaluated from

{Nu}_{H} = \frac{{- k}_{f} {\frac{\partial T_{F}}{\partial y^{*}}|}_{y^{*} = 0} H}{k_{f} (T_{w} - {T_{F}}_{B} (0))} = \frac{H}{δ} (2 - c) = \frac{10 P e}{{Pe}_{δ} ({Pe}_{δ} + 6)}

(35)

3 Results and Discussion

3.1 Uniform Blowing.

When the blowing mass flux

G (x)

is constant over the porous wall (i.e.,

Pe (x) = c_{p_{f}} G (x) H / k_{f} = const .

⁠) Eq. reduces to an algebraic equation

\frac{{Pe}_{δ}^{3} (7 {Pe}_{δ} + 32)}{{Pe}_{δ}^{2} + 6 {Pe}_{δ} + 10} = \frac{20 {Pe}^{3} ξ}{1 + Pe ξ}

(36)

For the limiting case of no blowing,

G \propto F \to 0

⁠, both

{Pe}_{δ}

and

Pe (> {Pe}_{δ})

become much less than unity. Thus, Eq. reduces to

ξ ≅ \frac{4}{25} {(\frac{{Pe}_{δ}}{Pe})}^{3}

(37)

Substituting this into Eq. , one finds

{Nu}_{H} = \frac{10 P e}{{Pe}_{δ} ({Pe}_{δ} + 6)} ≅ \frac{5 P e}{3 {Pe}_{δ}} = \frac{5}{3} \times {(\frac{4}{25 ξ})}^{\frac{1}{3}} = \frac{0.905}{ξ^{1 / 3}} : G \propto F \to 0

(38)

which is fairly close to the exact solution for the Graetz channel flow without blowing subject to constant heat flux [24], namely, ${Nu}_{H} = 0.939 / ξ^{1 / 3}$ ⁠, indicating the soundness of the present approximate analysis. Equation (38) based on the Leveque approximation holds over the entrance region from $0 < ξ < 0.1$ ⁠.

For the general cases of the dimensionless mass flux, $Pe$ = 0.1, 0.5, and 1.0, the dimensionless thermal boundary layer thickness ${Pe}_{δ} (ξ)$ ⁠, the overall cooling effectiveness $θ_{w} (ξ)$ ⁠, and the local Nusselt number ${Nu}_{H} (ξ)$ are plotted against the Graetz coordinate $ξ$ in Figs. 3(a)–3(c), respectively, using the corresponding algebraic equations (36a), (30), and (35). As Fig. 3(a) shows, the thermal boundary layer thickness $δ \propto {Pe}_{δ} = c_{p_{f}} G δ / k_{f}$ grows downstream in a manner similar to the case of no blowing, namely, $δ \propto ξ^{1 / 3}$ ⁠. Naturally, ${Pe}_{δ} = Pe (δ / H)$ is higher for higher $Pe$ ⁠.

Fig. 3

View large Download slide

Development of thermal boundary layer for the case of constant blowing mass flux: (a) dimensionless thermal boundary layer thickness, (b) overall cooling effectiveness, and (c)local Nusselt number

As illustrated in Fig. 3(b), the overall cooling effectiveness $θ_{w} (ξ)$ can be translated into the dimensionless wall temperature by $T_{w} (ξ) - T_{c} / {{T_{F}}_{B} (0) - T}_{c} = 1 - θ_{w} (ξ)$ ⁠. $θ_{w} (ξ)$ approaches unity when $T_{w} (ξ)$ reaches the coolant temperature $T_{c}$ ⁠, while it vanishes when there is no cooling effectiveness with $T_{w} (ξ)$ being the same as the hot fluid temperature ${T_{F}}_{B} (0)$ ⁠. Naturally, when the blowing mass flux $G \propto Pe$ is higher, the level of the overall cooling effectiveness $θ_{w} (ξ)$ is higher. These figures clearly indicate that the thermal resistance increases downstream as the thermal boundary layer thickens, enhancing the overall cooling effectiveness $θ_{w} (ξ)$ downstream. As a result, the wall temperature $T_{w} (ξ)$ is higher near the entrance and lower downstream. Thus, possible thermal damage is expected at the thermal entrance region around $ξ = 0$ ⁠.

The local Nusselt number variation is presented in Fig. 3(c), which shows that ${Nu}_{H}$ depends almost solely on the Graetz coordinate $ξ$ and is very much insensitive to the blowing mass flux $G \propto Pe$ ⁠. Since the Nusselt number is high at the entrance, the convective heat transfer rate from the hot fluid to the porous wall is high, heating the wall surface up. On the contrary, it becomes low downstream as the coolant fluid film covers over the wall thick, protecting the wall surface thermally.

3.2 Nonuniform Blowing.

For thermal management of heat generating stacks such in EV batteries and PEMFC, it is essential to control not only its average temperature (i.e., the average $θ_{w}$ over the wall) to a desired working temperature as close as possible but also the temperature distribution as uniform as possible, within a certain allowable deviation from its average. Thus, the task in this study is to find the relationship between the overall cooling effectiveness and the intensity of the blowing mass flux, and also its corresponding blowing mass flux distribution over the wall, which eventually accomplish these requirements for proper thermal managements.

To this end, one may revisit the energy balance relationship as given by Eq. , which may be translated in terms of the local thermal boundary layer thickness

δ (x)

as follows:

\frac{({T_{F}}_{B} (0) - T_{w} (x))}{\frac{δ (x)}{k_{f}}} ≅ \frac{(T_{w} (x) - T_{c})}{\frac{1}{G (x) c_{p_{f}}}}

(39)

This relationship is illustrated in Fig. 4, in which the highest potential of the hot fluid

{T_{F}}_{B} (0)

and the lowest potential of the coolant

T_{c}

are jointed at the potential

T_{w} (x)

at the wall by the thermal resistances

δ (x) / k_{f}

and

1 / G (x) c_{p_{f}}

in series. Obviously, the ratio of these two resistances must be kept constant in order to make

T_{w} (x)

unchanged over the wall surface such that

{Pe}_{δ} (x) = \frac{c_{p_{f}} G (x) δ (x)}{k_{f}} = const

(40)

Fig. 4

View large Download slide

Thermal circuit for blowing

This preliminary consideration prompts one to grade the blowing mass flux to follow a 1/3 power function of the distance as follows:

G (ξ) = G_{0} / ξ^{1 / 3}

which provides a higher blowing mass flux near the entrance, thickening the coolant film to reduce the wall temperature more there. For this mass flux distribution, the blowing Peclet number varies over the wall as

Pe (ξ) = \frac{c_{p_{f}} G (ξ) H}{k_{f}} = (\frac{c_{p_{f}} G_{0} H}{k_{f}}) / ξ^{1 / 3}

(41)

which reduces Eq. to

\frac{{Pe}_{δ}^{3} (7 {Pe}_{δ} + 32)}{{Pe}_{δ}^{2} + 6 {Pe}_{δ} + 10} = \frac{30 {(\frac{c_{p_{f}} G_{0} H}{k_{f}})}^{3}}{1 + \frac{3}{2} (\frac{c_{p_{f}} G_{0} H}{k_{f}}) ξ^{2 / 3}}

(42)

Note the denominator in the right-hand side term,

u_{B} (ξ) / u_{B} (0) = 1 + (3 / 2) (c_{p_{f}} G_{0} H / k_{f}) ξ^{2 / 3}

⁠, stays almost unity, since the blowing ratio is usually less than a few percent. Thus, as expected,

{Pe}_{δ} (ξ) = c_{p_{f}} G (ξ) δ (ξ) / k_{f}

in the left-hand side term remains nearly constant over the wall

(i . e ., δ (x) \propto ξ^{\frac{1}{3}})

such that, for the first approximation, one obtains

{Pe}_{δ} (ξ) = \frac{c_{p_{f}} G (ξ) δ (ξ)}{k_{f}} ≅ {Pe}_{δ} (0)

(43)

θ_{w} (ξ) ≅ \frac{{Pe}_{δ} (0) ({Pe}_{δ} (0) + 6)}{{{Pe}_{δ} (0)}^{2} + 6 {Pe}_{δ} (0) + 10}

(44)

from Eq. and

{Nu}_{H} (ξ) = \frac{10 (\frac{c_{p_{f}} G_{0} H}{k_{f}})}{{Pe}_{δ} (0) ({Pe}_{δ} (0) + 6) ξ^{1 / 3}}

(45)

from Eq. where

\frac{{{Pe}_{δ} (0)}^{3} (7 {Pe}_{δ} (0) + 32)}{{{Pe}_{δ} (0)}^{2} + 6 {Pe}_{δ} (0) + 10} = 30 {(\frac{c_{p_{f}} G_{0} H}{k_{f}})}^{3}

(46)

Once ${Pe}_{δ} (0)$ is determined from the implicit algebraic equation (46) for the given intensity of the blowing mass flux $(c_{p_{f}} G_{0} H / k_{f})$ ⁠, one can readily evaluate the dimensionless thermal boundary layer thickness $δ (ξ) / H = {Pe}_{δ} (0) / (c_{p_{f}} G_{0} H / k_{f}) ξ^{1 / 3} .$ Then, the overall cooling effectiveness $θ_{w} (ξ) ≅ const .$ is given by Eq. (44), while the local Nusselt number ${Nu}_{H} (ξ) \propto {1 / ξ}^{\frac{1}{3}}$ is determined from Eq. (45).

Alternatively, when the desired level of the overall cooling effectiveness

θ_{w}

(i.e., the wall temperature

T_{w} = T_{c} + ({{T_{F}}_{B} (0) - T}_{c}) (1 - θ_{w})

⁠) is specified,

{Pe}_{δ} (0)

is readily found by solving the quadratic equation as

{Pe}_{δ} (0) = 3 (\sqrt{\frac{1 + \frac{θ_{w}}{9}}{1 - θ_{w}}} - 1)

(47)

Then, the required level of the blowing mass flux $(c_{p_{f}} G_{0} H / k_{f})$ can be evaluated from Eq. (46). With ${Pe}_{δ} (0)$ and $(c_{p_{f}} G_{0} H / k_{f})$ thus evaluated, the local Nusselt number ${Nu}_{H} (ξ)$ can finally be estimated from Eq. (45).

In order to confirm the validity of the foregoing arguments based on $G (ξ) = G_{0} / ξ^{1 / 3}$ ⁠, Eq. (42) is iteratively solved without approximations (with retaining the denominator in the right-hand side term as it is) over the range of $ξ$ = 0.005–0.10 for the cases of the blowing mass flux intensity $(c_{p_{f}} G_{0} H / k_{f}) =$ 0.1, 0.5, and 1. The resulting streamwise variations of ${Pe}_{δ} (ξ)$ ⁠, $θ_{w} (ξ)$ ⁠, and ${Nu}_{H} (ξ)$ are presented in Figs. 5(a)–5(c), respectively. In these figures, the approximate results based on ${Pe}_{δ} (0)$ are indicated by dashed lines, while those based on the local value ${Pe}_{δ} (ξ)$ iteratively obtained are indicated by solid lines. As expected, ${Pe}_{δ} (ξ) \propto G (ξ) δ (ξ)$ in Fig. 5(a) remains almost constant over the wall, despite of slight decrease in ${Pe}_{δ} (ξ)$ downstream. As a result, the wall temperature (i.e., ${1 - θ}_{w} (ξ))$ is kept nearly constant, as can be confirmed in Fig. 5(b), while ${Nu}_{H} (ξ)$ decreases as increasing the blowing mass flux intensity, as indicated in Fig. 5(c), which is consistent with the well-known fact reported by Kays and Crawford [25] for transpiration cooling.

Fig. 5

View large Download slide

Development of thermal boundary layer for the case of nonuniform: (a) local blowing Peclet number based thermal boundary layer thickness, (b) overall cooling effectiveness, and (c) local Nusselt number

As demonstrated in Fig. 5(b), the blowing mass flux graded according to $G (ξ) = G_{0} / ξ^{1 / 3}$ is found quite satisfactory to keep the wall temperature uniform. Furthermore, the figure indicates that one can control the wall temperature freely to a desired level by adjusting the blowing mass flux intensity $(c_{p_{f}} G_{0} H / k_{f})$ ⁠. Naturally, the higher $(c_{p_{f}} G_{0} H / k_{f})$ brings out the lower wall temperature, as illustrated in Fig. 6, in which $(c_{p_{f}} G_{0} H / k_{f})$ is plotted against the wall temperature (i.e., ${1 - θ}_{w} (ξ)),$ using Eqs. (47) and (46). The figure shows that the blowing mass flux intensity $(c_{p_{f}} G_{0} H / k_{f})$ increases steeply from unity, as one tries to reduce the wall temperature to $T_{w} (ξ) - T_{c} / {{T_{F}}_{B} (0) - T}_{c} < 0.3$ ⁠. Thus, in engineering applications, $(c_{p_{f}} G_{0} H / k_{f})$ should be set below $(c_{p_{f}} G_{0} H / k_{f}) < 1$ ⁠, since 70% (or less) overall thermal effectiveness (i.e., $(T_{w} (ξ) - T_{c})$ > 0.3 $({{T_{F}}_{B} (0) - T}_{c})$ ⁠) is usually quite sufficient for most cases of thermal managements. The figure also indicates that, in this range, the sensitivity $- d θ_{w} / d (c_{p_{f}} G_{0} H / k_{f})$ is high so that the wall temperature can easily be set to a desired level by adjusting the blowing mass flux intensity $(c_{p_{f}} G_{0} H / k_{f})$ ⁠.

Fig. 6

View large Download slide

Blowing mass flux intensity against dimensionless wall temperature

3.3 Comparison of Numerical and Analytical Results.

Finally, numerical simulations based on the finite volume method [20] were carried out for a limited number of cases to verify the validity of the present analysis. The in-house computer code SAINTS [20] based on staggered grids was used for numerical calculations. Appropriate velocity boundary treatments such as no-slip streamwise velocity component and variable vertical velocity component are applied on the wall boundaries, while the thermal boundary condition was set according to the heat balance relationship on the permeable wall, as given by Eq. (20). In actual numerical computations, the variables normalized by the references, such as $ρ_{f,} u_{B} (0)$ ⁠, H and $({T_{F}}_{B} (0) - T_{c}),$ were used so that the convergence was measured in terms of the maximum change in each normalized variable during an iteration. The maximum change allowed for the convergence check was set to an arbitrarily small value such as 10⁻⁵.

The wall temperature variations $(1 - θ_{w})$ for the case of $G_{H} H / μ_{f} = 100$ and $c_{p_{f}} μ_{f} / k_{f} = 1$ (i.e., ${c_{p_{f}} G}_{H} H / k_{f} = 100$ ⁠) were obtained for the variable mass flux distribution $G (ξ) = G_{0} / ξ^{1 / 3}$ with $(c_{p_{f}} G_{0} H / k_{f})$ = 0.1, 0.5, and 1.0, and compared with those based on the present analysis in Fig. 7. The channel length is set to 10H such that the whole channel stays within the entrance region, namely, $0 \leq ξ \leq 0.1 .$ As will be shown shortly, the blowing parameter $(c_{p_{f}} G_{0} H / k_{f})$ = 0.1, 0.5, and 1.0 correspond with the cases of uniform blowing ratio of F = 0.32%, 1.6%, and 3.2%, respectively. As predicted in the analysis, the wall temperature in fact stays nearly constant over the channel wall. The levels predicted by both numerical and analytical methods agree well with each other. The temperature rise toward the entrance indicated by the numerical method is due to the fact that the limiting condition ${G_{0} / ξ^{1 / 3}|}_{ξ \to 0} \to \infty$ can never faithfully be described in the finite volume method.

Fig. 7

View large Download slide

Axial variation of dimensionless wall temperature for ${c_{p_{f}} G}_{H} H / k_{f} = 100$

In Fig. 8, the temperature field numerically obtained for the case of graded mass flux distribution $(c_{p_{f}} G_{0} H / k_{f})$ = 1.0 is presented along with the case of constant blowing ratio $F$ in which the mass flowrate at $ξ = 0.1$ coincides with that of graded mass flux, such that

Fig. 8

View large Download slide

Velocity and temperature fields for constant and variable blowing mass fluxes: (a) constant mass flux ${G / G}_{H} = 0.0323$ and (b) variable mass flux ${G / G}_{H} = (1 / 100) (1 / ξ^{1 / 3})$

F = \frac{1}{0.1} \int_{0}^{0.1} \frac{G}{G_{H}} d ξ = \frac{1}{0.1} \int_{0}^{0.1} \frac{(\frac{c_{p_{f}} G_{0} H}{k_{f}})}{(\frac{c_{p_{f}} G_{H} H}{k_{f}})} \frac{1}{ξ^{1 / 3}} d ξ = \frac{1}{0.1} \int_{0}^{0.1} \frac{1}{100} \frac{1}{ξ^{1 / 3}} d ξ = 0.032

The velocity and temperature fields for the case of constant blowing mass flux and that of graded blowing mass flux are presented in Figs. 8(a) and 8(b), respectively. The higher blowing mass flux near the entrance for the case of variable blowing mass flux suppresses the wall temperature rise there, resulting the wall temperature much lower there than the corresponding temperature for the case of constant blowing mass flux. On the other hand, the wall temperature downstream for the case of variable blowing mass flux in Fig. 8(b) stays comparatively higher than that for the case of constant blowing mass flux in Fig. 8(a). It can be confirmed from the temperature field presented in Fig. 8(b) that the temperature can be kept nearly constant over the wall by controlling the blowing mass flux locally.

4 Conclusions

Thermal entry problems associated with transpiration cooling were analytically treated assuming hydrodynamically fully developed forced convective flow in a channel in which the coolant fluid is supplied through the porous wall at a locally variable blowing mass flux. First, the hydrodynamically fully developed velocity fields in a channel were investigated allowing an arbitrary distribution of local blowing mass flux. The velocity profiles were found fairly insensitive to the blowing mass flux for a practical range of the blowing intensity. Second, the energy balance equation in the entrance region of the channel was considered along with the heat balance equation between the enthalpy increase of the coolant and the convective heat flux over the porous wall, in order to elucidate the effects of the blowing mass flux distribution on the wall temperature distribution. The analysis clearly shows that the uniform blowing results in a thin coolant film near the inlet and a thick one downstream, yielding a high wall temperature at the inlet, whereas the variable blowing rate inversely proportional to one-third power of the distance remedies the situation keeping the wall temperature uniform everywhere.

Conflict of Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Nomenclature

$c$ =: shape parameter for the temperature profile
${c_{p}}_{f}$ =: specific heat of fluid at constant pressure (J/kg K)
$f$ =: $u / u_{B}$ velocity profile function
$F$ =: $G / G_{H}$ blowing ratio
$f^{'} (0)$ =: ${d f / d η^{*}|}_{η^{*} = 0}$ velocity gradient at the wall
$G$ =: blowing mass flux (kg/m²s)
$G_{H}$ =: hot fluid mass flow rate per unit area at the inlet (kg/m²s)
$G_{0}$ =: blowing mass flux intensity $G (ξ) = G_{0} / ξ^{1 / 3}$
H =: half of the channel height
$k_{f}$ =: thermal conductivity of the fluid (W/m K)
${Nu}_{H}$ =: ${- k}_{f} {(\partial T_{F} / \partial y^{*})|}_{y^{*} = 0} H / k_{f} (T_{w} - {T_{F}}_{B} (0))$ Nusselt number
$p$ =: pressure (Pa)
$Pe$ =: ${c_{p}}_{f} G H / {k_{F}}_{e}$ Peclet number based on $G$
${Pe}_{δ}$ =: $c_{p_{f}} G δ / k_{f}$ Peclet number based on $G$ and $δ$
$\Pr$ =: ${{c_{p}}_{f} μ}_{f} / k_{f}$ Prandtl number
$Re$ =: $G H / μ_{f}$ blowing Reynolds number
$T_{c}$ =: coolant reservoir temperature (K)
$T_{F}$ =: hot fluid temperature (K)
${T_{F}}_{B}$ =: hot fluid bulk mean temperature (K)
$T_{w}$ =: channel wall temperature (K)
$u$ ⁠, v =: x and y velocity components (m/s)
$u_{B}$ =: bulk mean velocity of hot fluid (m/s)
$x$ =: streamwise coordinate (m)
$y$ =: vertical coordinate (m) from the channel center
$y^{*}$ =: vertical coordinate (m) from the inner wall

Greek Symbols

$δ$ =: thermal boundary layer thickness $ζ = y^{*} / δ$
$η^{*}$ =: $y^{*} / H$ dimensionless vertical coordinate
$θ$ =: $θ = T_{F} - {T_{F}}_{B} (0) / T_{c} - {T_{F}}_{B} (0)$ dimensionless temperature profile of the hot fluid
$θ_{w}$ =: $T_{w} - {T_{F}}_{B} (0) / T_{c} - {T_{F}}_{B} (0)$ overall cooling effectiveness
$Λ$ =: $Λ = (- H^{2} / μ_{f} u_{B} (x)) d p / d x$ dimensionless pressure gradient
$μ_{f}$ =: fluid viscosity (Pa·s)
$ξ$ =: $(x / H) / ({c_{p_{f}} G}_{H} H / k_{f})$ Graetz number
$ρ_{f}$ =: fluid density (kg/m³)

Subscripts

B =: bulk mean
c =: cold reservoir side
f =: fluid
F =: hot fluid
$w$ =: wall

References

1.

Graetz

,

L.

,

1885

, “

Ueber Die Wärmeleitungsfähigkeit Von Flüssigkeiten

,”

Annu. Rev. Phys. Chem.

,

25

, pp.

337

–

357

.10.1002/andp.18822540106

2.

Shah

,

R. K.

, and

London

,

A. L.

,

1978

,

Laminar Flow Forced Convection in Ducts: A Source Book for Compact Heat Exchanger Analytical Data

,

Academic Press

,

New York

, pp.

78

–

138

.

Google Scholar

Crossref

3.

White

,

F. M.

,

1974

,

Viscous Fluid Flow

,

McGraw-Hill

,

New York

, pp.

132

–

138

.

4.

Bejan

,

A.

,

2013

,

Convective Heat Transfer

, 4th ed.,

Wiley

, Hoboken, NJ, pp.

121

–

129

.

Google Scholar

Crossref

5.

Lienhard

,

J. H.

, IV , and

Lienhard

,

V. J. H.

,

2011

,

A Heat Transfer Textbook

, 3rd ed.,

Dover Publications

, Cambridge, MA, pp.

327

–

330

.

6.

White

,

F. M.

,

1974

,

Viscous Fluid Flow

,

McGraw-Hill

,

New York

, pp.

153

–

159

.

7.

Berman

,

R. S.

,

1953

, “

Laminar Flow in Channels With Porous Walls

,”

J. Appl. Phys.

,

24

(

9

), pp.

1232

–

1235

.10.1063/1.1721476

Google Scholar

Crossref

8.

Sellars

,

J. R.

,

1955

, “

Laminar Flow in Channels With Porous Walls at High Suction Reynolds Numbers

,”

J. Appl. Phys.

,

26

(

4

), pp.

489

–

490

.10.1063/1.1722024

Google Scholar

Crossref

9.

Yuan

,

S. W.

,

1956

, “

Further Investigation of Laminar Flow in Channels With Porous Walls

,”

J. Appl. Phys.

,

27

(

3

), pp.

267

–

269

.10.1063/1.1722355

Google Scholar

Crossref

10.

Donoughe

,

P. L.

,

1954

, “

Exact Solutions of Laminar-Boundary-Layer Equations With Constant Property Values for Porous Wall With Variable Temperature

,” NACA, Washington, DC, Report No.

NACA-TN

-

3579

.

11.

Eckert

,

E. R. G.

,

Donoughe

,

P. L.

, and

Moore

,

B. J.

,

1957

, “Velocity and Friction Characteristics of Laminar Viscous Boundary-Layer and Channel Flow Over Surfaces With Ejection or Suction,” NACA, Washington, DC, Report No. NACA-TN-4102.

12.

White

,

F. M.

,

1962

, “

Laminar Flow in a Uniformly Porous Tube

,”

ASME J. Appl. Mech.

,

29

(

1

), pp.

201

–

204

.10.1115/1.3636459

Google Scholar

Crossref

13.

Zhou

,

F.

, and

Catton

,

I.

,

2012

, “

Volume Averaging Theory (VAT) Based Modeling and Closure Evaluation for Fin-and-Tube Heat Exchangers

,”

Heat Mass Transfer

,

48

(

10

), pp.

1813

–

1823

. 10.1007/s00231-012-1025-7

Google Scholar

Crossref

14.

Zhou

,

F.

,

Hansen

,

N. E.

,

Geb

,

D. J.

, and

Catton

,

I.

,

2011

, “

Obtaining Closure for Fin-and-Tube Heat Exchanger Modeling Based on Volume Averaging Theory (VAT)

,”

ASME J. Heat Transfer

,

133

(

11

), p.

111802

.10.1115/1.4004393

Google Scholar

Crossref

15.

Zhang

,

W.

,

Yi

,

Y.

,

Bai

,

X.

, and

Nakayama

,

A.

,

2020

, “

A Local Thermal Non-Equilibrium Analysis for Convective and Radiative Heat Transfer in Gaseous Transpiration Cooling Through a Porous Wall

,”

Int. J. Heat Mass Transfer

,

162

, p.

120389

.10.1016/j.ijheatmasstransfer.2020.120389

Google Scholar

Crossref

16.

Bai

,

X.

,

Yi

,

Y.

,

Liu

,

C.

,

Zhang

,

W.

, and

Nakayama

,

A.

,

2021

, “

A Simple Conjugate Analysis and Its Comparison With Experiment for Heat Transfer Problems Associated With Hot Gas Flows in a Partially Transpiration-Cooled Channel

,”

Int. J. Heat Mass Transfer

,

165

, p.

120729

.10.1016/j.ijheatmasstransfer.2020.120729

Google Scholar

Crossref

17.

Nakayama

,

A.

,

Yi

,

Y.

,

Bai

,

X.

, and

Zhang

,

W.

,

2022

, “

A Two-Domain Analytical Approach for Laminar and Turbulent Forced Convection in a Transpiration-Cooled Channel

,”

Int. J. Heat Mass Transfer

,

184

, p.

122247

.10.1016/j.ijheatmasstransfer.2021.122247

Google Scholar

Crossref

18.

Afzal

,

A.

,

Mohammed Samee

,

A. D.

,

Abdul Razak

,

R. K.

, and

Ramis

,

M. K.

,

2021

, “

Thermal Management of Modern Electric Vehicle Battery Systems (MEVBS)

,”

J. Therm. Anal. Calorim.

,

144

(

4

), pp.

1271

–

1285

.10.1007/s10973-020-09606-x

Google Scholar

Crossref

19.

Kandlikar

,

S. G.

, and

Lu

,

Z.

,

2009

, “

Thermal Management Issues in a PEMFC Stack—A Brief Review of Current Status

,”

Appl. Therm. Eng.

,

29

(

7

), pp.

1276

–

1280

.10.1016/j.applthermaleng.2008.05.009

Google Scholar

Crossref

20.

Nakayama

,

A.

,

1995

,

PC-Aided Numerical Heat Transfer and Convective Flow

,

CRC Press

,

Boca Raton, FL

.

21.

Yi

,

Y.

,

2021

, “General Local Thermal Non-Equilibrium Treatments for External and Internal Flows in Porous Media,” Ph.D. dissertation, Shizuoka University, Shizuoka, Japan.

22.

von Wolfersdorf

,

J.

,

2005

, “

Effect of Coolant Side Heat Transfer on Transpiration Cooling

,”

Heat Mass Transfer

,

41

(

4

), pp.

327

–

337

.10.1007/s00231-004-0549-x

Google Scholar

Crossref

23.

Bai

,

X.

, and

Nakayama

,

A.

,

2019

, “

A Local Thermal Non-Equilibrium Integral Analysis for Forced Convective Thermal Boundary Development in a Channel Filled With in a Fluid Saturated Porous Medium

,”

Int. J. Heat Mass Transfer

,

142

, p.

118446

.10.1016/j.ijheatmasstransfer.2019.118446

Google Scholar

Crossref

24.

Nakayama

,

A.

,

Sano

,

Y.

,

Bai

,

X.

, and

Tado

,

K.

,

2017

, “

A Boundary Layer Analysis for Determination of the Limiting Current Density in an Electrodialysis Desalination

,”

Desalination

,

404

(

17

), pp.

41

–

49

.10.1016/j.desal.2016.10.013

25.

Kays

,

W. M.

, and

Crawford

,

M. E.

,

1993

,

Convective Heat and Mass Transfer

, 3rd ed.,

McGraw-Hill

,

New York

.

A Novel Way to Control Wall Temperature Distribution by Grading Blowing Rate

Abstract

1 Introduction

2 Physical Model for Forced Convective Flow in a Channel Subject to Blowing

2.1 Momentum Balance.

2.2 Energy Balance.

3 Results and Discussion

3.1 Uniform Blowing.

3.2 Nonuniform Blowing.

3.3 Comparison of Numerical and Analytical Results.

4 Conclusions

Conflict of Interest

Nomenclature

Greek Symbols

Subscripts

References

Get Email Alerts

Cited By

Journals

Conference Proceedings

eBooks

Resources

Opportunities

A Novel Way to Control Wall Temperature Distribution by Grading Blowing Rate

Abstract

1 Introduction

2 Physical Model for Forced Convective Flow in a Channel Subject to Blowing

2.1 Momentum Balance.

2.2 Energy Balance.

3 Results and Discussion

3.1 Uniform Blowing.

3.2 Nonuniform Blowing.

3.3 Comparison of Numerical and Analytical Results.

4 Conclusions

Conflict of Interest

Nomenclature

Greek Symbols

Subscripts

References

Get Email Alerts

Cited By

Related Articles

Related Proceedings Papers

Related Chapters

Journals

Conference Proceedings

eBooks

Resources

Opportunities

This Feature Is Available To Subscribers Only