• Journal of the Korean Society of Visualization
• /
• v.20 no.3
• /
• pp.74-85
• /
• 2022

We propose a Cubic-Interpolated Pseudo-Particle Lattice Boltzmann method (CIP-LBM) for the convection-diffusion equation (CDE) based on the Bhatnagar-Gross-Krook (BGK) scheme equation. The CIP-LBM relies on an accurate numerical lattice equilibrium particle distribution function on the advection term and the use of a splitting technique to solve the Lattice Boltzmann equation. Different schemes of lattice spaces such as D1Q3, D2Q5, and D2Q9 have been used for simulating a variety of problems described by the CDE. All simulations were carried out using the BGK model, although another LB scheme based on a collision term like two-relation time or multi-relaxation time can be easily applied. To show quantitative agreement, the results of the proposed model are compared with an analytical solution.

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

• 한국전산유체공학회:학술대회논문집
• /
• 2010.05a
• /
• pp.109-114
• /
• 2010

Interface tracking of two phase is significant to analyze multi-phase phenomena. The VOF(Volume of Fluid) and level set are well known interface tracking method. However, they have limitations to solve compressible flow and incompressible flow at the same time. CIP(Cubic Interpolate Propagation) method is appropriate for considering compressible and incompressible flow at once by solving the governing equation which is divided up into advection and non-advection term. In this article, we analyze the droplet impingement according to various We number using improved CIP method which treats nonlinear term once more comparison with original CIP method. Furthermore, we compare spread radius after droplet impingement on the wall with the experimental data and original CIP original CIP method, and it reduces the mass conservation error which is generated in the numerical analysis comparison with original CIP method.

• Journal of Environmental Health Sciences
• /
• v.31 no.6
• /
• pp.489-497
• /
• 2005

A numerical method for the mathematical water modeling in 2-dimensional flow has been developed. The model based on a split operator technique, in which, the advection term is calculated using the upwind scheme. The diffusion term is one- dimensionalized and calculated using Crank-Nicholson's implicit finite difference scheme to reduce the numerical errors from large time steps and variable spacings. It also provides a relatively simple and economic method for more accurate simulation of pollutant dispersion. Water depths and flow velocities in the Boreyong reservoir during the normal water periods were predicted by numerical experiments with a 2-dimensional flow model so as to provide current field data for the study of advection and diffusion of pollutants. Developed 2-dimensional water quality model is applied to Boreyong reservoir to simulate a spatial and periodical changes of water quality.

• Journal of Korea Water Resources Association
• /
• v.45 no.2
• /
• pp.179-189
• /
• 2012

To simulate the transport of pollutant, a numeric model for the advection-diffusion equation with the decay term is developed. This is finite-difference model using the implicit method (with the weight factor ${\alpha}$) and Gauss-Seidel SOR(successive over-relaxation). This model is compared to the analytical solutions (of simpler dimensional or boundary conditions), and in the condition of Peclet number < 5~20, the result shows stable condition, and Crank-Nicolson method (${\alpha}$=0.5) shows the more accurate results than fully-implicit method (${\alpha}$=1). The mass of advection, diffusion and decay is calculated and the error of mass balance is less than 3%. This model can evaluate the 3-D concentrations of the advection-diffusion and decay problems, but this model uses only the finite-difference method with the fixd grid system, so it can be effectively used in the problems with small Peclet numbers like the pollutant transport in groundwater.

• Journal of the Korea Computer Graphics Society
• /
• v.15 no.1
• /
• pp.1-6
• /
• 2009

This paper presents a new CIP method for unstructured mesh to reduce the numerical dissipation. To reflect precise physical characteristics, CIP method updates both the physical quantity and the derivative information. The proposed method uses the Finite Volume Method(FVM) to solve the non-advection term of CIP equation. And we performed several experiments to improve the accuracy of third-order interpolation. Our result shows that our algorithm has less numerical dissipation than that of linear advection solver.

• Proceedings of the Korean Environmental Sciences Society Conference
• /
• 2003.11a
• /
• pp.165-170
• /
• 2003

A significant surface net heat loss appears around the Kuroshio and the Tsushima Warm Current regions. The area where the surface heat loss occurs should require heat to be supplied by the current to maintain the long-term annual heat balance. Oceanic heat advection in these regions plays an important role in the heat budget. The spatial distribution of the heat supply by the Tsushima Warm Current near the surface was examined by calculating the horizontal heat supply in the surface layer of the East/Japan Sea, directly from historical sea surface temperature and current data. We have also found a simple estimation of the effective vertical scale of heat supply by the current to compensate net heat loss using the heat supplied by the current in the surface 10m layer. The heat supplied by the current for the annual heat balance was large in the Korea/Tsushima Strait and along the Japanese Coast, and was small in the northwestern part of the East/Japan Sea. The amount of heat supplied by the current was large in the northwestern part and small in the southeastern part of the East/Japan Sea. These features suggest that the heat supplied by the Tsushima Warm Current is restricted to near the surface around the northeastern part and extends to a deeper layer around the southeastern part of the East/Japan Sea.

• Journal of computational fluids engineering
• /
• v.15 no.4
• /
• pp.25-31
• /
• 2010

Interface tracking of two phase is significant to analyze multi-phase phenomena. The VOF(Volume of Fluid) and level set are well known interface tracking method. However, they have limitations to solve compressible flow and incompressible flow at the same time. CIP(Cubic Interpolate Propagation) method is appropriate for considering compressible and incompressible flow at once by solving the governing equation which is divided up into advection and non-advection term. In this article, we analyze the droplet impingement according to various We number using improved CIP method which treats nonlinear term once more comparison with original CIP method. Furthermore, we compare spread radius after droplet impingement on the wall with the experimental data and original CIP method. The result using improved CIP method shows the better result of the experiments, comparison with result of original CIP method, and it reduces the mass conservation error which is generated in the numerical analysis comparison with original CIP method.

• 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 구현을 통해 벡터 연산 등의 효율성을 높이며 시뮬레이션의 속도를 향상시킨다.

• Journal of Astronomy and Space Sciences
• /
• v.25 no.4
• /
• pp.415-424
• /
• 2008

Kwak et al. (2008) found that the mean neutral wind pattern in the high-latitude lower thermosphere is dominated by rotational flow than by divergent flow. As an extension of the our previous work (Kwak et al. 2008), we performed a term analysis of vorticity equation that describes the driving forces for the rotational component of the horizontal wind in order to determine key processes that causes strong rotational flow in the high-latitude lower thermospheric winds. For this study the National Center for Atmospheric Research Thermosphere-Ionosphere Electrodynamics General Circulation Model (NCAR-TIEGCM) is used. The primary forces that determine variations of the vorticity are the ion drag term and the horizontal advection term. Significant contributions, however, can be made by the stretching term. The effects of IMF on the vorticity forces are seen down to around 105-110km.