A predição de estruturas de proteínas (Protein Structure Prediction PSP) é um problema computacionalmente complexo. Para tratar esse problema, modelos simplificados de proteínas, como o Modelo HP, têm sido empregados para representar as conformações e Algoritmos Evolutivos (AEs) são utilizados na busca por soluções adequadas para PSP. Entretanto, abordagens utilizando AEs muitas vezes não tratam adequadamente as soluções geradas, prejudicando o desempenho da busca. Neste trabalho, é apresentada uma formulação multiobjetivo para PSP em Modelo HP, de modo a avaliar de forma mais robusta as conformações produzidas combinando uma avaliação baseada no número de contatos hidrofóbicos com a distância entre os monômeros. Foi adotado o Algoritmo Evolutivo Multiobjetivo em Tabelas (AEMT) a fim de otimizar essas métricas. O algoritmo pode adequadamente explorar o espaço de busca com pequeno número de indivíduos. Como consequência, o total de avaliações da função objetivo é significativamente reduzido, gerando um método para PSP utilizando Modelo HP mais rápido e robusto / Protein Structure Prediction (PSP) is a computationally complex problem. To overcome this drawback, simplified models of protein structures, such as the HP Model, together with Evolutionary Algorithms (EAs) have been investigated in order to find appropriate solutions for PSP. EAs with the HP Model have shown interesting results, however, they do not adequately evaluate potential solutions by using only the usual metric of hydrophobic contacts, hamming the performance of the algorithm. In this work, we present a multi-objective approach for PSP using HP Model that performs a better evaluation of the solutions by combining the evaluation based on the number of hydrophobic contacts with the distance among the hydrophobic amino acids. We employ a Multi-objective Evolutionary Algorithm based on Sub-population Tables (MEAT) to deal with these two metrics. MEAT can adequately explore the search space with relatively low number of individuals. As a consequence, the total assessments of the objective function is significantly reduced generating a method for PSP using HP Model that is faster and more robust

Submitted by JOSIANE SANTOS DE OLIVEIRA (josianeso) on 2018-02-07T12:55:00Z No. of bitstreams: 1 Vladimir da Costa Alves_.pdf: 1430205 bytes, checksum: a7c55af3297cf449b13a07261f5a2586 (MD5) / Made available in DSpace on 2018-02-07T12:55:00Z (GMT). No. of bitstreams: 1 Vladimir da Costa Alves_.pdf: 1430205 bytes, checksum: a7c55af3297cf449b13a07261f5a2586 (MD5) Previous issue date: 2017-11-30 / Nenhuma / Este trabalho tem o objetivo de verificar a existência de ciclos na indústria da construção civil no Brasil e analisar seus efeitos em outras atividades da economia brasileira, bem como verificar a sua relevância no que se refere ao PIB deste país. O período analisado corresponde aos anos entre 1996 e 2016. A motivação para pesquisar este segmento reside na relevância da indústria da construção civil na economia brasileira, tanto pela sua magnitude na participação do PIB nacional quanto pelo seu grau de encadeamento na cadeia produtiva e, ainda, pelo posicionamento na atividade econômica como um todo. A metodologia utilizada para atingir os objetivos propostos foi o uso do filtro HP e do modelo VAR a partir dos dados publicados pelo IBGE. O filtro HP permitiu verificar o nível de correlação dos ciclos da construção civil com outros segmentos selecionados. Já o modelo VAR permitiu constatar os impactos específicos mais significativos de um choque de 1% no nível de atividade na própria construção civil e indústria extrativa mineral. De um modo geral, os resultados mostram que existe uma correlação mais forte com a indústria extrativa mineral, a indústria da transformação e o PIB total no Brasil. / This paper aims to verify the existence of cycles in the Brazilian construction industry and to analyze its effects on other activities of the Brazilian economy, as well as to verify its relevance regarding the GDP of this country. The period analyzed corresponds to the years between 1996 and 2016. The motivation to research this segment lies in the relevance of the civil construction industry in the Brazilian economy, both for its magnitude in the participation of the national GDP and for its degree of linkage in the productive chain and, still, by positioning in the economic activity as a whole. The methodology used to reach the proposed objectives was the use of the HP filter and the VAR model from the data published by IBGE. The HP filter allowed to verify the level of correlation of the construction cycles with other selected segments. On the other hand, the VAR model showed the most significant specific impacts of a 1% shock on the level of activity in the civil construction and mining industry. Overall, the results show that there is a stronger correlation with mineral extractive industry, manufacturing industry and total GDP in Brazil.

