• Proceedings of the Korean Society of Computer Information Conference
• /
• 2024.01a
• /
• pp.407-410
• /
• 2024

본 논문에서는 Semi-Lagrangian 이류 과정에서 역추적(Backward tracing)한 위치의 주변 속도를 Divergence-constrained MLS(Moving least squares)를 이용하여 보간하고 그 결과를 이류된 속도 데이터의 외력으로 적용해 연기 시뮬레이션의 난류 표현을 개선할 수 있는 새로운 프레임워크를 제안한다. 일반적인 MLS는 고차보간법이기 때문에 시간에 따른 연속성 보장이 안 되기 때문에 그 결과가 노이즈한 형태로 나타난다. 본 논문에서는 연기의 원본 속도와 제안하는 기법을 통해 생성된 속도 간의 각도 변화를 통해 난류를 생성하며 이를 통해 안정적이고 연기의 밀도를 이류시킨다.

• The KIPS Transactions:PartA
• /
• v.14A no.4
• /
• pp.191-196
• /
• 2007

In the field of computer graphics, Navier-Stokes equations would be used for realistic simulations of smokes and currents. However, implementations derived from these equations are hard to achieve for real-time simulations, mainly due to its massive and complex calculations. Thus, there have been various attempts to approximate these equations for real-time simulation of smokes and others. When the advection terms of the equations are approximated by the Semi-Lagrange methods, the fluid density can be rapidly reduced and small-scale vorticity phenomena are easy to be missed, mainly due to the numerical losses over time. In this paper, we propose an improved numerical method to approximately calculate the advection terms, and thus eliminate these problems. To calculate the advection terms, our method starts to set critical regions around the target grid points. Then, among the grid points in a specific critical region, we search for a grid point which will be advected to the target grid point, and use the velocity of this grid point as its advection vector. This method would reduce the numerical losses in the calculation of densities and vorticity phenomena, and finally can implement more realistic smoke simulations. We also improve the overall efficiency of vector calculations and related operations through GPU-based implementation techniques, and thus finally achieve the real-time simulation.

• KSII Transactions on Internet and Information Systems (TIIS)
• /
• v.12 no.6
• /
• pp.2595-2618
• /
• 2018

In this paper, we investigate secure communication with the presence of multiple eavesdroppers (Eves) in a two-tier downlink dense heterogeneous network, wherein there is a macrocell base station (MBS) and multiple femtocell base stations (FBSs). Each base station (BS) has multiple users. And Eves attempt to wiretap a macrocell user (MU). To keep Eves ignorant of the confidential message, we propose a physical-layer security scheme based on cross-layer cooperation to exploit interference in the considered network. Under the constraints on the quality of service (QoS) of other legitimate users and transmit power, the secrecy rate of system can be maximized through jointly optimizing the beamforming vectors of MBS and cooperative FBSs. We explore the problem of maximizing secrecy rate in both non-colluding and colluding Eves scenarios, respectively. Firstly, in non-colluding Eves scenario, we approximate the original non-convex problem into a few semi-definite programs (SDPs) by employing the semi-definite relaxation (SDR) technique and conservative convex approximation under perfect channel state information (CSI) case. Furthermore, we extend the frame to imperfect CSI case and use the Lagrangian dual theory to cope with uncertain constraints on CSI. Secondly, in colluding Eves scenario, we transform the original problem into a two-tier optimization problem equivalently. Among them, the outer layer problem is a single variable optimization problem and can be solved by one-dimensional linear search. While the inner-layer optimization problem is transformed into a convex SDP problem with SDR technique and Charnes-Cooper transformation. In the perfect CSI case of both non-colluding and colluding Eves scenarios, we prove that the relaxation of SDR is tight and analyze the complexity of proposed algorithms. Finally, simulation results validate the effectiveness and robustness of proposed scheme.

• Journal of the Korea Concrete Institute
• /
• v.20 no.3
• /
• pp.317-323
• /
• 2008

In the current study, A model for the flow analysis of fresh and highly flowable concrete has been developed using a particle method, the moving particle semi-implicit (MPS) method. The phenomena on the flow of concrete has been considered as a visco-plastic flow problem, and the basic governing equation of concrete particle dynamics has been based on the Navier-Stokes equation in Lagrangian form and the conservation of mass. In order to formulate a visco-plastic flow constitutive law of fresh concrete, concrete is modeled as a highly viscous material in the state of non-flow and as a visco-plastic material in the state of flow after reaching the yield stress of fresh concrete. A flow test of fresh concrete in the L-box was simulated and the predicted flow was well matched with the experimental result. The developed method was well showed the flow motion of concrete particles because it was formulated to be based on the motion of visco-plastic fluid dynamics.

• Journal of Ocean Engineering and Technology
• /
• v.27 no.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 Ocean Engineering and Technology
• /
• v.31 no.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.

• Journal of computational fluids engineering
• /
• v.17 no.2
• /
• pp.35-41
• /
• 2012

This paper presents incompressible Navier-Stokes solution algorithm for 2D Free-surface flow problems on the Cartesian mesh, which was implemented to run on Graphics Processing Units(GPU). The INS solver utilizes the variable arrangement on the Cartesian mesh, Finite Volume discretization along Constrained Interpolation Profile-Conservative Semi-Lagrangian(CIP-CSL). Solution procedure of incompressible Navier-Stokes equations for free-surface flow takes considerable amount of computation time and memory space even in modern multi-core computing architecture based on Central Processing Units(CPUs). By the recent development of computer architecture technology, Graphics Processing Unit(GPU)'s scientific computing performance outperforms that of CPU's. This paper focus on the utilization of GPU's high performance computing capability, and presents an efficient solution algorithm for free surface flow simulation. The performance of the GPU implementations with double precision accuracy is compared to that of the CPU code using an representative free-surface flow problem, namely. dam-break problem.

• Journal of Mechanical Science and Technology
• /
• v.19 no.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.

• Bulletin of the Society of Naval Architects of Korea
• /
• v.23 no.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.

• Proceedings of the Korean Information Science Society Conference
• /
• 2006.10a
• /
• pp.143-147
• /
• 2006

자연 현상에서 나타나는 연기나 난류의 움직임을 사실적으로 시뮬레이션을 할 때 Navier-Stokes 방정식을 이용한다. 이 방정식을 이용한 구현은 방대한 연산량과 계산의 복잡성으로 인하여 실시간 시뮬레이션이 어렵다. 이 때문에 실시간 처리를 위하여 복잡한 수식을 근사화한다. 유체 시뮬레이션의 이류(advect) 과정에서 근사화를 위해 Semi-Lagrangian 방법을 이용할 때, 연기 시뮬레이션은 시간이 지남에 따라 밀도가 현저히 줄어들고 소규모의 소용돌이(small-scale vorticity) 현상이 급격히 감소하는 등의 수치적 소실이 발생한다. 본 논문에서는 이 문제를 해결하기 위해 이류항(advection term)을 계산할 때 새로운 수치적 방법을 제안한다. 본 논문에서는 이류항의 값을 구할 때, 현재 격자 주변의 값 중에서 다음 단계에 현재 격자의 위치로 오는 속도를 가진 격자를 찾아, 그 격자의 속도를 이류 속도 벡터로 활용한다. 이는 밀도와 소용돌이 현상의 수치적 소실을 줄여서 사실성을 높이고 실시간 처리도 가능하게 한다. 또한 본 논문에서는 GPU 구현을 통해 벡터 연산 등의 효율성을 높이며 시뮬레이션의 속도를 향상시킨다.