• 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 applied mathematics & informatics
• /
• v.41 no.1
• /
• pp.181-192
• /
• 2023

This work deals with the numerical modeling of dam-break flow over erodible bed. The mathematical model consists of the shallow water equations, the transport diffusion and the bed morphology change equations. The system is solved by central upwind scheme. The obtained results of the resolution of dam-beak problem is presented in order to show the performance of the numerical scheme. Also a comparison of central upwind and Roe schemes is presented.

• Journal of computational fluids engineering
• /
• v.8 no.3
• /
• pp.21-31
• /
• 2003

A non-dissipative and very accurate one-dimensional upwind leapfrog method is extended to higher-order and multi-dimensional advection and acoustic equations. The higher-order versions are developed by extending the stencils in space and time. The schemes are then successfully applied to the classical test cases for advection and acoustics.

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

This paper describes the numerical solution to the hyperbolic system of magnetohydrodynamic (MHD) equations. First, by pointing out the approximations involved, the deal MHD equations are presented. Next, the MHD waves as well as the associated shocks and discontinuities, are presented. Then, based on the hyperbolicity of the ideal MHD equations, the application of upwind schemes, which have been developed for hydrodynamics, is discussed to solve the equations numerically. As an definite example, one and multi-dimensional codes based on the Total Variation Diminishing scheme are presented. The treatment in the multi-dimensional code, which maintains ${\nabla}{\cdot}$B = 0, is described. Through tests, the robustness of the upwind schemes for MHDs is demonstrated.

• Journal of the korean Society of Automotive Engineers
• /
• v.11 no.3
• /
• pp.29-36
• /
• 1989

A multidimensional numerical simulation for flows in an engine with axisymmetric geometry was performed. Three kinds of differencing schemes, namely, skew upwind differencing scheme (SUDS), interpolated upwind differencing scheme (IUDS), upwind differencing scheme (UDS), are used in a comparative study. Simultaneously, the effects of the artificial dampings and the grids on numerical results are estimated. Compared with the measurements, the calculations with SUDS and proper artificial damping show very similar qualitative tendency with observed results. But there are some discrepancies due to numerical errors and unclear boundary conditions.

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

• Journal of Mechanical Science and Technology
• /
• v.15 no.5
• /
• pp.671-680
• /
• 2001

Boundary condition is one of the major factors to influence the numerical stability and solution accuracy in numerical analysis. One of the most important physical boundary conditions in the flowfield analysis is the wall boundary condition imposed on the body surface. To solve a two-dimensional Euler equation, totally four numerical wall boundary conditions should be prescribed. Two of them are supplied by the flow tangency condition. The other two conditions, therefore, should be prepared additionally in a suitable way. In this paper, four different sets of wall boundary conditions are proposed and then applied to solve high-speed flowfields around a quarter circle geometry. A two-dimensional compressible Euler solver is prepared based on the finite volume method. This solver hires three different upwind schemes; Steger-Warmings flux vector splitting, Roes flux difference splitting, and Lious advection upstream splitting method. It is found that the way to specify the additional numerical wall boundary conditions strongly affects the overall stability and accuracy of the upwind schemes in high-speed flow calculation. The optimal wall boundary conditions should be also chosen very carefully depending on the numerical schemes used to solve the problem.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.27 no.8
• /
• pp.1122-1133
• /
• 2003

Numerical characteristics of implicit upwind schemes, such as upwind ADI, line Gauss-Seidel (LGS) and point Gauss-Seidel (LU) algorithms, for Navier-Stokes equations have been investigated. Time-derivative preconditioning method was applied for efficient convergence at low Mach/Reynolds number regime as well as at large grid aspect ratios. All the algorithms were expressed in approximate factorization form and von Neumann stability analysis was performed to identify stability characteristics of the above algorithms in the presence of high grid aspect ratios. Stability analysis showed that for high aspect ratio computations, the ADI and LGS algorithms showed efficient damping effect up to moderate aspect ratio if we adopt viscous preconditioning based on min-CFL/max-VNN time-step definition. The LU algorithm, on the other hand, showed serious deterioration in stability characteristics as the grid aspect ratio increases. Computations for several practical applications also verified these results.

• Journal of Ocean Engineering and Technology
• /
• v.20 no.4 s.71
• /
• pp.31-36
• /
• 2006

Numerical schemes for spacing and time are tested to accurately simulate the wave propagation. The tested numerical schemesinclude 2nd-order central differencing, l-order upwind scheme, 2nd-order Leith scheme, 3rd-order MUSCLE, QUICK and QUICKEST schemes in spacing and the Euler and 4th-order Runge-Kutta(R-K) schemes in time. It is seen that more accurate results are expected when the higher-order schemes, especially the schemes combined with a TVD control limiter, are used for solving the wave equation. The 3rd-order upwind scheme with limiter and the 4th-order R-K scheme in tim￡ are finally applied to the wave-making simulation in a digital wave tank.

• Proceedings of the KSME Conference
• /
• 2004.11a
• /
• pp.1294-1299
• /
• 2004

Transient radiative heat transfer is analyzed in a one-dimensional slab using finite volume method (FVM). In this study, the step, $2^{nd}$ order upwind, and QUICK schemes are used for incident diffuse radiation and collimated beam, respectively. The results for diffuse radiation show that all schemes applied in this study give good agreements with available published results. In case of collimated beam, however, the results show deviations from the analytical solutions. To successfully describe the propagations of collimated beam, shock capturing schemes such as TVD scheme are need to be developed.