Cette thèse met l'accent sur la recherche scientifique pour développer une classe de procédés à hautesvitesses de déformation des tôles métalliques en alliages d'aluminium et en acier à haute résistance (AHR).Ces technologies emploient une impulsion de pression de courte durée qui propulse la tôle dans une matrice.Ces procédés sont généralement décrits comme procédés de formage par impulsion de pression (PPF). Letravail proposé dans ce mémoire de thèse a permis de surmonter trois obstacles pour l'utilisation desprocédés PPF et la fabrication à moindre coût de structures légères. Le premier obstacle a été le manque decorrélation entre formabilité et vitesses de déformation qui se développent lors d'un procédé PPF. Nous avonsproposé d'analyser la formabilité et la rupture des tôles, et de caractériser les vitesses de déformation et leurl'hétérogénéité pendant le procédé PPF. Le deuxième obstacle a été le manque de lois constitutives validéespour les métaux déformés par le procédé PPF. Nous avons étudié la microstructure et l'évolution despropriétés mécaniques durant le procédé PPF. Le troisième obstacle est le manque de modèles de formabilitéprédictifs validés pour le procédé PPF. Nous avons utilisé la méthode Marciniak-Kuczynski pour la prédictionde la formabilité de l'alliage AA5182 et de l'acier DP600 sous un large éventail de vitesses de déformation etsous différentes directions de ces vitesses. La combinaison de ces résultats de recherche permet une plusgrande capacité prédictive pour concevoir et développer des procédés PPF pour composants d'automobiledésirés à partir d'aluminium et d'acier AHR.

Cette thèse traite de l'étude de phénomènes ultra-rapides, et plus particulièrement de dynamiques de relaxation initiées par une impulsion femtoseconde (10-15 s). Dans la première partie, l'imagerie de vecteurs vitesses a permis d'étudier deux systèmes : l'iodure de méthyle (CH3I), excité dans son premier état de Rydberg et le tetrathiafulvalène (TTF - C6H4S4). Dans le TTF, le but est de caractériser les deux bandes majoritaires d'absorption vers 300 nm. Durant l'étude de la prédissociation de CH3I, la distribution vibrationnelle du fragment CH3 a été caractérisée, alors que pour le TTF une force de liaison du dimère a été extraite. L'interprétation des spectres de photoélectrons de CH3I permet une bonne compréhension des dynamiques depuis et/ou en résonance avec des états proches du potentiel d'ionisation. Dans la seconde partie, la structure électronique d'atomes ou de molécules est sondée par la génération d'harmoniques d'ordre élevé (GHOE). Le minimum de la section efficace totale de photoionisation de l'argon (minimum de Cooper) a été observé et étudié dans le spectre harmonique. Cette étude a permis, entre autre, la mise au point d'un modèle théorique fiable et complet. Dans une autre étude, la GHOE a montré la sensibilité du rayonnement harmonique à la chiralité lors de interaction des 2 énantiomères de la fenchone (C10H16O), avec un rayonnement polarisé elliptiquement. Pour finir, l'utilisation de la GHOE pour sonder une dynamique moléculaire, amorcée par un réseau transitoire d'excitation à 400 nm, a été réalisée sur le dioxyde d'azote (NO2). Ceci à l'échelle picoseconde : photodissociation, et femtoseconde : transfert de population par le biais d'une intersection conique.

Cette thèse présente la construction de schémas distribuant le résidu (RD) d'ordre très élevés, pour la discrétisation d'équations d'advection-diffusion multidimensionnelles et stationnaires sur maillages non structurés. Des schémas linéaires ainsi que des schémas non linéaires sont considérés. Une approximation de la solution polynomiale par morceaux et continue sur chaque élément est adoptée, de plus une procédure de reconstruction du gradient que celle de la solution numérique est utilisée afin d'avoir une représentation continue de la solution numérique et de son gradient. Il est montré que le gradient doit être reconstruit avec la même précision de la solution, sans quoi la précision formel du schéma numérique est perdue dans les cas où les effets de diffusion prévalent sur les effets d'advection, et aussi quand l'advection et la diffusion sont également importants. Ensuite, la méthode est étendue à des systèmes d'équations, en particulier aux équations de Navier-Stokes et aux équations RANS. La précision, l'efficacité et la robustesse du solveur RD implicite sont démontrées sur plusieurs cas tests.

