Particular solutions of the linearized Boltzmann equation for a binary mixture of rigid spheres

Particular solutions that correspond to inhomogeneous driving terms in the linearized Boltzmann equation for the case of a binary mixture of rigid spheres are reported. For flow problems (in a plane channel) driven by pressure, temperature, and density gradients, inhomogeneous terms appear in the Boltzmann equation, and it is for these inhomogeneous terms that the particular solutions are developed. The required solutions for temperature and density driven problems are expressed in terms of previously reported generalized (vector-valued) Chapman–Enskog functions. However, for the pressure-driven problem (Poiseuille flow) the required particular solution is expressed in terms of two generalized Burnett functions defined by linear integral equations in which the driving terms are given in terms of the Chapman–Enskog functions. To complete this work, expansions in terms of Hermite cubic splines and a collocation scheme are used to establish numerical solutions for the generalized (vector-valued) Burnett functions.     

Particular solutions of the linearized Boltzmann equation for a binary mixture of rigid spheres

Particular solutions that correspond to inhomogeneous driving terms in the linearized Boltzmann equation for the case of a binary mixture of rigid spheres are reported. For flow problems (in a plane channel) driven by pressure, temperature, and density gradients, inhomogeneous terms appear in the Boltzmann equation, and it is for these inhomogeneous terms that the particular solutions are developed. The required solutions for temperature and density driven problems are expressed in terms of previously reported generalized (vector-valued) Chapman–Enskog functions. However, for the pressure-driven problem (Poiseuille flow) the required particular solution is expressed in terms of two generalized Burnett functions defined by linear integral equations in which the driving terms are given in terms of the Chapman–Enskog functions. To complete this work, expansions in terms of Hermite cubic splines and a collocation scheme are used to establish numerical solutions for the generalized (vector-valued) Burnett functions.     

Numerical resolution of a mono-disperse model of bubble growth in magmas

Growth of gas bubbles in magmas may be modeled by a system of differential equations that account for the evolution of bubble radius and internal pressure and that are coupled with an advection–diffusion equation defining the gas flux going from magma to bubble. This system of equations is characterized by two relaxation parameters linked to the viscosity of the magma and to the diffusivity of the dissolved gas, respectively. Here, we propose a numerical scheme preserving, by construction, the total mass of water of the system. We also study the asymptotic behavior of the system of equations by letting the relaxation parameters vary from 0 to ∞, and show the numerical convergence of the solutions obtained by means of the general numerical scheme to the simplified asymptotic limits. Finally, we validate and compare our numerical results with those obtained in experiments.     

A higher-order moment method of the lattice Boltzmann model for the conservation law equation

In this paper, we proposed a higher-order moment method in the lattice Boltzmann model for the conservation law equation. In contrast to the lattice Bhatnagar–Gross–Krook (BGK) model, the higher-order moment method has a wide flexibility to select equilibrium distribution function. This method is based on so-called a series of partial differential equations obtained by using multi-scale technique and Chapman–Enskog expansion. According to Hirt’s heuristic stability theory, the stability of the scheme can be controlled by modulating some special moments to design the third-order dispersion term and the fourth-order dissipation term. As results, the conservation law equation is recovered with higher-order truncation error. The numerical examples show the higher-order moment method can be used to raise the accuracy of the truncation error of the lattice Boltzmann scheme for the conservation law equation.     

Solution to the Boltzmann kinetic equation for high-speed flows

The Boltzmann kinetic equation is solved by a finite-difference method on a fixed coordinate-velocity grid. The projection method is applied that was developed previously by the author for evaluating the Boltzmann collision integral. The method ensures that the mass, momentum, and energy conservation laws are strictly satisfied and that the collision integral vanishes in thermodynamic equilibrium. The last property prevents the emergence of the numerical error when the collision integral of the principal part of the solution is evaluated outside Knudsen layers or shock waves, which considerably improves the accuracy and efficiency of the method. The differential part is approximated by a second-order accurate explicit conservative scheme. The resulting system of difference equations is solved by applying symmetric splitting into collision relaxation and free molecular flow. The steady-state solution is found by the relaxation method.     

Lattice Boltzmann model for two-dimensional generalized sine-Gordon equation

