In this article, we study the smoothing effect of the Cauchy problem for the spatially homogeneous non-cutoff Boltzmann equation for hard potentials. It has long been suspected that the non-cutoff Boltzmann equation e...In this article, we study the smoothing effect of the Cauchy problem for the spatially homogeneous non-cutoff Boltzmann equation for hard potentials. It has long been suspected that the non-cutoff Boltzmann equation enjoys similar regularity properties as to whose of the fractional heat equation. We prove that any solution with mild regularity will become smooth in Gevrey class at positive time, with a sharp Gevrey index, depending on the angular singularity. Our proof relies on the elementary L^(2) weighted estimates.展开更多
The author generalized the propagator function theory introduced first by Sigmund, and gave a explicitly proof of a equivalence between forward and backward Boltzmann equations in a multi component medium by using the...The author generalized the propagator function theory introduced first by Sigmund, and gave a explicitly proof of a equivalence between forward and backward Boltzmann equations in a multi component medium by using the generalized propagator function theory.展开更多
We study the acoustomagnetoelectric (AME) effect in two-dimensional graphene with an energy bandgap using the semiclassical Boltzmann transport equation within the hypersound regime, (where represents the acoustic wav...We study the acoustomagnetoelectric (AME) effect in two-dimensional graphene with an energy bandgap using the semiclassical Boltzmann transport equation within the hypersound regime, (where represents the acoustic wavenumber and is the mean free path of the electron). The Boltzmann transport equation and other relevant equations were solved analytically to obtain an expression for the AME current density, consisting of longitudinal and Hall components. Our numerical results indicate that both components of the AME current densities display oscillatory behaviour. Furthermore, geometric resonances and Weiss oscillations were each defined using the relationship between the current density and Surface Acoustic Wave (SAW) frequency and the inverse of the applied magnetic field, respectively. Our results show that the AME current density of bandgap graphene, which can be controlled to suit a particular electronic device application, is smaller than that of (gapless) graphene and is therefore, more suited for nanophotonic device applications.展开更多
In the present work,both computational and experimentalmethods are employed to study the two-phase flow occurring in a model pump sump.The twofluid model of the two-phase flow has been applied to the simulation of the...In the present work,both computational and experimentalmethods are employed to study the two-phase flow occurring in a model pump sump.The twofluid model of the two-phase flow has been applied to the simulation of the threedimensional cavitating flow.The governing equations of the two-phase cavitating flow are derived from the kinetic theory based on the Boltzmann equation.The isotropic RNG k−ε−kca turbulence model of two-phase flows in the form of cavity number instead of the formof cavity phase volume fraction is developed.The RNG k−ε−kca turbulence model,that is the RNG k−εturbulence model for the liquid phase combined with the kca model for the cavity phase,is employed to close the governing turbulent equations of the two-phase flow.The computation of the cavitating flow through a model pump sump has been carried out with this model in three-dimensional spaces.The calculated results have been compared with the data of the PIV experiment.Good qualitative agreement has been achievedwhich exhibits the reliability of the numerical simulation model.展开更多
The existence and uniqueness of the solutions for the Boltzmann equations with measures as initial value are still an open problem which is posed by P. L. Lions (2000). The aim of this paper is to discuss the Cauchy...The existence and uniqueness of the solutions for the Boltzmann equations with measures as initial value are still an open problem which is posed by P. L. Lions (2000). The aim of this paper is to discuss the Cauchy problem of the system of discrete Boltzmann equations of the formδtfi+(fi^mi)=Qi(f1,f2...,fn),(mi〉1,i-1,...n)with non-negative finite Radon measures as initial conditions. In particular, the existence and uniqueness of BV solutions for the above problem are obtained.展开更多
In this article we study the Gevrey regularization effect for the spatially inhomogeneous Boltzmann equation without angular cut-off.This equation is partially elliptic in the velocity direction and degenerates in the...In this article we study the Gevrey regularization effect for the spatially inhomogeneous Boltzmann equation without angular cut-off.This equation is partially elliptic in the velocity direction and degenerates in the spatial variable.We consider the nonlinear Cauchy problem for the fluctuation around the Maxwellian distribution and prove that any solution with mild regularity will become smooth in the Gevrey class at positive time with the Gevrey index depending on the angular singularity.Our proof relies on the symbolic calculus for the collision operator and the global subelliptic estimate for the Cauchy problem of the linearized Boltzmann operator.展开更多
A previous study is continued by investigating the Boltzmann equation for particles with quantum effects (BQE). First, the corresponding entropy identity is proved, then if the initial data f(x,v,0) satisfies 0...A previous study is continued by investigating the Boltzmann equation for particles with quantum effects (BQE). First, the corresponding entropy identity is proved, then if the initial data f(x,v,0) satisfies 0≤f(x,v,0)≤CΦ(x,v,0) for a constant 0<C<∞ and function Φ(x,v,t), we prove the existence and uniqueness of spatial decay solutions of the BQE within a given function space B(Φ) using fixed point theory. Moreover, if there is a continuous function F(x,v) which belongs to a function set, then there exists a mild solution f(x,v,t) of the BQE such that f ∞(x,v)= limt→∞f(x+vt,v,t)=F(x,v).展开更多
The Boltzmann equations for Fermi-Dirac particles and Bose-Einstein particles, both in the absence of external force fields, are combined into a more general form called the Boltzmann equation with quantum effects (BQ...The Boltzmann equations for Fermi-Dirac particles and Bose-Einstein particles, both in the absence of external force fields, are combined into a more general form called the Boltzmann equation with quantum effects (BQE). It is assumed that the initial data f(x,v,0) satisfies 0≤f(x,v,0)≤cΦ(x,v,0) for a positive constant c and certain types of control functions Φ(x,v,t). Then within a given function space B(Φ), we prove that f(x+tv,v,t) uniformly converges to f ∞(x,v) in a certain norm where f ∞(x,v)= limt→∞f(x+tv,v,t) and different initial data determines different long time limits.展开更多
A class of high resolution positivity preserving Boltzmann type difference schemes for one and two dimensional Euler equations is studied. First, the relation between Boltzmann and Euler equations is analyzed. By usi...A class of high resolution positivity preserving Boltzmann type difference schemes for one and two dimensional Euler equations is studied. First, the relation between Boltzmann and Euler equations is analyzed. By using a kind of special interpolation, the high resolution Boltzmann type difference scheme is constructed. Finally, numerical tests show that the schemes are effective and useful.展开更多
A lattice Boltzmann flux solver(LBFS)is presented for simulation of fluid flows.Like the conventional computational fluid dynamics(CFD)solvers,the new solver also applies the finite volume method to discretize the gov...A lattice Boltzmann flux solver(LBFS)is presented for simulation of fluid flows.Like the conventional computational fluid dynamics(CFD)solvers,the new solver also applies the finite volume method to discretize the governing differential equations,but the numerical flux at the cell interface is not evaluated by the smooth function approximation or Riemann solvers.Instead,it is evaluated from local solution of lattice Boltzmann equation(LBE)at cell interface.Two versions of LBFS are presented in this paper.One is to locally apply one-dimensional compressible lattice Boltzmann(LB)model along the normal direction to the cell interface for simulation of compressible inviscid flows with shock waves.The other is to locally apply multi-dimensional LB model at cell interface for simulation of incompressible viscous and inviscid flows.The present solver removes the drawbacks of conventional lattice Boltzmann method(LBM)such as limitation to uniform mesh,tie-up of mesh spacing and time interval,limitation to viscous flows.Numerical examples show that the present solver can be well applied to simulate fluid flows with non-uniform mesh and curved boundary.展开更多
In this paper, we propose a novel incompressible finite-difference lattice Boltzmann Equation (FDLBE). Because source terms that reflect the interaction between phases can be accurately described, the new model is s...In this paper, we propose a novel incompressible finite-difference lattice Boltzmann Equation (FDLBE). Because source terms that reflect the interaction between phases can be accurately described, the new model is suitable for simulating two-way coupling incompressible multiphase flow The 2-D particle-laden flow over a backward-facing step is chosen as a test case to validate the present method. Favorable results are obtained and the present scheme is shown to have good prospects in practical applications.展开更多
In this paper, the dielectric properties of CO2, CO2/air, CO2/O2, CO2/N2, CO2/CF4, CO2/CH4, CO2/He, C02/H2, CO2/NH3 and CO2/CO were investigated based on the Boltzmann equation analysis, in which the reduced critical ...In this paper, the dielectric properties of CO2, CO2/air, CO2/O2, CO2/N2, CO2/CF4, CO2/CH4, CO2/He, C02/H2, CO2/NH3 and CO2/CO were investigated based on the Boltzmann equation analysis, in which the reduced critical electric field strength (E/N)cr of the gases was derived from the calculated electron energy distribution function (EEDF) by solv- ing the Boltzmann transport equation. In this work, it should be noted that the fundamental data were carefully selected by the published experimental results and calculations to ensure the validity of the calculation. The results indicate that if He, H2, N2 and CH4, in which there axe high ionization coefficients or a lack of attachment reactions, are added into CO2, the dielectric properties will decrease. On the other hand, air, O2, NH3 and CFa (ranked in terms of (E/N)cr value in increasing order) have the potential to improve the dielectric property of CO2 at room temperature.展开更多
In this article, we study the Cauchy problem for the linearized spatially homogeneous non-cutoff Boltzamnn equation with Maxwellian molecules. By using the spectral decomposition, we solve the Cauchy problem with init...In this article, we study the Cauchy problem for the linearized spatially homogeneous non-cutoff Boltzamnn equation with Maxwellian molecules. By using the spectral decomposition, we solve the Cauchy problem with initial datum in the sense of distribution, which contains the dual space of a Gelfand-Shilov class. We also prove that this solution belongs to the Gelfand-Shilov space for any positive time.展开更多
A method combining the immersed boundary technique and a multi- relaxation-time (MRT) lattice Boltzmann flux solver (LBFS) is presented for numerical simulation of incompressible flows over circular and elliptic c...A method combining the immersed boundary technique and a multi- relaxation-time (MRT) lattice Boltzmann flux solver (LBFS) is presented for numerical simulation of incompressible flows over circular and elliptic cylinders and NACA 0012 Airfoil. The method uses a simple Cartesian mesh to simulate flows past immersed complicated bodies. With the Chapman-Enskog expansion analysis, a transform is performed between the Navier-Stokes and lattice Boltzmann equations (LBEs). The LBFS is used to discretize the macroscopic differential equations with a finite volume method and evaluate the interface fluxes through local reconstruction of the lattice Boltzmann solution. The immersed boundary technique is used to correct the intermediate velocity around the solid boundary to satisfy the no-slip boundary condition. Agreement of simulation results with the data found in the literature shows reliability of the proposed method in simulating laminar flows on a Cartesian mesh.展开更多
The regularity of solutions to the Boltzmann equation is a fundamental problem in the kinetic theory. In this paper, the case with angular cut-off is investigated. It is shown that the macroscopic parts of solutions t...The regularity of solutions to the Boltzmann equation is a fundamental problem in the kinetic theory. In this paper, the case with angular cut-off is investigated. It is shown that the macroscopic parts of solutions to the Boltzmann equation, i.e., the density, momentum and total energy are continuous functions of (x, t) in the region R3 × (0, +∞). More precisely, these macroscopic quantities immediately become continuous in any positive time even though they are initially discontinuous and the discontinuities of solutions propagate only in the microscopic level. It should be noted that such kind of phenomenon can not hap- pen for the compressible Navier-Stokes equations in which the initial discontinuities of the density never vanish in any finite time, see [22]. This hints that the Boltzmann equation has better regularity effect in the macroscopic level than compressible Navier-Stokes equations.展开更多
The Boltzmann-Bhatnagar-Gross-Krook(BGK)model is investigated for its validity regarding the collision term approximation through relaxation evaluation. The evaluation is based on theoretical analysis and numerical ...The Boltzmann-Bhatnagar-Gross-Krook(BGK)model is investigated for its validity regarding the collision term approximation through relaxation evaluation. The evaluation is based on theoretical analysis and numerical comparison between the BGK and direct simulation Monte Carlo(DSMC) results for three specifically designed relaxation problems. In these problems, one or half component of the velocity distribution is characterized by another Maxwellian distribution with a different temperature. It is analyzed that the relaxation time in the BGK model is unequal to the molecular mean collision time. Relaxation of component distribution fails to involve enough contribution from other component distributions, which makes the BGK model unable to capture details of velocity distribution, especially when discontinuity exists in distribution. The BGK model,however, predicts satisfactory results including fluxes during relaxation when the temperature difference is small. Particularly, the model-induced error in the BGK model increases with the temperature difference, thus the model is more reliable for low-speed rarefied flows than for hypersonic flows.展开更多
Interaction between turbulence and particles is investigated in a channel flow. The fluid motion is calculated using direct numerical simulation(DNS) with a lattice Boltzmann(LB) method, and particles are tracked in a...Interaction between turbulence and particles is investigated in a channel flow. The fluid motion is calculated using direct numerical simulation(DNS) with a lattice Boltzmann(LB) method, and particles are tracked in a Lagrangian framework through the action of force imposed by the fluid. The particle diameter is smaller than the Kolmogorov length scale, and the point force is used to represent the feedback force of particles on the turbulence. The effects of particles on the turbulence and skin friction coefficient are examined with different particle inertias and mass loadings. Inertial particles suppress intensities of the spanwise and wall-normal components of velocity, and the Reynolds shear stress. It is also found that, relative to the reference particle-free flow,the overall mean skin-friction coefficient is reduced by particles. Changes of near wall turbulent structures such as longer and more regular streamwise low-speed streaks and less ejections and sweeps are the manifestation of drag reduction.展开更多
In either a periodic box Td or Rd (1 ≤ d ≤ 3), we obtain global bounded solution of the relativistic Boltzmann equation near global relativistic Maxwellian, in terms of natural mass, energy conservation and the en...In either a periodic box Td or Rd (1 ≤ d ≤ 3), we obtain global bounded solution of the relativistic Boltzmann equation near global relativistic Maxwellian, in terms of natural mass, energy conservation and the entropy inequality.展开更多
Aim.The well known JST(Jameson-Schmidt-Turkel) scheme requires the use of a dissipation term.We propose using gas-kinetic BGK(Bhatnagar-Gross-Krook) method,which is based on the more fundamental Boltzmann equation,in ...Aim.The well known JST(Jameson-Schmidt-Turkel) scheme requires the use of a dissipation term.We propose using gas-kinetic BGK(Bhatnagar-Gross-Krook) method,which is based on the more fundamental Boltzmann equation,in order to obviate the use of dissipation term and obtain,we believe,an improved solution.Section 1 deals essentially with three things:(1) as analytical solution of molecular probability density function at the cell interface has been obtained by the Boltzmann equation with BGK model,we can compute the flux term by integrating the density function in the phase space;eqs.(8) and(11) require careful attention;(2) the integrations can be expressed as the moments of Maxwellian distribution with different limits according to the analytical solution;eqs.(9) and(10) require careful attention;(3) the discrete equation by finite volume method can be solved using the time marching method.Computations are performed by the BGK method for the Sod′s shock tube problem and a two-dimensional shock reflection problem.The results are compared with those of the conventional JST scheme in Figs.1 and 2.The BGK method provides better resolution of shock waves and other features of the flow fields.展开更多
It is commonly known that the hydrodynamic equations can be derived from the Boltzmann equation. In this paper, we derive similar spin-dependent balance equations based on the spinor Boltzmann equation. Besides the us...It is commonly known that the hydrodynamic equations can be derived from the Boltzmann equation. In this paper, we derive similar spin-dependent balance equations based on the spinor Boltzmann equation. Besides the usual charge current, heat current, and pressure tensor, we also explore the characteristic spin accumulation and spin current as well as the spin-dependent pressure tensor and heat current in spintronics. The numerical results of these physical quantities are demonstrated using an example of spin-polarized transport through a mesoscopic ferromagnet.展开更多
基金supported by the NSFC(12101012)the PhD Scientific Research Start-up Foundation of Anhui Normal University.Zeng’s research was supported by the NSFC(11961160716,11871054,12131017).
文摘In this article, we study the smoothing effect of the Cauchy problem for the spatially homogeneous non-cutoff Boltzmann equation for hard potentials. It has long been suspected that the non-cutoff Boltzmann equation enjoys similar regularity properties as to whose of the fractional heat equation. We prove that any solution with mild regularity will become smooth in Gevrey class at positive time, with a sharp Gevrey index, depending on the angular singularity. Our proof relies on the elementary L^(2) weighted estimates.
基金The scientific research fund from the Educational Bureau of Anhui ProvinceP.R.C.(2 0 0 1Kj2 2 6)
文摘The author generalized the propagator function theory introduced first by Sigmund, and gave a explicitly proof of a equivalence between forward and backward Boltzmann equations in a multi component medium by using the generalized propagator function theory.
文摘We study the acoustomagnetoelectric (AME) effect in two-dimensional graphene with an energy bandgap using the semiclassical Boltzmann transport equation within the hypersound regime, (where represents the acoustic wavenumber and is the mean free path of the electron). The Boltzmann transport equation and other relevant equations were solved analytically to obtain an expression for the AME current density, consisting of longitudinal and Hall components. Our numerical results indicate that both components of the AME current densities display oscillatory behaviour. Furthermore, geometric resonances and Weiss oscillations were each defined using the relationship between the current density and Surface Acoustic Wave (SAW) frequency and the inverse of the applied magnetic field, respectively. Our results show that the AME current density of bandgap graphene, which can be controlled to suit a particular electronic device application, is smaller than that of (gapless) graphene and is therefore, more suited for nanophotonic device applications.
基金This research work was funded by the Chinese National Foundation of Natural Science(Nos.51076077,51176168,51249003 and 51076144)This work is also supported by Science Foundation of Zhejiang Sci-Tech University(ZSTU) under Grant No.11130032241201.
文摘In the present work,both computational and experimentalmethods are employed to study the two-phase flow occurring in a model pump sump.The twofluid model of the two-phase flow has been applied to the simulation of the threedimensional cavitating flow.The governing equations of the two-phase cavitating flow are derived from the kinetic theory based on the Boltzmann equation.The isotropic RNG k−ε−kca turbulence model of two-phase flows in the form of cavity number instead of the formof cavity phase volume fraction is developed.The RNG k−ε−kca turbulence model,that is the RNG k−εturbulence model for the liquid phase combined with the kca model for the cavity phase,is employed to close the governing turbulent equations of the two-phase flow.The computation of the cavitating flow through a model pump sump has been carried out with this model in three-dimensional spaces.The calculated results have been compared with the data of the PIV experiment.Good qualitative agreement has been achievedwhich exhibits the reliability of the numerical simulation model.
文摘The existence and uniqueness of the solutions for the Boltzmann equations with measures as initial value are still an open problem which is posed by P. L. Lions (2000). The aim of this paper is to discuss the Cauchy problem of the system of discrete Boltzmann equations of the formδtfi+(fi^mi)=Qi(f1,f2...,fn),(mi〉1,i-1,...n)with non-negative finite Radon measures as initial conditions. In particular, the existence and uniqueness of BV solutions for the above problem are obtained.
基金supported by National Natural Science Foundation of China(Grant No.11631011)supported by National Natural Science Foundation of China(Grant Nos.11961160716,11871054 and 11771342)+1 种基金the Natural Science Foundation of Hubei Province(Grant No.2019CFA007)the Fundamental Research Funds for the Central Universities(Grant No.2042020kf0210)。
文摘In this article we study the Gevrey regularization effect for the spatially inhomogeneous Boltzmann equation without angular cut-off.This equation is partially elliptic in the velocity direction and degenerates in the spatial variable.We consider the nonlinear Cauchy problem for the fluctuation around the Maxwellian distribution and prove that any solution with mild regularity will become smooth in the Gevrey class at positive time with the Gevrey index depending on the angular singularity.Our proof relies on the symbolic calculus for the collision operator and the global subelliptic estimate for the Cauchy problem of the linearized Boltzmann operator.
基金Supported by the Tsinghua U niversity Science Fund
文摘A previous study is continued by investigating the Boltzmann equation for particles with quantum effects (BQE). First, the corresponding entropy identity is proved, then if the initial data f(x,v,0) satisfies 0≤f(x,v,0)≤CΦ(x,v,0) for a constant 0<C<∞ and function Φ(x,v,t), we prove the existence and uniqueness of spatial decay solutions of the BQE within a given function space B(Φ) using fixed point theory. Moreover, if there is a continuous function F(x,v) which belongs to a function set, then there exists a mild solution f(x,v,t) of the BQE such that f ∞(x,v)= limt→∞f(x+vt,v,t)=F(x,v).
基金Supported by the Tsinghua U niversity Science Fund
文摘The Boltzmann equations for Fermi-Dirac particles and Bose-Einstein particles, both in the absence of external force fields, are combined into a more general form called the Boltzmann equation with quantum effects (BQE). It is assumed that the initial data f(x,v,0) satisfies 0≤f(x,v,0)≤cΦ(x,v,0) for a positive constant c and certain types of control functions Φ(x,v,t). Then within a given function space B(Φ), we prove that f(x+tv,v,t) uniformly converges to f ∞(x,v) in a certain norm where f ∞(x,v)= limt→∞f(x+tv,v,t) and different initial data determines different long time limits.
文摘A class of high resolution positivity preserving Boltzmann type difference schemes for one and two dimensional Euler equations is studied. First, the relation between Boltzmann and Euler equations is analyzed. By using a kind of special interpolation, the high resolution Boltzmann type difference scheme is constructed. Finally, numerical tests show that the schemes are effective and useful.
基金Supported by the National Natural Science Foundation of China(11272153)
文摘A lattice Boltzmann flux solver(LBFS)is presented for simulation of fluid flows.Like the conventional computational fluid dynamics(CFD)solvers,the new solver also applies the finite volume method to discretize the governing differential equations,but the numerical flux at the cell interface is not evaluated by the smooth function approximation or Riemann solvers.Instead,it is evaluated from local solution of lattice Boltzmann equation(LBE)at cell interface.Two versions of LBFS are presented in this paper.One is to locally apply one-dimensional compressible lattice Boltzmann(LB)model along the normal direction to the cell interface for simulation of compressible inviscid flows with shock waves.The other is to locally apply multi-dimensional LB model at cell interface for simulation of incompressible viscous and inviscid flows.The present solver removes the drawbacks of conventional lattice Boltzmann method(LBM)such as limitation to uniform mesh,tie-up of mesh spacing and time interval,limitation to viscous flows.Numerical examples show that the present solver can be well applied to simulate fluid flows with non-uniform mesh and curved boundary.
基金The project supported by the National Natural Science Foundation of China(60073044)the State Key Development Programme for Basic Research of China(G1990022207).
文摘In this paper, we propose a novel incompressible finite-difference lattice Boltzmann Equation (FDLBE). Because source terms that reflect the interaction between phases can be accurately described, the new model is suitable for simulating two-way coupling incompressible multiphase flow The 2-D particle-laden flow over a backward-facing step is chosen as a test case to validate the present method. Favorable results are obtained and the present scheme is shown to have good prospects in practical applications.
基金supported in part by the National Key Basic Research Program of China(973 Program)(No.2015CB251002)the Science and Technology Project Funds of the Grid State Corporation of China(No.SGSNK00KJJS1501564)+2 种基金National Natural Science Foundation of China(Nos.51221005,51577145)the Fundamental Research Funds for the Central Universities of Chinathe Program for New Century Excellent Talents in University,China
文摘In this paper, the dielectric properties of CO2, CO2/air, CO2/O2, CO2/N2, CO2/CF4, CO2/CH4, CO2/He, C02/H2, CO2/NH3 and CO2/CO were investigated based on the Boltzmann equation analysis, in which the reduced critical electric field strength (E/N)cr of the gases was derived from the calculated electron energy distribution function (EEDF) by solv- ing the Boltzmann transport equation. In this work, it should be noted that the fundamental data were carefully selected by the published experimental results and calculations to ensure the validity of the calculation. The results indicate that if He, H2, N2 and CH4, in which there axe high ionization coefficients or a lack of attachment reactions, are added into CO2, the dielectric properties will decrease. On the other hand, air, O2, NH3 and CFa (ranked in terms of (E/N)cr value in increasing order) have the potential to improve the dielectric property of CO2 at room temperature.
基金supported by the Fundamental Research Funds for the Central Unversities and National Science Foundation of China(11171261and 11422106)
文摘In this article, we study the Cauchy problem for the linearized spatially homogeneous non-cutoff Boltzamnn equation with Maxwellian molecules. By using the spectral decomposition, we solve the Cauchy problem with initial datum in the sense of distribution, which contains the dual space of a Gelfand-Shilov class. We also prove that this solution belongs to the Gelfand-Shilov space for any positive time.
文摘A method combining the immersed boundary technique and a multi- relaxation-time (MRT) lattice Boltzmann flux solver (LBFS) is presented for numerical simulation of incompressible flows over circular and elliptic cylinders and NACA 0012 Airfoil. The method uses a simple Cartesian mesh to simulate flows past immersed complicated bodies. With the Chapman-Enskog expansion analysis, a transform is performed between the Navier-Stokes and lattice Boltzmann equations (LBEs). The LBFS is used to discretize the macroscopic differential equations with a finite volume method and evaluate the interface fluxes through local reconstruction of the lattice Boltzmann solution. The immersed boundary technique is used to correct the intermediate velocity around the solid boundary to satisfy the no-slip boundary condition. Agreement of simulation results with the data found in the literature shows reliability of the proposed method in simulating laminar flows on a Cartesian mesh.
基金partially supported by by National Center for Mathematics and Interdisciplinary Sciences,AMSS,CAS and NSFC(11371349 and 11688101)partially supported by NSFC(11688101,11771429)
文摘The regularity of solutions to the Boltzmann equation is a fundamental problem in the kinetic theory. In this paper, the case with angular cut-off is investigated. It is shown that the macroscopic parts of solutions to the Boltzmann equation, i.e., the density, momentum and total energy are continuous functions of (x, t) in the region R3 × (0, +∞). More precisely, these macroscopic quantities immediately become continuous in any positive time even though they are initially discontinuous and the discontinuities of solutions propagate only in the microscopic level. It should be noted that such kind of phenomenon can not hap- pen for the compressible Navier-Stokes equations in which the initial discontinuities of the density never vanish in any finite time, see [22]. This hints that the Boltzmann equation has better regularity effect in the macroscopic level than compressible Navier-Stokes equations.
基金supported by the National Natural Science Foundation of China(91116013,11372325,and 11111120080)
文摘The Boltzmann-Bhatnagar-Gross-Krook(BGK)model is investigated for its validity regarding the collision term approximation through relaxation evaluation. The evaluation is based on theoretical analysis and numerical comparison between the BGK and direct simulation Monte Carlo(DSMC) results for three specifically designed relaxation problems. In these problems, one or half component of the velocity distribution is characterized by another Maxwellian distribution with a different temperature. It is analyzed that the relaxation time in the BGK model is unequal to the molecular mean collision time. Relaxation of component distribution fails to involve enough contribution from other component distributions, which makes the BGK model unable to capture details of velocity distribution, especially when discontinuity exists in distribution. The BGK model,however, predicts satisfactory results including fluxes during relaxation when the temperature difference is small. Particularly, the model-induced error in the BGK model increases with the temperature difference, thus the model is more reliable for low-speed rarefied flows than for hypersonic flows.
基金Project supported by the National Natural Science Foundation of China(Nos.11572183 and 11272198)
文摘Interaction between turbulence and particles is investigated in a channel flow. The fluid motion is calculated using direct numerical simulation(DNS) with a lattice Boltzmann(LB) method, and particles are tracked in a Lagrangian framework through the action of force imposed by the fluid. The particle diameter is smaller than the Kolmogorov length scale, and the point force is used to represent the feedback force of particles on the turbulence. The effects of particles on the turbulence and skin friction coefficient are examined with different particle inertias and mass loadings. Inertial particles suppress intensities of the spanwise and wall-normal components of velocity, and the Reynolds shear stress. It is also found that, relative to the reference particle-free flow,the overall mean skin-friction coefficient is reduced by particles. Changes of near wall turbulent structures such as longer and more regular streamwise low-speed streaks and less ejections and sweeps are the manifestation of drag reduction.
文摘In either a periodic box Td or Rd (1 ≤ d ≤ 3), we obtain global bounded solution of the relativistic Boltzmann equation near global relativistic Maxwellian, in terms of natural mass, energy conservation and the entropy inequality.
文摘Aim.The well known JST(Jameson-Schmidt-Turkel) scheme requires the use of a dissipation term.We propose using gas-kinetic BGK(Bhatnagar-Gross-Krook) method,which is based on the more fundamental Boltzmann equation,in order to obviate the use of dissipation term and obtain,we believe,an improved solution.Section 1 deals essentially with three things:(1) as analytical solution of molecular probability density function at the cell interface has been obtained by the Boltzmann equation with BGK model,we can compute the flux term by integrating the density function in the phase space;eqs.(8) and(11) require careful attention;(2) the integrations can be expressed as the moments of Maxwellian distribution with different limits according to the analytical solution;eqs.(9) and(10) require careful attention;(3) the discrete equation by finite volume method can be solved using the time marching method.Computations are performed by the BGK method for the Sod′s shock tube problem and a two-dimensional shock reflection problem.The results are compared with those of the conventional JST scheme in Figs.1 and 2.The BGK method provides better resolution of shock waves and other features of the flow fields.
基金Project supported by the National Natural Science Foundation of China(Grant No.11274378)the Key Research Program of the Chinese Academy of Sciences(Grant No.XDPB08-3)the MOST of China(Grant No.2013CB933401)
文摘It is commonly known that the hydrodynamic equations can be derived from the Boltzmann equation. In this paper, we derive similar spin-dependent balance equations based on the spinor Boltzmann equation. Besides the usual charge current, heat current, and pressure tensor, we also explore the characteristic spin accumulation and spin current as well as the spin-dependent pressure tensor and heat current in spintronics. The numerical results of these physical quantities are demonstrated using an example of spin-polarized transport through a mesoscopic ferromagnet.