• 대한화학회지
• /
• 제42권4호
• /
• pp.398-410
• /
• 1998

Boltzmann 방정식의 비선형 해법으로서 cumulant 모멘트 방법을 연구하였으며, Maxwell 분자모형 단원자분자 기체계의 정상충격파 문제에 대하여 적용하였다. 모멘트 방정식의 해는 Maxwell-Ikenberry-Truesdell(MIT) 반복법을 사용하였다. 원래의 MIT 반복법은 초기값을 평형분포함수로부터 구하지만, 본 연구에서는 반복계산의 초기값을 Mott-Smith의 두방식(bimodal)함수로부터 구하였다. 모멘트 계산은 2차 반복단계까지 수행하였으며, 강한 충격파에 대한 밀도, 온도, stress, heat flux 등의 윤곽과 충격파의 두께, 그리고 마하수 1.4 미만의 약한 충격파의 두께를 계산하였다. 1차 반복계산에서 충격파 윤곽에 대한 간단한 형태의 해석적 표현을 얻었으며, 이로부터 도출한 약한 충격파 두께에 대한 극한법칙은 Navier-Stokes 이론과 정확히 일치한다. 2차 반복계산에 의한 결과는 강한 충격파의 윤곽곡선 및 충격파 두께가 Monte Carlo 문헌값과 정량적으로 일치함을 보인다.

• 한국해안·해양공학회논문집
• /
• 제23권3호
• /
• pp.248-257
• /
• 2011

본 연구에서는 규칙파동장에 있어서 인장계류된 원형부유체의 파랑응답해석에 2차원 Navier-Stokes solver에 기초한 새로운 수치파동수조모델을 제안하였다. 본 수치파동수조모델에서는 이동구조물과 유체와의 상호작용을 해석하기 위하여 직각좌표계에서 임의형상의 불투과경계를 갖는 구조물과 유체와의 연성해석이 가능한 IBM(Immersed Boundary Method)과 자유수면 추적을 위한 VOF(Volume of Fluid)법을 결합하였다. 부유체운동에 대한 수치결과를 기존의 FAVOR(Fractional Area Volume Obstacle Representation)법에 의한 계산결과 및 수리실험과 비교하였다. 게다가, 수치모델의 보다 자세한 검증을 위하여 원형부유체의 동요 및 자유수면변동에 관한 수리모형을 추가로 실시하였고, 제안된 수치모델의 범용성과 타당성을 검증하기 위하여 직사각형부유체에 대한 수치 및 수리실험도 병행하였다. 이로부터 추정된 수치계산결과는 실험결과를 잘 재현하고 있는 것으로 판단되었다.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2011년 춘계학술대회논문집
• /
• pp.65-67
• /
• 2011

The moment-of-fluid (MOF) method is a new volume-tracking method that accurately treats evolving material interfaces. The MOF method uses moment data, namely the material volume fraction, as well as the centroid, for a more accurate representation of the material configuration, interfaces and concomitant volume advection. In this paper, unstructured mesh extension of the MOF method is to be presented. The MOF method is coupled with a stabilized finite element incompressible Navier-Stokes solver for two materials. The effectiveness of the MOF method is demonstrated with a free-surface dam-break problem.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2009년 추계학술대회논문집
• /
• pp.223-227
• /
• 2009

In the present study, a least square weighted residual method and Taylor-Galerkin method were formulated and tested for the discretization of the two hyperbolic type equations of level set method; advection and reinitialization equations. The two approaches were compared by solving a time reversed vortex flow and three-dimensional broken dam flow by employing a four-step splitting finite element method for the solution of the incompressible Navier-Stokes equations. From the numerical experiments, it was shown that the least square method is more accurate and conservative than Taylor-Galerkin method and both methods are approximately first order accurate when both advection and reinitialization phase are involved in the evolution of free surface.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제20권1호
• /
• pp.83-106
• /
• 2016

This paper reviews and compares three different methods for modeling incompressible and immiscible ternary fluid flows: the immersed boundary, level set, and phase-field methods. The immersed boundary method represents the moving interface by tracking the Lagrangian particles. In the level set method, an interface is defined implicitly by using the signed distance function, and its evolution is governed by a transport equation. In the phase-field method, the advective Cahn-Hilliard equation is used as the evolution equation, and its order parameter also implicitly defines an interface. Each method has its merits and demerits. We perform the several simulations under different conditions to examine the merits and demerits of each method. Based on the results, we determine the most suitable method depending on the specific modeling needs of different situations.

