La distillation dite « charentaise » est une méthode de distillation discontinue permettant la production d’eaux-de-vie de cognac à partir de vin. La distillation est effectuée dans des alambics en cuivre avec un chauffage à feu nu qui peut être favorable à la génération de composés volatils. Or, les précurseurs et les mécanismes de réactions formant ces composés, pendant la distillation, restent mal connus. La première partie de l’étude consiste à identifier les composés volatils formés au cours de la distillation charentaise. Le bilan matière effectué sur de nombreux composés lors du procédé de distillation a révélé que 2 esters, 3 aldéhydes, 3 terpènes et 12 norisoprénoïdes étaient générés. Par la suite, deux distillations utilisant un mode de chauffe différent (feu nu et vapeur) ont été menées sur un alambic pilote. Le but était d’évaluer l’impact du mode de chauffe sur la génération en composés volatils. L’étude a montré que le mode de chauffe a peu d’effet sur la génération en composés volatils. De plus, les composés formés sont similaires à ceux formés lors de la distillation en alambic traditionnel. Enfin, le troisième axe de l’étude porte sur la caractérisation de l’hydrolyse acide de l’α- terpenyl-O-β-glucopyranoside, précurseur impliqué dans la formation de l’α-terpineol, identifié comme étant généré lors de la distillation charentaise ainsi qu’en alambic pilote. Le suivi de la dégradation dans un réacteur hermétiquement fermé de l’α-terpenyl-O-β-glucopyranoside et dans des conditions représentatives du vin durant la distillation montre l’hydrolyse pour former l’α- terpineol, le trans-terpin et son isomère. Cette étude a révélé que l’hydrolyse du précurseur est favorisée en milieu aqueux et suit une cinétique d’ordre 1 / The « charentaise » distillation is a batch distillation method that allows the production of cognac spirits from wine. The distillation is performed in copper alembics through a direct open flame heating that could favor the formation of volatile compounds. The first part of this study consists in identifying the volatile compounds formed during the distillation of cognac spirits. The mass balance performed on volatile compounds revealed that 2 esters, 3 aldehydes, 3 terpenes and 12 norisoprenoids were generated. Thereafter, two distillations using a different heating mode (direct open flame and steam) were conducted on a small-scale alembic. The goal was to assess the impact of the heating mode on the formation of volatile compounds. The study showed that the mode of heating has little effect on volatile compounds’ generation. Moreover, the compounds formed are similar to the ones during the traditional distillation of cognac spirit. Finally, the third part of the study focuses on the characterization of the acid hydrolysis of α-terpenyl-O-β-glucopyranoside which is the suspected precursor to be involved in the formation of α-terpineol, identified as generated during the charentaise distillation and in small-scale distillations. α-terpenyl-O-β- glucopyranoside was placed in two representative model solutions corresponding to the initial wine and the stillage. Both solutions were exposed to 100 °C in a closed reactor system. Results showed that the hydrolysis of the precurseur formed α- terpineol, trans-terpin and its isomere that seems to be 4-(2-hydroxypropan-2-yl)-1- methylcyclohexan-2-ol). Data also revealed that the hydrolysis of the precursor follows a first order reaction model ant that an aqueous media promotes the formation of trans-terpin.

En 1972, l’article 732 du Code criminel canadien a été introduit, permettant aux personnes condamnées à un emprisonnement maximal de 90 jours de purger leur peine de façon discontinue. Cette sentence revêt un statut particulier puisqu’elle s’inscrit à la fois dans deux types de mesures prévues par le système pénal, soit les peines carcérales et alternatives. À partir du cadre conceptuel de la rationalité pénale moderne (RPM), la présente étude vise à identifier, au-delà des considérations législatives prévues à l’article 732 du Code criminel, de quelle façon le système judiciaire motive l’octroi de la peine de prison discontinue. Plus précisément, nous cherchons à distinguer les motifs juridiques mobilisés par les juges pour justifier la nécessité d’incarcérer un individu sur certains jours de la semaine seulement. Par l’analyse de contenu, 77 décisions de la Cour d’appel du Québec ont été codifiées pour faire ressortir les motifs juridiques se rapportant à la rationalité pénale moderne et ceux qui, au contraire, présentent des innovations cognitives. Les vocabulaires de motifs présents dans les arrêts peuvent être classés en quatre grandes catégories. Tout d’abord, nous retrouvons plusieurs décisions dans lesquelles les juges mobilisent des motifs juridiques s’inscrivant directement dans les théories de la peine constitutives de la RPM, soit pour justifier la nécessité de révoquer la peine discontinue au profit d’une sentence d’incarcération ferme ou pour démontrer le caractère suffisamment sévère de la prison intermittente vis-à-vis des objectifs généraux de dissuasion et de dénonciation. Ensuite, plusieurs juges motivent leurs décisions en s’appuyant sur les principes et objectifs de détermination de la peine codifiés à l’article 718 du Code criminel. Ces motifs, en s’appuyant sur des considérations plus individualisées, contribuent parfois à prolonger explicitement les principes de la RPM et parfois à diluer leur rigueur, sans pour autant proposer des raisonnements novateurs. En troisième lieu, nous retrouvons certaines décisions dans lesquelles les juges parviennent à sortir du cadre conceptuel de la RPM pour innover et proposer des vocabulaires de motifs alternatifs, mais en demeurant toujours encadrés par l’obligation de punir qui caractérise le droit criminel moderne. Enfin, nous avons identifié quelques cas non usuels, c’est-à-dire des décisions dans lesquelles les juges s’appuient sur des considérations législatives ou constitutionnelles qui ne traduisent ni des idées de la RPM, ni des innovations cognitives.

