Le stockage géologique du CO2 est l'une des diverses technologies étant explorées afin de réduire les émissions de carbone atmosphérique des processus industriels (i.e. combustion de l'énergie fossile). L'une des spécifiques caractéristiques de l'injection du CO2 en profondeur reste la possibilité de réactions géochimiques (dissolution-précipitation) entre la saumure réactive mobile (e.g. eau de formation enrichie en CO2) et la roche encaissante durant l'évolution spatiale et temporelle du CO2, conduisant à des modifications dans la structure des pores et par conséquent dans les propriétés d'écoulement du réservoir (e.g. la perméabilité k). Donc, ces changements structuraux peuvent largement contrôler l'injectivité, ainsi que le champ de pression dans le réservoir et aussi la propagation du CO2. Il demeure ainsi crucial d'explorer les changement dans les propriétés de réservoirs (e.g. structurales et hydrodynamiques) induits durant une injection de CO2 et explicitement les relations existantes entre eux (e.g. k ou surface réactive-Sr versus porosité- , k versus hétérogénéité de la roche), afin de développer des outils de modélisation prédictive des processus de transport et réactionnels se produisant durant une injection de CO2 et d'évaluer de façon fiable les risques. Dans le cas des réservoirs carbonatés, l'application des modèles prédictifs de transport réactif demeure toujours un enjeu, car contrainte par la forte hétérogénéité en leur sein ainsi que par l'incertitude dans la cinétique de réactions des minéraux carbonatés dans ce contexte. Dans cette optique, nous avons réalisé des expériences de percolation à travers des échantillons de roches carbonatées dans les conditions thermodynamiques de stockage en profondeur (T = 100°C et P =12 MPa). L'évolution de la perméabilité est suivie au cours des expériences ; et la variation de la porosité est calculée à partir des résultats d'analyses chimiques au ICP-AES des fluides de sortie échantillonnés. L'investigation des modifications apportées à la structure des pores est réalisée par le biais de la Micro-Tomographie haute résolution à rayon X, acquise au synchrotron de Grenoble (e.g. ESRF). Dépendant du régime de dissolution, contrôlé par la fabrique de la roche réservoir et la composition chimique de la saumuré chargée en CO2 (e.g. PCO2 engagée), on a observé qu'une modification de la structure de la roche peut soit améliorer soit détériorer (résultat atypique en contexte de dissolution) la valeur de la perméabilité k. Mots clés : Stockage géologique du CO2, transport, réactions géochimiques, structure des pores, propriétés hydrodynamiques, expériences de percolation de CO2, micro-tomographie à rayon X. / Geological storage of CO2 is one of diverse technologies being explored to reduce atmospheric carbon from industrial processes (i.e. fossil fuel combustion). One of the specific features of CO2 injection is the possibility of geochemical reactions (dissolution – precipitation) between mobile reactive brine (e.g. formation water enriched in CO2) and the host rock during the spatial and temporal evolution of CO2. That leads to modifications in the pore structure which in turn change the flow dynamics of the reservoir (e.g. the permeability k). Then, theses structural modifications can largely control the injectivity, so that the pressure field in the reservoir and also the CO2 propagation. Accordingly, it is crucial to explore the changes in the reservoir properties (e.g. structural and hydrodynamic) induced during a CO2 injection and specially the relationships between them (e.g. k or reactive surface-Sr versus porosity- , k versus rock heterogeneity), for developing predictive modelling tools of the transport and reaction processes occurring during a CO2 injection and reliable risk assessment. In the case of carbonate rocks, the application of the predictive models of transport and reaction is still challenging, because of their high heterogeneity so that the incertitude in the reaction kinetics of carbonate minerals. From this perspective, we realized brine-enriched in CO2 percolation experiments through carbonate rock samples in thermodynamic conditions expected during CO2 injection in deep reservoirs (T = 100°C et P =12 MPa). The permeability changes k(t) is monitored during the experiments and the porosity variation is calculated from chemical analyses of the sampled outlet fluids, using ICP-EAS. The pore structure modifications are investigated from high resolution X ray micro tomography images acquired from the synchrotron of Grenoble (ESRF). Depending to the dissolution regime, controlled by the reservoir rock fabric and the chemical composition of the brine (e.g. PCO2), we observed that a modification of pore structure can either improve (atypical result in dissolution context) or impair the value of the permeability k. Keywords: CO2 geological storage, transport, geochemical reactions, pore structure, hydrodynamic properties, brine enriched in CO2 percolation experiments, X ray microtomography.

