Numerical simulations of magma chamber dynamics at Campi Flegrei, and associated seismicity, deformation and gravity changes

Understanding the complex relationships between quantities measured by volcanic monitoring network and shallow magma processes is a crucial headway for the comprehension of volcanic processes and a more realistic evaluation of the associated hazard. This question is very relevant at Campi Flegrei, a volcanic quiescent caldera immediately north-west of Napoli (Italy). The system activity shows a high fumarole release and periodic ground slow movement (bradyseism) with high seismicity. This activity, with the high people density and the presence of military and industrial buildings, makes Campi Flegrei one of the areas with higher volcanic hazard in the world. In such a context my thesis has been focused on magma dynamics due to the refilling of shallow magma chambers, and on the geophysical signals detectable by seismic, deformative and gravimetric monitoring networks that are associated with this phenomenologies. Indeed, the refilling of magma chambers is a process frequently occurring just before a volcanic eruption; therefore, the faculty of identifying this dynamics by means of recorded signal analysis is important to evaluate the short term volcanic hazard. The space-time evolution of dynamics due to injection of new magma in the magma chamber has been studied performing numerical simulations with, and implementing additional features in, the code GALES (Longo et al., 2006), recently developed and still on the upgrade at the Istituto Nazionale di Geofisica e Vulcanologia in Pisa (Italy). GALES is a finite element code based on a physico-mathematical two dimensional, transient model able to treat fluids as multiphase homogeneous mixtures, compressible to incompressible. The fundamental equations of mass, momentum and energy balance are discretised both in time and space using the Galerkin Least-Squares and discontinuity-capturing stabilisation technique. The physical properties of the mixture are computed as a function of local conditions of magma composition, pressure and temperature.The model features enable to study a broad range of phenomenologies characterizing pre and sin-eruptive magma dynamics in a wide domain from the volcanic crater to deep magma feeding zones. The study of displacement field associated with the simulated fluid dynamics has been carried out with a numerical code developed by the Geophysical group at the University College Dublin (O’Brien and Bean, 2004b), with whom we started a very profitable collaboration. In this code, the seismic wave propagation in heterogeneous media with free surface (e.g. the Earth’s surface) is simulated using a discrete elastic lattice where particle interactions are controlled by the Hooke’s law. This method allows to consider medium heterogeneities and complex topography. The initial and boundary conditions for the simulations have been defined within a coordinate project (INGV-DPC 2004-06 V3_2 “Research on active volcanoes, precursors, scenarios, hazard and risk - Campi Flegrei”), to which this thesis contributes, and many researchers experienced on Campi Flegrei in volcanological, seismic, petrological, geochemical fields, etc. collaborate. Numerical simulations of magma and rock dynamis have been coupled as described in the thesis. The first part of the thesis consists of a parametric study aimed at understanding the eect of the presence in magma of carbon dioxide in magma in the convection dynamics. Indeed, the presence of this volatile was relevant in many Campi Flegrei eruptions, including some eruptions commonly considered as reference for a future activity of this volcano. A set of simulations considering an elliptical magma chamber, compositionally uniform, refilled from below by a magma with volatile content equal or dierent from that of the resident magma has been performed. To do this, a multicomponent non-ideal magma saturation model (Papale et al., 2006) that considers the simultaneous presence of CO2 and H2O, has been implemented in GALES. Results show that the presence of CO2 in the incoming magma increases its buoyancy force promoting convection ad mixing. The simulated dynamics produce pressure transients with frequency and amplitude in the sensitivity range of modern geophysical monitoring networks such as the one installed at Campi Flegrei . In the second part, simulations more related with the Campi Flegrei volcanic system have been performed. The simulated system has been defined on the basis of conditions consistent with the bulk of knowledge of Campi Flegrei and in particular of the Agnano-Monte Spina eruption (4100 B.P.), commonly considered as reference for a future high intensity eruption in this area. The magmatic system has been modelled as a long dyke refilling a small shallow magma chamber; magmas with trachytic and phonolitic composition and variable volatile content of H2O and CO2 have been considered. The simulations have been carried out changing the condition of magma injection, the system configuration (magma chamber geometry, dyke size) and the resident and refilling magma composition and volatile content, in order to study the influence of these factors on the simulated dynamics. Simulation results allow to follow each step of the gas-rich magma ascent in the denser magma, highlighting the details of magma convection and mixing. In particular, the presence of more CO2 in the deep magma results in more ecient and faster dynamics. Through this simulations the variation of the gravimetric field has been determined. Afterward, the space-time distribution of stress resulting from numerical simulations have been used as boundary conditions for the simulations of the displacement field imposed by the magmatic dynamics on rocks. The properties of the simulated domain (rock density, P and S wave velocities) have been based on data from literature on active and passive tomographic experiments, obtained through a collaboration with A. Zollo at the Dept. of Physics of the Federici II Univeristy in Napoli. The elasto-dynamics simulations allow to determine the variations of the space-time distribution of deformation and the seismic signal associated with the studied magmatic dynamics. In particular, results show that these dynamics induce deformations similar to those measured at Campi Flegrei and seismic signals with energies concentrated on the typical frequency bands observed in volcanic areas. The present work shows that an approach based on the solution of equations describing the physics of processes within a magmatic fluid and the surrounding rock system is able to recognise and describe the relationships between geophysical signals detectable on the surface and deep magma dynamics. Therefore, the results suggest that the combined study of geophysical data and informations from numerical simulations can allow in a near future a more ecient evaluation of the short term volcanic hazard.

