• 스마트미디어저널
• /
• 제4권2호
• /
• pp.9-16
• /
• 2015

Fluids appear in innumerable phenomena; therefore, it is interesting to reproduce those phenomena by computer graphics techniques. However, this process is not trivial. We work with a fluid simulation that uses Navier-Stokes equations to model the fluid, a semi-Lagrangian approach to solve it and the level set method to track the surface of the fluid. Modified versions of the Navier-Stokes equations for computer graphics allow us to create a wide diversity of effects. In this paper, we propose a technique that allows us to integrate a force inspired by surface tension into the model. We describe which information we need and how to modify the model with this new approach. We end up with a modified simulation that has additional effects that might be suitable for computer graphics purposes. The effects that we are able to recreate are small waves and droplet-like formations close to the surface of the fluid. This model preserves the overall behavior governed by the Navier-Stokes equations.

• 한국경영과학회지
• /
• 제25권1호
• /
• pp.27-36
• /
• 2000

We perform a computational complexity analysis of a heuristic algotithm proposed in the literature for the combinatorial optimization problems extended with a single side-constraint. This algorithm, although such a view was not given in the original work, is a disguised version of an optimal Lagrangian dual solution technique. It also has been observed to be a very efficient heuristic producing near-optimal solutions for the primal problems in some experiments. Especially, the number of iterations grows sublinearly in terms of the network node size so that the heuristic seems to be particularly suitable for the applicatons such as routing with semi-real time requirements. The goal of this paper is to establish a polynomal worst-case complexity of the algorithm. In particular, the obtained complexity bound suports the sublinear growth of the required iterations.

• 한국해양공학회지
• /
• 제31권6호
• /
• pp.388-396
• /
• 2017

A Lagrangian approach based computational fluid dynamics (CFD) was used to simulate large and/or sharp deformations and fragmentations of interfaces, including free surfaces, through tracing each particle with physical quantities. According to the concept of the particle-based CFD method, it is possible to apply it to both fluid particles and solid particles such as sand, gravel, and rock. However, the presence of more than two different phases in the same domain can make it complicated to calculate the interaction between different phases. In order to solve multiphase problems, particle interaction models for multiphase problems, including surface tension, buoyancy-correction, and interface boundary condition models, were newly adopted into the moving particle semi-implicit (MPS) method. The newly developed MPS method was used to simulate a typical validation problem involving dam breaking. Because the soil and other particles, excluding the water, may have different viscosities, various viscosity coefficients were applied in the simulations for validation. The newly developed and validated MPS method was used to simulate the mobile beds induced by broken dam flows. The effects of the viscosity on soil particles were also investigated.

• 한국해양공학회지
• /
• 제27권4호
• /
• pp.90-97
• /
• 2013

Oil booms are one of the most widely used types of equipment for the protection of coastal areas against oil spills. In some situations, however, there are several types of oil leaks from the oil boom. Important factors regarding these phenomena include the surrounding ocean environment, such as waves, the density and viscosity of oil, the length of the oil boom skirt, etc. To estimate the performance of the oil boom, it is necessary to predict the behavior of the spilled oil and oil boom. In the present study, the prediction of oil boom performance in waves was carried out using the Pusan-National-University-modified Moving Particle Semi-implicit (PNU-MPS) method, which is an improved version of the original MPS proposed by Koshizuka and Oka (1996). The governing equations, which consist of continuity and Navier-Stokes equations, are solved by Lagrangian moving particles, and all terms expressed by differential operators in the governing equations are replaced by the particle interaction models based on a kernel function. The simulation results were validated through a comparison with the results of Violeau et al. (2007)..

• Journal of Mechanical Science and Technology
• /
• 제19권spc1호
• /
• pp.388-394
• /
• 2005