Ce travail entre dans le cadre de l’étude des interactions eau-roche dans le cas du stockage du CO2 en milieu géologique. Un intérêt particulier est accordé aux hétérogénéités des paramètres associés aux phénomènes géochimiques. Ces hétérogénéités peuvent s’observer à différentes échelles: celle des grains (les minéraux présentent des défauts de cristallinité et des impuretés), et l’échelle centimétrique/pluri-décamétrique. En particulier, les paramètres thermodynamiques (logK) et de cinétique chimique (dans ce travail nous avons considéré le produit de la constante cinétique k par la surface spécifique S soit kS comme "paramètre de cinétique chimique") sont connus à partir des expériences de laboratoire pour des échantillons de quelques centimètres de dimension, alors que l’on s’intéresse aux réactions minéralogiques à l’échelle des réservoirs.Nous avons évalué les caractéristiques géostatistiques de la variabilité spatiale après réaction à travers des simulations de transport réactif dans lesquelles différents paramètres (logK et kS) sont perturbés avec une première variabilité imposée. Une combinaison de deux approches est ainsi abordée : déterministe et géostatistique. Le code du transport-réactif COORES (IFP-EN et Ecole nationale supérieure des mines de Saint-Etienne) a été utilisé pour les simulations déterministes et le système géochimique étudié concerne la dissolution du diopside avec précipitation de minéraux secondaires comme la calcite et la magnésite.Après analyse par la méthode des plans d’expériences, les résultats montrent qu’une corrélation spatiale élevée combinée avec une grande variance de dispersion des minéraux favorise une réactivité importante des minéraux lorsqu’on perturbe le paramètre de cinétique chimique kS. Par ailleurs une vitesse d’injection élevée accélère le processus de dissolution du minéral étudié. La variabilité spatiale du paramètre thermodynamique n’a cependant pas d’effet significatif sur les résultats, le système se comporte comme dans le cas homogène. Du point de vue de l’homogénéisation du paramètre kS, on retrouve l’influence de l’historique de dissolution. / The present work is based on the study of water-rock interactions in the case of CO2 storage in geological media. Particular attention is devoted to heterogeneities at different observation scales geochemical phenomena. These heterogeneities can be observed at different scales: the grain (mineral crystallinity present defects and impurities), and the centimeter scale / multi- decametric (rocks are heterogeneous at different scales). In particular, the thermodynamic parameters logK and chemical kinetics kS (in this work we considered the product of the rate constant k by the specific surface area S is kS as "chemical kinetics parameter") are known from laboratory experiments to a few centimeters in size, while we are interested in mineralogical reactions across tanks.We propose to evaluate the geostatistical characteristics of the local variability after reaction through simulations of reactive transport on a small scale in which various parameters (logK and kS) are perturbed with a first spatial variability imposed. A combination of both approaches is discussed: deterministic and geostatistical for the study of geochemical problems at different scales. The reactive transport code - COORES (IFP - EN and Ecole nationale supérieure des mines de Saint -Etienne) was used for deterministic simulations and the geochemical system studied concerns the dissolution of diopside with precipitation of secondary minerals such as calcite and magnesite.After analysis by the method of design of experiments, the results show that high spatial correlation variance combined with high dispersion of minerals promotes a high reactivity when minerals chemically disturbing is the kinetic parameter kS. In addition, a high velocity injection accelerates the dissolution of the mineral studied. However, the effect of spatial variability of the thermodynamic parameter, did not significantly affect the results, the system behaves as in the homogeneous case. From the standpoint of homogenizing the parameter kS, include the influence of the history of dissolution.