La magnéto-tellurique (MT) (Cagniard 1953, Tikhonov 1950) est une technique d'exploration de la Terre basée sur des mesures de champs électromagnétiques (EM). Une source naturelle (non artificielle) harmonique en temps et située dans l'ionosphère (Weaver 1994) produit un champ EM régi par les équations de Maxwell. Les champs électromagnétiques sont enregistrés par plusieurs récepteurs placés sur la surface de la Terre. Ces mesures sont utilisées pour produire une image du sous-sol à partir d'un procédé d'inversion utilisant des méthodes numériques. Nous utilisons la méthode hp-FEM résultant d'une extension du travail de Demkowicz 2005. Nous avons développé un logiciel qui résout, pour la première fois, le problème MT avec des éléments finis auto-adaptatifs. La méthode hp-FEM permet des raffinements locaux, à la fois en taille h et en ordre p sur les éléments, ce qui est un avantage notoire puisque la combinaison de ces deux types de critères permet de mieux capter la présence de singularités, fournissant ainsi des erreurs de discrétisation faible. C'est donc une méthode très précise dont la convergence est exponentielle (Gui and Babuska 1986, Babuska and Guo 1996). En raison des défis d'implémentation encore non résolus (Demkowicz et al. 2002) et de la complexité technique des calculs hp-FEM en 3D, nous nous limitons, dans ce travail, à des calculs en 1D et 2D.Le domaine de calcul est tronqué par un matériau absorbant (Perfectly Matched Layer PML, Berenger 1994), qui est conçu pour s'adapter automatiquement aux propriétés physiques des matériaux. En particulier, il s'ajuste efficacement à l'interface air-sol, où le contraste entre la conductivité des matériaux atteint jusqu'à seize ordres de grandeur. Dans cette thèse, nous présentons également des résultats préliminaires pour la mise en place d'une technique dimensionnelle adaptative plus connue sous le nom de DAM (Dimensionally Adaptive Method (DAM)). Lorsque la distribution de la résistivité du sous-sol dépend de multiples variables spatiales, une analyse correcte de la dimensionnalité (Ledo 2005, Martí et al. 2009, Weaver and Agarwal 2000) rend parfois possible de considérer les différentes régions avec des dimensions spatiales différentes. Par exemple, il est parfois possible d’interpréter la distribution comme une formation unidimensionnelle plus quelques hétérogénéités en 2D (ou 3D). Basée sur cette interprétation, la DAM tire profit d’une telle situation. Ainsi, l'idée principale de cette méthode est d'effectuer l'adaptativité sur la dimension spatiale en commençant par un problème de faible dimension et en utilisant les résultats obtenus pour minimiser le coût des problèmes de dimension supérieure. Nous commençons l'inversion avec un modèle 1D. Les résultats de ce problème d'inversion 1D sont utilisés comme information a priori sur les modèles de dimension supérieure. Un avantage fondamental de cette approche est que nous pouvons utiliser les solutions des problèmes de dimension inférieure précédemment calculées comme composantes du terme de régularisation associé à un problème de dimension supérieure afin d'augmenter la robustesse de l'inversion. Cette thèse propose également une analyse numérique rigoureuse de divers aspects des problèmes MT. En particulier, nous avons: (a) étudié l'effet de la source, (b) effectué une analyse fréquentielle de sensibilité, (c) illustré l'augmentation du taux de convergence lorsque l'adaptativité hp est employée, (d) séparé les effets 1D et 2D dans la solution numérique et (e) exploré l'intérêt de considérer différentes variables pour effectuer l'inversion. / The magnetotelluric (MT) method is a passive exploration technique that aims at estimating the resistivity distribution of the Earth's subsurface, and therefore at providing an image of it. This process is divided into two different steps. The first one consists in recording the data. In a second step, recorded measurements are analyzed by employing numerical methods. This dissertation focuses in this second task. We provide a rigorous mathematical setting in the context of the Finite Element Method (FEM) that helps to understand the MT problem and its inversion process. In order to recover a map of the subsurface based on 2D MT measurements, we employ for the first time in Mts a multi-goal oriented self adaptive hp-Finite Element Method (FEM). We accurately solve both the full formulation as well as a secondary field formulation where the primary field is given by the solution of a 1D layered media. To truncate the computational domain, we design a Perfectly Matched Layer (PML) that automatically adapts to high-contrast material properties that appear within the subsurface and on the air-ground interface. For the inversion process, we develop a first step of a Dimensionally Adaptive Method (DAM) by considering the dimension of the problem as a variable in the inversion. Additionally, this dissertation supplies a rigorous numerical analysis for the forward and inverse problems. Regarding the forward modelization, we perform a frequency sensitivity analysis, we study the effect of the source, the convergence of the hp-adaptivity, or the effect of the PML in the computation of the electromagnetic fields and impedance. As far as the inversion is concerned, we study the impact of the selected variable for the inversion process, the different information that each mode provides,and the gains of the DAM approach.