No description available.

Cette thèse est principalement dédiée à la simulation des écoulements massivement décollés dans le domaine spatial. Nous avons restreint notre étude aux écoulements d'arrière-corps, pour lesquels ces décollements sont imposés par des changements brutaux de la géométrie. Dans le domaine spatial, le caractère fortement compressible des écoulements rencontrés impose l'utilisation de schémas numériques robustes. D'un autre coté, la simulation fine de la turbulence impose des schémas d'ordre élevé et peu dissipatifs. Ces deux spécifications, apparemment contradictoires, doivent pourtant coexister au sein d'une même simulation. Les modèles de turbulence ainsi que les schémas de discrétisation sont indissociables et leur couplage doit impérativement être considéré. Les schémas numériques doivent garder leur précision formelle dans des géométries complexes et des maillages très irréguliers imposés par le contexte industriel. Cette étude analyse le schéma de discrétisation utilisé dans le code de calcul FLUSEPA développé par Airbus Defence & Space. Ce schéma est robuste et précis pour des écoulements avec chocs et il présente une faible sensibilité au maillage (l'ordre 3 étant conservé même sur des maillages fortement perturbés). Malheureusement, le schéma possède une trop faible résolvabilité liée à un niveau de dissipation trop élevé pour envisager des simulations hybrides RANS/LES. Pour pallier à cet inconvénient, nous nous sommes penchés vers une solution basée sur un recentrage conditionnel et local : dans les zones dominées par des structures tourbillonnaires, une fonction analytique assure un recentrage local lorsque la stabilité numérique le permet. Cette condition de stabilité assure le couplage entre le schéma et le modèle. De cette manière, les viscosités laminaire et tourbillonnaire sont les seules à jouer un rôle dans les régions dominées par la vorticité et servent aussi à stabiliser le schéma numérique. Cette étude présente de plus une comparaison qualitative et quantitative de plusieurs modèles hybrides RANS/LES, à égalité de maillage et de schéma utilisés Pour cela, un certain nombre d'améliorations (notamment de leur capacité à résoudre les instabilités de Kelvin-Helmohlotz sans retard), proposées dans la littérature ou bien introduites dans cette thèse, sont prises en compte. Les applications numériques étudiées concernent des géométries allant de la marche descendante au lanceur spatial complet à échelle réduite. / This thesis is mainly dedicated to the simulation of massively separated flows in the space domain. We restricted our study to afterbody flows, where the separation is imposed by abrupt geometry changes. In the space domain, highly compressible flows require the use of robust numerical schemes. On the other hand, the simulation of turbulence imposes high-order low dissipative numerical schemes. These two specifications, apparently contradictory, must coexist within the same simulation. The coupling between turbulence models and discretization schemes is of the utmost importance and must be considered. Numerical schemes should keep their formal accuracy on complex geometries and on very irregular meshes imposed by the industrial context. In this research, we analyze the discretization scheme implemented in the FLUSEPA solver, developed by Airbus Defence & Space. Such a scheme is robust and accurate for flows with shocks and exhibits a low sensitivity to the grid (the third order of accuracy being ensured, even on highly irregular grids). Unfortunately, the scheme possesses a too low resolvability related to a too high numerical dissipation for RANS/LES simulations. To circumvent this problem, we considered a conditional and local re-centering strategy: in regions dominated by vortical structures, an analytic function provides local re-centering when a numerical stability condition is satisfied. This stability condition ensures the coupling between the numerical scheme and the model. In this way, only the turbulent and the laminar viscosities play a role in regions dominated by vorticity, and also allow to stabilize the numerical scheme. This study provides also a qualitative and quantitative assessment of several hybrid RANS/LES models, using the same grids and discretization scheme. For this purpose some recent improvements (improving their ability to trigger the Kelvin-Helmohlotz instabilities without delay), proposed in the litterature or suggested in this work, are taken into account. Numerical applications include geometrical configurations ranging from a backward facing step to realistic launcher configurations.

