• Journal of computational fluids engineering
• /
• v.15 no.2
• /
• pp.55-61
• /
• 2010

A dependence of weighting parameter in preconditioning method for solving low Mach number flow with incompressible flow nature is investigated. The present preconditioning method employs a finite-difference method applied Roe‘s flux difference splitting approximation with the MUSCL-TVD scheme and 4th-order Runge-Kutta method in curvilinear coordinates. From the computational results of benchmark flows through a 2-D backward-facing step duct it is confirmed that there exists a suitable value of the weighting parameter for accurate and stable computation. A useful method to determine the weighting parameter is introduced. With this method, high accuracy and stable computational results were obtained for the flow with low Mach number in the range of Mach number less than 0.3.

• Journal of computational fluids engineering
• /
• v.12 no.2
• /
• pp.53-62
• /
• 2007

A high resolution numerical method aimed at solving cavitating flow is proposed and applied to gas-liquid two-phase shock tube problem. The present method employs a finite-difference 4th-order Runge-Kutta method and Roe's flux difference splitting approximation with the MUSCL TVD scheme. By applying the homogeneous equilibrium cavitation model, the present density-based numerical method permits simple treatment of the whole gas-liquid two-phase flow field, including wave propagation and large density changes. The speed of sound for gas-liquid two-phase media is derived on the basis of thermodynamic relations and compared with that by eigenvalues. By this method, a Riemann problem for Euler equations of one dimensional shock tube was computed. Numerical results such as detailed observations of shock and expansion wave propagations through the gas-liquid two-phase media at isothermal condition and some data related to computational efficiency are made. Comparisons of predicted results and exact solutions are provided and discussed.

• 한국전산유체공학회:학술대회논문집
• /
• 2007.10a
• /
• pp.249-257
• /
• 2007

A high resolution numerical method aimed at solving gas-liquid two-phase flow is proposed and applied to gas-liquid two-phase shock tube problem. The present method employs a finite-difference 4th-order Runge-Kutta method and Roe's flux difference splitting approximation with the MUSCL TVD scheme. By applying the homogeneous equilibrium cavitation model, the present density-based numerical method permits simple treatment of the whole gas-liquid two-phase flow field, including wave propagation and large density changes. The speed of sound for gas-liquid two-phase media is derived on the basis of thermodynamic relations and compared with that by eigenvalues. By this method, a Riemann problem for Euler equations of one dimensional shock tube was computed. Numerical results 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 exact solutions are provided and discussed.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.38 no.12
• /
• pp.1153-1161
• /
• 2010

An efficient solution method for harmonic balance techniques with Fourier transform is presented for periodic unsteady flow problems. The present partially-implicit harmonic balance treats the flux terms implicitly and the harmonic source term is solved explicitly. The convergence of the partially Implicit method is much faster than the explicit Runge-Kutta harmonic balance method. The method does not need to compute the additional flux Jacobian matrix from the implicit harmonic source term. Compared with fully implicit harmonic balance method, this partial approach turns out to have good convergence property. Oscillating flows over NACA0012 airfoil are considered to verify the method and to compare with results of explicit R-K(Runge-Kutta) and dual time stepping methods.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• v.y2005m4
• /
• pp.343-346
• /
• 2005

Three dimensional structures of detonation wave propagating through a square-shaped duct were investigated using computational method and parallel processing. Inviscid fluid dynamics equations coupled with $variable-{\gamma}$ formulation and simplified one-step Arrhenius chemical reaction model were analysed by MUSCL-type TVD scheme and four stage Runge-Kutta time integration. The unsteady computational results in three dimension show the detailed mechanism of transverse mode and diagonal mode of detonation wave instabilities resulting same cell length but different cell width in smoked-foil record.

• Computational Structural Engineering
• /
• v.9 no.1
• /
• pp.133-139
• /
• 1996

In this paper, an approximate method is developed to obtain the natural frequencies of the out of plane vibration of circular curved beams. The governing differential equations are derived using the dynamic equilibrium equations with the Timoshenko theory, and solved numerically. The Runge-Kutta method and Regula-Falsi method are used to integrate the differential equations and to determine the natural frequencies, respectively. In numerical examples, the hinged-hinged and clamped-clamped end constraints are considered. For each case, the four lowest natural frequencies are reported as functions of four non-dimensional system parameters.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.34 no.5
• /
• pp.1383-1393
• /
• 2014

Recently, with rapid improvement in computer hardware and theoretical development in the field of computational fluid dynamics, high-order accurate schemes also have been applied in the realm of computational hydraulics. In this study, numerical solutions of 1D shallow water equations are presented with TVD Runge-Kutta discontinuous Galerkin (RKDG) finite element method. The transcritical flows such as dam-break flows due to instant dam failure and transcritical flow with bottom elevation change were studied. As a formulation of approximate Riemann solver, the local Lax-Friedrichs (LLF), Roe, HLL flux schemes were employed and MUSCL slope limiter was used to eliminate unnecessary numerical oscillations. The developed model was applied to 1D dam break and transcritical flow. The results were compared to the exact solutions and experimental data.

• 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.

• Proceedings of the Korean Nuclear Society Conference
• /
• 1997.05a
• /
• pp.96-101
• /
• 1997

양자가속기를 외부 중성자 공급원으로 하여 미임계 운전을 가능하게 하고 토륨을 핵연료로 사용하므로 장주기 핵종과 핵무기 재료물질의 발생량을 현저히 줄일 수 있는 새로운 노형인 energy amplifier에 대한 연구가 CERN을 중심으로 활발히 진행되고 있다. 본 연구에서는 토륨주기에 대하여 고정 중성자속 조사에 의한 핵분열 및 방사붕괴에 관한 모델을 정립하여 다수의 연립선형 미분방정식으로 구성하여 Runge Kutta 5-6차 자동시간 간격 수치해법을 이용하여 계산하였다. 결과는 1014의 고정 중성자속에 대하여 충분한 U233의 생산이 평형상태에 도달하고 장주기 핵종도 우라늄 주기에 비하여 현저히 줄어듬을 보이므로 가속기를 이용한 토륨 핵연료 주기의 타당성을 확인하였다.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.5 no.4
• /
• pp.266-273
• /
• 1981

갑자기 작용하는 외압을 받는 직교이방성원통쉘의 동적좌굴을 해석하였다. Donnell-Karnam 형의 비선형방정식을 유도하였으며 쉘의 초기불완전성도 고려하였다. Galerkin의 방법을 사용하여 운동방정식을 구하고 Runge-Kutta 수직해법으로 비선형방정식을 풀었다. 쉘의 직교이방성특성이 처짐-하중 관계식의 비선형성에 미치는 영향을 검토하였으며 동적산출하중의 판별법을 정의하였다. 본 연구의 결과, 직교이방성원통쉘은 쉘의 초기불완전성에 그리 민감하지 않음을 보여주었다.