• 한국수자원학회:학술대회논문집
• /
• 한국수자원학회 2007년도 학술발표회 논문집
• /
• pp.1392-1396
• /
• 2007

The governing equations in generalized curvilinear coordinates for 3D laminar flow are the Incompressible Navier-Stokes (INS) equations with the artificial dissipative terms. and continuity equation discretized using a second-order accurate, finite volume method on the nonstaggered computational grid. This method adopts a dual or pseudo time-stepping Artificial Compressibility (AC) method integrated in pseudo-time. Multigrid methods are also applied because solving the equations on the coarse grids requires much less computational effort per iteration than on the fine grid. The algorithm yields practically identical velocity profiles and secondary flows that are in excellent overall agreement with an experimental measurement (Humphrey et al., 1977).

• 한국전산유체공학회지
• /
• 제19권4호
• /
• pp.45-51
• /
• 2014

We develop a numerical method for solving the radiation hydrodynamic equations in one-dimensional spherical coordinates. The present method is validated through simulations of shock tube, thermal radiative diffusion and point explosion problems. The transient growth of the fireball is investigated by varying explosion yields. The present study clearly captures well-known breakaway phenomena related to the shock separation between pressure waves and thermal shock front. The fireball radius at the breakaway point is roughly increased by the yield to power of 0.4.

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

This paper presents the results of the application of a computational fluid dynamics algorithm for the simulation of plasma flows of arc-heated jet. The underlying physical model is based on the axisymmetric form of the conservation equations that are coupled with an arc model including Ohm heating, electromagnetic forces. The arc model given as a source term in fluid dynamic equations is determined by a solution of electric potential field governed by an elliptic partial differential equation. The governing equation of electric field is loosely coupled with fluid dynamic equations by an electric conductivity that is a function of state variables. However, the electric fields and flow fields cannot be solved In fully coupled manner, but should be solved iteratively due to the different characteristics of governing equations. With this solution approach, several applications of arc flow analysis will be presented including Arc Thruster and Circuit Breaker.

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

Present study examines the numerical issues of cell structure simulation for various regimes of detonation phenomena ranging from weakly unstable to highly unstable detonations. Inviscid fluid dynamics equations with $variable-{\gamma}$ formulation and one-step Arrhenius reaction model are solved by a MUSCL-type TVD scheme and 4th order accurate Runge-Kutta time integration scheme. A series of numerical studies are carried out for the different regimes of the detonation phenomena to investigate the computational requirements for the simulation of the detonation wave cell structure by varying the reaction constants and grid resolutions. The computational results are investigated by comparing the solution of steady ZND structure to draw out the minimum grid resolutions and the size of the computational domain for the capturing cell structures of the different regimes of the detonation phenomena.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2006년도 PARALLEL CFD 2006
• /
• pp.185-187
• /
• 2006

• 한국전산구조공학회:학술대회논문집
• /
• 한국전산구조공학회 2006년도 정기 학술대회 논문집
• /
• pp.389-394
• /
• 2006

Fully flexible cell preserves Hamiltonian in structure, so the symplectic time integrator is applied to the equations of motion. Primarily, generalized leapfrog time integration (GLF) is applicable, but the equations of motion by GLF have some of implicit formulas. The implicit formulas give rise to a complicate calculation for coding and need an iteration process. In this paper, the time integration formulas are obtained for the fully flexible cell molecular dynamics simulation by using the splitting time integration. It separates flexible cell Hamiltonian into terms corresponding to each of Hamiltonian term, so the simple and completely explicit recursion formula was obtained. The explicit formulas are easy to implementation for coding and may be reduced the integration time because they are not need iteration process. We are going to compare the resulting splitting time integration with the implicit generalized leapfrog time integration.

• 대한수학회보
• /
• 제53권1호
• /
• pp.247-262
• /
• 2016

We propose and analyze spectral and pseudo-spectral Jacobi-Galerkin approaches for weakly singular Volterra integral equations (VIEs). We provide a rigorous error analysis for spectral and pseudo-spectral Jacobi-Galerkin methods, which show that the errors of the approximate solution decay exponentially in $L^{\infty}$ norm and weighted $L^2$-norm. The numerical examples are given to illustrate the theoretical results.

• 한국CDE학회논문집
• /
• 제7권4호
• /
• pp.219-232
• /
• 2002

In numerically controlled(NC) machining simulation, a Z-map has been used frequently for representing a workpiece. Since the Z-map is usually represented by a set of Z-axis aligned vectors, the machining process can be simulated through calculating the intersection points between the vectors and the surface swept by a machining tool. In this paper, we present an efficient method to calculate those intersection points when an APT-type tool moves along a linear tool path. Each of the intersection points can be expressed as the solution of a system of non-linear equations. We transform this system of equations into a single-variable equation, and calculate the candidate interval in which the unique solution exists. We prove the existence of a solution and its uniqueness in this candidate interval. Based on these characteristics, we can effectively apply numerical methods to finally calculate the solution of the non-linear equations within a given precision. The whole process of NC simulation can be achieved by updating the Z-map properly. Our method can provide more accurate results with a little more processing time, in comparison with the previous closed-form solution.

• 한국전산유체공학회지
• /
• 제11권4호
• /
• pp.90-100
• /
• 2006

A digital wave tank (DWT) simulation technique has been developed by authors to investigate the interactions of fully nonlinear waves with 3D marine structures. A finite-difference/volume method and a modified marker-and-cell (MAC) algorithm have been used, which are based on the Navier-Stokes (NS) and continuity equations. The fully nonlinear kinematic free-surface condition is implemented by the marker-density function (MDF) technique or the Level-Set (LS) technique developed for one or two fluid layers. In this paper, some applications for various engineering problems with free-surface are introduced and discussed. It includes numerical simulation of marine environments by simulation equipments, fully nonlinear wave motions around offshore structures, nonlinear ship waves, ship motions in waves and marine flow simulation with free-surface. From the presented simulations, it seems that the developed DWT simulation technique can handle various engineering problems with free-surface and reliably predict hydrodynamic features due to the fully-nonlinear wave motions interacting with such marine structures.

• 한국전산유체공학회지
• /
• 제15권2호
• /
• pp.35-40
• /
• 2010

In the present work, LES with new variational multiscale method is conducted on the fully developed channel flow with Reynolds number, 180 based on the friction velocity and the channel half width. Incompressible Navier-Stokes equations are integrated using finite element method with the basis function of NURBS. To solve space-time equations, Newton's method with two stage predictor multicorrector algorithm is employed. The code is parallelized using MPI. The computational domain is a rectangular box of size $2{\pi}{\times}2{\times}4/3{\pi}$ in the streamwise, wall normal and spanwise direction. Mean velocity profiles and velocity fluctuations are compared with the data of DNS. The results agree well with those of DNS and other traditional LES.