Cette thèse vise à développer et évaluer la méthode de Galerkin discontinue (DG) pour la simulationdes grandes échelles (LES) des écoulements turbulents. L'approche DG présente un nombre d'avantages intéressants pour la LES : ordre élevé, stencil compact, prise en compte des maillages non structurés et expression de la solution numérique dans une base de polynômes permettant l'utilisation de modèles de turbulence multi-échelle. Parmi ce type de modèles, nous nous sommes intéressés ici à la méthode Variational Multiscale (VMS) qui consiste à séparer les échelles résolues dans la base de polynômes pour restreindre l'influence du modèle à une gamme réduite d'échelles. Les modèles considérés ont été paramétrés en prenant en compte les fonctions de transfert spécifiques aux discrétisations DG. La précision de la méthode pour la représentation de phénomènes turbulents variés a été évaluée à travers la réalisation de DNS de configurations académiques. Enfin, l'approche VMS/DGa été éprouvée sur des configurations simples à haut nombre de Reynolds. Il apparaît que cette méthodologie permet la représentation précise des phénomènes turbulents pour un coût réduit en terme de degrés de liberté.

L'objet de cette thèse est de modéliser et simuler des écoulements diphasique/triphasique à phase non miscibles. L'impact de plusieurs gouttes d' émulsion (eau/huile) sur une plaque mouillante l'huile est simulée. Ainsi, une méthode pour lisser l'interface (SVOF), afin d'obtenir une courbure précise, basée sur une méthode eulérienne de type "Volume Of Fluid" (VOF) spécifique au caractère multiphasique de l' écoulement est développée. Un modèle de ligne triple et un modèle de mouillabilité sont développés pour calculer les forces capillaires. Ces modèles et ces méthodes sont validées partir de données expérimentales puis utilisées pour simuler le Plate-Outet les écoulements macroscopiques au voisinage du cylindre de laminage / The scope of this dissertation is to model and simulate non-miscible multiphase °ows. Theimpact of several emulsion droplet on the wetting steel strip is simulated. So, the method ofsmooth VOF based on Eulerian "Volume Of Fluid" approach which is particulary adapted tointerfacial °ows is developed. The new method SVOF allows to calculate the curvature with abetter precison than other method. A wettability model and a triple line model are developedto calculate the capillary forces. This models and this methods validated on the experimentaldata and used to simulate the Plate-Out and the macroscopic °ows in neighbourhood of coldrolling system.

Mémoire numérisé par la Division de la gestion de documents et des archives de l'Université de Montréal.

La présence d'inclusions d'oxyde au sein d'aubes de turbine haute pression monocristallines en alliage AM1 est suspectée d'être à l'origine de la germination de grains parasites responsables d'une diminution critique des performances mécaniques des aubes. Afin de caractériser l'évolution de la propreté inclusionnaire de l'alliage AM1 au cours des différentes étapes du procédé fonderie par induction sous vide, une étude expérimentale et numérique a été réalisée. La partie expérimentale combine d'une part des visualisations in situ du radeau inclusionnaire à la surface du bain de métal en cours de fusion et d'autre part la réalisation du test du bouton sur un ensemble d'échantillons de matière prélevés aux différentes étapes du procédé de fonderie, dans le but de caractériser la population d'inclusions (nature chimique, distribution de taille, fraction volumique). Les résultats mettent en évidence une modification de la population inclusionnaire majoritaire au cours du procédé en raison d'une alimentation régulière du bain de métal en particules d'oxyde issues du creuset de fusion. La partie numérique s'appuie sur un ensemble de trois modèles permettant la modélisation successive des forces de brassage électromagnétiques, de l'écoulement turbulent du métal liquide et de la trajectographie des particules d'oxydes au sein du bain métallique. La surface du bain de métal est considérée déformée et non évolutive au cours du temps. Le modèle de trajectographie tient compte de la turbulence de l'écoulement et de différentes formes de particules. Les résultats numériques montrent une bonne concordance avec les tendances dégagées expérimentalement / The presence of oxide inclusions within single-crystal turbine blade high pressure in AM1 involves the germination of parasite grains succeeding the critical reduction of their mechanical performances.In order to evaluate the cleanliness of material AM1 and to follow its evolution during various stages of the process foundry (VIM), an experimental and a numerical study were carried out. The experimental part consists in following the evolution of the inclusion population of the raw material until the liquid bath by combining visualizations of the surface of the liquid bath and the application of the electron beam button melting test on a set of material extracted during the process of foundry. The results highlight a modification of the major inclusion population during the process and a regular supply of particles within the liquid bath which come from the crucible. The digital part is based on a set of three models allowing the modeling of the electromagnetic, hydrodynamic phenomena and the trajectory calculation of the particles within the liquid bath. The value of the intensity of the electric current was given from a calibration of the electromagnetic model based to experimental measures of the axial magnetic induction within the industrial furnace. For the whole of the mathematical models, the surface of the liquid bath selected is deformed and non-evolutive. The model of trajectory calculation is based on a Lagrangian approach by taking account of the turbulence of the flow. The digital results show a good agreement with in experiments released tendencies

