Aplicación del método de los volúmenes finitos al cálculo del parámetro de masa añadida de un cilindro vibrante en el seno de un fluido

Máster Universitario en Sistemas Inteligentes y Aplicaciones Numéricas en Ingeniería (SIANI)

Fluidodinamica biomedica computazionale del flusso sanguigno in presenza di affezioni strumentali e/o patologiche

La tesi di Dottorato studia il flusso sanguigno tramite un codice agli elementi finiti (COMSOL Multiphysics). Nell’arteria è presente un catetere Doppler (in posizione concentrica o decentrata rispetto all’asse di simmetria) o di stenosi di varia forma ed estensione. Le arterie sono solidi cilindrici rigidi, elastici o iperelastici. Le arterie hanno diametri di 6 mm, 5 mm, 4 mm e 2 mm. Il flusso ematico è in regime laminare stazionario e transitorio, ed il sangue è un fluido non-Newtoniano di Casson, modificato secondo la formulazione di Gonzales & Moraga. Le analisi numeriche sono realizzate in domini tridimensionali e bidimensionali, in quest’ultimo caso analizzando l’interazione fluido-strutturale. Nei casi tridimensionali, le arterie (simulazioni fluidodinamiche) sono infinitamente rigide: ricavato il campo di pressione si procede quindi all’analisi strutturale, per determinare le variazioni di sezione e la permanenza del disturbo sul flusso. La portata sanguigna è determinata nei casi tridimensionali con catetere individuando tre valori (massimo, minimo e medio); mentre per i casi 2D e tridimensionali con arterie stenotiche la legge di pressione riproduce l’impulso ematico. La mesh è triangolare (2D) o tetraedrica (3D), infittita alla parete ed a valle dell’ostacolo, per catturare le ricircolazioni. Alla tesi sono allegate due appendici, che studiano con codici CFD la trasmissione del calore in microcanali e l’ evaporazione di gocce d’acqua in sistemi non confinati. La fluidodinamica nei microcanali è analoga all’emodinamica nei capillari. Il metodo Euleriano-Lagrangiano (simulazioni dell’evaporazione) schematizza la natura mista del sangue. La parte inerente ai microcanali analizza il transitorio a seguito dell’applicazione di un flusso termico variabile nel tempo, variando velocità in ingresso e dimensioni del microcanale. L’indagine sull’evaporazione di gocce è un’analisi parametrica in 3D, che esamina il peso del singolo parametro (temperatura esterna, diametro iniziale, umidità relativa, velocità iniziale, coefficiente di diffusione) per individuare quello che influenza maggiormente il fenomeno.

Carbonate conical mounds of the eastern Anti-Atlas(Morocco) and their relationship with potential analogues on Mars

Longstanding debates concerning the origin of the Kess Kess Emsian carbonate mounds exposed at Hamar Laghdad Ridge (eastern Anti-Atlas, Morocco) centre around the processes that induced precipitation of carbonate mud and the preservation of steep morphologies. Although in the last years an origin related to hydrothermalism seemed to be more likely, to date the Kess Kess are still considered controversial vent deposits. This study combines in updated research review information from previous work and new detailed field observations coupled with new analytical results to define a consistent framework and some new insight of current knowledge about the origin of these mounds. We obtain a complete minero-petrographic and palaeobiological data set and a detailed geochemical characterization of the different lithologies and facies of the Hamar Laghdad stratigraphic succession, including mounds, and we compared the results with the data from Maïder Basin mounds (Anti-Atlas, Morocco). Our data support the hydrothermal model proposed for the genesis and development of the Kess Kess mounds. The mechanisms linked to the mounds formation and growth are discussed in the light of the new finding of fluid-sediment interaction within a scenario driven by late magmatic fluids circulation. Conical mounds and other fluids related morphologies were also reported from Crommelin crater area (Arabia Terra, Mars). These mounds consist in meter-sized conical buildups hosted in the Equatorial Layered Deposits (ELDs) deposed during a regional groundwater fluid upwelling. Geometries and geological conditions that might have controlled the development of such morphologies were discussed. According to our data the morphological and stratigraphical characteristics of Crommelin area mounds are most consistent with a formation by fluids advection. Then we compare terrestrial and Martian data and examine the geological settings of hydrothermal mound occurrences on Earth in order to describe potential target areas for hydrothermal structures on Mars.

The simultaneous role of an alveolus as flow mixer and flow feeder for the deposition of inhaled submicron particles

