Recent technological advances are based on effectively using complex multiphysics concepts. Therefore, there is an ever increasing need for accurate numerical al-gorithms of reduced complexity for solving multiphysics problems. Traditional mesh-based simulation methods depend on a neighbor connectivity information for formulation of operators like derivatives. In large deformation problems, de-pendence on a mesh could prove a limitation in terms of accuracy and cost of preprocessing. Meshless methods obviate the need to construct meshes thus al-lowing simulations involving severe geometric deformations such as breakup of a contiguous domain into multiple fragments. Smoothed Particle Hydrodynamics (SPH) is a meshless particle based Lagrangian numerical method that has the longest continuous history of development ever since it was introduced in 1977. Commensurate with the significant growth in computational power, SPH has been increasingly applied to solve problems of greater complexity in fluid mechanics, solid mechanics, interfacial flows and astrophysics to name a few. The SPH approximation of the continuity and momentum equations govern-ing fluid flow traditionally involves a stiff equation of state relating pressure and density, when applied to incompressible flow problems. Incompressible Smoothed Particle Hydrodynamics (ISPH) is a variant of SPH that replaces this weak com-pressibility approach with a pressure equation that gives a hydrostatic pressure field which ensures a divergence-free (or density invariant) velocity field. The present study explains the development of an ISPH algorithm and its implementa-tion with focus on application to free surface flows, interaction of fluid with rigid bodies and coupling of incompressible fluids with a compressible second phase. Several improvements to the exiting ISPH algorithm are proposed in this study. A semi-analytic free surface model which is more accurate and robust compared to existing algorithms used in ISPH methods is introduced, validated against experi-ments and grid based CFD results. A surface tension model with specific applica-bility to free surfaces is presented and tested using 2D and 3D simulations. Using theoretical arguments, a volume conservation error in existing particle methods in general is demonstrated. A deformation gradient based approach is used to derive a new pressure equation which reduces these errors. The method is ap-plied to both free surface and internal flow problems and is shown to have better volume conservation and therefore reduced density fluctuations. Also, comments on instabilities arising from particle distributions are made and the role of the smoothing functions in such instabilities is discussed. The challenges in imple-menting the ISPH algorithm in a computer code are discussed and the experience of developing an in-house ISPH code is described. A parametric study on water entry of cylinders of different shapes, angular velocity and density is performed and aspects such as surface profiles, impact pressures and penetration velocities are compared. An analysis on the energy transfer between the solid and the fluid is also performed. Low Froude number water entry of a sphere is studied and the impact pressure is compared with the theoretical estimates. The Incompressible SPH formulation, employing the proposed improvements from the study is then coupled with a compressible SPH formulation to perform two phase flow simulations interacting compressible and incompressible fluids. To gain confidence in its applicability, the simulations are compared against the theoretical predication given by the Rayleigh-Plesset equation for the problem of compressible drop in an incompressible fluid.

This thesis presents a method based on Smoothed Particle Hydrodynamics to simulate sparse particle systems with fluid like properties in real-time. The simulation supports interactions with terrain and objects and is scaled depending on activity of the fluid. We use a carpet method on the GPU to visualize the water surface with translucency, reflection, refraction and added topology. Splash effects and foam are imitated and added as a last step.

This thesis focuses on fluid simulation, particularly on coupling between particle based simulation and grid based simulation and thus modeling evaporation. Mentioned coupling is based on the article Evaporation and Condensation of SPH-based Fluids of authors Hendrik Hochstetter a Andreas Kolb. The goal of this thesis is not purely implementing ideas of the mentioned article, but also study of different methods used for fluid simulation.

Particle simulations in real-time become reality only a few years before, when in computer science occured the idea of GPGPU. This new technology allows use the massive force of graphics card for general purposes. Today, the trend is to accelerate existing algorithms by rewriting into parallel form. On this priciple operate the particle systems too. An interesting area of particle systems are fluid simulations. The simulations are based on the theory of Navier-Stokes equations and their numerical solutions with SPH (Smoothed particle hydrodynamics). Liquids are part of everyday life, and therefore it is important to render them realistically. They are used in modern computer games and different visualizations that run in real time, therefore they must be quickly displayed.

