检索结果-内蒙古大学图书馆

Multiple-relaxation-time lattice Boltzmann model-based four-level finite-difference scheme for one-dimensional diffusion equations

学校读者我要写书评

暂无评论

Physical Review E 2021年第1期104卷 015312-015312页

作者： Yuxin Lin Ning Hong Baochang Shi Zhenhua Chai School of Mathematics and Statistics Huazhong University of Science and Technology Wuhan 430074 China School of General Education Wuchang University of Technology Wuhan 430223 China Hubei Key Laboratory of Engineering Modeling and Scientific Computing Huazhong University of Science and Technology Wuhan 430074 China

In this paper, we first present a multiple-relaxation-time lattice Boltzmann (MRT-LB) model for one-dimensional diffusion equation where the D1Q3 (three discrete velocities in one-dimensional space) lattice structure is considered. Then through the theoretical analysis, we derive an explicit four-level finite-difference scheme from this MRT-LB model. The results show that the four-level finite-difference scheme is unconditionally stable, and through adjusting the weight coefficient ω0 and the relaxation parameters s1 and s2 corresponding to the first and second moments, it can also have a sixth-order accuracy in space. Finally, we also test the four-level finite-difference scheme through some numerical simulations and find that the numerical results are consistent with our theoretical analysis.

关键词： Lattice-Boltzmann methods

A lattice boltzmann model for the coupled cross-diffusion-fluid system

学校读者我要写书评

暂无评论

arXiv 2021年

作者： Zhan, Chengjie Chai, Zhenhua Shi, Baochang School of Mathematics and Statistics Huazhong University of Science and Technology Wuhan430074 China Hubei Key Laboratory of Engineering Modeling and Scientific Computing Huazhong University of Science and Technology Wuhan430074 China

In this paper, we propose a lattice Boltzmann (LB) model for the generalized coupled cross-diffusion-fluid system. Through the direct Taylor expansion method, the proposed LB model can correctly recover the macroscopic equations. The cross diffusion terms in the coupled system are modeled by introducing additional collision operators, which can be used to avoid special treatments for the gradient terms. In addition, the auxiliary source terms are constructed properly such that the numerical diffusion caused by the convection can be eliminated. We adopt the developed LB model to study two important systems, i.e., the coupled chemotaxis-fluid system and the double-diffusive convection system with Soret and Dufour effects. We first test the present LB model through considering a steady-state case of coupled chemotaxis-fluid system, then we analyze the influences of some physical parameters on the formation of sinking plumes. Finally, the double-diffusive natural convection system with Soret and Dufour effects is also studied, and the numerical results agree well with some previous works. © 2021, CC BY.

关键词： Biochemistry

Multiple-distribution-function lattice Boltzmann method for convection-diffusion-system based incompressible Navier-Stokes equations

学校读者我要写书评

暂无评论

arXiv 2021年

作者： Chai, Zhenhua Shi, Baochang Zhan, Chengjie School of Mathematics and Statistics Huazhong University of Science and Technology Wuhan430074 China Hubei Key Laboratory of Engineering Modeling and Scientific Computing Huazhong University of Science and Technology Wuhan430074 China

In this paper, a multiple-distribution-function lattice Boltzmann method (MDF-LBM) with multiple-relaxation-time model is proposed for incompressible Navier-Stokes equations (NSEs) which are considered as the coupled convection-diffusion equations (CDEs). Through direct Taylor expansion analysis, we show that the Navier-Stokes equations can be recovered correctly from the present MDF-LBM, and additionally, it is also found that the velocity and pressure can be directly computed through the zero and first-order moments of distribution function. Then in the framework of present MDF-LBM, we develop a locally computational scheme for the velocity gradient where the first-order moment of the non-equilibrium distribution is used, this scheme is also extended to calculate the velocity divergence, strain rate tensor, shear stress and vorticity. Finally, we also conduct some simulations to test the MDF-LBM, and find that the numerical results not only agree with some available analytical and numerical solutions, but also have a second-order convergence rate in space. © 2021, CC BY.

关键词： Kinetic theory

The optimal displacement of immiscible two-phase fluids in a pore doublet

学校读者我要写书评

暂无评论

arXiv 2022年

作者： Shan, Fang Chai, Zhenhua Shi, Baochang Zhao, Meng School of Mathematics and Statistics Huazhong University of Science and Technology Wuhan430074 China Institute of Interdisciplinary Research for Mathematics and Applied Science Huazhong University of Science and Technology Wuhan430074 China Hubei Key Laboratory of Engineering Modeling and Scientific Computing Huazhong University of Science and Technology Wuhan430074 China