Studio del metodo dell'assorbimento diretto della CO2 per l'analisi del radiocarbonio

L’attuale rilevanza rappresentata dalla stretta relazione tra cambiamenti climatici e influenza antropogenica ha da tempo posto l’attenzione sull’effetto serra e sul surriscaldamento planetario così come sull’aumento delle concentrazioni atmosferiche dei gas climaticamente attivi, in primo luogo la CO2. Il radiocarbonio è attualmente il tracciante ambientale per eccellenza in grado di fornire mediante un approccio “top-down” un valido strumento di controllo per discriminare e quantificare il diossido di carbonio presente in atmosfera di provenienza fossile o biogenica. Ecco allora che ai settori applicativi tradizionali del 14C, quali le datazioni archeometriche, si affiancano nuovi ambiti legati da un lato al settore energetico per quanto riguarda le problematiche associate alle emissioni di impianti, ai combustibili, allo stoccaggio geologico della CO2, dall’altro al mercato in forte crescita dei cosiddetti prodotti biobased costituiti da materie prime rinnovabili. Nell’ambito del presente lavoro di tesi è stato quindi esplorato il mondo del radiocarbonio sia dal punto di vista strettamente tecnico e metodologico che dal punto di vista applicativo relativamente ai molteplici e diversificati campi d’indagine. E’ stato realizzato e validato un impianto di analisi basato sul metodo radiometrico mediante assorbimento diretto della CO2 ed analisi in scintillazione liquida apportando miglioramenti tecnologici ed accorgimenti procedurali volti a migliorare le performance del metodo in termini di semplicità, sensibilità e riproducibilità. Il metodo, pur rappresentando generalmente un buon compromesso rispetto alle metodologie tradizionalmente usate per l’analisi del 14C, risulta allo stato attuale ancora inadeguato a quei settori applicativi laddove è richiesta una precisione molto puntuale, ma competitivo per l’analisi di campioni moderni ad elevata concentrazione di 14C. La sperimentazione condotta su alcuni liquidi ionici, seppur preliminare e non conclusiva, apre infine nuove linee di ricerca sulla possibilità di utilizzare questa nuova classe di composti come mezzi per la cattura della CO2 e l’analisi del 14C in LSC.

Hydrogen production by enhanced methane reforming with membrane reactors or with CO2 capture materials

Pure hydrogen production from methane is a multi-step process run on a large scale for economic reasons. However, hydrogen can be produced in a one-pot continuous process for small scale applications, namely Low Temperature Steam Reforming. Here, Steam Reforming is carried out in a reactor whose walls are composed by a membrane selective toward hydrogen. Pd is the most used membrane material due to its high permeability and selectivity. However, Pd deteriorates at temperatures higher than 500°C, thus the operative temperature of the reaction has to be lowered. However, the employment of a membrane reactor may allow to give high yields thanks to hydrogen removal, which shifts the reaction toward the products. Moreover, pure hydrogen is produced. This work is concentrated on the synthesis of a catalytic system and the investigation of its performances in different processes, namely oxy-reforming, steam reforming and water gas shift, to find appropriate conditions for hydrogen production in a catalytic membrane reactor. The catalyst supports were CeZr and Zr oxides synthesized by microemulsion, impregnated with different noble metals. Pt, Rh and PtRh based catalysts were tested in the oxy reforming process at 500°C, where Rh on CeZr gave the most interesting results. On the opposite, the best performances in low temperature steam reforming were obtained with Rh impregnated on Zr oxide. This catalyst was selected to perform low temperature steam reforming in a Pd membrane reactor. The hydrogen removal given by the membrane allowed to increase the methane conversion over the equilibrium of a classical fixed bed reactor thanks to an equilibrium shift effect. High hydrogen production and recoveries were also obtained, and no other compound permeated through the membrane which proved to be hydrogen selective.