• Journal of Korean Society for Atmospheric Environment
• /
• 제9권E호
• /
• pp.347-357
• /
• 1993

The Lagrangian grid-free numerical method, the discrete vortex method, was applied to solve the Navier-Stokes euqations. This method avoids the introduction of numerical viscosity swamping the real physical viscosity at high Reynolds number, unlike Eulerian method, e.g. finite difference and element methods. The boundary integral equation method for the potential flow solution was included to make the discrete vortex method more feasible for complex geometries. The fast adaptive multipole expansion method was incorporated to reduce the computational time from $O(N^2)$ to O(N) for the computations of vortex-vortex interactions. The test problems were air flow around one circular cylinder and two circular cylinders in tandem with various gaps. The numerical results were in excellent gareement with the experimental and other computational results. The applicabilty of the method was discussed with the indoor and the outdoor air pollution problems, especially the contaminant transport in the recirculation regions.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2008년도 추계학술대회B
• /
• pp.2558-2561
• /
• 2008

In the present study, a parallel Taylor-Galerkin/level set based two-phase flow code was developed using finite element discretization and domain decomposition method based on MPI (Message Passing Interface). The proposed method can be utilized for the analysis of a large scale free surface problem in a complex geometry due to the feature of FEM and domain decomposition method. Four-step fractional step method was used for the solution of the incompressible Navier-Stokes equations and Taylor-Galerkin method was adopted for the discretization of hyperbolic type redistancing and advection equations. A Parallel ILU(0) type preconditioner was chosen to accelerate the convergence of a conjugate gradient type iterative solvers. From the present parallel numerical experiments, it has been shown that the proposed method is applicable to the simulation of large scale free surface flows.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2000년도 춘계학술대회논문집B
• /
• pp.704-708
• /
• 2000

A new efficient numerical method for computing unsteady, incompressible flows, DRANS (Decoupled Reynolds-Averaged Navier-Stokes), is presented. To eliminate the restriction of CFL condition, a fully-implicit time advancement in which the Crank-Nicolson method is used fer both the diffusion and convection terms. is adopted. Based on decomposition method, the velocity-turbulent quantity decoupling is achieved. The additional decoupling of the intermediate velocity components in the convection term is made for the fully-implicit time advancement scheme. Since the iterative procedures for the momentum, ${\kappa}\;and\;{\varepsilon}$ equations are not required, the components decouplings bring fourth the reduction of computational cost. The second-order accuracy in time of the present numerical algorithm is ascertained by computing decaying vortices. The present decoupling method is applied to turbulent boundary layer with local forcing.

• 한국전산구조공학회논문집
• /
• 제23권6호
• /
• pp.675-682
• /
• 2010

The flood impact pressure acting on a vertical wall resulting from a dam-breaking problem is simulated using a navier-Stokes(N-S) solver. The N-S solver uses Eulerian Finite Volume Method(FVM) along with Volume Of Fluid(VOF) method for 2-D incompressible free surface flows. A Split Lagrangian Advection(SLA) scheme for VOF method is implemented in this paper. The SLA scheme is developed based on an algorithm of Piecewise Linear Interface Calculation(PLIC). The coupling between the continuity and momentum equations is affected by using a well-known Semi-Implicit Method for Pressure-Linked Equations (SIMPLE) algorithm. Several two-dimensional numerical simulations of the dam-breaking problem are presented to validate the accuracy and demonstrate the capability of the present algorithm. The significance of the time step and grid resolution are also discussed. The computational results are compared with experimental data and with computations by other numerical methods. The results showed a favorable agreement of water impact pressure as well as the global fluid motion.

• 한국추진공학회:학술대회논문집
• /
• 한국추진공학회 2008년 영문 학술대회
• /
• pp.260-265
• /
• 2008

A numerical scheme for solid propellant rocket has been studied using preconditioning method to research unsteady combustion processes for the double-base propellant with a converging-diverging nozzle. The Navier-Stokes equation is solved by dualtime stepping method with finite volume method. The turbulence model uses a shear stress transport modeling. The species equation follows up the method of Xinping WI, Mridul Kumar and Kenneth K. Kuo. A preconditioned algorithm is applied to solve incompressible regime inside the combustor and compressible flow at nozzle. Mass flux was evaluated using modified advective upwind splitting method. The simulated result the comparison a fully coupled implicit method and a semi implicit method in terms of accuracy and efficiency. This report shows the result of solid rocket propellant combustion.