Theory manual to OctVCE – a cartesian cell CFD code with special application to blast wave problems, Department of Mechanical Engineering Report 2007/12

OctVCE is a cartesian cell CFD code produced especially for numerical simulations of shock and blast wave interactions with complex geometries, in particular, from explosions. Virtual Cell Embedding (VCE) was chosen as its cartesian cell kernel for its simplicity and sufficiency for practical engineering design problems. The code uses a finite-volume formulation of the unsteady Euler equations with a second order explicit Runge-Kutta Godonov (MUSCL) scheme. Gradients are calculated using a least-squares method with a minmod limiter. Flux solvers used are AUSM, AUSMDV and EFM. No fluid-structure coupling or chemical reactions are allowed, but gas models can be perfect gas and JWL or JWLB for the explosive products. This report also describes the code’s ‘octree’ mesh adaptive capability and point-inclusion query procedures for the VCE geometry engine. Finally, some space will also be devoted to describing code parallelization using the shared-memory OpenMP paradigm. The user manual to the code is to be found in the companion report 2007/13.

Finite element modelling of rock alteration and metamorphic process in hydrothermal systems

We use the finite element method to simulate the rock alteration and metamorphic process in hydrothermal systems. In particular, we consider the fluid-rock interaction problems in pore-fluid saturated porous rocks. Since the fluid rock interaction takes place at the contact interface between the pore-fluid and solid minerals, it is governed by the chemical reaction which usually takes place very slowly at this contact interface, from the geochemical point of view. Due to the relative slowness of the rate of the chemical reaction to the velocity of the pore-fluid flow in the hydrothermal system to be considered, there exists a retardation zone, in which the conventional static theory in geochemistry does not hold true. Since this issue is often overlooked by some purely numerical modellers, it is emphasized in this paper. The related results from a typical rock alteration and metamorphic problem in a hydrothermal system have shown not only the detailed rock alteration and metamorphic process, but also the size of the retardation zone in the hydrothermal system. Copyright (C) 2001 John Wiley & Sons, Ltd.

A modified pore-filling isotherm for liquid-phase adsorption in activated carbon

This article modifies the usual form of the Dubinin-Radushkevich pore-filling model for application to liquid-phase adsorption data, where large molecules are often involved. In such cases it is necessary to include the repulsive part of the energy in the micropores, which is accomplished here by relating the pore potential to the fluid-solid interaction potential. The model also considers the nonideality of the bulk liquid phase through the UNIFAC activity coefficient model, as well as structural heterogeneity of the carbon. For the latter the generalized adsorption integral is used while incorporating the pore-size distribution obtained by density functional theory analysis of argon adsorption data. The model is applied here to the interpretation of aqueous phase adsorption isotherms of three different esters on three commercial activated carbons. Excellent agreement between the model and experimental data is observed, and the fitted Lennard-Jones size parameter for the adsorbate-adsorbate interactions compares well with that estimated from known critical properties, supporting the modified approach. On the other hand, the model without consideration of bulk nonideality, or when using classical models of the characteristic energy, gives much poorer bts of the data and unrealistic parameter values.

Multicomponent supercritical adsorption in microporous activated carbon materials

In this paper, a new technique for predicting multicomponent adsorption equilibria of supercritical fluids in microporous carbons is presented. In difference from adsorption on a surface, which is a function of the fluid-solid interaction only, adsorption in porous media is influenced by the proximity of the pore walls, resulting in the enhancement in adsorption affinity. The degree of this enhancement is different for different adsorbates, and it increases with a decrease in pore size. The theory is applied to a number of carbonaceous systems with good success.

Determination of stresses around a cylindrical single pile caused by horizontal movements of soils with the finitive element method