In an effort to understand the fate of inhaled submicron particles in the small sacs, or alveoli, comprising the gas-exchange region of the lung, we calculated the flow in three-dimensional (3D) rhythmically expanding models of alveolated ducts. Since convection toward the alveolar walls is a precursor to particle deposition, it was the goal of this paper to investigate the streamline maps' dependence upon alveoli location along the acinar tree. On the alveolar midplane, the recirculating flow pattern exhibited closed streamlines with a stagnation saddle point. Off the midplane we found no closed streamlines but nested, funnel-like, spiral, structures (reminiscent of Russian nesting dolls) that were directed towards the expanding walls in inspiration, and away from the contracting walls in expiration. These nested, funnel-like, structures were surrounded by air that flowed into the cavity from the central channel over inspiration and flowed from the cavity to the central channel over expiration. We also found that fluid particle tracks exhibited similar nested funnel-like spiral structures. We conclude that these unique alveolar flow structures may be of importance in enhancing deposition. In addition, due to inertia, the nested, funnel-like, structures change shape and position slightly during a breathing cycle, resulting in flow mixing. Also, each inspiration feeds a fresh supply of particle-laden air from the central channel to the region surrounding the mixing region. Thus, this combination of flow mixer and flow feeder makes each individual alveolus an effective mixing unit, which is likely to play an important role in determining the overall efficiency of convective mixing in the acinus.

Effect of mesh distortion on the accuracy of high order vorticity-velocity CFD approaches

The numerical solution of the incompressible Navier-Stokes equations offers an alternative to experimental analysis of fluid-structure interaction (FSI). We would save a lot of time and effort and help cut back on costs, if we are able to accurately model systems by these numerical solutions. These advantages are even more obvious when considering huge structures like bridges, high rise buildings or even wind turbine blades with diameters as large as 200 meters. The modeling of such processes, however, involves complex multiphysics problems along with complex geometries. This thesis focuses on a novel vorticity-velocity formulation called the Kinematic Laplacian Equation (KLE) to solve the incompressible Navier-stokes equations for such FSI problems. This scheme allows for the implementation of robust adaptive ordinary differential equations (ODE) time integration schemes, allowing us to tackle each problem as a separate module. The current algortihm for the KLE uses an unstructured quadrilateral mesh, formed by dividing each triangle of an unstructured triangular mesh into three quadrilaterals for spatial discretization. This research deals with determining a suitable measure of mesh quality based on the physics of the problems being tackled. This is followed by exploring methods to improve the quality of quadrilateral elements obtained from the triangles and thereby improving the overall mesh quality. A series of numerical experiments were designed and conducted for this purpose and the results obtained were tested on different geometries with varying degrees of mesh density.

Effect of high order interpolation in the stability and efficiency of the time-integration process in vorticity-velocity CFD algorithms

The numerical solution of the incompressible Navier-Stokes Equations offers an effective alternative to the experimental analysis of Fluid-Structure interaction i.e. dynamical coupling between a fluid and a solid which otherwise is very complex, time consuming and very expensive. To have a method which can accurately model these types of mechanical systems by numerical solutions becomes a great option, since these advantages are even more obvious when considering huge structures like bridges, high rise buildings, or even wind turbine blades with diameters as large as 200 meters. The modeling of such processes, however, involves complex multiphysics problems along with complex geometries. This thesis focuses on a novel vorticity-velocity formulation called the KLE to solve the incompressible Navier-stokes equations for such FSI problems. This scheme allows for the implementation of robust adaptive ODE time integration schemes and thus allows us to tackle the various multiphysics problems as separate modules. The current algorithm for KLE employs a structured or unstructured mesh for spatial discretization and it allows the use of a self-adaptive or fixed time step ODE solver while dealing with unsteady problems. This research deals with the analysis of the effects of the Courant-Friedrichs-Lewy (CFL) condition for KLE when applied to unsteady Stoke’s problem. The objective is to conduct a numerical analysis for stability and, hence, for convergence. Our results confirmthat the time step ∆t is constrained by the CFL-like condition ∆t ≤ const. hα, where h denotes the variable that represents spatial discretization.

Dual-Phase Liquid Xenon Detectors for Dark Matter Searches

Dual-phase time projection chambers (TPCs) filled with the liquid noble gas xenon (LXe) are currently the most sensitive detectors searching for interactions of WIMP dark matter in a laboratory-based experiment. This is achieved by combining a large, monolithic dark matter target of a very low background with the capability to localize the interaction vertex in three dimensions, allowing for target fiducialization and multiple-scatter rejection. The background in dual-phase LXe TPCs is further reduced by the simultaneous measurement of the scintillation and ionization signal from a particle interaction, which is used to distinguish signal from background signatures. This article reviews the principle of dual-phase LXe TPCs, and provides an overview about running as well as future experimental efforts.

FEM/FD Immersed Boundary FSI Simulations

