On the numerical solution to the truncated discrete SPH formulation of the hydrostatic problem
- PDF / 2,052,435 Bytes
- 11 Pages / 595.22 x 842 pts (A4) Page_size
- 78 Downloads / 216 Views
On the numerical solution to the truncated discrete SPH formulation of the hydrostatic problem * Pablo Eleazar Merino-Alonso1, Fabricio Macià1, Antonio Souto-Iglesias2 1. Mathematical Modeling, Analysis, and Simulation Applied to Engineering Research Group (M2ASAI), Universidad politécnica de Madrid (UPM), School of marine engineering (ETSIN), Madrid, Spain 2. Marine and Hydrodynamical Model Basin Research Group (CEHINAV), Universidad Politécnica de Madrid (UPM), School of Marine Engineering (ETSIN), Madrid, Spain (Received May 30, 2020, Revised June 7, 2020, Accepted June 8, 2020, Published August 6, 2020) ©China Ship Scientific Research Center 2020 Abstract: The aim of this work is to study the solution of the smoothed particle hydrodynamics (SPH) discrete formulation of the hydrostatic problem with a free surface. This problem, in which no time dependency is considered, takes the form of a system of linear equations. In particular, the problem in one dimension is addressed. The focus is set on the convergence when both the particle spacing and the smoothing length tend to zero by keeping constant their ratio. Values of this ratio of the order of one, corresponding to a limited number of neighbors, are of practical interest. First, the problem in which each particle has one single neighbor at each side is studied. The explicit expressions of the numerical solution and the quadratic error are provided in this case. The expression of the quadratic error demonstrates that the SPH solution does not converge to the exact one in general under the specified conditions. In this case, the error converges to a residue, which is in general large compared to the norm of the exact solution. The cases with two and three neighbors are also studied. An analytical study in the case of two neighbors is performed, showing how the kernel influences the accuracy of the solution through modifying the condition number of the referred system of linear equations. In addition to that, a numerical investigation is carried out using several Wendland kernel formulas. When two and three neighbors are involved it is found that the error tends in most cases to a small limiting value, different from zero, while divergent solutions are also found in the case of two neighbors with the Wendland Kernel C 2 . Other cases with more neighbors are also considered. In general, the Wendland Kernel
C 2 turns out to be the worst choice, as the solution is divergent for certain values of the ratio between the particle spacing and the smoothing length, associated with an ill-conditioned matrix. Key words: Smoothed particle hydrodynamics (SPH), convergence, water-at-rest problem, hydrostatic tank, free surface, kernel truncation
Introduction The important topic of convergence in smoothed particle hydrodynamics (SPH) has been scarcely treated in the literature, with few theoretical studies[1-4]. The convergence of the SPH semidiscrete scheme for the Euler equation is treated in detail by Franz and Wendland[5]. In that work, an interesting constructive kerne
Data Loading...
