Algorithm refinement for stochastic partial differential equations.

A hybrid particle/continuum algorithm is formulated for Fickian diffusion in the fluctuating hydrodynamic limit. The particles are taken as independent random walkers; the fluctuating diffusion equation is solved by finite differences with deterministic and white-noise fluxes. At the interface between the particle and continuum computations the coupling is by flux matching, giving exact mass conservation. This methodology is an extension of Adaptive Mesh and Algorithm Refinement to stochastic partial differential equations. A variety of numerical experiments were performed for both steady and time-dependent scenarios. In all cases the mean and variance of density are captured correctly by the stochastic hybrid algorithm. For a non-stochastic version (i.e., using only deterministic continuum fluxes) the mean density is correct, but the variance is reduced except within the particle region, far from the interface. Extensions of the methodology to fluid mechanics applications are discussed.

[1]  Robert E. Rudd,et al.  Concurrent Coupling of Length Scales in Solid State Systems , 2000 .

[2]  I. Oppenheim,et al.  Spatial correlations in bounded nonequilibrium fluid systems , 1997 .

[3]  Moulay D. Tidriri,et al.  Coupling Boltzmann and Navier-Stokes Equations by Friction , 1996 .

[4]  A. Patera,et al.  Heterogeneous Atomistic-Continuum Representations for Dense Fluid Systems , 1997 .

[5]  Alejandro L. Garcia,et al.  A simple model for nonequilibrium fluctuations in a fluid , 1996 .

[6]  J. Dufty.  Long‐Range Correlations in Nonequilibrium Fluids a , 1998 .

[7]  Noam Bernstein,et al.  Spanning the length scales in dynamic simulation , 1998 .

[8]  N. Kampen,et al.  Stochastic processes in physics and chemistry , 1981 .

[9]  M. Ortiz,et al.  An adaptive finite element approach to atomic-scale mechanics—the quasicontinuum method , 1997, cond-mat/9710027.

[10]  Axel Klar,et al.  An adaptive domain decomposition procedure for Boltzmann and Euler equations , 1998 .

[11]  Hassan Hassan,et al.  A decoupled DSMC/Navier-Stokes analysis of a transitional flow experiment , 1996 .

[12]  W. Ebeling.  Stochastic Processes in Physics and Chemistry , 1995 .

[13]  Alejandro L. Garcia,et al.  Numerical integration of the fluctuating hydrodynamic equations , 1987 .

[14]  D. C. Wadsworth,et al.  One-dimensional hybrid continuum/particle simulation approach for rarefied hypersonic flows , 1990 .

[15]  Alejandro L. Garcia.  Numerical methods for physics , 1994 .

[16]  C. Gardiner.  Handbook of Stochastic Methods , 1983 .

[17]  R Schmitz,et al.  FLUCTUATIONS IN NONEQUILIBRIUM FLUIDS , 1988 .

[18]  Eirik Grude Flekkøy,et al.  Hybrid model for combined particle and continuum dynamics , 2000 .

[19]  Berni J. Alder.  Highly discretized dynamics , 1997 .

[20]  Patrick Le Tallec,et al.  Coupling Boltzmann and Navier-Stokes Equations by Half Fluxes , 1997 .

[21]  Alejandro L. Garcia,et al.  Adaptive Mesh and Algorithm Refinement Using Direct Simulation Monte Carlo , 1999 .

[22]  J Feder,et al.  Coupling particles and fields in a diffusive hybrid model. , 2001, Physical review. E, Statistical, nonlinear, and soft matter physics.