A Computational Fluid Dynamics (CFD) tool is developed in this thesis to solve complex fluid-structure interaction (FSI) problems. The fluid domain is based on Smoothed Particle Hydro-dynamics (SPH) and the structural domain employs large-deformation Finite Element Method (FEM). Validation tests of SPH and FEM are first performed individually. A loosely-coupled SPH-FEM model is then proposed for solving FSI problems. Validation results of two benchmark FSI problems are illustrated (Antoci et al., 2007; Souto-Iglesias et al., 2008). The first test case is flow in a sloshing tank interacting with an elastic body and the second one is dam-break flow through an elastic gate. The results obtained with the SPH-FEM model show good agreement with published results and suggest that the SPH-FEM model is a viable and effective numerical tool for FSI problems. This research is then applied to simulate a two-dimensional free-stream flow interacting with a deformable, pressurized surface, such as an ACV/SES bow seal. The dynamics of deformable surfaces such as the skirt/seal systems of the ACV/SES utilize the large-deformation FEM model. The fluid part including the air inside the chamber and water are simulated by SPH. A validation case is performed to investigate the application of SPH-FEM model in ACV/SES via comparison with experimental data (Zalek and Doctors, 2010). The thesis provides the theory of the SPH and FEM models incorporated and the derivation of the loosely-coupled SPH-FEM model. The validation results have suggested that this SPH-FEM model can be readily applied to skirt/seal dynamics of ACV/SES interacting with free-surface flow. / Ph. D.

Nesse trabalho, foi desenvolvido um simulador numérico (C/C++) para a resolução de escoamentos de fluidos newtonianos incompressíveis, baseado no método de partículas Lagrangiano, livre de malhas, Smoothed Particle Hydrodynamics (SPH). Tradicionalmente, duas estratégias são utilizadas na determinação do campo de pressões de forma a garantir-se a condição de incompressibilidade do fluido. A primeira delas é a formulação chamada Weak Compressible Smoothed Particle Hydrodynamics (WCSPH), onde uma equação de estado para um fluido quase-incompressível é utilizada na determinação do campo de pressões. A segunda, emprega o Método da Projeção e o campo de pressões é obtido mediante a resolução de uma equação de Poisson. No estudo aqui desenvolvido, propõe-se três métodos iterativos, baseados noMétodo da Projeção, para o cálculo do campo de pressões, Incompressible Smoothed Particle Hydrodynamics (ISPH). A fim de validar os métodos iterativos e o código computacional, foram simulados dois problemas unidimensionais: os escoamentos de Couette entre duas placas planas paralelas infinitas e de Poiseuille em um duto infinito e foram usadas condições de contorno do tipo periódicas e partículas fantasmas. Um problema bidimensional, o escoamento no interior de uma cavidade com a parede superior posta em movimento, também foi considerado. Na resolução deste problema foi utilizado o reposicionamento periódico de partículas e partículas fantasmas. / In this work, we have developed a numerical simulator (C/C++) to solve incompressible Newtonian fluid flows, based on the meshfree Lagrangian Smoothed Particle Hydrodynamics (SPH) Method. Traditionally, two methods have been used to determine the pressure field to ensure the incompressibility of the fluid flow. The first is calledWeak Compressible Smoothed Particle Hydrodynamics (WCSPH) Method, in which an equation of state for a quasi-incompressible fluid is used to determine the pressure field. The second employs the Projection Method and the pressure field is obtained by solving a Poissons equation. In the study developed here, we have proposed three iterative methods based on the Projection Method to calculate the pressure field, Incompressible Smoothed Particle Hydrodynamics (ISPH) Method. In order to validate the iterative methods and the computational code we have simulated two one-dimensional problems: the Couette flow between two infinite parallel flat plates and the Poiseuille flow in a infinite duct, and periodic boundary conditions and ghost particles have been used. A two-dimensional problem, the lid-driven cavity flow, has also been considered. In solving this problem we have used a periodic repositioning technique and ghost particles.

Nesse trabalho, foi desenvolvido um simulador numérico (C/C++) para a resolução de escoamentos de fluidos newtonianos incompressíveis, baseado no método de partículas Lagrangiano, livre de malhas, Smoothed Particle Hydrodynamics (SPH). Tradicionalmente, duas estratégias são utilizadas na determinação do campo de pressões de forma a garantir-se a condição de incompressibilidade do fluido. A primeira delas é a formulação chamada Weak Compressible Smoothed Particle Hydrodynamics (WCSPH), onde uma equação de estado para um fluido quase-incompressível é utilizada na determinação do campo de pressões. A segunda, emprega o Método da Projeção e o campo de pressões é obtido mediante a resolução de uma equação de Poisson. No estudo aqui desenvolvido, propõe-se três métodos iterativos, baseados noMétodo da Projeção, para o cálculo do campo de pressões, Incompressible Smoothed Particle Hydrodynamics (ISPH). A fim de validar os métodos iterativos e o código computacional, foram simulados dois problemas unidimensionais: os escoamentos de Couette entre duas placas planas paralelas infinitas e de Poiseuille em um duto infinito e foram usadas condições de contorno do tipo periódicas e partículas fantasmas. Um problema bidimensional, o escoamento no interior de uma cavidade com a parede superior posta em movimento, também foi considerado. Na resolução deste problema foi utilizado o reposicionamento periódico de partículas e partículas fantasmas. / In this work, we have developed a numerical simulator (C/C++) to solve incompressible Newtonian fluid flows, based on the meshfree Lagrangian Smoothed Particle Hydrodynamics (SPH) Method. Traditionally, two methods have been used to determine the pressure field to ensure the incompressibility of the fluid flow. The first is calledWeak Compressible Smoothed Particle Hydrodynamics (WCSPH) Method, in which an equation of state for a quasi-incompressible fluid is used to determine the pressure field. The second employs the Projection Method and the pressure field is obtained by solving a Poissons equation. In the study developed here, we have proposed three iterative methods based on the Projection Method to calculate the pressure field, Incompressible Smoothed Particle Hydrodynamics (ISPH) Method. In order to validate the iterative methods and the computational code we have simulated two one-dimensional problems: the Couette flow between two infinite parallel flat plates and the Poiseuille flow in a infinite duct, and periodic boundary conditions and ghost particles have been used. A two-dimensional problem, the lid-driven cavity flow, has also been considered. In solving this problem we have used a periodic repositioning technique and ghost particles.

