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

Title: A semi-analytical method for simulating matrix diffusion in numerical transport models

Publication Date:
Sponsoring Org.:
OSTI Identifier:
Grant/Contract Number:
Resource Type:
Journal Article: Publisher's Accepted Manuscript
Journal Name:
Journal of Contaminant Hydrology
Additional Journal Information:
Journal Volume: 197; Journal Issue: C; Related Information: CHORUS Timestamp: 2017-12-18 18:51:48; Journal ID: ISSN 0169-7722
Country of Publication:

Citation Formats

Falta, Ronald W., and Wang, Wenwen. A semi-analytical method for simulating matrix diffusion in numerical transport models. Netherlands: N. p., 2017. Web. doi:10.1016/j.jconhyd.2016.12.007.
Falta, Ronald W., & Wang, Wenwen. A semi-analytical method for simulating matrix diffusion in numerical transport models. Netherlands. doi:10.1016/j.jconhyd.2016.12.007.
Falta, Ronald W., and Wang, Wenwen. Wed . "A semi-analytical method for simulating matrix diffusion in numerical transport models". Netherlands. doi:10.1016/j.jconhyd.2016.12.007.
title = {A semi-analytical method for simulating matrix diffusion in numerical transport models},
author = {Falta, Ronald W. and Wang, Wenwen},
abstractNote = {},
doi = {10.1016/j.jconhyd.2016.12.007},
journal = {Journal of Contaminant Hydrology},
number = C,
volume = 197,
place = {Netherlands},
year = {Wed Feb 01 00:00:00 EST 2017},
month = {Wed Feb 01 00:00:00 EST 2017}

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record at 10.1016/j.jconhyd.2016.12.007

Citation Metrics:
Cited by: 1work
Citation information provided by
Web of Science

Save / Share:
  • We study a system with bimolecular irreversible kinetic reaction A+B→∅ where the underlying transport of reactants is governed by diffusion, and the local reaction term is given by the law of mass action. We consider the case where the initial concentrations are given in terms of an average and a white noise perturbation. Our goal is to solve the diffusion–reaction equation which governs the system, and we tackle it with both analytical and numerical approaches. To obtain an analytical solution, we develop the equations of moments and solve them approximately. To obtain a numerical solution, we develop a grid-less Montemore » Carlo particle tracking approach, where diffusion is modeled by a random walk of the particles, and reaction is modeled by annihilation of particles. The probability of annihilation is derived analytically from the particles' co-location probability. We rigorously derive the relationship between the initial number of particles in the system and the amplitude of white noise represented by that number. This enables us to compare the particle simulations and the approximate analytical solution and offer an explanation of the late time discrepancies. - Graphical abstract:.« less
  • The computational modeling of the biodegradation of contaminated groundwater systems accounting for biochemical reactions coupled to contaminant transport is a valuable tool for both the field engineer/planner with limited computational resources and the expert computational researcher less constrained by time and computer power. There exists several analytical and numerical computer models that have been and are being developed to cover the practical needs put forth by users to fulfill this spectrum of computational demands. Generally, analytical models provide rapid and convenient screening tools running on very limited computational power, while numerical models can provide more detailed information with consequent requirementsmore » of greater computational time and effort. While these analytical and numerical computer models can provide accurate and adequate information to produce defensible remediation strategies, decisions based on inadequate modeling output or on over-analysis can have costly and risky consequences. In this chapter we consider both analytical and numerical modeling approaches to biodegradation and reactive transport. Both approaches are discussed and analyzed in terms of achieving bioremediation goals, recognizing that there is always a tradeoff between computational cost and the resolution of simulated systems.« less
  • Perturbative experiments in magnetically confined fusion plasmas have shown that edge cold pulses travel to the centre of the device on a time scale much faster than expected on the basis of diffusive transport. An open issue is whether the observed fast pulse propagation is due to non-local transport mechanisms or if it could be explained on the basis of local transport models. To elucidate this distinction, perturbative experiments involving ICRH power modulation in addition to cold pulses have been conducted in JET for the same plasma. Local transport models have found problematic the reconciliation of the fast propagation ofmore » cold pulses with the comparatively slower propagation of heat waves generated by power modulation. In this paper, a non-local model based on the use of fractional diffusion operators is used to describe these experiments. A numerical study of the parameter dependence of the pulse speed and the amplitude and phase of the heat wave is also presented.« less
  • Research on solving inverse problems governed by the linear integrodifferential transport equation has been around for >30 yr. Inverse radiative transfer problems in clouds is a very important topic in environmental sciences. Moreover, in cancer therapy, it is usually desirable for the absorbed dose to be as uniform as possible within the tumor; i.e., the existence of a hot spot inside the tumor usually indicates that a more optimal treatment plan exists-one in which a lesser dose is deposited outside the tumor. In the field of nuclear geophysics, inverse transport methods play an important role in nuclear oil well logging.more » In connection with these applications, the authors describe an analytical numerical method applied to two different types of inverse S{sub N} transport problems: (a) boundary condition estimation and (b) interior source estimation.« less