The transport equation for the Nambu-Jona-Lasinio (NJL) model isderived phenomenologically. The finite-temperature effective mass for the quarkis analysed, By means of the Chapman-Enskog method and hydrodynamic ap-pro...The transport equation for the Nambu-Jona-Lasinio (NJL) model isderived phenomenologically. The finite-temperature effective mass for the quarkis analysed, By means of the Chapman-Enskog method and hydrodynamic ap-proach the different transport coefficients for the NJL plasma are calculated tothe first order in the relaxation time.展开更多
The resolution by Chen and Sun of divergent Chapman-Enskog expansion problem will not only build a unified foundation for non-equilibrium dynamics modeling at all Mach number and Knudsen number, but also shed light to...The resolution by Chen and Sun of divergent Chapman-Enskog expansion problem will not only build a unified foundation for non-equilibrium dynamics modeling at all Mach number and Knudsen number, but also shed light to a large class of difficult theoretical problems involving divergent expansion on strong nonlinearity.展开更多
This work develops asymptotic expansions of systems of partial differential equations associated with multi-scale switching diffusions. The switching process is modeled by using an inhomogeneous continuous- time Marko...This work develops asymptotic expansions of systems of partial differential equations associated with multi-scale switching diffusions. The switching process is modeled by using an inhomogeneous continuous- time Markov chain. In the model, there are two small parameters ε and δ. The first one highlights the fast switching, whereas the other delineates the slow diffusion. Assuming that ε and δ are related in that ε = δγ, our results reveal that different values of γ lead to different behaviors of the underlying systems, resulting in different asymptotic expansions. Although our motivation comes from stochastic problems, the approach is mainly analytic and is constructive. The asymptotic series are rigorously justified with error bounds provided. An example is provided to demonstrate the results.展开更多
In [16] a visco-elastic relaxation system, called the relaxed Burnett system, was proposed by Jin and Slemrod as a moment approximation to the Boltzmann equation. The relaxed Burnett system is weakly parabolic, has a ...In [16] a visco-elastic relaxation system, called the relaxed Burnett system, was proposed by Jin and Slemrod as a moment approximation to the Boltzmann equation. The relaxed Burnett system is weakly parabolic, has a linearly hyperbolic convection part, and is endowed with a generalized entropy inequality. It agrees with the solution of the Boltzmann equation up to the Burnett order via the Chapman-Enskog expansion. We develop a one-dimensional non-oscillatory numerical scheme based on the relaxed Burnett system for the Boltzmann equation. We compare numerical results for stationary shocks based on this relaxation scheme, and those obtained by the DSMC (Direct Simulation Monte Carlo), by the Navier-Stokes equations and by the extended thermodynamics with thirteen moments (the Grad equations). Our numerical experiments show that the relaxed Burnett gives more accurate approximations to the shock profiles of the Boltzmann equation obtained by the DSMC, for a range of Mach numbers for hypersonic flows, than those obtained by the other hydrodynamic systems.展开更多
This paper proposes a lattice Boltzmann model with an amending function for one-dimensional nonlinear partial differential equations (NPDEs) in the form ut +αuux +βu^nuz +γuxx +δuzxx +ζxxxx = 0. This model...This paper proposes a lattice Boltzmann model with an amending function for one-dimensional nonlinear partial differential equations (NPDEs) in the form ut +αuux +βu^nuz +γuxx +δuzxx +ζxxxx = 0. This model is different from existing models because it lets the time step be equivalent to the square of the space step and derives higher accuracy and nonlinear terms in NPDEs. With the Chapman-Enskog expansion, the governing evolution equation is recovered correctly from the continuous Boltzmann equation. The numerical results agree well with the analytical solutions.展开更多
In this paper, in order to extend the lattice Boltzmann method to deal with more nonlinear equations, a one-dimensional (1D) lattice Boltzmann scheme with an amending function for the nonlinear Klein-Gordon equation i...In this paper, in order to extend the lattice Boltzmann method to deal with more nonlinear equations, a one-dimensional (1D) lattice Boltzmann scheme with an amending function for the nonlinear Klein-Gordon equation is proposed. With the Taylor and Chapman-Enskog expansion, the nonlinear Klein-Gordon equation is recovered correctly from the lattice Boltzmann equation. The method is applied on some test examples, and the numerical results have been compared with the analytical solutions or the numerical solutions reported in previous studies. The L2, L∞ and Root-Mean-Square (RMS) errors in the solutions show the efficiency of the method computationally.展开更多
By combining Chapman-Enskog expansion with the BGK approximation to Baltzmann equation and Navier-Stokes equation was obtained. And an expression of (Darcy's) law was obtained through taking variable average over ...By combining Chapman-Enskog expansion with the BGK approximation to Baltzmann equation and Navier-Stokes equation was obtained. And an expression of (Darcy's) law was obtained through taking variable average over Navier-Stokes equation on some representative space in porous media,and finally an example was taken to prove its reliability.展开更多
In this paper, we consider a one-dimensional nonlinear partial differential equation that has the form ut + αuux + βunux - γuxx + δuxxx = F(u). A higher order lattice Bhatnager-Gross-Krook (BGK) model with an amen...In this paper, we consider a one-dimensional nonlinear partial differential equation that has the form ut + αuux + βunux - γuxx + δuxxx = F(u). A higher order lattice Bhatnager-Gross-Krook (BGK) model with an amending-function is proposed. With the Chapman-Enskog expansion, different kinds of nonlinear partial differential equations are recovered correctly from the continuous Boltzmann equation. The numerical results show that this method is very effective.展开更多
In this paper,we develop a lattice Boltzmann model for a class ofone-dimensional nonlinear wave equations,including the second-order hyperbolictelegraph equation,the nonlinear Klein-Gordon equation,the damped and unda...In this paper,we develop a lattice Boltzmann model for a class ofone-dimensional nonlinear wave equations,including the second-order hyperbolictelegraph equation,the nonlinear Klein-Gordon equation,the damped and undampedsine-Gordon equation and double sine-Gordon equation.By choosing properly theconservation condition between the macroscopic quantity u,and the distributionfunctions and applying the Chapman-Enskog expansion,the governing equation isrecovered correctly from the lattice Boltzmann equation.Moreover,the local equilib-rium distribution function is obtained.The results of numerical examples have beencompared with the analytical solutions to confirm the good accuracy and the applica-bility of our scheme.展开更多
We develop a two-relaxation-time (TRT) Lattice Boltzmann model for hydrodynamicequations with variable source terms based on equivalent equilibriumfunctions. A special parametrization of the free relaxation parameter ...We develop a two-relaxation-time (TRT) Lattice Boltzmann model for hydrodynamicequations with variable source terms based on equivalent equilibriumfunctions. A special parametrization of the free relaxation parameter is derived. Itcontrols, in addition to the non-dimensional hydrodynamic numbers, any TRT macroscopicsteady solution and governs the spatial discretization of transient flows. Inthis framework, the multi-reflection approach [16, 18] is generalized and extended forDirichlet velocity, pressure and mixed (pressure/tangential velocity) boundary conditions.We propose second and third-order accurate boundary schemes and adapt themfor corners. The boundary schemes are analyzed for exactness of the parametrization,uniqueness of their steady solutions, support of staggered invariants and for the effectiveaccuracy in case of time dependent boundary conditions and transient flow.When the boundary scheme obeys the parametrization properly, the derived permeabilityvalues become independent of the selected viscosity for any porous structureand can be computed efficiently. The linear interpolations [5, 46] are improved withrespect to this property.展开更多
In this paper,the authors study the 1 D steady Boltzmann flow in a channel.The walls of the channel are assumed to have vanishing velocity and given temperaturesθ0andθ1.This problem was studied by Esposito-Lebowitz-...In this paper,the authors study the 1 D steady Boltzmann flow in a channel.The walls of the channel are assumed to have vanishing velocity and given temperaturesθ0andθ1.This problem was studied by Esposito-Lebowitz-Marra(1994,1995)where they showed that the solution tends to a local Maxwellian with parameters satisfying the compressible Navier-Stokes equation with no-slip boundary condition.However,a lot of numerical experiments reveal that the fluid layer does not entirely stick to the boundary.In the regime where the Knudsen number is reasonably small,the slip phenomenon is significant near the boundary.Thus,they revisit this problem by taking into account the slip boundary conditions.Following the lines of[Coron,F.,Derivation of slip boundary conditions for the Navier-Stokes system from the Boltzmann equation,J.Stat.Phys.,54(3-4),1989,829-857],the authors will first give a formal asymptotic analysis to see that the flow governed by the Boltzmann equation is accurately approximated by a superposition of a steady CNS equation with a temperature jump condition and two Knudsen layers located at end points.Then they will establish a uniform L∞estimate on the remainder and derive the slip boundary condition for compressible Navier-Stokes equations rigorously.展开更多
In this paper,a simplified lattice Boltzmann method(SLBM)without evolution of the distribution function is developed for simulating incompressible viscous flows.This method is developed from the application of fractio...In this paper,a simplified lattice Boltzmann method(SLBM)without evolution of the distribution function is developed for simulating incompressible viscous flows.This method is developed from the application of fractional step technique to the macroscopic Navier-Stokes(N-S)equations recovered from lattice Boltzmann equation by using Chapman-Enskog expansion analysis.In SLBM,the equilibrium distribution function is calculated from the macroscopic variables,while the non-equilibrium distribution function is simply evaluated from the difference of two equilibrium distribution functions.Therefore,SLBM tracks the evolution of the macroscopic variables rather than the distribution function.As a result,lower virtual memories are required and physical boundary conditions could be directly implemented.Through numerical test at high Reynolds number,the method shows very nice performance in numerical stability.An accuracy test for the 2D Taylor-Green flow shows that SLBM has the second-order of accuracy in space.More benchmark tests,including the Couette flow,the Poiseuille flow as well as the 2D lid-driven cavity flow,are conducted to further validate the present method;and the simulation results are in good agreement with available data in literatures.展开更多
In this paper,a gas kinetic scheme for the compressible multicomponent flows is presented by making use of two-species BGK model in[A.D.Kotelnikov and D.C.Montgomery,A Kinetic Method for Computing Inhomogeneous Fluid ...In this paper,a gas kinetic scheme for the compressible multicomponent flows is presented by making use of two-species BGK model in[A.D.Kotelnikov and D.C.Montgomery,A Kinetic Method for Computing Inhomogeneous Fluid Behavior,J.Comput.Phys.134(1997)364-388].Different from the conventional BGK model,the collisions between different species are taken into consideration.Based on the Chapman-Enskog expansion,the corresponding macroscopic equations are derived from this two-species model.Because of the relaxation terms in the governing equations,the method of operator splitting is applied.In the hyperbolic part,the integral solutions of the BGK equations are used to construct the numerical fluxes at the cell interface in the framework of finite volume method.Numerical tests are presented in this paper to validate the current approach for the compressible multicomponent flows.The theoretical analysis on the spurious oscillations at the interface is also presented.展开更多
Concepts of the lattice Boltzmann method are discussed in detail for the one-dimensional kinetic model.Various techniques of constructing lattice Boltzmann models are discussed,and novel collision integrals are derive...Concepts of the lattice Boltzmann method are discussed in detail for the one-dimensional kinetic model.Various techniques of constructing lattice Boltzmann models are discussed,and novel collision integrals are derived.Geometry of the ki-netic space and the role of the thermodynamic projector is elucidated.展开更多
A lattice Boltzmann method is developed for modeling viscous elementary flows.An adjustable source term is added to the lattice Boltzmann equation,which can be tuned to model different elementary flow features like a ...A lattice Boltzmann method is developed for modeling viscous elementary flows.An adjustable source term is added to the lattice Boltzmann equation,which can be tuned to model different elementary flow features like a doublet or a point source of any strength,including a negative source(sink).The added source term is dimensionally consistent with the lattice Boltzmann equation.The proposed model has many practical applications,as it can be used in the framework of the potential flow theory of viscous and viscoelastic fluids.The model can be easily extended to the three dimensional case.The model is verified by comparing its results with the analytical solution for some benchmark problems.The results are in good agreement with the analytical solution of the potential flow theory.展开更多
基金The project supported by the Postdoctor Science Foundation the Nuclear Industry Foundation of China
文摘The transport equation for the Nambu-Jona-Lasinio (NJL) model isderived phenomenologically. The finite-temperature effective mass for the quarkis analysed, By means of the Chapman-Enskog method and hydrodynamic ap-proach the different transport coefficients for the NJL plasma are calculated tothe first order in the relaxation time.
文摘The resolution by Chen and Sun of divergent Chapman-Enskog expansion problem will not only build a unified foundation for non-equilibrium dynamics modeling at all Mach number and Knudsen number, but also shed light to a large class of difficult theoretical problems involving divergent expansion on strong nonlinearity.
基金supported in part by the Air Force Office of Scientific Research under FA9550-15-1-0131
文摘This work develops asymptotic expansions of systems of partial differential equations associated with multi-scale switching diffusions. The switching process is modeled by using an inhomogeneous continuous- time Markov chain. In the model, there are two small parameters ε and δ. The first one highlights the fast switching, whereas the other delineates the slow diffusion. Assuming that ε and δ are related in that ε = δγ, our results reveal that different values of γ lead to different behaviors of the underlying systems, resulting in different asymptotic expansions. Although our motivation comes from stochastic problems, the approach is mainly analytic and is constructive. The asymptotic series are rigorously justified with error bounds provided. An example is provided to demonstrate the results.
基金Supported by NSF grant DMS-0196106 Supported by NSF grant DMS-9803223 and DMS-00711463.
文摘In [16] a visco-elastic relaxation system, called the relaxed Burnett system, was proposed by Jin and Slemrod as a moment approximation to the Boltzmann equation. The relaxed Burnett system is weakly parabolic, has a linearly hyperbolic convection part, and is endowed with a generalized entropy inequality. It agrees with the solution of the Boltzmann equation up to the Burnett order via the Chapman-Enskog expansion. We develop a one-dimensional non-oscillatory numerical scheme based on the relaxed Burnett system for the Boltzmann equation. We compare numerical results for stationary shocks based on this relaxation scheme, and those obtained by the DSMC (Direct Simulation Monte Carlo), by the Navier-Stokes equations and by the extended thermodynamics with thirteen moments (the Grad equations). Our numerical experiments show that the relaxed Burnett gives more accurate approximations to the shock profiles of the Boltzmann equation obtained by the DSMC, for a range of Mach numbers for hypersonic flows, than those obtained by the other hydrodynamic systems.
基金Project supported by the National Natural Science Foundation of China (Grant No. 10661005)Fujian Province Science and Technology Plan Item (Grant No. 2008F5019)
文摘This paper proposes a lattice Boltzmann model with an amending function for one-dimensional nonlinear partial differential equations (NPDEs) in the form ut +αuux +βu^nuz +γuxx +δuzxx +ζxxxx = 0. This model is different from existing models because it lets the time step be equivalent to the square of the space step and derives higher accuracy and nonlinear terms in NPDEs. With the Chapman-Enskog expansion, the governing evolution equation is recovered correctly from the continuous Boltzmann equation. The numerical results agree well with the analytical solutions.
文摘In this paper, in order to extend the lattice Boltzmann method to deal with more nonlinear equations, a one-dimensional (1D) lattice Boltzmann scheme with an amending function for the nonlinear Klein-Gordon equation is proposed. With the Taylor and Chapman-Enskog expansion, the nonlinear Klein-Gordon equation is recovered correctly from the lattice Boltzmann equation. The method is applied on some test examples, and the numerical results have been compared with the analytical solutions or the numerical solutions reported in previous studies. The L2, L∞ and Root-Mean-Square (RMS) errors in the solutions show the efficiency of the method computationally.
文摘By combining Chapman-Enskog expansion with the BGK approximation to Baltzmann equation and Navier-Stokes equation was obtained. And an expression of (Darcy's) law was obtained through taking variable average over Navier-Stokes equation on some representative space in porous media,and finally an example was taken to prove its reliability.
基金Supported by the National Natural Science Foundation of China (Grant No 10661005) the Science and Technology Plan Item of Fujian Province (Grant No 2008F5019)
文摘In this paper, we consider a one-dimensional nonlinear partial differential equation that has the form ut + αuux + βunux - γuxx + δuxxx = F(u). A higher order lattice Bhatnager-Gross-Krook (BGK) model with an amending-function is proposed. With the Chapman-Enskog expansion, different kinds of nonlinear partial differential equations are recovered correctly from the continuous Boltzmann equation. The numerical results show that this method is very effective.
基金The authors are very thankful to the reviewers for their valuable suggestions toimprove the quality of the paper.This work is supported by National Natural Science Foundation of China(Nos.11101399,11271171,11301234)the Provincial Natural Science Foundation of Jiangxi(Nos.20161ACB20006,20142BCB23009,20151BAB201012).
文摘In this paper,we develop a lattice Boltzmann model for a class ofone-dimensional nonlinear wave equations,including the second-order hyperbolictelegraph equation,the nonlinear Klein-Gordon equation,the damped and undampedsine-Gordon equation and double sine-Gordon equation.By choosing properly theconservation condition between the macroscopic quantity u,and the distributionfunctions and applying the Chapman-Enskog expansion,the governing equation isrecovered correctly from the lattice Boltzmann equation.Moreover,the local equilib-rium distribution function is obtained.The results of numerical examples have beencompared with the analytical solutions to confirm the good accuracy and the applica-bility of our scheme.
文摘We develop a two-relaxation-time (TRT) Lattice Boltzmann model for hydrodynamicequations with variable source terms based on equivalent equilibriumfunctions. A special parametrization of the free relaxation parameter is derived. Itcontrols, in addition to the non-dimensional hydrodynamic numbers, any TRT macroscopicsteady solution and governs the spatial discretization of transient flows. Inthis framework, the multi-reflection approach [16, 18] is generalized and extended forDirichlet velocity, pressure and mixed (pressure/tangential velocity) boundary conditions.We propose second and third-order accurate boundary schemes and adapt themfor corners. The boundary schemes are analyzed for exactness of the parametrization,uniqueness of their steady solutions, support of staggered invariants and for the effectiveaccuracy in case of time dependent boundary conditions and transient flow.When the boundary scheme obeys the parametrization properly, the derived permeabilityvalues become independent of the selected viscosity for any porous structureand can be computed efficiently. The linear interpolations [5, 46] are improved withrespect to this property.
基金supported by the National Natural Science Foundation of China(Nos.11971201,11731008)the General Research Fund from RGC of Hong Kong(No.14301719)+1 种基金a Direct Grant from CUHK(No.4053397)the Fundamental Research Funds for the Central Universities and a fellowship award from the Research Grants Council of the Hong Kong Special Administrative Region,China(No.SRF2021-1S01)。
文摘In this paper,the authors study the 1 D steady Boltzmann flow in a channel.The walls of the channel are assumed to have vanishing velocity and given temperaturesθ0andθ1.This problem was studied by Esposito-Lebowitz-Marra(1994,1995)where they showed that the solution tends to a local Maxwellian with parameters satisfying the compressible Navier-Stokes equation with no-slip boundary condition.However,a lot of numerical experiments reveal that the fluid layer does not entirely stick to the boundary.In the regime where the Knudsen number is reasonably small,the slip phenomenon is significant near the boundary.Thus,they revisit this problem by taking into account the slip boundary conditions.Following the lines of[Coron,F.,Derivation of slip boundary conditions for the Navier-Stokes system from the Boltzmann equation,J.Stat.Phys.,54(3-4),1989,829-857],the authors will first give a formal asymptotic analysis to see that the flow governed by the Boltzmann equation is accurately approximated by a superposition of a steady CNS equation with a temperature jump condition and two Knudsen layers located at end points.Then they will establish a uniform L∞estimate on the remainder and derive the slip boundary condition for compressible Navier-Stokes equations rigorously.
文摘In this paper,a simplified lattice Boltzmann method(SLBM)without evolution of the distribution function is developed for simulating incompressible viscous flows.This method is developed from the application of fractional step technique to the macroscopic Navier-Stokes(N-S)equations recovered from lattice Boltzmann equation by using Chapman-Enskog expansion analysis.In SLBM,the equilibrium distribution function is calculated from the macroscopic variables,while the non-equilibrium distribution function is simply evaluated from the difference of two equilibrium distribution functions.Therefore,SLBM tracks the evolution of the macroscopic variables rather than the distribution function.As a result,lower virtual memories are required and physical boundary conditions could be directly implemented.Through numerical test at high Reynolds number,the method shows very nice performance in numerical stability.An accuracy test for the 2D Taylor-Green flow shows that SLBM has the second-order of accuracy in space.More benchmark tests,including the Couette flow,the Poiseuille flow as well as the 2D lid-driven cavity flow,are conducted to further validate the present method;and the simulation results are in good agreement with available data in literatures.
基金Natural Science Foundation of China(NSFC)No.10931004,No.11171037 and No.91130021.
文摘In this paper,a gas kinetic scheme for the compressible multicomponent flows is presented by making use of two-species BGK model in[A.D.Kotelnikov and D.C.Montgomery,A Kinetic Method for Computing Inhomogeneous Fluid Behavior,J.Comput.Phys.134(1997)364-388].Different from the conventional BGK model,the collisions between different species are taken into consideration.Based on the Chapman-Enskog expansion,the corresponding macroscopic equations are derived from this two-species model.Because of the relaxation terms in the governing equations,the method of operator splitting is applied.In the hyperbolic part,the integral solutions of the BGK equations are used to construct the numerical fluxes at the cell interface in the framework of finite volume method.Numerical tests are presented in this paper to validate the current approach for the compressible multicomponent flows.The theoretical analysis on the spurious oscillations at the interface is also presented.
基金support by the BFE Project 100862.S.S.C.was supported by the ETH Project 0-20280-05.
文摘Concepts of the lattice Boltzmann method are discussed in detail for the one-dimensional kinetic model.Various techniques of constructing lattice Boltzmann models are discussed,and novel collision integrals are derived.Geometry of the ki-netic space and the role of the thermodynamic projector is elucidated.
文摘A lattice Boltzmann method is developed for modeling viscous elementary flows.An adjustable source term is added to the lattice Boltzmann equation,which can be tuned to model different elementary flow features like a doublet or a point source of any strength,including a negative source(sink).The added source term is dimensionally consistent with the lattice Boltzmann equation.The proposed model has many practical applications,as it can be used in the framework of the potential flow theory of viscous and viscoelastic fluids.The model can be easily extended to the three dimensional case.The model is verified by comparing its results with the analytical solution for some benchmark problems.The results are in good agreement with the analytical solution of the potential flow theory.