• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.33 no.3
• /
• pp.164-169
• /
• 2009

The slip velocity and the temperature jumps for low-speed flow in microchannels are investigated using Langmuir slip boundary condition. This slip boundary condition is suggested to simulate micro flow. The current study analyzes Langmuir slip boundary condition theoretically and it analyzed numerically micro-Couette flow, micro-Poiseuille flow and grooved microchannel flow. First, to prove validity for Langmuir slip condition, an analytical solution for micro-Couette flow is derived from Navier-Stokes equations with Langmuir slip conditions and is compared with DSMC and an analytical solution with Maxwell slip boundary condition. Second, the numerical analysis is performed for micro-Poiseuille flow and grooved microchannel flow. The slip velocity and temperature distribution are compared with results of DSMC or Maxwell slip condition and those are shown in good agreement.

• 한국전산유체공학회:학술대회논문집
• /
• 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.

• 한국전산유체공학회:학술대회논문집
• /
• 2006.05a
• /
• pp.306-308
• /
• 2006

Numerical simulations were carried out of the effects of momentum and heat produced by a plasma actuator on neutral flow. Momentum and heat generated during plasma discharge were modeled as a body force and heat source using results of experiments and DSMC of particle. These force and heat model were inserted into a Navier-Stokes equation and the flow around the plasma actuator could be explored by solving fluid dynamics only. Fluid simulation showed that force produced in DSMC generated a jet flow in the vicinity of the plasma actuator and heat accounted for density change.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.28 no.12
• /
• pp.1494-1500
• /
• 2004

A zero-dimensional direct simulation Monte Carlo(DSMC) model is developed for simulating diatomic gas including vibrational kinetics. The method is applied to the simulation of two systems: vibrational relaxation of a simple harmonic oscillator and translational-rotational-vibrational energy exchange process under heating and cooling. 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.

• Transactions of the Korean Society of Machine Tool Engineers
• /
• v.17 no.2
• /
• pp.45-50
• /
• 2008

Molecular flows inside a guide block in the OLED(organic luminescent emitting device) deposition process have been simulated using DSMC(direct simulation Monte Carlo) method. Because the organic materials are evaporated under vacuum, molecules flow at a high Knudsen number of the free molecular regime, where the continuum mechanics is not valid. A guide block is designed as a part of the linear cell source to transport the evaporated materials to a deposition chamber, When solving the flows, the inlet boundary condition is proved to affect significantly the whole flow pattern. Thus, it is proposed that the pressure should be specified at the inlet. From the analysis of the density distributions at the nozzle exit of the guide block, it is shown that the longer nozzle can emit molecules more straightly. Finally, a nondimensionalized mass flow profile is obtained by numerical experiments, where various nozzle widths and inlet pressures are tested.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2008.03a
• /
• pp.602-605
• /
• 2008

The anode sheath structure in the hollow anode of an anode-layer type Hall thruster was numerically computed using a fully kinetic 2D3V Particle-in-Cell and Direct Simulation Monte Carlo(PIC-DSMC) code. By treating both ions and electrons as particles, anode surface region, which is electrically non-neutral, was analyzed. In order to analyze in detail, the calculation code was parallelized using Message Passing Interface (MPI). The code successfully simulated the discharge current oscillation. In the low magnetic induction case, ion sheath appears in the anode surface because ionization is enough to maintain the plasma occurs in the anode hollow. As the magnetic induction increases, main ionization region move to outside of the anode. At the same time, anode sheath voltage decreases. In the high magnetic induction case, electron sheath appears on the anode surface periodically because the ionization occurs mainly in the discharge channel. This anode sheath condition shift can be explained using the simple sheath model.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.12 no.1
• /
• pp.83-94
• /
• 2000

The performance of a turbomolecular pump(TMP) in both molecular and transition flow regions is predicted by the numerical solutions of the Boltzmann equation obtained by the direct simulation Monte Carlo method. The compression characteristics of the TMP are investigated for a wide range of the Knudsen number( Kn ). The maximum compression ratios strongly depend on Kn in transition region, while do they weakly on Kn in free molecular flow region. The present numerical results of the single blade row in both molecular and transition regions are used to predict the overall performance of a TMP, which has three kinds of blade with 24-rows.

• Journal of the KSME
• /
• v.44 no.3
• /
• pp.55-63
• /
• 2004

이 글에서는 마이크로와 나노스케일의 해석에 사용하는 수치모사 방법인 직접모사 몬테 카를로 (Direct Simulation Monte Carlo : DSMC)방법과 분자동역학(Molecular Dynamics: MD)과이 관계에 대하여 설명한다.

• Geophysics and Geophysical Exploration
• /
• v.15 no.2
• /
• pp.85-91
• /
• 2012

Multiple Lapse Time Window (MLTW) analysis for obtaining intrinsic attenuation value require numerous data without directional bias to compensate focal mechanism. The first window of MLTW, therefore, shows large deviation in fitting smoothed theoretical curve. The information on the focal mechanism may reduce burdens of number and distribution. This study combined algorithm of computing focal mechanism to DSMC method by Yoshimoto (2000). However, the MLTW method based on the numerous data was not applicable to this study, because of the limited data to the almost same fault plane solution. This study showed that the available data was too insufficient to construct smoothed theoretical curve, although the deviation of the first window was improved. Instead of conventional solution by more data, the study seems to be needed for new constraints to obtain smoothed curve.

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