Les méthodes Direct Simulation Monte Carlo (DSMC) sont utilisées par Ariane group pour calculer les torseurs d'efforts aérodynamiques et les flux thermiques sur les engins spatiaux pour des écoulements hypersoniques en milieu raréfié.Afin de pouvoir caractériser la dislocation des étages de lanceurs et donc l'empreinte de retombée de débris, une modélisation précise des mécanismes générateurs de flux thermiques est nécessaire. Les réactions chimiques étant dimensionnantes pour le calcul du flux thermique, l'objectif de cette thèse est de développer l'outil de calcul avec la méthode DSMC nommé IEMC de manière à pouvoir prendre en compte les écoulements réactifs.Deux modèles de chimie sont mis en place pour pouvoir prendre en compte la totalité des réactions. Après leur vérifications sur des cas élémentaires, ils sont appliqués et validés sur des cas tests de rentrée pour différentes atmosphères. Les différents modèles considérés sont testés afin d'évaluer leur influence. Les modèles de chimie dépendent de nouveaux paramètres d'entrée, dont les valeurs numériques sont incertaines. Une étude de quantification de leur incertitude est menée et a permis de vérifier que les grandeurs de sorties de la simulation avec un écoulement réactif, notamment le flux thermique, n'est que peu impacté par ces paramètres incertains. / Direct Simulation Monte Carlo (DSMC) methods are used in Ariane group to compute aerodynamic forces and moments and heat fluxes on space objects for hypersonic flows in rarefied regimes.To caracterise the dislocation of the stages and the debris footprints, a precise modelisation of the mechanism that contribute to the heat flux is necessary. The contribution of the chemical reactions is important for the determination of the heat flux. The purpose of this thesis is to develop the in house IEMC tool using the DSMC method so that it can compute reactive flows.The different steps of the developments are presented in this work. The first step is the presentation, implementation and verification of two different chemistry models. They are validated for simulations on real test cases. Different models are tested in order to evaluate their effect. Chemical models implemented in the code depend on new input parameters, whose numerical data are uncertain. Using a uncertainty quantification study, it is shown that the output data of the reactive simulation, especially the heat flux, is weakly impacted by the tested uncertain parameters.

L’eau des circuits d’une centrale nucléaire est purifiée à l’aide de résines échangeuses d’ions. La prédiction de leurs performances constitue une aide importante pour l’exploitation de ces réacteurs. Les méthodes du transport réactif sont particulièrement adaptées pour cela et constituent la base du code OPTIPUR, dédié à la modélisation de ces résines. Le travail présenté comporte trois axes principaux. Le premier est l’intégration d’une limitation au transfert de masse dans une colonne de résines échangeuses d’ions, avec une mobilité spécifique à chaque espèce chimique, dans le cadre d’un découplage des calculs de chimie et de transport. Ce modèle permet, sans paramètre ajustable, de reproduire assez fidèlement une série d’expériences réalisées précédemment par le CEA. Le second axe concerne les aspects numériques du transport réactif, avec l’utilisation de la méthode d’Anderson pour accélérer la convergence du couplage chimie-transport dans un schéma itératif. En utilisant les informations issues des itérations précédentes et sans modification majeure du code, la robustesse et les temps de calcul ont pu être nettement améliorés. La troisième thématique abordée est celle de l’équilibre d’échange d’ions. Les bases d’un modèle prenant en considération l’évolution de l’humidité de la résine, ainsi que son élasticité sont proposées ; les interactions entre groupes fonctionnels, contre-ions et eau sont considérées comme des équilibres chimiques. Les constantes d’équilibre sont ajustées à partir de mesures de la teneur en eau de la résine à différentes pressions de vapeur d’eau. Finalement, des coefficients de sélectivité apparents peuvent être calculés et comparés aux mesures disponibles. / In nuclear power plants, the water contained in different circuits is purified by passing through ion exchange resins. Prediction of the performance of these resins is an important help to the plant operators. To this end, the method of reactive transport modeling are well suited and is the basis of the OPTIPUR code that was designed to model the resins. The work presented in this manuscript covers three main aspects. The first one is the integration of a limitation to mass transfer in a ion exchange deep bed, taking into account a specific mobility for each chemical species, in the context of separated calculations for chemistry and transport. This model was shown to reproduce experimental data, without adjustable parameters. The second part of this work deals with the numerical aspects of reactive transport modelling. A method developped by Anderson was used to accelerate the convergence of the chemistry transport coupling in an iterative scheme. Using the information from previous iterations, and without major changes in the code, calculation times were largely decreased, as well as the number of calculations failures. The third topic is ion exchange equilibrium. The basis of a model that takes into account the change in the water content of the resin and its elasticity are described. The interactions between the fonctional groups, the counterions and water are considered as chemical reactions. The corresponding equilibrium constants are fitted to measurements of the water content of the resin at different relative humidity. Finaly, the selectivity coefficients can be calculated and compared to litterature values.

