Study of Tsunami-Induced Fluid and Debris Load on Bridges using the Material Point Method

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

Dynamic fluid-structure interactions using finite volume unstructured mesh procedures

A three-dimensional finite volume, unstructured mesh (FV-UM) method for dynamic fluid–structure interaction (DFSI) is described. Fluid structure interaction, as applied to flexible structures, has wide application in diverse areas such as flutter in aircraft, wind response of buildings, flows in elastic pipes and blood vessels. It involves the coupling of fluid flow and structural mechanics, two fields that are conventionally modelled using two dissimilar methods, thus a single comprehensive computational model of both phenomena is a considerable challenge. Until recently work in this area focused on one phenomenon and represented the behaviour of the other more simply. More recently, strategies for solving the full coupling between the fluid and solid mechanics behaviour have been developed. A key contribution has been made by Farhat et al. [Int. J. Numer. Meth. Fluids 21 (1995) 807] employing FV-UM methods for solving the Euler flow equations and a conventional finite element method for the elastic solid mechanics and the spring based mesh procedure of Batina [AIAA paper 0115, 1989] for mesh movement. In this paper, we describe an approach which broadly exploits the three field strategy described by Farhat for fluid flow, structural dynamics and mesh movement but, in the context of DFSI, contains a number of novel features: a single mesh covering the entire domain, a Navier–Stokes flow, a single FV-UM discretisation approach for both the flow and solid mechanics procedures, an implicit predictor–corrector version of the Newmark algorithm, a single code embedding the whole strategy.

Development of a model to predict the life of pneumatic conveyor bends subject to erosive wear

