Publication Cover
Numerical Heat Transfer, Part A: Applications
An International Journal of Computation and Methodology
Volume 69, 2016 - Issue 12
581
Views
7
CrossRef citations to date
0
Altmetric
Original Articles

Modeling of compressible phase-change heat transfer in a Taylor-Bubble with application to pulsating heat pipe (PHP)

&
Pages 1355-1375 | Received 04 Sep 2015, Accepted 02 Dec 2015, Published online: 13 May 2016
 

ABSTRACT

The present work deals with the development of a compressible phase-change solver and implementation toward the numerical modeling and investigation of a part-unit cell of a pulsating heat pipe (PHP). The fundamental understanding of the working of the part-unit cell is imperative in the development of a complete Computational Fluid Dynamics (CFD) model of a PHP. The compressible model developed in the present work is based on the Volume-of-Fluid solver of the open source CFD software, OpenFOAM, in which the contour-based interface reconstruction algorithm and the contact-line evaporation model have been incorporated. Owing to the lack of a single standard benchmark validation case for a compressible phase-change solver, a huge emphasis in the present work is laid on the solver development and validation, the latter part of which is conducted in stages. Furthermore, simulations for the formation of a Taylor-Bubble through a constrained bubble growth are performed and the fallacy of an incompressible solver is shown distinctly. The validated solver is used to model a part-unit cell of a PHP and a parametric study is performed on the part-unit cell. The effect of variation of evaporator length, evaporator superheat, and liquid fill ratio on the performance of the PHP is discussed.

Nomenclature

a=

major-axis length of ellipse, m

b=

minor-axis length of ellipse, m

c=

specific heat capacity, J/(kg K)

D=

diameter, m

F=

volume fraction of the liquid phase, -

Fs=

surface tension force, N

f1=

objective function 1, -

g=

acceleration due to gravity, m/s2

=

phase-change mass flux, kg/(m2 s)

k=

thermal conductivity, W/(m K)

L=

length, m

p=

pressure, Pa

R=

heat resistance, K m2/W

=

radius, m

s=

second, s

T=

temperature, K

t=

time, s

U=

velocity vector, m/s

u=

magnitude of velocity, m/s

V=

volume, m3

Δh=

change in enthalpy, J/kg

Δt=

time step, s

Δτ=

artificial time step, s

Δx=

length interval, m

ρ=

density, kg/m3

=

sharp mass source term field, kg/(m3s)

=

smeared mass source term field, kg/(m3s)

γ=

polytropic constant, -

μ=

dynamic viscosity, Pa s

σ=

surface tension, N/m

κ=

interface curvature, 1/m

ϕ=

any arbitrary vector field, -

Subscripts=
A=

adiabatic section

E=

evaporator

F=

liquid fill region

I=

component/phase

L=

liquid

V=

vapor

Int=

interface

Lv=

difference of the parameter between liquid and vapor phases

Max=

maximum

Sat=

saturation

=

infinity

Nomenclature

a=

major-axis length of ellipse, m

b=

minor-axis length of ellipse, m

c=

specific heat capacity, J/(kg K)

D=

diameter, m

F=

volume fraction of the liquid phase, -

Fs=

surface tension force, N

f1=

objective function 1, -

g=

acceleration due to gravity, m/s2

=

phase-change mass flux, kg/(m2 s)

k=

thermal conductivity, W/(m K)

L=

length, m

p=

pressure, Pa

R=

heat resistance, K m2/W

=

radius, m

s=

second, s

T=

temperature, K

t=

time, s

U=

velocity vector, m/s

u=

magnitude of velocity, m/s

V=

volume, m3

Δh=

change in enthalpy, J/kg

Δt=

time step, s

Δτ=

artificial time step, s

Δx=

length interval, m

ρ=

density, kg/m3

=

sharp mass source term field, kg/(m3s)

=

smeared mass source term field, kg/(m3s)

γ=

polytropic constant, -

μ=

dynamic viscosity, Pa s

σ=

surface tension, N/m

κ=

interface curvature, 1/m

ϕ=

any arbitrary vector field, -

Subscripts=
A=

adiabatic section

E=

evaporator

F=

liquid fill region

I=

component/phase

L=

liquid

V=

vapor

Int=

interface

Lv=

difference of the parameter between liquid and vapor phases

Max=

maximum

Sat=

saturation

=

infinity

Acknowledgments

The authors would like to acknowledge the technical support provided by Prof. Dr.-Ing. Peter Stephan and Dr.-Ing. Stefan Batzdorf of Technische Universitat Darmstadt toward successful completion of the present work.

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.