Il est probable qu'à terme les émissions de nanoparticules soient réglementées et ce sont donc les concentrations en nombre qui seront considérées. Il convient donc d'adapter les modèles afin de pouvoir simuler correctement les concentrations en nombre, dans les ambiances confinées comme dans l'atmosphère. Un modèle de dynamique des particules capable de suivre avec autant de précision la concentration en nombre que la concentration en masse, avec un temps de calcul optimal, a été développé. La dynamique des particules dépend de divers processus, les plus importants étant la condensation/évaporation, suivie par la nucléation, la coagulation, et les phénomènes de dépôts. Ces processus sont bien connus pour les particules fines et grossières, mais dans le cas des nanoparticules, certains phénomènes additionnels doivent être pris en compte, notamment l'effet Kelvin pour la condensation/ évaporation et les forces de van der Waals pour la coagulation. Le travail a tout d'abord porté sur le processus de condensation/évaporation, qui s'avère être le plus compliqué numériquement. Les particules sont présumées sphériques. L'effet Kelvin est pris en compte car il devient considérable pour les particules de diamètre inférieur à 50 nm. Les schémas numériques utilisés reposent sur une approche sectionnelle : l'échelle granulométrique des particules est discrétisée en sections, caractérisées par un diamètre représentatif. Un algorithme de répartition des particules est utilisé, après condensation/évaporation, afin de conserver les diamètres représentatifs à l'intérieur de leurs sections respectives. Cette redistribution peut se faire en terme de masse ou de nombre. Un des points clé de l'algorithme est de savoir quelle quantité, de la masse ou du nombre, doit être redistribuée. Une approche hybride consistant à répartir la quantité dominante dans la section de taille considérée (le nombre pour les nanoparticules et la masse pour les particules fines et grossières) a été mise en place et a permis d'obtenir une amélioration de la précision du modèle par rapport aux algorithmes existants, pour un large choix de conditions. Le processus de coagulation pour les nanoparticules a aussi été résolu avec une approche sectionnelle. La coagulation est régie par le mouvement brownien des nanoparticules. Pour cette approche, il a été constaté qu'il est plus efficace de calculer le noyau de coagulation en utilisant le diamètre représentatif de la section plutôt que de l'intégrer sur la section entière. Les simulations ont aussi pu montrer que les interactions de van der Waals amplifient fortement le taux de coagulation pour les nanoparticules. La nucléation a été intégrée au modèle nouvellement développé en incorporant un terme source de nanoparticules dans la première section, commençant à un nanomètre. La formulation de ce taux de nucléation correspond à celle de l'acide sulfurique mais le traitement des interactions numériques entre nucléation, coagulation et condensation/évaporation est générique. Différentes stratégies de couplage visant à résoudre séparément ou en même temps les trois processus sont discutées. Afin de pouvoir proposer des recommandations, différentes méthodes numériques de couplage ont été développées puis évaluées par rapport au temps de calcul et à la précision obtenue en terme de concentration massique et numérique / It is necessary to adapt existing models in order to simulate the number concentration, and correctly account for nanoparticles, in both free and confined atmospheres. A model of particle dynamics capable of following accurately the number as well as the mass concentration of particles, with an optimal calculation time, has been developed. The dynamics of particles depends on various processes, the most important ones being condensation/evaporation, followed by nucleation, coagulation, and deposition phenomena. These processes are well-known for fine and coarse particles, but some additional phenomena must be taken into account when applied to nanoparticles, such as the Kelvin effect for condensation/evaporation and the van der Waals forces for coagulation. This work focused first on condensation/evaporation, which is the most numerically challenging process. Particles were assumed to be of spherical shape. The Kelvin effect has been taken into account as it becomes significant for particles with diameter below 50 nm. The numerical schemes are based on a sectional approach : the particle size range is discretized in sections characterized by a representative diameter. A redistribution algorithm is used, after condensation/ evaporation occurred, in order to keep the representative diameter between the boundaries of the section. The redistribution can be conducted in terms of mass or number. The key point in such algorithms is to choose which quantity has to be redistributed over the fixed sections. We have developed a hybrid algorithm that redistributes the relevant quantity for each section. This new approach has been tested and shows significant improvements with respect to most existing models over a wide range of conditions. The process of coagulation for nanoparticles has also been solved with a sectional approach. Coagulation is monitored by the brownian motion of nanoparticles. This approach is shown to be more efficient if the coagulation rate is evaluated using the representative diameter of the section, rather than being integrated over the whole section. Simulations also reveal that the van derWaals interactions greatly enhance coagulation of nanoparticles. Nucleation has been incorporated into the newly developed model through a direct source of nanoparticles in the first size section, beginning at one nanometer. The formulation of this rate of nucleation corresponds to that of sulfuric acid but the treatment of the numerical interactions between nucleation, coagulation and condensation/evaporation is generic. Various strategies aiming to solve separately or jointly these three processes are discussed. In order to provide recommendations, several numerical splitting methods have been implemented and evaluated regarding their CPU times and their accuracy in terms of number and mass concentrations