Synthesis, Modification and Characterization of Polymeric Membranes for CO2 Separation

This doctorate focused on the development of dense polymeric membranes for carbon capture, mostly in post combustion applications, and for natural gas sweetening. The work was supported by the European Project NANOMEMC2 funded under H2020 program. Different materials have been investigated, that rely on two main transport mechanisms: the solution-diffusion and the facilitated transport. In both cases, proper nano-fillers have been added to the matrix, in order to boost the mechanical and permselective properties of the membranes. Facilitated transport membranes were based on the use of was polyvinylamine (PVAm), as main matrix with fixed-site carriers, and L-Arginine as mobile carrier; the filler, used mostly as reinforcer, was carboxymethylated nanocellulose (cNFC). Humid test showed interesting results, and especially the blend made of PVAm/cNFC/Arg in weight ratio 27,5/27,5/45 crossed the Robeson CO2/N2 upper bound, representing current state of the art membranes, with a CO2 permeability of 271 Barrer and CO2/N2 selectivity of 70. Solution diffusion membranes were based on Pebax®2533 matrix which was added with three different graphene oxide (GO)-based materials, namely pristine GO, Porous Graphene Oxide (PGO) and a GO functionalized with polyetheramine (PEAGO). All of them provided a modest but clear increment of permeability of the Pebax matrix, from plus 2% (GO) to plus 8% (PGO), with no change in selectivity. The gas tested with this type of composites were CO2 and N2, for Post combustion capture applications. Pebax®2533 was also chemically modified, obtaining the product called “Benzoyl-P2533”, that was fully characterized, and tested in term of permeation using five gas: CO2, N2, CH4, O2, and He. Modified material showed an increment of the overall permeability of the material of a fair 10% for all gases tested, apart from helium, that increased of almost 50%.

Development and Detailed Characterization of Innovative, High-Performance Membrane Materials for CO2 Capture Processes

The scope of this dissertation is to study the transport phenomena of small molecules in polymers and membranes for gas separation applications, with particular attention to energy efficiency and environmental sustainability. This work seeks to contribute to the development of new competitive selective materials through the characterization of novel organic polymers such as CANALs and ROMPs, as well as through the combination of selective materials obtaining mixed matrix membranes (MMMs), to make membrane technologies competitive with the traditional ones. Kinetic and thermodynamic aspects of the transport properties were investigated in ideal and non-ideal scenarios, such as mixed-gas experiments. The information we gathered contributed to the development of the fundamental understanding related to phenomenon like CO2-induced plasticization and physical aging. Among the most significant results, ZIF-8/PPO MMMs provided materials whose permeability and selectivity were higher than those of the pure materials for He/CO2 separation. The CANALs featured norbornyl benzocyclobutene backbone and thereby introduced a third typology of ladder polymers in the gas separation field, expanding the structural diversity of microporous materials. CANALs have a completely hydrocarbon-based and non-polar rigid backbone, which makes them an ideal model system to investigate structure-property correlations. ROMPs were synthesized by means of the ring opening metathesis living polymerization, which allowed the formation of bottlebrush polymers. CF3-ROMP reveled to be ultrapermeable to CO2, with unprecedented plasticization resistance properties. Mixed-gas experiments in glassy polymer showed that solubility-selectivity controls the separation efficiency of materials in multicomponent conditions. Finally, it was determined that plasticization pressure in not an intrinsic property of a material and does not represent a state of the system, but rather comes from the contribution of solubility coefficient and diffusivity coefficient in the framework of the solution-diffusion model.

Coupled wellbore-reservoir simulation of thermal effects during multiphase CO2 injection

