DOCX

# Wall Boiling Models

7 views 0
Wall Boiling Models

17.5.16. Wall Boiling Models

17.5.16.1. Overview

The term “subcooled boiling” is used to describe the physical situation where the wall temperature is high enough to cause boiling to occur at the wall even though the bulk volume averaged liquid temperature is less than the saturation value. In such cases, the energy is transferred directly from the wall to the liquid. Part of this energy will cause the temperature of the liquid to increase and part will generate vapor. Interphase heat transfer will also cause the average liquid temperature to increase, however, the saturated vapor will condense. Additionally, some of the energy may be transferred directly from the wall to the vapor. These basic mechanisms are the foundations of the so called Rensselaer Polytechnic Institute (RPI) models.

In ANSYS FLUENT, the wall boiling models are developed in the context of the Eulerian multiphase model. The multiphase flows are governed by the conservation equations for phase continuity (Equation 17119), momentum

(Equation 17120), and energy (Equation 17126). The wall boiling

phenomenon is modeled by the RPI nucleate boiling model of Kurual and Podowski [190] and an extended formulation for the departed nucleate boiling regime (DNB) by Lavieville et al [200].

The wall boiling models are compatible with three different wall boundaries: isothermal wall, specified heat flux, and specified heat transfer coefficient (coupled wall boundary).

Specific submodels have been considered to account for the interfacial transfers of momentum, mass, and heat, as well as turbulence models in boiling flows, as described below.

To learn how to set up the boiling model, please refer to Including the

Boiling Model.

17.5.16.2. RPI Model

According to the basic RPI model, the total heat flux from the wall to the liquid is partitioned into three components, namely the convective heat flux, the quenching heat flux, and the evaporative heat flux:

(17274)

The heated wall surface is subdivided into area , which is covered by nucleating bubbles and a portion , which is covered by the fluid.

; The convective heat flux is expressed as

(17275)

where is the single phase heat transfer coefficient, and and

are the wall and liquid temperatures, respectively.

; The quenching heat flux models the cyclic averaged transient

energy transfer related to liquid filling the wall vicinity after

bubble detachment, and is expressed as

(17276)

Where is the conductivity, is the periodic time, and is the diffusivity.

; The evaporative flux is given by

(17277)

Where is the volume of the bubble based on the bubble departure diameter, is the active nucleate site density, is the vapor density, and is the latent heat of evaporation, and is the bubble departure frequency. These equations need closure for the following parameters:

; Area of Influence

Its definition is based on the departure diameter and the nucleate site density:

(17278)

Note that in order to avoid numerical instabilities due to unbound empirical correlations for the nucleate site density, the area of influence has to be restricted. The area of influence is limited as follows:

(17279)

The value of the empirical constant is usually set to 4, however it has been found that this value is not universal and may vary between 1.8 and 5. The following relation for this constant has also been implemented based on Del Valle and Kenning's findings [77]:

(17280)

and is the subcooled Jacob number defined as

(17281)

; Frequency of Bubble Departure

Implementation of the RPI model normally uses the frequency of bubble departure as the one based on inertia controlled growth (not really applicable to subcooled boiling) [64]

(17282)

; Nucleate Site Density

The nucleate site density is usually represented by a correlation based on the wall superheat. The general expression is of the form

(17283)

Here the empirical parameters from Lemmert and Chawla [205] are used,

where and . Other formulations are also available,

such as Kocamustafaogullari and Ishii [184] where

(17284)

Here

Where is the bubble departure diameter and the density function is defined as

(17285)

; Bubble Departure Diameter

The default bubble departure diameter (mm) for the RPI model is based on empirical correlations [190] and is calculated as

(17286)

while Kocamustafaogullari and Ishii [184] use

(17287)

with being the contact angle in degrees.

17.5.16.3. Non-equilibrium Subcooled Boiling

When using the basic RPI model, the temperature of the vapor is not calculated, instead it is fixed at the saturation temperature. In order to model different boiling regimes like DNB and critical heat flux, it is necessary to include the vapor temperature in the solution process. The wall heat partition is now modified as follows:

(17288)