Submitted by Jaqueline Silva (jtas29@gmail.com) on 2014-09-30T20:52:20Z No. of bitstreams: 2 dissertação - Júnio César Mendes Simão - 2013.pdf: 1594402 bytes, checksum: fada126ea8a1b1a7b56dfeab5037ff06 (MD5) license_rdf: 23148 bytes, checksum: 9da0b6dfac957114c6a7714714b86306 (MD5) / Approved for entry into archive by Jaqueline Silva (jtas29@gmail.com) on 2014-09-30T21:42:39Z (GMT) No. of bitstreams: 2 dissertação - Júnio César Mendes Simão - 2013.pdf: 1594402 bytes, checksum: fada126ea8a1b1a7b56dfeab5037ff06 (MD5) license_rdf: 23148 bytes, checksum: 9da0b6dfac957114c6a7714714b86306 (MD5) / Made available in DSpace on 2014-09-30T21:42:39Z (GMT). No. of bitstreams: 2 dissertação - Júnio César Mendes Simão - 2013.pdf: 1594402 bytes, checksum: fada126ea8a1b1a7b56dfeab5037ff06 (MD5) license_rdf: 23148 bytes, checksum: 9da0b6dfac957114c6a7714714b86306 (MD5) Previous issue date: 2013-03-01 / Coordenação de Aperfeiçoamento de Pessoal de Nível Superior - CAPES / This paper deals with the learning of some content from Financial Mathematics using the HP-12C calculator, to prepare citizens for both managing its own nances and for the labor market. The matter will be addressed through practical applications in a pleasant and e ective, with the intention to make students able to resolve nancial situations that present themselves in a critical way, not being carried away by misleading proposals. At the end we present some practical applications where advertisements will be analyzed. / Este trabalho trata da aprendizagem de alguns conteúdos de Matemática Financeira com a utilização da calculadora HP-12C, de modo a preparar cidadãos tanto para administração das suas próprias nanças quanto para o mercado de trabalho. O assunto será tratado através de aplicações práticas de um modo agradável e e caz, com a intenção de tornar os alunos aptos a resolver situações nanceiras que se apresentam de modo crítico, não sendo levados por propostas enganosas. Ao nal iremos apresentar algumas aplicações práticas onde anúncios publicitários serão analisados.

La modélisation de la nature (minéralogie, composition) et de la dynamique du manteau terrestre nécessite la connaissance des propriétés de ses minéraux constitutifs à haute pression (HP, 1-135 GPa) et haute température (HT, 1000-3000K). En particulier, les grandeurs thermodynamiques (enthalpies de formation, entropies, volume, ... ) contrôlent la stabilité des minéraux dans l'espace pression-température. Dans une première approche dite "macroscopique", les capacités calorifiques de la forstérite et de ses analogues, du pyrope et du spinelle sont mesurées par calorimétrie à HT. Les données de dilatation thermique et de compressibilité sont tirées de mesures de volume à HP et HT. On dispose alors d'un ensemble de données qui permet de mettre en évidence, pour la plupart de ces minéraux, un comportement anharmonique à HT. Dans une seconde approche "microscopique", les propriétés thermodynamiques d'un minéral sont calculées par modélisation vibrationnelle à partir de ses spectres infrarouge et Raman. Une méthode originale de mesure de l'anharmonicité est développée, fondée sur l'enregistrement des spectres vibrationnels à HP et HT. On montre alors pourquoi et comment les valeurs des capacités calorifiques s'écartent à HT de la limite prédite par une théorie harmonique. Cette méthode "microscopique" originale, qui ne nécessite que quelques milligrammes de matière, est ainsi développée et testée sur des matériaux pour lesquels elle peut être confrontée aux données "macroscopiques" (e.g. calorimétrie). Elle pourra dans l'avenir être appliquée aux phases de très haute pression du manteau terrestre, dont la métastabilité et les faibles quantités synthétisées empêchent la caractérisation par les méthodes classiques "macroscopiques".