T2Well-ECO2M is a coupled wellbore reservoir simulator still under development at Lawrence Berkeley National Laboratory (USA) with the ability to deal with a mixture of H2O-CO2-NaCl and includes the simulation of CO2 phase transition and multiphase flow. The code was originally developed for the simulation of CO2 injection into deep saline aquifers and the modelling of enhanced geothermal systems; however, the focus of this research was to modify and test T2Well-ECO2M to simulate CO2 injection into depleted gas reservoirs. To this end, the original code was properly changed in a few parts and a dedicated injection case was developed to study CO2 phase transition inside of a wellbore and the corresponding thermal effects. In the first scenario, the injection case was run applying the fully numerical approach of wellbore to formation heat exchange calculation. Results were analysed in terms of wellbore pressure and temperature vertical profiles, wellhead and bottomhole conditions, and characteristic reservoir displacement fronts. Special attention was given to the thorough analysis of bottomhole temperature as the critical parameter for hydrate formation. Besides the expected direct effect of wellbore temperature changes on reservoir conditions, the simulation results indicated also the effect of CO2 phase change in the near wellbore zone on BH pressure distribution. To test the implemented software changes, in a second scenario, the same injection case was reproduced using the improved semi-analytical time-convolution approach for wellbore to formation heat exchange calculation. The comparison of the two scenarios showed that the simulation of wellbore and reservoir parameters after one year of continuous CO2 injection are in good agreement with the computation time to solve the time-convolution semi-analytical reduced. The new updated T2Well-ECO2M version has shown to be a robust and performing wellbore-reservoir simulator that can be also used to simulate the CO2 injection into depleted gas reservoirs.

Strategies for lowering the environmental impact of different electrochemical energy storage system

The continuous growth of global population brings an exponential increase on energy consumption and greenhouse gas emission in the atmosphere contributing to the increase of the planet temperature. Therefore, it is mandatory to adopt renewable energy production systems like photovoltaic or wind power: unfortunately, the main limit of these technologies is the natural intermittence of the energy sources that limits their applicability. The key enabling technology for a widespread usage of clean power sources are electrochemical energy storage systems, most commonly known as batteries. Batteries will enable the storage of energy during overproduction period and the release during low production period stabilizing the power outcome, allowing the connection to the main grid and increasing the applicability of renewable energy sources. Despite the high number of benefits that the widespread use of batteries will bring, starting from the reduction of CO2 emitted in the atmosphere, it is necessary also to take care of the environmental impact of processes and materials used for the production of electrochemical storage systems. In addition, there are many different battery systems, with different chemistries and designs that require specific strategies. Nowadays, the most part of the materials and chemicals used for battery production are toxic for humans and the environment. For this reason, this Ph.D. thesis addresses the challenging scope of lowering the environmental impact of manufacturing processes of different electrochemical energy storage systems using natural derived or low carbon footprint materials while increasing the performances with respect to commercial devices. The activities carried out during my Ph.D. cover a high number of different electrochemical storage systems involving a wide range of electrochemical processes from capacitive to faradic. New materials, different production processes and new battery design, all in view of sustainability and low environmental impact, increased the innovative and challenging aspects of this work.

Design and testing of innovative, artificial intelligence-based techniques to model and control the combustion process and to reduce the Co2 emissions in modern, high-performance, gasoline direct injection engines

This work deals with the development of calibration procedures and control systems to improve the performance and efficiency of modern spark ignition turbocharged engines. The algorithms developed are used to optimize and manage the spark advance and the air-to-fuel ratio to control the knock and the exhaust gas temperature at the turbine inlet. The described work falls within the activity that the research group started in the previous years with the industrial partner Ferrari S.p.a. . The first chapter deals with the development of a control-oriented engine simulator based on a neural network approach, with which the main combustion indexes can be simulated. The second chapter deals with the development of a procedure to calibrate offline the spark advance and the air-to-fuel ratio to run the engine under knock-limited conditions and with the maximum admissible exhaust gas temperature at the turbine inlet. This procedure is then converted into a model-based control system and validated with a Software in the Loop approach using the engine simulator developed in the first chapter. Finally, it is implemented in a rapid control prototyping hardware to manage the combustion in steady-state and transient operating conditions at the test bench. The third chapter deals with the study of an innovative and cheap sensor for the in-cylinder pressure measurement, which is a piezoelectric washer that can be installed between the spark plug and the engine head. The signal generated by this kind of sensor is studied, developing a specific algorithm to adjust the value of the knock index in real-time. Finally, with the engine simulator developed in the first chapter, it is demonstrated that the innovative sensor can be coupled with the control system described in the second chapter and that the performance obtained could be the same reachable with the standard in-cylinder pressure sensors.