Le moulage RTM à haute cadence est un procédé de fabrication composite prometteur qui satisfait les exigences de l’industrie automobile pour produire des pièces structurelles complexes avec un temps de cycle court. Cependant, les réductions de temps de cycle sont un véritable défi. Dans ce procédé, une résine est injectée avec mélange en dans la cavité d’un moule contenant un renfort fibreux. Ce flux de résine réactive génère des schémas d’écoulement complexes et des couplages thermo-chimio-rhéologiques forts. En raison de la grande sensibilité de la résine et des temps de cycle serrés, la prédiction de la stratégie d’injection optimale est très difficile et très coûteuse à mener expérimentalement. Le travail réalisé a donc poursuivi deux objectifs: 1. Identifier et quantifier expérimentalement les mécanismes influençant le procédé RTM réactif avec mélange en tête et 2. Développer une méthode de simulation numérique en vue d’introduire les mécanismes identifiés dans le logiciel industriel PAM COMPOSITE développé par ESI Group. L’identification et la quantification des mécanismes ont été réalisées grâce à des investigations expérimentales et numériques. Un nouveau montage expérimental a été développé pour l’étude du mécanisme de stockage de résines intra-mèche grâce à des observations aux échelles macro- et microscopiques. De plus, une méthode numérique a été développée pour simuler l’écoulement réactif de la résine dans des matériaux à simple et à double échelle de porosité. Cette méthode a permis d’étudier les mécanismes locaux difficiles à mesurer expérimentalement et de préparer le transfert vers le logiciel industriel d’ESI. / High Speed Resin Transfer Molding (RTM) is a promising composite manufacturing process fitting automotive industry requirements to produce complex structural parts with a perspective of short cycle times. However, cycle time reductions are a real challenge. In this process, a resin mixed on-line with curing agents is injected in the cavity of a mold containing a fibrous reinforcement. This flow of reactive resin generates acomplex flow pattern and strong thermo-chemo rheological couplings. Due to the high sensitivity of the resin cure, and the tight cycle times, prediction of the optimal injection strategy is very difficult and very expensive to conduct experimentally. In this context, two goals where followed in this work: 1. Identify and quantify experimentally the mechanisms, related to the process or to the reinforcement, influencing the reactive RTMprocess with on-line mixing and 2. Develop a numerical simulation method in a view of introducing the identified mechanisms in the industrial software PAMCOMPOSITE developed by ESI Group. Identification and quantification of the mechanisms were realized thanks to experimental investigations and numerical simulations. A new experimental setup has been developed for the investigation of the mechanism of intra-tow resin storage through macro-scale and micro-scale observations. Additionally, a numerical method has been developed to simulate the reactive flow of a resin in single and dual scale porous materials. This method allowed both to investigate local mechanisms difficult to study experimentally and prepare the transfer to the industrial software of ESI.

