Metodología de simulación del almacenamiento geológico de CO2 como contribución a la solución al problema del Cambio Climático, aplicada a la estructura BG-GE-08 (oeste de la Comunidad de Murcia)

Una de las medidas de mitigación del Cambio Climático propuestas al amparo de la ONU por el IPCC (Intergovermental Panel on Climate Change) en su ‘Informe de Síntesis 2007’ consiste en la puesta en marcha de acciones para la captación y almacenamiento de dióxido de carbono, existiendo tres tipos de formaciones geológicas idóneas para el almacenamiento geológico de este gas: yacimientos de petróleo y gas agotados, capas de carbón no explotables y formaciones salinas profundas. En el caso de las formaciones salinas profundas, el problema fundamental para llevar a cabo un estudio de almacenamiento de CO2, reside en la dificultad de obtención de datos geológicos del subsuelo en una cierta estructura seleccionada, cuyas características pueden ser a priori idóneas para la inyección y almacenamiento del gas. Por este motivo la solución para poder analizar la viabilidad de un proyecto de almacenamiento en una estructura geológica pasa por la simulación numérica a partir de la modelización 3D del yacimiento. Los métodos numéricos permiten simular la inyección de un caudal determinado de dióxido de carbono desde un pozo de inyección localizado en una formación salina. En la presente tesis se ha definido una metodología de simulación de almacenamiento geológico de CO2, como contribución a la solución al problema del Cambio Climático, aplicada de forma concreta a la estructura BG-GE-08 (oeste de la Comunidad de Murcia). Esta estructura geológica ha sido catalogada por el IGME (Instituto Geológico y Minero de España) como idónea para el almacenamiento de dióxido de carbono, dada la existencia de una capa almacén confinada entre dos capas sello. ABSTRACT One of the climate change mitigation proposals suggested by the IPCC (Intergovermental Panel on Climate Change) in its ‘Synthesis Report 2007’ involves the launch of actions for capturing and storing carbon dioxide, existing three different geological structures suitable for gas storage: oil and gas reservoirs already drained, useless coal layers and deep saline structures. In case of deep saline structures, the main problem to develop a study of CO2 storage is the difficulty of obtaining geological data for some selected structure with characteristics that could be suitable for injection and gas storage. According to this situation, the solution to analyze the feasibility of a storage project in a geological structure will need numerical simulation from a 3D model. Numerical methods allow the simulation of the carbon dioxide filling in saline structures from a well, used to inject gas with a particular flow. In this document a simulation methodology has been defined for geological CO2 storage, as a contribution to solve the Climatic Change problem, applied to the structure BG-GE-08 (west of Murcia region). This geological structure has been classified by the IGME (Geological and Mining Institute of Spain) as suitable for the storage of carbon dioxide given the existence of a storage layer confined between two seal layers.

Simulation and geometrical design of multi-section tapered semiconductor optical amplifiers at 1.57 µm

Fully integrated semiconductor master-oscillator power-amplifiers (MOPA) with a tapered power amplifier are attractive sources for applications requiring high brightness. The geometrical design of the tapered amplifier is crucial to achieve the required power and beam quality. In this work we investigate by numerical simulation the role of the geometrical design in the beam quality and in the maximum achievable power. The simulations were performed with a Quasi-3D model which solves the complete steady-state semiconductor and thermal equations combined with a beam propagation method. The results indicate that large devices with wide taper angles produce higher power with better beam quality than smaller area designs, but at expenses of a higher injection current and lower conversion efficiency.

Nonlinear seismic behaviour of concrete arch bridges

