• Bulletin of the Korean Mathematical Society
• /
• v.55 no.3
• /
• pp.881-898
• /
• 2018

This article deals with implementation of a high-order finite difference scheme for numerical solution of the incompressible Navier-Stokes equations on curvilinear grids. The numerical scheme is based on pseudo-compressibility approach. A fifth-order upwind compact scheme is used to approximate the inviscid fluxes while the discretization of metric and viscous terms is accomplished using sixth-order central compact scheme. An implicit Euler method is used for discretization of the pseudo-time derivative to obtain the steady-state solution. The resulting block tridiagonal matrix system is solved by approximate factorization based alternating direction implicit scheme (AF-ADI) which consists of an alternate sweep in each direction for every pseudo-time step. The convergence and efficiency of the method are evaluated by solving some 2D benchmark problems. Finally, computed results are compared with numerical results in the literature and a good agreement is observed.

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

• 한국방재학회:학술대회논문집
• /
• 2008.02a
• /
• pp.431-434
• /
• 2008

Pratical dispersion-correction scheme is applicated to simulate the distant propagation of tsunami. This scheme is based on the leap-frog finite difference scheme for the linear shallow-water equations. The new scheme has the advantage of using the constant spatial grid size and time step size even in area of variable depths. And this new model constructed by using the 2nd upwind scheme, dynamic linking method, and staggered grid system. This model is simulated to near Sokcho harbor about The Central East Sea Tsunami in 1983. And this result is compared to tide gage and result of former model.

• 한국전산유체공학회:학술대회논문집
• /
• 2000.10a
• /
• pp.54-59
• /
• 2000

A numerical study of axisymmetric rocket main nozzle flow has been accomplished. The CSCM upwind flux difference splitting method with an iterative time marching scheme having second order accuracy in time and space has been used to simulate unsteady flow characteristics in an axisymmetric rocket main nozzle. Though the pressure vary at nozzle inlet with the lapse of time, Mach No. and the density were not changed significontly compared with the temperature. Specific heat ratio $\gamma$=1.134 predicted higher temperature at nozzle throat and exit and nondimensional thrust coefficients at exit than specific heat ratio $\gamma$=1.4 did.

• Journal of Wetlands Research
• /
• v.5 no.1
• /
• pp.67-78
• /
• 2003

The unsteady 2D convection and diffusion equation is solved numerically for the real-time simulation of suspended load propagation. The streamlined upwind scheme efficiently reduces numerical oscillations due to the high Peclet number in the convection dominant flow. By using the mixed boundary condition to express the external source terms or externally induced suspended load as a function of time in the algorithm, the model is capable of handling not only continuous load cases but also non-continuous suspended load influx. The suspended load transport modelwas verified using a case study for which an analytical exact solution is available and was applied to the real-time simulation of a suspended load influx case on the Mississippi River. The model algorithm can provide a framework upon which water quality as well as contaminant transport models can be built.

• Proceedings of the Korea Committee for Ocean Resources and Engineering Conference
• /
• 2004.11a
• /
• pp.154-159
• /
• 2004

A new finite difference scheme is studied fur the simulation of free surface, surface the third derivative tenn for the wave elevation is artificially added in the the free-surface boundary condition. This study presents a comparative analysis with simulations performed by using the classical MAC method. More systematic computations are carried out by changing the submergence depth and angle of attack. Through the numerical simulation, it is found that a new numerical method becomes more efficient for the reason that the free surface elevation is reasonably developed at tire rear of trailing edge.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.9
• /
• pp.1148-1156
• /
• 2000

The QUICK scheme for convection terms is modified for unstructured finite volume method by using linear reconstruction technique and validated through the computation of two well defined laminar flows. It uses two upstream grid points and one downstream grid point in approximating the convection terms. The most upstream grid point is generated by considering both the direction of flow and local grid line. Its value is calculated from surrounding grid points by using a linear construction method. Numerical error by the modified QUICK scheme is shown to decrease about 2.5 times faster than first order upwind scheme as grid size decreases. Computations are also carried out to see effects of the skewness and irregularity of grid on numerical solution. All numerical solutions show that the modified QUICK scheme is insensitive to both the skewness and irregularity of grid in terms of the accuracy of solution.

• 한국방재학회:학술대회논문집
• /
• 2008.02a
• /
• pp.233-236
• /
• 2008

To analyze the saline intrusion in the place, such as an estuary, the three-dimensional numerical model is developed. In this study, the advection terms of the governing equations are discretized by upwind scheme. By using an explicit scheme for the longitudinal direction and an implicit scheme for the vertical direction, the numerical model is free from the restriction of temporal step size caused by a relatively small grid ratio. The equation of state is used to consider the density, and the scalar transport equation for salinity is employed the third order TVD to scheme to prevent unphysical oscillation near discontinuity. In order to verify saline intrusion, the numerical model is conducted to compare the previous model in the lock exchange. The present model generally show a good agreement with the previous one.

• Journal of the Korean Institute of Telematics and Electronics
• /
• v.26 no.1
• /
• pp.48-61
• /
• 1989

In this paper, two-dimensional numerical simulation of GaAs MESFFT with 0.7${\mu}m$ gate length is perfomed. Drift-diffusion model which consider that mobility is a function of local electric field, is used. As a discretization method, instead of FDM (finite difference method) and FEM (finite element method), the Control-Volume Formulation (CVF) is used and as a numerical scheme current hybrid scheme or upwind scheme is replaced by power-law scheme which is very approximate to exponential scheme. In the process of numerical analysis, Peclet number which represents the velocity ratio of drift and diffusion, is introduced. And using this concept a current equation which consider numerical scheme at the interface of control volume, is proposed. The I-V characteristics using the model and numerical method has a good agreement with that of previous paper by others. Therefore, it is confined that it may be useful as a simulator for GaAs MESFET. Besides I-V characteristics, the mechanism of both velocity saturation in drift-diffusion model is described from the view of velocity and electric field distribution at the bottom of the channel. In addition, the relationship between the mechanism and position of dipole and drain current, are described.

• 유체기계공업학회:학술대회논문집
• /
• 2001.11a
• /
• pp.23-29
• /
• 2001

Three-dimensional flow through a tubular centrifugal fan with airfoil type blades are analyzed using CFX-TASCflow. Standard k - $\epsilon$ model and k - $\omega$ model are used as turbulence closures. The numerical schemes for convetion terms, i.e., Upwind Differencing Scheme(UDS), Mass Weighted Skewed Upstream Differencing Scheme(MWS), Linear Profile Skewed Upstream Differencing Scheme(LPS), and Modified Linear Profile Skewed Upstream Differencing Scheme(MLPS) are also tested. And, the performance of these schemes coupled with two turbulence models are evaluated. Computational static pressure distributions are compared with experimental data obtained in this work.