Human health is multidimensional: apart from physical, mental, and social aspects, it also incorporates subjective perceptions of health, and functional status (FS). Given that elderly persons have very distinctive health and social needs, multidimensional assessment (MA) of health proves particularly useful in this age group.¶ Aboriginal populations suffer poor health, and there are relatively few studies addressing the health problems of older Aboriginal Australians, mainly because of their distinctive demographic structure, and the low proportion of their elderly. Also, there is no prior information available on MA of health in this Australian population group.¶ This thesis offers a MA of health in older Aboriginal persons from two, urban and rural/isolated, locations in the NT, Katherine and Lajamanu (the NT survey).¶ This thesis specifically addresses the following questions: - what is the physical health, FS, subjective perception of health, and social functioning amongst the NT survey participants? - what are the possible similarities and differences in various dimensions of health between the two major survey locations, what age and gender patterns are observed, and what are the reasons for these patterns, similarities and differences? - how do various dimensions of health relate to each other, and why? - how do current findings relate to broader Aboriginal and non-Aboriginal populations, and why? - what can MA add to a better understanding of various aspects of morbidity and health care use? - what are its possible implications for health planning?¶ Findings from this work indicate poor physical health amongst participants in almost all investigated aspects, comparable to information available from other Aboriginal populations. These are accompanied by low levels of ability for physical functioning. Despite this, subjective perception of health is rather optimistic amongst participants, and levels of social functioning high. Use of health services is mainly related to available health infrastructure. Important health differences exist between Katherine and Lajamanu, and they became particularly visible when all dimensions of health are considered together.¶ The Main conclusions from the current work are that 1) poor physical health is not necessarily accompanied by similar level of deterioration in other dimensions of health: even though participants from the isolated community of Lajamanu experience most chronic diseases, their ability for physical functioning is better, self-perceived health (SPH) more optimistic and levels of social functioning highest 2) institutionalised participants from Katherine suffer by far the worst health of all sample segments in this study; at least some of the poor health outcomes are potentially avoidable, and could be improved by more appropriate residential choices for Aboriginal elderly 3) better health infrastructure does not necessarily bring better health in all its dimensions, suggesting that other factors (primarily socio-economic and cultural) should be addressed in conjunction with this in solving complex health problems of Aboriginal Australians, and 4) it provides strong support that MA can become a useful tool in comprehensive health assessment of older Aboriginals.