La thèse porte sur l'étude de la réduction de traînée par injection de bulles. La réduction de traînée présente un intérêt pour les applications navales puisqu'elle est issue d'une modification des structures cohérentes qui contribuent le plus à la résistance à l'avancement. Le but de cette étude est d'analyser les mécanismes à l'origine de cette diminution du frottement pariétal. L'approche utilisée dans le cadre de cette étude est numérique, elle emploie le code JADIM par une approche Euler-Lagrange : la phase continue est simulée par Simulation Numérique Directe et la phase dispersée est simulée en suivant individuellement chaque bulle. La configuration retenue dans le cadre de cette étude est celle de l'écoulement de Taylor-Couette (écoulement compris entre deux cylindres en rotation). La première partie de la thèse vise à adapter l'outil numérique employé, afin de pouvoir prendre en compte le retour de la phase dispersée via des termes de forçage dans les équations bilan de matière et de quantité de mouvement. La deuxième partie de la thèse vise à étudier l'écoulement porteur en configuration monophasique, afin de disposer d'une référence sur l'écoulement non perturbé. La troisième partie de la thèse a pour objectif d'étudier la dispersion passive des bulles dans le système, afin d'analyser les mécanismes de migrations mis en jeu. Enfin la dernière partie de la thèse vise à étudier l'influence des bulles sur l'écoulement porteur en analysant l'effet de certains paramètres, notamment le taux de vide et la flottabilité. / The study deals with drag reduction induced by bubble injection, its application concerns naval transport. The aim of the study is to shed more light on mechanisms that are involved in this wall friction reduction. The study is based on a numerical approach, and use the JADIM code with an Euler-Lagrange approach: the continuous phase is solved by Direct Numerical Simulation, and the disperse phase by a tracking of each bubble. Within the framework of this study we consider the Taylor-Couette flow configuration (flow between two concentric cylinders in rotation). The first part of the study deals with the modification of the numerical tool, in order to take into account the influence of the disperse phase on the continuous one through forcing terms in the mass and momentum balance equations. The aim of the second part is to study de Taylor-Couette flow in its monophasic configuration, for the purpose of providing a reference of the undisturbed flow. The third part deals with the passive dispersion of bubble in Taylor-Couette flow, in order to analyze the migration mechanisms involved. And the aim of the last part is to study the effects of the disperse phase on the continuous one, by analyzing the influence of bubbly phase parameters (like void fraction and buoyancy).