The nonlinear sine-Gordon equation arises in various problems in science and engineering. In this paper, we propose a numerical model based on lattice Boltmann method to obtain the numerical solutions of two-dimensional generalized sine-Gordon equation, including damped and undamped sine-Gordon equation. By choosing properly the conservation condition between the macroscopic quantity $u_t$ and the distribution functions and applying the Chapman-Enskog expansion, the governing equation is recovered correctly from the lattice Boltzmann equation. Moreover, the local equilibrium distribution function is obtained. The numerical results of the first three examples agree well with the analytic solutions, which indicates the lattice Boltzmann model is satisfactory and efficient. Numerical solutions for cases involving the most known from the bibliography line and ring solitons are given. Numerical experiments also show that the present scheme has a good long-time numerical behavior for the generalized sine-Gordon equation. Moreover, the model can also be applied to other two-dimensional nonlinear wave equations, such as nonlinear hyperbolic telegraph equation and Klein-Gordon equation.     

A lattice Boltzmann model for the eddy–stream equations in two-dimensional incompressible flows

A lattice Boltzmann model for two-dimensional incompressible flows with eddy–stream equations is proposed. By using two kinds of distribution functions and employing several higher-order moments of equilibrium distribution functions, the eddy equation and stream function equation with the second-order truncation error are obtained. In the numerical examples, we compared the numerical results of this scheme with those obtained by other classical method. The numerical results agree well with the classical ones.     

Relaxation approximation to bed-load sediment transport

In this work we propose and apply a numerical method based on finite volume relaxation approximation for computing the bed-load sediment transport in shallow water flows, in one and two space dimensions. The water flow is modeled by the well-known nonlinear shallow water equations which are coupled with a bed updating equation. Using a relaxation approximation, the nonlinear set of equations (and for two different formulations) is transformed to a semilinear diagonalizable problem with linear characteristic variables. A second order MUSCL-TVD method is used for the advection stage while an implicit–explicit Runge–Kutta scheme solves the relaxation stage. The main advantages of this approach are that neither Riemann problem solvers nor nonlinear iterations are required during the solution process. For the two different formulations, the applicability and effectiveness of the presented scheme is verified by comparing numerical results obtained for several benchmark test problems.     

Conservative numerical method for solving the averaged Boltzmann equation

A method is proposed for averaging the Boltzmann kinetic equation with respect to transverse velocities. A system of two integro-differential equations for two desired functions depending only on the longitudinal velocity is derived. The gas particles are assumed to interact as absolutely hard spheres. The integrals in the equations are double. The reduction in the number of variables in the desired functions and the low multiplicity of the integrals ensure the computational efficiency of the averaged equations. A numerical method of discrete ordinates is developed that effectively solves the gas relaxation problem based on the averaged equations. The method is conservative, and the number of particles, momentum, and energy are conserved automatically at every time step.     

Numerical simulations for the contraction flow using grid generation

We study the incomprssible Navier Stokes equations for the flow inside contraction geometry. The governing equations are expressed in the vorticity-stream function formulations. A rectangular computational domain is arised by elliptic grid generation technique. The numerical solution is based on a technique of automatic numerical generation of acurvilinear coordinate system by transforming the governing equation into computational plane. The transformed equations are approximated using central differences and solved simultaneously by successive over relaxation iteration. The time dependent of the vorticity equation solved by using explicit marching procedure. We will apply the technique on several irregularshapes.     

Application of a direct method for solving the Boltzmann equation in supersonic flow computation

The Boltzmann kinetic equation is used to numerically study the evolution of separated flows over a backward-facing step at low Knudsen numbers. The Boltzmann equation is solved by applying an explicit–implicit scheme. To improve the efficiency of the solution algorithm, it is parallelized with the help of MPI. The solution obtained with the kinetic equation is compared with those based on continuous medium equations. It is shown that the kinetic approach makes it possible to reproduce the distributions of surface pressure, friction coefficient, and heat transfer, as well as to obtain a flow structure close to experimental data.     

Study of the contraction flow using grid generation technique

We study a model procedure to solve the incompressible Navier-Stokes equations on the flow inside contraction geometry. The governing equations are expressed in the primitive variable formulation. A rectangular computational plane is arises by elliptic grid generation technique. The numerical solution is based on a technique of automatic numerical generation of a curvilinear coordinate system. By transformed the governing equation into computational plane. The time dependent momentum equations are solved explicitly for the velocity field using the explicit marching procedure, the continuity equation is applied at each grid point in the solution of pressure equation, while the successive over relaxation (SOR) method is used for the Neumann problem for pressure. We will apply the technique on several irregular-shape.     

Relaxation time simulation method with internal energy exchange for perfect gas flow at near-continuum conditions

This paper presents an internal energy exchange scheme for the relaxation time simulation method (RTSM) which solves the BGK equation for the perfect gas flow at near-continuum region discrete rotational energies are introduced to model the relaxation of internal energy modes. This development improved the agreements between RTSM and DSMC with little additional computational cost. The result shows a possibility of an improved hybrid RTSM/DSMC code for the continuum/rarefied gas flow.     

