• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.27 no.7
• /
• pp.995-1006
• /
• 2003

The suitability of high-order accurate, centered and upwind-biased compact difference schemes for large eddy simulation is evaluated by a dynamic analysis. Large eddy simulation of isotropic turbulence is performed with various dissipative and non-dissipative schemes to investigate the effect of numerical dissipation on the resolved solutions. It is shown by the present dynamic analysis that upwind schemes reduce the aliasing error and increase the finite differencing error. The existence of optimal upwind scheme that minimizes total numerical error is verified. It is also shown that the finite differencing error from numerical dissipation is the leading source of numerical errors by upwind schemes. Simulations of a turbulent channel flow are conducted to show the existence of the optimal upwind scheme.

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

RoeM and AUSMPW+ schemes are two of the most accurate and efficient schemes which are recently developed for the analysis of single phase gas dynamics. In this paper, we developed two-phase versions of these schemes for the analysis of gas-liquid large density ratio two-phase flow. We adopt homogeneous equilibrium model (HEM) using mass fraction to describe different two phases. In the Eulerian-Eulerian framework, HEM assumes dynamic and thermal equilibrium of the two phases in the same computational mesh. From the mixture equation of state (EOS), we derived new shock-discontinuity sensing term (SDST), which is commonly used in RoeM and AUSMPW+ for the stable numerical flux calculation. The proposed two-phase versions of RoeM and AUSMPW+ schemes are applied on several air-water two-phase test problems. In spite of the large discrepancy of material properties such as density, enthalpy, and speed of sound, the numerical results show that both schemes provide very satisfactory solutions.

• Journal of computational fluids engineering
• /
• v.11 no.1 s.32
• /
• pp.36-42
• /
• 2006

The existing numerical approximations of convection flux, especially the spatial higher-order difference schemes, in unstructured cell-centered finite volume methods are examined in detail with each other and evaluated with respect to the accuracy through their application to a 2-D benchmark problem. Six higher-order schemes are examined, which include two second-order upwind schemes, two central difference schemes and two hybrid schemes. It is found that the 2nd-order upwind scheme by Mathur and Murthy(1997) and the central difference scheme by Demirdzic and Muzaferija(1995) have more accurate prediction performance than the other higher-order schemes used in unstructured cell-centered finite volume methods.

• Journal of Ocean Engineering and Technology
• /
• v.12 no.2 s.28
• /
• pp.97-110
• /
• 1998

The numerical damping and dispersion error characteristics associated with difference schemes and a panel shift method used for the calculation of steady free surface flows by a panel method are an analysed in this paper. First, 12 finite difference operators used for the double model flow by Letcher are applied to a two dimensional cylinder with the Kelvin free surface condition and the numerical errors with these schemes are compared with those by the panel shift method. Then, 3-D waves due to a submerged source are calculated by the difference schemes, the panel shift method and also by a higher order boundary element method(HOBEM). Finally, the waves and wave resistance for Wigley's hull are calculated with these three schemes. It is shown that the panel shift method is free of numerical damping and dispersion error and performs better than the difference schemes. However, it can be concluded that the HOBEM also free of the numerical damping and dispersion error is the most stable, accurate and efficient.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.16 no.3
• /
• pp.193-204
• /
• 2012

This paper is comparing numerical schemes for a differential equation with convection and fourth-order diffusion. Our model equation is $h_t+(h^2-h^3)_x=-(h^3h_{xxx})_x$, which arises in the context of thin film flow driven the competing effects of an induced surface tension gradient and gravity. These films arise in thin coating flows and are of great technical and scientific interest. Here we focus on the several numerical methods to apply the model equation and the comparison and analysis of the numerical results. The convection terms are treated with well known WENO methods and the diffusion term is treated implicitly. The diffusion and convection schemes are combined using a fractional step-splitting method.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2006.05a
• /
• pp.22-27
• /
• 2006

