Boltzmann type schemes for gas dynamics and the entropy property
SIAM Journal on Numerical Analysis
Rarefied flow computations using nonlinear model Boltzmann equations
Journal of Computational Physics
Coupling Boltzmann and Navier-Stokes equations by half fluxes
Journal of Computational Physics
A hybrid continuum/particle approach for modeling subsonic, rarefied gas flows
Journal of Computational Physics
A hybrid kinetic/fluid model for solving the gas dynamics Boltzmann-BGK equation
Journal of Computational Physics
A Smooth Transition Model between Kinetic and Diffusion Equations
SIAM Journal on Numerical Analysis
A smooth transition model between kinetic and hydrodynamic equations
Journal of Computational Physics
A hybrid particle-continuum method applied to shock waves
Journal of Computational Physics
Development and verification of a coupled DSMC-NS scheme using unstructured mesh
Journal of Computational Physics
A moving interface method for dynamic kinetic-fluid coupling
Journal of Computational Physics
Numerical Regularized Moment Method of Arbitrary Order for Boltzmann-BGK Equation
SIAM Journal on Scientific Computing
Hi-index | 31.46 |
We introduce a coupled method for hydrodynamic and kinetic equations on 2-dimensional h-adaptive meshes. We adopt the Euler equations with a fast kinetic solver in the region near thermodynamical equilibrium, while use the Boltzmann-BGK equation in kinetic regions where fluids are far from equilibrium. A buffer zone is created around the kinetic regions, on which a gradually varying numerical flux is adopted. Based on the property of a continuously discretized cut-off function which describes how the flux varies, the coupling will be conservative. In order for the conservative 2-dimensional specularly reflective boundary condition to be implemented conveniently, the discrete Maxwellian is approximated by a high order continuous formula with improved accuracy on a disc instead of on a square domain. The h-adaptive method can work smoothly with a time-split numerical scheme. Through h-adaptation, the cell number is greatly reduced. This method is particularly suitable for problems with hydrodynamics breakdown on only a small part of the whole domain, so that the total efficiency of the algorithm can be greatly improved. Three numerical examples are presented to validate the proposed method and demonstrate its efficiency.