• 한국전산유체공학회:학술대회논문집
• /
• 2006.10a
• /
• pp.36-40
• /
• 2006

A comparative study on flux functions for the 2-dimensional Euler equations has been conducted. Explicit 4-stage Runge-Kutta method is used to integrate the equations. Flux functions used in the study are Steger-Warming's, van Leer's. Godunov's, Osher's(physical order and natural order), Roe's, HILE, AUSM, AUSM+ and AUSMPW+. The performance of MUSCL limiters and MLP limiters in conjunction with flux functions are compared extensively for steady and unsteady problems.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.10a
• /
• pp.378-385
• /
• 2008

The present paper deals with the accurate and robust limiting procedure for the multi-dimensional flow analysis on unstructured meshes. The multi-dimensional limiting process (MLP) which was successfully proposed on structured grid system is extended to unstructured meshes. Based on MUSCL-type framework on unstructured meshes, the new slope limiter is devised to satisfy the MLP condition, which is quite effective to regulate the unwanted oscillations, especially on multiple dimensions. Considering the neighborhood based on the vertex of the cell, as well as the edge, this limiting strategy captures the multi-dimensional flow features very accurately with the proper stencils. From the various numerical results, these desirable characteristics of the proposed limiting strategy are clearly shown.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2010.11a
• /
• pp.760-763
• /
• 2010

The expansion of 1D numerical code to 2D or 3D is needed in order to improve the analysis accuracy of the interior ballistics. The cut cell method has been imposed for the code expansion to multi dimensions. The MUSCL-Hancock scheme as a high resolution method has been selected. A feasibility of the cut cell method has been verified by analyzing the free piston problem.

• Journal of computational fluids engineering
• /
• v.19 no.3
• /
• pp.91-97
• /
• 2014

A high resolution numerical method aimed at solving cavitating flow was proposed and applied to gas-liquid two-phase shock tube problem with arbitrary void fraction. The present method with compressibility effects employs a finite-difference 4th-order Runge-Kutta method and Roe's flux difference splitting approximation with the MUSCL TVD scheme. The Jacobian matrix from the inviscid flux of constitute equation is diagonalized analytically and the speed of sound for the two-phase media is derived by eigenvalues. So that the present method is appropriate for the extension of high order upwind schemes based on the characteristic theory. By this method, a Riemann problem for Euler equations of one dimensional shock tube was computed. Numerical results of high speed flow phenomena such as detailed observations of shock and expansion wave propagations through the gas-liquid two-phase media and some data related to computational efficiency are made. Comparisons of predicted results and solutions at isothermal condition are provided and discussed.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.27 no.6
• /
• pp.406-418
• /
• 2015

The present modified FUNWAVE-TVD model, which is a modification to its previous version 2.1, is applied to solitary wave propagation and is tested against the experiments of Vincent and Briggs(1989) and Luth et al.(1994). The eddy viscosity breaking scheme is used for comparison with the existing study in the case of breaking experiment. The symmetry of wave-induced current is maintained when the modified model is employed to Vincent and Briggs(1989) breaking experiment, but the symmetry of wave-induced current in previous model is not maintained. A better agreement with the breaking experimental data is obtained in the modified model using eddy viscosity breaking scheme than the shock capturing breaking scheme using nonlinear shallow water equation. For comparison with the schemes in the model, the fourth order MUSCL-TVD scheme by Erduran et al.(2005) and the third order MUSCL-TVD scheme using minmod limiter is applied, and the numerical solutions of solitary wave are compared.

• Journal of Korea Water Resources Association
• /
• v.42 no.12
• /
• pp.1079-1089
• /
• 2009

