Multiscale lattice Boltzmann approach to modeling gas flows

For multiscale gas flows, kinetic-continuum hybrid method is usually used to balance the computational accuracy and efficiency. However, the kinetic-continuum coupling is not straightforward since the coupled methods are based on different theoretical frameworks. In particular, it is not easy to recover the non-equilibrium information required by the kinetic method which is lost by the continuum model at the coupling interface. Therefore, we present a multiscale lattice Boltzmann (LB) method which deploys high-order LB models in highly rarefied flow regions and low-order ones in less rarefied regions. Since this multiscale approach is based on the same theoretical framework, the coupling precess becomes simple. The non-equilibrium information will not be lost at the interface as low-order LB models can also retain this information. The simulation results confirm that the present method can achieve model accuracy with reduced computational cost

The Moment Guided Monte Carlo method for the Boltzmann equation

In this work we propose a generalization of the Moment Guided Monte Carlo method developed in [11]. This approach permits to reduce the variance of the particle methods through a matching with a set of suitable macroscopic moment equations. In order to guarantee that the moment equations provide the correct solutions, they are coupled to the kinetic equation through a non equilibrium term. Here, at the contrary to the previous work in which we considered the simplified BGK operator, we deal with the full Boltzmann operator. Moreover, we introduce an hybrid setting which permits to entirely remove the resolution of the kinetic equation in the limit of infinite number of collisions and to consider only the solution of the compressible Euler equation. This modification additionally reduce the statistical error with respect to our previous work and permits to perform simulations of non equilibrium gases using only a few number of particles. We show at the end of the paper several numerical tests which prove the efficiency and the low level of numerical noise of the method.Comment: arXiv admin note: text overlap with arXiv:0908.026

Finite element modeling of non-equilibrium fluid-wall interaction beyond the continuum regime

The numerical modeling of the aerodynamic interactions at high-altitudes and high-Mach numbers is considered in view of its importance when studying problems where the continuum hypothesis at the foundation of the Navier-Stokes equations becomes invalid. One of the difficulties associated with these flight conditions is that both the velocity and the temperature of the fluid do not abide by the no - slip conditions at the wall. A weak-Galerkin Finite Element formulation of the Maxwell-Smoluchowsky model is introduced to discretize the velocity slip and temperature jump conditions with better accuracy than the standard Finite Element approximation. The methodology is assessed on configurations such as cylinders and spheres for flow conditions ranging from quasi-equilibrium to non-equilibrium. Improvements are observed in the slip regime compared to available data. Nonetheless, the results in the transition regime highlight the need for more sophisticated physical modeling to address non-equilibrium at the wall

Finite element modeling of nonequilibrium fluid-wall interaction at high-Mach regime

The numerical modeling of the aerodynamic interactions at high-Altitudes and high-Mach numbers is considered in view of its importance when studying problems where the continuum hypothesis at the foundation of the Navier- Stokes equations becomes invalid. One of the difficulties associated with these flight conditions is that both the velocity and the temperature of the fluid do not abide by the no-slip conditions at the wall. A weak Galerkin finite element formulation of the Maxwell-Smoluchowki model is introduced to discretize the velocity slip and temperature jump conditions with better accuracy than the standard finite element approximation. The methodology is assessed on configurations such as cylinders and spheres for flow conditions ranging from quasi-equilibrium to nonequilibrium. Improvements are observed in the slip regime compared with available data. Nonetheless, the results in the transition regime highlight the need for more sophisticated physical modeling to address nonequilibrium at the wall

Comparison of DAC and MONACO DSMC Codes with Flat Plate Simulation

Various implementations of the direct simulation Monte Carlo (DSMC) method exist in academia, government and industry. By comparing implementations, deficiencies and merits of each can be discovered. This document reports comparisons between DSMC Analysis Code (DAC) and MONACO. DAC is NASA's standard DSMC production code and MONACO is a research DSMC code developed in academia. These codes have various differences; in particular, they employ distinct computational grid definitions. In this study, DAC and MONACO are compared by having each simulate a blunted flat plate wind tunnel test, using an identical volume mesh. Simulation expense and DSMC metrics are compared. In addition, flow results are compared with available laboratory data. Overall, this study revealed that both codes, excluding grid adaptation, performed similarly. For parallel processing, DAC was generally more efficient. As expected, code accuracy was mainly dependent on physical models employed

Fluid Simulations with Localized Boltzmann Upscaling by Direct Simulation Monte-Carlo

In the present work, we present a novel numerical algorithm to couple the Direct Simulation Monte Carlo method (DSMC) for the solution of the Boltzmann equation with a finite volume like method for the solution of the Euler equations. Recently we presented in [14],[16],[17] different methodologies which permit to solve fluid dynamics problems with localized regions of departure from thermodynamical equilibrium. The methods rely on the introduction of buffer zones which realize a smooth transition between the kinetic and the fluid regions. In this paper we extend the idea of buffer zones and dynamic coupling to the case of the Monte Carlo methods. To facilitate the coupling and avoid the onset of spurious oscillations in the fluid regions which are consequences of the coupling with a stochastic numerical scheme, we use a new technique which permits to reduce the variance of the particle methods [11]. In addition, the use of this method permits to obtain estimations of the breakdowns of the fluid models less affected by fluctuations and consequently to reduce the kinetic regions and optimize the coupling. In the last part of the paper several numerical examples are presented to validate the method and measure its computational performances

A Hybrid Particle Scheme for Simulating Multiscale Gas Flows with Internal Energy Nonequilibrium

Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/83650/1/AIAA-2010-820-828.pd

Assessment of an All-Particle Hybrid Method for Hypersonic Rarefied Flow

Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/106450/1/AIAA2013-1203.pd

Extension of a Multiscale Particle Scheme to Near-Equilibrium Viscous Flows

Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/76880/1/AIAA-40262-641.pd

A fast spectral method for the Boltzmann equation for monatomic gas mixtures

Although the fast spectral method has been established for solving the Boltzmann equation for single-species monatomic gases, its extension to gas mixtures is not easy because of the non-unitary mass ratio between the di↵erent molecular species. The conventional spectral method can solve the Boltzmann collision operator for binary gas mixtures but with a computational cost of the order m3rN6, where mr is the mass ratio of the heavier to the lighter species, and N is the number of frequency nodes in each frequency direction. In this paper, we propose a fast spectral method for binary mixtures of monatomic gases that has a computational cost O(pmrM2N4 logN), where M2 is the number of discrete solid angles. The algorithm is validated by comparing numerical results with analytical Bobylev- Krook-Wu solutions for the spatially-homogeneous relaxation problem, for mr up to 36. In spatially-inhomogeneous problems, such as normal shock waves and planar Fourier/Couette flows, our results compare well with those of both the numerical kernel and the direct simulation Monte Carlo methods. As an application, a two-dimensional temperature-driven flow is investigated, for which other numerical methods find it difficult to resolve the flow field at large Knudsen numbers. The fast spectral method is accurate and elective in simulating highly rarefied gas flows, i.e. it captures the discontinuities and fine structures in the velocity distribution functions