Arch bridge structural solution has been known for centuries, in fact the simple nature of arch that require low tension and shear strength was an advantage as the simple materials like stone and brick were the only option back in ancient centuries. By the pass of time especially after industrial revolution, the new materials were adopted in construction of arch bridges to reach longer spans. Nowadays one long span arch bridge is made of steel, concrete or combination of these two as "CFST", as the result of using these high strength materials, very long spans can be achieved. The current record for longest arch belongs to Chaotianmen bridge over Yangtze river in China with 552 meters span made of steel and the longest reinforced concrete type is Wanxian bridge which also cross the Yangtze river through a 420 meters span. Today the designer is no longer limited by span length as long as arch bridge is the most applicable solution among other approaches, i.e. cable stayed and suspended bridges are more reasonable if very long span is desired. Like any super structure, the economical and architectural aspects in construction of a bridge is extremely important, in other words, as a narrower bridge has better appearance, it also require smaller volume of material which make the design more economical. Design of such bridge, beside the high strength materials, requires precise structural analysis approaches capable of integrating the combination of material behaviour and complex geometry of structure and various types of loads which may be applied to bridge during its service life. Depend on the design strategy, analysis may only evaluates the linear elastic behaviour of structure or consider the nonlinear properties as well. Although most of structures in the past were designed to act in their elastic range, the rapid increase in computational capacity allow us to consider different sources of nonlinearities in order to achieve a more realistic evaluations where the dynamic behaviour of bridge is important especially in seismic zones where large movements may occur or structure experience P - _ effect during the earthquake. The above mentioned type of analysis is computationally expensive and very time consuming. In recent years, several methods were proposed in order to resolve this problem. Discussion of recent developments on these methods and their application on long span concrete arch bridges is the main goal of this research. Accordingly available long span concrete arch bridges have been studied to gather the critical information about their geometrical aspects and properties of their materials. Based on concluded information, several concrete arch bridges were designed for further studies. The main span of these bridges range from 100 to 400 meters. The Structural analysis methods implemented in in this study are as following: Elastic Analysis: Direct Response History Analysis (DRHA): This method solves the direct equation of motion over time history of applied acceleration or imposed load in linear elastic range. Modal Response History Analysis (MRHA): Similar to DRHA, this method is also based on time history, but the equation of motion is simplified to single degree of freedom system and calculates the response of each mode independently. Performing this analysis require less time than DRHA. Modal Response Spectrum Analysis (MRSA): As it is obvious from its name, this method calculates the peak response of structure for each mode and combine them using modal combination rules based on the introduced spectra of ground motion. This method is expected to be fastest among Elastic analysis. Inelastic Analysis: Nonlinear Response History Analysis (NL-RHA): The most accurate strategy to address significant nonlinearities in structural dynamics is undoubtedly the nonlinear response history analysis which is similar to DRHA but extended to inelastic range by updating the stiffness matrix for every iteration. This onerous task, clearly increase the computational cost especially for unsymmetrical buildings that requires to be analyzed in a full 3D model for taking the torsional effects in to consideration. Modal Pushover Analysis (MPA): The Modal Pushover Analysis is basically the MRHA but extended to inelastic stage. After all, the MRHA cannot solve the system of dynamics because the resisting force fs(u; u_ ) is unknown for inelastic stage. The solution of MPA for this obstacle is using the previously recorded fs to evaluate system of dynamics. Extended Modal Pushover Analysis (EMPA): Expanded Modal pushover is a one of very recent proposed methods which evaluates response of structure under multi-directional excitation using the modal pushover analysis strategy. In one specific mode,the original pushover neglect the contribution of the directions different than characteristic one, this is reasonable in regular symmetric building but a structure with complex shape like long span arch bridges may go through strong modal coupling. This method intend to consider modal coupling while it take same time of computation as MPA. Coupled Nonlinear Static Pushover Analysis (CNSP): The EMPA includes the contribution of non-characteristic direction to the formal MPA procedure. However the static pushovers in EMPA are performed individually for every mode, accordingly the resulted values from different modes can be combined but this is only valid in elastic phase; as soon as any element in structure starts yielding the neutral axis of that section is no longer fixed for both response during the earthquake, meaning the longitudinal deflection unavoidably affect the transverse one or vice versa. To overcome this drawback, the CNSP suggests executing pushover analysis for governing modes of each direction at the same time. This strategy is estimated to be more accurate than MPA and EMPA, moreover the calculation time is reduced because only one pushover analysis is required. Regardless of the strategy, the accuracy of structural analysis is highly dependent on modelling and numerical integration approaches used in evaluation of each method. Therefore the widely used Finite Element Method is implemented in process of all analysis performed in this research. In order to address the study, chapter 2, starts with gathered information about constructed long span arch bridges, this chapter continuous with geometrical and material definition of new models. Chapter 3 provides the detailed information about structural analysis strategies; furthermore the step by step description of procedure of all methods is available in Appendix A. The document ends with the description of results and conclusion of chapter 4.

