• 한국전산유체공학회:학술대회논문집
• /
• 1999.11a
• /
• pp.42-47
• /
• 1999

In this study, we adopt the point symmetric Gauss-Seidel relaxation algorithm to obtain the steady state solution of the Navier-Stokes equations for the thermal and chemical nonequilibrium flow of air. All of the inviscid, viscous flux Jacobians and thermochemical source Jacobians are included in the implicit part Numerical simulation is performed for the thermal and chemical nonequilibrium flow over blunt body and computational results are presented. The convergence history and CPU time of the present computation are compared with the LU-SGS scheme which employs the approximate Jacobians.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 1997.11a
• /
• pp.17-17
• /
• 1997

극초음속 여객기와 군사용 항공기에 대한 수요가 증가함에 따라서 새로운 개념의 다양한 추진기관이 연구가 진행되고 개발되어 왔다. 초음속 항공기의 속도 영역은 마하 10-20 정도가 되는데 이 속도 한계를 극복하기 위하여 초음속 연소 램제트 엔진(SCRamjet; Supersonic Combustion Ramjet)이 제안되었다. 스크램 제트를 개발하기 위해서는 연료와 산화제의 혼합 효율 문제, 화염의 안정화 문제, 벽면의 냉각에 관한 문제 등 몇 가지 기본적인 문제들을 해결해야 한다. Univ of Michigan에서 실험한 연소기를 모델로 본 연구에서는 연료와 공기의 혼합에 관한 수치 연구를 수행하였다. 다원 혼합기체에 관한 축대칭 Navier-Stokes 방정식을 지배 방정식을 이용하였고 비평형 화학반응식을 고려하였다. 공간 차분에는 유한 체적법을 이용하였다. 대류 플럭스 항은 Roe의 Upwind FDS 기법을 사용하여 차분하였고 점성항에는 중심 차분법을 이용하였다. 시간 적분법으로는 근사 자코비안과 LU분할 기법을 이용한 완전 내재적 방법이 쓰였다. 난류 모델로는 Mentor에 의해 제안된 2 방정식 k-$\varepsilon$/k-$\omega$ 혼합모델을 사용하였다. 유동장이 실험에서의 찍은 사진과 유사한 모습의 충격파 간섭을 수치 모사하였고 수소가 확산되는 모습과 함께 노즐 lip 주위의 재순환 영역에 대해서 살펴볼 수 있었다.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2000.11a
• /
• pp.20-20
• /
• 2000

2개의 경사 충격파와 하나의 수직 충격파로 초음속 유동을 압축하는 초음속 공기 흡입구의 수치적 연구를 수행하였다. 지배방정식으로는 Navier-Stokes방정식을 사용하였고 난류모델로는 SST 모델을 사용하였다. 지배방정식의 점성항 계산에는 중심차분법을 사용하였고 대류항 계산에는 풍상차분법인 Roe의 FDS기법을 MUSCL기법과 결합하여 이용하였다. 유한 체적법을 이용하여 차분된 방정식은 LU분할 기법을 이용한 완전 내재적 방법으로 2차 정확도 시간 적분으로 비정상 과정의 연구를 수행하였다. 흡입구 배압을 정해주어야 하는 어려움을 해결하기 위해 흡입구 후면에 노즐을 달고 노즐의 면적을 조절하여 배압이 형성되도록 하였다.(중략)

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2000.04a
• /
• pp.16-16
• /
• 2000

