Cette thèse porte sur le développement d’une nouvelle technique de modélisation des problèmes IFS utilisant les méthodes particulaires. Ce travail s’inscrit dans la continuité des travaux de recherche de l’équipe biomécanique du LAMIH, concernant la compréhension du comportement de l’os humain dans son environnement de moelle osseuse. La méthode SPH a été utilisée pour la modélisation des travées osseuses, supposées dans une première approche comme des solides élastiques. La méthode LB a été développée pour la modélisation des écoulements de moelle considérée comme un fluide visqueux incompressible. L’efficacité et la performance de ces deux méthodes ont été démontrées grâce aux benchmarks académiques évalués et les résultats comparés à ceux de la littérature ou ceux obtenus par des logiciels commerciaux. A l’issue d’une revue de l’état de l’art des techniques de couplage fluide-structure, une approche partitionnée en temps a été choisie, permettant d’utiliser deux codes distincts basés sur des algorithmes de résolution de type dynamique explicite. La discrétisation spatiale est faite par une technique spécifique basée sur les domaines fictifs, cette technique est très efficace car elle ne nécessite pas de rediscrétisation des domaines. L’approche de couplage développée a été appliquée à des benchmarks académiques ainsi qu’à une application en biomécanique, ayant permis d’aboutir à des résultats numériques satisfaisants. Plusieurs pistes d’amélioration sont maintenant nécessaires afin d’aller vers des modélisations plus biofidèles telles que la prise en compte du contact et de l’endommagement. / The objective of this thesis is the development of a new technique for the FSI problems modelling using particulars methods. This work is in the continuity of the LAMIH biomechanics team research works, regarding the comprehension of behavior of bone in its environment of marrow. The SPH method was used for the trabeculae modelling, supposed in a first attempt as an elastic solid. The LB method was developed for the marrow flow modelling considered as a viscous incompressible liquid. The efficacy and performance of these two methods were demonstrated using academics benchmarks which were evaluated and the results were compared of those of literature and of those obtained from commercials softwares. Following a bibliographic review of FSI coupling techniques, a partitioned approach in time was chosen, allowing the use of two separates codes, both based on a dynamic explicit algorithm resolution scheme. The special discretization was done based on a specific technique of fictional domain, this technique is very efficient because it doesn’t require an additional domain discretization. The coupling approach developed was applied on academic benchmarks and on a biomechanical application, leading to satisfactory numerical results. Many Improvement track are now necessary to go towards more biofidelic modeling as taking into account the contact and the damage.

Desde a década de 70, há um crescente interesse em simulações em computador de fenômenos físicos visto sua diversidade de aplicações. Dentre esses fenômenos, podem ser destacados a interação entre corpos rígidos, elásticos, plásticos, quebráveis e também fluidos. Neste trabalho realizamos a simulação de um desses fenômenos, o escoamento de fluidos, por um método conhecido como Smoothed Particles Hydrodynamics, uma abordagem lagrangeana baseada em partículas para resolução das equações que modelam o movimento do fluido. Várias são as vantagens de métodos lagrangeanos usando partículas sobre os que usam malhas, por exemplo, as propriedades do material transladam com as partículas como função do tempo, além da capacidade de lidar com grandes deformações. Dentre as desvantagem, destacamos uma deficiência relacionada ao ganho de energia total do sistema e estabilidade das partículas. Para lidar com isso, utilizamos uma abordagem baseada na lei da conservação da energia: em um sistema isolado a energia total se mantém constante e ela não pode ser criada ou destruida. Dessa forma, alterando o integrador temporal nós restringimos o aumento arbitrário de energia, tornando a simulação mais tolerante às condições iniciais / Since the late 70s, there is a growing interest in physically-based simulations due to its increasing range of application. Among these simulations, we may highlight interaction between rigid, elastic, plastic and breakable bodies and also fluids. In this work, one of these phenomena, fluid flow, is simulated using a technique known as Smoothed Particle Hydrodynamics, a meshless lagrangean method that solves the equations of the flow behavior of fluids. There are several advantages of meshless methods over mesh-based methods, for instance, the material properties are translated along with particles as a function of time and the ability to handle arbitrary deformations. Among the disadvantages, we may highlight a problem related to the gain of energy by the system and stability issues. In order to handle this, we used an approach based on the law of conservation of energy: in an isolated system the total energy remains constant and cannot be created or destroyed. Based on this, we used a technique that bounds the total energy and the simulation becomes less sensitive to initial conditions