Immersed boundary simulations have been under development for physiological flows, allowing for elegant handling of fluid-structure interaction modelling with large deformations due to retained domain-specific meshing. We couple a structural system in Lagrangian representation that is formulated in a weak form with a Navier-Stokes system discretized through a finite differences scheme. We build upon a proven highly scalable imcompressible flow solver that we extend to handle FSI. We aim at applying our method to investigating the hemodynamics of Aortic Valves. The code is going to be extended to conform to the new hybrid-node supercomputers.

Microscale mapping of alteration conditions and potential biosignatures in basaltic-ultramafic rocks on early Earth and beyond

Subseafloor environments preserved in Archean greenstone belts provide an analogue for investigating potential subsurface habitats on Mars. The c. 3.5-3.4 Ga pillow lava metabasalts of the mid-Archean Barberton greenstone belt, South Africa, have been argued to contain the earliest evidence for microbial subseafloor life. This includes candidate trace fossils in the form of titanite microtextures, and sulfur isotopic signatures of pyrite preserved in metabasaltic glass of the c. 3.472 Ga Hooggenoeg Formation. It has been contended that similar microtextures in altered martian basalts may represent potential extraterrestrial biosignatures of microbe-fluid-rock interaction. But despite numerous studies describing these putative early traces of life, a detailed metamorphic characterization of the microtextures and their host alteration conditions in the ancient pillow lava metabasites is lacking. Here, we present a new nondestructive technique with which to study the in situ metamorphic alteration conditions associated with potential biosignatures in mafic-ultramafic rocks of the Hooggenoeg Formation. Our approach combines quantitative microscale compositional mapping by electron microprobe with inverse thermodynamic modeling to derive low-temperature chlorite crystallization conditions. We found that the titanite microtextures formed under subgreenschist to greenschist facies conditions. Two chlorite temperature groups were identified in the maps surrounding the titanite microtextures and record peak metamorphic conditions at 315 ± 40°C (XFe3+(chlorite) = 25-34%) and lower-temperature chlorite veins/microdomains at T = 210 ± 40°C (lower XFe3+(chlorite) = 40-45%). These results provide the first metamorphic constraints in textural context on the Barberton titanite microtextures and thereby improve our understanding of the local preservation conditions of these potential biosignatures. We suggest that this approach may prove to be an important tool in future studies to assess the biogenicity of these earliest candidate traces of life on Earth. Furthermore, we propose that this mapping approach could also be used to investigate altered mafic-ultramafic extraterrestrial samples containing candidate biosignatures.

Igneous geochemistry of OP Leg 126 samples

Major element, trace element, and radiogenic isotope compositions of samples collected from Ocean Drilling Program Leg 126 in the Izu-Bonin forearc basin are presented. Lavas from the center of the basin (Site 793) are high-MgO, low-Ti, two-pyroxene basaltic andesites, and represent the products of synrift volcanism in the forearc region. These synrift lavas share many of the geochemical and petrographic characteristics of boninites. In terms of their element abundances, ratios, and isotope systematics they are intermediate between low-Ti arc tholeiites from the active arc and boninites of the outer-arc high. These features suggest a systematic geochemical gradation between volcanics related to trench distance and a variably depleted source. A basement high drilled on the western flank of the basin (Site 792) comprises a series of plagioclase-rich two-pyroxene andesites with calc-alkaline affinities. These lavas are similar to calc-alkaline volcanics from Japan, but have lower contents of Ti, Zr, and low-field-strength elements (LFSE). Lavas from Site 793 show inter-element variations between Zr, Ti, Sr, Ni, and Cr that are consistent with those predicted during crystallization and melting processes. In comparison, concentrations of P, Y, LFSE, and the rare-earth elements (REE) are anomalous. These elements have been redistributed within the lava pile, concentrating particularly in sections of massive and pillowed flows. Relative movement of these two-element groupings can be related to the alteration of interstitial basaltic andesite glass to a clay mineral assemblage by a post-eruptive process. Fluid-rock interaction has produced similar effects in the basement lavas of Site 792. In this sequence, andesites and dacites have undergone a volume change related to silica mobility. As a result of this process, some lithologies have the major element characteristics of basaltic andesite and rhyolite, but can be related to andesitic or dacitic precursors by silica removal or addition.

Chemical and isotopic compositions of rocks and minerals from the Palenyi Island, Por'ya Guba (Bay), White Sea

