• Wind and Structures
• /
• v.20 no.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
• /
• v.1 no.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
• /
• v.11 no.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
• /
• v.16 no.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.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2008.03a
• /
• 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.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.23 no.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.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.14 no.5
• /
• pp.1145-1158
• /
• 1994

A three-dimensional, finite difference, numerical model with free surface was developed on ${\sigma}$-coordinate. A semi-implicit numerical scheme in time has been adopted for computational efficiency. The scheme is essentially independent of the stringent stability criteria (CFL condition) for explicit schemes of external surface gravity wave. Implicit algorithm was applied for vertical shear stress, Coriolis force and pressure gradient terms. The reliability of the model with vertically variable grid system was checked by the comparison of simulation results with analytic solution of wind-driven currents in a one-dimensional channel. Sensitivity analysis of differencing parameters was carried out by applying the model to the calculation of wind-driven currents in a square lake.

• Interaction and multiscale mechanics
• /
• v.5 no.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.

• Journal of computational fluids engineering
• /
• v.13 no.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.12a
• /
• 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.