• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.30 no.5
• /
• pp.32-40
• /
• 2002

The study of nonlinear gas transport in rarefied condition or associated with the microscale length of the geometry has emerged as an interesting topic in recent years. Along with the DSMC method, several fluid dynamic models that come under the general category of the moment method or the Chapman-Enskog method have been used for this type of problem. In the present study, on the basis of Eu's generalized hydrodynamics, computational models for diatomic gases are developed. The rotational nonequilibrium effect is included by introducing excess normal stress associated with the bulk viscosity of the gas. The new models are applied to study the one-dimensional shock structure and the multi-dimensional rarefied hypersonic flow about a blunt body. The results indicate that the bulk viscosity plays a considerable role in fundamental flow problems such as the shock structure and shear flow. An excellent agreement with experiment is observed for the inverse shock density thickness.

• 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.

• 유체기계공업학회:학술대회논문집
• /
• 2006.08a
• /
• pp.493-496
• /
• 2006

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.

• Journal of Mechanical Science and Technology
• /
• v.19 no.2
• /
• pp.664-673
• /
• 2005

The heat transfer characteristics of rarefied flows in a micro-actuator are studied numerically. The effect of Knudsen number (Kn) on the heat transfer of the micro-actuator flows is also examined. The Kn based on gas density and characteristic dimension is varied from near-continuum to highly rarefied conditions. Direct simulation Monte Carlo calculations have been performed to estimate the performance of the micro-actuator. The results show that the magnitude of the temperature jump at the wall increases with Kn. Also, the heat transfer to the isothermal wall is found to increase significantly with Kn.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.4
• /
• pp.509-517
• /
• 1997

Results of flat plate compressible boundary layer calculation, based on discrete formulation of DSMC method, are presented in low Mach number and low Knudsen number range. The free stream is a uniform flow of pure nitrogen at various Mach numbers in low pressures (i.e. rarefied gas). Complete thermal accommodation and diffuse molecular reflections are used as the wall boundary condition, replacing unreal no-slip condition used in continuum calculations. In the discrete formulation of DSMC method, there is no need to use ad hoc assumptions on transport properties like viscosity and thermal conductivity, instead viscosity is calculated from values of other field variables (velocity and shear stress). Also the results are compared with existing self-similar continuum solutions. In all Mach number cases computed, velocity slip is most pronounced in regions near the leading edge where continuum formulation renders the solution singular. As the boundary layer develops further downstream, velocity slips asymptote to values that are between 10 to 20% of the magnitude of free stream velocity. When the free stream number density is reduced, so the gas more rarefied, the velocity slip increases as expected.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.31 no.2
• /
• pp.1-9
• /
• 2003

In the present study, a 3-D Parallel DSMC method in developed on unstructured meshes for the efficient simulation of rarefied gas flows. Particle tracing between cells in achieved based on a linear shape function extended to three dimensions. For high parallel efficiency, successive domain decomposition is applied to achieve load balancing between processors by accounting for the number of particles. A particle weighting technique is also adopted to handle flows containing gases of significantly dirrerent number densities in the same flow domain. Application is made for flow past a 3-D delta wing and the result is compared with that from experiment and other calculation. Flow around a rocket payload at 100km altitude is also solved and the effect of plume back flow from the nozzle in studied.

• Aerospace Engineering and Technology
• /
• v.8 no.1
• /
• pp.144-153
• /
• 2009

To study the plume effects in the rarefied region, the Direct Simulation Monte Carlo(DSMC) method is usually adopted because the plume field usually contains the entire range of flow regime from the near-continuum in the vicinity of nozzle exit through transitional state to free molecular at far field region from the nozzle. The objective of this study is to investigate the behaviors of a small monopropellant thruster plume in the rarefied region numerically using DSMC method. To deduce accurate results efficiently, the preconditioned scheme is introduced to calculate continuum flow fields inside thruster to predict nozzle exit properties used for inlet conditions of DSMC method. By combining these two methods, the rarefied flow characteristics of plume such as strong nonequilibrium near nozzle exit, large back flow region, etc, can be investigated.

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

The direct simulation Monte Carlo(DSMC) method is employed to calculate the etch rate on Al wafer. The etchant is assumed to be Cl$_2$. The etching process of an Al wafer in a helicon plasma etcher is examined by simulating molecular collisions of reactant and product. The flow field inside a plasma etch reactor is also simulated by the DSMC method fur a chlorine feed gas flow. The surface reaction on the Al wafer is simply modelled by one-step reaction: 3C1$_2$+2Allongrightarrow1 2AIC1$_3$. The gas flow inside the reactor is compared for six different nozzle locations. It is found that the flow field inside the reactor is affected by the nozzle locations. The Cl$_2$ number density on the wafer decreases as the nozzle location moves toward the side of the reactor. Also, the present numerical results show that the nozzle location 1, which is at the top of the reactor chamber, produces a higher etch rate.

• Journal of computational fluids engineering
• /
• v.1 no.1
• /
• pp.142-149
• /
• 1996

이차원 노즐을 통하여 저밀도 환경으로 팽창하는 희박류의 분석에 있어서 불연속좌표법과 결합된 유한차분법(finite-difference method coupled with the discrete-ordinate method, FDDO)과 직접모사법(direct-simulation Monte-Carlo method, DSMC)이 비교되었다. FDDO를 이용한 분석에서는 충돌적분모델을 도입하여 간단해진 볼츠만식(Boltzmann equation)이 불연속좌표법을 이용하여 물리적 공간에서는 연속이나 분자속도 공간에서는 불연속좌표로 표시되는 편미분방정식군으로 변환되어 유한차분법에의하여 수치해석 되었다. 직접모사법에서는 분자모델로 가변강구모델(variable hard sphere model, VHS)이, 충돌샘플링모델로는 비시계수법(no time counter method, NTC)이 채택되었다. 전혀 다른 두 가지 방법에 의한 노즐 내부에서의 유체흐름 해석결과는 매우 잘 일치하였으며, 노즐 외부의 plume 영역에서는 FDDO에 의한 해석결과가 직접모사법에 의한 해석결과에 비하여 약간 느린 팽창을 보였다.

• Advances in aircraft and spacecraft science
• /
• v.2 no.1
• /
• pp.45-56
• /
• 2015

The attitude aerodynamic control is an important subject in the design of an aerospace plane. Usually, at high altitudes, this control is fulfilled by thrusters so that the implementation of an aerodynamic control of the vehicle has the advantage of reducing the amount of thrusters fuel to be loaded on board. In the present paper, the efficiency of a wing-flap has been evaluated considering a NACA 0010 airfoil with a trailing edge flap of length equal to 35% of the chord. Computational tests have been carried out in hypersonic, rarefied flow by a direct simulation Monte Carlo code at the altitudes of 65 and 85 km, in the range of angle of attack 0-40 deg. and with flap deflection equal to 0, 15 and 30 deg.. Effects of the flap deflection have been quantified by the variations of the aerodynamic force and of the longitudinal moment. The shock wave-boundary layer interaction and the shock wave-shock wave interaction have been also considered. A possible interaction of the leading edge shock wave and of the shock wave arising from the vertex of the convex corner, produced on the lower surface of the airfoil when the flap is deflected, generates a shock wave whose intensity is stronger than those of the two interacting shock waves. This produces a consistent increment of pressure and heat flux on the lower surface of the flap, where a thermal protection system is required.