Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userSemi-Lagrangian implicit Bhatnagar-Gross-Krook collision model for the finite-volume discrete Boltzmann method
67
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
A new implicit BGK collision model using a semi-Lagrangian approach is proposed in this paper. Unlike existing models, in which the implicit BGK collision is resolved either by a temporal extrapolation or by a variable transformation, the new model removes the implicitness by tracing the particle distribution functions (PDFs) back in time along their characteristic paths during the collision process. An interpolation scheme is needed to evaluate the PDFs at the traced-back locations. By using the first-order interpolation, the resulting model allows for the straightforward replacement of ${f_{alpha}}^{eq,n+1}$ by ${f_{alpha}}^{eq,n}$ no matter where it appears. After comparing the new model with the existing models under different numerical conditions (e.g. different flux schemes and time marching schemes) and using the new model to successfully modify the variable transformation technique, three conclusions can be drawn. First, the new model can improve the accuracy by almost an order of magnitude. Second, it can slightly reduce the computational cost. Therefore, the new scheme improves accuracy without extra cost. Finally, the new model can significantly improve the ${Delta}t/{tau}$ limit compared to the temporal interpolation model while having the same ${Delta}t/{tau}$ limit as the variable transformation approach. The new scheme with a second-order interpolation is also developed and tested; however, that technique displays no advantage over the simple first-order interpolation approach. Both numerical and theoretical analyses are also provided to explain why the new implicit scheme with simple first-order interpolation can outperform the same scheme with second-order interpolation, as well as the existing temporal extrapolation and variable transformation schemes.
rate research
Read More
A two-fluid Discrete Boltzmann Model(DBM) for compressible flows based on Ellipsoidal Statistical Bhatnagar-Gross-Krook(ES-BGK) is presented. The model has flexible Prandtl number or specific heat ratio. Mathematically, the model is composed of two coupled Discrete Boltzmann Equations(DBE). Each DBE describes one component of the fluid. Physically, the model is equivalent to a macroscopic fluid model based on Navier-Stokes(NS) equations, and supplemented by a coarse-grained model for thermodynamic non-equilibrium behaviors. To obtain a flexible Prandtl number, a coefficient is introduced in the ellipsoidal statistical distribution function to control the viscosity. To obtain a flexible specific heat ratio, a parameter is introduced in the energy kinetic moments to control the extra degree of freedom. For binary mixture, the correspondence between the macroscopic fluid model and the DBM may be several-to-one. Five typical benchmark tests are used to verify and validate the model. Some interesting non-equilibrium results, which are not available in the NS model or the single-fluid DBM, are presented.
The discrete effect on the boundary condition has been a fundamental topic for the lattice Boltzmann method in simulating heat and mass transfer problems. In previous works based on the halfway anti-bounce-back (ABB) boundary condition for convection-diffusion equations (CDEs), it is reported that the discrete effect cannot be commonly removed in the Bhatnagar-Gross-Krook (BGK) model except for a special value of relaxation time. Targeting this point in the present paper, we still proceed within the framework of BGK model for two-dimensional CDEs, and analyze the discrete effect on a non-halfway ABB boundary condition which incorporates the effect of the distance ratio. By analyzing an unidirectional diffusion problem with a parabolic distribution, the theoretical derivations with three different discrete velocity models show that the numerical slip is a combined function of the relaxation time and the distance ratio. Different from previous works, we definitely find that the relaxation time can be freely adjusted by the distance ratio in a proper range to eliminate the numerical slip. Some numerical simulations are carried out to validate the theoretical derivations, and the numerical results for the cases of straight and curved boundaries confirm our theoretical analysis. Finally, it should be noted that the present analysis can be extended from the BGK model to other lattice Boltzmann (LB) collision models for CDEs, which can broaden the parameter range of the relaxation time to approach 0.5.
We present a multi-scale lattice Boltzmann scheme, which adaptively refines particles velocity space. Different velocity sets, i.e., higher- and lower-order lattices, are consistently and efficiently coupled, allowing us to use the higher-order lattice only when and where needed. This includes regions of either high Mach number or high Knudsen number. The coupling procedure of different lattices consists of either projection of the moments of the higher-order lattice onto the lower-order lattice or lifting of the lower-order lattice to the higher-order velocity space. Both lifting and projection are local operations, which enable a flexible adaptive velocity set. The proposed scheme can be formulated both in a static and an optimal, co-moving reference frame, in the spirit of the recently introduced Particles on Demand method. The multi-scale scheme is first validated through a convected athermal vortex and also studied in a jet flow setup. The performance of the proposed scheme is further investigated through the shock structure problem and a high Knudsen Couette flow, typical examples of highly non-equilibrium flows in which the order of the velocity set plays a decisive role. The results demonstrate that the proposed multi-scale scheme can operate accurately, with flexibility in terms of the underlying models and with reduced computational requirements.
We present a finite-difference integration algorithm for solution of a system of differential equations containing a diffusion equation with nonlinear terms. The approach is based on Crank-Nicolson method with predictor-corrector algorithm and provides high stability and precision. Using a specific example of short-pulse laser interaction with semiconductors, we give a detailed description of the method and apply it for the solution of the corresponding system of differential equations, one of which is a nonlinear diffusion equation. The calculated dynamics of the energy density and the number density of photoexcited free carriers upon the absorption of laser energy are presented for the irradiated thin silicon film. The energy conservation within 0.2% has been achieved for the time step $10^4$ times larger than that in case of the explicit scheme, for the chosen numerical setup. We also present a few examples of successful application of the method demonstrating its benefits for the theoretical studies of laser-matter interaction problems.
A monolithic coupling between the material point method (MPM) and the finite element method (FEM) is presented. The MPM formulation described is implicit, and the exchange of information between particles and background grid is minimized. The reduced information transfer from the particles to the grid improves the stability of the method. Once the residual is assembled, the system matrix is obtained by means of automatic differentiation. In such a way, no explicit computation is required and the implementation is considerably simplified. When MPM is coupled with FEM, the MPM background grid is attached to the FEM body and the coupling is monolithic. With this strategy, no MPM particle can penetrate a FEM element, and the need for computationally expensive contact search algorithms used by existing coupling procedures is eliminated. The coupled system can be assembled with a single assembly procedure carried out element by element in a FEM fashion. Numerical results are reported to display the performances and advantages of the methods here discussed.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
