• 한국추진공학회:학술대회논문집
• /
• 한국추진공학회 2005년도 제24회 춘계학술대회논문집
• /
• pp.343-346
• /
• 2005

전산 해석과 병렬처리를 이용하여 정사각 관 내부를 전파하는 데토네이션파의 삼차원 구조를 살펴보았다. 가변 비열비 공식과 간략화 된 일단계 Arrhenius 반응 모델과 연계된 비점성 유체 방정식을 MUSCL 기반 TVD 해법과 4단계 Runge-Kutta 적분 방법을 이용하여 해석하였다. 삼차원에서의 비정상 해석 결과로부터 그을음 막 기록(smoked-foil record)에서 같은 길이와 다른 폭을 가지는 수평 및 대각 방향 불안정에 의한 상세한 파면 구조를 파악할 수 있었다.

• 한국추진공학회:학술대회논문집
• /
• 한국추진공학회 2008년도 제30회 춘계학술대회논문집
• /
• pp.332-335
• /
• 2008

본 연구의 목적은 수소 및 탄화수소를 이용한 bluff-body 화염 안정화 및 hyper-mixer와 같은 초음속 연소기 연료 분사기 주변의 유동해석 이다. 대규모 박리영역을 갖는 후방단 형상의 초음속 주변 유동을 해석하기 위한 검증 단계로 다양한 수치기법을 평가해 보았다. 고차 공간 정확도로 확장된 5차 TVD-MUSCL기법, Roe FDS, S-A DES/DDES 조합이 기저유동과 기저유출 유동에 좋은 결과를 보여 주었다.

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

A high resolution scheme for solving gas-liquid two-phase flows with cavitation is described. This scheme uses the curvilinear coordinate grid and solves the density based momentum equations for mixture of gas-liquid medium with a preconditioning method to treat both compressible and incompressible flow characteristics. The present preconditioned method is based on the Runge-Kutta explicit finite-difference scheme, and is improved by using the diagonalization, the flux difference splitting and the MUSCL-TVD schemes to save computational effort and to increase stability and resolvability, especially at gas-liquid contact surfaces. A homogeneous equilibrium cavitation model is used to treat the gas-liquid two-phase medium in cavitating flow as a locally homogeneous pseudo-single-phase medium. Therefore, it is easy to solve cavitating flow, including wave propagation, large density changes and incompressible flow characteristic at low Mach number. Some numerical results obtained by the present scheme are shown.

• 한국추진공학회지
• /
• 제10권2호
• /
• pp.1-14
• /
• 2006

본 논문은 약한 불안정 데토네이션 영역부터 강한 불안정 데토네이션 영역까지 여러 영역에 걸친 데토네이션 파 셀 구조 모사에 대한 수치적 문제점들을 살펴보았다. 비열 비 값이 변하는 점성 유체 역학 방정식 및 1단계 Arrhenius 반응 모델 해석을 위하여 MUSCL-type TVD 기법을 이용한 공간 차분과 4차 정확도의 Runge-Kutta 시간 적분을 이용하였다. 일련의 수치해석 연구는 여러 반응 상수 및 격자 해상도에 따른 데토네이션 셀 구조를 해석하기 위하여 요구되는 계산 조건을 구하기 위하여 다양한 데토네이션 현상 영역에서 수행되었다. 다른 영역의 데토네이션 현상에서 셀 구조를 포착하기 위한 계산 영역의 크기와 최소 격자 해상도를 찾아내기 위하여 정상 1차원 ZND 해석 결과와 전산 해석 결과를 비교 검토하였다.

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

In this paper, the preconditioned multistage time stepping methods which are popular multigrid smoothers is implemented for the compressible Navier-Stokes calculation with full-coarsening multigrid method. The convergence characteristic of the point-Jacobi and Alternating direction line Jacobi(DDADI) preconditioners are studied. The performance of 2nd order upwind numerical fluxes such as 2nd order upwind TVD scheme and MUSCL-type linear reconstruction scheme are compared in the inviscid and viscous turbulent flow caculations.

• 한국수자원학회논문집
• /
• 제42권12호
• /
• pp.1079-1089
• /
• 2009

