skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Revisiting the Fundamental Analytical Solutions of Heat and Mass Transfer: The Kernel of Multirate and Multidimensional Diffusion

Abstract

There are two types of analytical solutions of temperature/concentration in and heat/mass transfer through boundaries of regularly shaped 1D, 2D, and 3D blocks. These infinite-series solutions with either error functions or exponentials exhibit highly irregular but complementary convergence at different dimensionless times, t d0. In this paper, approximate solutions were developed by combining the error-function-series solutions for early times and the exponential-series solutions for late times and by using time partitioning at the switchover time, t d0. The combined solutions contain either the leading term of both series for normal-accuracy approximations (with less than 0.003 relative error) or the first two terms for high-accuracy approximations (with less than 10-7 relative error) for 1D isotropic (spheres, cylinders, slabs) and 2D/3D rectangular blocks (squares, cubes, rectangles, and rectangular parallelepipeds). This rapid and uniform convergence for rectangular blocks was achieved by employing the same time partitioning with individual dimensionless times for different directions and the product of their combined 1D slab solutions. The switchover dimensionless time was determined to minimize the maximum approximation errors. Furthermore, the analytical solutions of first-order heat/mass flux for 2D/3D rectangular blocks were derived for normal-accuracy approximations. These flux equations contain the early-time solution with a three-term polynomial inmore » √td and the late-time solution with the limited-term exponentials for rectangular blocks. The heat/mass flux equations and the combined temperature/concentration solutions form the ultimate kernel for fast simulations of multirate and multidimensional heat/mass transfer in porous/fractured media with millions of low-permeability blocks of varying shapes and sizes.« less

Authors:
ORCiD logo [1];  [1];  [1];  [1]
  1. Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States)
Publication Date:
Research Org.:
Montana State Univ., Bozeman, MT (United States)
Sponsoring Org.:
USDOE Office of Fossil Energy (FE)
OSTI Identifier:
1430907
Alternate Identifier(s):
OSTI ID: 1409860
Grant/Contract Number:
FC26-05NT42587; AC02-05CH11231; FE0023323
Resource Type:
Journal Article: Accepted Manuscript
Journal Name:
Water Resources Research
Additional Journal Information:
Journal Volume: 53; Journal Issue: 11; Journal ID: ISSN 0043-1397
Publisher:
American Geophysical Union (AGU)
Country of Publication:
United States
Language:
English
Subject:
42 ENGINEERING

Citation Formats

Zhou, Quanlin, Oldenburg, Curtis M., Rutqvist, Jonny, and Birkholzer, Jens T. Revisiting the Fundamental Analytical Solutions of Heat and Mass Transfer: The Kernel of Multirate and Multidimensional Diffusion. United States: N. p., 2017. Web. doi:10.1002/2017WR021040.
Zhou, Quanlin, Oldenburg, Curtis M., Rutqvist, Jonny, & Birkholzer, Jens T. Revisiting the Fundamental Analytical Solutions of Heat and Mass Transfer: The Kernel of Multirate and Multidimensional Diffusion. United States. doi:10.1002/2017WR021040.
Zhou, Quanlin, Oldenburg, Curtis M., Rutqvist, Jonny, and Birkholzer, Jens T. Tue . "Revisiting the Fundamental Analytical Solutions of Heat and Mass Transfer: The Kernel of Multirate and Multidimensional Diffusion". United States. doi:10.1002/2017WR021040.
@article{osti_1430907,
title = {Revisiting the Fundamental Analytical Solutions of Heat and Mass Transfer: The Kernel of Multirate and Multidimensional Diffusion},
author = {Zhou, Quanlin and Oldenburg, Curtis M. and Rutqvist, Jonny and Birkholzer, Jens T.},
abstractNote = {There are two types of analytical solutions of temperature/concentration in and heat/mass transfer through boundaries of regularly shaped 1D, 2D, and 3D blocks. These infinite-series solutions with either error functions or exponentials exhibit highly irregular but complementary convergence at different dimensionless times, td0. In this paper, approximate solutions were developed by combining the error-function-series solutions for early times and the exponential-series solutions for late times and by using time partitioning at the switchover time, td0. The combined solutions contain either the leading term of both series for normal-accuracy approximations (with less than 0.003 relative error) or the first two terms for high-accuracy approximations (with less than 10-7 relative error) for 1D isotropic (spheres, cylinders, slabs) and 2D/3D rectangular blocks (squares, cubes, rectangles, and rectangular parallelepipeds). This rapid and uniform convergence for rectangular blocks was achieved by employing the same time partitioning with individual dimensionless times for different directions and the product of their combined 1D slab solutions. The switchover dimensionless time was determined to minimize the maximum approximation errors. Furthermore, the analytical solutions of first-order heat/mass flux for 2D/3D rectangular blocks were derived for normal-accuracy approximations. These flux equations contain the early-time solution with a three-term polynomial in √td and the late-time solution with the limited-term exponentials for rectangular blocks. The heat/mass flux equations and the combined temperature/concentration solutions form the ultimate kernel for fast simulations of multirate and multidimensional heat/mass transfer in porous/fractured media with millions of low-permeability blocks of varying shapes and sizes.},
doi = {10.1002/2017WR021040},
journal = {Water Resources Research},
number = 11,
volume = 53,
place = {United States},
year = {Tue Oct 24 00:00:00 EDT 2017},
month = {Tue Oct 24 00:00:00 EDT 2017}
}