The objective of this study is to develop the accurate, robust and high resolution two-dimensional numerical model that solves the computationally difficult hydraulic problems, including the wave front propagation over dry bed and abrupt change in bathymetry. The developed model in this study solves the conservative form of the two-dimensional shallow water equations using an unsplit finite volume scheme and HLLC approximate Riemann solvers to compute the interface fluxes. Bed-slope term is discretized by the divergence theorem in the framework of FVM for application of unsplit scheme. Accurate and stable SGM, in conjunction with the MUSCL which is second-order-accurate both in space and time, is adopted to balance with fluxes and source terms. The exact C-property is shown to be satisfied for balancing the fluxes and the source terms. Since the spurious oscillations in second-order schemes are inherent, an efficient slope limiting technique is used to supply TVD property. The accuracy, conservation property and application of developed model are verified by comparing numerical solution with analytical solution and experimental data through the simulations of one-dimensional dam break flow without bed slope, steady transcritical flow over a hump and two-dimensional dam break flow with a constriction.

• 한국방재학회:학술대회논문집
• /
• 2009.02b
• /
• pp.152.1-152.1
• /
• 2009

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

It is well known that preconditioning methods are efficient for convergence acceleration in the compressible low Mach number flows. In this study, the original Euler equations and three differently nondimensionalized preconditioning methods are implemented in two dimensional inviscid bump flows using the 3rd order MUSCL and DADI schemes as numerical flux discretization and time integration, respectively. The multigrid and local time stepping methods are also used to accelerate the convergence. The test case indicates that a properly modified local preconditioning technique involving concepts of a global preconditioning allows Mach number independent convergence. Besides, an asymptotic analysis for properties of preconditioning methods is added.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.10 no.2
• /
• pp.1-14
• /
• 2006

Present study examines the numerical issues of cell structure simulation for various regimes of detonation phenomena ranging from weakly unstable to highly unstable detonations. Inviscid fluid dynamics equations with $variable-{\gamma} $ formulation and one-step Arrhenius reaction model are solved by a MUSCL-type TVD scheme and 4th order accurate Runge-Kutta time integration scheme. A series of numerical studies are carried out for the different regimes of the detonation phenomena to investigate the computational requirements for the simulation of the detonation wave cell structure by varying the reaction constants and grid resolutions. The computational results are investigated by comparing the solution of steady ZND structure to draw out the minimum grid resolutions and the size of the computational domain for the capturing cell structures of the different regimes of the detonation phenomena.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2012.05a
• /
• pp.581-581
• /
• 2012

경사면의 잠김과 드러남(wetting and drying)에 대한 천수방정식의 수치 실험으로서 그 전선 근처에서 수심이 매우 작은 경우, 수심이 작을 때 2차 정확도를 위한 MUSCL의 적용 여부, 그리고 마찰 적용의 한계 수심 등에 대해 수심경사법(DGM)과 수면경사법(SGM)을 비교 검토하였다. 먼저, 포물선형 융기에서 정상상태의 기준해와 비교하고 하류에서 배수(drain)가 되는 경우를 상정하였다. 그리고 삼각형 바닥 턱을 지나는 댐 붕괴 흐름에 대한 실험의 결과와 비교하였다. 마지막으로 유사난수 발생에 의한 불규칙 바닥 구간을 지나는 가상 댐 붕괴 흐름에 대해 모의하였다. 정상 상태의 잠긴 바닥(wet bed)에 대해 정수 상태에서 SGM의 선평형성(well-balancedness)이 만족되었으나, 흐름이 상류와 사류일 때 각각 DGM과 SGM에서 유량 변동이 있음을 확인하였다. 드러난 바닥(dry bed)으로 흐름이 진행될 때 수심이 음의 값이 되지 않도록 잠김 드러남 전선근처에서 수심이 매우 작아지는 경우($<1{\times}10^{-6}m{\sim}1{\times}10^{-5}m$) 수심을 0으로 강제하고, n이 0.011과 0.04인 경우 수심이 각각 $1{\times}10^{-4}m$와 $1{\times}10^{-2}m$보다 작은 경우 마찰의 적용을 배제하는 것이 불가피하였다. SGM에 의한 MUSCL에서 수심이 $5{\times}10^{-5}m{\sim}5{\times}10^{-2}m$보다 작을 때 1차 정확도로 회귀가 불가피하였으나, DGM에서는 그와 같은 정확도의 훼손이 나타나지 않았다.