This study was carried out with the aim of modeling in 2D, in plain strain, the movement of a soft cohesive soil around a pile, in order to enable the determination of stresses resulting along the pile, per unit length. The problem in study fits into the large deformations problem and can be due to landslide, be close of depth excavations, to be near of zones where big loads are applied in the soil, etc. In this study is used an constitutive Elasto-Plastic model with the failure criterion of Mohr-Coulomb to model the soil behavior. The analysis is developed considering the soil in undrained conditions. To the modeling is used the finite element program PLAXIS, which use the Updated Lagrangian - Finite Element Method (UL-FEM). In this work, special attention is given to the soil-pile interaction, where is presented with some detail the formulation of the interface elements and some studies for a better understand of his behavior. It is developed a 2-D model that simulates the effect of depth allowing the study of his influence in the stress distribution around the pile. The results obtained give an important base about how behaves the movement of the soil around a pile, about how work the finite element program PLAXIS and how is the stress distribution around the pile. The analysis demonstrate that the soil-structure interaction modeled with the UL-FEM and interface elements is more appropriate to small deformations problems.

Pontes integrais: estado da arte

Dissertação para obtenção do grau de mestre em Engenharia Civil

Avaliação do comportamento de uma ponte metálica centenária sujeita à ação de tráfego ferroviário ligeiro

A presente dissertação tem como objetivo principal a análise numérica do comportamento dinâmico de uma ponte ferroviária, sob ação de tráfego ligeiro ferroviário. Neste contexto são apresentados alguns fundamentos teóricos a ter em conta nestes domínios, visando uma melhor compreensão dos fenómenos existentes no comportamento dinâmico de pontes ferroviárias quando sujeitas ao tráfego. O caso de estudo teve como foco a ponte Luiz I, uma ponte metálica situada sobre o rio Douro, que liga as cidades do Porto e Vila Nova de Gaia, sob ação de tráfego ligeiro ferroviário no seu tabuleiro superior para a condição anterior aos trabalhos de reabilitação e reforço realizados entre 2004 e 2005. Para o efeito foi desenvolvido um modelo numérico de elementos finitos da ponte realizado com recurso ao programa ANSYS, assim como um modelo numérico do veículo do Metro de Lisboa. Com base nestes modelos foram obtidos os parâmetros modais, nomeadamente as frequências naturais e os modos de vibração de toda a estrutura e do veículo. O estudo do comportamento dinâmico da ponte foi realizado por intermédio de uma metodologia de cargas móveis e de interação veículo-estrutura, através da ferramenta computacional Train-Bridge Interaction (TBI). As análises dinâmicas foram efetuadas para a passagem dos veículos de passageiros das redes de Metros do Porto e Lisboa. Nestas análises é estudada a resposta da estrutura em função da variabilidade ao nível da secção transversal, dependência do tramo, influência do veículo, da sua velocidade de circulação e impacto das frequências de vibração estimadas pelo modelo numérico.

Dynamic stiffness of the piled foundation of a onshore wind turbine

One of today's biggest concerns is the increase of energetic needs, especially in the developed countries. Among various clean energies, wind energy is one of the technologies that assume greater importance on the sustainable development of humanity. Despite wind turbines had been developed and studied over the years, there are phenomena that haven't been yet fully understood. This work studies the soil-structure interaction that occurs on a wind turbine's foundation composed by a group of piles that is under dynamic loads caused by wind. This problem assumes special importance when the foundation is implemented on locations where safety criteria are very demanding, like the case of a foundation mounted on a dike. To the phenomenon of interaction between two piles and the soil between them it's given the name of pile-soil-pile interaction. It is known that such behavior is frequency dependent, and therefore, on this work evaluation of relevant frequencies for the intended analysis is held. During the development of this thesis, two methods were selected in order to assess pile-soil-pile interaction, being one of analytical nature and the other of numerical origin. The analytical solution was recently developed and its called Generalized pile-soil-pile theory, while for the numerical method the commercial nite element software PLAXIS 3D was used. A study of applicability of the numerical method is also done comparing the given solution by the nite element methods with a rigorous solution widely accepted by the majority of the authors.

An 40Ar/39Ar, Rb/Sr, and stable isotope study of micas in low-grade fold-and-thrust belt: An example from the Swiss Helvetic Alps

