Smoothed Particle Hydrodynamics method and flow dynamics: the case of lava numerical modeling and simulation

Smoothed Particle Hydrodynamics is a Lagrangian mesh-free method that has been gaining momentum in the field of Computational Fluid Dynamics. Thanks to its nature, SPH is able to deal with complex flows and their features, such as non-Newtonian rheologies, free surface, thermal dependency, phase transition, large deformations, etc. SPH shows an intrinsically parallel nature, that allows its execution on high-performance parallel computing hardware, such as modern Graphics Processing Units (GPUs), gaining advantages in terms of simulation time. In this thesis we will work on GPUSPH, an implementation of the SPH method that runs on GPUs. We will study the simulation of a very complex fluid: lava. The combination of free surface, natural topography, phase transition and the formation of structures such as levees and tunnels makes the modeling and simulation of lava flows an extremely challenging task for CFD that has an important impact in numerous fields of engineering and scientific research. We will see the introduction in GPUSPH of models and strategies that deal with the features characterizing lava flows, including the development of a semi-implicit scheme, that allows to simulate very high-viscosity fluids ensuring robustness and reducing simulation times. The new implementation will be tested to verify its correctness and study the accuracy and the performance achieved.

Lo Smoothed Particle Hydrodynamics è un metodo Lagrangiano mesh-free che sta riscuotendo interesse nel campo della fluidodinamica computazionale. Grazie alla sua natura, il metodo SPH è in grado di gestire fluidi complessi caratterizzati da reolgie non-Newtoniane, superfici libere, dipendenze dalla temperatura, transizioni di fase, grandi deformazioni, e così via. Il metodo SPH è caratterizzato da una natura intrinsecamente parallela, che ne rende possibile l'esecuzione su hardware per il calcolo parallelo ad alte prestazioni, come ad esempio le moderne schede grafiche (GPU), ottenendo così tempi di simulazione vantaggiosi. In questa tesi lavoreremo su GPUSPH, un'implementazione su GPU del metodo SPH. Studieremo la simulazione di un fluido molto complesso: la lava. La combinazione di superficie libera, topografia naturale, transizioni di fase e la formazione di strutture come argini e tunnel, rende la modellazione e la simulazione della lava un aspetto di interesse nel campo della fluidodinamica computazionale, con un notevole impatto in numerosi campi dell'Ingegneria e della ricerca scientifica. Vedremo l'introduzione in GPUSPH di modelli e strategie che permetteranno di trattare le diverse caratteristiche delle colate laviche, includendo lo sviluppo di uno schema di integrazione semi-implicito, che permetterà la simulazione di fluidi altamente viscosi, assicurando robustezza e riducendo i tempi di simulazione. La nuova implementazione sarà testata per verificarne la correttezza e studiarne i livelli accuratezza e performance raggiunti.