scispace - formally typeset

Journal ArticleDOI

Double Diffusive Natural Convection in a Nuclear Waste Repository

01 Jul 2008-Nuclear Technology (American Nuclear Society)-Vol. 163, Iss: 1, pp 38-46

Abstract: In this study, we conduct a two-dimensional numerical analysis of double diffusive natural convection in an emplacement drift for a nuclear waste repository. In-drift heat and moisture transport is driven by combined thermal- and compositional-induced buoyancy forces. Numerical results demonstrate buoyancy-driven convective flow patterns and configurations during both repository heat-up and cool-down phases. It is also shown that boundary conditions, particularly on the drip-shield surface, have strong impacts on the in-drift convective flow and transport.
Topics: Combined forced and natural convection (61%), Natural convection (59%), Convection (56%), Heat transfer (53%), Thermal (50%)

Content maybe subject to copyright    Report

Double Diffusive Natural Convection in a Nuclear Waste Repository
Yue Hao, John
J.
Nitao, Thomas A. Buscheck, and Yunwei Sun
L-631, Lawrence Livermore National Laboratory, Livermore,
CA
94551
Abstract
-
In this study, we conduct a two-dimensional numerical analysis of double diffusive natural convection in an
er?lplacement drift for a nuclear waste repository. In-drift heat and moisture transport is driven by combined thermal- and
conlpositional-induced buoyancy forces. Numerical results demonstrate buoyancy-driven convective flow patterns and
configurations during both repository heat-up and cool-down phases. It is also shown that bounda ry conditions, particularly
on the drip-shield
surfnce, have strong impacts on the in-drift convective flow and transport.
I. INTRODUCTION Stokes equations. The governing equations are expressed
in the following mass, momentum, energy, and vapor
It is recognized that natural convection plays an
concentration conservation forms:
important role in
many heat and mass transfer vrocesses.
In
ionisothermal appfications with binary fluid mixtures,
the interaction of natural heat convection
with.mass
transport of two components results in a complex flow
and transport phenomenon called double
diffusive
convection. Double diffusive convection, resulting from
buoyancy forces caused by temperature and
compositional gradients, is found in many natural-system
and engineering applications. For this reason it has
attracted considerable attention
[l, 21. One potential
example is related to the transport of heat and moisture
inside emplacement drifts in a nuclear-waste repository in
the unsaturated zone. After the emplacement of waste
packages, the radioactive heat of decay generates water
vapor due. to the evaporation of water in the adjoining
host rock, which migrates into the drift. Within the drift,
natural convection contributes to transport of water vapor
from hotter to cooler locations, where it may condense.
Continuity:
Momentum:
where
t,
p,
u,
p,
and
g
represent time, binary mixture
density, velocity, pressure, and body force, and
z
is the
viscous stress tensor,
in which
p
is fluid viscosity and the superscript
T
denotes
the transpose.
Energy:
In addition to affecting heat and mass transport, the
complicated flow patterns and structures induced by
(4)
combined thermaland compositional buoyancy effects
will influence evaporation and condensation on
engineered material surfaces within the drift. Because of
the potential for influencing the corrosion of drip shields
and waste packages, moisture condensation within drifts
is of concern for total system performance assessment of
the repository. Moreover, analyses of in-drift heat and
moisture convection can
ljrovide a better understanding of
the basic physics of flow and transport phenomena inside
engineered tunnels. Recent CFD models that apply the
FLUENT code have been used to describe in-drift flows,
with a focus on thermal-induced natural convection and
determining effective dispersion coefficients for models
predicting moisture transport and condensation
[3,4].
The .goal of this work is to investigate double-
diffusive convection inside an emplacement drift, using a
Navier-Stokes model approach to examine and capture
the predominant convection modes.
with
T,
C,,,
k,
and
Q
as temperature, specific heat, thermal
conductivity, and heat generation term.
Vapor concentration (expressed as mass fraction) in a
I
vapor-air mixture:
C, D,
and
R
are vapor mass fraction, binary diffusion
coefficient, and vapor source term, respectively.
For this study, the total pressure of the
airlvapor
mixture
p,
is assumed spatially uniform, and the binary
mixture obeys the ideal gas law. Based on these
conditions, the fluid density can be expressed as
,
.
11.
MATHEMATICAL FORMULATION
with
T
as temperature,
R
as universal gas constant,
My
as
In this study, natural heat and mass convection in an
molecular mass of water vapor and
MA
as molecular mass
.air/vapor binary system is considered, using Navier-
of air.

If both temperature and vapor-concentration
variations are small, then the Boussinesq approximation is
applied with the fluid density simplified as
in which subscript
oo
denotes the reference state, and
BT
and
/Ic
are the thermal and concentration expansion
coefficients, respectively,
For thermal natural convection, we use the dimensionless
Rayleigh number
Ra,
=
TL3
and Prandtl number
For mass transfer, we use the corresponding Rayleigh
number
Ra,
=
g&AcL3
VD
and Schmidt number
Here AT, AC,
u,
a,
and L are temperature difference,
vapor-concentration difference, fluid kinematics viscosity,
thermal diffusivity, and length scale, respectively.
In order to compare the magnitudes of thermal- and
compositional-induced natural convections, the so called
"buoyancy ratio" is introduced by
[5],
The dominant driving buoyant force is determined by the
buoyancy ratio N. It is obvious that as N decreases,
thermal buoyancy effects will become dominant over
compositional effects. The
k-w
turbulent model [6] is
used to account for turbulent flow effects for large
Rayleigh-number problems. The above governing
equations are solved by the Navier-Stokes module
implemented in the NUFT code
[7], which has been
validated against benchmark
problemS.
Ill.
RESULTS
AND
DISCUSSION
This study addresses both heat and vapor transport
within an emplacement drift (Fig la). As depicted in Fig.
I
b, the drip shield and waste package are lumped together
as a monolithic heat source that is impermeable with
respect to mass transport. We model a two-dimensional
in-drift flow domain (Fig. la), which is bounded by the
outer drip-shield surface, upper invert surface, and
drift-
wall surface above the invert. Due to the symmetry of the
problem, only half of the drift needs to be represented in
the numerical model (Fig. lb).
For the problem shown in Fig. 1, in-drift flow and
transport processes are highly dependent on
thermal-
hydrological conditions in the adjoining porous host rock.
Hence the key aspect of modeling of in-drift flow is the
specification of boundary conditions on the interfaces
between the fieelopen-flow system in the drift and the
porous-flow system in the host rock. Ideally one would
determine these conditions by directly coupling
freelopen
flow with porous flow. However, this is beyond the scope
of the current study. For the purpose of this study, we
assign appropriate boundary conditions on the surfaces of
the drift wall, invert, and drip shield on the basis of results
fiom the Multiscale Thermohydrologic Model (MSTHM)
[8,9, 101. The MSTHM, which is based on porous-
medium Darcy-flow approximations, predicts the coupled
thermal-hydrological conditions at both drift-scale and
mountain-scale. Fig. 2 plots the time history of the
representative drift-wall temperature and vapor
concentration obtained fiom the MSTHM
LDTH-
submodel simulations.
During the initial 50-year preclosure ventilation
0
period, the heat generated by waste packages is removed
by forced convective cooling. After drift ventilation
ceases, the postclosure period begins and the drift-wall
temperature abruptly rises (Fig. 2a) to well above the
local boiling point of water
(96°C). The vapor
concentration at the drift-wall surfaces also sharply
increases along with temperatures (Fig. 2b). The
dryout
phase lasts over 1000 years until the drift wall cools down
to
96OC (Fig. 2b).
For simplicity, the temperature
Tdw, Tbvert, and Tds are
uniformly imposed along the drift-wall surface, upper
invert surface, and drip-shield surface (Fig. lb). In
addition, the surfaces of the drift wall and invert are
assumed to be permeable with constant vapor
concentrations
Cdw and ChVe,. In order to explore all the
possible flow patterns and convection modes inside the
drift, we investigate two possible conditions on the outer
drip-shield surface with respect to mass transfer.
The simulated transient behavior discussed in the
following sections is not intended to be exactly
representative of a real repository system. The primary
purpose of this study is to illustrate the influence of
boundary conditions on double diffusive natural
convection within emplacement drifts, with a focus on the
interactionbetween the thermal and compositional
buoyancy forces. These simulations neglect phase change
(evaporation and condensation) within the
dnft.
Moreover, the transient aspects of the simulated behavior
are presented to illustrate the interaction between the
thermal and compositional buoyancy forces.




Citations
More filters

Journal ArticleDOI
Abstract: In this article, convection driven by combined thermal and solutal concentration buoyancy effects in a lid-driven square cavity is examined using velocity-vorticity form of Navier-Stokes equations. The governing equations consist of vorticity transport equation, velocity Poisson equations, energy equation, and concentration equation. Validation results are discussed for convection due to heat and mass transfer in a lid-driven square cavity at Re = 500, Le = 2, and GRT = GRS = 100. These results indicate that the present velocity-vorticity formulation could predict the characteristic parameters of flow, temperature, and solutal concentration fields using a much coarser mesh compared to the mesh used in a stream function-vorticity formulation. The capability of the proposed algorithm to handle complex geometry is demonstrated by application to mixed convection in a lid-driven square cavity with a square blockage. The effect of buoyancy ratio on the convection phenomenon is discussed for buoyancy ratio var...

30 citations


Journal ArticleDOI
Antonio Soto-Meca1, J. Serna1, F.J.S. Velasco1Institutions (1)
TL;DR: Results show that the pulsating regime enhances heat/mass transfer within a square cavity up to a 14%/38% respectively with respect to the non-pulsating case, due to the promotion of additional shear stress fields.
Abstract: This work assess the effectiveness of the pulsating flow regime for the enhancement of heat and mass transfer in double diffusive problems by studying a two-dimensional, heated, lid cavity. The research characterizes the influence of pulsating parameters such as temporal frequency, wave number and amplitude in the process. Results show that the pulsating regime enhances heat/mass transfer within a square cavity up to a 14%/38% respectively with respect to the non-pulsating case, due to the promotion of additional shear stress fields. As Richardson number, Brownian diffusion or the solute/thermal buoyancy ratio increase or the cavity becomes narrower, heat/mass transfer increases.

11 citations


Journal ArticleDOI
Abstract: In this paper, detailed uncertainty propagation analysis (UPA) and variance-based global sensitivity analysis (GSA) are performed on the widely adopted double-diffuse convection (DDC) benchmark problem of a square porous cavity with horizontal temperature and concentration gradients. The objective is to understand the impact of uncertainties related to model parameters on metrics characterizing flow, heat and mass transfer processes, and to derive spatial maps of uncertainty and sensitivity indices which can provide physical insights and a better understanding of DDC processes in porous media. DDC simulations are computationally expensive and UPA and GSA require large number of simulations, so an appropriate strategy is developed to reduce the computational burden. The approach is built on two pillars: (a) an efficient numerical simulator based on the Fourier series method that generates training data, and (b) polynomial chaos expansion (PCE) meta-models that are trained using the simulator data, and then replace the numerical model in UPA and GSA. Assuming that the Rayleigh number (Ra), the solutal to thermal buoyancy ratio (Nb) and the Lewis number (Le) are the uncertain input variables, the results of UPA show that the zones of high temperature and concentration variability are located in the regions where the flow is mainly driven by the buoyancy effects. GSA indicates that Nb is the most influential parameter affecting the temperature and concentration fields, followed respectively by Ra and Le. For the heat-driven flow case ( N b > − 1 ), the concentration field is more influenced by Le than Ra. For deeper understanding of uncertainty propagation, we estimate the bias introduced by replacing uncertain parameters by deterministic values. The resulting spatial maps of the difference between deterministic output and stochastic mean show that a deterministic approach leads to different zones where the temperature, concentration and velocity fields can be either overestimated or underestimated. The conclusions drawn in this work are likely to be helpful in different applications involving DDC in porous enclosures leading to convective circulation cells.

8 citations


Journal ArticleDOI
Abstract: We investigate natural convection driven by a horizontal temperature gradient and a vertical concentration gradient in fluid-filled enclosures with obstructions inside it. Within the domain, nine adiabatic and impermeable cylinders are placed, occupying 30% of the domain volume. The Boussinesq approximation is used to account for density variations within the fluid and the flow is fully resolved. The solutal Rayleigh number has been fixed at R a C = 10 6 and the Prandtl number at P r = 5.4 . The Lewis number has been varied in the range of 1 ≤ Le ≤ 100 and the buoyancy ratio in the range of 0.1 ≤ |N| ≤ 10. The rate of heat and mass transfer are compared to those found in single-scalar natural convection, i.e solely thermal or concentration driven convection. Besides, the obtained heat and mass transfer rate in the cylinder-packed enclosure have been compared to those found in a fluid-only domain. We observe that the addition of a destabilizing concentration gradient to a side-heated enclosure results in heat transfer enhancement, which decreases with Lewis number and thermal Rayleigh number. Similarly, the temperature gradient increases the mass transfer, especially at high Lewis numbers and lower concentration buoyancy force over its thermal counterpart. Although the presence of the cylindrical obstacles reduced the flow velocity, the mass transfer was enhanced at lower buoyancy ratio.

5 citations


Journal ArticleDOI
Abstract: The nuclear waste storage concept according to the baseline design of the proposed high-level nuclear waste repository at Yucca Mountain is analyzed. The hightemperature storage concept, in which the emplacement area is heated above the boiling temperature of water, is subject to criticism on the basis of uncertainties due to nonlinear multiphysics processes in the rock mass and in the storage airspace. The storage environment around the nuclear waste containers is reexamined using a new thermal-hydrologic airflow model. The complex nature ofthethermal-hydraulicbehaviorinasuperheatedwaste repository is described with fewer simplifying assumptions than those used in the baseline design. The emplacement area in the mountain is described as an open system, in which the air pressure is connected to the barometric pressure through fractures, faults, and partially sealed drifts. The cyclic variation of the atmospheric pressure that affects the heat and mass transport processes in the near-field rock mass is also modeled. The implications of evaporation into the drift airspace are discussed, and a hypothesis of salt accumulation in the near-field rock mass is established. Model calculation is also presented for a below–boiling temperature storage concept that is easier to predict and has fewer anomalies. The price for a below–boiling temperature storage is the extended preclosure ventilation time period. However, as demonstrated for a trade-off, it is possible to design a repository with below-boiling temperatures and doubled waste inventory at the same time.

5 citations


References
More filters

Book
01 Jan 1993
Abstract: Keywords: ecoulement : compressible Note: + disquette Reference Record created on 2005-11-18, modified on 2016-08-08

6,811 citations


Journal ArticleDOI
Benjamin Gebhart1, L Pera1Institutions (1)
Abstract: This paper concerns the laminar flows which arise in fluids due to the interaction of the force of gravity and density differences caused by the simultaneous diffusion of thermal energy and of chemical species. Species concentration levels are assumed small, as is typical for many processes in water and in atmospheric air. The usual Boussinesq approximations yield a set of equations which are shown to have solutions of similarity form for combined buoyancy effects, for vertical flows adjacent to surfaces and in plumes. This similarity is of the same form as that found for single buoyancy mechanism flows. The resulting equations were integrated for air and water for various practical values of the Schmidt number and for multiple buoyancy effects aiding and opposing. The results show many interesting effects on velocity, heat and mass transfer, and on laminar stability. A comparison of the results with those of integral method analysis shows the limits and reasons for failure of these approximate calculations in the more complicated of such combined buoyancy mechanism flows.

525 citations


Journal ArticleDOI
C. Béghein1, Fariborz Haghighat2, F. Allard1Institutions (2)
Abstract: Steady-state thermosolutal convection in a square cavity filled with air, submitted to horizontal temperature and concentration gradients, is studied numerically. In the first series of numerical simulations, the influence of solutal buoyancy force on heat or mass transfer rate is investigated: Lewis and thermal Rayleigh numbers are kept constant (Le = 1, RaT = 107), solutal Rayleigh number is varied (RaS = 105−5 × 107). The second series deals with the influence of Lewis number on fluid motion for heat transfer driven flow (RaT = 107, RaS = 0) and mass transfer driven flow (RaT = 0, RaS = 107) configurations. Lewis number is varied from 0.3 to 5. Correlations are obtained between heat and mass transfer rates and the non-dimensional numbers characterizing both phenomena.

208 citations


S. Ostrach1Institutions (1)
01 Jan 1980
Abstract: The problem of free and natural convection with combined driving forces is considered in general and all possible configurations are identified. Dimensionless parameters are discussed in order to help categorize the various problems, and existing work is critically evaluated. Four distinct cases are considered for conventional convection and for the situation when the body force and the density gradient are parallel but opposed. Considerable emphasis is given to unstable convection in horizontal layers.

130 citations


Journal ArticleDOI
Abstract: The effect of buoyancy ratio on the flow structure is investigated numerically for a binary mixture gas in a rectangular enclosure subject to opposing horizontal thermal and compositional buoyancies. The following conditions were considered: RaT = 105, Pr = 1, Le = 2 and N = 0.0–2.0 for A = 2. The numerical solution predicts that oscillatory double-diffusive convection with the secondary cell flow structure occurs for a certain range of buoyancy ratio. The key mechanism for oscillatory flow is that the unstably stratified region of species shifts from the central part of the enclosure to the upper and lower parts, and vice versa in a time-periodic sense, due to the interaction of heat and mass transfer with different diffusivities near the vertical walls. Bifurcation structures of the oscillatory flow in the present system are discussed.

103 citations