Computational Bayesian maximum entropy solution of a stochastic advection‐reaction equation in the light of site‐specific information

This work presents a computational formulation of the Bayesian maximum entropy (BME) approach to solve a stochastic partial differential equation (PDE) representing the advection‐reaction process across space and time. The solution approach provided by BME has some important features that distinguish it from most standard stochastic PDE techniques. In addition to the physical law, the BME solution can assimilate other sources of general and site‐specific knowledge, including multiple‐point nonlinear space/time statistics, hard measurements, and various forms of uncertain (soft) information. There is no need to explicitly solve the moment equations of the advection‐reaction law since BME allows the information contained in them to consolidate within the general knowledge base at the structural (prior) stage of the analysis. No restrictions are posed on the shape of the underlying probability distributions or the space/time pattern of the contaminant process. Solutions of nonlinear systems of equations are obtained in four space/time dimensions and efficient computational schemes are introduced to cope with complexity. The BME solution at the prior stage is in excellent agreement with the exact analytical solution obtained in a controlled environment for comparison purposes. The prior solution is further improved at the integration (posterior) BME stage by assimilating uncertain information at the data points as well as at the solution grid nodes themselves, thus leading to the final solution of the advection‐reaction law in the form of the probability distribution of possible concentration values at each space/time grid node. This is the most complete way of describing a stochastic solution and provides considerable flexibility concerning the choice of the concentration realization that is more representative of the physical situation. Numerical experiments demonstrated a high solution accuracy of the computational BME approach. The BME approach can benefit from the use of parallel processing (the relevant systems of equations can be processed simultaneously at each grid node and multiple integrals calculations can be accelerated significantly, etc.).

[1]  Athanasios Papoulis,et al.  Probability, Random Variables and Stochastic Processes , 1965 .

[2]  James M. Ortega,et al.  Iterative solution of nonlinear equations in several variables , 2014, Computer science and applied mathematics.

[3]  D. Brillinger Introduction to random differential equations and their applications , 1971 .

[4]  Samuel D. Conte,et al.  Elementary Numerical Analysis: An Algorithmic Approach , 1975 .

[5]  M. Golberg Solution Methods for Integral Equations , 1979 .

[6]  D. Jordan,et al.  Nonlinear Ordinary Differential Equations: An Introduction for Scientists and Engineers , 1979 .

[7]  Bruno O. Shubert,et al.  Random variables and stochastic processes , 1979 .

[8]  Iraj Javandel,et al.  Groundwater Transport: Handbook of Mathematical Models , 1984 .

[9]  John G. Proakis,et al.  Probability, random variables and stochastic processes , 1985, IEEE Trans. Acoust. Speech Signal Process..

[10]  P. Kitanidis Parameter Uncertainty in Estimation of Spatial Functions: Bayesian Analysis , 1986 .

[11]  Joseph E. Flaherty,et al.  Adaptive Methods for Partial Differential Equations , 1989 .

[12]  J. Strikwerda Finite Difference Schemes and Partial Differential Equations , 1989 .

[13]  G. Christakos A Bayesian/maximum-entropy view to the spatial estimation problem , 1990 .

[14]  George Christakos,et al.  Some Applications of the Bayesian, Maximum-Entropy Concept in Geostatistics , 1991 .

[15]  Pierre M. Adler,et al.  Porous media : geometry and transports , 1992 .

[16]  George Christakos,et al.  Random Field Models in Earth Sciences , 1992 .

[17]  R. Knox Subsurface Transport and Fate Processes , 1993 .

[18]  Kenneth H. Reckhow,et al.  Bayesian Model for Fate and Transport of Polychlorinated Biphenyl in Upper Hudson River , 1995 .

[19]  Francis A. DiGiano,et al.  Process Dynamics in Environmental Systems , 1996 .

[20]  Karlheinz Spitz,et al.  A Practical Guide to Groundwater and Solute Transport Modeling , 1996 .

[21]  William J. Mitsch,et al.  Environmental Modeling: Fate and Transport of Pollutants in Water, Air, and Soil , 1997 .

[22]  George Christakos,et al.  Spatiotemporal information systems in soil and environmental sciences , 1998 .

[23]  Dionissios T. Hristopulos,et al.  Spatiotemporal Environmental Health Modelling: A Tractatus Stochasticus , 1998, Springer US.

[24]  Arthur E. Mynett Hydroinformatics and its applications at Delft Hydraulics , 1999 .

[25]  G. Christakos,et al.  BME STUDIES OF STOCHASTIC DIFFERENTIAL EQUATIONS REPRESENTING PHYSICAL LAWS -PART I , 1999 .

[26]  Matthew W. Farthing,et al.  A comparison of high-resolution, finite-volume, adaptive-stencil schemes for simulating advective-dispersive transport , 2000 .

[27]  Dionissios T. Hristopulos,et al.  Practical Calculation of Non-Gaussian Multivariate Moments in Spatiotemporal Bayesian Maximum Entropy Analysis , 2001 .

[28]  You‐Kuan Zhang Stochastic Methods for Flow in Porous Media: Coping with Uncertainties , 2001 .

[29]  Arnold Neumaier,et al.  Introduction to Numerical Analysis , 2001 .