White micas in carbonate-rich tectonites and a few other rock types of large thrusts in the Swiss Helvetic fold-and-thrust belt have been analyzed by Ar-40/Ar-39 and Rb/Sr techniques to better constrain the timing of Alpine deformation for this region. Incremental Ar-40/Ar-39 heating experiments of 25 weakly metamorphosed (anchizone to low greenschist) samples yield plateau and staircase spectra. We interpret most of the staircase release spectra result from variable mixtures of syntectonic (neoformed) and detrital micas. The range in dates obtained within individual spectra depends primarily on the duration of mica nucleation and growth, and relative proportions of neoformed and detrital mica. Rb/Sr analyses of 12 samples yield dates of ca. 10-39 Ma (excluding one anomalously young sample). These dates are slightly younger than the Ar-40/Ar-39 total gas dates obtained for the same samples. The Rb/ Sr dates were calculated using initial Sr-87/Sr-86 ratios obtained from the carbonate-dominated host rocks, which are higher than normal Mesozoic carbonate values due to exchange with fluids of higher Sr-87/Sr-86 ratios (and lower O-18/O-16 ratios). Model dates calculated using Sr-87/Sr-86 values typical of Mesozoic marine carbonates more closely approximate the Ar-40/Ar-39 total gas dates for most of the samples. The similarities of Rb/Sr and Ar-40/Ar-39 total gas dates are consistent with limited amounts of detrital mica in the samples. The delta(18)O values range from 24-15%. (VSMOW) for 2-6 mum micas and 27-16parts per thousand for the carbonate host rocks. The carbonate values are significantly lower than their protolith values due to localized fluid-rock interaction and fluid flow along most thrust surfaces. Although most calcite-mica pairs are not in oxygen isotope equilibrium at temperatures of ca. 200-400 degreesC, their isotopic fractionations are indicative of either 1) partial exchange between the minerals and a common external fluid, or 2) growth or isotopic exchange of the mica with the carbonate after the carbonate had isotopically exchanged with an external fluid. The geological significance of these results is not easily or uniquely determined, and exemplifies the difficulties inherent in dating very fine-grained micas of highly deformed tectonites in low-grade metamorphic terranes. Two generalizations can be made regarding the dates obtained from the Helvetic thrusts: 1) samples from the two highest thrusts (Mt. Gond and Sublage) have all of their Ar-40/Ar-39 steps above 20 Ma, and 2) most samples from the deepest Helvetic thrusts have steps (often accounting for more than 80% of Ar-39 release) between 15 and 25 Ma. These dates are consistent with the order of thrusting in the foreland-imbricating system and increase proportions of neoformed to detrital mica in the more metamorphosed hinterland and deeply buried portions of the nappe pile. Individual thrusts accommodated the majority of their displacement during their initial incorporation into the foreland-imbricating system, and some thrusts remained active or were reactivated down to 15 Ma.

Sr, C and O isotope systematics in the Pucara Basin, central Peru - Comparison between Mississippi Valley-type deposits and barren areas

