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

• Nuclear Engineering and Technology
• /
• v.55 no.12
• /
• pp.4395-4407
• /
• 2023

Bubble columns are widely encountered in several industries, especially in the field of nuclear safety. The Eulerian-Eulerian and the Eulerian-Lagrangian methods are commonly used to investigate bubble columns. Eulerian approaches require additional tasks such as strict volume conservation at the interface and a predefined well-structured grid. In contrast, the Lagrangian approach can be easily implemented. Hence, we introduce a fully Lagrangian approach for the simulation of bubble columns using the discrete bubble model (DBM) and moving particle semi-implicit (MPS) methods. Additionally, we propose a rigorous method to estimate the volume fraction accurately, and verified it through experimental data and analytical results. The MPS method was compared with the experimental data of Dambreak. The DBM was verified by analyzing the terminal velocity of a single bubble for each bubble size. It agreed with the analytical results for each of the four drag correlations. Additionally, the improved method for calculating the volume fraction showed agreement with the Ergun equation for the pressure drop in a packed bed. The implemented MPS-DBM was used to simulate the bubble column, and the results were compared with the experimental results. We demonstrated that the MPS-DBM was in quantitative agreement with the experimental data.

• East Asian mathematical journal
• /
• v.35 no.1
• /
• pp.59-66
• /
• 2019

In this paper we study implicit-explicit (IMEX) methods combined with a semi-Lagrangian scheme to evaluate the prices of fixed strike arithmetic Asian options under jump-diffusion models. An Asian option is described by a two-dimensional partial integro-differential equation (PIDE) that has no diffusion term in the arithmetic average direction. The IMEX methods with the semi-Lagrangian scheme to solve the PIDE are discretized along characteristic curves and performed without any fixed point iteration techniques at each time step. We implement numerical simulations for the prices of a European fixed strike arithmetic Asian put option under the Merton model to demonstrate the second-order convergence rate.

• Proceedings of the Korean Society of Precision Engineering Conference
• /
• 2008.11a
• /
• pp.185-186
• /
• 2008

• Journal of applied mathematics & informatics
• /
• v.31 no.5_6
• /
• pp.795-812
• /
• 2013

The merit function arises from the development of the solution methods for the complementarity problems defined over the cone of non negative real vectors and has been well extended to the complementarity problems defined over the symmetric cones. In this paper, we focus on the extension of the merit functions including the gap function, the regularized gap function, the implicit Lagrangian and others to the complementarity problems defined over the nonsymmetric matrix cone. These theoretical results of this paper suggest new solution methods based on unconstrained and/or simply constrained methods to solve the matrix cone complementarity problems (MCCP).

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

• Journal of the Korean Society of Fisheries and Ocean Technology
• /
• v.32 no.2
• /
• pp.164-176
• /
• 1996

The main intention of this paper is to develop and compare the algorithm based on finite element procedures for nonlinear transient dynamic analysis which has combined effects of material and geometric nonlinearities. Incremental equilibrium equations based on the principle of virtual work are derived by the finite element approach. For the elasto - plastic large deformation analysis of shells and the determination of the displacement-time configuration under time-varying loads, the explicit, implicit and combined explicit-implicit time integration algorithm is adopted. In the time structure is selected and the results are compared with each others. Isoparametric 8-noded quadrilateral curved elements are used for shell structure in the analysis and for geometrically nonlinear elastic behaviour, a total Lagrangian coordinate system was adopted. On the other hands, material nonlinearity is based on elasto-plastic models with Von-Mises yield criteria. Thus, the combined explicit-implicit time integration algorithm is benefit in general case of shell structure, which is the result of this paper.

• Ocean Systems Engineering
• /
• v.7 no.3
• /
• pp.299-318
• /
• 2017

The paper aims at illustrating several key issues and ongoing efforts for development of a reliable fully-Lagrangian particle-based solver for simulation of hydroelastic slamming. Fluid model is founded on the solution of Navier-Stokes along with continuity equations via an enhanced version of a projection-based particle method, namely, Moving Particle Semi-implicit (MPS) method. The fluid model is carefully coupled with a structure model on the basis of conservation of linear and angular momenta for an elastic solid. The developed coupled FSI (Fluid-Structure Interaction) solver is applied to simulations of high velocity impact of an elastic aluminum wedge and hydroelastic slammings of marine panels. Validations are made both qualitatively and quantitatively in terms of reproduced pressure as well as structure deformation. Several remaining challenges as well as important key issues are highlighted. At last, a recently developed multi-scale MPS method is incorporated in the developed FSI solver towards enhancement of its adaptivity.

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

• Bulletin of the Society of Naval Architects of Korea
• /
• v.27 no.3
• /
• pp.107-116
• /
• 1990

The liquid sloshing in an elastic tank is a fluid-structure interaction problem. It requires nonlinear analysis to solve the complicated physics involved in the large interaction of fluid-structure, the variation of dynamic characteristics of structure due to hydrodynamic loading, and the distorsion of fluid flow due to structural vibration. In this paper a Lagrangian FEM is introduced to analyze the liquid sloshing in an elastic tank assuming that the elastic wall is one degree of freedom rigid wall. Numerical integration is performed using an implicit-explicit algorithm, which is formed by mixing the predictor-corrector method and the Runge-Kutta 4th order method. The influence of dynamic characteristics of the sloshing tank on the fluid flow is discussed. The numerical method is also applied for the simulation of the wall generated wave in the tank.