La première partie de ce travail est consacré à la démonstration d'un théorème d'existence et d'unicité de la solution du système de Maxwell dans le cas général, où les coefficients sont des tenseurs symétriques définis positifs, qui dépendent d'un façon non régulière de la variable d'espace. Dans ces conditions, le milieu de propagation pourrait être aussi bien isotrope qu'anisotrope. Dans la seconde partie, nous nous sommes intéressés à l'étude et au développpement de plusieurs méthodes numériques dans un domaine isotrope où les coefficients peuvent être discontinus; nous avons étudié deux méthodes de type volumes finies, une basé sur un calcul de flux décentrés, et l'autre basée sur un calcul de flux centrés. Nous avons également adapté une méthode d'éléments finis dite Galerkin Discontinue, et enfin une méthode hybride volumes finies / différences finies avec une étude de stabilité de cette dernière. Pour des raisons géométriques, nous avons choisi les éléments du maillage comme volumes d'intégration. De nombreuses validations et comparisons numériques ont montré que ces méthodes sont bien adaptées au cas hétérogène. Néanmoins, il semble que la méthode volumes finis avec flux centrés et une discrétisation temporelle de type saute-mouton est la plus optimale en terme de compromis entre qualité des résultats et le coût de calcul.

Le développement du trafic aérien autour des agglomérations est limité notamment par les nuisances sonores qu'il engendre. Il est donc indispensable de parvenir à comprendre l'origine de ces phénomènes afin d'en réduire les effets. Les hélicoptères, qui sont amenés à survoler les agglomérations à basse altitude, font ainsi l'objet de nombreuses études. Le but principal de notre étude est de développer des méthodes numériques permettant d'évaluer le rôle des instabilités elliptiques de tourbillons dans la production du BWI. La première partie de cette thèse s'intéresse aux caractéristiques du bruit produit par l'interaction d'une pale avec un tourbillon instable. Pour ce faire, une simulation des instabilités elliptiques se développant dans une paire de tourbillons analytiques corotatifs parallèles rectilignes uniformes est d'abord réalisée. Ensuite, un outil numérique, basé sur la théorie d'Amiet, est développé pour obtenir les perturbations de pression sur pales issus d'une interaction entre une pale et le tourbillon instable. Un bon accord est trouvé entre ces résultats et les données expérimentales. Dans une deuxième partie nous nous intéressons à l'influence des caractéristiques géométriques du tourbillon d'extrémité de pale jeune isolé (coeur du tourbillon elliptique, trajectoire courbe) sur le développement des instabilités elliptiques. Pour cela le profil d'un tourbillon d'extrémité de pale a été calculé au moyen du code volumes finis. Une étude linéaire a ensuite été menée en considérant cette solution comme champ porteur, au moyen d'une méthode Euler linéarisée harmonique développée et validée au cours de cette thèse.

