Table 1. Pressure, temperature and isotope ratios between H2O and CO2 for all experiments

Traditionally, the application of stable isotopes in Carbon Capture and Storage (CCS) projects has focused on d13C values of CO2 to trace the migration of injected CO2 in the subsurface. More recently the use of d18O values of both CO2 and reservoir fluids has been proposed as a method for quantifying in situ CO2 reservoir saturations due to O isotope exchange between CO2 and H2O and subsequent changes in d18OH2O values in the presence of high concentrations of CO2. To verify that O isotope exchange between CO2 and H2O reaches equilibrium within days, and that d18OH2O values indeed change predictably due to the presence of CO2, a laboratory study was conducted during which the isotope composition of H2O, CO2, and dissolved inorganic C (DIC) was determined at representative reservoir conditions (50°C and up to 19 MPa) and varying CO2 pressures. Conditions typical for the Pembina Cardium CO2 Monitoring Pilot in Alberta (Canada) were chosen for the experiments. Results obtained showed that d18O values of CO2 were on average 36.4±2.2 per mil (1 sigma, n=15) higher than those of water at all pressures up to and including reservoir pressure (19 MPa), in excellent agreement with the theoretically predicted isotope enrichment factor of 35.5 per mil for the experimental temperatures of 50°C. By using 18O enriched water for the experiments it was demonstrated that changes in the d18O values of water were predictably related to the fraction of O in the system sourced from CO2 in excellent agreement with theoretical predictions. Since the fraction of O sourced from CO2 is related to the total volumetric saturation of CO2 and water as a fraction of the total volume of the system, it is concluded that changes in d18O values of reservoir fluids can be used to calculate reservoir saturations of CO2 in CCS settings given that the d18O values of CO2 and water are sufficiently distinct.

A second order in time modified Lagrange-Galerkin finite element method for the incompressible Navier-Stokes equations

We introduce a second order in time modified Lagrange--Galerkin (MLG) method for the time dependent incompressible Navier--Stokes equations. The main ingredient of the new method is the scheme proposed to calculate in a more efficient manner the Galerkin projection of the functions transported along the characteristic curves of the transport operator. We present error estimates for velocity and pressure in the framework of mixed finite elements when either the mini-element or the $P2/P1$ Taylor--Hood element are used.

A method for incorporating residual stresses into patient-specific finite element simulations of arteries with an example on AAAs

Through progress in medical imaging, image analysis and finite element (FE) meshing tools it is now possible to extract patient-specific geometries from medical images of abdominal aortic aneurysms(AAAs), and thus to study clinically-relevant problems via FE simulations. Such simulations allow additional insight into human physiology in both healthy and diseased states. Medical imaging is most often performed in vivo, and hence the reconstructed model geometry in the problem of interest will represent the in vivo state, e.g., the AAA at physiological blood pressure. However, classical continuum mechanics and FE methods assume that constitutive models and the corresponding simulations begin from an unloaded, stress-free reference condition.

Automatic detection of crop rows in maize fields with high weeds pressure

This paper proposes a new method, oriented to crop row detection in images from maize fields with high weed pressure. The vision system is designed to be installed onboard a mobile agricultural vehicle, i.e. submitted to gyros, vibrations and undesired movements. The images are captured under image perspective, being affected by the above undesired effects. The image processing consists of three main processes: image segmentation, double thresholding, based on the Otsu’s method, and crop row detection. Image segmentation is based on the application of a vegetation index, the double thresholding achieves the separation between weeds and crops and the crop row detection applies least squares linear regression for line adjustment. Crop and weed separation becomes effective and the crop row detection can be favorably compared against the classical approach based on the Hough transform. Both gain effectiveness and accuracy thanks to the double thresholding that makes the main finding of the paper.

Effect of substrate-target distance and sputtering pressure in the synthesis of AlN thin films