The origin and evolution of granites has been widely studied because granitoid rocks constitute a major portion of the Earth ́s crust. The formation of granitic magma is, besides temperature mainly triggered by the water content of these rocks. The presence of water in magmas plays an important role due to the ability of aqueous fluids to change the chemical composition of the magma. The exsolution of aqueous fluids from melts is closely linked to a fractionation of elements between the two phases. Then, aqueous fluids migrate to shallower parts of the Earth ́s crust because of it ́s lower density compared to that of melts and adjacent rocks. This process separates fluids and melts, and furthermore, during the ascent, aqueous fluids can react with the adjacent rocks and alter their chemical signature. This is particularly impor- tant during the formation of magmatic-hydrothermal ore deposits or in the late stages of the evolution of magmatic complexes. For a deeper insight to these processes, it is essential to improve our knowledge on element behavior in such systems. In particular, trace elements are used for these studies and petrogenetic interpretations because, unlike major elements, they are not essential for the stability of the phases involved and often reflect magmatic processes with less ambiguity. However, for the majority of important trace elements, the dependence of the geochemical behavior on temperature, pressure, and in particular on the composition of the system are only incompletely or not at all experimentally studied. Former studies often fo- cus on the determination of fluid−melt partition coefficients (Df/m=cfluid/cmelt) of economically interesting elements, e.g., Mo, Sn, Cu, and there are some partitioning data available for ele- ments that are also commonly used for petrological interpretations. At present, no systematic experimental data on trace element behavior in fluid−melt systems as function of pressure, temperature, and chemical composition are available. Additionally, almost all existing data are based on the analysis of quenched phases. This results in substantial uncertainties, particularly for the quenched aqueous fluid because trace element concentrations may change upon cooling. The objective of this PhD thesis consisted in the study of fluid−melt partition coefficients between aqueous solutions and granitic melts for different trace elements (Rb, Sr, Ba, La, Y, and Yb) as a function of temperature, pressure, salinity of the fluid, composition of the melt, and experimental and analytical approach. The latter included the refinement of an existing method to measure trace element concentrations in fluids equilibrated with silicate melts di- rectly at elevated pressures and temperatures using a hydrothermal diamond-anvil cell and synchrotron radiation X-ray fluorescence microanalysis. The application of this in-situ method enables to avoid the main source of error in data from quench experiments, i.e., trace element concentration in the fluid. A comparison of the in-situ results to data of conventional quench experiments allows a critical evaluation of quench data from this study and literature data. In detail, starting materials consisted of a suite of trace element doped haplogranitic glasses with ASI varying between 0.8 and 1.4 and H2O or a chloridic solution with m NaCl/KCl=1 and different salinities (1.16 to 3.56 m (NaCl+KCl)). Experiments were performed at 750 to 950◦C and 0.2 or 0.5 GPa using conventional quench devices (externally and internally heated pressure vessels) with different quench rates, and at 750◦C and 0.2 to 1.4 GPa with in-situ analysis of the trace element concentration in the fluids. The fluid−melt partitioning data of all studied trace elements show 1. a preference for the melt (Df/m < 1) at all studied conditions, 2. one to two orders of magnitude higher Df/m using chloridic solutions compared to experiments with H2O, 3. a clear dependence on the melt composition for fluid−melt partitioning of Sr, Ba, La, Y, and Yb in experiments using chloridic solutions, 4. quench rate−related differences of fluid−melt partition coefficients of Rb and Sr, and 5. distinctly higher fluid−melt partitioning data obtained from in-situ experiments than from comparable quench runs, particularly in the case of H2O as starting solution. The data point to a preference of all studied trace elements for the melt even at fairly high salinities, which contrasts with other experimental studies, but is supported by data from studies of natural co-genetically trapped fluid and melt inclusions. The in-situ measurements of trace element concentrations in the fluid verify that aqueous fluids will change their composition upon cooling, which is in particular important for Cl free systems. The distinct differences of the in-situ results to quench data of this study as well as to data from the literature signify the im- portance of a careful fluid sampling and analysis. Therefore, the direct measurement of trace element contents in fluids equilibrated with silicate melts at elevated PT conditions represents an important development to obtain more reliable fluid−melt partition coefficients. For further improvement, both the aqueous fluid and the silicate melt need to be analyzed in-situ because partitioning data that are based on the direct measurement of the trace element content in the fluid and analysis of a quenched melt are still not completely free of quench effects. At present, all available data on element complexation in aqueous fluids in equilibrium with silicate melts at high PT are indirectly derived from partitioning data, which involves in these experiments assumptions on the species present in the fluid. However, the activities of chemical components in these partitioning experiments are not well constrained, which is required for the definition of exchange equilibria between melt and fluid species. For example, the melt-dependent variation of partition coefficient observed for Sr imply that this element can not only be complexed by Cl− as suggested previously. The data indicate a more complicated complexation of Sr in the aqueous fluid. To verify this hypothesis, the in-situ setup was also used to determine strontium complexation in fluids equilibrated with silicate melts at desired PT conditions by the application of X-ray absorption near edge structure (XANES) spectroscopy. First results show a strong effect of both fluid and melt composition on the resulting XANES spectra, which indicates different complexation environments for Sr. / Die Entstehung und Entwicklung von Graniten steht seit Jahrzehnten im Fokus vieler geologischer Studien, da sich die Erdkruste zu großen Teilen aus granitoiden Gesteinen zusammensetzt. Von besonderer Bedeutung für die Bildung von granitischen Schmelzen ist neben der Temperatur, der Wassergehalt der Schmelze, da dieser Parameter die chemische Zusammensetzung der Schmelze entscheidend verändern kann. Die Entmischung wässriger Fluide aus Schmelzen führt zur Neuverteilung von Elementen zwischen diesen Phasen. Bedingt durch die geringere Dichte des wässrigen Fluids im Vergleich zur Schmelze und dem Nebengestein, beginnt dieses aus tieferen Erdschichten aufzusteigen. Damit verknüpft ist nicht nur eine räumliche Trennung von Schmelze und Fluid, sondern auch die Alterierung des Nebengestein. Dieser Prozess ist insbesondere bei der Bildung von magmatisch-hydrothermalen Lagerstätten und in späten Entwicklungsstadien magmatischer Komplexe wichtig. Für ein detailliertes Verständnis dieser Prozesse ist es notwendig, das Elementverhalten in solchen Systemen in Abhängigkeit von Parametern wie Temperatur, Druck und chemischer Zusammensetzung des Systems experimentell zu untersuchen, und Elementverteilungskoeffizienten als Funktion dieser Variablen zu bestimmen. Für die Untersuchungen sind insbesondere Spurenelemente geeignet, da diese im Gegensatz zu Hauptelementen nicht essentiell für die Stabilität weiterer auftretender Phasen sind, aber sehr sensibel auf Änderungen intensiver Variablen reagieren können. Zudem werden bei geochemischen Mineral- und Gesteinsanalysen viele Spurenelemente, Spurenelementverhältnisse, und Spurenelementisotope als petrogenetische Indikatoren verwendet, d.h. diese Daten liefern Informationen darüber, wann und in welcher Tiefe und bei welchen chemischen Bedingungen ein Gestein gebildet worden ist, und welche weiteren Prozesse es auf dem Weg zur Erdoberfläche durchlaufen hat. Allerdings sind für vie- le Spurenelemente die Abhängigkeiten der Verteilung zwischen Fluiden und Schmelzen von intensiven Variablen nicht, oder nur unzureichend experimentell untersucht worden. Zusätzlich dazu basiert die Mehrheit der experimentell gewonnenen Verteilungskoeffizienten und deren Interpretation, insbesondere hinsichtlich der Elementkomplexierung im Fluid, auf der Analyse von schnell abgekühlten Phasen. Bisher ist nicht geklärt, ob solche Analysen repräsentativ sind für die Zusammensetzungen der Phasen bei hohen Drücken und Temperaturen. Das Ziel dieser Studie war die Erarbeitung eines experimentellen Datensatzes zur Spu- renelementverteilung zwischen granitischen Schmelzen und wässrigen Fluiden in Abhängigkeit von der Schmelzzusammensetzung, der Salinität des Fluids, des Drucks und der Temperatur. Ein Hauptanliegen der Arbeit bestand in der Weiterentwicklung einer experimentellen Methode bei welcher der Spurenelementgehalt im Fluid in-situ, d.h. unter hohen Drücken und Temperaturen, und im Gleichgewicht mit einer silikatischen Schmelze bestimmt wird. Die so gewonnenen Daten können anschließend mit den Resultaten von Abkühlexperimenten vergli- chen werden, um diese und auch Literaturdaten kritisch zu bewerten. Die Daten aller unter- suchten Spurenelemente dieser Arbeit (Rb, Sr, Ba, La, Y und Yb) zeigen: 1. unter den untersuchten Bedingungen eine Präferenz für die Schmelze unabhängig von der chemischen Zusammensetzung von Schmelze und Fluid, Druck oder Temperatur, 2. die Verwendung von chloridhaltigen Fluiden kann die Verteilungskoeffizienten um 1 bis 2 Größenordnungen anheben und 3. für die Verteilungskoeffizienten von Sr, Ba, La, Y und Yb eine starke Abhängigkeit von der Schmelzzusammensetzung im chloridischen System. Der Vergleich der Daten der verschiedenen Methoden zeigt, dass insbesondere für chloridfreie Fluide große Diskrepanzen zwischen den in-situ Daten und Analysen von abgeschreckten Proben bestehen. Dieses Ergebnis beweist eindeutig, dass beim Abschrecken der Proben Rückreaktionen stattfinden, und dass Daten, welche auf Analysen abgeschreckter Fluide basieren, nur eingeschränkt verwendet werden sollten. Die Variation der Verteilungskoeffizienten von Sr, Ba, La, Yb, und Y als Funktion der Schmelzzusammensetzung ist entweder auf eine Änderung der Komplexierung im Fluid und/oder einen anderen veränderten Einbau dieser Elemente in die Schmelze zurückzuführen. Daher wurde im Rahmen dieser Arbeit erstmals versucht, die Elementkomplexierung in silikatischen Fluiden direkt bei hohen Temperaturen und Drücken zu bestimmen. Die Daten für Sr zeigen, dass abhängig von der Schmelzzusammensetzung unterschiedliche Komplexe stabil sein müssen.

