• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 1997.04a
• /
• pp.608-614
• /
• 1997

2차원 압축성 오일러 방정식을 지배 방정식으로 하여 음향학적 문제를 직접 수치 계산하였다. Roe의 개략적인 리만해(Approximate Riemann Solver)에 의거하여 대류항을 구성하였고 높은 정확도를 유지하기 위해 Vanleer의 MUSCL 방법을 사용하여 공간을 이산화하였다. 시간 적분으로는 2차의 Runge-Kutta 방법을 사용하였다. Resonance tube 내의 음향학적 문제 해석을 위해 최근까지 제시된 음향 경계 조건들을 비교하였고 각 경계 조건들이 갖는 물리적 의미를 고찰하였다. 이들을 통하여 압력파의 물리적인 행태를 모사하였고 실험과 비교하였다. 계산결과들을 통하여 사용된 방법은 전형적인 압축성 유동 해석 문제 뿐만 아니라 선형, 비선형의 음향학적 문제들을 적절한 경계 조건과 결합하여 해석할 수 있음을 알 수 있었고 resonance tube와 같은 공학적으로 실제적인 문제들에 대한 응용이 기대된다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2007.05a
• /
• pp.428-432
• /
• 2007

댐 제방 등의 붕괴로 인하여 발생하는 급격한 유량의 변화와 흐름영역의 변화로 인한 천이류 및 도수의 발생, 불규칙한 하천단면에서 갈수기 저수기의 흐름해석은 기존의 수치해법의 한계로 인하여 수리모형실험 및 경험식 또는 단면의 단순화 등에 의존하고 있는 실정이다. 본 연구에서는 자연하천에서 비선형 흐름율 계산에 불연속초기조건의 해석해인 Riemann 근사해법을 사용하여 수치적으로 안정되고 정확한 1차원 모형을 개발하고자 한다. 이를 위하여 유한체적법을 사용하였고, 수위와 유량의 계산을 위하여 요구되는 유한체적을 유출입하는 흐름율의 계산에 HLL Riemann 해법을 사용하였으며, MUSCL 기법으로 2차 정확도기법으로 확장하였다. Riemann 해법을 통하여 계산된 비선형의 흐름율과 보존 특성을 만족시켜줄 수 있는 하상 및 하폭변화로 인한 생성항을 처리하는 기법을 제안함으로서 새로운 1차원 수치해석모형을 개발하였다. 개발된 모형의 실제하천의 적용성을 확인하기 위하여 하상과 하폭이 변화하는 부정류 흐름에 적용하여 모형의 적용성 및 정확성을 검증하였다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2007.05a
• /
• pp.1332-1336
• /
• 2007

댐 제방 붕괴파는 갑작스러운 유량의 증가가 발생하여 불연속적인 흐름특성을 가지는 충격파(shock wave)가 전파되며, 갈수기 저수기에는 중소하천의 상류, 여울과 소에서의 흐름 또는 낙차공이나 보, 댐 여수로 등의 수공구조물에서 부분적인 사류 흐름이 발생된다. 이 때 흐름은 한계수위를 통과하게 되므로 기존 수치해법의 적용에 어려움이 존재한다. 본 연구에서는 실제하천에 적용될 수 있는 1차원 HLL, Roe Riemann 근사해법들을 간단히 소개하고, 시간공간적으로 2차의 고정확도 기법으로 확장하는 방법에 대하여 소개하였다. 각 기법을 정확해가 존재하는 댐붕괴 및 마른하도 전파의 경우에 적용하여 각 기법의 적용성 및 정확성을 비교하였다. 그리고 기존 Lax-Friedrichs 기법과 Lax-Wendroff 기법의 적용결과를 비교하였다. 적용결과 Lax-Friedrichs 기법을 제외한 모든 기법이 정확해와 잘 일치하였으며 특히 HLL 기법을 2차 정확도로 확장한 WAF 기법이 가장 높은 정확도로 계산되었다. 그러나 비선형 생성항이 존재하는 실제하천에 있어서 MUSCL 기법을 이용한 2차 정확도 기법이 합리적일 것으로 판단된다.

• 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
• /
• 2008.03a
• /
• pp.597-601
• /
• 2008

Three-dimensional structures of detonation wave propagating in circular tube were investigated. Inviscid fluid dynamics equations coupled with a conservation equation of reaction progress variable were analyzed by a MUSCL-type TVD scheme and four stage Runge-Kutta time integration. Variable-$\gamma$ formulation was used to account for the variable properties between unburned and burned states and the chemical reaction was modeled by using a simplified one-step irreversible kinetics model. The computational code was parallelized based on domain decomposition technique using MPI-II message passing library. The computations were carried out using a home made Windows based PC cluster having 160 AMD AthloxXP and Athlon64 processor. The computational domain consisted of through a roundshaped tube with wall conditions. As an initial condition, analytical ZND solution was distributed over the computational domain with disturbances. The disturbances has circumferential large gradient. The unsteady computational results in three-dimension show the detailed mechanisms of multi-cell mode of detonation wave instabilities resulting diamond shape in smoked-foil record.

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

