Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userBurnett Spectral Method for High-Speed Rarefied Gas Flows
118
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We introduce a numerical solver for the spatially inhomogeneous Boltzmann equation using the Burnett spectral method. The modelling and discretization of the collision operator are based on the previous work [Z. Cai, Y. Fan, and Y. Wang, Burnett spectral method for the spatially homogeneous Boltzmann equation, arXiv:1810.07804], which is the hybridization of the BGK operator for higher moments and the quadratic collision operator for lower moments. To ensure the preservation of the equilibrium state, we introduce an additional term to the discrete collision operator, which equals zero when the number of degrees of freedom tends to infinity. Compared with the previous work [Z. Hu, Z. Cai, and Y. Wang,Numerical simulation of microflows using Hermite spectral methods, arXiv:1807.06236], the computational cost is reduced by one order. Numerical experiments such as shock structure calculation and Fourier flows are carried out to show the efficiency and accuracy of our numerical method.
rate research
Read More
A high-performance gas kinetic solver using multi-level parallelization is developed to enable pore-scale simulations of rarefied flows in porous media. The Boltzmann model equation is solved by the discrete velocity method with an iterative scheme. The multi-level MPI/OpenMP parallelization is implemented with the aim to efficiently utilise the computational resources to allow direct simulation of rarefied gas flows in porous media based on digital rock images for the first time. The multi-level parallel approach is analyzed in details confirming its better performance than the commonly-used MPI processing alone for an iterative scheme. With high communication efficiency and appropriate load balancing among CPU processes, parallel efficiency of 94% is achieved for 1536 cores in the 2D simulations, and 81% for 12288 cores in the 3D simulations. While decomposition in the spatial space does not affect the simulation results, one additional benefit of this approach is that the number of subdomains can be kept minimal to avoid deterioration of the convergence rate of the iteration process. This multi-level parallel approach can be readily extended to solve other Boltzmann model equations.
The reduced speed of sound technique (RSST) has been used for efficient simulation of low Mach number flows in solar and stellar convection zones. The basic RSST equations are hyperbolic, and are suitable for parallel computation by domain decomposition. The application of RSST is limited to cases where density perturbations are much smaller than the background density. In addition, non-conservative variables are required to be evolved using this method, which is not suitable in cases where discontinuities like shock waves co-exist in a single numerical domain. In this study, we suggest a new semi-conservative formulation of the RSST that can be applied to low Mach number flows with large density variations. We derive the wave speed of the original and newly suggested methods to clarify that these methods can reduce the speed of sound without affecting the entropy wave. The equations are implemented using the finite volume method. Several numerical tests are carried out to verify the suggested methods. The analysis and numerical results show that the original RSST is not applicable when mass density variations are large. In contrast, the newly suggested methods are found to be efficient in such cases. We also suggest variants of the RSST that conserve momentum in the machine precision. The newly suggested variants are formulated as semi-conservative equations, which reduce to the conservative form of the Euler equations when the speed of sound is not reduced. This property is advantageous when both high and low Mach number regions are included in the numerical domain. The newly suggested forms of RSST can be applied to a wider range of low Mach number flows.
The general synthetic iteration scheme (GSIS) is extended to find the steady-state solution of nonlinear gas kinetic equation, removing the long-standing problems of slow convergence and requirement of ultra-fine grids in near-continuum flows. The key ingredients of GSIS are that the gas kinetic equation and macroscopic synthetic equations are tightly coupled, and the constitutive relations in macroscopic synthetic equations explicitly contain Newtons law of shear stress and Fouriers law of heat conduction. The higher-order constitutive relations describing rarefaction effects are calculated from the velocity distribution function, however, their constructions are simpler than our previous work (Su et al. Journal of Computational Physics 407 (2020) 109245) for linearized gas kinetic equations. On the other hand, solutions of macroscopic synthetic equations are used to inform the evolution of gas kinetic equation at the next iteration step. A rigorous linear Fourier stability analysis in periodic system shows that the error decay rate of GSIS can be smaller than 0.5, which means that the deviation to steady-state solution can be reduced by 3 orders of magnitude in 10 iterations. Other important advantages of the GSIS are (i) it does not rely on the specific form of Boltzmann collision operator and (ii) it can be solved by sophisticated techniques in computational fluid dynamics, making it amenable to large scale engineering applications. In this paper, the efficiency and accuracy of GSIS is demonstrated by a number of canonical test cases in rarefied gas dynamics.
Fluid motion driven by thermal effects, such as that due to buoyancy in differentially heated three-dimensional (3D) enclosures, arise in several natural settings and engineering applications. It is represented by the solutions of the Navier-Stokes equations (NSE) in conjunction with the thermal energy transport equation represented as a convection-diffusion equation (CDE) for the temperature field. In this study, we develop new 3D lattice Boltzmann (LB) methods based on central moments and using multiple relaxation times for the three-dimensional, fifteen velocity (D3Q15) lattice, as well as it subset, i.e. the three-dimensional, seven velocity (D3Q7) lattice to solve the 3D CDE for the temperature field in a double distribution function framework. Their collision operators lead to a cascaded structure involving higher order terms resulting in improved stability. In this approach, the fluid motion is solved by another 3D cascaded LB model from prior work. Owing to the differences in the number of collision invariants to represent the dynamics of flow and the transport of the temperature field, the structure of the collision operator for the 3D cascaded LB formulation for the CDE is found to be markedly different from that for the NSE. The new 3D cascaded (LB) models for thermal convective flows are validated for natural convection of air driven thermally on two vertically opposite faces in a cubic cavity enclosure at different Rayleigh numbers against prior numerical benchmark solutions. Results show good quantitative agreement of the profiles of the flow and thermal fields, and the magnitudes of the peak convection velocities as well as the heat transfer rates given in terms of the Nusselt number.
The capability to simulate a two-way coupled interaction between a rarefied gas and an arbitrary-shaped colloidal particle is important for many practical applications, such as aerospace engineering, lung drug deliver and semiconductor manufacturing. By means of numerical simulations based on the Direct Simulation Monte Carlo (DSMC) method, we investigate the influence of the orientation of the particle and rarefaction on the drag and lift coefficients, in the case of prolate and oblate ellipsoidal particles immersed in a uniform ambient flow. This is done by modelling the solid particles using a cut-cell algorithm embedded within our DSMC solver. In this approach, the surface of the particle is described by its analytical expression and the microscopic gas-solid interactions are computed exactly using a ray-tracing technique. The measured drag and lift coefficients are used to extend the correlations available in the continuum regime to the rarefied regime, focusing on the transitional and free-molecular regimes. The functional forms for the correlations for the ellipsoidal particles are chosen as a generalisation from the spherical case. We show that the fits over the data from numerical simulations can be extended to regimes outside the simulated range of $Kn$ by testing the obtained predictive model on values of $Kn$ that where not included in the fitting process, allowing to achieve an higher precision when compared with existing predictive models from literature. Finally, we underline the importance of this work in providing new correlations for non-spherical particles that can be used for point-particle Euler-Lagrangian simulations to address the problem of contamination from finite-size particles in high-tech mechanical systems.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