A combined Sr, O and C isotope study has been carried out in the Pucara basin, central Peru, to compare local isotopic trends of the San Vicente and Shalipayco Zn-Pb Mississippi Valley-type (MVT) deposits with regional geochemical patterns of the sedimentary host basin. Gypsum, limestone and regional replacement dolomite yield Sr-87/Sr-86 ratios that fall within or slightly below the published range of seawater Sr-87/Sr-86 values for the Lower Jurassic and the Upper Triassic. Our data indicate that the Sr isotopic composition of seawater between the Hettangian and the Toarcian may extend to lower Sr-87/Sr-86 ratios than previously published values. An Sr-87-enrichment is noted in (1) carbonate rocks from the lowermost part of the Pucara basin, and (2) different carbonate generations at the MVT deposits. This indicates that host rocks at MVT deposits and in the lowermost part of the carbonate sequence interacted with Sr-87-enriched fluids. The fluids acquired their radiogenic nature by interaction with lithologies underlying the carbonate rocks of the Pucara basin. The San Ramon granite, similar Permo-Triassic intrusions and their elastic derivatives in the Mitu Group are likely sources of radiogenic Sr-87. The Brazilian shield and its erosion products are an additional potential source of radiogenic Sr-87. Volcanic rocks of the Mitu Group are not a significant source for radiogenic Sr-87; however, molasse-type sedimentary rocks and volcaniclastic rocks cannot be ruled out as a possible source of radiogenic Sr-87. The marked enrichment in Sr-87 of carbonates toward the lower part of the Pucara Group is accompanied by only a slight decrease in delta(18)O values and essentially no change in delta(13)C values, whereas replacement dolomite and sparry carbonates at the MVT deposits display a coherent trend of progressive Sr-87-enrichment, and O-18- and C-13-depletion. The depletion in O-18 in carbonates from the MVT deposits are likely related to a temperature increase, possibly coupled with a O-18-enrichment of the ore-forming fluids. Progressively lower delta(13)C values throughout the paragenetic sequence at the MVT deposits are interpreted as a gradually more important contribution from organically derived carbon. Quantitative calculations show that a single fluid-rock interaction model satisfactorily reproduces the marked Sr-87-enrichment and the slight decrease in delta(18)O values in carbonate rocks from the lower part of the Pucara Group. By contrast, the isotopic covariation trends of the MVT deposits are better reproduced by a model combining fluid mixing and fluid-rock interaction. The modelled ore-bearing fluids have a range of compositions between a hot, saline, radiogenic brine that had interacted with lithologies underlying the Pucara sequence and cooler, dilute brines possibly representing local fluids within the Pucara sequence. The composition of the local fluids varies according to the nature of the lithologies present in the neighborhood of the different MVT deposits. The proportion of the radiogenic fluid in the modelled fluid mixtures interacting with the carbonate host rocks at the MVT deposits decreases as one moves up in the stratigraphic sequence of the Pucara Group.

Segmentation and grid generation for numerical simulations of vad connections with patient-specific data

Objectives: We are interested in the numerical simulation of the anastomotic region comprised between outflow canula of LVAD and the aorta. Segmenta¬tion, geometry reconstruction and grid generation from patient-specific data remain an issue because of the variable quality of DICOM images, in particular CT-scan (e.g. metallic noise of the device, non-aortic contrast phase). We pro¬pose a general framework to overcome this problem and create suitable grids for numerical simulations.Methods: Preliminary treatment of images is performed by reducing the level window and enhancing the contrast of the greyscale image using contrast-limited adaptive histogram equalization. A gradient anisotropic diffusion filter is applied to reduce the noise. Then, watershed segmentation algorithms and mathematical morphology filters allow reconstructing the patient geometry. This is done using the InsightToolKit library (www.itk.org). Finally the Vascular Model¬ing ToolKit (www.vmtk.org) and gmsh (www.geuz.org/gmsh) are used to create the meshes for the fluid (blood) and structure (arterial wall, outflow canula) and to a priori identify the boundary layers. The method is tested on five different patients with left ventricular assistance and who underwent a CT-scan exam.Results: This method produced good results in four patients. The anastomosis area is recovered and the generated grids are suitable for numerical simulations. In one patient the method failed to produce a good segmentation because of the small dimension of the aortic arch with respect to the image resolution.Conclusions: The described framework allows the use of data that could not be otherwise segmented by standard automatic segmentation tools. In particular the computational grids that have been generated are suitable for simulations that take into account fluid-structure interactions. Finally the presented method features a good reproducibility and fast application.

Textural, chemical, and isotopic effects of late-magmatic carbonatitic fluids in the carbonatite-syenite Tamazeght complex, High Atlas Mountains, Morocco