Lorsqu’un matériau solide est soumis à un chargement dynamique (par l’impact d’un projectile, la détonation d’un explosif ou l’irradiation par un laser intense), il se forme une onde de choc, qui se propage dans le matériau depuis la surface chargée. Si cette onde débouche sur une surface libre comportant des défauts géométriques tels que des rugosités, des rayures ou des cavités, son interaction avec ces défauts conduit à l’éjection, sous forme de jets de matière, de débris dont la taille caractéristique est de l’ordre du micromètre et dont la vitesse est typiquement de quelques km/s. La maîtrise de ce processus, appelé microjetting ou micro-éjection, est essentielle pour de nombreuses applications (conception de blindage, découpe pyrotechnique, usinage à très haute vitesse, expériences de Fusion par Confinement Inertiel…). Dans le cadre de cette thèse, menée en collaboration avec le centre CEA de Bruyères-le-Châtel, ce phénomène est étudié dans quatre métaux (Aluminium, Etain, Cuivre et Plomb) à partir de rainures calibrées de deux types : triangulaires isolées de demi-angles d’ouverture contrôlés (20°, 30° et 45°) ou sinusoïdales périodiques. Les influences du matériau, de la forme et de l’ouverture des défauts, de la pression de choc et de l’état du milieu (solide ou fondu sous choc ou en détente) sur les propriétés balistiques des éjectas (vitesses de jet, distribution en taille et densité surfacique des débris constituant les jets) sont évaluées via trois approches complémentaires : expérimentale, théorique et numérique.L’étude expérimentale comporte plusieurs campagnes de chocs laser, effectuées sur l’installation LULI2000 du Laboratoire pour l’Utilisation des Lasers Intenses (Ecole Polytechnique, Palaiseau), avec plusieurs techniques de diagnostic : Ombroscopie Transverse, Vélocimétrie Hétérodyne, radiographie X rapide in-situ, récupération d’éjectas dans des gels (analysés ensuite en microtomographie). Les résultats sont confrontés à des prédictions théoriques (hydrodynamique des chocs obliques et des charges creuses pour les rainures triangulaires, instabilités de Richtmyer-Meshkov pour les rainures sinusoïdales). Enfin, les simulations numériques réalisées avec le code Radioss utilisent deux approches complémentaires : les Eléments Finis Lagrangiens et la formulation SPH (Smoothed Particles Hydrodynamics), encore très peu appliquée au microjetting, plus empirique que la première mais mieux adaptée aux grandes déformations dans les jets et permettant d’accéder à des distributions de tailles de fragments / When a dense material is subjected to a dynamic load (such as projectile impact, explosive detonation or irradiation by a high energy laser beam), a shock wave propagates from the loaded surface. If this shock wave interacts with a free surface with geometrical defects such as grooves, scratches or cavities, it can lead to the ejection of micrometric debris with typical velocities of a few km/s. Understanding this microjetting process is a key issue for many applications, including shielding design, pyrotechnics, high-speed machining and Inertial Confinement Fusion experiments.In this work in collaboration with the CEA-DIF at Bruyères-le-Châtel, this phenomenon is studied under laser-driven shock loading in four materials (Aluminum, Tin, Copper and Lead) with calibrated grooves of two types: isolated triangular profile with controlled aperture half-angles (20°, 30° and 45°) or periodic sinusoidal shape. The influences of the material, of the geometry of the defects, of the shock pressure and of the state of matter (solid or melted under shock or release wave) on the ballistic properties of the ejecta (jet velocity, size distribution and areal mass of the debris constituting the jet) are investigated with three complementary approaches: experimental, theoretical and numerical.The experimental study involves several campaigns performed at the LULI2000 facility of the Laboratoire pour l’Utilisation des Lasers Intenses (Ecole Polytechnique, Palaiseau) and complementary diagnostic techniques: Transverse Shadowgraphy, Heterodyne Velocimetry, fast in situ X-ray radiography, recovery of the ejecta in a gel followed by microtomography. The results are compared with theoretical predictions (2D shocks and shaped charges hydrodynamics for the triangular grooves, Richtmyer-Meshkov Instabilities for the sinusoidal grooves). Then, numerical simulations are performed with the Radioss code with two complementary approaches: the Lagrangian Finite Elements and the SPH (Smoothed Particles Hydrodynamics) formulation, still very scarcely applied to microjetting, more empirical than the first approach but more suitable to the high strains in the jets and allowing access to size distributions of the debris.