The low-speed preconditioning technique is applied to solve the compressible Reynolds averaged Navier-Stokes equations for low-speed flows. The space discretization is based on Roe's flux-difference splitting with third-order-accurate MUSCL extrapolation. Time integration is performed employing a diagonal approximate factorization algorithm. The dual-time stepping has been incorporated to solve the unsteady flows. Full multigrid method is implemented to accelerate the convergence rate. To verify the algorithms several cases have been tested. Demonstrated the improvement on convergence and quality of the solution.

• Journal of Korea Water Resources Association
• /
• v.44 no.7
• /
• pp.511-522
• /
• 2011

Two-dimensional shallow water model based on the cut cell and the adaptive mesh refinement techniques is presented in this paper. These two mesh generation methods are combined to facilitate modeling of complex geometries. By using dynamically adaptive mesh, the model can achieve high resolution efficiently at the interface where flow changes rapidly. The HLLC Reimann solver and the MUSCL method are employed to calculate advection fluxes with numerical stability and precision. The model was applied to simulate the extreme urban flooding experiments performed by the IMPACT (Investigation of Extreme Flood Processes and Uncertainty) project. Simulation results were in good agreement with observed data, and transient flows as well as the impact of building structures on flood waves were calculated with accuracy. The cut cell method eased the model sensitivity to refinement. It can be concluded that the model is applicable to the urban flood simulation in case the effects of sewer and stormwater drainage system on flooding are relatively small like the dam brake.

• Journal of Korea Water Resources Association
• /
• v.41 no.8
• /
• pp.761-772
• /
• 2008

The object of this study is to develop the model that solves the numerically difficult problems in hydraulic engineering and to demonstrate the applicability of this model by means of various test examples, such as, verification in the gradually varied unsteady condition, three steady flow problems with the change of bottom slope with exact solution, and frictional bed with analytical solution. The governing equation of this model is the integral form of the Saint-Venant equation satisfying the conservation laws, and finite volume method with the Riemann solver is used. The evaluation of the mass and momentum flux with the HLL Riemann approximate solver is executed. MUSCL-Hancock scheme is used to achieve the second order accuracy in space and time. This study introduce the new and simple technique to discretize the source terms of gravity and hydrostatic pressure force due to longitudinal width variation for the balance of quantity between nonlinear flux and source terms. The results show that the developed model's implementation is accurate, robust and highly stable in various flow conditions with source terms, and this model is reliable for one-dimensional applications in hydraulic engineering.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.29 no.5B
• /
• pp.429-439
• /
• 2009

When Catastrophic extreme flood occurs due to dam break, the response time for flood warning is much shorter than for natural floods. Numerical models can be powerful tools to predict behaviors in flood wave propagation and to provide the information about the flooded area, wave front arrival time and water depth and so on. But flood wave propagation due to dam break can be a process of difficult mathematical characterization since the flood wave includes discontinuous flow and dry bed propagation. Nevertheless, a lot of numerical models using finite volume method have been recently developed to simulate flood inundation due to dam break. As Finite volume methods are based on the integral form of the conservation equations, finite volume model can easily capture discontinuous flows and shock wave. In this study the numerical model using Riemann approximate solvers and finite volume method applied to the conservative form for two-dimensional shallow water equation was developed. The MUSCL scheme with surface gradient method for reconstruction of conservation variables in continuity and momentum equations is used in the predictor-corrector procedure and the scheme is second order accurate both in space and time. The developed finite volume model is applied to 2D partial dam break flows and dam break flows with triangular bump and validated by comparing numerical solution with laboratory measurements data and other researcher's data.

• Journal of Mechanical Science and Technology
• /
• v.16 no.11
• /
• pp.1511-1521
• /
• 2002

High-speed moist air or steam flow has long been of important subject in engineering and industrial applications. Of many complicated gas dynamics problems involved in moist air flows, the most challenging task is to understand the nonequilibrium condensation phenomenon when the moist air rapidly expands through a flow device. Many theoretical and experimental studies using supersonic wind tunnels have devoted to the understanding of the nonequilibrium condensation flow physics so far. However, the nonequilibrium condensation can be also generated in the subsonic flows induced by the unsteady expansion waves in shock tube. The major flow physics of the nonequilibrium condensation in this application may be different from those obtained in the supersonic wind tunnels. In the current study, the nonequilibrium condensation phenomenon caused by the unsteady expansion waves in a shock tube is analyzed by using the two-dimensional, unsteady, Navier-Stokes equations, which are fully coupled with a droplet growth equation. The third-order TVD MUSCL scheme is applied to solve the governing equation systems. The computational results are compared with the previous experimental data. The time-dependent behavior of nonequilibrium condensation of moist air in shock tube is investigated in details. The results show that the major characteristics of the nonequilibrium condensation phenomenon in shock tube are very different from those in the supersonic wind tunnels.