Nos travaux de recherches sont des travaux supports pour les gestionnaires de l’hôpital Robert Pax de Sarreguemines, et plus généralement de tout centre hospitalier pour développer des approches centrées sur le patient. Nous nous sommes inspirés des approches centrées sur le produit issues du domaine industriel qui ont été proposées pour répondre aux exigences croissantes de gestion des produits dans un environnement de plus en plus incertain. Par analogie, les systèmes de production de soins centrés patient peuvent être assimilés aux systèmes de production de biens centrés produit. Cependant, il ne faut pas perdre de vue la spécificité des systèmes de production de soins : le facteur humain. Cette spécificité est à l’origine de leur caractère complexe, aléatoire et imprévisible. Par ailleurs, les approches de pilotage, que ce soit dans le milieu industriel ou hospitalier, sont majoritairement des modèles mathématiques et des modèles de simulation utilisant une approche de gestion centrée sur une ou plusieurs ressources considérées comme critiques. C’est pourquoi il nous a paru judicieux d’opter pour une approche centrée patient basée sur le paradigme multi-agents. Nous proposons alors, un système multi-agents de pilotage réactif dynamique et distribué centré patient du parcours patient au sein des systèmes hospitaliers. L’alternative que nous proposons consiste à utiliser une approche centrée patient et basée sur les agents permettant de minimiser les délais d’attente, ainsi que la durée de séjour, et par conséquent les coûts des soins, tout en assurant un soin de qualité pour l’ensemble des patients et une meilleure gestion des ressources hospitalières / Through our research, we offer a support tool for the managers in Robert Pax hospital in Sarreguemine (France), and more generally any hospital center wishing to develop a patient centered approach. We were inspired by « product centered » approaches emerging from industrial domain which were proposed to answer the increasing requirements of products management. By analogy, health care systems patient centered can be assimilated to production systems product centered. However, it is important not to lose sight of the specificities of health care systems which is the human factor. This specificity makes this system complex, random and unpredictable. Besides, the approaches used in industrial or hospital environment, are mainly mathematical models and simulation approaches centered on one ore several resources categorized as critical. In that sense, it seemed judicious to choose a patient centered approach based on a multi-agent paradigm. We proposed a Reactive multi-agent control system of the patient flow in healthcare system. The proposed alternative is a patient centered approach allowing to minimize the patients waiting time and the length of their stay and consequently the care costs. The proposed approach also ensures the care quality and an optimal use of the hospital resources

Les travaux présentés dans cette thèse visent à produire automatiquement des prototypes de circuits de communication et de contrôle à partir de spécifications temporelles déclaratives. Partant d'un ensemble de propriétés écrites en langage PSL, nous produisons un modèle RTL synthétisable automatiquement. La méthode proposée est modulaire, contrairement aux méthodes publiées antérieurement qui étaient fondées sur la théorie des automates. Pour chaque propriété, nous produisons un composant qui observe certains opérandes et génère des chronogrammes pour les autres opérandes : le module réactif. Tout d'abord, une bibliothèque des modules réactifs primitifs a été développée pour les opérateurs FL et SERE. Pour ce faire, une relation de dépendance a été définie pour chaque opérateur : fondée sur la sémantique de l'opérateur, elle exprime la dépendance entre ses opérandes. Ensuite, la relation de dépendance de chaque opérateur est interprétée comme un composant matériel qui met en œuvre l'opérateur : c'est le module réactif primitif de l'opérateur. À l'aide de cette formalisation, nous proposons une méthode pour déterminer automatiquement quels signaux d'une propriété sont observés et lesquels sont générés. Dans le cas où il n'est pas possible de déterminer le sens du signal, un solveur est ajouté pour identifier la valeur du signal. Le solveur sert aussi à déterminer la valeur d'un signal généré par plusieurs propriétés. Le circuit final est l'interconnexion des modules réactifs et des solveurs pour l'ensemble des propriétés. Un outil prototype, SyntHorus2, qui est une extension d'HORUS, a été mis développé. Il prend les propriétés PSL comme entrées et génère le code VHDL synthétisable du circuit. En outre, il génère des propriétés complémentaires pour vérifier si l'ensemble des spécifications est cohérent et complet. La méthode est efficace et synthétise des circuits de commande en quelques secondes. Les résultats que nous avons obtenus sur des jeux d'essais classiques montrent que notre technique compile les propriétés plus efficacement que les outils prototypes qui l'ont précédée. / The work presented in this thesis aims at automatically prototype communication and control designs from declarative temporal specifications. From a set of PSL properties, we produce a synthesizable RTL design automatically. The proposed method is modular, in contrast to previously published methods that were based on automata theory. From each property, we produce a component that observes some operands and generates waveforms for the other operands: the reactant. First, a library of primitive reactants has been provided for FL and SERE operators. To this goal, a dependency relation is defined for each operator that expresses the dependency among its operands using the operator's semantics. Then, the dependency relation of each operator is interpreted as a hardware component that implements the operator: the operator's primitive reactant. Using this formalization, a method is proposed to automatically decide which signals of a property are observed and which are generated. In the cases when specifying the signal direction is not possible, a solver is implemented to identify the signal value. In addition, the way of identifying the value of the signal that is generated in several properties is addressed. The final circuit is the interconnection of the properties' reactants and solvers. A prototype tool SyntHorus2, which is an extension to HORUS, has been developed. It takes PSL properties as its inputs, and generates the synthesizable VHDL code of the circuit. In addition, it generates some complementary properties to verify if the set of specification is coherent and complete. The method is efficient, and synthesizes control circuits in a few seconds. Results obtained on classical benchmarks show that our technique compiles properties more efficiently than previous prototype tools.

