A structure-preserving scheme for computing effective diffusivity and anomalous diffusion phenomena of random flows

Tan Zhang,Zhongjian Wang,Jack Xin,Zhiwen Zhang
2024-05-29
Abstract:This paper aims to investigate the diffusion behavior of particles moving in stochastic flows under a structure-preserving scheme. We compute the effective diffusivity for normal diffusive random flows and establish the power law between spatial and temporal variables for cases with anomalous diffusion phenomena. From a Lagrangian approach, we separate the corresponding stochastic differential equations (SDEs) into sub-problems and construct a one-step structure-preserving method to solve them. Then by modified equation systems, the convergence analysis in calculating the effective diffusivity is provided and compared between the structure-preserving scheme and the Euler-Maruyama scheme. Also, we provide the error estimate for the structure-preserving scheme in calculating the power law for a series of super-diffusive random flows. Finally, we calculate the effective diffusivity and anomalous diffusion phenomena for a series of 2D and 3D random fields.
Numerical Analysis
What problem does this paper attempt to address?
### What problems does this paper attempt to solve? This paper aims to study the diffusion behavior of particles in random flows, especially the calculation of the effective diffusion rate and the anomalous diffusion phenomenon under the structure - preserving scheme. Specifically, the article focuses on the following core issues: 1. **Calculation of the effective diffusion rate**: - For the random flow of normal diffusion, the author proposes a structure - preserving method to calculate the effective diffusion rate \( D_E \). - Specifically, the corresponding stochastic differential equations (SDEs) are decomposed into sub - problems by the Lagrangian method, and a one - step structure - preserving method is constructed to solve these sub - problems. 2. **Research on the anomalous diffusion phenomenon**: - For the random flow with the anomalous diffusion phenomenon (such as super - diffusion and sub - diffusion), the author establishes a power - law relationship between the spatial variable and the time variable. - The anomalous diffusion phenomena in different types of random flows are studied, and the validity and accuracy of the calculation by the structure - preserving method in this case are verified. 3. **Comparison and analysis of numerical methods**: - The convergence analysis of the structure - preserving method in calculating the effective diffusion rate is provided, and it is compared with the traditional Euler - Maruyama method. - For a series of super - diffusive random flows, the error estimate of the power - law calculated by the structure - preserving method is provided. 4. **Treatment of high - dimensional random flows**: - The effective diffusion rate and the anomalous diffusion phenomenon in two - dimensional and three - dimensional random fields are studied. - Different numerical scheme construction methods for two - dimensional and three - dimensional random flows are proposed to ensure that the calculation can be effectively carried out in high - dimensional cases as well. ### Summary of mathematical formulas - **Stochastic differential equations (SDEs)**: \[ dx(t)=v(x, t, \omega)dt+\sigma dw(t),\quad x\in\mathbb{R}^d,\quad t\geq0 \] where \( x(t)=(x_1(t),\ldots,x_d(t))^T\in\mathbb{R}^d \) is the spatial position of the passive tracer particle, \( v(x, t, \omega) \) is the incompressible velocity field, \( \sigma>0 \) is the molecular diffusion rate, and \( w(t) \) is the standard \( d \)-dimensional Brownian motion. - **Effective diffusion rate matrix \( D_E \)**: \[ e^T D_E e = D_0+(-L_0 G_e, G_e)_{L^2(\mathbb{R}^d)} \] where \( D_0=\frac{\sigma^2}{2} \), \( L_0 f=\frac{\sigma^2}{2}\Delta f + v\cdot\nabla f \) is the generator, and \( G \) satisfies the continuous correction problem \( L_0 G=-v \). - **Volume - preserving property**: \[ \det\left(\frac{\partial x_n^*}{\partial x_n}\right)=1 \] Through these studies, the author not only solves the key challenge of accurately characterizing and quantifying the large - scale effective diffusion in complex fluid flows, but also provides theoretical support for applications in fields such as atmospheric science, oceanography, chemical engineering, and combustion.