Here is the diffusive heat flux of the vapor bubble phase,

is the vapor heat transfer coefficient based on turbulent wall functions. The function depends on the local liquid volume fraction with similar limiting values as the liquid volume fraction. Lavieville et al [200] proposed the following expression:

(17289)

Here, the critical value for the vapor fraction is

17.5.16.4. Interfacial Momentum Transfer

The interfacial momentum transfer may include four parts: drag, lift,

mass and turbulent drift forces (all described in Conservation virtual

Equations, Interphase Exchange Coefficients, and Turbulence Models. In

the wall boiling models, the virtual mass force is modeled using the standard correlation implemented in the Eulerian multiphase model within ANSYS FLUENT, while specific sub-models have been implemented for drag, lift, and turbulent drift forces. Also, user-defined options are available for both drag and lift forces.

17.5.16.4.1. Interfacial Area

The interfacial area is an important parameter for the drag and the heat transfer process. For dispersed boiling, the interfacial area, based on the diameter of the bubble, would be enough. However, as bubble coalescence takes place, this needs to be modified. The following options are included:

(17290)

(17291)

(17292)

17.5.16.4.2. Interfacial Drag Force

The interfacial drag force is calculated using the standard model described in Interphase Exchange Coefficients (and defined in the context

of the interfacial area in Equation 17290) is of the general form

(17293)

Where the drag coefficient is determined by choosing the minimum of the viscous regime and the distorted regime , defined as follows:

(17294)

The bubble diameter can be a constant value, a UDF, or a correlation function of local subcooling [190]:

(17295)

17.5.16.4.3. Interfacial Lift Force

The coefficient for the interfacial lift force is calculated using the correlation proposed by Moraga et al. [262]:

(17296)

Where . The lift coefficient combines the opposing action of two lift forces: the classical aerodynamics lift force resulting from interaction between bubble and liquid shear, and the lateral force resulting from interaction between bubble and vortexes shed by bubble wakes. Here is the bubble Reynolds number, and

is the bubble shear Reynolds number.

The formulation proposed by Tomiyama et al. [441] is also available with

the lift coefficient expressed as , where

(17297)

and

(17298)

where is the bubble Reynolds number and

(17299)

is the Eötvos number, with as the gravitational acceleration and

the surface tension number.

17.5.16.4.4. Turbulence Drift Force

In the ANSYS FLUENT Eulerian multiphase model, the general correlation for turbulence drift force (turbulent dispersion) is based on Simonin [351]. Due to numerical instabilities this force is now included in the Rhie & Chow interpolation [323] for the volume flux calculations.

Simonin’s [351] approach can be also used for the boiling model. However, for completeness of the RPI model, the default for the turbulent drift force is given by

(17300)

Where the turbulent dispersion coefficient is, by default, set to 1.0

17.5.16.5. Interfacial Heat Transfer

17.5.16.5.1. Vapor to Liquid Heat Transfer

As the bubbles depart from the wall and move towards the subcooled region, there is heat transfer from the bubble to the liquid, that is defined as

(17301)

Where is the interfacial area defined by Equation 17291 and

is the heat transfer coefficient based on the Ranz-Marshall correlation [315]

(17302)

17.5.16.5.2. Superheated Liquid to Vapor Heat Transfer The interface to vapor heat transfer is calculated using the constant time

method [200]. It is assumed that the vapor scale return to saturation

retains the saturation temperature by rapid evaporation/condensation. The formulation is as follows:

(17303)

Where is the time scale set to a default value of 0.05 and is the isobaric heat capacity.

17.5.16.6. Mass Transfer

17.5.16.6.1. Mass Transfer From the Wall to Vapor The evaporation mass flow is applied at the cell near the wall and it is

Equation 17303 derived from the evaporation heat flux,

(17304)

17.5.16.6.2. Interfacial Mass Transfer

The interfacial mass transfer depends directly on the interfacial heat transfer. Assuming that all the heat transferred to the interface is used in mass transfer (i.e. evaporation or condensation), the interfacial mass transfer rate can be written as:

(17305)

Report this document

For any questions or suggestions please email
cust-service@docsford.com