• 한국해양공학회:학술대회논문집
• /
• 한국해양공학회 2004년도 학술대회지
• /
• 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년도 추계 학술대회논문집
• /
• 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년도 정기총회 및 학술발표대회
• /
• 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.

• 대한전자공학회논문지
• /
• 제26권1호
• /
• pp.48-61
• /
• 1989

본 논문에서는 게이트의 길이가 0.7${\mu}m$인 n형 GaAs MESFET를 2차원적으로 수치 해석하였으며, 이동도를 국부 전계의 함수로 취하는 드리프트 -확산 모델을 사용하였다. 이산화 방법으로는 종래에 사용되던 FDM(finite difference method), FEM(finite element method)을 사용치 아낳고 Control-Volume Formulation을 사용하였으며, numerical scheme으로는 기존의 hybrid scheme이나 upwind scheme 대신에 exponential scheme과 거의 근사한 power-law scheme을 사용하였다. 이때 드리프트 속도와 확산 속도의 비율을 나타내는 Peclet number의 개념을 사용하였으며, 이 개념을 사용하여 control volume의 경계에서 numerical scheme을 고려한 전류식을 제안하였다. 앞에서 고려한 모델들과 수치해석 방법을 사용하여 시뮬레이션한 I-V 특성은 기존 노문의 결과와 일치하였다. 따라서 본 논문의 결과가 GaAs MESFET를 위한 유용한 2차원 시뮬레이터가 될 수 있음을 확인하였다. 또한 I-V 특성외에 채널 밑바닥에서이 속도 및 전계 분포를 통해 드리프트-확산 모델을 고려한 경우에 발생하는 속도 포화의 메카니즘을 제시했고, Dipole의 발생위치 및 발생 원인과 드레인 전류와의 관계 등에 대해서도 제시했다.

• 한국항해항만학회지
• /
• 제27권5호
• /
• pp.473-478
• /
• 2003

평행한 벽 사이에 원형배관(circular pipeline)을 놓고 그 주위의 유동특성에 대한 수치연구를 수행하였다. 비압축성 유체를 가지고, Navier-Stokes 방정식을 풀었고 3차 풍상(upwind) 차분의 수치해법을 이용하였다. 한쪽 벽과의 거리가 매우 작아질 때, 볼텍스 떨어짐이 상당히 억압되는데 이것은 벽 경계와의 상호 박리 작용 때문으로 간수된다. 본 연구는 레이놀드 수의 변함과 물체가 벽에 접근함에 따른 볼텍스 떨어짐의 특성을 규명하는데 있다. 원형배관 후류와 평행벽내 유기된 박리의 상호작용을 집중적으로 다루며 서로 다른 조건에서 박리와 와역한(vorticity dynamics)의 특성을 해석하였다.

• 대한수학회보
• /
• 제55권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년도 정기총회 및 학술발표대회
• /
• 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
• /
• 제15권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.

• 천문학회지
• /
• 제28권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.

• 대한전기학회논문지:전기물성ㆍ응용부문C
• /
• 제53권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.