Carbonatites of the Eocene Tamazeght complex, High Atlas Mountains, Morocco, consist of calciocarbonatites (alvikite and sovite dykes) and magnesiocarbonatites (diatreme breccias and dykes rocks). These are associated with ultramafic, shonkinitic, gabbroic to monzonitic and various foid syenitic silicate units. Stable and radiogenic isotope compositions for carbonatites and silicate rocks indicate that they share a common source in the mantle, although for some carbonatitic samples contamination with sedimentary rocks seems important. The observed isotopic heterogeneity is mainly attributed to source characteristics, fractional crystallization (accompanied by various degrees of assimilation), and late- to post-magmatic fluid-rock interaction. During the late fluid-rock interaction, Sr, Mn, and possibly also Fe were mobilized and redistributed to form secondary carbonate minerals in carbonatites. These fluids also penetrated into the adjacent syenitic rocks, causing enrichment in the same elements.

Volatile transport during the crystallization of anatectic melts: Oxygen, boron and hydrogen stable isotope study on the metamorphic complex of Naxos, Greece

Crystallization of anatectic melts in high-temperature metamorphic terrains releases volatile-rich magmas that can be transported into adjacent lithologies. This study addresses the variations in the oxygen, boron and hydrogen isotopic composition of aplite-pegmatite dikes that formed during the crystallization of anatectic melts in regional high-temperature metamorphism on the island of Naxos, Greece, and propagated upward into the overlying sequences of metamorphic schist. The transport distance of these dikes was increased through a significant horizontal component of travel that was imposed by contemporaneous low-angle extensional shearing. Laser fluorination oxygen isotope analyses of quartz, tourmaline, garnet, and biotite mineral separates from the aplite-pegmatite dikes show a progressive rise in delta(18)O values with increasing distance from the core. Oxygen isotope fractionations among quartz, tourmaline, and garnet show temperature variations from > 700degreesC down to similar to400degreesC. This range is considered to reflect isotopic fractionation beginning with crystallization at high temperatures in water-undersaturated conditions and then evolving through lower temperature crystallization and retrograde sub-solidus exchange. Two processes are examined for the cause of the progressive increase in delta(18)O values: (1) heterogeneous delta(18)O sources and (2) fluid-rock exchange between the aplite/pegmatite magmas and their host rock. Although the former process cannot be ruled out, there is as yet no evidence in the exposed sequences on Naxos for the presence of a suitable high delta(18)O magma source. In contrast, a tendency for the delta(18)O of quartz in the aplite/pegmatite dikes to approach that of the quartz in the metamorphic rock suggests that fluid-rock exchange with the host rock may potentially be an important process. Advection of fluid into the magma is examined based on Darcian fluid flow into an initially water-undersaturated buoyantly propagating aplitic dike magma. It is shown that such advective flow could only account for part of the O-18-enrichment, unless it were amplified by repeated injection of magma pulses, fluid recycling, and deformation-assisted post-crystallization exchange. The mechanism is, however, adequate to account for hydrogen isotope equilibration between dike and host rock. In contrast, variations in the delta(11)B values of tourmalines suggest that B-11/B-10 fractionation during crystallization and/or magma degassing was the major control of boron geochemistry rather than fluid-rock interaction and that the boron isotopic system was decoupled from that of oxygen. Copyright (C) 2003 Elsevier Ltd.

Origin of mineralizing fluids of the sediment-hosted Navachab gold mine, Namibia: Constraints from stable (O, H, C, S) isotopes

