• Journal of computational fluids engineering
• /
• v.4 no.3
• /
• pp.63-70
• /
• 1999

A Direct simulation Monte-Carlo (DSMC) code is developed, which employs the Monte-Carlo statistical sampling technique to investigate hypersonic rarefied gas flows accompanying chemical reactions. The DSMC method is a numerical simulation technique for analyzing the Boltzmann equation by modeling a real gas flow using a representative set of molecules. Due to the limitations in computational requirements. the present method is applied to a flow around a simple two-dimensional object in exit velocity of 7.6 km/sec at an altitude of 90 km. For the calculation of chemical reactions an air model with five species (O₂, N₂, O, N, NO) and 19 chemical reactions is employed. The simulated result showed various rarefaction effects in the hypersonic flow with chemical reactions.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.03a
• /
• pp.84-91
• /
• 2008

Generalized hydrodynamic (GH) theory for multi-species gas and the computational models are developed for the numerical simulation of hypersonic rarefied gas flow on the basis of Eu's GH theory. The rotational non-equilibrium effect of diatomic molecules is taken into account by introducing excess normal stress associated with the bulk viscosity. The numerical model for the diatomic GH theory is developed and tested. Moreover, with the experience of developing the dia-tomic GH computational model, the GH theory is extended to a multi-species gas including 5 species; O$_2$, N$_2$, NO, O, N. The multi-species GH model includes diffusion relation due to the molecular collision and thermal phenomena. Two kinds of GH models are developed for an axisymmetric flow solver. By compar-ing the computed results of diatomic and multi-species GH theories with those of the Navier-Stokes equations and the DSMC results, the accuracy and physical consistency of the GH computational models are examined.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.10a
• /
• pp.84-91
• /
• 2008

Generalized hydrodynamic (GH) theory for multi-species gas and the computational models are developed for the numerical simulation of hypersonic rarefied gas flow on the basis of Eu's GH theory. The rotational non-equilibrium effect of diatomic molecules is taken into account by introducing excess normal stress associated with the bulk viscosity. The numerical model for the diatomic GH theory is developed and tested. Moreover, with the experience of developing the dia-tomic GH computational model, the GH theory is extended to a multi-species gas including 5 species; $O_2,\;N_2$, NO, O, N. The multi-species GH model includes diffusion relation due to the molecular collision and thermal phenomena. Two kinds of GH models are developed for an axisymmetric flow solver. By compar-ing the computed results of diatomic and multi-species GH theories with those of the Navier-Stokes equations and the DSMC results, the accuracy and physical consistency of the GH computational models are examined.

• Journal of computational fluids engineering
• /
• v.19 no.2
• /
• pp.40-48
• /
• 2014

For rarefied gas flow regimes, physical phenomena such as velocity slip and temperature jump occur on the solid body surface. To predict these phenomena accurately, either the Navier-Stokes solver with a slip boundary condition or the direct simulation Monte Carlo method should be used. In the present study, flow simulations of a wedge were conducted in Mach-10 flow of argon gas for several different flow regimes using a two-dimensional Navier-Stokes solver with the Maxwell slip boundary condition. The results of the simulations were compared with those of the direct simulation Monte Carlo method to assess the present method. It was found that the values of the velocity slip and the temperature jump predicted increase as the Knudsen number increases. Also, the results are comparatively reasonable up to the Knudsen number of 0.05.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.26 no.10
• /
• pp.1394-1405
• /
• 2002

The performance of micro-actuators utilizing radiometric forces are studied numerically. The Knudsen number based on gas density and characteristic dimension is varied from near-continuum to highly rarefied conditions. Direct simulation Monte Carlo(DSMC) calculations have been performed to estimate the performance of the micro-actuators. In the present DSMC method, the variable hard sphere molecular model and no time counter technique are used to simulate the molecular collision kinetics. For simulation of diatomic gas flows, the Borgnakke-Larsen phenomenological model is adopted to redistribute the translational and internal energies.

• Proceedings of the KSME Conference
• /
• 2000.04b
• /
• pp.625-630
• /
• 2000