Les simulations numériques réalistes en électrophysiologie cardiaque ont un coût de calcul extrêmement élevé. Ce coût s’explique en grande partie par la raideur, à la fois en temps et en espace, d’une onde de « potentiel d’action » (PA). Par ailleurs, les phénomènes observés sont très instationnaires et s’étudient en temps long. Une description précise de la dynamique des PA est cruciale pour construire des modèles numériques pertinents d’un point de vue médical ou clinique. Cet aspect fondamental ne peut être contourné dans les études numériques réalistes.La raideur de l’onde de PA ne peut être captée numériquement qu’en ayant recours à des maillages très fins. Ces maillages très fins induisent un coût de calcul très important, et introduisent aussi des erreurs supplémentaires : les systèmes linéaires à résoudre deviennent très mal conditionnés. Au final, les erreurs numériques peuvent être particulièrement grandes dans les simulations alors que leur contrôle est évidemment essentiel pour assurer la fiabilité des résultats. Jusqu’à présent, très peu de résultats sont disponibles pour assurer cette fiabilité. Dans les faits, les erreurs sont la plupart du temps contrôlées par des procédés empiriques. Il existe quelques résultats théoriques étudiant la convergence et la stabilité des schémas numériques associés. En pratique, en plus d'avoir un contrôle de l'erreur sur le potentiel, il est aussi nécessaire d'avoir un contrôle de l’erreur sur des quantités macroscopiques décrivant la dynamique de l’onde de PA : temps d’activation, durée du PA, propriétés de restitution... Ces quantités ont en effet une interprétation physiologique qui permet de caractériser le caractère arythmogène des tissus.Les modèles sont des systèmes d’EDP de réaction-diffusion couplés avec des systèmes d’équations différentielles pouvant être très raides, les modèles ioniques. Ils sont actuellement discrétisés par éléments finis conforme (Lagrange) et par des schémas en temps d’ordre un ou deux. Dans ce travail, nous concevons et évaluons l’intérêt d'utiliser des méthodes d’ordre supérieure pour ces systèmes. Parallèlement nous introduisons d'une part une nouvelle classe de schémas appelé schémas exponentiel Adams Bashforth intégral (IEAB), et d'autre part des schémas Rush Larsen (RL) d'ordre élevé. Ces nouveaux schémas sont des schémas multipas de type exponentiels. Nous montrons qu'ils possèdent des bonnes propriétés de stabilité et permettent de faire face efficacement à la raideur des modèles ioniques. Les schémas que nous proposons sont comparés numériquement (en terme de précision, coût en temps de calcul et stabilité) à plusieurs schémas classiques, ainsi qu'aux schémas exponentiels (RL1, RL2) communément utilisés pour des simulations en électrophysiologie cardiaque. Nous proposons des techniques permettant de calculer avec précision les quantités d’intérêts cliniques (temps d’activation, de récupération, durée du potentiel d’action). Des résultats théoriques de convergence en temps et de convergence globale (espace et temps) sont énoncés et prouvés. Ces résultats sont ensuite illustrés numériquement à travers le modèle monodomaine et les modèles ioniques de Beeler Reuter, de Ten Tusscher et al. L’intérêt d'utiliser des schémas d'ordre élevés est aussi évalué sur des ondes spirales en 2D et 3D. / Realistic numerical simulations in cardiac electrophysiology have a computational cost of extremely high. This cost is largely explained by the stiffness both in time and space, of the action potential (AP) wave. Moreover, the observed phenomena are very unsteady and are studied in long time. A precise description of the dynamic of AP is crucial for constructing relevant numerical models, from a medical or clinical perspective. This fundamental aspect can not be circumvented in realistic numerical studies.The stiffness of AP wave can only be captured numerically, by using very fine meshes. In addition to the high computational cost, these very fine meshes also introduce additional errors : the linear systems to solve become very badly conditioned. In the end, the numerical errors can be particularly large whereas their control is obviously essential to ensure the reliability of the results. So far very few results are available to ensure this reliability. In practice, the errors are mostly controlled by empirical processes. In practice, in addition of having a control of the error on the potential, it is also necessary to have an error control on macroscopic quantities describing the dynamics of the AP wave : activation time, AP duration, properties of restitution ... These quantities have indeed a physiological interpretation which allows to characterize the arrhythmogenic character of the tissues.The models are systems of reaction diffusion PDE coupled with systems of differential equations that can be very stiffs (ionic models). They are currently discretized by conforming finite elements (Lagrange finite elements methods) and by schemes in time of order one or two. In this work, we design and evaluate the interest of using higher order methods for these systems. At the same time, we introduce on the one hand, a new class of schemes called Integral Exponential Adams Bashforth (IEAB) schemes and, on the other hand, high order Rush Larsen (RL) schemes. These new schemes are exponential time-stepping schemes. We show that they have good stability properties and can efficiently cope with the stiffness of ionic models. The schemes we propose are numerically compared (in terms of accuracy, CPU time and stability) with several classical schemes, as well as with the exponential schemes (RL1, RL2), commonly used for cardiac electrophysiology simulations. We propose good techniques for accurately calculating quantities of clinical interest (activation time, recovery time, duration of action potential). Theoretical results of convergence in time and global convergence (in space and time) are stated and proved. These results are then illustrated numerically through the monodomain model and the ionic models of Beeler Reuter, Ten Tusscher et al. The advantage of using high order schemes is also evaluated on spiral waves in 2D and 3D.

