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

• Nuclear Engineering and Technology
• /
• 제55권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
• /
• 제35권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.

• 한국정밀공학회:학술대회논문집
• /
• 한국정밀공학회 2008년도 추계학술대회 논문집
• /
• pp.185-186
• /
• 2008

• Journal of applied mathematics & informatics
• /
• 제31권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).

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

• 수산해양기술연구
• /
• 제32권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
• /
• 제7권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
• /
• 제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.

• 대한조선학회지
• /
• 제27권3호
• /
• pp.107-116
• /
• 1990

구조물의 탄성을 고려한 탱크내 유동은 자유수면을 갖는 유체와 탄성변형하는 구조물이 연성된 시스템으로서 유체유동으로 인한 과도한 구조물변형, 유체의 부가질량 및 부가감쇠력에 의한 구조물의 동특성변화, 구조물 진동으로 인한 유체유동의 왜곡 등이 복합된 비선형 해석이 요구된다. 본 논문에서는 탱크 벽을 1자유도 수평운동하는 강체로 가정하였으며 Lagrangian 유한요소법을 이용하여 유동해석을 수행하였고 유체-구조물 연성문제의 수치적분을 위하여 조합된 implicit-explicit 알고리듬을 도입하였다. 탱크벽의 동특성 변화에 따른 유체-구조물연성 탱크의 동특성변화를 관찰하였으며 파도생성 문제에 관한 수치계산을 수행하였다.