Ce travail de thèse, menée dans le cadre d'une convention CIFRE entre TURBOMECA et le CERFACS, s'inscrit dans un contexte d'amélioration des performances des turbines de type axial équipant les turboréacteurs d'hélicoptère. L'une des principales difficultés rencontrée dans cette démarche concerne la maîtrise de la température que voient les pales de ce composant, notamment la roue haute pression. Les travaux de cette thèse s'articulent autour de deux axes principaux: - Le premier traite l'analyse de la Simulations aux Grandes Echelles (SGE) autour de pales. Une approche numérique SGE sur des maillages non-structurés est comparée aux résultats Reynolds Averaged Navier-Stokes (RANS) sur des maillages structurés, usuels dans ce type de configuration, ainsi qu'à une approche SGE sur maillages structurés. La SGE sur maillage non-structuré démontre sa capacité à prendre en compte les phénomènes qui ont un impact sur les flux de chaleur pariétaux. - Le second axe de recherche a pour objectif de développer un outil numérique de couplage pour assurer le transfert d'information entre un code SGE réactif sur maillage non-structuré, employé dans les chambres de combustion, et un code non-réactif en RANS, utilisé par les industriels pour modéliser l'étage turbine. Cet outil a été validé sur plusieurs cas tests qui montrent le potentiel de cette méthodologie pour le couplage multi-composant. / This PhD dissertation, conducted as part of a CIFRE research project between TURBOMECA and CERFACS, deals with improving performance of axial turbines from helicopter engines. One of the main difficulties with such an objective is the control of the temperature prediction around the blades, especially the temperature of the high pressure rotor. The work of this thesis focusses on two axes: - First concerns the analysis of Large Eddy Simulation (LES) predictions around blades: a numerical LES approach on unstructured meshes is compared to Reynolds Averaged Navier-Stokes (RANS) results on structured meshes as well as to LES on structured meshes. LES on unstructured meshes demonstrates its capacity of taking into account the phenomena which have an impact on wall heat flux around blades. - The second axis deals with the development of a numerical tool for coupling and transferring information between a reactive LES code, used in combustion chambers, and a non-reactive RANS solver, employed by industrial actors for modeling the turbine stage. This tool is validated on a number of test cases which show the potential of this methodology for multi-component predictions.