본 연구에서는 마른하도 및 복잡한 지형에서의 파의 전파와 같은 수공학 분야에서 해결하기 어려운 문제를 해석하기 위한 고정확도 2차원 수치모형을 개발하기 위해, unsplit 유한체적기법과 HLLC Riemann 해법을 이용한 흐름율 계산으로 쌍곡선형 적분 보존형의 2차원 천수방정식을 해석하였다. Unsplit 기법의 적용을 위해 하상경사항은 발산정리를 이용하여 이산화한 형태를 적용하였으며, 흐름율과 생성항의 균형을 이루기 위해 수면경사법을 시간과 공간에 대해 2차정확도를 가지는 MUSCL 기법과 연계하였다. 그리고 적용한 생성항 처리기법과 흐름율과의 보존특성이 만족함을 보였다. 2차정확도의 사용으로 불연속 지점에서 발생할 수 있는 수치진동을 제거하기 위해서 경사제한자를 사용한 TVD 기법을 적용하였다. 개발모형을 정확해가 존재하는 생성항이 없는 1차원 댐 붕괴 흐름에 적용하여 흐름율 계산의 정확성을 검증하였고, 하상융기를 가진 하도의 정상류 및 천이류 모의를 통해 개발모형의 보존특성을 검증하였으며, 하상경사 및 단면의 확대/축소구간이 존재하는 2차원 댐 붕괴 흐름에 적용하여 개발모형의 적용성을 검증하였다.

• 한국수자원학회논문집
• /
• 제44권2호
• /
• pp.145-156
• /
• 2011

본 연구에서는 자연하천의 흐름에서 흔히 발생하는 천이류, 불연속류, 마른하도로의 파의 전파 등을 포함하는 복잡한 흐름을 해석하기 위한 고정확도 2차원 수치모형을 개발하였다. 하상경사항을 효율적으로 처리하기 위해 quasi-steady wave propagation 기법을 적용하여 해당 격자에 대한 생성항의 영향을 효율적으로 반영함으로써 쌍곡선형 적분 보존형의 2차원 천수방정식을 해석하였다. Fractional Step Method를 적용한 유한체적기법의 사용을 위해 HLL Riemann 해법을 이용하여 흐름률을 계산하였고, 시간 및 공간에 대한 2차 정확도를 만족하기 위해 MUSCL 기법을 적용하였다. 2차 정확도의 사용으로 불연속지점에서 발생하는 수치진동은 TVD 기법 적용을통해 제어하였다. 개발된모형은 2차원 제방 붕괴 및 댐하류부에 구조물이 존재하는 경우의댐 붕괴 모의를 통해실측치와의 검증을 실시하였다. 또한 하류부에 역경사가 존재하는 경우의 댐 붕괴 모의를 통해 실측치와 비교함으로써 생성항의 영향에 대한 모형의 적용성을 검증하였다.

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

In this paper, some numerical methods recently developed for gas-liquid two-phase flows are reviewed. And then, a preconditioning method to solve cavitating flow by the author is introduced. This method employs a finite-difference Runge-Kutta method combined with MUSCL TVD scheme, and a homogeneous equilibrium cavitation model. So that it permits to treat simply the whole gas-liquid two-phase flow field including wave propagation, large density changes and incompressible flow characteristic at low Mach number. Finally, numerical results such as detailed observations of the unsteady cavity flows, a sheet cavitation break-off phenomena and some data related to performance characteristics of hydrofoils are shown.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2003년도 The Fifth Asian Computational Fluid Dynamics Conference
• /
• pp.181-182
• /
• 2003

A preconditioned numerical method for gas-liquid to-phase flow is applied to solve cavitating flow. The present method employs a density based finite-difference method of dual time-stepping integration procedure and Roe's flux difference splitting approximation with MUSCL-TVD scheme. A homogeneous equilibrium cavitation model is used. The method permits simple treatment of the whole gas-liquid two-phase flow field including wave propagation, large density changes and incompressible flow characteristics at low Mach number. By this method, two-dimensional internal flows through a venturi tuve and decelerating cascades are computed and discussed.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2008년도 학술대회
• /
• pp.215-220
• /
• 2008

A numerical method for gas-liquid two-phase flow is applied to solve shock-bubble interaction problems. The present method employs a finite-difference Runge-Kutta method and Roe's flux difference splitting approximation with the MUSCL-TVD scheme. A homogeneous equilibrium cavitation model is used. By this method, a Riemann problem for shock tube was computed for validation. Then, shock-bubble interaction problems between cylindrical bubbles located in the liquid and incident liquid shock wave are computed.