Numerical and experimental investigations are peformed for the rarefied gas flows in pumping channels of a helical-type drag pump. Modern turbomolecular pumps include a drag stage in the discharge side, operating roughly in $10^{-2}{\sim}10Torr$. The flow occurring in the pumping channel develops from the molecular transition to slip flow traveling downstream. Two different numerical methods are used in this analysis: the first one is a continuum approach in solving the Navier-Stokes equations with slip boundary conditions, and the second one is a stochastic particle approach through the use of the direct simulation Monte Carlo(DSMC) method. The flow in a pumping channel is three-dimensional(3D), and the main difficulty in modeling a 3D case comes from the rotating frame of reference. Thus, trajectories of particles are no longer straight lines. In the Present DSMC method, trajectories of particles are calculated by integrating a system of differential equations including the Coriolis and centrifugal forces. Our study is the first instance to analyze the rarefied gas flows in rotating frame in the presence of noninertial effects.

• 한국전산유체공학회:학술대회논문집
• /
• 1998.05a
• /
• pp.138-144
• /
• 1998

The BGK equation, which is the kinetic model equation of Boltzmann equation, is solved using FDDO(finite difference with the discrete-ordinate method) to compute the rarefied flow of monatomic gas. Using reduced velocity distribution and discrete ordinate method, the scalar equation is transformed into a system of hyperbolic equations. High resolution ENO(Essentially Non-Oscillatory) scheme based on Harten-Yee's MFA(Modified Flux Approach) method with Strang-type explicit time integration is applied to solve the system equations. The calculated results are well compared with the experimental density field of NACA0012 airfoil, validating the developed computer code. Next. the computed results of circular cylinder flow for various Knudsen numbers are compared with the DSMC(Direct Simulation Monte Carlo) results by Vogenitz et al. The present scheme is found to be useful and efficient far the analysis of two-dimensional rarefied gas flows, especially in the transitional flow regime, when compared with the DSMC method.

• Journal of computational fluids engineering
• /
• v.9 no.3
• /
• pp.1-7
• /
• 2004

A kinetic theory analysis is made of low-speed rarefied gas flows around a flat plate. The Boltzmann equation simplified by a collision model is solved by means of a finite difference approximation with the discrete ordinate method. The method does not suffer from statistical noise which is common in particle based methods and requires much less amount of computational effort. Calculations are made for flows around a micro-scale flat plate with a finite length of 20 microns. The method is assessed by comparing the results with those from several different methods and available experimental data.

• The KSFM Journal of Fluid Machinery
• /
• v.10 no.5
• /
• pp.27-33
• /
• 2007

Rarefied gas flows through two-dimensional micro channels are studied numerically for the performance optimization of a nanomembrane-based Knudsen compressor. The effects of the wall temperature distributions on the thermal transpiration flow patterns are examined. The flow has a pumping effect, and the mass flow rates through the channel are calculated. The results show that a steady one-way flow is induced for a wide range of the Knudsen number. The DSMC(direct simulation Monte Carlo) method with VHS(variable hard sphere) model and NTC(no time counter) techniques has been applied in this work to obtain numerical solutions. A critical element that drives Knudsen compressor Is the thermal transpiration membrane. The membranes are based on aerosol or machined aerogel. The aerogel is modeled as a single micro flow channel.

• Journal of computational fluids engineering
• /
• v.19 no.3
• /
• pp.77-84
• /
• 2014

The lattice Boltzmann (LB) method has been used to simulate rarefied gas flows in a micro-system as an alternative tool. However, previous results were mainly focused on a simple geometry with flat walls because the LB method is modeled on uniform Cartesian lattices. When previous boundary conditions for the microflows are applied to curved walls, the use of them requires approximation of the curved boundary by a series of stair steps, and introduces additional errors. For macroflows, no-slip curved wall boundary treatments have been developed remarkably in order to overcome these limits. However, the investigations for the slip curved wall boundary have rarely been performed for microflows. In this work, a curved boundary treatment of the LB method for a slip flow has been introduced. The results of the LB method for 2D microchannel and 3D microtube flows are in excellent agreement with the analytical solutions.