Construção tridimensional de atlas dos núcleos anteriores do tálamo a partir de cortes de um bloco de encéfalo

Tese de mestrado integrado em Engenharia Biomédica e Biofísica, apresentada à Universidade de Lisboa, através da Faculdade de Ciências, 2016

Functional characterisation of an engineered multidomain human P450 E1 by molecular Lego

The human cytochrome P450s constitute an important family of monooxygenase enzymes that carry out essential roles in the metabolism of endogenous compounds and foreign chemicals. We present here results of a fusion between a human P450 enzyme and a bacterial reductase that for the first time is shown does not require the addition of lipids or detergents to achieve wild-type-like activities. The fusion enzyme, P450 2E1-BMR, contains the N-terminally modified residues 22-493 of the human P450 2E1 fused at the C-terminus to residues 473-1049 of the P450 BM3 reductase (BMR). The P450 2E1-BMR enzyme is active, self-sufficient and presents the typical marker activities of the native human P450 2E1: the hydroxylation of p-nitrophenol (K (M)=1.84 +/- 0.09 mM and k (cat) of 2.98 +/- 0.04 nmol of p-nitrocatechol formed per minute per nanomole of P450) and chlorzoxazone (K (M)=0.65 +/- 0.08 mM and k (cat) of 0.95 +/- 0.10 nmol of 6-hydroxychlorzoxazone formed per minute per nanomole of P450). A 3D model of human P450 2E1 was generated to rationalise the functional data and to allow an analysis of the surface potentials. The distribution of charges on the model of P450 2E1 compared with that of the FMN domain of BMR provides the ground for the understanding of the interaction between the fused domains. The results point the way to successfully engineer a variety of catalytically self-sufficient human P450 enzymes for drug metabolism studies in solution.

Estimating parameters in stochastic systems:a variational Bayesian approach

This work is concerned with approximate inference in dynamical systems, from a variational Bayesian perspective. When modelling real world dynamical systems, stochastic differential equations appear as a natural choice, mainly because of their ability to model the noise of the system by adding a variation of some stochastic process to the deterministic dynamics. Hence, inference in such processes has drawn much attention. Here a new extended framework is derived that is based on a local polynomial approximation of a recently proposed variational Bayesian algorithm. The paper begins by showing that the new extension of this variational algorithm can be used for state estimation (smoothing) and converges to the original algorithm. However, the main focus is on estimating the (hyper-) parameters of these systems (i.e. drift parameters and diffusion coefficients). The new approach is validated on a range of different systems which vary in dimensionality and non-linearity. These are the Ornstein–Uhlenbeck process, the exact likelihood of which can be computed analytically, the univariate and highly non-linear, stochastic double well and the multivariate chaotic stochastic Lorenz ’63 (3D model). As a special case the algorithm is also applied to the 40 dimensional stochastic Lorenz ’96 system. In our investigation we compare this new approach with a variety of other well known methods, such as the hybrid Monte Carlo, dual unscented Kalman filter, full weak-constraint 4D-Var algorithm and analyse empirically their asymptotic behaviour as a function of observation density or length of time window increases. In particular we show that we are able to estimate parameters in both the drift (deterministic) and the diffusion (stochastic) part of the model evolution equations using our new methods.

ThreeDimensional Epigraphic Recording at Stobi (Former Yugoslav Republic of Macedonia): Creating a Virtual Lapidarium