Desde a década de 70, há um crescente interesse em simulações em computador de fenômenos físicos visto sua diversidade de aplicações. Dentre esses fenômenos, podem ser destacados a interação entre corpos rígidos, elásticos, plásticos, quebráveis e também fluidos. Neste trabalho realizamos a simulação de um desses fenômenos, o escoamento de fluidos, por um método conhecido como Smoothed Particles Hydrodynamics, uma abordagem lagrangeana baseada em partículas para resolução das equações que modelam o movimento do fluido. Várias são as vantagens de métodos lagrangeanos usando partículas sobre os que usam malhas, por exemplo, as propriedades do material transladam com as partículas como função do tempo, além da capacidade de lidar com grandes deformações. Dentre as desvantagem, destacamos uma deficiência relacionada ao ganho de energia total do sistema e estabilidade das partículas. Para lidar com isso, utilizamos uma abordagem baseada na lei da conservação da energia: em um sistema isolado a energia total se mantém constante e ela não pode ser criada ou destruida. Dessa forma, alterando o integrador temporal nós restringimos o aumento arbitrário de energia, tornando a simulação mais tolerante às condições iniciais / Since the late 70s, there is a growing interest in physically-based simulations due to its increasing range of application. Among these simulations, we may highlight interaction between rigid, elastic, plastic and breakable bodies and also fluids. In this work, one of these phenomena, fluid flow, is simulated using a technique known as Smoothed Particle Hydrodynamics, a meshless lagrangean method that solves the equations of the flow behavior of fluids. There are several advantages of meshless methods over mesh-based methods, for instance, the material properties are translated along with particles as a function of time and the ability to handle arbitrary deformations. Among the disadvantages, we may highlight a problem related to the gain of energy by the system and stability issues. In order to handle this, we used an approach based on the law of conservation of energy: in an isolated system the total energy remains constant and cannot be created or destroyed. Based on this, we used a technique that bounds the total energy and the simulation becomes less sensitive to initial conditions

Made available in DSpace on 2016-12-23T14:04:30Z (GMT). No. of bitstreams: 1 Wesley Jose Nunes Pinto.pdf: 2090367 bytes, checksum: e676fde8423a3a2cfeac61da24020ea8 (MD5) Previous issue date: 2013-08-19 / Neste estudo foi aplicado do método numérico, sem malhas, baseado em partículas, denominado SPH (Smoothed Particles Hydrodynamics). E um código numérico na linguagem computacional FORTRAN foi utilizado para solucionar as equações de Navier-Stokes. O clássico problema da literatura da dinâmica dos fluidos Computacional, denotado como problema da cavidade quadrada bidimensional (Shear-Driven Cavity Flow) , foi estudado com a intenção de verificar o comportamento do código numérico em relação a resultados específicos já existentes do assunto. O citado problema físico das cavidades abertas é amplamente empregado como benchmark, visando a validação do método numérico utilizado no trabalho desenvolvido na pesquisa. O trabalho de análise e validação do código numérico foi dividido em três seções: a primeira lista as localizações dos centros dos vórtices principais gerados pelo escoamento na aresta superior das cavidades; a segunda plota os perfis das componentes das velocidades centrais das cavidades; e a terceira: lista os desvios absolutos dos perfis das velocidades centrais do presente trabalho, comparados com dados de outros estudos. Constata-se que o método SPH apresentou boa acurácia nas simulações realizadas, obtendo boa concordância entre os resultados das simulações dinâmicas com os dados de referências, validando-se o modelo numérico proposto, tendo melhores resultados para baixos números de Reynolds / In this study, it was applied the numerical method, grid-free, based on particles named SPH (Smoothed Particles Hydrodynamics). Also, a numerical code in the computer language FORTRAN was used to solve the Navier-Stokes Equations. This classic problem of the literature related to Computational Fluid Dynamics indicated as Shear-Driven Cavity Flow was studied to check the behavior of the numerical code regarding specific existing results. Such problem is highly used as Benchmark, aiming the validation of the numerical method used to develop the research. The analysis and validation of the numerical code was divided into three sections: the first one lists the location of the centre of the main vortex generated by the flow of the upper edge of the cavities; the second one plots the profiles of the components of the central speed of the cavities; the third one lists the absolute deviation of the profiles of the central speed of this study compared with other cases data. It is established that the SPH Method presented accuracy in the performed simulations, in a consonance between the results of the dynamic simulations and the reference data, thus the proposed numerical model was validated with better results for low Reynolds numbers