In this work, we analyze the influence of the processing pressure and the substrate–target distance on the synthesis by reactive sputtering of c-axis oriented polycrystalline aluminum nitride thin films deposited on Si(100) wafers. The crystalline quality of AlN has been characterized by high-resolution X-ray diffraction (HR-XRD). The films exhibited a very high degree of c-axis orientation especially when a low process pressure was used. After growth, residual stress measurements obtained indirectly from radius of curvature measurements of the wafer prior and after deposition are also provided. Two different techniques are used to determine the curvature—an optically levered laser beam and a method based on X-ray diffraction. There is a transition from compressive to tensile stress at a processing pressure around 2 mTorr. The transition occurs at different pressures for thin films of different thickness. The degree of c-axis orientation was not affected by the target–substrate distance as it was varied in between 30 and 70 mm.

Modelling of landslides propagation with SPH : effects of rheology and pore water pressure

El estudio desarrollado en este trabajo de tesis se centra en la modelización numérica de la fase de propagación de los deslizamientos rápidos de ladera a través del método sin malla Smoothed Particle Hydrodynamics (SPH). Este método tiene la gran ventaja de permitir el análisis de problemas de grandes deformaciones evitando operaciones costosas de remallado como en el caso de métodos numéricos con mallas tal como el método de los Elementos Finitos. En esta tesis, particular atención viene dada al rol que la reología y la presión de poros desempeñan durante estos eventos. El modelo matemático utilizado se basa en la formulación de Biot-Zienkiewicz v - pw, que representa el comportamiento, expresado en términos de velocidad del esqueleto sólido y presiones de poros, de la mezcla de partículas sólidas en un medio saturado. Las ecuaciones que gobiernan el problema son: • la ecuación de balance de masa de la fase del fluido intersticial, • la ecuación de balance de momento de la fase del fluido intersticial y de la mezcla, • la ecuación constitutiva y • una ecuación cinemática. Debido a sus propiedades geométricas, los deslizamientos de ladera se caracterizan por tener una profundidad muy pequeña frente a su longitud y a su anchura, y, consecuentemente, el modelo matemático mencionado anteriormente se puede simplificar integrando en profundidad las ecuaciones, pasando de un modelo 3D a 2D, el cual presenta una combinación excelente de precisión, sencillez y costes computacionales. El modelo propuesto en este trabajo se diferencia de los modelos integrados en profundidad existentes por incorporar un ulterior modelo capaz de proveer información sobre la presión del fluido intersticial a cada paso computacional de la propagación del deslizamiento. En una manera muy eficaz, la evolución de los perfiles de la presión de poros está numéricamente resuelta a través de un esquema explicito de Diferencias Finitas a cada nodo SPH. Este nuevo enfoque es capaz de tener en cuenta la variación de presión de poros debida a cambios de altura, de consolidación vertical o de cambios en las tensiones totales. Con respecto al comportamiento constitutivo, uno de los problemas principales al modelizar numéricamente deslizamientos rápidos de ladera está en la dificultad de simular con la misma ley constitutiva o reológica la transición de la fase de iniciación, donde el material se comporta como un sólido, a la fase de propagación donde el material se comporta como un fluido. En este trabajo de tesis, se propone un nuevo modelo reológico basado en el modelo viscoplástico de Perzyna, pensando a la viscoplasticidad como a la llave para poder simular tanto la fase de iniciación como la de propagación con el mismo modelo constitutivo. Con el fin de validar el modelo matemático y numérico se reproducen tanto ejemplos de referencia con solución analítica como experimentos de laboratorio. Finalmente, el modelo se aplica a casos reales, con especial atención al caso del deslizamiento de 1966 en Aberfan, mostrando como los resultados obtenidos simulan con éxito estos tipos de riesgos naturales. The study developed in this thesis focuses on the modelling of landslides propagation with the Smoothed Particle Hydrodynamics (SPH) meshless method which has the great advantage of allowing to deal with large deformation problems by avoiding expensive remeshing operations as happens for mesh methods such as, for example, the Finite Element Method. In this thesis, special attention is given to the role played by rheology and pore water pressure during these natural hazards. The mathematical framework used is based on the v - pw Biot-Zienkiewicz formulation, which represents the behaviour, formulated in terms of soil skeleton velocity and pore water pressure, of the mixture of solid particles and pore water in a saturated media. The governing equations are: • the mass balance equation for the pore water phase, • the momentum balance equation for the pore water phase and the mixture, • the constitutive equation and • a kinematic equation. Landslides, due to their shape and geometrical properties, have small depths in comparison with their length or width, therefore, the mathematical model aforementioned can then be simplified by depth integrating the equations, switching from a 3D to a 2D model, which presents an excellent combination of accuracy, computational costs and simplicity. The proposed model differs from previous depth integrated models by including a sub-model able to provide information on pore water pressure profiles at each computational step of the landslide's propagation. In an effective way, the evolution of the pore water pressure profiles is numerically solved through a set of 1D Finite Differences explicit scheme at each SPH node. This new approach is able to take into account the variation of the pore water pressure due to changes of height, vertical consolidation or changes of total stress. Concerning the constitutive behaviour, one of the main issues when modelling fast landslides is the difficulty to simulate with the same constitutive or rheological model the transition from the triggering phase, where the landslide behaves like a solid, to the propagation phase, where the landslide behaves in a fluid-like manner. In this work thesis, a new rheological model is proposed, based on the Perzyna viscoplastic model, thinking of viscoplasticity as the key to close the gap between the triggering and the propagation phase. In order to validate the mathematical model and the numerical approach, benchmarks and laboratory experiments are reproduced and compared to analytical solutions when possible. Finally, applications to real cases are studied, with particular attention paid to the Aberfan flowslide of 1966, showing how the mathematical model accurately and successfully simulate these kind of natural hazards.