Dans le cadre d’une thématique générale du Laboratoire Thermique Energétique et Procédés de Pau consacré à l’étude et à la valorisation des déchets tels que la pyrolyse, la combustion et la gazéification, la thèse se focalise sur la modélisation des phénomènes de transports de masse, de quantité de mouvement et d’énergie en milieu multiphasique multiconstituant réactif dans le cas d’un réacteur de gazéification à lit fixe. La première étape de la modélisation consiste à décrire classiquement le mouvement de chacune des phases continues par les équations de conservation de la masse, de la quantité de mouvement et de la chaleur. Bien qu’à cette échelle les mécanismes soient parfaitement décrits, le passage à la simulation impose une étape d’homogénéisation par prise de moyenne. Ce changement d’échelle, décrit dans ce travail, conduit à un système d’équations homogène à l’échelle locale. Ce modèle est appliqué au cas du procédé de gazéification de déchets de la société Europlasma nommé CHO-Power. Les simulations bidimensionnelles instationnaires du problème ont permis de mettre en lumière les différents mécanismes en présence au cours du procédé ainsi qu’un certain nombre de verrous dans l’obtention du chemin de convergence conduisant au régime permanent. / One of the thematics of the “Laboratoire Thermique Energétique et Procédés de Pau” is the study and the valorization of waste as pyrolysis, combustion and gasification for example. In this context, this work deals with modeling of heat, mass and momentum transport in a multiphase multi components reactive medium for a fixed bed gasifier. Conservation equations are first written for each phase. The macroscopic partial differential equations are expressed by integrating these microscopic conservation laws over a representative volume. This change of scale, described in this work, leads to a homogeneous system of equations. This model is applied to the case of the gasification of waste process of the Europlasma company named CHO-Power. The unsteady two-dimensional simulations of the problem allowed to highlight the different mechanisms present during the process and a number of locks in obtaining convergence path leading to the steady state.

La protection et la remédiation des ressources en eau sont un enjeu sociétal majeur, ainsi il est nécessaire de comprendre l’évolution de solutés, tels les polluants, au sein de la zone saturée et non saturée. Dans ce but, de nombreux travaux ont été consacrés à la modélisation du transport réactif en milieux poreux. Son déroulement à l’échelle de Darcy dépend des hétérogénéités microscopiques du milieu. Les modèles de réseau de pores qui simplifient la géométrie en un ensemble de pores reliés par des liens de sections constantes, permettent de se placer à une échelle mésoscopique, faisant le lien entre l’échelle porale et l’échelle de Darcy. Sur de telles formes géométriques, l’écoulement admet un traitement analytique. En ce qui concerne le transport réactif des solutés, nous proposons une solution analytique dans les liens qui permet de calculer le débit de masse entre pores. Le modèle de transport se formule alors comme un système d’équations de Volterra de secondes espèces dont les noyaux de convolution sont des séries d’exponentielles décroissantes (hormis le premier terme qui est constant). Leurs temps de relaxation sont pilotés essentiellement par le temps de dispersion td. Dans la limite où td tend vers 0 à Péclet constant, les termes transitoires des noyaux se réduisent à un Dirac, débouchant sur un premier modèle simplifié à réponse instantanée c'est-à-dire un modèle de transport quasi-statique. Dans le cas où les volumes des pores sont suffisamment grands, les noyaux se réduisent à leur premier terme. Ces formulations du transport généralisent celles de la littérature. En particulier pour des Péclet petit ou grand on retrouve respectivement les modèles usuels en régime dispersif et convectif. Numériquement, la décroissance exponentielle des noyaux permet d’optimiser le calcul des convolutions avec une précision arbitrairement fixée, réduisant drastiquement le temps de résolution. / Protection and remediation of ground water resources are a major societal challenge. It implies to understand the evolution of solutes as pollutant in the saturated and non-saturated zones. For that purpose numerous studies have been conducted for modeling the reactive transport in a porous media. At Darcy scale, the behavior of solutes depends on microscopic heterogeneity for the media. The Pore Network Models (PNM) simplifies drastically its geometry and considers pores linked by straight throats the section of which is constant. They give a description which is in between the macroscopic and the pore descriptions. With such geometry it is possible to use a Poiseuille flow modeling the flux. With respect to the reactiontransport equation, we seek the analytical solution of the CDE in throats, which in turn allows computing the mass flux in pores. The transport solution consists of a Volterra equation system. Its convolution kernels result in a summation of time function which is decreasing exponentially with time (except the first term which still constant). The time constant is driven by the diffusion time td. As td goes to zero, keeping the Peclet number fixed, each term of the summation reduces to a Dirac. The response of the system is then instantaneous. When the volume of the pore is large enough it is possible to neglect all the term of the kernel except the constant one. In the limit where the Peclet number goes to zero, usual models are recovered. Numerically, the exponential time decreasing of the kernel allow to optimize their computational time up to an arbitrary fixed precision.