Le stockage géologique du CO2 dans les aquifères salins représente une stratégie prometteuse pour la réduction des émissions de CO2 anthropiques. Ce type de stockage requiert des connaissances fondamentales afin d'évaluer les scénarios d'injection, d'estimer la capacité des réservoirs et les risques de fuite. C'est dans ce contexte que des outils microfluidiques haute pression/haute température ont été développés afin d'étudier différents mécanismes liés aux technologies de stockage de CO2. Les laboratoires géologiques sur puce (GLoCs) permettent de mener des expériences à des conditions de pression et de température typiques des réservoirs (25 < T < 50°C, 50 < p < 10 MPa) et d'en mimer des propriétés géologiques.Ce manuscrit présente dans un premier temps les stratégies de fabrication des GLoCs et l'accès à leurs caractéristiques (porosité et perméabilité). La détection du CO2 directement au sein de GLoCs grâce à l'implémentation de fibres optiques par spectroscopie infrarouge est ensuite présentée, ainsi que la mise en œuvre de la technique de laminographie X menées à l'ESRF pour le suivi de dissolution de carbonates dans des microcanaux. Le manuscrit explicite ensuite les investigations, menées avec des GLoCs concernant les différents mécanismes de piégeage du CO2 à l'échelle du pore. La visualisation directe et le traitement d'image ont permis de suivre l'évolution des phases de CO2 et des phases aqueuses au sein du GLoC et les mécanismes de déplacement et de saturation. Enfin, les travaux en cours sont exposés, comme les expériences de drainage avec des saumures réactives ou la formation d'hydrates au sein de milieux poreux. / CO2 geological storage in deep saline aquifers represents a mediation solution for reducing the anthropogenic CO2 emissions. Consequently, this kind of storage requires adequate scientific knowledge to evaluate injection scenarios, estimate reservoir capacity and assess leakage risks. In this context, we have developed and used high pressure/high temperature microfluidic tools to investigate the different mechanisms associated with CO2 geological storage in deep saline aquifers. The silicon-Pyrex 2D porous networks (Geological Labs On Chips) can replicate the reservoir p,T conditions (25 < T < 50°C, 50 < p < 10 MPa), geological and topological properties. This thesis manuscript first highlights the strategies developed during this work to fabricate the GLoCs and to access to global characteristics of our porous media such as porosity and permeability, which are later compared to numerical modelling results. The carbon dioxide detection in GLoCs mimicking p,T conditions of geological reservoirs by using the direct integration of optical fiber for IR spectroscopy is presented. I then detail the strategies for following the dissolution of carbonates in GLoCs with X-rays laminography experiments.Then, the manuscript focuses on the use of GLoCs to investigate each CO2 trapping mechanism at the pore scale. The direct optical visualization and image processing allow us to follow the evolution of the injected CO2/aqueous phase within the reservoir, including displacement mechanisms and pore saturation levels.Eventually, I present the ongoing works such as experiments with reactive brines and hydrates formations in porous media.