There is an increased need for 3D recording of archaeological sites and digital preservation of their artifacts. Digital photogrammetry with prosumer DSLR cameras is a suitable tool for recording epigraphy in particular, as it allows for the recording of inscribed surfaces with very high accuracy, often better than 2 mm and with only a short time spent in the field. When photogrammetry is fused with other computational photography techniques like panoramic tours and Reflectance Transformation Imaging, a workflow exists to rival traditional LiDAR­based methods. The difficulty however, arises in the presentation of 3D data. It requires an enormous amount of storage and end­user sophistication. The proposed solution is to use game­engine technology and high definition virtual tours to provide not only scholars, but also the general public with an uncomplicated interface to interact with the detailed 3D epigraphic data. The site of Stobi, located near Gradsko, in the Former Yugoslav Republic of Macedonia (FYROM) was used as a case study to demonstrate the effectiveness of RTI, photogrammetry and virtual tour imaging working in combination. A selection of nine sets of inscriptions from the archaeological site were chosen to demonstrate the range of application for the techniques. The chosen marble, sandstone and breccia inscriptions are representative of the varying levels of deterioration and degradation of the epigraphy at Stobi, in which both their rates of decay and resulting legibility is varied. This selection includes those which are treated and untreated stones as well as those in situ and those in storage. The selection consists of both Latin and Greek inscriptions with content ranging from temple dedication inscriptions to statue dedications. This combination of 3D modeling techniques presents a cost and time efficient solution to both increase the legibility of severely damaged stones and to digitally preserve the current state of the inscriptions.

CAD Files and European Design Law

Three-dimensional printing (“3DP”) is an additive manufacturing technology that starts with a virtual 3D model of the object to be printed, the so-called Computer-Aided-Design (“CAD”) file. This file, when sent to the printer, gives instructions to the device on how to build the object layer-by-layer. This paper explores whether design protection is available under the current European regulatory framework for designs that are computer-created by means of CAD software, and, if so, under what circumstances. The key point is whether the appearance of a product, embedded in a CAD file, could be regarded as a protectable element under existing legislation. To this end, it begins with an inquiry into the concepts of “design” and “product”, set forth in Article 3 of the Community Design Regulation No. 6/2002 (“CDR”). Then, it considers the EUIPO’s practice of accepting 3D digital representations of designs. The enquiry goes on to illustrate the implications that the making of a CAD file available online might have. It suggests that the act of uploading a CAD file onto a 3D printing platform may be tantamount to a disclosure for the purposes of triggering unregistered design protection, and for appraising the state of the prior art. It also argues that, when measuring the individual character requirement, the notion of “informed user” and “the designer’s degree of freedom” may need to be reconsidered in the future. The following part touches on the exceptions to design protection, with a special focus on the repairs clause set forth in Article 110 CDR. The concluding part explores different measures that may be implemented to prohibit the unauthorised creation and sharing of CAD files embedding design-protected products.

La generazione di modelli tridimensionali densi in ambito architettonico mediante la fotogrammetria digitale

The Italian territory offers a wide range of treasures in the field of Cultural Assets. This is a highly relevant property, which needs an accurate management and preservation performed by appropriate tools, also giving attention to the maintenance and safeguard from risk factors. Nowadays the increasing development of new digital technologies, added by remarkable steps forward got by the subject of Geomatic makes possible an efficient integration among different techniques, helped also by spread of solutions to improve the data import-export and transmission between different devices. The main objective of this thesis is to experience the photogrammetric restitution implemented in a commercial software of digital photogrammetry, in order to generate a dense 3D model of the facade of the Basilica Sant'Apollinare Nuovo in Ravenna. The 1st Chapter, after a general introduction regarding the 3D survey of Cultural Heritage and some considerations linked to the use of digital photogrammetry in this field, is focused to analyze the case of stereoscopic and the monoscopic approach. In particular, it develops the theme of close-range photogrammetry. The 2nd Chapter, exposes the theme of digital images, from color theory until their appearing on the monitor. The 3rd Chapter, develops the case study of the Basilica di Sant'Apollinare Nuovo, the historical, architectural and religious of the same. Also, it is examined the issue of photogrammetry and laser scanning of the case study. The final part of the same chapter, treats the processing of data processing the software Agisoft PhotoScan, in order to generate, by means of Structure from Motion technique, a digital geometric 3D model of the Basilica Facade. The digital model has been scaled on the basis of measurements made on the field. With the software it was possible to accomplish the three phases of the photogrammetric data processing: internal orientation, exterior orientation and restitution.

Core Reactor Calculation Using the Adaptive Remeshing with a Current Error Estimator