These thesis review some recent results on the construction of very high order multidimensional upwind schemes for the solution of steady and unsteady conservation laws on unstructured triangular grids.<p>We also consider the extension to the approximation of solutions to conservation laws containing second order dissipative terms. To build this high order schemes we use a subtriangulation of the triangular Pk elements where we apply the distribution used for a P1 element.<p>This manuscript is divided in two parts. The first part is dedicated to the design of the high order schemes for scalar equations and focus more on the theoretical design of the schemes. The second part deals with the extension to system of equations, in particular we will compare the performances of 2nd, 3rd and 4th order schemes.<p><p>The first part is subdivided in four chapters:<p>The aim of the second chapter is to present the multidimensional upwind residual distributive schemes and to explain what was the status of their development at the beginning of this work.<p>The third chapter is dedicated to the first contribution: the design of 3rd and 4th order quasi non-oscillatory schemes.<p>The fourth chapter is composed of two parts: we start by understanding the non-uniformity of the accuracy of the 2nd order schemes for advection-diffusion problem. To solve this issue we use a Finite Element hybridisation.<p>This deep study of the 2nd order scheme is used as a basis to design a 3rd order scheme for advection-diffusion.<p>Finally, in the fifth chapter we extend the high order quasi non-oscillatory schemes to unsteady problems.<p>In the second part, we extend the schemes of the first part to systems of equations as follows:<p>The sixth chapter deals with the extension to steady systems of hyperbolic equations. In particular, we discuss how to solve some issues such as boundary conditions and the discretisation of curved geometries.<p>Then, we look at the performance of 2nd and 3rd order schemes on viscous flow.<p>Finally, we test the space-time schemes on several test cases. In particular, we will test the monotonicity of the space-time non-oscillatory schemes and we apply residual distributive schemes to acoustic problems. / Doctorat en Sciences de l'ingénieur / info:eu-repo/semantics/nonPublished