Le travail présenté ici avait pour but le développement d'un schéma de type Galerkin discontinu (GD) d'ordre élevé pour la résolution des équations de la dynamique des gaz écrites dans un formalisme Lagrangien total, sur des maillages bi-dimensionnels totalement déstructurés. A cette fin, une méthode progressive a été utilisée afin d'étudier étape par étape les difficultés numériques inhérentes à la discrétisation Galerkin discontinue ainsi qu'aux équations de la dynamique des gaz Lagrangienne. Par conséquent, nous avons développé dans un premier temps des schémas de type Galerkin discontinu jusqu'à l'ordre trois pour la résolution des lois de conservation scalaires mono-dimensionnelles et bi-dimensionnelles sur des maillages déstructurés. La particularité principale de la discrétisation GD présentée est l'utilisation des bases polynomiales de Taylor. Ces dernières permettent, dans le cadre de maillages bi-dimensionnels déstructurés, une prise en compte globale et unifiée des différentes géométries. Une procédure de limitation hiérarchique, basée aux noeuds et préservant les extrema réguliers a été mise en place, ainsi qu'une forme générale des flux numériques assurant une stabilité globale L2 de la solution. Ensuite, nous avons tâché d'appliquer la discrétisation Galerkin discontinue développée aux systèmes mono-dimensionnels de lois de conservation comme celui de l'acoustique, de Saint-Venant et de la dynamique des gaz Lagrangienne. Nous avons noté au cours de cette étude que l'application directe de la limitation mise en place dans le cadre des lois de conservation scalaires, aux variables physiques des systèmes mono-dimensionnels étudiés provoquait l'apparition d'oscillations parasites. En conséquence, une procédure de limitation basée sur les variables caractéristiques a été développée. Dans le cas de la dynamique des gaz, les flux numériques ont été construits afin que le système satisfasse une inégalité entropique globale. Fort de l'expérience acquise, nous avons appliqué la discrétisation GD mise en place aux équations bi-dimensionnelles de la dynamique des gaz, écrites dans un formalisme Lagrangien total. Dans ce cadre, le domaine de référence est fixe. Cependant, il est nécessaire de suivre l'évolution temporelle de la matrice jacobienne associée à la transformation Lagrange-Euler de l'écoulement, à savoir le tenseur gradient de déformation. Dans le travail présent, la transformation résultant de l'écoulement est discrétisée de manière continue à l'aide d'une base Éléments Finis. Cela permet une approximation du tenseur gradient de déformation vérifiant l'identité essentielle de Piola. La discrétisation des lois de conservation physiques sur le volume spécifique, le moment et l'énergie totale repose sur une méthode Galerkin discontinu. Le schéma est construit de sorte à satisfaire de manière exacte la loi de conservation géométrique (GCL). Dans le cas du schéma d'ordre trois, le champ de vitesse étant quadratique, la géométrie doit pouvoir se courber. Pour ce faire, des courbes de Bézier sont utilisées pour la paramétrisation des bords des cellules du maillage. Nous illustrons la robustesse et la précision des schémas mis en place à l'aide d'un grand nombre de cas tests pertinents, ainsi que par une étude de taux de convergence.

Cette thèse s'inscrit dans un projet d'innovation duale RAPID financé par DGA/DS/MRIS et appelé GREAT faisant intervenir la société Axessim, l'ONERA, INRIA, l'IRMA et le CEA. Ce projet a pour but la mise en place d'une solution industrielle de simulation électromagnétique basée sur une méthode Galerkin Discontinue (GD) parallèle sur maillage hexaédrique. Dans un premier temps, nous établissons un schéma numérique adapté à un système de loi de conservation. Nous pouvons ainsi appliquer cette approche aux équations de Maxwell, mais également à tout système hyperbolique. Dans un second temps, nous mettons en place une parallélisation à deux niveaux de ce schéma. D'une part, les calculs sont parallélisés sur carte graphique au moyen de la bibliothèque OpenCL. D'autre part, plusieurs cartes graphiques peuvent être utilisées, chacune étant pilotée par un processus MPI. De plus, les communications MPI et les calculs OpenCL sont asynchronisés permettant d'obtenir une forte accélération. / This thesis is part of a dual innovation project funded by RAPID DGA/DS/MRIS and called GREAT involving Axessim company, ONERA, INRIA, IRMA and the CEA. This project aims at the establishment of an industrial solution of electromagnetic simulation based on a method Discontinuous Galerkin (DG) on parallel hexahedral mesh. First, we establish a numerical scheme adapted to a conservation law system. We can apply this approach to the Maxwell equations but also to any hyperbolic system. In a second step, we set up a two-level parallelization of this scheme. On the one hand, the calculations are parallelized on graphics card using the OpenCL library. On the other hand, multiple graphics cards can be used, each driven by a MPI process. In addition, MPI communications and OpenCL computations are launched asynchronously in order to obtain a strong acceleration.

