• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 1999년도 춘계 학술대회논문집
• /
• pp.73-78
• /
• 1999

극초음속 유동장의 정확한 해석을 위해 AUSMPW+ 수치기법과 충격파 포착시 생기는 수치오차를 제거하기 위해 충격파 정렬 기법(Shock-Aligned Grid Technique)을 개발하였다. AUSMPW+ 수치기법은 자체 수치점성이 적은 수치기법으로 점성 경계층 계산시 정확한 계산결과를 보여주며 기존의 AUSM 계열이 가지는 문제점인 물성치의 진동 현상을 제거한 수치기법이다. 원통형과 무딘 물체 주위의 극초음속 유동장 해석을 통해 공력이 진동현상 없이 정확하게 계산됨을 확인하였다. 그리고 충격파 정렬 기법의 특성을 파악하기 위해 충격파 반사문제와 충격파-충격파 상호작용 문제를 해석하여 수치오차 없이 충격파를 포착할 수 있음을 보였다. 또한 화학적 평형 비평형 유동 영역까지 충격파 정렬 격자 기법을 확장하였다.

• 대한기계학회논문집B
• /
• 제40권12호
• /
• pp.785-792
• /
• 2016

본 연구에서는 NACA0012 천음속 익형 유동에 있어 비평형 응축이 항력 발산 마하수 $M_D$에 미치는 영향을 TVD 유한 차분법을 사용하여 연구하였다. 받음각 ${\alpha}$가 동일한 경우, 정체점 상대습도 ${\phi}_0$가 높을수록 충격파 직전의 최대 마하수 $M_{max}$는 작게 되고 초음속 영역의 크기도 적게 된다. 주류 마하수 $M_{\infty}$, 정체점 상대습도 ${\phi}_0$ 및 정체온도 $T_0$가 동일한 경우 받음각 ${\alpha}$가 클수록 비평형 응축영역 길이 ${\Delta}_z$은 짧게 된다. 한편, 주류 마하수 $M_{\infty}$와 받음각 ${\alpha}$가 동일한 경우 정체점 상대습도 ${\phi}_0$가 높을수록 조파저항의 감소 때문에 항력계수 $C_D$는 적어진다. $M_D$는 ${\alpha}$가 동일한 경우 ${\phi}_0$가 클수록 크게 되며, ${\phi}_0$가 동일한 경우는 ${\alpha}$가 클수록 $M_D$는 적게 된다.

• 한국전산유체공학회지
• /
• 제5권2호
• /
• pp.47-54
• /
• 2000

A robust Navier-Stokes code has been developed to efficiently predict hypersonic flows in chemical nonequilibrium. The HLLE+ flux discretization scheme is used to improve accuracy and robustness of hypersonic flow analysis. An efficient LU approximate factorization method is also used to solve the flow equations and species continuity equations in fully coupled fashion to implicitly treat stiff source terms of chemical reactions. The HLLE+ scheme shows lower grid dependency for the wall heating rates than other schemes. The developed code has been used to compute chemical nonequilibrium air flow through expanding hypersonic nozzle and past two and three dimensional blunt-nosed bodies. The results are in good agreement with existing numerical and experimental results.

• 한국전산유체공학회지
• /
• 제20권3호
• /
• pp.1-7
• /
• 2015

When the steam is used as working fluid in fluid machinery, different from other gases as air, phase transition (steam condensation) can occur and it affects not only the flow fields, but also machine performance & efficiency. Therefore, considering phase transition phenomena in CFD calculation is required to achieve accurate prediction of steam flow and non-equilibrium wet-steam model is needed to simulate realistic steam condensing flow. In this research, non-equilibrium wet-steam model is implemented on in-house code(T-Flow), the flow fields including phase transition phenomena in convergent-divergent nozzle are studied and compared to results of advance researches.

• 대한기계학회논문집B
• /
• 제38권12호
• /
• pp.1009-1015
• /
• 2014

본 연구는 NACA0012 천음속 에어포일 유동에 있어서 비평형 응축이 Force 계수(압력, 양력 및 항력계수)에 미치는 영향을 TVD 수치해석을 통하여 연구하였다. 정체점 온도 298 K, 받음각 ${\alpha}=3^{\circ}$인 경우, 주류 마하수 0.78~0.81에서는 정체점 상대습도의 증가함에 따라 양력은 단순 감소한다. 반면 Lift force break 마하수 영역의 주류 마하수에서는 정체점 상대습도의 증가에 따라 양력은 오히려 증가한다. 받음 각 ${\alpha}=3^{\circ}$, 정체점 상대습도가 0%인 경우, 주류 마하수의 증가에 따라 항력은 급격하게 증가하지만, 응축의 영향이 큰 60%인 경우에는 주류 마하수의 증가에 조금 증가할 뿐이다. 동일한 주류 마하수인 경우 비평형 응축에 따른 전 항력의 감소는 받음각과 정체점 상대습도가 증가할수록 크게 된다. 응축이 없는 ${\Phi}_0=0%$인 경우는 주류 마하수가 크고 받음각이 클수록 Wave drag은 크게 되나 응축의 영향이 비교적 큰 ${\Phi}_0=50%$ 이상인 경우는 오히려 Wave drag이 작아지는 것으로 나타났다. 한편, 정체점 상대습도가 낮고, 주류 마하수가 클수록 충격파 직전의 최대 마하수는 커지는 것으로 나타났다.

