• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.17 no.3
• /
• pp.209-219
• /
• 2013

The present paper deals with the extension of AUSMPW+ scheme into two-fluid model for multiphase flow. AUSMPW+ scheme is the improvement of a single-phase AUSM+ scheme by designing pressure-based weighting functions to prevent oscillations near a wall and shock instability after a strong shock. Recently, Kitamura and Liou assessed a family of AUSM-type schemes with two-fluid model governing equations [K. Kitamura and M.-S. Liou, Comparative study of AUSM-Family schemes in compressible multi-phase flow simulations, ICCFD7-3702 (2012)]. It was observed that the direct application of the single-phase AUSMPW+ did not provide satisfactory results for most of numerical test cases, which motivates the current study. It turns out that, by designing pressure-based weighting functions, which play a key role in controlling numerical diffusion for two-fluid model, problems reported in can be overcome. Various numerical experiments validate the proposed modification of AUSMPW+ scheme is accurate and robust to solve multiphase flow within the framework of two-fluid model.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.11 no.1
• /
• pp.73-80
• /
• 1999

This paper describes a numerical study of three-dimensional buoyant turbulent flow in a stairwell model with three convective differencing schemes, which include the upwind differencing scheme, the hybrid scheme and QUICK scheme. The Reynolds-averaged Navier-Stokes and energy equations are solved with a two-equation turbulence model. The Boussinesq approximation is used to model buoyancy terms in the governing equations. Three-dimensional predictions of the velocity and temperature fields are presented and are compared with experimental data. Three-dimensional simulations with each scheme have predicted the overall features of the flow fairly satisfactorily. A better agreement with experimental is achieved with QUICK scheme.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.6
• /
• pp.1640-1647
• /
• 1993

The numerical diffusion arising from streamline-to-grid skewness produces a deteriorating effect on the numerical accuracy. The QUICKER scheme to reduce the numerical diffusion requires more computational effort than the HYBRID scheme. This paper deals with the relative computational efficiencies of adopting QUICKER scheme with a coarser grid system and of adopting HYBRID scheme with a denser grid system. Laminar driven cavity flow with Re=400, 1000 is used as a test problem. It is found that QUICKER scheme with a coarser grid system is more efficient than the HYBRID scheme with a denser grid system.

• Journal of The Korean Astronomical Society
• /
• v.34 no.4
• /
• pp.271-273
• /
• 2001

We present a 4-parameter implicit Lagrangean code which satisfies conservation of mass, linear and angular momenta, energy and entropy simultaneously. The primary advantage of this scheme is possibility to control dissipative properties of the scheme avoiding the effects of numerical viscosity.

• Journal of computational fluids engineering
• /
• v.14 no.4
• /
• pp.13-22
• /
• 2009

A two-phase (gas and liquid) flow analysis solver, named CUPID, has been developed for a realistic simulation of transient two-phase flows in light water nuclear reactor components. In the CUPID solver, a two-fluid three-field model is adopted and the governing equations are solved on unstructured grids for flow analyses in complicated geometries. For the numerical solution scheme, the semi-implicit method of the RELAP5 code, which has been proved to be very stable and accurate for most practical applications of nuclear thermal hydraulics, was used with some modifications for an application to unstructured non-staggered grids. This paper is concerned with the effects of interpolation schemes on the simulation of two-phase flows. In order to stabilize a numerical solution and assure a high numerical accuracy, the second-order upwind scheme is implemented into the CUPID code in the present paper. Some numerical tests have been performed with the implemented scheme and the comparison results between the second-order and first-order upwind schemes are introduced in the present paper. The comparison results among the two interpolation schemes and either the exact solutions or the mesh convergence studies showed the reduced numerical diffusion with the second-order scheme.

• 한국방재학회:학술대회논문집
• /
• 2007.02a
• /
• pp.172-176
• /
• 2007

In this study, the numerical model for simulation of tsunamis is constructed by using the dispersion-correction scheme, 2nd upwind scheme, dynamic linking method, and so forth. The composed numerical model is used to simulate a hitorical tsunami event. The target tsunami event is the 1983 Central East Sea Tsunami. And, the predicted run-up heights of the tsunami at Imwon port are very reasonable compared to available observed data.

• Communications of the Korean Mathematical Society
• /
• v.22 no.3
• /
• pp.453-464
• /
• 2007

A computationally efficient numerical scheme is presented for the phase-field model of two-phase systems for anisotropic interfacial energy. The scheme is solved by using a nonlinear multigrid method. When the coefficient for the anisotropic interfacial energy is sufficiently high, the interface of the system shows corners or missing crystallographic orientations. Numerical simulations with high and low anisotropic coefficients show excellent agreement with exact equilibrium shapes. We also present spinodal decomposition, which shows the robustness of the pro-posed scheme.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.17 no.1
• /
• pp.29-45
• /
• 2013

Central schemes offer a simple and versatile approach for computing approximate solutions of nonlinear systems of hyperbolic conservation laws. However, there are large numerical dissipation in case of contact discontinuity. We study semi-discrete central upwind scheme by changing flux functions to reduce the numerical dissipation and we perform numerical computations for various problems in case of contact discontinuity.

• Communications of the Korean Mathematical Society
• /
• v.38 no.1
• /
• pp.299-318
• /
• 2023

This article devises an exponentially fitted method for the numerical solution of two parameter singularly perturbed parabolic boundary value problems. The proposed scheme is able to resolve the two lateral boundary layers of the solution. Error estimates show that the constructed scheme is parameter-uniformly convergent with a quadratic numerical rate of convergence. Some numerical test examples are taken from recently published articles to confirm the theoretical results and demonstrate a good performance of the current scheme.

• Water for future
• /
• v.28 no.4
• /
• pp.125-136
• /
• 1995

ELM, a characteristic line based method, was applied to advection-dispersion equation, and the results obtained were compared with those of Eulerian schemes(Stone-Brian and QUICKEST). The calculation methods consisted of Lagrangian interpolation scheme and cubic spline interpolation scheme for the advection calculation, and the Crank-Nicholson scheme for the dispersion calculation. The results of numerical methods were as follows: (1) for Gaussian hill: ELM, using Lagrangian interpolation scheme, gave the most accurate computational result, ELM, using cubic spline interpolation scheme, and QUICKEST scheme gave numerical damping for Peclet number 50. Stone-Brian scheme gave phase shift introduced in the numerical solution for Peclet number 10 and 50. (2) for advanced front: All schemes gave accurate computational results for Peclet number 1 and 4. ELM, Lagrangian interpolation scheme, and Stone,Brian scheme gave dissipation error and ELM, using cubic spline interpolation scheme, and QUICKEST scheme gave numerical oscillation for Peclet number 50.