• Wind and Structures
• /
• 제20권3호
• /
• pp.423-448
• /
• 2015

In this paper the unsteady fluid-structure interaction (FSI) problems with large structural displacement are solved by partitioned solution approaches in the arbitrary Lagrangian-Eulerian finite element framework. The incompressible Navier-Stokes equations are solved by the characteristic-based split (CBS) scheme. Both a rigid body and a geometrically nonlinear solid are considered as the structural models. The latter is solved by Newton-Raphson procedure. The equation governing the structural motion is advanced by Newmark-${\beta}$ method in time. The dynamic mesh is updated by using moving submesh approach that cooperates with the ortho-semi-torsional spring analogy method. A mass source term (MST) is introduced into the CBS scheme to satisfy geometric conservation law. Three partitioned coupling strategies are developed to take FSI into account, involving the explicit, implicit and semi-implicit schemes. The semi-implicit scheme is a mixture of the explicit and implicit coupling schemes due to the fluid projection splitting. In this scheme MST is renewed for interfacial elements. Fixed-point algorithm with Aitken's ${\Delta}^2$ method is carried out to couple different solvers within the implicit and semi-implicit schemes. Flow-induced vibrations of a bridge deck and a flexible cantilever behind an obstacle are analyzed to test the performance of the proposed methods. The overall numerical results agree well with the existing data, demonstrating the validity and applicability of the present approaches.

• Coupled systems mechanics
• /
• 제1권1호
• /
• pp.99-113
• /
• 2012

We present a monolithic semi-implicit algorithm for solving fluid-structure interaction problem at small structural displacements. The algorithm uses one global mesh for the fluid-structure domain obtained by gluing the fluid and structure meshes which are matching on the interface. The continuity of velocity at the interface is automatically satisfied and the continuity of stress does not appear explicitly in the global weak form due to the action and reaction principle. At each time step, we have to solve a monolithic system of unknowns velocity and pressure defined on the global fluid-structure domain. Numerical results are presented.

• International Journal of Naval Architecture and Ocean Engineering
• /
• 제11권2호
• /
• pp.639-647
• /
• 2019

A realistic numerical simulation technology using a Lagrangian Fluid-Structure Interaction (FSI) model was combined with a fracture algorithm to predict the fluid-ice-structure interaction. The failure of ice was modeled as the tensile fracture of elastic material by applying a novel FSI model based on the Moving Particle Semi-implicit (MPS) method. To verify the developed fracture algorithm, a series of numerical simulations for 3-point bending tests with an ice beam were performed and compared with the experiments carried out in an ice room. For application of the developed FSI model, a dropping water droplet hitting a cantilever ice beam was simulated with and without the fracture algorithm. The simulation showed that the effects of fracture which can occur in the process of a FSI simulation can be studied.

• 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.

• 한국추진공학회:학술대회논문집
• /
• 한국추진공학회 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.

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

• 대한토목학회논문집
• /
• 제14권5호
• /
• pp.1145-1158
• /
• 1994

물의 유동(流動)에 관한 계산상 효율성(效率性)을 갖는 자유수면(自由水面)을 고려한 3차원(次元) 유한차분(有限差分) 수직모형(數値模型)을 개발하였다. 수직모형(數値模型)은 연직방향(鉛直方向)에 대해 정규화(正規化)한 좌표(座標)(${\sigma}$-coordinate)를 사용하며, 시간(時間) 적분방법(積分方法)으로는 반음해법(半陰解法)(semi-implicit)을 사용하여 계산시간(計算時間)의 효율성(效率性)을 도모하였다. 모드분리(mode-splitting)개념을 도입하여 내부모드(internal mode)에 대해서는 양해법(陽解法)을 사용하였으며, 외부모드(external mode)는 수평방향(水平方向) 운동량방정식(運動量方程式)들과 연속방정식(連續方程式)의 차분식(差分式)으로부터 구한 타원형(楕圓型) 차분방정식(差分方程式)을 SOR방법에 의하여 해석하였다. 이와 같은 방법은 계산(計算) 시간간격(詩間間隔)이 표면(表面) 중력파(重力波)에 대한 CFL(Courant-Fredrich-Lewy)조건에 의해 제약을 받지 않아 계산시간(計算時間)의 효율성을 도모할 수 있다. 개발된 모형은 1차원(次元) 수로(水路)에서 취송유(吹送流)의 연직분포(鉛直分布)에 대한 해석해(解析解)와 비교(比較) 및 연직(鉛直) 가변격자(可變格子)의 도입에 따른 오차분석(誤差分析) 정사각형(正四角形) 호수(湖水)에서 취송유(吹送流) 계산(計算) 및 차분화(差分化) 상수(常數)들의 민감도(敏感度) 분석(分析)을 수행하였다.

• Interaction and multiscale mechanics
• /
• 제5권1호
• /
• pp.27-40
• /
• 2012

The two-dimensional incompressible flow of a linear viscoelastic fluid we considered in this research has rapidly oscillating initial conditions which contain both the large scale and small scale information. In order to grasp this double-scale phenomenon of the complex flow, a multiscale analysis method is developed based on the mathematical homogenization theory. For the incompressible flow of a linear viscoelastic Maxwell fluid, a well-posed multiscale system, including averaged equations and cell problems, is derived by employing the appropriate multiple scale asymptotic expansions to approximate the velocity, pressure and stress fields. And then, this multiscale system is solved numerically using the pseudospectral algorithm based on a time-splitting semi-implicit influence matrix method. The comparisons between the multiscale solutions and the direct numerical simulations demonstrate that the multiscale model not only captures large scale features accurately, but also reflects kinetic interactions between the large and small scale of the incompressible flow of a linear viscoelastic fluid.

• 한국전산유체공학회지
• /
• 제13권3호
• /
• pp.51-61
• /
• 2008

The SMAC (Simplified Marker And Cell) algorithm is extended for an application to thermal non-equilibrium two-phase flows in light water nuclear reactors (LWRs). A two-fluid three-field model is adopted and a multi-dimensional unstructured grid is used for complicated geometries. The phase change and the time derivative terms appearing in the continuity equations are implemented implicitly in a pressure correction equation. The energy equations are decoupled from the momentum equations for faster convergence. The verification of the present numerical method was carried out against a set of test problems which includes the single and the two-phase flows. The results are also compared to those of the semi-implicit ICE method, where the energy equations are coupled with the momentum equation for pressure correction.

• 유체기계공업학회:학술대회논문집
• /
• 유체기계공업학회 1998년도 유체기계 연구개발 발표회 논문집
• /
• pp.27-32
• /
• 1998

This paper presents the development of a two-dimensional model for investigating the fluid flow in water jet and calculating the velocity and pressure distributions. The mathematical formulation as a standard k-$\epsilon$ model was solved employing a general thermofluid-mechanics computer program, PHOENICS code, which is based on the Semi-Implicit Method Pressure Linked Equations(SIMPLE) algorithm. The developed code was applied to water jet design to determine the nozzle size, and investigated the effect of the change of nozzle location. Calculated results showed that the flow pattern is not changed as the change of nozzle location.