Cette thèse traite des méthodes d’éléments finis Galerkin discontinus d’ordre élevé pour la résolution d’équations aux dérivées partielles, avec un intérêt particulier pour l’équation de transport des neutrons. Nous nous intéressons tout d’abord à une méthode de pré-traitement de matrices creuses par blocs, qu’on retrouve dans les méthodes Galerkin discontinues, avant factorisation par un solveur multifrontal. Des expériences numériques conduites sur de grandes matrices bi- et tri-dimensionnelles montrent que cette méthode de pré-traitement permet une réduction significative du ’fill-in’, par rapport aux méthodes n’exploitant pas la structure par blocs. Ensuite, nous proposons une méthode d’éléments finis Galerkin discontinus, employant des éléments d’ordre élevé en espace comme en angle, pour résoudre l’équation de transport des neutrons. Nous considérons des solveurs parallèles basés sur les sous-espaces de Krylov à la fois pour des problèmes ’source’ et des problèmes aux valeur propre multiplicatif. Dans cet algorithme, l’erreur est décomposée par projection(s) afin d’équilibrer les contraintes numériques entre les parties spatiales et angulaires du domaine de calcul. Enfin, un algorithme HP-adaptatif est présenté ; les résultats obtenus démontrent une nette supériorité par rapport aux algorithmes h-adaptatifs, à la fois en terme de réduction de coût de calcul et d’amélioration de la précision. Les valeurs propres et effectivités sont présentées pour un panel de cas test industriels. Une estimation précise de l’erreur (avec effectivité de 1) est atteinte pour un ensemble de problèmes aux domaines inhomogènes et de formes irrégulières ainsi que des groupes d’énergie multiples. Nous montrons numériquement que l’algorithme HP-adaptatif atteint une convergence exponentielle par rapport au nombre de degrés de liberté de l’espace éléments finis. / We consider high order discontinuous-Galerkin finite element methods for partial differential equations, with a focus on the neutron transport equation. We begin by examining a method for preprocessing block-sparse matrices, of the type that arise from discontinuous-Galerkin methods, prior to factorisation by a multifrontal solver. Numerical experiments on large two and three dimensional matrices show that this pre-processing method achieves a significant reduction in fill-in, when compared to methods that fail to exploit block structures. A discontinuous-Galerkin finite element method for the neutron transport equation is derived that employs high order finite elements in both space and angle. Parallel Krylov subspace based solvers are considered for both source problems and $k_{eff}$-eigenvalue problems. An a-posteriori error estimator is derived and implemented as part of an h-adaptive mesh refinement algorithm for neutron transport $k_{eff}$-eigenvalue problems. This algorithm employs a projection-based error splitting in order to balance the computational requirements between the spatial and angular parts of the computational domain. An hp-adaptive algorithm is presented and results are collected that demonstrate greatly improved efficiency compared to the h-adaptive algorithm, both in terms of reduced computational expense and enhanced accuracy. Computed eigenvalues and effectivities are presented for a variety of challenging industrial benchmarks. Accurate error estimation (with effectivities of 1) is demonstrated for a collection of problems with inhomogeneous, irregularly shaped spatial domains as well as multiple energy groups. Numerical results are presented showing that the hp-refinement algorithm can achieve exponential convergence with respect to the number of degrees of freedom in the finite element space