L'objet du présent travail est la modélisation, par simulation numérique, de la tenue mécanique de flasques de sièges automobile obtenus par semi-découpe de tôle d'acier à haute limite élastique. Des essais de compression, de traction et de cisaillement sont mis en place pour étudier l'anisotropie plastique, l'écrouissage, l'endommagement et la rupture. Des essais non proportionnels de cisaillement suivi de traction et de semi-découpe suivie de poinçonnement permettent d'investiguer l'effet de prédéformations, à bas taux de triaxialité des contraintes, sur le comportement du matériau à nouveau sollicité à un taux de triaxialité des contraintes plus élevé. Basé sur les travaux de Xue et Wierzbicki, un critère de rupture permettant de rendre compte des observations expérimentales est proposé et comparé à un modèle d'endommagement de Lemaitre. Le critère de rupture proposé est implanté dans deux logiciels de calcul par éléments finis : le logiciel Forge® dédié à la simulation de la mise en forme et le logiciel LS-Dyna® dédié à la simulation de crash. Dans un premier temps, la mise en forme par semi-découpe est simulée avec Forge®. Dans un second temps, les champs de variables d'état sont transportés du maillage résultat obtenu avec Forge®, vers le maillage d'entrée pour LS-Dyna®. Ensuite les calculs de tenue mécanique sont réalisés avec LS-Dyna®. Les résultats de calculs de tenue mécanique sont enfin comparés à des résultats expérimentaux. Ces résultats expérimentaux sont des essais de mise en forme par semi-découpe suivis d'essais d'arrachement qui reproduisent, en laboratoire, les sollicitations observées sur les pièces industrielles. / The purpose of this work is to model by digital simulation the mechanical strength of automotive seat flanges obtained by half-blanking of high yield strength steel sheets. Compression, tensile and shear tests are done in order to investigate plastic anisotropy, strain hardening, ductile damage and fracture. Non-proportional tests of shear followed by traction and of half-blanking followed by punching allow investigating the effect of predeformations at low stress triaxiality, on the behavior of the material loaded again at higher stress triaxiality. Based on the work of Xue and Wierzbicki, a fracture criterion that allows taking into account the experimental observations is proposed and compared to a Lemaitre damage model. The proposed fracture criterion is implemented in two finite element software programms: Forge®, dedicated to forming simulation and LS-Dyna®, dedicated to crash simulation. In a first step, forming by half-blanking is simulated with Forge®. In a second step, state variable fields are mapped from outcome mesh from Forge® to input mesh for LS-Dyna®. Then mechanical strength computations are performed using LS-Dyna®. Finally the mechanical strength computation results are compared to experimental results. These experimental results are forming tests by half-blanking followed by wrenching tests. These tests reproduce in the laboratory loading observed on industrial parts.