The displacement of multiphase fluid flow in a pore doublet is a fundamental problem, and is also of importance in understanding of the transport mechanisms of multiphase flows in the porous media. During the displacement of immiscible two-phase fluids in the pore doublet, the transport process is not only influenced by the capillary and viscous forces, but also affected by the channel geometry. In this paper, we first present a mathematical model to describe the two-phase fluid displacement in the pore doublet where the effects of capillary force, viscous force and the geometric structure are included. Then we derive an analytical solution of the model for the first time, and find that the displacement process is dominated by the capillary number, the viscosity ratio and the radius ratio. Furthermore, we define the optimal displacement that the wetting fluids in two daughter channels break through the branches simultaneously (both of them have the same breakthrough time), and also obtain the critical capillary number corresponding to the optimal displacement, which is related to the radius ratio of two daughter channels and viscosity ratio of two immiscible fluids. Finally, it is worthy noting that the present analytical results on the displacement in the pore doublet can be used to explain and understand the phenomenon of preferential imbibition or preferential flow in porous media. © 2022, CC BY.

关键词： Porous materials

GAUSSIAN QUADRATURE RULES FOR COMPOSITE HIGHLY OSCILLATORY INTEGRALS

学校读者我要写书评

暂无评论

arXiv 2021年

作者： Wu, Menghan Wang, Haiyong School of Mathematics and Statistics Huazhong University of Science and Technology Wuhan430074 China Hubei Key Laboratory od Engineering Modeling and Scientific Computing Huazhong University of Science and Technology Wuhan430074 China

Highly oscillatory integrals of composite type arise in electronic engineering and their calculations are a challenging problem. In this paper, we propose two Gaussian quadrature rules for computing such integrals. The first one is constructed based on the classical theory of orthogonal polynomials and its nodes and weights can be computed efficiently by using tools of numerical linear algebra. An interesting connection between the quadrature nodes and the Legendre points is proved and it is shown that the rate of convergence of this rule depends solely on the regularity of the non-oscillatory part of the integrand. The second one is constructed with respect to a sign-changing function and the classical theory of Gaussian quadrature can not be used anymore. We explore theoretical properties of this Gaussian quadrature, including the trajectories of the quadrature nodes and the convergence rate of these nodes to the endpoints of the integration interval, and prove its asymptotic error estimate under suitable hypotheses. Numerical experiments are presented to demonstrate the performance of the proposed *** Codes 65D30 Copyright © 2021, The Authors. All rights reserved.

关键词： Computation theory

Rectangular multiple-relaxation-time lattice Boltzmann method for the Navier-Stokes and nonlinear convection-diffusion equations: General equilibrium and some important issues

学校读者我要写书评

暂无评论

arXiv 2022年

作者： Chai, Zhenhua Yuan, Xiaolei Shi, Baochang Institute of Interdisciplinary Research for Mathematics and Applied Science School of Mathematics and Statistics Huazhong University of Science and Technology Wuhan430074 China Hubei Key Laboratory of Engineering Modeling and Scientific Computing Huazhong University of Science and Technology Wuhan430074 China College of Mathematics and Information Science Hebei University Baoding071002 China

In this paper, we develop a general rectangular multiple-relaxation-time lattice Boltzmann (RMRT-LB) method for the Navier-Stokes equations (NSEs) and nonlinear convection-diffusion equation (NCDE) by extending our recent unified framework of MRT-LB method [Z.H. Chai and B.C. Shi, Phys. Rev. E 102, 023306 (2020)], where a rectangular equilibrium distribution function (REDF) [J.H. Lu et al, Phil. Trans. R. Soc. A 369, 2311-2319 (2011)] on a rectangular lattice is utilized. Due to the anisotropy of discrete velocities on a rectangular lattice, the third-order moment of REDF is inconsistent with that of popular LB method, and thus the previous unified framework of MRT-LB method cannot be directly applied to the NSEs using the REDF on the rectangular rDdQq (q discrete velocities in d-dimensional space, d ≥ 1) lattice. The macroscopic NSEs can be recovered from the RMRT-LB method through the direct Taylor expansion method by properly selecting the relaxation sub-matrix S2 which is related to kinetic viscosity and bulk viscosity. While the rectangular lattice does not lead to the change of the zero-th, first and second-order moments of REDF, thus the unified framework of MRT-LB method can be directly applied to the NCDE. It should be noted that the RMRT-LB model for NSEs can be derived on the rDdQq lattice, including rD2Q9, rD3Q19, and rD3Q27 lattices, while there are no rectangular D3Q13 and D3Q15 lattices within this framework of RMRT-LB method. Thanks to the block-lower triangular-relaxation matrix introduced in the unified framework, the RMRT-LB versions (if existed) of the previous MRT-LB models can be obtained, including those based on raw (natural)-moment, central-moment, Hermite-moment, and central Hermite-moment, respectively. It is also found that when the standard or regular lattice is used, and the sound speed is taken as an adjustable parameter, the present RMRT-LB method becomes a new class of MRT-LB method for the NSEs and NCDE, and the commonly used MRT-LB