A novel method for measuring the dynamic fracture-initiation toughness under impulsive loadings

The design and development of a new method for performing fracture toughness tests under impulsive loadings using explosives is presented. The experimental set-up was complemented with pressure transducers and strain gauges in order to measure, respectively, the blast wave that reached the specimen and the loading history. Fracture toughness tests on a 7017-T73 aluminium alloy were carried out by using this device under impulsive loadings. Previous studies reported that such aluminium alloy had very little strain rate sensitivity, which made it an ideal candidate for comparison at different loading rates. The fracture-initiation toughness values of the 7017-T73 aluminium alloy obtained at impulsive loadings did not exhibit a significant variation from the cases studied at lower loading rates. Therefore, the method and device developed for measuring the dynamic fracture-initiation toughness under impulsive loadings was considered suitable for such a purpose

Pulsatile blood flow interaction with arterial walls of aorta : autoregulation and impedance pressure boundary condition and its biomedical applications

Para las decisiones urgentes sobre intervenciones quirúrgicas en el sistema cardiovascular se necesitan simulaciones computacionales con resultados fiables y que consuman un tiempo de cálculo razonable. Durante años los investigadores han trabajado en diversos métodos numéricos de cálculo que resulten atractivos para los cirujanos. Estos métodos, precisos pero costosos desde el punto de vista del coste computacional, crean un desajuste entre la oferta de los ingenieros que realizan las simulaciones y los médicos que operan en el quirófano. Por otra parte, los métodos de cálculo más simplificados reducen el tiempo de cálculo pero pueden proporcionar resultados no realistas. El objetivo de esta tesis es combinar los conceptos de autorregulación e impedancia del sistema circulatorio, la interacción flujo sanguíneo-pared arterial y modelos geométricos idealizados tridimensionales de las arterias pero sin pérdida de realismo, con objeto de proponer una metodología de simulación que proporcione resultados correctos y completos, con tiempos de cálculo moderados. En las simulaciones numéricas, las condiciones de contorno basadas en historias de presión presentan inconvenientes por ser difícil conocerlas con detalle, y porque los resultados son muy sensibles ante pequeñas variaciones de dichas historias. La metodología propuesta se basa en los conceptos de autorregulación, para imponer la demanda de flujo aguas abajo del modelo en el ciclo cardiaco, y la impedancia, para representar el efecto que ejerce el flujo en el resto del sistema circulatorio sobre las arterias modeladas. De este modo las historias de presión en el contorno son resultados del cálculo, que se obtienen de manera iterativa. El método propuesto se aplica en una geometría idealizada del arco aórtico sin patologías y en otra geometría correspondiente a una disección Stanford de tipo A, considerando la interacción del flujo pulsátil con las paredes arteriales. El efecto de los tejidos circundantes también se incorpora en los modelos. También se hacen aplicaciones considerando la interacción en una geometría especifica de un paciente anciano que proviene de una tomografía computarizada. Finalmente se analiza una disección Stanford tipo B con tres modelos que incluyen la fenestración del saco. Clinicians demand fast and reliable numerical results of cardiovascular biomechanic simulations for their urgent pre-surgery decissions. Researchers during many years have work on different numerical methods in order to attract the clinicians' confidence to their colorful contours. Though precise but expensive and time-consuming methodologies create a gap between numerical biomechanics and hospital personnel. On the other hand, simulation simplifications with the aim of reduction in computational time may cause in production of unrealistic outcomes. The main objective of the current investigation is to combine ideas such as autoregulation, impedance, fluid-solid interaction and idealized geometries in order to propose a computationally cheap methodology without excessive or unrealistic simplifications. The pressure boundary conditions are critical and polemic in numerical simulations of cardiovascular system, in which a specific arterial site is of interest and the rest of the netwrok is neglected but represented by a boundary condition. The proposed methodology is a pressure boundary condition which takes advantage of numerical simplicity of application of an imposed pressure boundary condition on outlets, while it includes more sophisticated concepts such as autoregulation and impedance to gain more realistic results. Incorporation of autoregulation and impedance converts the pressure boundary conditions to an active and dynamic boundary conditions, receiving feedback from the results during the numerical calculations and comparing them with the physiological requirements. On the other hand, the impedance boundary condition defines the shapes of the pressure history curves applied at outlets. The applications of the proposed method are seen on idealized geometry of the healthy arotic arch as well as idealized Stanford type A dissection, considering the interaction of the arterial walls with the pulsatile blood flow. The effect of surrounding tissues is incorporated and studied in the models. The simulations continue with FSI analysis of a patient-specific CT scanned geometry of an old individual. Finally, inspiring of the statistic results of mortality rates in Stanford type B dissection, three models of fenestrated dissection sac is studied and discussed. Applying the developed boundary condition, an alternative hypothesis is proposed by the author with respect to the decrease in mortality rates in patients with fenestrations.