• 대한기계학회논문집
• /
• 제13권5호
• /
• pp.999-1011
• /
• 1989

본 연구의 목적은 첫째, 재발달 경계층에서의 난류에서지 및 난류전단응력에 대한 전달방정식들의 각 항의 측정자료들을 보다 정확히 제시하고 항들간의 균형을 비교 평가함으로써 비평형 유동으로부터 평형유동으로 회복되는 과정을 검토하고, 둘째, 대표적인 난류 모델들로써 표존 k-.epsilon.모델 및 레이놀즈 응력 모델을 사용한 수치계산을 수행함으로써 이와같은 모델들이 비평형 유동을 서술함에 있어 발생될 수 있는 문제점들을 고찰하는데 있다.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2004년도 추계 학술대회논문집
• /
• pp.161-164
• /
• 2004

A finite volume numerical model is developed for simulating non-equilibrium electroosmotic flow in micro- and nanochannels. The Guoy-Chapman model is adopted to compute the flow and electric potential. The Nernst-Planck equation is employed to trace unsteady transports of ionic species, i.e., time-dependent net charge density. A new set of boundary conditions based on surface charge density are designed rather than using the conventionally-employed zeta potential. A few issues for an efficient computation of electroosmotic flows are discussed. Representative computational examples are given to illustrate the robustness of the numerical model.

• 대한기계학회논문집B
• /
• 제31권7호
• /
• pp.620-627
• /
• 2007

In this paper, we adapt a modified internal energy non-equilibrium first-order extrapolation thermal boundary condition to the thermal lattice Boltzmann model (TLBM). This model is the double populations approach to simulate hydrodynamic and thermal fields. The bounce-back boundary condition which is a traditional boundary condition of lattice Boltzmann method has only a first order in numerical accuracy at the boundary and numerical instability. A non-equilibrium first-order extrapolation boundary condition has been verified to be of better numerical stability than the bounce-back boundary condition and this boundary condition is proved to be of second-order accuracy for the flat boundaries. The two-dimensional natural convection flow in a square cavity with Pr=0.71 and various Rayleigh numbers are simulated. The results are found to be in good agreement with those of previous studies.

• 한국추진공학회:학술대회논문집
• /
• 한국추진공학회 1998년도 제11회 학술강연회논문집
• /
• pp.32-32
• /
• 1998

고체 추진제 로켓의 연소시에 발생되는 산화 알루미늄(A1$_2$O$_3$) 입자는 로켓 추진 노즐에서 팽창과정의 효율을 저하시키는 요소가 되며, 이러한 비효율성은 연소 가스와 입자간의 비평형 상태 효과와 기본적인 속도와 열적 차이에 의해서 발생된다고 보고되었다. 또한 연소시 발생된 산화 알루미늄 입자는 높은 열과 큰 운동량을 가지고 로켓 노즐 내부를 유동하게 되며, 매우 많은 량이 짧은 시간에 고온 고속으로 노즐 벽면이나 기타 구조물에 충돌 및 점착하기 때문에 로켓 노즐내의 표면이 손상을 입게 되고, 로켓의 방향 제어 및 조정 안정성이 저하되며, 구조적인 강도가 약화 될 수 있다. 또한 산화 알루미늄 액적들의 경우 노즐 벽면에 고착되게 되면 로켓의 중량 증가로 인해서 추력의 손실을 초래할 수 있다. 따라서 이러한 연소 부산물들의 운동 경로와 충돌 위치 및 표면에서의 충돌량과 그리고 충돌에 따른 마모량 및 점착 그리고 열전달 특성을 예측하는 것이 필수적이다.

• 대한기계학회논문집
• /
• 제18권6호
• /
• pp.1582-1589
• /
• 1994

The characteristics Prandt1-Meyer expansion of supersonic flow with condensation through a wavy wall in a channel are investigated by experiment and numerical direct marching method of characteristics. In the present study, for the case of moist air flow in the type of indraft supersonic wind tunnel, the dependency of location of formation and reflection of the oblique shock wave generated by the wavy wall and the distribution of flow properties, on the specific humidity and temperature at the entrance of wavy wall and the attack angle of the wavy wall to the main stream is clarified by schlieren photograph, distribution of static pressure and Mach number, and plots of numerical results. Also, we confirm that the wavy wall plays an important key role in the formation of oblique shock wave, and that the effect of condensation on the flow field appears apparently.