结合气体分子动力学方法的RKDG有限元格式求解一维Euler可压方程组

In this paper, two numerical methods are developed for solving one-dimensionl compressible ELder equations by the RKDG finite element method.The schemesare obtained based on an important relation between the Boltzmann equation andthe ELder equations.The schemes have the TVD-like property under the uniform meshes.Several numerical results also present the performance of the schemes.     

On a hyperbolic conservation law of electron transport in solid materials for electron probe microanalysis

The prediction of X-ray intensities based on the distribution of electrons throughout solid materials is essential to solve the inverse problem of quantifying the composition of materials in electron probe microanalysis (EPMA) [3]. We present a hyperbolic conservation law for electron transport in solid materials and investigate its validity under conditions typical for EPMA experiments. The conservation law is based on the time-stationary Boltzmann equation for binary electron-atom scattering. We model the energy loss of the electrons with a continuous slowing-down approximation. A first order moment approximation with respect to the angular variable is discussed. We propose to use a minimum entropy closure to derive a system of hyperbolic conservation laws, known as the M1 model [11]. A finite volume scheme for the numerical solution of the resulting equations is presented. Important numerical aspects of the scheme are discussed, such as bounds for the finite propagation speeds, as well as difficulties arising fromspatial discontinuities in thematerial coefficients and the scaling of the characteristic velocities with the stopping power of the electrons.We compare the accuracy and performance of the numerical solution of the hyperbolic conservation law to Monte Carlo simulations. The results indicate a reasonable accuracy of the proposed method and showthat compared to the MonteCarlo simulation the finite volume scheme is computationally less expensive.     

An iterative splitting approach for linear integro-differential equations

The motivation for this paper is to solve a model based on the dynamics of electrons in a plasma using a simplified Boltzmann equation. Such problems have arisen in active plasma resonance spectroscopy, which is used for plasma diagnostic techniques; see Braithwaite and Franklin (2009) [1]. We propose a modified iterative splitting approach to solve the Boltzmann equations as a system of integro-differential equations. To enable solution by fast and iterative computations, we first transform the integro-differential equations into second order differential equations. Second, we split each second order differential equations into two first order differential equations via a splitting approach. We carry out an error analysis of the higher order iterative approach. Numerical experiments with a simplified Boltzmann equation will be discussed, along with the benefits of computing with this splitting approach.     

Relaxation methods for mixed-integer optimal control of partial differential equations

We consider integer-restricted optimal control of systems governed by abstract semilinear evolution equations. This includes the problem of optimal control design for certain distributed parameter systems endowed with multiple actuators, where the task is to minimize costs associated with the dynamics of the system by choosing, for each instant in time, one of the actuators together with ordinary controls. We consider relaxation techniques that are already used successfully for mixed-integer optimal control of ordinary differential equations. Our analysis yields sufficient conditions such that the optimal value and the optimal state of the relaxed problem can be approximated with arbitrary precision by a control satisfying the integer restrictions. The results are obtained by semigroup theory methods. The approach is constructive and gives rise to a numerical method. We supplement the analysis with numerical experiments.     

Strong convergence of a stochastic Rosenbrock-type scheme for the finite element discretization of semilinear SPDEs driven by multiplicative and additive noise

This paper aims to investigate the numerical approximation of a general second order parabolic stochastic partial differential equation(SPDE) driven by multiplicative and additive noise. Our main interest is on such SPDEs where the nonlinear part is stronger than the linear part, usually called stochastic dominated transport equations. Most standard numerical schemes lose their good stability properties on such equations, including the current linear implicit Euler method. We discretize the SPDE in space by the finite element method and propose a novel scheme called stochastic Rosenbrock-type scheme for temporal discretization. Our scheme is based on the local linearization of the semi-discrete problem obtained after space discretization and is more appropriate for such equations. We provide a strong convergence of the new fully discrete scheme toward the exact solution for multiplicative and additive noise and obtain optimal rates of convergence. Numerical experiments to sustain our theoretical results are provided.     

二相流计算的一种差分算法

考虑两相流的力学行为,忽略相间的耗散作用,建立了Euler型的基本控制方程．状态方程采用刚性状态方程．基于Abgrall提出的准则,在流动区域内,对可压两相流提出了一个精度较高的Euler型数值方法,数值格式是Godunov型格式,对守恒型和非守恒型方程采用HLLC型和Lax-Friedrichs型近似Riemann解算器,引入了速度驰豫和压强驰豫过程来代替两相间的相互作用．在一维情形下给出数值算例,并且和Saurel的算例进行了比较,结果表明该算法不但精确而且稳定,且在间断处没有数值振荡．