• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 2006.05a
• /
• pp.99-104
• /
• 2006

A semi-Lagrangian finite element scheme with objective time stepping algorithm for solving viscoelastic flow problem is presented. The convection terms in the momentum and constitutive equations are treated using a quasi-monotone semi-Lagrangian scheme, in which characteristic feet on a regular grid are traced backwards over a single time-step. Concerned with the generalized midpoint rule type of algorithms formulated to exactly preserve objectivity, we use the geometric transformation such as pull-back, push-forward operation. The method is applied to the 4:1 planar contraction problem for an Oldroyd B fluid for both creeping and inertial flow conditions.

• Journal of computational fluids engineering
• /
• v.9 no.3
• /
• pp.8-17
• /
• 2004

A Semi-Lagrangian method based on CIP(Cubic Interpolated Pseudoparticle)method is proposed and it is applied to solve the two dimensional advection equation. Especially the attentions are given to settle the pole problem and to enhance the accuracy in solving the advection equation on the spherical coordinate system. Tn this algorithm, the CU method is employed as the Semi-Lagrangian method and extended to the spherical coordinate system. To enhance the accuracy of the solution, the spatial discretization is made by CIP method. The mathematical formulation and numerical results are also described. To verify the efficiency, accuracy and capability of proposed algorithm, two dimensional rotating cosine bell problem and the frontogenesis problem are simulated by the present scheme. As results, it is confirmed that the present scheme gives an accurate solution and settles the pole problem in the advection equation on the sphere.

• Journal of the Chungcheong Mathematical Society
• /
• v.28 no.4
• /
• pp.623-633
• /
• 2015

In this paper, a general Weighted Essentially Non-Oscillatory (WENO) interpolation is proposed and applied to a semi-Lagrangian method. The proposed method is based on the conservation law, and characteristic curves are used to complete the semi-Lagrangian method. Therefore, the proposed method satisfies conservation of mass and is free of the CFL condition which is a necessary condition for convergence. Using a several standard examples, the proposed method is compared with the third order Strong Stability Preserving (SSP) Runge-Kutta method to verify the high-order accuracy.

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

• Journal of Ocean Engineering and Technology
• /
• v.16 no.5
• /
• pp.41-48
• /
• 2002

Numerical experiments are conducted using a three-dimensional baroclinic equation model and a Lagrangian method for clarifying the effect of th settling velocity on the suspended solids diffusion caused by the dredging and the reclamination works. Diffusion characteristics of the neutral particles and the weighting particles is experimented by the Lagrangian particles trajectory model, The results show that the diffusion characteristics of the suspended solids is effected by the settling velocity classified by the particles size in the density layered semi-closed bay. To estimate exactly the diffusion characteristics of the suspended solids and the contaminant with weight the three-dimensional baroclinic equation model and the three-dimensional Lagrangian particles trajectory model considering the settling velocity of the particle in the density layered semi-closed bay must be used.

• The KIPS Transactions:PartA
• /
• v.16A no.6
• /
• pp.419-426
• /
• 2009

To realistically simulate fluid, the Navier-Stokes equations are generally used. Solving these Navier-Stokes equations on the Eulerian framework, the non-linear advection terms invoke heavy computation and thus Semi-Lagrangian methods are used as an approximated way of solving them. In the Semi-Lagrangian methods, the locations of advection sources are traced and the physical values at the traced locations are interpolated. In the case of Stam's method, there are relatively many chances of numerical losses, and thus there have been efforts to correct these numerical errors. In most cases, they have focused on the numerical interpolation processes, even simultaneously using particle-based methods. In this paper, we propose a new approach to reduce the numerical losses, through improving the tracing method during the advection calculations, without any modifications on the Eulerian framework itself. In our method, we trace the grids with the velocities which will let themselves to be moved to the current target position, differently from the previous approaches, where velocities of the current target positions are used. From the intuitive point of view, we adopted the simple physical observation: the physical quantities at a specific position will be moved to the new location due to the current velocity. Our method shows reasonable reduction on the numerical losses during the smoke simulations, finally to achieve real-time processing even with enhanced realities.

• Proceedings of the Korean Society of Computer Information Conference
• /
• 2023.07a
• /
• pp.581-583
• /
• 2023

기존의 불 영역 검출을 위한 이미지 전처리 기법은 임계값에 영향을 많이 받기 때문에 사진의 색상 분포에 의존하는 결과가 나온다. 이를 보완하고자 본 논문에서는 화염이 가지고 있는 특유의 비선형적이고 불규칙한 움직임을 추출할 수 있는 방법을 제안한다. 유체 시뮬레이션에서 활용되는 준 라그랑지안(Semi-Lagrangian) 활용하여 유동을 분석하고 이를 기반으로 영상에서 화염의 움직임을 추출한다.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.10 no.2
• /
• pp.73-82
• /
• 1998

In order to simulate surface discharged heat dispersion in costal area, a 2-dimensional Eulerian-Lagrangian model for far field and semi-active particle tracking random walk model in near field has been combined. The mass of discharged heat water in near field has treated as particles with buoyancy and this is eventually converted to horizontal additive dispersion in random walk equations. This model is applied to both a simplified coastal geometry and a real site. In simple application it can simulate plume-like characteristics around discharging point than a near field-model, CORMIX/3. Actual application in the Chonsu Bay shows farther spreading of heat water in near field comparing the observed data, and this shows that the developed model might be applied with satisfaction.

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