This paper describes the development of the stream-tube based dispersion model for modeling contaminant transport in open channels. The operator-splitting technique is employed to separate the 2D contaminant transport equation into the pure advection and pure dispersion equations. Then the total variation diminishing (TVD) schemes are combined with the second-order Lax-Wendroff and third-order QUICKEST explicit finite difference schemes respectively to solve the pure advection equation in order to prevent the occurrence of numerical oscillations. Due to various limiters owning different features, the numerical tests for 1D pure advection and 2D dispersion are conducted to evaluate the performance of different TVD schemes firstly, then the TVD schemes are applied to experimental data for simulating the 2D mixing in a straight trapezoidal channel to test the model capability. Both the numerical tests and model application show that the TVD schemes are very competent for solving the advection-dominated transport problems.

• Advances in Energy Research
• /
• v.5 no.3
• /
• pp.255-275
• /
• 2017

The advantages of using natural circulation (NC) as a cooling system, has prompted the worldwide development to investigate this phenomenon more than before. The interesting application of the NC in low power experimental facilities and research reactors, highlights the obligation of study in these laminar flows. The inherent oscillations of NC between hot source and cold sink in low Grashof numbers necessitates stability analysis of cooling flow with experimental or numerical schemes. For this type of analysis, numerical methods could be implemented to desired mass, momentum and energy equations as an efficient instrument for predicting the behavior of the flow field. In this work, using the explicit, implicit and Crank-Nicolson methods, the fluid flow parameters in a natural circulation experimental test loop are obtained and the sensitivity of solving approaches are discussed. In this way, at first, the steady state and transient results from explicit are obtained and compared with experimental data. The implicit and crank-Nicolson scheme is investigated in next steps and in subsequent this research is focused on the numerical aspects of instability prediction for these schemes. In the following, the assessment of the flow behavior with coarse and fine mesh sizes and time-steps has been reported and the numerical schemes convergence are compared. For more detail research, the natural circulation of fluid was modeled by ANSYS-CFX software and results for the experimental loop are shown. Finally, the stability map for rectangular closed loop was obtained with employing the Nyquist criterion.

• The Bulletin of The Korean Astronomical Society
• /
• v.39 no.2
• /
• pp.119.2-119.2
• /
• 2014

Many astrophysical phenomena involve processes of magnetohydrodynamics (MHD) and relativistic magnetohydrodynamics (RMHD). A number of numerical schemes have been developed to solve the equations of ideal MHD and RMHD. Recent codes are based on upwind schemes which solve hyperbolic systems of equations following the characteristics of the systems. Upwind schemes stand out by their robustness, clarity of the underlying physical model, and ability of achieving high resolution. We present MHD and RMHD codes based on the total variation diminishing (TVD) and weighted essentially non-oscillatory (WENO) schemes, which are second and higher order accurate extensions of upwind schemes. We demonstrate the ability and limitation of codes based on upwind schemes through a series of tests.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.27 no.4
• /
• pp.272-280
• /
• 2023

We study the numerical analysis for the Cahn-Hilliard (CH) equation using the decoupled projection (DP) method. The CH equation is a fourth order nonlinear partial differential equation that is hard to solve. Therefore, various of numerical schemes have been proposed to solve the CH equation. To verify the relation of each existing scheme for the CH equation, we consider the DP method for linear convex splitting schemes. We present the numerical experiments to demonstrate our analysis. Throughout this study, it is expected to construct a novel numerical scheme using the relation with existing numerical schemes.

• 한국전산유체공학회:학술대회논문집
• /
• 2005.04a
• /
• pp.128-134
• /
• 2005

In the CFD area, the numerical analysis of high Mach number flow over a blunt-body poses many issues. Various numerical schemes have been developed to cover the issues, but the traditional schemes are still used widely due to the complexities of new schemes and intricacy of modifying the established codes. In the present study, the well-known Roe's FDS based on TVD-MUSCL scheme is used for the solution of very high Mach number three-dimensional flows posing carbuncle and non-physical phenomena in numerical analysis. A parametric study was carried out to account for the effects of the entropy fixing, grid configurations and initial condition. The carbuncle phenomena could be easily overcome by the entropy fixing, and the non-physical solution could be eliminated by the use of the modified initial condition regardless of entropy fixing and grid configurations.