충격파를 포함하는 초음속 유동을 해석하는 수치해법 중에서 많이 사용되어진 것은 엄밀 및 근사 리만 해법과 플럭스 분할 기법들로서 이들은 Euler 방정식에 기반을 두고 선형 또는 비선형파의 상호작용을 풍상 차분법으로 기술하는 방법들이다. 이러한 수치기법들은 과거 광범위하게 사용되어 왔으나 최근 여러 가지 단점이 발견되었다. 이와 같은 문제점을 극복하고자 입자의 통계적인 운동을 기술하는 기체 운동론에 근거하여 BGK 수치기법이 제시되었다. 이는 비충돌 볼츠만 방정식으로부터 입자의 수준에서 플럭스 분할 기법 형태의 풍상차분법을 구현하는 것으로 볼츠만 방정식의 충돌항을 BGK 모델로 대치하고 이것의 적분해로부터 수치 플럭스를 구한다. 이 수치기법은 기존의 리만해법에 비하여 수치적으로나 물리적으로 매우 타당한 성질인 강건성, 정확성, 엔트로피 조건, 양수보존성 등을 가지고 있음이 밝혀졌다. 이와 같은 수치기법을 사용하여 로켓 노즐 내의 아음속, 천이음속, 초음속에서의 유동장 해석을 위한 프로그램을 작성하였다. 시간 적분에 대하여는 정상 상태의 계산을 위하여 내재적 시간 적분 방법을 사용하였으며, 공간 이산화 방법으로는 임의의 제어체적에 대하여 적분형 보존 방정식을 적용하는 유한 체적법을 사용하였다. 초음속 입구 유동과 출구에서 초음속과 저음속 유동의 두가지 경우를 고려하여 얻은 결과를 기존의 연구 결과와 비교하여 본 결과 잘 일치하였다. 입구 유동이 저음속이고 출구 유동이 초음속인 경우에 대하여도 해석결과가 실험결과와 잘 일치하였다. 상대적으로 낮은 온도, 압력 조건과 높은 온도, 압력 조건을 가지는 고체 로켓 모터 노즐 내의 유동을 해석하였다. 이들 해석 결과를 전압, 전온도로 표준화시킨 결과 서로 일치하였으며, 파라서 저온, 저압에서 얻은 결과도 표준화시킬 경우, 고온, 고압에서도 사용될 수 있음을 알 수 있었다.의 영향에 초점을 맞추었다.다고 판단되며 배기 가스 자체에 대기 공기중에 함유되어 있던 습기가 얼어붙는(Icing화) 문제가 발생하기 때문에 배기가스의 Icing을 방지하기 위하여 압축기 끝단에서 공기를 추출하여 배기부분에 송출할 필요성이 있는 것으로 판단되었다. 출구가스의 기체 유동속도가 매우 빠르므로 (100-l10m.sec) 이를 완화하기 위한 디퓨저의 설계가 요구된다고 판단된다. 또 연소기 후방에 물을 주입하는 경우 열교환기 및 기타 부분품에 발생할 수 있는 부식 및 열교환 효율 저하도 간과할 수 없는 문제로 파악되었다. 이러한 기술적 문제가 적절히 해결되는 경우 비활성 가스 제너레이터는 민수용으로는 대형 빌딩, 산림, 유조선 등의 화재에 매우 적절히 사용되어 질 수 있을 뿐 아니라 군사적으로도 군사작전 중 및 공군 기지의 화재 그리고 지하벙커에 설치되어 있는 고급 첨단 군사 장비 등의 화재 뿐 아니라 대간첩작전 등에 효과적으로 활용될 수 있을 것으로 판단된다.가 작으며, 본 연소관에 충전된 RDX/AP계 추진제의 경우 추진제의 습기투과에 의한 추진제 물성 변화는 미미한 것으로 나타났다.의 향상으로, 음성개선에 효과적이라고 사료되었으며, 이 방법이 편측 성대마비 환자의 효과적인 음성개선의 치료방법의 하나로 응용될 수 있으리라 생각된다..7%), 혈액투석, 식도부분절제술 및 위루술·위회장문합술을 시행한 경우가 각 1례(2.9%)씩이었다. 13) 심각한 합병증은 9례(26.5%)에서 보였는데 그중 식도협착증이 6례(17.6%), 급성신부전증 1례(2.9%), 종격동기흉과 폐염이 병발한 경우와 폐염이 각 1례(2.9%)였다. 14) 식도경 시행회수는 1회가 17례(54.8%), 2회가 9례(29.0%), 3회 이상이 5례(16.1%)였다.EX>$IC_{50}$/ 값이 210 $\mu\textrm{g}$<

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.30 no.8
• /
• pp.37-45
• /
• 2002

