Tesis:
A novel generalized diffusive SPH model: Theoretical analysis and 3D HPC implementation
- Autor: CERCÓS PITA, Jose Luis
- Título: A novel generalized diffusive SPH model: Theoretical analysis and 3D HPC implementation
- Fecha: 2016
- Materia: Sin materia definida
- Escuela: E.T.S. DE INGENIEROS AERONAUTICOS
- Departamentos: MECANICA DE FLUIDOS Y PROPULSIÓN AEROESPACIAL
- Acceso electrónico: http://oa.upm.es/43276/
- Director/a 1º: SOUTO IGLESIAS, Antonio
- Resumen: The subject of the present thesis is the development of a new free Smoothed Particle Hydrodynamics (SPH) code oriented to the researchers, conversely to the already available free/open source SPH codes. The new software should be competitive in performance terms, allowing at the same time to be modified and extended with a minimum amount of effort. Such code should be later applied to analyze SPH models where a diffusive term is added in the mass conserving equation, which are usually known as δ-SPH models. Therefore, the thesis can be split in 3 main parts: The theoretical aspects of the SPH models considered, the numerical model implementation, and their verification and validation. During the theoretical analysis, the possibility of adding a diffusive term inside the continuity equation in order to significantly reduce the characteristic numerical instabilities of the SPH model is investigated, compiling a list of the already existing models, discussing their features and relationships, and pointing out their main benefits and drawbacks. A close relation between the δ-SPH models and the Riemann solver based ones have been unveiled. Also, a new diffusive term not requiring neither a tuning parameter nor a Courant condition is presented. Such term is resulting from a slightly modification of an already existing one. Later, as part of the same theoretical investigation of the model, an energy conservation analysis has been carried out, considering the interactions with the boundary conditions, such that several extra energy terms are resulting, some of them representing extra energy dissipations if a set of conditions are fulfilled. All this theoretical aspects have been implemented in a new free SPH solver, accelerated with OpenCL, modular, and extensible by Python scripts. Implementation details regarding the software itself, as well as the SPH models considered have been provided. The main benefits of the new code, compared with other already existing free/open source alternatives, have been discussed. Finally, a set of 4 verification and validation test cases is described, covering all the theoretical and implementation aspects discussed along the thesis. In such verification and validation cases the implementation details, performance, and results quality are discussed. The verification and validation cases are complemented with two practical applications, where the capabilities of the new software to scale to significantly complex problems are shown. RESUMEN En esta tesis se documenta el desarrollo de un nuevo software de SPH (Smoothed Particle Hydrodynamics), que tiene la particularidad de estar orientado al trabajo de investigación, en lugar de enfocarse a aplicaciones industriales como los otros códigos libres/abiertos disponibles en la actualidad. El nuevo código debe ser competitivo en términos de coste computacional, al tiempo que provee una herramienta modular y extensible. Dicho código se empleará más tarde para analizar los modelos conocidos como δ-SPH, que se caracterizan por la presencia de un término difusivo en la equación de conservación de masa. Por tanto, el trabajo llevado a cabo durante el desarrollo de la tesis se puede dividir en tres partes: Los aspectos teóricos, la implementación del modelo numérico, y la verificación y validación de los mismos. A lo largo del análisis teórico se investigará la posibilidad de añadir un término difusivo a la ecuación de continuidad, con el objetivo de reducir significativamente las inestabilidades numéricas que vienen caracterizando al modelo de SPH. Para ello se recopilarán todos los modelos existentes actualmente, discutiendo sus propiedades principales, las relaciones que guardan entre ellos, y las principales ventajas a destacar de cada uno de ellos. Ello nos servirá para establecer una relación entre los modelos δ-SPH, y los modelos basados en Riemann solvers. Además se presentará un nuevo término, resultante de la modificación de uno ya existente. Dicho término tendrá la ventaja de no requerir ni un parámetro de tuneado, ni una condición de Courant adicional sobre el paso de tiempo. Como parte del mismo análisis teórico, se investigarán las propiedades de conservación de los modelos, en donde se considerará la interacción con los contornos. Como resultado, se demostrará la existencia de unos términos de energía adicionales, algunos de los cuales pueden ser considerados disipaciones extra si se cumplen unas determinadas condiciones. Todos estos aspectos teóricos se implementan en un nuevo código libre de SPH, el cual está acelerado con OpenCL, además de ser modular y extensible con códigos de Python. Se documentan algunos detalles sobre la implementación del software en si mismo, y sobre la implementación de los modelos SPH a considerar. El nuevo software se compara con otras alternativas ya existentes, resaltando sus ventajas. Finalmente, se describen 4 casos de verificación y validación que permiten cubrir todos los aspectos teóricos y computacionales abordados durante la tesis. En ellos se proveen algunos detalles sobre cómo se han desarrollado, el coste computacional, y la calidad de los resultados. Estos casos de verificación y validación se complementan con dos aplicaciones prácticas, donde se muestran las capacidades del software para afrontar casos significativamente complejos.