• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.9 no.2
• /
• pp.286-290
• /
• 2008

This paper summarizes the components of an explicit aeroelastic solver developed especially for the simulation of dynamic fracture events occurring during the flight of solid propellant rockets. The numerical method combines an explicit Arbitrary Lagrangian Eulerian (ALE) version of the Cohesive Volumetric Finite Element (CVFE) scheme, used to simulate the spontaneous motion of one or more cracks propagating dynamically through a domain with regressing boundaries, and an explicit unstructured finite volume Euler code to follow the flow field during the failure event. A key feature of the algorithm is the ability to adaptively repair and expand the fluid mesh to handle the large geometrical changes associated with grain deformation and crack motion.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2009.04a
• /
• pp.59-62
• /
• 2009

유체-고체 상호작용(FSI)은 산업전반에서 꼭 필요한 분야이면서도 쉽게 접근하기가 어려운 전산역학 분야의 난제 중 하나이다. 유체-고체 상호작용의 전산해석에서 유체와 고체 사이의 불일치 격자망을 어떻게 처리하는가는 매우 어렵고 민감한 부분이 된다. 운동학적 연속성과 계면을 따른 응력의 평형을 추적하기 위해 유체와 고체의 계면에서는 각각의 영역에서 해석된 물리량들을 다른 영역으로 정확히 전달해야 하는데 대부분의 유체-고체 상호작용의 문제들은 불일치 격자를 가지고 있기 때문에 불일치 격자망을 효과적으로 처리하는 수단이 필요하다. 그래서 넓은 분야에 걸쳐 적용 가능한 유체 고체 상호작용 문제에 대한 효과적인 해석방법의 제안이 큰 의미를 갖는다고 생각한다. 따라서 본 연구에서는 유체-고체 계면의 운동을 이동최소제곱 기반의 변절점 요소를 사용하여 모사함으로써 2차원 유체-고체의 상호작용(FSI)을 위한 새로운 접근방법을 제시하였다.

• 한국전산유체공학회:학술대회논문집
• /
• 2009.04a
• /
• pp.116-120
• /
• 2009

A two-dimensional unsteady inviscid flow solver has been developed for the simulation of complex geometric configurations on adaptive Cartesian meshes. Embedded condition was used for boundary condition and a predictor-corrector explicit time marching scheme was used for time-accurate numerical simulation. The Cartesian mesh generator, which was previously developed for steady problem, was used grid generation for unsteady flow. The solver was based on ALE formulation for body motion. For diminishing the effects of cut-cells, the cell merging method was used. Using cell merging method, it was eliminated the CFL constraints. The conservation problem, which is caused cell-type variation around region swept by solid boundary, was also solved using cell merging method. The results are presented for 2D circular cylinder and missile launching problem.

• Journal of computational fluids engineering
• /
• v.10 no.3 s.30
• /
• pp.1-8
• /
• 2005

In this study, unsteady aerodynamic analyses of a wedge-type airfoil based on nonlinear piston theory and Euler equations have been performed in supersonic and hypersonic flows. The third-order nonlinear piston theory (NPT) to calculate unsteady lift and moment coefficients is derived and applied in the time-domain. Also, unsteady flow quantities are obtained from the two-dimensional time-dependent Euler equations. For the CFD based unsteady aerodynamic analyses, an arbitrary Lagrangean-Eulerian (ALE) formulation for the Euler equations is used to calculate flow fluxes in the computational flow field with moving boundaries. Numerical comparisons for unsteady lift and moment coefficients are presented between NPT and Euler approaches. The results show very good agreements in the high supersonic and hypersonic flows. It means that the present NPT can be efficiently used to predict unsteady aerodynamic forces ol wedge type airfoils with dynamic motions in the high supersonic and hypersonic flow regimes.

• Wind and Structures
• /
• v.18 no.3
• /
• pp.215-233
• /
• 2014

This paper presents a method to extract flutter derivatives of bridge decks based on a combination of the computational fluid dynamics (CFD), system simulations and system identifications. The incompressible solver adopts an Arbitrary Lagrangian-Eulerian (ALE) formulation with the finite volume discretization in space. The imposed sectional motion in heaving or pitching relies on exponential time series as input, with aerodynamic forces time histories acting on the section evaluated as output. System identifications are carried out to fit coefficients of the inputs and outputs of ARMA models, as to establish discrete-time aerodynamic models. System simulations of the established models are then performed as to obtain the lift and moment exerting on the sections to a sinusoidal displacement. It follows that flutter derivatives are identified. The present approaches are applied to a hexagon thin plate and a real bridge deck. The results are compared to the Theodorsen closed-form solution and those from wind tunnel tests. Satisfactory agreements are observed.

