• 한국전산유체공학회지
• /
• 제9권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.

• 충청수학회지
• /
• 제28권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.

• 한국소성가공학회:학술대회논문집
• /
• 한국소성가공학회 2006년도 춘계학술대회 논문집
• /
• 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.

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

• 한국해양공학회지
• /
• 제16권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.

• 정보처리학회논문지A
• /
• 제16A권6호
• /
• pp.419-426
• /
• 2009

일반적으로 사실성 있는 유체를 시뮬레이션하기 위해 Navier-Stokes 방정식을 사용한다. Euler 구조에서 Navier-Stokes 방정식을 풀 때, 이 류항은 비선형이어서 계산이 복잡하기 때문에 근사화한 모델로 Semi-Lagrangian 방법을 사용한다. Semi-Lagrangian 방법에서는 먼저 이류하 는 위치를추적하고, 추적한 위치에서 값을 보간해서 사용한다. Stam이 제안한 방법으로 계산할 경우, 이 과정에서 수치적 소실이 많이 발생하 기 때문에 수치적 소실을 보정하려는 노력들이 있어 왔다. 그러나 대부분의 경우에 보간하는 과정에서의 소실을 줄이려는 노력이거나, 입자를 같이 사용하는 방법이었다. 따라서 본 논문에서는 Euler 구조에서 다른 추가나 변형을 가하지 않고 이류항의 연산에서 추적법을 개선함으로 수 치적 소실을 줄이는 방법을 제안한다. 우리의 방법에서는 현재 격자의 속도로 역추적하는 기존의 방법이 아니라, 현재의 격자로 오게 될 속도 를 가진 격자를 찾아서, 그 격자의 물리량들을 선형 보간하여 사용한다. 이는 직관적으로 생각할 때, 어느 지점의 물리량은 그 지점의 속도로 인해 다음 단계에 다른 지점에 있게 된다는 사실을 그대로 적용한 것이다. 본 논문에서 제안한 방법으로 기체를 시뮬레이션 했을 때 수치적 소 실이 줄었으며, 그로 인해 사실성을 높이면서도 실시간 처리가 가능했다.

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

• 대한기계학회논문집A
• /
• 제29권11호
• /
• pp.1455-1462
• /
• 2005

In this paper, we intend to introduce a nonlinear finite element method based on the fully nonlinear potential flow theory in order to simulate the large amplitude sloshing flow in two-dimensional baffled tank subject to horizontally forced excitation. The free surface is tracked by a direct time differentiation scheme with the four-step predictor-corrector time integration method. The flow velocity is accurately recovered from the velocity potential by second-order least square method. In order to maintain the finite element mesh regularity and total mass, the semi-Lagrangian surface tracking method with area conservation is applied. According to the numerical formulae, we perform the parametric experiments by varying the installation height and the opening width of baffles, in order to examine the effects of baffle on the nonlinear liquid sloshing. From the numerical results, the hydrodynamic characteristics of the large amplitude sloshing are investigated.

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

• 한국전산구조공학회:학술대회논문집
• /
• 한국전산구조공학회 2002년도 봄 학술발표회 논문집
• /
• pp.369-376
• /
• 2002

The stochastic analysis of semi-infinite domain is presented using the weighted integral method, which is improved to include the higher order terms in expanding the displacement vector. To improve the weighted integral method, the Lagrangian remainder is taken into account in the expansion of the status variable with respect to the mean value of the random variables. In the resulting formulae only the 'proportionality coefficients' are introduced in the resulting equation, therefore no additional computation time and memory requirement is needed. The equations are applied in analyzing the semi-infinite domain. The results obtained by the improved weighted integral method are reasonable and are in good agreement with those of the Monte Carlo simulation. To model the semi-infinite domain, the Bettess's infinite element is adopted, where the theoretical decomposition of the strain-displacement matrix to calculate the deviatoric stiffness of the semi-infinite domains is introduced. The calculated value of mean and the covariance of the displacement are revealed to be larger than those given by the finite domain assumptions which is thought to be rational and should be considered in the design of structures on semi-infinite domains.