L'objectif principal de cette thèse est l'étude des problèmes et des applications qu'ils se développent dans le domaine de la nanophotonique. Plus précisément, nous considérons les structures de métaux nobles où les modèles de dispersion locaux sont insuffisants et la non-localité doit être incluse dans le modèle. Ici, le système physique sous-jacent est typiquement modélisé comme des équations de Maxwell couplées à des lois de dispersion spatio-temporelles dans le régime des longueurs d'onde optiques. Bien que les solutions analytiques puissent être dérivées pour un petit nombre de problèmes, cela n'est généralement pas possible pour les dispositifs du monde réel, qui présentent souvent des géométries complexes et des compositions de matériaux. Suite à une analyse rigoureuse des propriétés physiques et mathématiques du modèle continu original, nous proposons une méthode de type à éléments finis d'ordre élevé pour discrétiser le modèle continu dans l'espace et le temps. Les méthodes discontinues Galerkin (DG) sont bien établies pour la discrétisation spatiale des équations de Maxwell. Cette thèse prolonge les travaux antérieurs sur les systèmes couplés des équations de Maxwell et les lois de dispersion spatiale. Nous utilisons des méthodes explicites de Runge-Kutta (RK) d'ordre élevé pour la discrétisation temporelle. L'intégration temporelle RK garantit un ordre de convergence espace-temps élevé du schéma entièrement discret, qui repose sur un schéma de preuve de convergence. Parallélisme MPI (Message Passing Interface), éléments curvilignes et PML (Perfectly Matched Layers) autour des aspects d'implémentation et d'évaluation des performances dans le cadre du logiciel développé à Inria Sophia Antipolis-Méditerannée (DIOGENES). La méthode développée est appliquée à de nombreuses simulations nanophotoniques réelles de dispositifs où des observables tels que la réflexion, la section transversale (CS) et la spectroscopie de perte d'énergie électronique (EELS) sont étudiés. Entre autres, nous élaborons une feuille de route pour un étalonnage expérimental robuste du modèle de dispersion non local linéarisé basé sur la solution de problèmes inverses et la quantification d'incertitude (UQ) des paramètres géométriques stochastiques. Nous avons également amélioré les accords de simulations numériques non locales et les résultats expérimentaux pour la résonance des plasmons d'espacement des nano-cubes d'argent. Cela démontre la pertinence de simulations non locales précises. / The main objective of this thesis is the study of problems and applications as they arise in the field of nanophotonics. More speci cally, we consider noble metal structures where local dispersion models are insu cient and nonlocality has to be included in the model. Here, the underlying physical system is typically modeled as Maxwell’s equations coupled to spatio- temporal dispersion laws in the regime of optical wavelengths. While analytical solutions can be derived for a small number of problems, this is typically not possible for real-world devices, which often feature complicated geometries and material compositions. Following a rigorous analysis of the physical and mathematical properties of the original continuous model, we propose a high order finite element type method for discretizing the continuous model in space and time. Discontinuous Galerkin (DG) methods are well established for the spatial discretization of Maxwell’s equations. This thesis extends previous work on the coupled systems of Maxwell’s equations and spatial dispersion laws. We use explicit high-order Runge-Kutta (RK) methods for the subsequent time discretiz- ation. RK time integration guarantees a high space-time convergence order of the fully-discrete scheme, which is underpinned by a sketch of a convergence proof. Message Passing Interface (MPI) parallelization, curvilinear elements and Perfectly Matched Layers (PMLs) round of implementation aspects and performance assessments in the scope of the Software developed at Inria Sophia Antipolis-Méditerannée (DIOGENeS). The developed method is applied to numerous real-world nanophotonics simulations of devices where observables like re ectance, Cross Section (CS) and Electron Energy Loss Spectroscopy (EELS) are studied. Inter alia, we elaborate a roadmap for a robust experimental calibration of the linearized nonlocal disper- sion model based on the solution of inverse problems and Uncertainty Quanti cation (UQ) of stochastic geometric parameters. We also find improved agreements of nonlocal numerical simulations and exper- imental results for the gap-plasmon resonance of silver nano-cubes. This demonstrates the relevance of accurate nonlocal simulations.