Use of intrinsic modes in biology: Examples of indicial response of pulmonary blood pressure to ± step hypoxia

Recently, a new method to analyze biological nonstationary stochastic variables has been presented. The method is especially suitable to analyze the variation of one biological variable with respect to changes of another variable. Here, it is illustrated by the change of the pulmonary blood pressure in response to a step change of oxygen concentration in the gas that an animal breathes. The pressure signal is resolved into the sum of a set of oscillatory intrinsic mode functions, which have zero “local mean,” and a final nonoscillatory mode. With this device, we obtain a set of “mean trends,” each of which represents a “mean” in a definitive sense, and together they represent the mean trend systematically with different degrees of oscillatory content. Correspondingly, the oscillatory content of the signal about any mean trend can be represented by a set of partial sums of intrinsic mode functions. When the concept of “indicial response function” is used to describe the change of one variable in response to a step change of another variable, we now have a set of indicial response functions of the mean trends and another set of indicial response functions to describe the energy or intensity of oscillations about each mean trend. Each of these can be represented by an analytic function whose coefficients can be determined by a least-squares curve-fitting procedure. In this way, experimental results are stated sharply by analytic functions.

Pressure-mediated oligonucleotide transfection of rat and human cardiovascular tissues

The application of gene therapy to human disease is currently restricted by the relatively low efficiency and potential hazards of methods of oligonucleotide or gene delivery. Antisense or transcription factor decoy oligonucleotides have been shown to be effective at altering gene expression in cell culture expreriments, but their in vivo application is limited by the efficiency of cellular delivery, the intracellular stability of the compounds, and their duration of activity. We report herein the development of a highly efficient method for naked oligodeoxynucleotide (ODN) transfection into cardiovascular tissues by using controlled, nondistending pressure without the use of viral vectors, lipid formulations, or exposure to other adjunctive, potentially hazardous substances. In this study, we have documented the ability of ex vivo, pressure-mediated transfection to achieve nuclear localization of fluorescent (FITC)-labeled ODN in approximately 90% and 50% of cells in intact human saphenous vein and rat myocardium, respectively. We have further documented that pressure-mediated delivery of antisense ODN can functionally inhibit target gene expression in both of these tissues in a sequence-specific manner at the mRNA and protein levels. This oligonucleotide transfection system may represent a safe means of achieving the intraoperative genetic engineering of failure-resistant human bypass grafts and may provide an avenue for the genetic manipultation of cardiac allograft rejection, allograft vasculopathy, or other transplant diseases.

