• 대한기계학회논문집A
• /
• 제30권10호
• /
• pp.1320-1327
• /
• 2006

Parallel procedure and performance of two-step Eulerian code have not been reported sufficiently yet even though it was developed and utilized widely in the impact simulation. In this study, parallel strategy of two-step Eulerian code was proposed and described in detail. The performance was evaluated in the self-made linux cluster computer. Compared with commercial code, a relatively good performance is achieved. Through the performance evaluation of each computation stage, remap is turned out to be the most time consuming part among the other part such as FE processing, communication, time marching etc.

• 한국항공우주학회지
• /
• 제34권8호
• /
• pp.47-54
• /
• 2006

본 연구에서는 Lagrangian 단계와 Remap 단계로 이루어진 Two-step Eulerian 기법에 적용되는 수치기법에 대해 살펴보고, 이를 구현하였다. 요소 값은 donor cell 과 van Leer 기법으로 이류처리 하였고, 절점 값은 Simplifed ALE [SALE] 기법으로 이류처리 하였으며, 재료 경계면은 Simple Line Interface Calculation [SLIC] 기법으로 결정하였다. 또한, 프로그램 구현 과정과 Lagrangian 단계와 Remap 단계의 연결 시 주의 깊게 고려해야할 사항들에 대해 자세히 기술하였다. 본 연구를 통해 구현된 Two-step Eulerian 코드의 검증을 위해, 순수 Lagrangian 기법, 그리고 상용코드의 결과와 비교하여 해석 결과의 정확성을 검증하였다.

• Journal of Mechanical Science and Technology
• /
• 제16권5호
• /
• pp.683-695
• /
• 2002

The objective of the present study is to analyze the fluid flow with moving boundary using a finite element method. The algorithm uses a fractional step approach that can be used to solve low-speed flow with large density changes due to intense temperature gradients. The explicit Lax-Wendroff scheme is applied to nonlinear convective terms in the momentum equations to prevent checkerboard pressure oscillations. The ALE (Arbitrary Lagrangian Eulerian) method is adopted for moving grids. The numerical algorithm in the present study is validated for two-dimensional unsteady flow in a driven cavity and a natural convection problem. To extend the present numerical method to engine simulations, a piston-driven intake flow with moving boundary is also simulated. The density, temperature and axial velocity profiles are calculated for the three-dimensional unsteady piston-driven intake flow with density changes due to high inlet fluid temperatures using the present algorithm. The calculated results are in good agreement with other numerical and experimental ones.

• Journal of Mechanical Science and Technology
• /
• 제16권5호
• /
• pp.720-731
• /
• 2002

Flow with moving free surfaces is analyzed with an the Eulerian coordinate system. This study proposes a semi-implicit filling algorithm using VOF in which the PLIC (Piecewise Linear Interface Calculation) -type interface reconstruction method and the donor-acceptor-type front advancing scheme are adopted. Also, a new scheme using extrapolation of the stream function is proposed to find the velocity of the node that newly enters the computational domain. The effect of wall boundary conditions on the flow field and temperature field is examined by numerically solving a two-dimensional casting process.

• 한국전산구조공학회논문집
• /
• 제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.

• Structural Engineering and Mechanics
• /
• 제47권5호
• /
• pp.599-622
• /
• 2013

Seismic response of two dimensional liquid tanks is numerically simulated using fully nonlinear velocity potential theory. Galerkin-weighted-residual based finite element method is used for solving the governing Laplace equation with fully nonlinear free surface boundary conditions and also for velocity recovery. Based on mixed Eulerian-Lagrangian (MEL) method, fourth order explicit Runge-Kutta scheme is used for time integration of free surface boundary conditions. A cubic-spline fitted regridding technique is used at every time step to eliminate possible numerical instabilities on account of Lagrangian node induced mesh distortion. An artificial surface damping term is used which mimics the viscosity induced damping and brings in numerical stability. Four earthquake motions have been suitably selected to study the effect of frequency content on the dynamic response of tank-liquid system. The nonlinear seismic response vis-a-vis linear response of rectangular liquid tank has been studied. The impulsive and convective components of hydrodynamic forces, e.g., base shear, overturning base moment and pressure distribution on tank-wall are quantified. It is observed that the convective response of tank-liquid system is very much sensitive to the frequency content of the ground motion. Such sensitivity is more pronounced in shallow tanks.

• Wind and Structures
• /
• 제17권3호
• /
• pp.263-274
• /
• 2013

A numerical technique for fluid-structure interaction, which is based on the finite element method (FEM) and computational fluid dynamics (CFD), was developed for application to an industrial chimney equipped with a pendulum tuned mass damper (TMD). In order to solve the structural problem, a one-dimensional beam model (Navier-Bernoulli) was considered and, for the dynamical problem, the standard second-order Newmark method was used. Navier-Stokes equations for incompressible flow are solved in several horizontal planes to determine the pressure in the boundary of the corresponding cross-section of the chimney. Forces per unit length were obtained by integrating the pressure and are introduced in the structure using standard FEM interpolation techniques. For the fluid problem, a fractional step scheme based on a second order pressure splitting has been used. In each fluid plane, the displacements have been taken into account considering an Arbitrary Lagrangian Eulerian approach. The stabilization of convection and diffusion terms is achieved by means of quasi-static orthogonal subscales. For each period of time, the fluid problem was solved and the geometry of the mesh of each fluid plane is updated according to the structure displacements. Using this technique, along-wind and across-wind effects have been properly explained. The method was applied to an industrial chimney in three scenarios (with or without TMD and for different damping values) and for two wind speeds, showing different responses.