Mesoscopic simulation of diffusive contaminant spreading in gas flows at low pressure

Many modern production and measurement facilities incorporate multiphase systems at low pressures. In this region of flows at small, non-zero Knudsen- and low Mach numbers the classical mesoscopic Monte Carlo methods become increasingly numerically costly. To increase the numerical efficiency of simulations hybrid models are promising. In this contribution, we propose a novel efficient simulation approach for the simulation of two phase flows with a large concentration imbalance in a low pressure environment in the low intermediate Knudsen regime. Our hybrid model comprises a lattice-Boltzmann method corrected for the lower intermediate Kn regime proposed by Zhang et al. for the simulation of an ambient flow field. A coupled event-driven Monte-Carlo-style Boltzmann solver is employed to describe particles of a second species of low concentration. In order to evaluate the model, standard diffusivity and diffusion advection systems are considered.Comment: 9 pages, 8 figure

Micro-scale cavities in the slip - and transition - flow regimes

Differences between Navier-Stokes-Fourier (NSF) slip/jump solutions and direct simulation Monte-Carlo (DSMC) computations are highlighted for a micro lid-driven cavity problem. The results indicate a need for better modelling techniques which at the same time retain low computational cost of NSF models. We also highlight the fact thatmany micro-flows that have been considered are simple planar flows and typical classification systems are defined on such flows. We show that for complex flows, such as thedriven cavity, non-equilibrium effects are more appreciable and their onset occurs at lower Knudsen numbers than expected

Effects of rarefaction on cavity flow in the slip regime

The Navier-Stokes-Fourier equations, with boundary conditions that account for the effects of velocity-slip and temperature-jump, are compared to the direct simulation Monte Carlo method for the case of a lid-driven micro-cavity. Results are presented for Knudsen numbers within the slip-flow regime where the onset of nonequilibrium effects are usually observed. Good agreement is found in predicting the general features of the velocity field and the recirculating flow. However, although the steady-state pressure distributions along the walls of the driven cavity are generally in good agreement with the Monte Carlo data, there is some indication that the results are starting to show noticeable differences, particularly at the separation and reattachment points. The modified Navier-Stokes-Fourier equations consistently overpredict the maximum and minimum pressure values throughout the slip regime. This highlights the need for alternative boundary formulations or modeling techniques that can provide accurate and computationally economic solutions over a wider range of Knudsen numbers

A Multiscale Kinetic-Fluid Solver with Dynamic Localization of Kinetic Effects

This paper collects the efforts done in our previous works [P. Degond, S. Jin, L. Mieussens, A Smooth Transition Between Kinetic and Hydrodynamic Equations, J. Comp. Phys., 209 (2005) 665--694.],[P.Degond, G. Dimarco, L. Mieussens, A Moving Interface Method for Dynamic Kinetic-fluid Coupling, J. Comp. Phys., Vol. 227, pp. 1176-1208, (2007).],[P. Degond, J.G. Liu, L. Mieussens, Macroscopic Fluid Model with Localized Kinetic Upscaling Effects, SIAM Multi. Model. Sim. 5(3), 940--979 (2006)] to build a robust multiscale kinetic-fluid solver. Our scope is to efficiently solve fluid dynamic problems which present non equilibrium localized regions that can move, merge, appear or disappear in time. The main ingredients of the present work are the followings ones: a fluid model is solved in the whole domain together with a localized kinetic upscaling term that corrects the fluid model wherever it is necessary; this multiscale description of the flow is obtained by using a micro-macro decomposition of the distribution function [P. Degond, J.G. Liu, L. Mieussens, Macroscopic Fluid Model with Localized Kinetic Upscaling Effects, SIAM Multi. Model. Sim. 5(3), 940--979 (2006)]; the dynamic transition between fluid and kinetic descriptions is obtained by using a time and space dependent transition function; to efficiently define the breakdown conditions of fluid models we propose a new criterion based on the distribution function itself. Several numerical examples are presented to validate the method and measure its computational efficiency.Comment: 34 page

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

Models for Metal Hydride Particle Shape, Packing, and Heat Transfer

A multiphysics modeling approach for heat conduction in metal hydride powders is presented, including particle shape distribution, size distribution, granular packing structure, and effective thermal conductivity. A statistical geometric model is presented that replicates features of particle size and shape distributions observed experimentally that result from cyclic hydride decreptitation. The quasi-static dense packing of a sample set of these particles is simulated via energy-based structural optimization methods. These particles jam (i.e., solidify) at a density (solid volume fraction) of 0.665+/-0.015 - higher than prior experimental estimates. Effective thermal conductivity of the jammed system is simulated and found to follow the behavior predicted by granular effective medium theory. Finally, a theory is presented that links the properties of bi-porous cohesive powders to the present systems based on recent experimental observations of jammed packings of fine powder. This theory produces quantitative experimental agreement with metal hydride powders of various compositions.Comment: 12 pages, 12 figures, 2 table

An immersed boundary method for the fluid--structure--thermal interaction in rarefied gas flow

An immersed boundary method for the fluid--structure--thermal interaction in rarefied gas flow is presented. In this method, the slip model is incorporated with the penalty immersed boundary method to address the velocity and temperature jump conditions at the fluid--structure interface in rarefied gas flow within slip regime. In this method, the compressible flow governed by Navier-Stokes equations are solved by using high-order finite difference method; the elastic solid is solved by using finite element method; the fluid and solid are solved independently and the fluid--structure--thermal interaction are achieved by using a penalty method in a partitioned way. Several validations are conducted including Poiseuille flow in a 2D pipe, flow around a 2D NACA airfoil, moving square cylinder in a 2D pipe, flow around a sphere and moving sphere in quiescent flow. The numerical results from present method show good agreement with the previous published data obtained by other methods, and it confirms the the good ability of the proposed method in handling fluid--structure--thermal interaction for both weakly compressible and highly compressible rarefied gas flow. To overcome the incapability of Navier-Stokes equations at high local Knudsen numbers in supersonic flow, an artificial viscosity is introduced to ease the sharp transition at the shock wave front. Inspired by Martian exploration, the application of proposed method to study the aerodynamics of flapping wing in rarefied gas flow is conducted in both 2D and 3D domains, to obtain some insights for the flapping-wing aerial vehicles operating in Martian environment