Publication Cover
Numerical Heat Transfer, Part A: Applications
An International Journal of Computation and Methodology
Volume 71, 2017 - Issue 7
174
Views
8
CrossRef citations to date
0
Altmetric
Original Articles

Numerical modeling of turbulent flow structure and heat transfer in a droplet-laden swirling flow in a pipe with a sudden expansion

&
Pages 721-736 | Received 27 Oct 2016, Accepted 16 Feb 2017, Published online: 27 Apr 2017
 

ABSTRACT

This contribution presents the mathematical model to simulate the swirling turbulent gas-droplet flow in a sudden pipe expansion. The set of axisymmetrical steady-state Reynolds averaged Navier–Stokes equations (RANS) for the two-phase flow is utilized. The dispersed phase is modeled by the Eulerian approach. The flow swirl causes an increase in the intensity of heat transfer (more than 1.5 times compared with the nonswirling mist flow at other identical inlet conditions). Evaporation of the droplets leads to a significant increase in the heat transfer intensity in the swirling two-phase flow (more than 2.5 times compared with the single-phase flow).

Nomenclature

Bi=

Biot number

b1D=

diffusion parameter of injection, determined with the use of the saturation curve

d=

droplet diameter, m

CD=

drag coefficient of the evaporating droplet, written by taking into account the deviation of the Stokes law

CP=

specific heat capacity, J kg−1 K−1

D=

diffusion coefficient, m2 s−1

G=

mass flow rate, kg s−1

J=

mass flux of steam from the surface of the evaporating droplet, kg m−2 s−1

ER=

heat transfer enhancement ratio

H=

step height, m

=

steam mass concentration at the “vapor-gas mixture–droplet” interface, corresponding to the saturation parameters at droplet temperature TL

KV0=

mass concentration of steam far from the droplet

k=

k = ⟨uiui′⟩/2 turbulent kinetic energy, m2 s−2

L=

latent heat of vaporization, J kg−1

ML=

droplets mass fraction

Nu=

Nusselt number

P=

pressure (Pa)

Pr=

Prandtl number

Pij=

stress production term, m s−2

2R1=

pipe diameter before sudden expansion, m

2R2=

pipe diameter after sudden expansion, m

Re=

Re = Um12R1/ν gas-phase Reynolds number

ReL=

ReL = |U − UL|d/ν Reynolds number of the dispersed phase

r=

radial coordinate, m

S=

swirl number

Sc=

Schmidt number

Sh=

Sherwood number

StD=

diffusional Stanton number

Stk=

Stk = τ/τf Stokes number in the mean motion

T=

temperature (K)

Ui=

components of mean velocity of the gas phase, m s−1

ULi=

components of droplets mean velocity, m s−1

u′⟩, ⟨v′⟩, ⟨w′⟩=

components of velocity pulsation, m s−1

uv′⟩=

turbulent Reynolds stresses, m2 s−2

uLiuLj′⟩=

kinetic stresses in the dispersed phase, m2 s−2

x=

axial coordinate, m

xR=

separation length, m

y=

distance normal to the wall, m

We=

Weber number

We*=

critical Weber number

Greek letters=
Φ=

Φ = MLρ/ρL volumetric concentration of the dispersed phase

ΓE=

turbulence macroscale, m

ΛT=

the turbulent macroscale time, s

ΩtL=

time of interaction with temperature pulsations of the carrying flow, s

ΩLag=

Lagrangian integral timescales, s

ΩE=

Eulerian integral timescales, s

Ωϵ=

time microscale of turbulence, s

αL=

heat transfer coefficient for the evaporating droplet, W m−2

β=

blending coefficient

ϵ=

dissipation of the turbulent kinetic energy, m2s−3

ϕij=

the velocity–pressure–gradient correlation (the pressure term), m s−2

μ=

dynamic viscosity, kg m−1 s−1

ν=

kinematic viscosity, m2 s−1

ρ=

density, kg m−3

σ=

surface tension, N m−1

τ=

dynamic relaxation time of droplets, s

τf=

turbulent timescale, s

Subscripts=
0=

parameter under conditions on the pipe axis

1=

parameter under inlet conditions

A=

air

L=

droplet

P=

solid particles (without the evaporation process)

S=

parameter in swirling flow

T=

turbulent parameter

V=

steam

W=

parameter on the wall condition

m=

mean-mass parameter

Superscripts=
*=

parameter on the droplet surface

