This paper studies the strong convergence of the quantum lattice Boltzmann(QLB)scheme for the nonlinear Dirac equations for Gross-Neveu model in 1+1 dimensions.The initial data for the scheme are assumed to be converg...This paper studies the strong convergence of the quantum lattice Boltzmann(QLB)scheme for the nonlinear Dirac equations for Gross-Neveu model in 1+1 dimensions.The initial data for the scheme are assumed to be convergent in L^(2).Then for any T≥0 the corresponding solutions for the quantum lattice Boltzmann scheme are shown to be convergent in C([0,T];L^(2)(R^(1)))to the strong solution to the nonlinear Dirac equations as the mesh sizes converge to zero.In the proof,at first a Glimm type functional is introduced to establish the stability estimates for the difference between two solutions for the corresponding quantum lattice Boltzmann scheme,which leads to the compactness of the set of the solutions for the quantum lattice Boltzmann scheme.Finally the limit of any convergent subsequence of the solutions for the quantum lattice Boltzmann scheme is shown to coincide with the strong solution to a Cauchy problem for the nonlinear Dirac equations.展开更多
With the help of the method of separation of variables and the Debye-Hüchel approximation, the Poisson-Boltzmann equation that describes the distribution of the potential in the electrical double layer of a cylin...With the help of the method of separation of variables and the Debye-Hüchel approximation, the Poisson-Boltzmann equation that describes the distribution of the potential in the electrical double layer of a cylindrical particle with a limited length has been firstly solved under a very low potential condition. Then with the help of the functional analysis theory this equation has been further analytically solved under general potential conditions and consequently, the corresponding surface charge densities have been obtained. Both the potential and the surface charge densities cointide with those results obtained from the Debye-Hüchel approximation when the very low potential of zeψ〈〈kT is introduced.展开更多
A Lagrangian lattice Boltzmann method for solving Euler equations is proposed. The key step in formulating this method is the introduction of the displacement distribution function. The equilibrium distribution functi...A Lagrangian lattice Boltzmann method for solving Euler equations is proposed. The key step in formulating this method is the introduction of the displacement distribution function. The equilibrium distribution function consists of macroscopic Lagrangian variables at time steps n and n + 1. It is different from the standard lattice Boltzmann method. In this method the element, instead of each particle, is required to satisfy the basic law. The element is considered as one large particle, which results in simpler version than the corresponding Eulerian one, because the advection term disappears here. Our numerical examples successfully reproduce the classical results.展开更多
In practical fluid dynamic simulations, the bou n dary condition should be treated carefully because it always has crucial influen ce on the numerical accuracy, stability and efficiency. Two types of boundary tr eatme...In practical fluid dynamic simulations, the bou n dary condition should be treated carefully because it always has crucial influen ce on the numerical accuracy, stability and efficiency. Two types of boundary tr eatment methods for lattice Boltzmann method (LBM) are proposed. One is for the treatment of boundaries situated at lattice nodes, and the other is for the appr oximation of boundaries that are not located at the regular lattice nodes. The f irst type of boundary treatment method can deal with various dynamic boundaries on complex geometries by using a general set of formulas, which can maintain sec ond\|order accuracy. Based on the fact that the fluid flows simulated by LBM are not far from equilibrium, the unknown distributions at a boundary node are expr essed as the analogous forms of their corresponding equilibrium distributions. T herefore, the number of unknowns can be reduced and an always\|closed set of equ ations can be obtained for the solutions to pressure, velocity and special bound ary conditions on various geometries. The second type of boundary treatment is a complete interpolation scheme to treat curved boundaries. It comes from careful analysis of the relations between distribution functions at boundary nodes and their neighboring lattice nodes. It is stable for all situations and of second\| order accuracy. Basic ideas, implementation procedures and verifications with ty pical examples for the both treatments are presented. Numerical simulations and analyses show that they are accurate, stable, general and efficient for practica l simulations.展开更多
This paper presents a self-contained description on the configuration of propagator method(PM)to calculate the electron velocity distribution function(EVDF) of electron swarms in gases under DC electric and magnetic f...This paper presents a self-contained description on the configuration of propagator method(PM)to calculate the electron velocity distribution function(EVDF) of electron swarms in gases under DC electric and magnetic fields crossed at a right angle. Velocity space is divided into cells with respect to three polar coordinates v,θ and f. The number of electrons in each cell is stored in three-dimensional arrays. The changes of electron velocity due to acceleration by the electric and magnetic fields and scattering by gas molecules are treated as intercellular electron transfers on the basis of the Boltzmann equation and are represented using operators called the propagators or Green’s functions. The collision propagator, assuming isotropic scattering, is basically unchanged from conventional PMs performed under electric fields without magnetic fields. On the other hand, the acceleration propagator is customized for rotational acceleration under the action of the Lorentz force. The acceleration propagator specific to the present cell configuration is analytically derived. The mean electron energy and average electron velocity vector in a model gas and SF6 were derived from the EVDF as a demonstration of the PM under the Hall deflection and they were in a fine agreement with those obtained by Monte Carlo simulations. A strategy for fast relaxation is discussed, and extension of the PM for the EVDF under AC electric and DC/AC magnetic fields is outlined as well.展开更多
基金partially supported by the NSFC(11421061,12271507)the Natural Science Foundation of Shanghai(15ZR1403900)。
文摘This paper studies the strong convergence of the quantum lattice Boltzmann(QLB)scheme for the nonlinear Dirac equations for Gross-Neveu model in 1+1 dimensions.The initial data for the scheme are assumed to be convergent in L^(2).Then for any T≥0 the corresponding solutions for the quantum lattice Boltzmann scheme are shown to be convergent in C([0,T];L^(2)(R^(1)))to the strong solution to the nonlinear Dirac equations as the mesh sizes converge to zero.In the proof,at first a Glimm type functional is introduced to establish the stability estimates for the difference between two solutions for the corresponding quantum lattice Boltzmann scheme,which leads to the compactness of the set of the solutions for the quantum lattice Boltzmann scheme.Finally the limit of any convergent subsequence of the solutions for the quantum lattice Boltzmann scheme is shown to coincide with the strong solution to a Cauchy problem for the nonlinear Dirac equations.
基金Supported by the National Natural Science Foundation of China(No.20473034) the Taihu Scholar Foundation of SouthernYangtze University(2003).
文摘With the help of the method of separation of variables and the Debye-Hüchel approximation, the Poisson-Boltzmann equation that describes the distribution of the potential in the electrical double layer of a cylindrical particle with a limited length has been firstly solved under a very low potential condition. Then with the help of the functional analysis theory this equation has been further analytically solved under general potential conditions and consequently, the corresponding surface charge densities have been obtained. Both the potential and the surface charge densities cointide with those results obtained from the Debye-Hüchel approximation when the very low potential of zeψ〈〈kT is introduced.
文摘A Lagrangian lattice Boltzmann method for solving Euler equations is proposed. The key step in formulating this method is the introduction of the displacement distribution function. The equilibrium distribution function consists of macroscopic Lagrangian variables at time steps n and n + 1. It is different from the standard lattice Boltzmann method. In this method the element, instead of each particle, is required to satisfy the basic law. The element is considered as one large particle, which results in simpler version than the corresponding Eulerian one, because the advection term disappears here. Our numerical examples successfully reproduce the classical results.
文摘In practical fluid dynamic simulations, the bou n dary condition should be treated carefully because it always has crucial influen ce on the numerical accuracy, stability and efficiency. Two types of boundary tr eatment methods for lattice Boltzmann method (LBM) are proposed. One is for the treatment of boundaries situated at lattice nodes, and the other is for the appr oximation of boundaries that are not located at the regular lattice nodes. The f irst type of boundary treatment method can deal with various dynamic boundaries on complex geometries by using a general set of formulas, which can maintain sec ond\|order accuracy. Based on the fact that the fluid flows simulated by LBM are not far from equilibrium, the unknown distributions at a boundary node are expr essed as the analogous forms of their corresponding equilibrium distributions. T herefore, the number of unknowns can be reduced and an always\|closed set of equ ations can be obtained for the solutions to pressure, velocity and special bound ary conditions on various geometries. The second type of boundary treatment is a complete interpolation scheme to treat curved boundaries. It comes from careful analysis of the relations between distribution functions at boundary nodes and their neighboring lattice nodes. It is stable for all situations and of second\| order accuracy. Basic ideas, implementation procedures and verifications with ty pical examples for the both treatments are presented. Numerical simulations and analyses show that they are accurate, stable, general and efficient for practica l simulations.
文摘This paper presents a self-contained description on the configuration of propagator method(PM)to calculate the electron velocity distribution function(EVDF) of electron swarms in gases under DC electric and magnetic fields crossed at a right angle. Velocity space is divided into cells with respect to three polar coordinates v,θ and f. The number of electrons in each cell is stored in three-dimensional arrays. The changes of electron velocity due to acceleration by the electric and magnetic fields and scattering by gas molecules are treated as intercellular electron transfers on the basis of the Boltzmann equation and are represented using operators called the propagators or Green’s functions. The collision propagator, assuming isotropic scattering, is basically unchanged from conventional PMs performed under electric fields without magnetic fields. On the other hand, the acceleration propagator is customized for rotational acceleration under the action of the Lorentz force. The acceleration propagator specific to the present cell configuration is analytically derived. The mean electron energy and average electron velocity vector in a model gas and SF6 were derived from the EVDF as a demonstration of the PM under the Hall deflection and they were in a fine agreement with those obtained by Monte Carlo simulations. A strategy for fast relaxation is discussed, and extension of the PM for the EVDF under AC electric and DC/AC magnetic fields is outlined as well.
