DSMC 및 IP 방법을 이용한 2차원 아음속 희박기체 유동 해석 = Numerical Simulations of Two Dimensional Subsonic Rarefied Gas Flows using DSMC and IP Methods

다국어 초록 (Multilingual Abstract)

In the present study, to investigate the characteristics of subsonic rarefied gas flows efficiently, the information preservation(IP) method was applied to the direct simulation Monte-Carlo(DSMC) method to reduce statistical scatters from the DSMC method for low speed rarefied gas flows. To find the effects of the IP method, the results of the IP method were compared to the results of the DSMC method for the same numbers of sampling. Before the comparison, it was necessary to verify the DSMC solver. The verification of the solver was conducted by simulations of the micro-channel flows by comparing to the other researchers’ DSMC results. The velocity profiles at the 2/3 section of the channel agreed well with the other researchers’ results. In addition, the normalized slip velocity distributions on the wall and the pressure distributions along the centerline also agreed well with the other researchers’ results. After the DSMC solvers were verified, two types of Couette flows were considered. One is the flows that have temperature differences between two plates, and the other is the flows that have moving plates. For the first Couette flows, the temperature of 373K was used for the upper plate, and of 173K for lower plate. For the second Couette flows, the upper plate is moving at 300m/s, and the temperature of 273K was used for both two plates. All simulations of Couette flows were conducted from near-continuum to free-molecular regimes. From both Couette flow simulations, it was found that the temperature jump and the velocity slip occurred on the plate surface. In addition, it was also shown that the temperature and the velocity differences between the plate and the gas became larger as the flow fields became more rarefied. Lastly, flow simulations around a NACA0012 airfoil were conducted to identify the effects of the IP method. The freestream Mach number of 0.8 was used, and the Knudsen number of 0.014 was considered. It was observed that more clear contours were obtained from the IP method than the DSMC method, since the statistical scatter error was reduced by the IP method. In addition, it was found that the velocity slip on the airfoil surface occurred due to the effects of the rarefied atmospheric environment.

참고문헌 (Reference)

1 Koura, K, "Variable Soft Sphere Molecular Model for Inverse-power-law of Lenard Jones Potential" 3 : 2459-2465, 1991

2 Koura, K, "Variable Soft Sphere Molecular Model for Air Species" 4 : 1083-1085, 1992

3 Fan, J, "Statitical Simulation of Low-speed Unidirectional Flow in Transition Regime" 245-252, 1998

4 Wang, M, "Simulations for Gas Flows in Microgeometries Using the Direct Simulation Monte Carlo Method" 25 : 975-985, 2004

5 Masters, N.D, "Octant Flux Splitting Information Preservation DSMC Method for Thermally Driven Flows" 226 : 2044-2062, 2007

6 Bird, G.A, "Molecular Gas Dynamics and the Direct Simulation of Gas Flows" Clarendon 1994

7 Otten, D.L, "Inportance Sampling Based Direct Simulation Monte Carlo Method" (5061) : 1-13, 2010

8 Yao, Z.H, "IP-DSMC Method for Micro-scale Flow with Temperature Variation" 35 : 2016-2023, 2011

9 Roohi, E, "Extending the Navier-Stokes Solution to Transition Regime in Two-dimensional Micro- and Nanochannel Flows Using Information Preservation Scheme" 21 (21): 1-12, 2009

10 Fan, J, "Computation of Rarifes Gas Flows Around a NACA0012 Airfoil" 39 (39): 618-625, 2001

11 Sun, Q, "A Direct Simulation Method for Subsonic, Mircoscale Gas Flows" 178 : 400-425, 2002

12 Sun, Q, "A Direct Simulation Method for Subsonic, Mircoscale Gas Flows" University of Michigan 2003

13 Fei, F, "A Diffusive Information Preservation Method for Small Knudsen Number Flows" 243 : 179-193, 2013