Les protéines extraites des feuilles de luzerne sont une source importante de protéines. La filtration membranaire, technologie de séparation respectueuse de l’environnement avec une productivité élevée et de faible coût a été utilisée pour séparer et concentrer les protéines des feuilles de luzerne à partir de leur jus. Cependant le phénomène du colmatage de la membrane qui réduit sérieusement le flux et la séparation des protéines est un facteur limitant important dans l'application de la filtration membranaire. Pour améliorer la récupération des protéines et amenuiser le phénomène du colmatage, la filtration membranaire associée à fort cisaillement a été utilisée pour la filtration du jus de luzerne. Toujours dans l’objectif d'optimiser le processus de la filtration, "le mode de la filtration" et "les paramètres de fonctionnement" ont été étudiés pour réduire le colmatage de la membrane et améliorer le rendement de la filtration. Puis, l’effluent du jus de luzerne a été filtré par des membranes dans des conditions de fort cisaillement afin de recycler les protéines. En outre, le mécanisme du colmatage a été étudié et a permis d’évaluer les stratégies de contrôle du colmatage. L'optimisation du procédé membranaire, via l’étude du "mode de filtration" et des "paramètres de fonctionnement" a été conduit dans le but d’améliorer la séparation et la concentration des protéines et de réduire le colmatage. Trois types de « mode de filtration » ont été testés : la filtration frontale sur le module de la cellule amicon (DA), la filtration tangentielle dynamique sur le module à disque rotatif (CRDM) et la filtration frontale sur le module à disque rotatif (DRDM)). Les « paramètres de fonctionnement » qui ont été étudiés sont les suivants : le type de membranes (ultrafiltration (UF) et microfiltration (MF)), la vitesse de rotation, la température et la pression transmembranaire (TMP). Le comportement du débit (évolution du flux du perméat au cours de la filtration), les performances de la séparation (taux de clarification et de concentration), l’efficacité du nettoyage de la membrane (récupération de la perméabilité membranaire) et la productivité lors des tests de recyclage et de concentration ont été étudiés dans le but de définir des stratégies dans le contrôle du colmatage. Puis, l’effluent de luzerne a été filtré par UF afin de séparer et purifier les protéines. Le mécanisme du phénomène du colmatage des membranes lors de la filtration du jus de luzerne a été étudié. Le processus du colmatage de la membrane a montré une tendance d’un colmatage multi-site progressif. Le modèle du colmatage multisite progressif selon la loi de Darcy (SMDM) a été proposé afin de mieux décrire et comprendre le processus du colmatage. Les effets de la composition du fluide d’alimentation, du choix de la membrane et des conditions hydraulique ont joué un rôle important dans le processus progressif du colmatage. De plus, les coefficients de résistance et de compressibilité dans les différentes étapes et sites ont été calculés afin d’expliquer le processus complexe du colmatage et d’évaluer l'efficacité des stratégies du contrôle du colmatage. Une série d'essais avec de longues durées de filtration a été réalisée pour étudier le déclin du flux et le colmatage de la membrane à diverses étapes du processus. Ces résultats présentent une utilité pour améliorer la récupération des protéines et contrôler le colmatage dans le processus de la filtration membranaire à fort cisaillement du jus de luzerne. Ces résultats sont aussi utiles pour la conception et la mise en place des technologies membranaires dans les processus industriels. / Alfalfa leaf proteins extracted from plants are an important protein source. As an environmentally friendly separation technology with high productivity and low cost, membrane filtration was used to separate and concentrate leaf protein from alfalfa juice. However membrane fouling seriously reduces flux and protein separation and is an important limitation in the application of membrane filtration. To improve protein recovery and fouling control, dynamic shear-enhanced membrane filtration with high shear rate on membrane surface and excellent anti-fouling capacity was used for alfalfa juice filtration in this work. In order to optimize filtration process, filtration mode and operation parameters were investigated to reduce membrane fouling and improve separation performance. Then, alfalfa wastewater was also treated by dynamic shear-enhanced membrane filtration to recycle proteins. Furthermore, the fouling mechanism was studied and served as a valuable evaluation for fouling strategies. In this study, process optimization including “Filtration mode” and “Operation parameters” was studied to improve protein recovery and fouling control. In “Filtration mode”, three types of filtration modules (dead end filtration using laboratory Amicon cell (DA), dynamic cross filtration using rotating disk module (CRDM) and dead end filtration using rotating disk module (DRDM)) were used to investigate the filtration performance. As for “Operation parameters”, the operation parameters including membranes (ultrafiltration (UF) and microfiltration (MF)), rotating speed, temperature and transmembrane pressure (TMP) were studied to optimize the filtration process. Flux behavior (permeate flux and flux decline), separation performance (clarification and concentration capacity), membrane cleaning efficiency (permeability recovery) and productivity in full recycling tests and concentration tests were utilized to evaluate the various operation strategies. In addition, alfalfa wastewater was treated by UF membrane, while waste proteins were recycled. Fouling mechanism for alfalfa juice filtration was investigated. The fouling process showed significantly stepwise multisite patterns. Based on Darcy’s law, the stepwise multisite Darcy’s law model (SMDM) was proposed to better describe and understand the fouling process. The effects of feed composition, membrane and hydraulic conditions played an important role in stepwise fouling process. Moreover, the resistance coefficient and compressibility for different steps and sites were calculated to explain the complex fouling process and estimate the efficiency of flux decline control strategies. Besides, a series of long tests were utilized to study flux decline and membrane fouling at various fouling step process. These results can be used to understand the protein recovery and fouling control during shear-enhanced membrane filtration process of alfalfa juice. They have important implications for process design of membrane technology in industrial scale.