Journal Article:
Free Publicly Available Full Text
This content will become publicly available on October 24, 2018
Publisher's Version of Record

Save / Share:
  • A field-scale reactive transport model was developed that incorporates laboratory-characterized U(VI) surface complexation reactions (SCR) and multi-rate mass transfer processes, and field-measured hydrogeochemical conditions at Department of Energy, Hanford 300A site, Washington, where an Integrated Field Research Challenge project is ongoing. The model was used to assess the importance of multi-rate mass transfer processes on sorption-retarded U(VI) reactive transport at the 300A site and to evaluate the effect of variable geochemical conditions on U(VI) plume migration caused by dynamic river stage fluctuations at the east side of the site. Model simulations revealed a complex spatio-temporal variations of groundwater geochemistry thatmore » affects U(VI) speciation, adsorption, and plume migration. In general, the river water intrusion enhances uranium adsorption and lowers groundwater aqueous uranium concentration as a result of river water dilution that decreases aqueous carbonate concentration, which subsequently weakens aqueous U(VI)-carbonate complexation and enhances U(VI)-surface complexation. The simulations also found that SCR-retarded U migration becomes more dynamic and more in sync with the groundwater flow field when multi-rate mass transfer processes are involved. Strong U(VI) adsorption was simulated at the 300A site based on the field-measured hydrogeochemical conditions, suggesting a slow dissipation of U(VI) plume, a phenomenon consistent with the observation at the site. Uranium breakthrough curves at selected observation points and the mass changes over time in the simulation indicate that uranium adsorption/desorption never attains steady state as a result of both the highly dynamic flow field and the chemistry variations caused by river water intrusion. Thus, the multi-rate SCM model appears to be a crucial feature for future reactive transport simulations of uranium at the 300A site.« less
  • An analytical solution is presented for the movement of chemicals through a sorbing porous medium with lateral or intra-aggregate diffusion. The liquid phase in the porous medium is divided into mobile and immobile regions. Diffusional transfer between the 2 liquid regions is assumed to be proportional to the concentration difference between the mobile and immobile liquids. Sorption processes in both the dynamic and stagnant regions of the medium are assumed to be instantaneous and adsorption isotherm is assumed to be linear. The analytical model derived describes the extensive tailing observed during the flow through an unsaturated, aggregated sorbing medium andmore » explains the often observed early breakthrough of chemicals in the effluent. (26 refs.)« less
  • The authors have presented an analytical description of the temperature and kinetic curves for adsorbent granules in the form of plates, and they have proposed methods for the determination of the mass- and heat-transfer parameters based on experimental temperature curves. In the case of the Xe-CaA system the adsorption process is controlled by the external heat transfer.
  • The physical processes that occur typically within an oil sand bed are considered when the bed is subjected to a hot gaseous stream. In this study, the extent of fluid volatilization was obtained from a consideration of the simultaneous heat and mass transfer processes within the oil sands. The resulting system of equations together with the boundary conditions were solved numerically using an implicit finite difference method. The transient fluid concentration and temperature distributions within the oil sand bed were then obtained under a wide range of operating conditions. The resulting theoretical rates of volatilization and temperatures show generally goodmore » agreement with corresponding experimental values that were obtained for the purpose.« less