The Navachab gold mine in the Damara belt of central Namibia is characterized by a polymetallic Au-Bi-As-Cu-Ag ore assemblage, including pyrrhotite, chalcopyrite, sphalerite, arsenopyrite, bismuth, gold, bismuthinite, and bismuth tellurides. Gold is hosted by quartz sulfide veins and semimassive sulfide lenses that are developed in a near-vertical sequence of shelf-type metasedimentary rocks, including marble, calcsilicate rock, and biotite schist. The sequence has been intruded by abundant syntectonic lamprophyre, aplite, and pegmatite dikes, documenting widespread igneous activity coeval with mineralization. The majority of quartz from the veins has delta(18)O values of 14 to 15 per mil (V-SMOW). The total variations in delta(18)O values of the biotite schist and calcsilicate rock are relatively small (12-14 parts per thousand), whereas the marble records steep gradients in delta(18)O values (17-21 parts per thousand), the lowest values being recorded at the vein margins. Despite this, there is no correlation between delta(18)O and delta(13)C values and the carbonate content of the rocks, indicating that fluid-rock interaction alone cannot explain the isotopic gradients. In addition, the marble records increased delta(13)C values at the contact to the veins, possibly related to a change in the physicochemical conditions during fluid-rock interaction. Gold is interpreted to have precipitated in equilibrium with metamorphic find (delta(18)O 12-14 parts per thousand; delta D = -40 to -60 parts per thousand) at peak metamorphic conditions of ca. 550 degrees C and 2 kbars, consistent with isotopic fractionations between coexisting calcite, garnet, and clinopyroxene in the alteration halos. The most likely source of the mineralizing fluid was a midcrustal fluid in equilibrium with the Damaran metapelites that underwent prograde metamorphism at amphibolite- to granulite-facies grades. Although there is no isotopic evidence for the contribution of magmatic fluids, they may have been important in contributing to the overall hydraulic regime and high apparent geothermal gradients (ca. 80 degrees C/km(-1)) in the mine area.

The magmatic to hydrothermal evolution of the intrusive Mont Saint-Hilaire Complex: Insights into the late-stage evolution of peralkaline rocks

The Cretaceous Mont Saint-Hilaire complex (Quebec, Canada) comprises three major rock units that were emplaced in the following sequence: (I) gabbros; (II) diorites; (III) diverse partly agpaitic foid syenites. The major element compositions of the rock-forming minerals, age-corrected Nd and oxygen isotope data for mineral separates and trace element data of Fe-Mg silicates from the various lithologies imply a common source for all units. The distribution of the rare earth elements in clinopyroxene from the gabbros indicates an ocean island basalt type composition for the parental magma. Gabbros record temperatures of 1200 to 800 degrees C, variable silica activities between 0 center dot 7 and 0 center dot 3, and f(O2) values between -0 center dot 5 and +0 center dot 7 (log delta FMQ, where FMQ is fayalite-magnetite-quartz). The diorites crystallized under uniform a(SiO2) (a(SiO2) = 0 center dot 4-0 center dot 5) and more reduced f(O2) conditions (log delta FMQ similar to-1) between similar to 1100 and similar to 800 degrees C. Phase equilibria in various foid syenites indicate that silica activities decrease from 0 center dot 6-0 center dot 3 at similar to 1000 degrees C to < 0 center dot 3 at similar to 550 degrees C. Release of an aqueous fluid during the transition to the hydrothermal stage caused a(SiO2) to drop to very low values, which results from reduced SiO(2) solubilities in aqueous fluids compared with silicate melts. During the hydrothermal stage, high water activities stabilized zeolite-group minerals. Fluid inclusions record a complex post-magmatic history, which includes trapping of an aqueous fluid that unmixed from the restitic foid syenitic magma. Cogenetic aqueous and carbonic fluid inclusions reflect heterogeneous trapping of coexisting immiscible external fluids in the latest evolutionary stage. The O and C isotope characteristics of fluid-inclusion hosted CO(2) and late-stage carbonates imply that the surrounding limestones were the source of the external fluids. The mineral-rich syenitic rocks at Mont Saint-Hilaire evolved as follows: first, alkalis, high field strength and large ion lithophile elements were pre-enriched in the (late) magmatic and subsequent hydrothermal stages; second, percolation of external fluids in equilibrium with the carbonate host-rocks and mixing processes with internal fluids as well as fluid-rock interaction governed dissolution of pre-existing minerals, element transport and precipitation of mineral assemblages determined by locally variable parameters. It is this hydrothermal interplay between internal and external fluids that is responsible for the mineral wealth found at Mont Saint-Hilaire.