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

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

Direct simulations of aerodynamic sound, especially sound emitted by rapidly rotating elliptic cylinder by the finite difference lattice Boltzmann method (FDLBM). Effect of pile-fabrics for noise reduction is also studied by the finite volume LBM (FVLBM) using an unstructured grid. Second order time integration and third order upwind scheme are shown to be enough for these simulations. Sound sources are detected to be doublets for both cases. For the elliptic cylinder, the doublet is generated in the interaction between the vortex and the edge. For the circular cylinders, they are generated synchronizing with the Karman vortex street, and it is also shown that the pile-fabrics covering the surface of the cylinder reduces the strength of the source.

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

For numerical analysis of a tidal flow, a two-dimensional hydrodynamic model is developed by solving the nonlinear shallow-water equations. The governing equations are discretized explicitly with a finite difference leap-frog scheme and a first-order upwind scheme on adaptive hierarchical quadtree grids. The developed model is verified by applying to prediction of tidal behaviors. The calculated tidal levels are compared to available field measurements. A very reasonable agreement is observed.

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

• Journal of Navigation and Port Research
• /
• v.27 no.5
• /
• pp.473-478
• /
• 2003

Numerical study was made on the flow characteristics around a circular pipeline between parallel walls. The incompressible Navier-Stokes equations were solved by using a third-order upwind differential scheme. When the distance near a wall is small enough, the vortex shedding is almost completely suppressed because of the interaction with the wall boundary layer separation. This study aims to clarify the characteristics of the vortex shedding regime as the body approaches a wall as Reynolds number varies. The feature of separated vorticity dynamics is analyzed at different conditions with particular attention to the interaction between the pipeline wake and the induced separation on the plane walls.

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

• 한국방재학회:학술대회논문집
• /
• 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.

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

• Journal of The Korean Astronomical Society
• /
• v.28 no.2
• /
• pp.223-243
• /
• 1995

We describe the implementation of a multi-dimensional numerical code to solve the equations for idea! magnetohydrodynamics (MHD) in cylindrical geometry. It is based on an explicit finite difference scheme on an Eulerian grid, called the Total Variation Diminishing (TVD) scheme, which is a second-order-accurate extension of the Roe-type upwind scheme. Multiple spatial dimensions are treated through a Strang-type operator splitting. Curvature and source terms are included in a way to insure the formal accuracy of the code to be second order. The constraint of a divergence-free magnetic field is enforced exactly by adding a correction, which involves solving a Poisson equation. The Fourier Analysis and Cyclic Reduction (FACR) method is employed to solve it. Results from a set of tests show that the code handles flows in cylindrical geometry successfully and resolves strong shocks within two to four computational cells. The advantages and limitations of the code are discussed.

• The Transactions of the Korean Institute of Electrical Engineers C
• /
• v.53 no.8
• /
• pp.433-442
• /
• 2004

Various models and various boundary conditions have been suggested for fluid transport simulations of high density plasma discharges such as the inductively coupled plasma discharge. In this work, we compare the various models using one-dimensional simulations based on the FDM(finite difference method), the upwind scheme, the power-law scheme, and the dielectric relaxation scheme［l］ Comparing the exactness, the numerical stability and the efficiency of the various models. the most adoptable model is suggested.