Cette thèse se concentre sur la contrôlabilité de l'indice de réfraction au niveau sub-micronique par changements d'indice induits par laser sur de longues dimensions i.e., avec des hauts rapports d'aspect élevés et des sections à l'échelle nanométrique. À cette fin, nous explorons les faisceaux ultracourts de Bessel non-diffractifs d'ordre zéro et les facteurs qui contribuent au confinement de l'énergie au-delà de la limite de diffraction. Le traitement par laser de matériaux transparents à l'aide de faisceaux non diffractifs offre un avantage important pour les structures sub-microniques en volume de haut rapport d'aspect à des fins applicatives en nanophotonique et en nanofluidique. Nous présentons l'effet de différentes conditions de focalisation et de paramètres laser sur la modification de la silice fondue, explorant ainsi les différents régimes d'interaction. Cette thèse aborde essentiellement des conditions modérées de focalisation car elles offrent un régime d'interaction stable sur une large gamme de paramètres laser, permettant l'ingénierie de la dispersion. La durée de l'impulsion laser s'est révélée être essentielle dans la définition du type de modification de l'indice de réfraction ou de modification structurale. Par exemple, l'usinage utilisant des impulsions laser femtosecondes entraîne une augmentation des structures d'indice de réfraction alors que les impulsions laser picosecondes engendrent une cavité uniforme i.e., des structures de faible indice. Pour acquérir un meilleur contrôle et une meilleure précision du dépôt d'énergie laser, un ensemble de mécanismes physiques responsables des dommages induits par laser dans des conditions d'excitation non-diffractives a été observé expérimentalement et examiné par des simulations indiquant le rôle essentiel de la diffusion de la lumière sur les électrons. Des mesures de microscopie pompe-sonde résolues en temps avec une résolution temporelle sub-picoseconde et spatiale sub-micronique donnent accès à l'excitation et à la relaxation dynamique instantanées. La transmission optique dynamique et le contraste de phase offrent des informations complémentaires sur la réponse électronique ou sur celle de la matrice vitreuse. La dynamique ultrarapide des porteurs libres a été particulièrement étudiée puisque le transfert d'énergie des électrons vers le réseau est la clé de transformation ultérieure du matériau. Le rôle de l'excitation instantanée pour différentes durées et énergie d'impulsion laser est exposé. Ainsi, la dynamique complète des porteurs de charge est présentée pour différents paramètres du laser. En particulier, la dynamique d'obtention de structures d'indice de réfraction positif et des cavités uniformes indique deux chemins différents de relaxation électronique et de dépôt de l'énergie: une relaxation rapide par l'intermédiaire de défauts pour les structures d'indice positif et une relaxation thermomécanique lente pour les cavités nanométriques. Enfin, en corrélant les résultats des études résolues en temps, les simulations et les résultats de photoluminescence après irradiation, nous formulons des scénarios potentiels de formation de l'indice de réfraction positif ainsi que des structures d'indice faible ou de vides uniformes / This thesis is focused on the controllability of laser-induced refractive index changes at sub-micron level over long dimensions i.e., with high aspect ratios and sections on the nanoscale. To this end, we explore non-diffractive zerothorder ultrafast Bessel beams and factors contributing to energy confinement beyond the diffraction limit. Laser processing of transparent materials using non-diffracting beams offers a strong advantage for high aspect ratio submicron structures inside the bulk in view of nanophotonics and nanouidics applications. We present the role of various focusing conditions and laser parameters on material modification in bulk fused silica and explore the different interaction regimes. This thesis tackles mostly the moderate focusing conditions as they offer a stable interaction regime backed up dispersion engineering over a large range of laser parameters. The laser pulse duration was found to be key in defining the type of laser induced refractive index or structural modification. For instance, machining using femtosecond laser pulses results in increased refractive index structures whereas picosecond laser pulses result in uniform void i.e., low index structures. To acquire better control over the laser energy deposition and precision, a range of physical mechanisms responsible for the laser induced damage in non-diffractive excitation conditions have been observed experimentally and further interrogated by simulations indicating a critical role of light scattering on carriers. Time-resolved pump-probe microscopy measurements with a sub-picosecond temporal and sub-micron spatial resolution allow access to the instantaneous excitation and relaxation dynamics. Dynamic optical transmission and phase contrast o_er complementary information of either electronic and glass matrix response. Primarily, ultrafast dynamics of free carriers was studied as the electron mediated energy transfer to the lattice is key to the subsequent material transformation. Role of instantaneous excitation at different laser pulse durations and energies is outlined. Then complete carrier dynamics is presented at different laser parameters. Particularly dynamics in conditions of positive refractive index structures and uniform voids is indicating two different paths of electronic relaxation and energy deposition: a fast defect mediated relaxation for positive index structures and slow thermomechanical relaxation for nanosize void structures. Finally, by correlating the results of time resolved studies, simulations and post-irradiated photoluminescence results, we formulate potential formation scenarios for the positive refractive index and low index or uniform void structures

