• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 1996년도 춘계 학술대회논문집
• /
• pp.131-136
• /
• 1996

The numerical method for the flow field of a gas atomization process is presented. For the analysis of the compressible supersonic jet flow of a gas. an axisymmetric Navier-Stokes equations are solved using a LU-factored upwind method. The MUSCL type TVD scheme is used for the discretization of inviscid flux, whereas Steger-Warming splitting and LU factorization is applied to the implicit operator. For the validation of the present method, we computed the flow field around the simple gas atomizer proposed by Issac. The numerical results has shown excellent agreement with the experimental data.

• 한국전산구조공학회논문집
• /
• 제26권1호
• /
• pp.39-48
• /
• 2013

본 논문은 기존의 1차원 Stefan 문제를 해석할 수 있는 이동최소제곱 차분법을 확장하여 복잡한 계면경계 형상을 갖는 2차원 문제에 적용할 수 있는 수치기법을 개발한다. 1차원 경우와 달리 2차원 영역에서 임의로 움직이는 이동경계의 위상변화를 효과적으로 모델링할 수 있는 기법을 제안했으며, 이동경계 모사시 절점만 사용하는 이동최소제곱 차분법의 강점을 그대로 살리면서 이동경계의 불연속 특이성과 kinetics 조건을 정확하게 만족시키는 이동최소제곱 미분근사식을 제시했다. 평형방정식은 implicit(음해)법으로 차분하여 수치 안정성을 확보했으며, 이동경계는 explicit(양해)법으로 update하여 계산효율성의 극대화했다. 몇 가지 수치예제를 통해 개발된 이동최소제곱 차분법이 다양한 계면경계 형상을 갖는 2차원 Stefan 문제를 정확하고 효율적으로 풀 수 있음을 검증했다.

• 한국소성가공학회:학술대회논문집
• /
• 한국소성가공학회 1996년도 자동차부품 제작기술의 진보
• /
• pp.97-109
• /
• 1996

In the present work computations are carried out for analysis of complicated sheet metal forming process such as forming of a rear hinge. Finite element formulation using dynamic explicit time integration scheme and step-wise combined Implicit/Explicit scheme are introduced for numerical analysis of sheet metal forming process. The rigid-plastic finite element method based on membrane elements has long been employed as a useful numerical technique for the analysis of sheet metal forming because of its time effectiveness. The explicit scheme in general use is based on the elastic-plastic modelling of material requiring large computation time. In finite element simulation of sheet metal forming processes, the robustness and stability of computation are important requirements since the computation time and convergency become major points of consideration besides the solution accuracy due to the complexity of geometry and boundary conditions. The implicit scheme employs a more reliable and rigorous scheme in considering the equilibrium at each step of deformation, while in the explicit scheme the problem of convergency is eliminated at the cost of solution accuracy. The explicit approach and the implicit approach have merits and demerits, respectively. In order to combine the merits of these two methods a step-wise combined implicit/explicit scheme has been developed.

• 대한기계학회논문집B
• /
• 제29권5호
• /
• pp.609-616
• /
• 2005

While the SIMPLE algorithm is most widely used for the simulations of flow phenomena that take place in the industrial equipments or the manufacturing processes, it is less adopted for the simulations of the free surface flow. Though the SIMPLE algorithm is free from the limitation of time step, the free surface behavior imposes the restriction on the time step. As a result, the explicit schemes are faster than the implicit scheme in terms of computation time when the same time step is applied to, since the implicit scheme includes the numerical method to solve the simultaneous equations in its procedure. If the computation time of SIMPLE algorithm can be reduced when it is applied to the unsteady free surface flow problems, the calculation can be carried out in the more stable way and, in the design process, the process variables can be controlled based on the more accurate data base. In this study, a modified SIMPLE algorithm is presented fur the free surface flow. The broken water column problem is adopted for the validation of the modified algorithm (MoSIMPLE) and for comparison to the conventional SIMPLE algorithm.

• Journal of applied mathematics & informatics
• /
• 제26권1_2호
• /
• pp.1-14
• /
• 2008

In this paper, A Riesz fractional diffusion equation with a nonlinear source term (RFDE-NST) is considered. This equation is commonly used to model the growth and spreading of biological species. According to the equivalent of the Riemann-Liouville(R-L) and $Gr\ddot{u}nwald$-Letnikov(G-L) fractional derivative definitions, an implicit difference approximation (IFDA) for the RFDE-NST is derived. We prove the IFDA is unconditionally stable and convergent. In order to evaluate the efficiency of the IFDA, a comparison with a fractional method of lines (FMOL) is used. Finally, two numerical examples are presented to show that the numerical results are in good agreement with our theoretical analysis.

