Variational Time Integration Approach for Smoothed Particle Hydrodynamics Simulation of Fluids
Leandro Tavares da Silva,Gilson Antonio Giraldi
DOI: https://doi.org/10.48550/arXiv.1512.04444
2018-05-03
Abstract:Variational time integrators are derived in the context of discrete mechanical systems. In this area, the governing equations for the motion of the mechanical system are built following two steps: (a) Postulating a discrete action; (b) Computing the stationary value for the discrete action. The former is formulated by considering Lagrangian (or Hamiltonian) systems with the discrete action being constructed through numerical approximations of the action integral. The latter derives the discrete Euler-Lagrange equations whose solutions give the variational time integrator. In this paper, we build variational time integrators in the context of smoothed particle hydrodynamics (SPH). So, we start with a variational formulation of SPH for fluids. Then, we apply the generalized midpoint rule, which depends on a parameter $\alpha$, in order to generate the discrete action. Then, the step (b) yields a variational time integration scheme that reduces to a known explicit one if $\alpha\in\{0,1\}$ but it is implicit otherwise. Hence, we design a fixed point iterative method to approximate the solution and prove its convergence condition. Besides, we show that the obtained discrete Euler-Lagrange equations preserve linear momentum. In the experimental results, we consider artificial viscous as well as boundary interaction effects and simulate a dam breaking set up. We compare the explicit and implicit SPH solutions and analyze momentum conservation of the dam breaking simulations.
Computational Physics,Fluid Dynamics