Los sistemas flexibles de alta resistencia anclados al terreno son una de las distintas técnicas existentes para la estabilización de taludes, ya sean de roca o de suelos. Están constituidos por una membrana (red de cables o malla de alambre) sujeta al terreno mediante placas de anclaje, cables de refuerzo y bulones. En la mayor parte de los modelos de cálculo existentes se supone un comportamiento activo del sistema, es decir, que evita que se produzcan deslizamientos a través de una supuesta pretensión del sistema y convexidad del terreno. El sistema ejercería una presión normal al terreno que incrementa la tensión tangencial en la superficie potencial de deslizamiento evitando que se alcance la rotura del terreno. En esta tesis se han medido las fuerzas en distintos componentes del sistema desde el momento de la instalación, y se ha comprobado que la pretensión es muy reducida. Por otro lado, la supuesta convexidad del terreno raramente tiene lugar. Todo esto conduce a demostrar la hipótesis de comportamiento pasivo del sistema, es decir, que éste contiene a la masa inestable una vez que se ha producido la rotura. Es por ello que se ha considerado el desarrollar un nueva metodología de cálculo basada en un comportamiento pasivo. La nueva metodología consiste en realizar una simulación numérica dinámica en 2D de la interacción sistema flexible - masa inestable – talud estable. Partiendo de las dimensiones de un círculo de rotura en suelos o una cuña en roca, se deja caer la masa inestable con la fuerza de la gravedad. En su caída, la masa inestable deformará a la membrana, cables de refuerzo y bulones desarrollándose en ellos tensiones que deberán considerarse para su correcto dimensionamiento. En el caso particular de los taludes de suelos, se ha recurrido a la modelización de la masa inestable mediante la discretización por puntos SPH (Smooth Particle Hydrodynamics). / Highly resistant flexible systems anchored to the ground are among the techniques for slope stabilisation, either soil or rock. The system is formed by a membrane (cable net or wire mesh) tightened to the ground through spike plates, reinforcement cables and bolts. In the majority of the existing design models, an active behaviour of the system is considered; which means, that it is able to avoid ground sliding through a pretension of the system and the convexity of the slope surface. The system would exert a normal pressure over the ground that increases the shear stress in the potential slip surface avoiding that failure takes place. In this thesis, forces on different system components have been measured, finding that the pretension force is very low. On the other hand, the supposed ground convexity rarely exists. All this demonstrates that actual system behaviour is passive; which means that it is able to contain the unstable mass once the failure has already occurred. Therefore, a new design methodology based on a passive behaviour has been developed. The new methodology consists in performing a dynamic numerical simulation in 2D of the interaction flexible system – unstable mass – stable slope. Starting from specific known dimensions of slip circle in soils or a wedge in rocks, unstable mass falls only under the action of gravity. During its falling, the unstable mass deforms the membrane, reinforcement cables and bolts. Maximum stresses developed in these components should be considered for their design. In the particular case of soil slopes, unstable mass has been discretised with the mesh free method SPH (Smooth Particle Hydrodynamics).

This thesis makes two noteworthy contributions in the are of mesh-free methods: a Pseudo-Divergence-Free (PDF) Element Free Galerkin (EFG) method which alleviates the volumetric locking and a Stabilized Updated Lagrangian formulation which allows to solve fast-transient dynamic problems involving large distortions. The thesis is organized in the following way. First of all, this thesis dedicates one chapter to the state of the art of mesh-free methods. The main reason is because there are many mesh-free methods that can be found in the literature which can be based on different ideas and with different properties. There is a real need of classifying, ordering and comparing these methods: in fact, the same or almost the same method can be found with different names in the literature. Secondly, a novel improved formulation of the (EFG) method is proposed in order to alleviate volumetric locking. It is based on a pseudo-divergence-free interpolation. Using the concept of diffuse derivatives an a convergence theorem of these derivatives to the ones of the exact solution, the new approximation proposed is obtained imposing a zero diffuse divergence. In this way is guaranteed that the method verifies asymptotically the incompressibility condition and in addition the imposition can be done a priori. This means that the main difference between standard EFG and the improved method is how is chosen the interpolation basis. Modal analysis and numerical results for two classical benchmark tests in solids corroborate that, as expected, diffuse derivatives converge to the derivatives of the exact solution when the discretization is refined (for a fixed dilation parameter) and, of course, that diffuse divergence converges to the exact divergence with the expected theoretical rate. For standard EFG the typical convergence rate is degrade as the incompressible limit is approached but with the improved method good results are obtained even for a nearly incompressible case and a moderately fine discretization. The improved method has also been used to solve the Stokes equations. In this case the LBB condition is not explicitly satisfied because the pseudo-divergence-free approximation is employed. Reasonable results are obtained in spite of the equal order interpolation for velocity and pressure. Finally, several techniques have been developed in the past to solve the well known tensile instability in the SPH (Smooth Particle Hydrodynamics) mesh-free method. It has been proved that a Lagrangian formulation removes completely the instability (but zero energy modes exist). In fact, Lagrangian SPH works even better than the Finite Element Method in problems involving distortions. Nevertheless, in problems with very large distortions a Lagrangian formulation will need of frequent updates of the reference configuration. When such updates are incorporated then zero energy modes are more likely to be activated. When few updates are carried out the error is small but when updates are performed frequently the solution is completely spoilt because of the zero energy modes. In this thesis an updated Lagrangian formulation is developed. It allows to carry out updates of the reference configuration without suffering the appearance of spurious modes. To update the Lagrangian formulation an incremental approach is used: an intermediate configuration will be the new reference configuration for the next time steps. It has been observed that this updated formulation suffers from similar numerical fracture to the Eulerian case. A modal analysis has proven that there exist zero energy modes. In the paper the updated Lagrangian method is exposed in detail, a stability analysis is performed and finally a stabilization technique is incorporated to preclude spurious modes.