Cette thèse aborde différents aspects de la résolution numérique des Equations aux Dérivées Partielles.Le premier chapitre est consacré à l'étude de la méthode Mixed High-Order (MHO). Il s'agit d'une méthode mixte de dernière génération permettant d'obtenir des approximations d'ordre arbitraire sur maillages généraux. Le principal résultat obtenu est l'équivalence entre la méthode MHO et une méthode primale de type Hybrid High-Order (HHO).Dans le deuxième chapitre, nous appliquons la méthode MHO/HHO à des problèmes issus de la mécanique des fluides. Nous considérons d'abord le problème de Stokes, pour lequel nous obtenons une discrétisation d'ordre arbitraire inf-sup stable sur maillages généraux. Des estimations d'erreur optimales en normes d'énergie et L2 sont proposées. Ensuite, nous étudions l'extension au problème d'Oseen, pour lequel on propose une estimation d'erreur en norme d'énergie où on trace explicitement la dépendance du nombre de Péclet local.Dans le troisième chapitre, nous analysons la version hp de la méthode HHO pour le problème de Darcy. Le schéma proposé permet de traiter des maillages généraux ainsi que de faire varier le degré polynomial d'un élément à l'autre. La dépendance de l'anisotropie locale du coefficient de diffusion est tracée explicitement dans l'analyse d'erreur en normes d'énergie et L2.La thèse se clôture par une ouverture sur la réduction de problèmes de diffusion à coefficients variables. L'objectif consiste à comprendre l'impact du choix de la formulation (mixte ou primale) utilisée pour la projection sur l'espace réduit sur la qualité du modèle réduit. / This Ph.D. thesis deals with different aspects of the numerical resolution of Partial Differential Equations.The first chapter focuses on the Mixed High-Order method (MHO). It is a last generation mixed scheme capable of arbitrary order approximations on general meshes. The main result of this chapter is the equivalence between the MHO method and a Hybrid High-Order (HHO) primal method.In the second chapter, we apply the MHO/HHO method to problems in fluid mechanics. We first address the Stokes problem, for which a novel inf-sup stable, arbitrary-order discretization on general meshes is obtained. Optimal error estimates in both energy- and L2-norms are proved. Next, an extension to the Oseen problem is considered, for which we prove an error estimate in the energy norm where the dependence on the local Péclet number is explicitly tracked.In the third chapter, we analyse a hp version of the HHO method applied to the Darcy problem. The resulting scheme enables the use of general meshes, as well as varying polynomial orders on each face.The dependence with respect to the local anisotropy of the diffusion coefficient is explicitly tracked in both the energy- and L2-norms error estimates.In the fourth and last chapter, we address a perspective topic linked to model order reduction of diffusion problems with a parametric dependence. Our goal is in this case to understand the impact of the choice of the variational formulation (primal or mixed) used for the projection on the reduced space on the quality of the reduced model.