A two-dimensional unsteady, inviscid flow solver has been developed for the simulation of airfoil-vortex interactions on unstructured dynamically adapted meshes. The Euler solver is based on a second-order accurate implicit time integration using a point Gauss-Seidel relaxation scheme and a dual time-step subiteration. A vertex-centered, finite-volume discretization is used in conjunction with the Roe's flux-difference splitting. An unsteady solution-adaptive dynamic mesh scheme is used by adding and deleting mesh points to take account of both spatial and temporal variations of the flow field. The effect of vortex interaction on the aeroelastic displacement of an airfoil attached to the idealized two degree-of-freedom spring system is investigated.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.16 no.12 s.103
• /
• pp.1246-1254
• /
• 2005

Theoretical investigation is made on the electromagnetic fields generated by an impulsive point current source, fur the VV, HV, and VH problems at the interface between an isotropic upper half-space medium and a normally uniaxial lower half-space medium. The electric fields of these problems are associated only with the extraordinary-wave components in the Fourier-Laplace domain. Applying the Cagniard-de Hoop method to each problem, the time-domain solutions of the wave fields are obtained. The fields of the VV case can be expressed in explicit(integral-free) forms. The fields of the HV and VH cases are not integral-free, but the major singularities implicit in the integral solutions can be analytically extracted. The interfacial far fields in the frequency domain are determined by the singularities in the time domain.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2003.05a
• /
• pp.221-224
• /
• 2003

A numerical study is carried out for the detonation wave propagation through a T-branch. The T-branch is a crucial part of the combustion wave igniter, a novel concept of rocket ignition system aimed for the simultaneous ignition of multiple combustion chambers by delivering detonation waves. Euler equation and induction parameter equation are used as governing equations with a reaction term modeled from the chemical kinetics database obtained from a detailed chemistry mechanism. Second-order accurate implicit time integration and third-order space accurate TVD algorithm were used for solution of the coupled equations. Over two-million grid points enabled the capture of the dynamics of the detonation wave propagation including the degeneration and re-initiation phenomena, and some of the design factors were be obtained for the CWI flame tubes.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.40 no.5
• /
• pp.405-412
• /
• 2012

A flow control on airfoil installed a vortex cell for high efficiency wind turbine blade in stationary and dynamic stall conditions have been numerically investigated by solving the compressible Navier-Stokes equations. The numerical scheme is based on a node-based finite-volume method with Roe's flux-difference splitting and an implicit time-integration method coupled with dual time step sub-iteration. The computed result for the airfoil in the stationary showed that lift-drag ratio increases due to low pressure by the vortex cell. The oscillating airfoil with the vortex cell showed that the magnitude of hysteresis loop is reduced due to the enhanced vortex in the cell.

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

A numerical simulation of an incompressible cavity flow is conducted. Two dimensional Navier-Stokes equations are integrated using staggered grid and a finite volume method with C-QUICK scheme for spatial derivatives and fully implicit scheme for the time derivatives. SIMPLE-C algorithm is employed to solve the pressure field. Computational results show that the third eddy is generated in the shear layer mode but not in the steady mode. This signifies that the third eddy plays an important role in cavity flow stability. As a means to control the flow, jet blowing is applied to a position below the cavity upstream edge. Effects of flow control parameters on the stability such as the frequency, the phase, and the velocity magnitude are reported.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.33 no.7
• /
• pp.1-8
• /
• 2005

For study on the unsteady wall interference effect, flows around a circular cylinder in closed test-section wind tunnels have been numerically investigated by solving compressible Navier-Stokes equations. The numerical scheme is based on a node-based finite-volume method with the Roe's flux-difference splitting and an implicit time-integration method coupled with dual time-step sub-iteration. The computed results showed that the unsteady pressure gradient over the cylinder is enhanced by the wall interference, and as a result the fluctuations of lift and drag are augmented. The drag is further increased because of the lower base pressure. The vortex shedding frequency is also increased by the wall interference. The pressure on the test section wall shows the harmonics having the shedding frequency contained in the wall effect.