Cette thèse propose une nouvelle approche du parallélisme et de la concurrence, posant les bases d'un langage de programmation à la fois sûr et "secure" (garantissant la sécurité des données), fondé sur une sémantique formelle claire et simple, tout en étant adapté aux architectures multi-cœur. Nous avons adopté le paradigme synchrone, dans sa variante réactive, qui fournit une alternative simple à la programmation concurrente standard en limitant l'impact des erreurs dépendant du temps ("data-races"). Dans un premier temps, nous avons considéré un langage réactif d'orchestration, DSL, dans lequel on fait abstraction de la mémoire (Partie 1). Dans le but de pouvoir traiter la mémoire et la sécurité, nous avons ensuite étudié (Partie 2) un noyau réactif, CRL, qui utilise un opérateur de parallélisme déterministe. Nous avons prouvé la réactivité bornée des programmes de CRL. Nous avons ensuite équipé CRL de mécanismes pour contrôler le flux d'information (Partie 3). Pour cela, nous avons d'abord étendu CRL avec des niveaux de sécurité pour les données, puis nous avons défini dans le langage étendu, SSL, un système de types permettant d'éviter les fuites d'information. Parallèlement (Partie 4), nous avons ajouté la mémoire à CRL, en proposant le modèle DSLM. En utilisant une notion d'agent, nous avons structuré la mémoire de telle sorte qu'il ne puisse y avoir de "data-races". Nous avons également étudié l'implémentation de DSLM sur les architectures multi-cœur, fondée sur la notion de site et de migration d'un agent entre les sites. L'unification de SSL et de DSLM est une piste pour un travail futur. / This thesis proposes a new approach to parallelism and concurrency, laying the basis for the design of a programming language with a clear and simple formal semantics, enjoying both safety and security properties, while lending itself to an implementation on multicore architectures. We adopted the synchronous programming paradigm, in its reactive variant, which provides a simple alternative to standard concurrent programming by limiting the impact of time-dependent errors ("data-races"). As a first step (Part 1), we considered a reactive orchestration language, DSL, which abstracts away from the memory. To set the basis for a formal treatment of memory and security, we then focussed on a reactive kernel, CRL, equipped with a deterministic parallel operator (Part 2). We proved bounded reactivity of CRL programs. Next, we enriched CRL with mechanisms for information flow control (Part 3). To this end, we first extended CRL with security levels for data. We then defined a type system on the extended language, SSL, which ensures the absence of information leaks. Finally, we added memory to CRL, as well as the notions of agent and site, thus obtaining the model DSLM (Part 4). We structured the memory in such a way that data-races cannot occur, neither within nor among agents. We also investigated the implementation of DSLM on multicore architectures, using the possibility of agent migration between sites. The unification of SSL and DSLM is left for future work.