Recently, the authors have developed a method for real-time dynamics of multibody systems, which combines a semi-recursive formulation to derive the equations of motion in dependent relative coordinates, along with an augmented Lagrangian technique to impose the loop closure conditions. The following numerical integration procedures, which can be grouped into the so-called structural integrators, were tested : trapezoidal rule, Newmark dissipative schemes, HHT rule, and the Generalized-${\alpha}$ family. It was shown that, for large multi body systems, Newmark dissipative was the best election since, provided that the adequate parameters were chosen, excellent behavior was achieved in terms of efficiency and robustness with acceptable levels of accuracy. In the present paper, the performance of the described method in combination with another group of integrators, the Implicit Runge-Kutta family (IRK), is analyzed. The purpose is to clarify which kind of IRK algorithms can be more suitable for real-time applications, and to see whether they can be competitive with the already tested structural family of integrators. The final objective of the work is to provide some practical criteria for those interested in achieving real-time performance for large and complex multibody systems.

• 대한조선학회지
• /
• 제23권2호
• /
• pp.1-13
• /
• 1986

The nonlinear hydrodynamic forces acting on a two-dimensional circular cylinder, oscillating with large amplitude in the free surface, are calculated by using the Semi-Lagrangian Time-Step-ping Method used by O.M. Faltinsen. In present calculation the position and the potential value of free surface are calculated using the exact kinematic and dynamic free surface boundary condition. At each time step an integral equation is solved to obtain the value of potential and normal velocity along the boundaries, consisting of both the body surface and the free surface. Some effort was devoted to the elimination of instability arising in the range of high frequency. Numerical simulations were performed up to the 3rd or 4th period which seems to be enough for the transient effect to die out. Each harmonic component and time-mean force are obtained by the Fourier transform of forces in time domain. The results are compared with others' experimental and theoretical results. Particularly, the calculation shows the tendency that the acceleration-phase 1st-harmonic component(added mass) increases as the motion amplitude increases and a reverse tendency in the velocity-phase 1st-harmonic component(damping coefficient). The Yamashita's experimental result also shows the same tendency. In general, the present result show relatively good agreement with the Yamashita's experimental result except for the time-mean force.

• Nuclear Engineering and Technology
• /
• 제52권6호
• /
• pp.1172-1179
• /
• 2020

As a Lagrangian particle method, Moving Particle Semi-implicit (MPS) method has great capability to capture interface/surface. In recent years, the multiphase flow simulation using MPS method has become one of the important directions of its developments. In this study, some key methods for multiphase flow have been introduced. The interface tension model in multiphase flow is modified to maintain the smooth of the interface and suitable for the three-phase flow. The mass transfer at immiscible liquid interface entrained by single bubble which could occur in Molten Core-Concrete Interaction (MCCI) has been investigated using this particle method. With the increase of bubble size, the height of entrainment column also increases, but the time of film rupture is slightly different. With the increase of density ratio between the two liquids, the height of entrained column decreases significantly due to the decreasing buoyancy of the denser liquid in the lighter liquid. In addition, the larger the interface tension coefficient is, the more rapidly the entrained denser liquid falls. This study validates that the MPS method has shown great performance for multiphase flow simulation. Besides, the influence of physical parameters on the mass transfer at immiscible interface has also been investigated in this study.

• 대한조선학회논문집
• /
• 제28권1호
• /
• pp.38-52
• /
• 1991

자유표면하에서 움직이는 임의의 형상의 3차원 물체로 인한 비선형 유체력을 준 Lagrangian 방법을 사용하여 해석하였다. 경계치 문제는 경계 적분 방법(Boundary Integral Method)을 이용하여 해결하였으며, 물체와 자유 표면의 형상은 곡면 Panel로 표현하였다. 이들 표면은 hi-cubic B-spline 방법을 사용하여 유한개의 작은 표면 요소로 나뉘어지게 되며, 또한 \phi와 (equation omitted) 표면 요소상에서 bi-linear하다고 가정한다. 특이점에 의한 유기 포텐시얼의 계산시 1/R에 비례하는 부분은 제거하고 해석적으로 처리하였다. 물체로부터 멀리 떨어진 곳에서의 유체 유통은 좌표계의 원점에 위치한 Dipole로 표현하였으며, Time Stepping시 Runge-Kutta의 4차 방법을 사용하였다. 3차원인 경우에 대한 적분 방정식과 포텐시얼의 시간 미분간에 대한 경계 조건이 유도되었으며, 이러한 식들을 사용하여 자유표면의 형상과 물체의 운동을 동시에 계산하였다. 대진폭을 가지고 규칙적으로 진동하는 물체에 작용하는 힘과 이때의 자유 표면 형상을 계산하고 기 발표된 자료와 비교하여 보았으며, 자유표면 근처에서 운동하는 물체에 작용하는 비선형 효과를 관찰하였다.