• 한국윤활학회:학술대회논문집
• /
• 한국윤활학회 1999년도 제30회 추계학술대회
• /
• pp.146-153
• /
• 1999

A numerical dynamic simulation is necessary to investigate the capacity of the HDD. The slider surface become more and more complicated to make the magnetized area smaller and readback signal stronger. So a numerical dynamic simulation must be preceded to develop a new slider in HDD. The dynamic simulations of air-lubricated slider bearing have been peformed using FIFD(Factored Implicit Finite Difference) method. The governing equation, Reynolds equation Is modified with Fukui and Kaneko model(FK model) which includes the first and the second-order slip. The equations of motion for the slider bearing are solved simultaneously with the modified Reynolds equation for the case of three degrees of freedom. The slider transient response for disk step bump and slider impulse force is given for various case and for iteration algorithm and new algorithm.

• 한국전산유체공학회지
• /
• 제2권1호
• /
• pp.8-12
• /
• 1997

The three-dimensional Euler equations are solved numerically for the analysis of contraction flows in wind or water tunnels. A second-order finite difference method is used for the spatial discretization on the nonstaggered grid system and the 4-stage Runge-Kutta scheme for the numerical integration in time. In order to speed up the convergence, the local time stepping and the implicit residual-averaging schemes are introduced. The pressure field is obtained by solving the pressure-Poisson equation with the Neumann boundary condition. For the evaluation of the present Euler solver, numerical computations are carried out for three contraction geometries, one of which was adopted in the Large Cavitation Channel for the U.S. Navy. The comparison of the computational results with the available experimental data shows good agreement.

• 한국해양공학회지
• /
• 제11권4호
• /
• pp.141-151
• /
• 1997

This paper deals with a numerical study of flow-guider applied to controlling current in a bay. Two dimensional numerical model for tidal currents based on the depth averaged equation is developed and standard k-.epsilon. model is adopted to determine the turbulence diffusion. Equations are described in a generalized coordinate system to be implemented by non-staggered grid system and discretized by using finite volume method. Unsteady flow is simulated by fully implicit scheme. Hybrid scheme and central differencing are used to compute the convective terms and source terms, respectively. The tidal current in a rectangular bay is simulated and it gives satisfactory results. The realistic and distinct models of a large structure placed in bay are also exemplified with or without flow-guiders. The simulation results show that the flow-guider gives the residual tidal current in the bay by the different flux with respect to the direction of tidal current.

• Wind and Structures
• /
• 제34권1호
• /
• pp.59-72
• /
• 2022

In this work the numerical results of the flow around a 5:1 rectangular cylinder at Reynolds numbers 3 000 and 40 000, zero angle of attack and smooth incoming flow condition are presented. Implicit Large Eddy Simulations (ILES) have been performed with a high-order accurate spatial scheme and an implicit high-order accurate time integration method. The spatial approximation is based on a discontinuous Galerkin (dG) method, while the time integration exploits a linearly-implicit Rosenbrock-type Runge-Kutta scheme. The aim of this work is to show the feasibility of high-fidelity flow simulations with a moderate number of DOFs and large time step sizes. Moreover, the effect of different parameters, i.e., dimension of the computational domain, mesh type, grid resolution, boundary conditions, time step size and polynomial approximation, on the results accuracy is investigated. Our best dG result at Re=3 000 perfectly agrees with a reference DNS obtained using Nek5000 and about 40 times more degrees of freedom. The Re=40 000 computations, which are strongly under-resolved, show a reasonable correspondence with the experimental data of Mannini et al. (2017) and the LES of Zhang and Xu (2020).

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2000년도 추계 학술대회논문집
• /
• pp.129-134
• /
• 2000

An efficient numerical method to solve the unsteady incompressible Navier-Stokes equations is developed. A fully implicit time advancement is employed to avoid the CFL(Courant-Friedrichs-Lewy) restriction, where the Crank-Nicholson discretization is used for both the diffusion and convection terms. Based on a block LU decomposition, velocity-pressure decoupling is achieved in conjunction with the approximate factorization. Main emphasis is placed on the additional decoupling of the intermediate velocity components with only n th time step velocity The temporal second-order accuracy is Preserved with the approximate factorization without any modification of boundary conditions. Since the decoupled momentum equations are solved without iteration, the computational time is reduced significantly. The present decoupling method is validated by solving the turbulent minimal channel flow unit.