Cette thèse a pour objet le développement d'un modèle numérique de simulation des fluides fondé sur la méthode Smoothed Particle Hydrodynamics (SPH). SPH est une méthode de simulation numérique sans maillage présentant un certain nombre d'avantages par rapport aux méthodes Eulériennes. Elle permet notamment de modéliser des écoulements à surface libre ou interfaces fortement déformées. Ce travail s'adresse principalement à quatre problématiques liées aux fondements de la méthode SPH : l'imposition des conditions aux limites, la prédiction précise des champs de pression, l'implémentation d'un modèle thermique et la réduction des temps de calcul. L'objectif est de modéliser des écoulements industriels complexes par la méthode SPH, en complément de ce qui peut se faire avec des méthodes à maillage. Typiquement, les problèmes visés sont des écoulements 3-D à surface libre ou confinés, pouvant interagir avec des structures mobiles et/ou transporter des scalaires, notamment des scalaires actifs (e.g. température). Dans ce but, on propose ici un modèle SPH incompressible (ISPH) basé sur une représentation semi-analytique des conditions aux limites. La technique des conditions aux limites semi-analytiques permet d'imposer des conditions sur la pression de manière précise et physique, contrairement à ce qui se fait avec des conditions aux limites classiques en SPH. Un modèle k-epsilon a été incorporé à ce nouveau modèle ISPH, à partir des travaux de Ferrand et al. (2013). Un modèle de flottabilité a également été ajouté, reposant sur l'approximation de Boussinesq. Les interactions entre flottabilité et turbulence sont prises en compte. Enfin, une formulation pour les frontières ouvertes dans le nouveau modèle est établie. La validation du modèle en 2-D a été réalisée sur un ensemble de cas-tests permettant d'estimer les capacités de prédiction du nouveau modèle en ce qui concerne les écoulements isothermes et non-isothermes, laminaires ou turbulents. Des cas confinés sont présentés, ainsi que des écoulements à surface libre (l'un d'eux incluant un corps solide mobile dans l'écoulement). La formulation pour les frontières ouvertes a été testée sur un canal de Poiseuille plan laminaire et sur deux cas de propagation d'une onde solitaire. Des comparaisons sont présentées avec des méthodes à maillage, ainsi qu'avec un modèle SPH quasi-incompressible (WCSPH) avec le même type de conditions aux limites. Les résultats montrent que le modèle permet de représenter des écoulements dans des domaines à géométrie complexe, tout en améliorant la prédiction des champs de pression par rapport à la méthode WCSPH. L'extension du modèle en trois dimensions a été réalisée dans un code massivement parallèle fonctionnant sur carte graphique (GPU). Deux cas de validation en 3-D sont proposés, ainsi que des résultats sur un cas simple d'application en 3-D / In this work a numerical model for fluid flow simulation was developed, based on the Smoothed Particle Hydrodynamics (SPH) method. SPH is a meshless Lagrangian Computational Fluid Dynamics (CFD) method that offers some advantages compared to mesh-based Eulerian methods. In particular, it is able to model flows presenting highly distorted free-surfaces or interfaces. This work tackles four issues concerning the SPH method : the imposition of boundary conditions, the accuracy of the pressure prediction, the modelling of buoyancy effects and the reduction of computational time. The aim is to model complex industrial flows with the SPH method, as a complement of what can be done with mesh-based methods. Typically, the targetted problems are 3-D free-surface or confined flows that may interact with moving solids and/or transport scalars, in particular active scalars (e.g. the temperature). To achieve this goal, a new incompressible SPH (ISPH) model is proposed, based on semi-analytical boundary conditions. This technique for the representation of boundary conditions in SPH makes it possible to accurately prescribe consistent pressure boundary conditions, contrary to what is done with classical boundary conditions in SPH. A k-epsilon turbulence closure is included in the new ISPH model. A buoyancy model was also added, based on the Boussinesq approximation. The interactions between buoyancy and turbulence are modelled. Finally, a formulation for open boundary conditions is proposed in this framework. The 2-D validation was performed on a set of test-cases that made it possible to assess the prediction capabilities of the new model regarding isothermal and non-isothermal flows, in laminar or turbulent regime. Confined cases are presented, as well as free-surface flows (one of them including a moving body in the flow). The open boundary formulation was tested on a laminar plane Poiseuille flow and on two cases of propagation of a solitary wave. Comparisons with mesh-based methods are provided with, as well as comparisons with a weakly-compressible SPH (WCSPH) model using the same kind of boundary conditions. The results show that the model is able to represent flows in complex boundary geometries, while improving the pressure prediction compared to the WCSPH method. The extension of the model to 3-D was done in a massively parallel code running on a Graphic Processing Unit (GPU). Two validation cases in 3-D are presented, as well as preliminary results on a simple 3-D application case

Avec les dernières augmentations de l'intensité maximale de laser réalisable grâce à de courtes impulsions à haute puissance (gamme femtoseconde) un intérêt a surgi dans les sources de plasma laser potentiels. Les lasers sont utilisés en radiographie proton, allumage rapide, hadronthérapie, la production de radioisotopes et de laboratoire astrophysique. Au cours de l'interaction laser-cible, les ions sont accélérés par des processus physiques différents, en fonction de la zone de la cible. Tous ces mécanismes ont un point commun: les ions sont accélérés par des champs électriques intenses, qui se produisent en raison de la séparation de forte charge induite par l'interaction de l'impulsion laser avec la cible, directement ou indirectement. Deux principales sources distinctes pour le déplacement de charge peuvent être mis en évidence. Le premier est le gradient de charge provoquée par l'action directe de la force ponderomotive de laser sur les électrons dans la surface avant de la cible, qui est la prémisse pour le processus d'accélération des radiations de pression (RPA). Une deuxième source peut être identifiée comme provenant du rayonnement laser qui est transformée en énergie cinétique d'une population d'électrons relativistes chaud (~ quelques MeV). Les électrons chauds se déplacent et font recirculer à travers la cible et forment un nuage d'électrons relativistes à la sortie de la cible dans le vide. Ce nuage, qui se prolonge pour plusieurs longueurs de Debye, crée un champ électrique extrêmement intense longitudinal, la plupart du temps dirigé le long de la surface normale, ce qui, par conséquent, est la cause de l'accélération d'ions efficace, qui conduit à l'accélération cible normale gaine (TNSA) processus . Le mécanisme TNSA permet d'utiliser des géométries différentes cibles afin de parvenir à une meilleure focalisation des faisceaux de particules de l'ordre de plusieurs dizaines de microns, avec des densités d'énergie élevées. Les électrons chauds sont produits par l'irradiation d'une feuille solide avec une impulsion laser intense; ces électrons sont transportés à travers la cible, la formation d'un champ électrostatique fort, normal à la surface cible. Protons et les ions chargés positivement de la surface arrière de la cible sont accélérés par ce domaine jusqu'à ce que la charge de l'électron est compensée. La densité d'électrons chauds et la température dans le vide arrière dépendent des propriétés géométriques et de composition cibles tels que la courbure de la cible, les structures de mise au point d'impulsion et de microstructure pour l'accélération de protons améliorée. Au cours de ma première année, j'ai étudié les effets de la géométrie de la cible sur le proton et l'ion énergie et la distribution angulaire afin d'optimiser les faisceaux de particules laser accéléré au moyen de deux dimensions (2D) particule-in-cell (PIC) simulations de l'interaction de l'ultra-court impulsions laser avec plusieurs cibles microstructurées. Également au cours de cette année, je l'ai étudié la théorie derrière les modèles utilisés. / With the latest increases in maximum laser intensity achievable through short pulses at high power (femtosecond range) an interest has arisen in potential laser plasma sources. Lasers are used in proton radiography, rapid ignition, hadrontherapy, production of radioisotopes and astrophysical laboratory. During the laser-target interaction, the ions are accelerated by different physical processes, depending on the area of ​​the target. All these mechanisms have one thing in common: the ions are accelerated by intense electric fields, which occur due to the separation of high charge induced by the interaction of the laser pulse with the target, directly or indirectly. Two main distinct sources for charge displacement can be identified. The first is the charge gradient caused by the direct action of the laser ponderomotive force on the electrons in the front surface of the target, which is the premise for the pressure ramping acceleration (RPA) process. A second source can be identified as coming from the laser radiation which is transformed into kinetic energy of a hot relativistic electron population (~ a few MeV). The hot electrons move and recirculate through the target and form a cloud of relativistic electrons at the exit of the target in a vacuum. This cloud, which extends for several lengths of Debye, creates an extremely intense longitudinal electric field, mostly directed along the normal surface, which is therefore the cause of effective ion acceleration, which leads to the normal target sheath acceleration (TNSA) process. The TNSA mechanism makes it possible to use different target geometries in order to obtain a better focusing of the beams of particles on the order of several tens of microns, with high energy densities. Hot electrons are produced by irradiating a solid sheet with an intense laser pulse; these electrons are transported through the target, forming a strong electrostatic field, normal to the target surface. Protons and positively charged ions from the back surface of the target are accelerated by this domain until the charge of the electron is compensated. The density of hot electrons and the temperature in the back vacuum depend on the target geometric and compositional properties such as target curvature, pulse and microstructure tuning structures for enhanced proton acceleration. In my first year I studied the effects of target geometry on the proton and energy ion and angular distribution in order to optimize the accelerated laser particle beams by means of two-dimensional (2D) particle -in-cell (PIC) simulations of the interaction of ultra-short laser pulses with several microstructured targets. Also during this year, I studied the theory behind the models used.