De nos jours, la simulation des écoulements diphasiques a de plus en plus d’importance dans les chambres de combustion aéronautiques en tant qu’un des éléments requis pour analyser et maîtriser le processus complet de combustion, afin d’améliorer la performance du moteur et de mieux prédire les émissions polluantes. Dans les applications industrielles, la modélisation du combustible liquide trouvé en aval de l’injecteur sous forme de brouillard de gouttes polydisperse, appelé spray, est de préférence faite à l’aide de méthodes Eulériennes. Ce choix s’explique par les avantages qu’offrent ces méthodes par rapport aux méthodes Lagrangiennes, notamment la convergence statistique intrinsèque, le couplage aisé avec la phase gazeuse ainsi que l’efficacité pour le calcul haute performance. Dans la présente thèse, on utilise une approche Eulérienne basée sur une fermeture au niveau cinétique de type distribution Gaussienne Anisotrope (AG). L’AG résout des moments de vitesse jusqu’au deuxième ordre et permet de capter les croisements des trajectoires (PTC) à petite échelle de manière statistique. Le système d’équations obtenu est hyperbolique, le problème est bien-posé et satisfait les conditions de réalisabilité. L’AG est comparé au modèle monocinétique (MK) d’ordre 1 en vitesse. Il est approprié pour la description des particules faiblement inertielles. Il mène à un système faiblement hyperbolique qui peut générer des singularités. Plusieurs schémas numériques, utilisés pour résoudre les systèmes hyperboliques et faible- ment hyperboliques, sont évalués. Ces schémas sont classifiés selon leur capacité à traiter les singularités naturellement présentes dans les modèles Eulériens, sans perdre l’ordre global de la méthode ni rompre les conditions de réalisabilité. L’AG est testé sur un champ turbulent 3D chargé de particules dans des simulations numériques directes. Le code ASPHODELE est utilisé pour la phase gazeuse et l’AG est implémenté dans le code MUSES3D pour le spray. Les résultats sont comparés aux de simulations Lagrangiennes de référence et aux modèle MK. L’AG est validé pour des gouttes modérément inertielles à travers des résultats qualitatifs et quantitatifs. Il s’avère prometteur pour les applications complexes comprenant des PTC à petite échelle. Finalement, l’AG est étendu à la simulation aux grandes échelles nécessaire dans les cas réels turbulents dans le domaine industriel en se basant sur un filtrage au niveau cinétique. Cette stratégie aide à garantir les conditions de réalisabilités. Des résultats préliminaires sont évalués en 2D pour tester la sensibilité des résultats LES sur les paramètres des modèles de fermetures de sous mailles. / In aeronautical combustion chambers, the ability to simulate two-phase flows gains increasing importance nowadays since it is one of the elements needed for the full understanding and prediction of the combustion process. This matter is motivated by the objective of improving the engine performance and better predicting the pollutant emissions. On the industrial scale, the description of the fuel spray found downstream of the injector is preferably done through Eulerian methods. This is due to the intrinsic statistical convergence of these methods, their natural coupling to the gas phase and their efficiency in terms of High Performance Computing compared to Lagrangian methods. In this thesis, the use of Kinetic-Based Moment Method with an Anisotropic Gaussian (AG) closure is investigated. By solving all velocity moments up to second order, this model reproduces statistically the main features of small scale Particles Trajectories Crossing (PTC). The resulting hyperbolic system of equations is mathematically well-posed and satisfies the realizability properties. This model is compared to the first order model in the KBMM hierarchy, the monokinetic model MK which is suitable of low inertia particles. The latter leads to a weakly hyperbolic system that can generate δ-shocks. Several schemes are compared for the resolution of the hyperbolic and weakly hyperbolic system of equations. These methods are assessed based on their ability to handle the naturally en- countered singularities due to the moment closures, especially without globally degenerating to lower order or violating the realizability constraints. The AG is evaluated for the Direct Numerical Simulation of 3D turbulent particle-laden flows by using ASPHODELE solver for the gas phase, and MUSES3D solver for the Eulerian spray in which the new model is implemented. The results are compared to the reference Lagrangian simulation as well as the MK results. Through the qualitative and quantitative results, the AG is found to be a predictive method for the description of moderately inertial particles and is a good candidate for complex simulations in realistic configurations where small scale PTC occurs. Finally, within the framework of industrial turbulence simulations a fully kinetic Large Eddy Simulation formalism is derived based on the AG model. This strategy of directly applying the filter on the kinetic level is helpful to devise realizability conditions. Preliminary results for the AG-LES model are evaluated in 2D, in order to investigate the sensitivity of the LES result on the subgrid closures.