The erosion processes resulting from flow of fluids (gas-solid or liquid-solid) are encountered in nature and many industrial processes. The common feature of these erosion processes is the interaction of the fluid (particle) with its boundary thus resulting in the loss of material from the surface. This type of erosion in detrimental to the equipment used in pneumatic conveying systems. The puncture of pneumatic conveyor bends in industry causes several problems. Some of which are: (1) Escape of the conveyed product causing health and dust hazard; (2) Repairing and cleaning up after punctures necessitates shutting down conveyors, which will affect the operation of the plant, thus reducing profitability. The most common occurrence of process failure in pneumatic conveying systems is when pipe sections at the bends wear away and puncture. The reason for this is particles of varying speed, shape, size and material properties strike the bend wall with greater intensity than in straight sections of the pipe. Currently available models for predicting the lifetime of bends are inaccurate (over predict by 80%. The provision of an accurate predictive method would lead to improvements in the structure of the planned maintenance programmes of processes, thus reducing unplanned shutdowns and ultimately the downtime costs associated with these unplanned shutdowns. This is the main motivation behind the current research. The paper reports on two aspects of the first phases of the study-undertaken for the current project. These are (1) Development and implementation; and (2) Testing of the modelling environment. The model framework encompasses Computational Fluid Dynamics (CFD) related engineering tools, based on Eulerian (gas) and Lagrangian (particle) approaches to represent the two distinct conveyed phases, to predict the lifetime of conveyor bends. The method attempts to account for the effect of erosion on the pipe wall via particle impacts, taking into account the angle of attack, impact velocity, shape/size and material properties of the wall and conveyed material, within a CFD framework. Only a handful of researchers use CFD as the basis of predicting the particle motion, see for example [1-4] . It is hoped that this would lead to more realistic predictions of the wear profile. Results, for two, three-dimensional test cases using the commercially available CFD PHOENICS are presented. These are reported in relation to the impact intensity and sensitivity to the inlet particle distributions.

A Numerical study of fluid flow in peristaltic pumps and polymeric elastic pipes

This thesis work deals with a mathematical description of flow in polymeric pipe and in a specific peristaltic pump. This study involves fluid-structure interaction analysis in presence of complex-turbulent flows treated in an arbitrary Lagrangian-Eulerian (ALE) framework. The flow simulations are performed in COMSOL 4.4, as 2D axial symmetric model, and ABAQUS 6.14.1, as 3D model with symmetric boundary conditions. In COMSOL, the fluid and structure problems are coupled by monolithic algorithm, while ABAQUS code links ABAQUS CFD and ABAQUS Standard solvers with single block-iterative partitioned algorithm. For the turbulent features of the flow, the fluid model in both codes is described by RNG k-ϵ. The structural model is described, on the basis of the pipe material, by Elastic models or Hyperelastic Neo-Hookean models with Rayleigh damping properties. In order to describe the pulsatile fluid flow after the pumping process, the available data are often defective for the fluid problem. Engineering measurements are normally able to provide average pressure or velocity at a cross-section. This problem has been analyzed by McDonald's and Womersley's work for average pressure at fixed cross section by Fourier analysis since '50, while nowadays sophisticated techniques including Finite Elements and Finite Volumes exist to study the flow. Finally, we set up peristaltic pipe simulations in ABAQUS code, by using the same model previously tested for the fl uid and the structure.

Finite-size particles, advection, and chaos: A collective phenomenon of intermittent bursting

We consider finite-size particles colliding elastically, advected by a chaotic flow. The collisionless dynamics has a quasiperiodic attractor and particles are advected towards this attractor. We show in this work that the collisions have dramatic effects in the system's dynamics, giving rise to collective phenomena not found in the one-particle dynamics. In particular, the collisions induce a kind of instability, in which particles abruptly spread out from the vicinity of the attractor, reaching the neighborhood of a coexisting chaotic saddle, in an autoexcitable regime. This saddle, not present in the dynamics of a single particle, emerges due to the collective particle interaction. We argue that this phenomenon is general for advected, interacting particles in chaotic flows.

Positioning of Deadeners for Vibration Reduction in Vehicle Roof Using Embedded Sensitivity

The noise, vibration and harshness (NVH) performance of passenger vehicles strongly depends on the fluid-structure interaction between the air in the vehicle cavity and the sheet metal structure of the vehicle. Most of the noise and vibration problems related to this interaction come from resonance peaks of the sheet metal, which are excited by external forces (road, engine, and wind). A reduction in these resonance peaks can be achieved by applying bitumen damping layers, also called deadeners, in the sheet metal. The problem is where these deadeners shall be fixed, which is usually done in a trial-and-error basis. In this work, one proposes the use of embedded sensitivity to locate the deadeners in the sheet metal of the vehicle, more specifically in the vehicle roof. Experimental frequency response functions (FRFs) of the roof are obtained and the data are processed by adopting the embedded sensitivity method, thus obtaining the sensitivity of the resonance peaks on the local increase in damping due to the deadeners. As a result, by examining the sensitivity functions, one can find the optimum location of the deadeners that maximize their effect in reducing the resonance peaks of interest. After locating the deadeners in the optimum positions, it was possible to verify a strong reduction in resonance peaks of the vehicle roof, thus showing the efficiency of the procedure. The main advantage of this procedure is that it only requires FRF measurements of the vehicle in its original state not needing any previous modification of the vehicle structure to find the sensitivity functions. [DOI: 10.1115/1.4000769]

Quantum dynamics of three coupled atomic Bose-Einstein condensates

The simplest model of three coupled Bose-Einstein condensates is investigated using a group theoretical method. The stationary solutions are determined using the SU(3) group under the mean-field approximation. This semiclassical analysis, using system symmetries, shows a transition in the dynamics of the system from self trapping to delocalization at a critical value for the coupling between the condensates. The global dynamics are investigated by examination of the stable points, and our analysis shows that the structure of the stable points depends on the ratio of the condensate coupling to the particle-particle interaction, and undergoes bifurcations as this ratio is varied. This semiclassical model is compared to a full quantum treatment, which also displays a dynamical transition. The quantum case has collapse and revival sequences superimposed on the semiclassical dynamics, reflecting the underlying discreteness of the spectrum. Nonzero circular current states are also demonstrated as one of the higher-dimensional effects displayed in this system.

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.

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.

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.