Les données de la littérature concernant le vieillissement et le survieillissement des alliages plomb-calcium sont bien souvent incomplètes et contradictoires. Ceci est certainement lié aux difficultés expérimentales rencontrées pour observer les transformations et le fait que celles-ci sont nombreuses, entraînant une certaine confusion dans les travaux des auteurs. De plus, de faibles écarts dans les conditions d'élaboration et dans la composition des alliages peuvent avoir de l'importance sur leur comportement. Le présent travail, réalisé pour le CEA-DAM, nous a permis d'obtenir un ensemble de diagrammes TTT plus réalistes et plus justes que ceux donnés dans la littérature. Les moyens techniques mis en place (en particulier, la préservation de la chaîne du froid, condition indispensable pour garantir la répétabilité des résultats), nous ont permis d'étudier plus particulièrement les premières transformations et de mieux maîtriser les cinq étapes de vieillissement et de survieillissement. Nos travaux ont permis de déterminer clairement les cinétiques et les mécanismes des transformations. Ce travail constitue un approfondissement de la compréhension des phénomènes de vieillissement et de survieillissement de ces alliages.

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'objectif général de la thèse est la simulation numérique des mouvements forts du sol dûs aux séismes. Les déformations importantes du sol engendrent des comportements nonlinéaires dans les couches superficielles. L'apport principal de la thèse est la prise en compte de la nonlinéarité des milieux dans un contexte éléments finis Galerkin discontinus. Différentes lois de comportement sont implémentées et analysées. Le cas particulier du modèle élastoplastique de Masing-Prandtl-Ishlinskii-Iwan (MPII) est approfondi. Cette étude est divisée en deux parties. Une première qui vise à poser la structure du problème en présentant les équations et modèles utilisés pour décrire les mouvements du sol. Dans cette partie nous présentons également la méthode d’approximation spatiale Galerkin Discontinue ainsi que les différents schémas temporels que nous avons considérés. Une attention particulière est portée sur la complexité algorithmique du modèle nonlinéaire élastoplastique MPII en vue de réduire le temps de calcul des simulations. La deuxième partie est dédiée aux applications numériques. Ces applications sont réparties en trois catégories distinctes. 1) Nous nous intéressons toutd’abord à la configuration unidimensionnelle où une seule onde de cisaillement est propagée. Dans ce contexte, un flux numérique décentré est établi et des applications aux cas nonlinéaire élastique et nonlinéaire élastoplastique sont étudiées. Une solution analytique concernant le cas nonlinéaire élastique est proposée, ce qui permet de réaliser une étude numérique de convergence. 2) Le problème unidimensionnel étendu aux trois composantes du mouvement est étudié et utilisé comme un premier pas vers le 3D compte tenu du couplage entre les ondes de cisaillement et de compression. Nous nous intéressons ici à des signaux synthétiques et réels. L’application d’une méthode permettant de réduire significativement le temps de calcul du modèle élastoplastique est détaillée. 3) Une configuration tridimensionnelle est examinée. Après différentes applications de vérification en milieu linéaire, deux cas d’étude élastoplastique sont analysés. Une première sur un mode propre d’un cube puis une seconde sur un milieu plus réaliste composé d’un bassin hémisphérique à couches sédimentaires ayant un comportement élastoplastique / The general objective of this thesis is the numerical simulation of strong ground motions due to earthquakes. Significant deformations of the soil generate nonlinear behaviors in the superficial layers. The main contribution of this work is to take into account the nonlinearity of the media in a discontinuous Galerkin finite elements context. Different constitutive laws are implemented and analyzed. The particular case of theMasing-Prandtl-Ishlinskii-Iwan (MPII) elastoplastic model is looked at in-depth. This study is divided into two parts. A first one that aims at defining the framework of the problem by presenting the equations and models used to describe the soil motion. In this part we also present the Galerkin Discontinuous spatial approximation method as well as the different temporal schemes that we considered. Particular attention is paid to the algorithmic complexity of the nonlinear elastoplastic MPII model in order to reduce the computation time of simulations. The second part is dedicated to numerical applications. These applications are divided into three distinct categories. 1) We are first interested in the one-dimensional configuration where a single shear wave is propagated. In this context, an upwind numerical flux is established and applications to nonlinear elastic and nonlinear elastoplastic cases are studied. Ananalytical solution concerning the nonlinear elastic case is proposed, which makes it possible to carry out a numerical study of convergence. 2) The one-dimensional problem extended to the three components of the motion is studied and used as a first step towards 3D applications considering the coupling between the shear and compression waves. We are interested here in synthetic and real input signals. The application of a method that significantly reduces the calculation time of the elastoplastic model