DOI QR코드

DOI QR Code

Semi-Lagrangian 이류항 계산의 추적법 개선

Improved Trajectory Calculation on the Semi-Lagrangian Advection Computation

  • 박수완 (경북대학교 컴퓨터공학과) ;
  • 백낙훈 (경북대학교 전자전기컴퓨터학부) ;
  • 유관우 (경북대학교 컴퓨터공학과)
  • 발행 : 2009.12.31

초록

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

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.

키워드

참고문헌

  1. E. Wu, Y. Liu, and X. Liu, “An improved study of real-time fluid simulation on GPU,” Comp. Anim. Virtual Worlds, pp.139-146, 2004
  2. N. Foster and D. Metaxas, “Modeling the Motion of Hot, Turbulent Gas,” SIGGRAPH'99, pp.181-188, 1997
  3. J. Stam, “Stable fluids,” SIGGRAPH'99 , pp.121-128, 1999
  4. P. Bartello and S. J. Thomas, “The Cost-Effectiveness of Semi-Lagrangian Advection,” Monthly Weather Review, pp.2883-2897, 1996
  5. R. Fedkiw, J. Stam, and H. W. Jensen, “Visual simulation of smoke,” SIGGRAPH'01, pp.15-22, 2001
  6. O.-Y. Song, H. Shin, and H.-S. KO, “Stable but Non-Dissipative Water,” ACM Transactions on Graphics, pp.81-97, 2005
  7. D. Kim, O.-Y. Song, and H.-S. KO, “A Semi-Lagrangian CIP Fluid Solver without Dimensional Splitting,” Proceedings of EUROGRAPHICS 2008, pp.467-475, 2008
  8. T. YABE and T. AOKI, “A universal solver for hyperbolic equations by cubic-polynomial interpolation i. one dimensional solver” Computer Physics Communications 66, 2-3, 219-232, 1991 https://doi.org/10.1016/0010-4655(91)90071-R
  9. T. YABE, T. ISHIKAWA, P. Y. WANG, T. AOKI, Y. KADOTA, and F. IKEDA, “A universal solver for hyperbolic equations by cubic-polynomial interpolation ii. Two and three-dimensional solvers” Computer Physics Communications 66, 2-3, 233-242. 1991 https://doi.org/10.1016/0010-4655(91)90072-S
  10. B.-M. Kim, Y. Liu, I. Llamas, and J. Rossignac, “Using BFECC for fluid simulation” In Eurographics Workshop on Natural Phenomena, 2005
  11. B.-M. Kim, Y. Liu, I. Llamas, and J. Rossignac, “Advections with significantly reduced dissipation and diffusion,” IEEE Trans. Vis. Comput. Graph. 13(1), pp.135-144, 2007 https://doi.org/10.1109/TVCG.2007.3
  12. A. Selle , R. Fedkiw, B.-M. Kim, Y. Liu and J.Rossignac, “An unconditionally stable MacCormack method,” Journal of Scientific Computing, pp.350-371, 2008
  13. D. Enright, S. Marschner, and R. Fedkiw, “Animation and rendering of complex water surfaces,” ACM Trans. Graph. 21, pp.736-744, 2002 https://doi.org/10.1145/566570.566645
  14. A. Selle, N. Rasmussen, and R. Fedkiw, “A vortex particle method for smoke, water and explosions,” ACM Trans. Graph. 24, 2005 https://doi.org/10.1145/1073204.1073282
  15. O.-Y. Song, D. KIM, H.-S. KO, “Derivative particles for simulating detailed movements of fluids,” IEEE Transactions on Visualization and Computer Graphics 13, pp.711-719, 2007 https://doi.org/10.1109/TVCG.2007.1022
  16. J. Stam, “Real-Time Fluid Dynamics for Games,” Proceedings of the Game Developer Conference, 2003