• Proceedings of the KSME Conference
• /
• 2008.11a
• /
• pp.442-447
• /
• 2008

This paper presents a new approach to simulate fluid-solid interaction problems involving non-matching interfaces. The coupling between fluid and solid domains with dissimilar finite element meshes consisting of 4-node quadrilateral elements is achieved by using the interface element method (IEM). Conditions of compatibility between fluid and solid meshes are satisfied exactly by introducing the interface elements defined on interfacing regions. Importantly, a consistent transfer of loads through matching interface element meshes guarantees the present method to be an efficient approach of the solution strategy to fluid-solid interaction problems. An arbitrary Lagrangian-Eulerian (ALE) description is adopted for the fluid domain, while for the solid domain an updated Lagrangian formulation is considered to accommodate finite deformations of an elastic structure. The stabilized equal order velocity-pressure elements for incompressible flows are used in the motion of fluids. Fully coupled equations are solved simultaneously in a single computational domain. Numerical results are presented for fluid-solid interaction problems involving nonmatching interfaces to demonstrate the effectiveness of the methodology.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.27 no.9
• /
• pp.1327-1334
• /
• 2003

Flow characteristics of the fluid laden with many particles in the two-dimensional channel are investigated using the Navier-Stokes equations coupled with the equation of motion of particles by direct numerical simulation. A four-step fractional step method with Crank-Nicolson scheme and ALE technique is used for P2P1 mixed finite element method. The motion and distribution of particles in the fluid is virtually described as a result of direct numerical simulation and the increase of viscosity is compared with theoretical equations. The effect of channel height on the relative viscosity and the tubular pinch effect are discussed.

• Proceedings of the Korean Institute of Navigation and Port Research Conference
• /
• v.1
• /
• pp.181-186
• /
• 2006

In this paper, generation mechanisms of ocean freak waves are briefly introduced in the context of wave-current interaction phenomena. The present model of the fluid motion is based on the Navier-Stokes equations incorporating velocity-pressure formulation because of need to model the nonlinear wave interaction with spatially non-uniform current field. In order to deal with the free surface motion, an Arbitrary Lagrangian-Eulerian (ALE) description is adopted. As an accurate and efficient numerical tool, the spectral element method is presented with general features and specific treatment for the wave-current interaction problem. As an intermediate stage of development, solution procedure and characteristics aspects of the present modeling and numerical method are addressed in detail, and preliminary numerical results prove its accuracy and convergence.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.32 no.5
• /
• pp.323-331
• /
• 2019

In this study, a computational model approach for the modeling of hydrodynamic pressures acting on radial gates during strong earthquakes is proposed. The use of the dynamic layering method with the Arbitrary Lagrangian Eulerian (ALE) algorithm and the SIMPLE method for simulating free reservoir surface flow in addition to moving boundary interfaces between the fluid domain and a structure due to earthquake excitation are suggested. The verification and validation of the proposed approach are realized by comparisons performed using the renowned formulation derived by the experimental results for vertical and inclined dam surfaces subjected to earthquake excitation. A parameter study for the truncated lengths of the two-dimensional fluid domain demonstrates that twice the water level leads to efficient and converged computational results. Finally, numerical simulations for large radial gates with different curvatures subjected to two strong earthquakes are successfully performed using the suggested computational model.

• Structural Engineering and Mechanics
• /
• v.77 no.4
• /
• pp.441-461
• /
• 2021

Reinforced concrete (RC) columns are crucial in building structures and they are of higher vulnerability to terrorist threat than any other structural elements. Thus it is of great interest and necessity to achieve a comprehensive understanding of the possible responses of RC columns when exposed to high intensive blast loads. The primary objective of this study is to derive analytical formulas to assess vulnerability of RC columns using an advanced numerical modelling approach. This investigation is necessary as the effect of blast loads would be minimal to the RC structure if the explosive charge is located at the safe standoff distance from the main columns in the building and therefore minimizes the chance of disastrous collapse of the RC columns. In the current research, finite element model is developed for RC columns using LS-DYNA program that includes a comprehensive discussion of the material models, element formulation, boundary condition and loading methods. Numerical model is validated to aid in the study of RC column testing against the explosion field test results. Residual capacity of RC column is selected as damage criteria. Intensive investigations using Arbitrary Lagrangian Eulerian (ALE) methodology are then implemented to evaluate the influence of scaled distance, column dimension, concrete and steel reinforcement properties and axial load index on the vulnerability of RC columns. The generated empirical formulae can be used by the designers to predict a damage degree of new column design when consider explosive loads. With an extensive knowledge on the vulnerability assessment of RC structures under blast explosion, advancement to the convention design of structural elements can be achieved to improve the column survivability, while reducing the lethality of explosive attack and in turn providing a safer environment for the public.