With the objective to improve the reactor physics calculation on a 2D and 3D nuclear reactor via the Diffusion Equation, an adaptive automatic finite element remeshing method, based on the elementary area (2D) or volume (3D) constraints, has been developed. The adaptive remeshing technique, guided by a posteriori error estimator, makes use of two external mesh generator programs: Triangle and TetGen. The use of these free external finite element mesh generators and an adaptive remeshing technique based on the current field continuity show that they are powerful tools to improve the neutron flux distribution calculation and by consequence the power solution of the reactor core even though they have a minor influence on the critical coefficient of the calculated reactor core examples. Two numerical examples are presented: the 2D IAEA reactor core numerical benchmark and the 3D model of the Argonauta research reactor, built in Brasil.

Spatio-temporal coverage optimization of sensor networks

Les réseaux de capteurs sont formés d’un ensemble de dispositifs capables de prendre individuellement des mesures d’un environnement particulier et d’échanger de l’information afin d’obtenir une représentation de haut niveau sur les activités en cours dans la zone d’intérêt. Une telle détection distribuée, avec de nombreux appareils situés à proximité des phénomènes d’intérêt, est pertinente dans des domaines tels que la surveillance, l’agriculture, l’observation environnementale, la surveillance industrielle, etc. Nous proposons dans cette thèse plusieurs approches pour effectuer l’optimisation des opérations spatio-temporelles de ces dispositifs, en déterminant où les placer dans l’environnement et comment les contrôler au fil du temps afin de détecter les cibles mobiles d’intérêt. La première nouveauté consiste en un modèle de détection réaliste représentant la couverture d’un réseau de capteurs dans son environnement. Nous proposons pour cela un modèle 3D probabiliste de la capacité de détection d’un capteur sur ses abords. Ce modèle inègre également de l’information sur l’environnement grâce à l’évaluation de la visibilité selon le champ de vision. À partir de ce modèle de détection, l’optimisation spatiale est effectuée par la recherche du meilleur emplacement et l’orientation de chaque capteur du réseau. Pour ce faire, nous proposons un nouvel algorithme basé sur la descente du gradient qui a été favorablement comparée avec d’autres méthodes génériques d’optimisation «boites noires» sous l’aspect de la couverture du terrain, tout en étant plus efficace en terme de calculs. Une fois que les capteurs placés dans l’environnement, l’optimisation temporelle consiste à bien couvrir un groupe de cibles mobiles dans l’environnement. D’abord, on effectue la prédiction de la position future des cibles mobiles détectées par les capteurs. La prédiction se fait soit à l’aide de l’historique des autres cibles qui ont traversé le même environnement (prédiction à long terme), ou seulement en utilisant les déplacements précédents de la même cible (prédiction à court terme). Nous proposons de nouveaux algorithmes dans chaque catégorie qui performent mieux ou produits des résultats comparables par rapport aux méthodes existantes. Une fois que les futurs emplacements de cibles sont prédits, les paramètres des capteurs sont optimisés afin que les cibles soient correctement couvertes pendant un certain temps, selon les prédictions. À cet effet, nous proposons une méthode heuristique pour faire un contrôle de capteurs, qui se base sur les prévisions probabilistes de trajectoire des cibles et également sur la couverture probabiliste des capteurs des cibles. Et pour terminer, les méthodes d’optimisation spatiales et temporelles proposées ont été intégrées et appliquées avec succès, ce qui démontre une approche complète et efficace pour l’optimisation spatio-temporelle des réseaux de capteurs.

Study of a beam-column joint of a precast system by nonlinear numerical analysis

The present work consists of a detailed numerical analysis of a 4-way joint made of a precast column and two partially precast beams. The structure has been previously built and experimentally analyzed through a series of cyclic loads at the Laboratory of Tests on Structures (Laboratorio di Prove su Strutture, La. P. S.) of the University of Bologna. The aim of this work is to design a 3D model of the joint and then apply the techniques of nonlinear finite element analysis (FEA) to computationally reproduce the behavior of the structure under cyclic loads. Once the model has been calibrated to correctly emulate the joint, it is possible to obtain new insights useful to understand and explain the physical phenomena observed in the laboratory and to describe the properties of the structure, such as the cracking patterns, the force-displacement and the moment-curvature relations, as well as the deformations and displacements of the various elements composing the joint.

Exploring the World's Visual History

Thesis (Ph.D.)--University of Washington, 2016-08

Improving the Accuracy of 3-D Reconstruction in Robotic Vision Applications

Thesis (Ph.D.)--University of Washington, 2016-08