• Nuclear Engineering and Technology
• /
• 제42권6호
• /
• pp.620-635
• /
• 2010

In this paper we describe current activities on the project Multi-Scale Modeling and Analysis of convective boiling (MSMA), conducted jointly by the Paul Scherrer Institute (PSI) and the Swiss Nuclear Utilities (Swissnuclear). The long-term aim of the MSMA project is to formulate improved closure laws for Computational Fluid Dynamics (CFD) simulations for prediction of convective boiling and eventually of the Critical Heat Flux (CHF). As boiling is controlled by the competition of numerous phenomena at various length and time scales, a multi-scale approach is employed to tackle the problem at different scales. In the MSMA project, the scales on which we focus range from the CFD scale (macro-scale), bubble size scale (meso-scale), liquid micro-layer and triple interline scale (micro-scale), and molecular scale (nano-scale). The current focus of the project is on micro- and meso-scales modeling. The numerical framework comprises a highly efficient, parallel DNS solver, the PSI-BOIL code. The code has incorporated an Immersed Boundary Method (IBM) to tackle complex geometries. For simulation of meso-scales (bubbles), we use the Constrained Interpolation Profile method: Conservative Semi-Lagrangian $2^{nd}$ order (CIP-CSL2). The phase change is described either by applying conventional jump conditions at the interface, or by using the Phase Field (PF) approach. In this work, we present selected results for flows in complex geometry using the IBM, selected bubbly flow simulations using the CIP-CSL2 method and results for phase change using the PF approach. In the subsequent stage of the project, the importance of effects of nano-scale processes on the global boiling heat transfer will be evaluated. To validate the models, more experimental information will be needed in the future, so it is expected that the MSMA project will become the seed for a long-term, combined theoretical and experimental program.

• Steel and Composite Structures
• /
• 제45권2호
• /
• pp.175-191
• /
• 2022

A review of previous studies shows that although there is a considerable difference between buckling loads of structures under follower and non-follower lateral loads, only the buckling load of FGM elliptical cylindrical shell under non-follower lateral load was investigated in the literature. This study is the first to obtain the buckling load of elliptical FGM cylindrical shells under follower lateral load and also make a comparison between buckling loads of elliptical FGM cylindrical shells under follower and non-follower lateral loads. Moreover, this research is the first one to derive the load potential function of elliptical cylindrical shell. In this regard, the FGM cylindrical elliptical shell was modeled using the semi-analytical finite strip method and based on the First Shear Deformation Theory (FSDT). The shell is discretized by strip elements aligned in the longitudinal direction. The Lagrangian and harmonic shape functions were considered in the circumference and longitudinal directions, respectively. The buckling pressure of the shell under follower and non-follower lateral loads was obtained from eigenvalue problem. The results obtained from the model were compared with those presented in the literature to evaluate the validity of the model. A comparison index was defined to compare the buckling loads of the shell under follower and non-follower lateral load. A parametric study was carried out to investigate the effects of material properties and shell geometry characteristics on the comparison index. For the elliptical cylindrical shells with length-to-radius ratio greater than 16 and major-to-minor axis ratio greater than 0.6, the comparison index reaches to more than 20 percent which is significant. Moreover, the maximum difference is about 30 percent in some cases. The results obtained from the parametric study indicate that the buckling load of long elliptical cylindrical shell under non-follower load is not reliable.