nÃo hà / Nesta dissertaÃÃo abordamos questÃes concernentes a MatemÃtica Financeira a partir de dois enfoques, o teÃrico metodolÃgico e o uso da tecnologia da mÃquina calculadora HP 12C. Nesse sentido, nosso objetivo foi de contribuir com o uso da calculadora HP 12C para os ensino fundamental, mÃdio e superior, proporcionando a esses alunos mais interesse e motivaÃÃo para compreender melhor um assunto tÃo importante para sua vida humana e profissional. Para tanto elaboramos um tutorial com operaÃÃes elementares sob o contexto da MatemÃtica Financeira proporcionando uma aprendizagem consistente para resolver situaÃÃes prÃticas e ter o discernimento entre comprar à vista e à prazo. Dessa forma o uso adequado do tutorial facilita um avanÃo na aprendizagem do contexto da MatemÃtica Financeira, auxiliando em uma avaliaÃÃo ou orÃamento financeiro sob o ponto de vista de prevà resultados de investimentos em um negÃcio. Abordamos tambÃm nesse contexto os sistemas de amortizaÃÃes de financiamentos atualmente em uso, capitalizaÃÃo simples e composta, fluxo de caixa, valor presente lÃquido, taxa interna de retorno e fundo de pensÃo. / The Financial Mathematics treated in this issue presents two points of view, the methodological theoretical and the use of the machine technology. Its purpose is to contribute with the use of the calculator HP 12C in the elementary, secondary and higher education, providing for the students more interest and motivation to learn subject like this in their professional lives. So, we elaborate a tutorial with elementary operations about the structure of Financial Mathematics that help to solve practice situations and obtain the discernment between buying the spot or forward. The adequate use of this tutorial, help to solve practice situations and understand how to work the Financial Mathematics and allowing to foresee the results of investments in business. We also treat in this tutorial the amortization of financing systems currently in use, simple and compound capitalization, cash flow, net present value, internal rate of return and pension fund.