M=

mist

S=

swirling

Acronym=
LRN=

low-Reynolds number

RANS=

Reynolds averaged Navier–Stokes equations

SMC=

second moment closure

TKE=

turbulent kinetic energy

Nomenclature

Bi=

Biot number

b1D=

diffusion parameter of injection, determined with the use of the saturation curve

d=

droplet diameter, m

CD=

drag coefficient of the evaporating droplet, written by taking into account the deviation of the Stokes law

CP=

specific heat capacity, J kg−1 K−1

D=

diffusion coefficient, m2 s−1

G=

mass flow rate, kg s−1

J=

mass flux of steam from the surface of the evaporating droplet, kg m−2 s−1

ER=

heat transfer enhancement ratio

H=

step height, m

=

steam mass concentration at the “vapor-gas mixture–droplet” interface, corresponding to the saturation parameters at droplet temperature TL

KV0=

mass concentration of steam far from the droplet

k=

k = ⟨uiui′⟩/2 turbulent kinetic energy, m2 s−2

L=

latent heat of vaporization, J kg−1

ML=

droplets mass fraction

Nu=

Nusselt number

P=

pressure (Pa)

Pr=

Prandtl number

Pij=

stress production term, m s−2

2R1=

pipe diameter before sudden expansion, m

2R2=

pipe diameter after sudden expansion, m

Re=

Re = Um12R1/ν gas-phase Reynolds number

ReL=

ReL = |U − UL|d/ν Reynolds number of the dispersed phase

r=

radial coordinate, m

S=

swirl number

Sc=

Schmidt number

Sh=

Sherwood number

StD=

diffusional Stanton number

Stk=

Stk = τ/τf Stokes number in the mean motion

T=

temperature (K)

Ui=

components of mean velocity of the gas phase, m s−1

ULi=

components of droplets mean velocity, m s−1

u′⟩, ⟨v′⟩, ⟨w′⟩=

components of velocity pulsation, m s−1

uv′⟩=

turbulent Reynolds stresses, m2 s−2

uLiuLj′⟩=

kinetic stresses in the dispersed phase, m2 s−2

x=

axial coordinate, m

xR=

separation length, m

y=

distance normal to the wall, m

We=

Weber number

We*=

critical Weber number

Greek letters=
Φ=

Φ = MLρ/ρL volumetric concentration of the dispersed phase

ΓE=

turbulence macroscale, m

ΛT=

the turbulent macroscale time, s

ΩtL=

time of interaction with temperature pulsations of the carrying flow, s

ΩLag=

Lagrangian integral timescales, s

ΩE=

Eulerian integral timescales, s

Ωϵ=

time microscale of turbulence, s

αL=

heat transfer coefficient for the evaporating droplet, W m−2

β=

blending coefficient

ϵ=

dissipation of the turbulent kinetic energy, m2s−3

ϕij=

the velocity–pressure–gradient correlation (the pressure term), m s−2

μ=

dynamic viscosity, kg m−1 s−1

ν=

kinematic viscosity, m2 s−1

ρ=

density, kg m−3

σ=

surface tension, N m−1

τ=

dynamic relaxation time of droplets, s

τf=

turbulent timescale, s

Subscripts=
0=

parameter under conditions on the pipe axis

1=

parameter under inlet conditions

A=

air

L=

droplet

P=

solid particles (without the evaporation process)

S=

parameter in swirling flow

T=

turbulent parameter

V=

steam

W=

parameter on the wall condition

m=

mean-mass parameter

Superscripts=
*=

parameter on the droplet surface

M=

mist

S=

swirling

Acronym=
LRN=

low-Reynolds number

RANS=

Reynolds averaged Navier–Stokes equations

SMC=

second moment closure

TKE=

turbulent kinetic energy

Additional information

Funding

This work was supported by the grant of the Russian Ministry of Science & Education through Contract № 14.Z50.31.0041 dated 13th of February 2017.

Reprints and Corporate Permissions

Please note: Selecting permissions does not provide access to the full text of the article, please see our help page How do I view content?

To request a reprint or corporate permissions for this article, please click on the relevant link below:

Academic Permissions

Please note: Selecting permissions does not provide access to the full text of the article, please see our help page How do I view content?

Obtain permissions instantly via Rightslink by clicking on the button below:

If you are unable to obtain permissions via Rightslink, please complete and submit this Permissions form. For more information, please visit our Permissions help page.