关键词： Navier Stokes equations

GLOBAL SOLVABILITY IN A TWO-SPECIES CHEMOTAXIS SYSTEM WITH SIGNAL PRODUCTION

学校读者我要写书评

暂无评论

arXiv 2021年

作者： Ren, Guoqiang Xiang, Tian School of Mathematics and Statistics Huazhong University of Science and Technology Hubei Wuhan430074 China Hubei Key Laboratory of Engineering Modeling and Scientific Computing Huazhong University of Science and Technology Hubei Wuhan430074 China Institute for Mathematical Sciences Renmin University of China Bejing100872 China

In this paper, we study the following prototypical two-species chemotaxis system with Lotka-Volterra competition and signal production: (Equation presented). We show that if (Equation presented). the associated Neumann initial-boundary value problem for (∗) admits a global generalized solutions in a bounded and smooth domain Ω ⊂ N (N ≥ 2) under merely integrable initial *** Codes 35A01, 35K57, 35Q92, 92C17 © 2021, CC BY.

关键词： Initial value problems

GLOBAL SOLVABILITY AND ASYMPTOTICAL BEHAVIOR IN A TWO-SPECIES CHEMOTAXIS MODEL WITH SIGNAL ABSORPTION

学校读者我要写书评

暂无评论

arXiv 2021年

In this work, we study global existence, eventual smoothness and asymptotical behavior of positive solutions for the following two-species chemotaxis consumption model: ut = ∆u - χ1∇ · (u∇w), x ∈ Ω, t > 0, {{{{{{{ vt = ∆v - χ2∇ · (u∇w), x ∈ Ω, t > 0, wt = ∆w - (αu + βv)w, x ∈ Ω, t > 0, in a bounded smooth but not necessarily convex domain Ω ⊂ Rn(n = 2, 3, 4, 5) with nonnegative initial data u0, v0, w0 and homogeneous Neumann boundary data. Here, the parameters χ1, χ2 are positive and α, β are nonnegative. Under a smallness condition max{χ1, χ2}kw0kL∞ MSC Codes 35K59, 35K57 © 2021, CC BY.

关键词： Biochemistry

DENSITY CONVERGENCE OF A FULLY DISCRETE FINITE DIFFERENCE METHOD FOR STOCHASTIC CAHN-HILLIARD EQUATION

学校读者我要写书评

暂无评论

arXiv 2022年

作者： Hong, Jialin Jin, Diancong Sheng, Derui LSEC ICMSEC Academy of Mathematics and Systems Science Chinese Academy of Sciences School of Mathematical Sciences University of Chinese Academy of Sciences Beijing100049 China School of Mathematics and Statistics Huazhong University of Science and Technology Hubei Key Laboratory of Engineering Modeling and Scientific Computing Huazhong University of Science and Technology Wuhan430074 China

This paper focuses on investigating the density convergence of a fully discrete finite difference method when applied to numerically solve the stochastic Cahn-Hilliard equation driven by multiplicative space-time white noises. The main difficulty lies in the control of the drift coefficient that is neither globally Lipschitz nor one-sided Lipschitz. To handle this difficulty, we propose a novel localization argument and derive the strong convergence rate of the numerical solution to estimate the total variation distance between the exact and numerical solutions. This along with the existence of the density of the numerical solution finally yields the convergence of density in L1() of the numerical solution. Our results partially answer positively to the open problem emerged in [J. Cui and J. Hong, J. Differential Equations (2020)] on computing the density of the exact solution numerically. Copyright © 2022, The Authors. All rights reserved.

关键词： Finite difference method