Hypersthene-garnet-sillimanite-quartz enclaves were studied in orthopyroxene-plagioclase and orthopyroxene-clinopyroxene crystalline schists and gneisses from shear zones exposed in the Palenyi Island within the Early Proterozoic Belomorian Mobile Belt. Qualitative analysis of mineral assemblages indicates that these rocks were metamorphosed to the granulite facies (approximately 900°C and 10-11 kbar). Oxygen isotopic composition was determined in rock-forming minerals composing zones of the enclaves of various mineral and chemical composition. Closure temperatures of the isotopic systems obtained by methods of oxygen isotopic thermometry are close to values obtained with mineralogical geothermometers (garnet-orthopyroxene and garnet-biotite) and correspond to the high-temperature granulite facies (860-900°C). Identified systematic variations in d18O values were determined in the same minerals from zones of different mineral composition. Inasmuch as these zones are practically in contact with one another, these variations in d18O cannot be explained by primary isotopic heterogeneity of the protolith. Model calculations of the extent and trend of d18O variations in minerals suggest that fluid-rock interaction at various integral fluid/rock ratios in discrete zones was the only mechanism that could generate the zoning. This demonstrates that focused fluid flux could occur in lower crustal shear zones. Preservation of high-temperature isotopic equilibria of minerals testifies that the episode of fluid activity at the peak of metamorphism was very brief.

Geochemistry of lavas and dikes from the upper oceanic crust of Hess Deep Rift, eastern Pacific Ocean

Strontium isotopes are useful tracers of fluid-rock interaction in marine hydrothermal systems and provide a potential way to quantify the amount of seawater that passes through these systems. We have determined the whole-rock Sr-isotopic compositions of a section of upper oceanic crust that formed at the fast-spreading East Pacific Rise, now exposed at Hess Deep. This dataset provides the first detailed comparison for the much-studied Ocean Drilling Program (ODP) drill core from Site 504B. Whole-rock and mineral Sr concentrations indicate that Sr-exchange between hydrothermal fluids and the oceanic crust is complex, being dependent on the mineralogical reactions occurring; in particular, epidote formation takes up Sr from the fluid increasing the 87Sr/86Sr of the bulk-rock. Calculating the fluid-flux required to shift the Sr-isotopic composition of the Hess Deep sheeted-dike complex, using the approach of Bickle and Teagle (1992, doi:10.1016/0012-821X(92)90221-G) gives a fluid-flux similar to that determined for ODP Hole 504B. This suggests that the level of isotopic exchange observed in these two regions is probably typical for modern oceanic crust. Unfortunately, uncertainties in the modeling approach do not allow us to determine a fluid-flux that is directly comparable to fluxes calculated by other methods.

Table 1. Mean geochemical values for Baseline (pre-injection) and Monitor-1 (post-injection) samples, and the percent change between the two weighted means

The Weyburn Oil Field, Saskatchewan is the site of a large (5000 tonnes/day of CO2) CO2-EOR injection project By EnCana Corporation. Pre- and post-injection samples (Baseline and Monitor-1, respectively) of produced fluids from approximately 45 vertical wells were taken and chemically analyzed to determine changes in the fluid chemistry and isotope composition between August 2000 and March 2001. After 6 months of CO2 injection, geochemical parameters including pH, [HCO3], [Ca], [Mg], and ?13CO2(g) point to areas in which injected CO2 dissolution and reservoir carbonate mineral dissolution have occurred. Pre-injection fluid compositions suggest that the reservoir brine in the injection area may be capable of storing as much as 100 million tonnes of dissolved CO2. Modeling of water-rock reactions show that clay minerals and feldspar, although volumetrically insignificant, may be capable of acting as pH buffers, allowing injected CO2 to be stored as bicarbonate in the formation water or as newly precipitated carbonate minerals, given favorable reaction kinetics.

Chemistry of phyllosilicates of ODP Hole 140-504B

Replacement minerals in olivine record the evolution of hydrothermal alteration between 1600 and 2000 mbsf in the sheeted dike complex in Hole 504B. 1. Talc (+ magnetite) rim on olivine represents the earliest alteration. Talc probably crystallized during initial cooling of the dikes. 2. The partial breakdown of talc to "deweylite", a chaotic mixture of serpentine and Al-free stevensite, was facilitated by further cooling and a somewhat increased fluid:rock interaction in the dikes. 3. The presence of chlorite veins and the replacement of unaltered olivine cores, talc, and deweylite and of other silicates by chlorite suggest fracturing of the rocks during cooling (shrinkage cracks) and local influx of seawater into the dikes. 4. Late amphibole veins and locally extensive amphibole alteration indicate increasing temperature and the development of new sets of fractures, possibly due to the injection of fresh magma. Several generations of chlorite and amphibole veins are present in the dikes. Offset veins and the crack-seal texture within veins in the dikes suggest that the alteration cycle was probably repeated with the injection of each set of new dikes. Presently measured temperatures (195°C) at 2000 m depth in Hole 504B indicate that deweylite, which was previously considered a low-temperature mineral, can form well above its previously estimated crystallization temperature of 50°C.