3D numerical modelling of acoustic horns using the method of fundamental solutions

In the present work, a three-dimensional (3D) formulation based on the method of fundamental solutions (MFS) is applied to the study of acoustic horns. The implemented model follows and extends previous works that only considered two-dimensional and axisymmetric horn configurations. The more realistic case of 3D acoustic horns with symmetry regarding two orthogonal planes is addressed. The use of the domain decomposition technique with two interconnected sub-regions along a continuity boundary is proposed, allowing for the computation of the sound pressure generated by an acoustic horn installed on a rigid screen. In order to reduce the model discretization requirements for these cases, Green’s functions derived with the image source methodology are adopted, automatically accounting for the presence of symmetry conditions. A strategy for the calculation of an optimal position of the virtual sources used by the MFS to define the solution is also used, leading to improved reliability and flexibility of the proposed method. The responses obtained by the developed model are compared to reference solutions, computed by well-established models based on the boundary element method. Additionally, numerically calculated acoustic parameters, such as directivity and beamwidth, are compared with those evaluated experimentally.

Refined method of potential enhancement in the equilibria characterization of activated carbon. Comparison with GCMC and DFT

In this paper, we present a technique for equilibria characterization of activated carbon having slit-shaped pores. This method was first developed by Do (Do, D. D. A new method for the characterisation of micro-mesoporous materials. Presented at the International Symposium on New Trends in Colloid and Interface Science, September 24-26, 1998 Chiba, Japan) and applied by his group and other groups for characterization of pore size distribution (PSD) as well as adsorption equilibria determination of a wide range of hydrocarbons. It is refined in this paper and compared with the grand canonical Monte Carlo (GCMG) simulation and density functional theory (DFT). The refined theory results in a good agreement between the pore filling pressure versus pore width and those obtained by GCMG and DFT. Furthermore, our local isotherms are qualitatively in good agreement with those obtained by the GCMC simulations. The main advantage of this method is that it is about 4 orders of magnitude faster than the GCMC simulations, making it suitable for optimization studies and design purposes. Finally, we apply our method and the GCMG in the derivation of the PSD of a commercial activated carbon. It was found that the PSD derived from our method is comparable to that derived from the GCMG simulations.

Corrosion behaviour of a pressure die cast magnesium alloy

High pressure die casting is the most important production method for casting magnesium alloy components, and uniformity of appearance is an important criterion for acceptance of a component by customers. This paper investigates the influence of uniformity in surface appearance of diecast AZ91D plates on their corrosion behaviour. Through immersion, hydrogen collection and weight loss measurements it was found that corrosion is more likely to occur on the areas of the plate that appear to be darker, leading to a non-uniformly corroded surface. Microstructural analysis showed that the non-uniformity in appearance is related to a difference in the morphology and distribution of porosity across the surface of a diecast AZ91D plate. The darker areas of the surface are high in porosity which breaks the continuity of the beta-phase network and provides shortcut paths for corrosion from the surface to the interior of the casting. The brighter shiny areas of the surface are much less porous, with isolated pores being confined by corrosion resistant beta-precipitates thus reducing the corrosion rate.