The mineralizing fluid in the Piaba gold deposit, São Luís cratonic fragment (NW-Maranhão, Brazil) based on fluid inclusion studies in quartz veins

O depósito aurífero de Piaba tornou-se a primeira mina em operação no fragmento cratônico São Luís, noroeste do Maranhão. Seu ambiente geológico compreende rochas metavulcanossedimentares do Grupo Aurizona e granitoides da Suíte Tromaí, entre outras unidades menores, formadas em ambiente de arcos de ilhas entre 2240 e 2150 Ma, juntamente com outras unidades menores. A mineralização é caracterizada por uma trama stockwork de veios e vênulas de quartzo com seus halos de alteração (clorita + muscovita + carbonato + pirita + calcopirita e ouro) hospedada em um granodiorito granofírico fino (Granófiro Piaba) e em rocha subvulcânica andesítica do Grupo Aurizona. O corpo mineralizado é espacialmente limitado à zona de cisalhamento ENE-WSW rúptil-dúctil (Falha Piaba). Estudos petrográficos, microtermométricos e por espectroscopia microRaman no quartzo definiram inclusões aquo-carbônicas bifásicas e trifásicas, produzidas por aprisionamento heterogêneo durante separação de fases, e fluidos aquosos tardios. A solução mineralizadora corresponde a um fluido aquo-carbônico composto por CO₂ (5 - 24 mol%, densidade de 0,96 - 0,99 g/cm³), H₂O (74 - 93 mol%), N₂ (< 1 mol%), CH₄ (<1mol%) e 5,5 % em peso NaCl equivalente. O minério depositou a 267 - 302ºC e 1,25 - 2,08 kbar, correspondendo a profundidades de 4 a 7 km, em consonância com o regime estrutural. A composição e o intervalo de P-T do fluido mineralizador, combinados com o caráter redutor (log ƒ_O2 -31,3 a -34,3) e a sulfetação da rocha hospedeira, sugerem que o ouro foi transportado como um complexo sulfetado. O minério foi depositado em consequência da separação de fase, redução da atividade de enxofre e da ƒO₂ pela interação fluido-rocha.

Application of reactive fluid transport modeling to hydrothermal systems

A one-dimensional multi-component reactive fluid transport algorithm, 1DREACT (Steefel, 1993) was used to investigate different fluid-rock interaction systems. A major short coming of mass transport calculations which include mineral reactions is that solid solutions occurring in many minerals are not treated adequately. Since many thermodynamic models of solid solutions are highly non-linear, this can seriously impact on the stability and efficiency of the solution algorithms used. Phase petrology community saw itself faced with a similar predicament 10 years ago. To improve performance and reliability, phase equilibrium calculations have been using pseudo compounds. The same approach is used here in the first, using the complex plagioclase solid solution as an example. Thermodynamic properties of a varying number of intermediate plagioclase phases were calculated using ideal molecular, Al-avoidance, and non-ideal mixing models. These different mixing models can easily be incorporated into the simulations without modification of the transport code. Simulation results show that as few as nine intermediate compositions are sufficient to characterize the diffusional profile between albite and anorthite. Hence this approach is very efficient, and can be used with little effort. A subsequent chapter reports the results of reactive fluid transport modeling designed to constrain the hydrothermal alteration of Paleoproterozoic sediments of the Southern Lake Superior region. Field observations reveal that quartz-pyrophyllite (or kaolinite) bearing assemblages have been transformed into muscovite-pyrophyllite-diaspore bearing assemblages due to action of fluids migrating along permeable flow channels. Fluid-rock interaction modeling with an initial qtz-prl assemblage and a K-rich fluid simulates the formation of observed mineralogical transformation. The bulk composition of the system evolves from an SiO2-rich one to an Al2O3+K2O-rich one. Simulations show that the fluid flow was up-temperature (e.g. recharge) and that fluid was K-rich. Pseudo compound approach to include solid solutions in reactive transport models was tested in modeling hydrothermal alteration of Icelandic basalts. Solid solutions of chlorites, amphiboles and plagioclase were included as the secondary mineral phases. Saline and fresh water compositions of geothermal fluids were used to investigate the effect of salinity on alteration. Fluid-rock interaction simulations produce the observed mineral transformations. They show that roughly the same alteration minerals are formed due to reactions with both types of fluid which is in agreement with the field observations. A final application is directed towards the remediation of nitrate rich groundwaters. Removal of excess nitrate from groundwater by pyrite oxidation was modeled using the reactive fluid transport algorithm. Model results show that, when a pyrite-bearing, permeable zone is placed in the flow path, nitrate concentration in infiltrating water can be significantly lowered, in agreement with proposals from the literature. This is due to nitrogen reduction. Several simulations investigate the efficiency of systems with different mineral reactive surface areas, reactive barrier zone widths, and flow rates to identify the optimum setup.

Investigating the Effects of Various Crowd Characteristics on the Dynamic Properties of an Occupied Structure

One of the challenges for structural engineers during design is considering how the structure will respond to crowd-induced dynamic loading. It has been shown that human occupants of a structure do not simply add mass to the system when considering the overall dynamic response of the system, but interact with it and may induce changes of the dynamic properties from those of the empty structure. This study presents an investigation into the human-structure interaction based on several crowd characteristics and their effect on the dynamic properties of an empty structure. The dynamic properties including frequency, damping, and mode shapes were estimated for a single test structure by means of experimental modal analysis techniques. The same techniques were utilized to estimate the dynamic properties when the test structure was occupied by a crowd with different combinations of size, posture, and distribution. The goal of this study is to isolate the occupant characteristics in order to determine the significance of each to be considered when designing new structures to avoid crowd serviceability issues. The results are presented and summarized based on the level of influence of each characteristic. The posture that produces the most significant effects based on the scope of this research is standing with bent knees with a maximum decrease in frequency of the first mode of the empty structure by 32 percent atthe highest mass ratio. The associated damping also increased 36 times the damping of the empty structure. In addition to the analysis of the experimental data, finite element models and a two degree-of-freedom model were created. These models were used to gain an understanding of the test structure, model a crowd as an equivalent mass, and also to develop a single degree-of-freedom (SDOF) model to best represent a crowd of occupants based on the experimental results. The SDOF models created had an averagefrequency of 5.0 Hz, within the range presented in existing biomechanics research, and combined SDOF systems of the test structure and crowd were able to reproduce the frequency and damping ratios associated with experimental tests. Results of this study confirmed the existence of human-structure interaction andthe inability to simply model a crowd as only additional mass. The two degree-offreedom model determined was able to predict the change in natural frequency and damping ratio for a structure occupied by multiple group sizes in a single posture. These results and model are the preliminary steps in the development of an appropriate methodfor modeling a crowd in combination with a more complex FE model of the empty structure.

Fluid-related inclusions in Alpine high-pressure peridotite reveal trace element recycling during subduction-zone dehydration of serpentinized mantle (Cima di Gagnone, Swiss Alps).

Serpentinites release at sub-arc depths volatiles and several fluid-mobile trace elements found in arc magmas. Constraining element uptake in these rocks and defining the trace element composition of fluids released upon serpentinite dehydration can improve our understanding of mass transfer across subduction zones and to volcanic arcs. The eclogite-facies garnet metaperidotite and chlorite harzburgite bodies embedded in paragneiss of the subduction melange from Cima di Gagnone derive from serpentinized peridotite protoliths and are unique examples of ultramafic rocks that experienced subduction metasomatism and devolatilization. In these rocks, metamorphic olivine and garnet trap polyphase inclusions representing the fluid released during high-pressure breakdown of antigorite and chlorite. Combining major element mapping and laser-ablation ICP-MS bulk inclusion analysis, we characterize the mineral content of polyphase inclusions and quantify the fluid composition. Silicates, Cl-bearing phases, sulphides, carbonates, and oxides document post-entrapment mineral growth in the inclusions starting immediately after fluid entrapment. Compositional data reveal the presence of two different fluid types. The first (type A) records a fluid prominently enriched in fluid-mobile elements, with Cl, Cs, Pb, As, Sb concentrations up to 10(3) PM (primitive mantle), similar to 10(2) PM Tit Ba, while Rb, B, Sr, Li, U concentrations are of the order of 10(1) PM, and alkalis are similar to 2 PM. The second fluid (type B) has considerably lower fluid-mobile element enrichments, but its enrichment patterns are comparable to type A fluid. Our data reveal multistage fluid uptake in these peridotite bodies, including selective element enrichment during seafloor alteration, followed by fluid-rock interaction along with subduction metamorphism in the plate interface melange. Here, infiltration of sediment-equilibrated fluid produced significant enrichment of the serpentinites in As, Sb, B, Pb, an enriched trace element pattern that was then transferred to the fluid released at greater depth upon serpentine dehydration (type A fluid). The type B fluid hosted by garnet may record the composition of the chlorite breakdown fluid released at even greater depth. The Gagnone study-case demonstrates that serpentinized peridotites acquire water and fluid-mobile elements during ocean floor hydration and through exchange with sediment-equilibrated fluids in the early subduction stages. Subsequent antigorite devolatilization at subarc depths delivers aqueous fluids to the mantle wedge that can be prominently enriched in sediment-derived components, potentially triggering arc magmatism without the need of concomitant dehydration/melting of metasediments or altered oceanic crust.

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.

Vortex-induced vibrations on bridges : the interaction mechanism in complex sections and a new proposed semi-empirical model

El presente trabajo de investigación se ocupa del estudio de las vibraciones verticales inducidas por vórtices (VIV) en aquellos puentes que, por sus características geométricas y propiedades dinámicas, muestran cierta sensibilidad este tipo de fenómeno aeroelástico. El objeto principal es el análisis del mecanismo de interacción viento-estructura sobre secciones no fuseladas de geometría simple, con objeto de realizar una adecuada caracterización del problema y poder abordar posteriormente el análisis de otras secciones de geometría más compleja, representativas de los principales elementos estructurales de los puentes, como arcos, tableros, torres y pilas. Este aspecto es fundamental durante la fase de diseño del puente, donde deberán tenerse en cuenta también una serie de detalles que pueden influir significativamente su sensibilidad ante problemas aerodinámicos, como la morfología y dimensiones principales de la sección transversal del tablero, la disposición de barreras de seguridad y barreras cortaviento, o las riostras que unen diferentes elementos estructurales. La configuración de dos elementos en tándem o la construcción de un puente en las inmediaciones de otro existente son otros aspectos a considerar respecto a la sensibilidad frente a efectos aeroelásticos. El estudio se ha llevado a cabo principalmente mediante la implementación de simulaciones numéricas que reproducen la interacción entre la corriente de aire y secciones representativas de modelos estructurales, a partir de un código CFD basado en el método de las partículas de vórtices (VPM), siguiendo por tanto un esquema Lagrangiano. Los resultados han sido validados con datos experimentales existentes, valores procedentes de ensayos en túnel de viento y registros reales a partir de diferentes casos de estudio: Alconétar (2006), Niterói (1980), Trans- Tokyo Bay (1995) y Volgogrado (2010). Finalmente, se propone un modelo semi-empírico para la estimación del rango de velocidades críticas y amplitudes de oscilación basado en la utilización de las derivadas de flameo de Scanlan, y la densidad espectral de las fuerzas aerodinámicas en el dominio de la frecuencia. The present research work concerns the study of vertical vortex-induced vibrations (VIV) in bridges which show certain sensitivity to this type of aeroelastic phenomenon. It focuses on the analysis of the wind-structure interaction mechanism on bluff sections, with the objective of making a good characterisation of the problem and subsequently addressing the analysis of sections with a complex geometry, which are representative of the bridge structural elements, such as arches, decks, towers and piers. This issue is of relative importance during the bridge design phase, since minor details of the aforementioned elements can significantly influence its sensitivity to aerodynamic problems. The shape and main dimensions of the deck cross section, the addition of safety barriers and windshields, the presence of braces to enhance the structure mechanical properties, the utilisation of cross sections in tandem arrangement, or the erection of a new bridge in the vicinity of another existing one are some of the aspects to be considered regarding the sensitivity to the aeroelastic effects. The study has been carried out mainly through the implementation of numerical simulations that reproduces the interaction between the airflow and the representative cross section of a structural bridge model, by the use of a CFD code based on the vortex particle method (VPM), thus following a Lagrangian scheme. The results have been validated with existing experimental data, values from wind tunnel tests and full scale observations from the different case studies: Alconétar (2006), Niterói (1980), Trans-Tokyo Bay (1995) and Volgograd (2010). Finally, a new semi-empirical model is proposed for the estimation of the critical wind velocity ranges and oscillation amplitudes based on the use of the Scanlan’s flutter derivatives and the power spectral density of aerodynamic force time history in the frequency domain.

On the development of advanced techniques for mixed-elastohydrodynamic lubrification modelling of journal and sliding bearing systems.

The present thesis is focused on the development of a thorough mathematical modelling and computational solution framework aimed at the numerical simulation of journal and sliding bearing systems operating under a wide range of lubrication regimes (mixed, elastohydrodynamic and full film lubrication regimes) and working conditions (static, quasi-static and transient conditions). The fluid flow effects have been considered in terms of the Isothermal Generalized Equation of the Mechanics of the Viscous Thin Films (Reynolds equation), along with the massconserving p-Ø Elrod-Adams cavitation model that accordingly ensures the so-called JFO complementary boundary conditions for fluid film rupture. The variation of the lubricant rheological properties due to the viscous-pressure (Barus and Roelands equations), viscous-shear-thinning (Eyring and Carreau-Yasuda equations) and density-pressure (Dowson-Higginson equation) relationships have also been taken into account in the overall modelling. Generic models have been derived for the aforementioned bearing components in order to enable their applications in general multibody dynamic systems (MDS), and by including the effects of angular misalignments, superficial geometric defects (form/waviness deviations, EHL deformations, etc.) and axial motion. The bearing exibility (conformal EHL) has been incorporated by means of FEM model reduction (or condensation) techniques. The macroscopic in fluence of the mixedlubrication phenomena have been included into the modelling by the stochastic Patir and Cheng average ow model and the Greenwood-Williamson/Greenwood-Tripp formulations for rough contacts. Furthermore, a deterministic mixed-lubrication model with inter-asperity cavitation has also been proposed for full-scale simulations in the microscopic (roughness) level. According to the extensive mathematical modelling background established, three significant contributions have been accomplished. Firstly, a general numerical solution for the Reynolds lubrication equation with the mass-conserving p - Ø cavitation model has been developed based on the hybridtype Element-Based Finite Volume Method (EbFVM). This new solution scheme allows solving lubrication problems with complex geometries to be discretized by unstructured grids. The numerical method was validated in agreement with several example cases from the literature, and further used in numerical experiments to explore its exibility in coping with irregular meshes for reducing the number of nodes required in the solution of textured sliding bearings. Secondly, novel robust partitioned techniques, namely: Fixed Point Gauss-Seidel Method (PGMF), Point Gauss-Seidel Method with Aitken Acceleration (PGMA) and Interface Quasi-Newton Method with Inverse Jacobian from Least-Squares approximation (IQN-ILS), commonly adopted for solving uid-structure interaction problems have been introduced in the context of tribological simulations, particularly for the coupled calculation of dynamic conformal EHL contacts. The performance of such partitioned methods was evaluated according to simulations of dynamically loaded connecting-rod big-end bearings of both heavy-duty and high-speed engines. Finally, the proposed deterministic mixed-lubrication modelling was applied to investigate the in fluence of the cylinder liner wear after a 100h dynamometer engine test on the hydrodynamic pressure generation and friction of Twin-Land Oil Control Rings.

Comparison of Lattice-Fluid Binary Parameters For Mixtures and Block Copolymers

The aim of this report is to discuss the method of determination of lattice-fluid binary interaction parameters by comparing well characterized immiscible blends and block copolymers of poly(methyl methacrylate) (PMMA) and poly(ϵ−caprolactone) (PCL). Experimental pressure-volume-temperature (PVT) data in the liquid state were correlated with the Sanchez—Lacombe (SL) equation of state with the scaling parameters for mixtures and copolymers obtained through combination rules of the characteristic parameters for the pure homopolymers. The lattice-fluid binary parameters for energy and volume were higher than those of block copolymers implying that the copolymers were more compatible due to the chemical links between the blocks. Therefore, a common parameter cannot account for both homopolymer blend and block copolymer phase behaviors based on current theory. As we were able to adjust all data of the mixtures with a single set of lattice-binary parameters and all data of the block copolymers with another single set we can conclude that both parameters did not depend on the composition for this system. This characteristic, plus the fact that the additivity law of specific volumes can be suitably applied for this system, allowed us to model the behavior of the immiscible blend with the SL equation of state. In addition, a discussion on the relationship between lattice-fluid binary parameters and the Flory–Huggins interaction parameter obtained from Leibler's theory is presented.

Lithologic, mineral, chemical and isotopic compositions of altered rocks with increasing distance across the eastern flank of the Juan de Fuca Ridge, ODP Leg 168 data

During ODP Leg 168, 10 sites were drilled across the eastern flank of the Juan de Fuca Ridge (JdFR), to examine the conditions of fluid-rock interaction in three distinct hydrothermal regimes (referred to as the Hydrothermal Transition (HT), Buried Basement (BB) and Rough Basement (RB) transects), extending over a ~120 km linear transect perpendicular to the spreading ridge. This was carried out in an attempt to constrain the conditions and processes that control the location, style and magnitude of low temperature (<150°C) fluid-rock interaction within this setting. This paper presents new data on the petrology, mineral chemistry and whole rock strontium and oxygen isotopic compositions of basalts from the eastern flank of the JdFR, in order to investigate the extent, style and sequence of low-temperature hydrothermal alteration and to establish how the hydrothermal regime evolved with time. Throughout the flank, a progressive sequence of low-temperature hydrothermal alteration has been identified, marked by changes in the dominant secondary mineral assemblage, changing from: chlorite+chlorite/smectite; to iron oyxhydroxide+celadonite; to saponite+/-pyrite; culminating at present with Ca- to CaMg(+/-Fe,Mn)-carbonate. The changes in secondary mineralogy have been used to infer a series of systematic shifts in the conditions of alteration that occurred as the basement moved off-axis and was progressively buried by sediment. In general, hydrothermal alteration of the uppermost oceanic crust commenced under open, oxidative conditions, with interaction between unmodified to slightly modified seawater and basaltic crust, to a regime in which circulation of a strongly modified seawater-derived fluid was more restricted, and alteration occurred under non-oxidative conditions. Across the flank, petrological observations and microprobe analyses indicate that the observed ranges in secondary mineral composition are directly related to changes in the geochemical and textural characteristics of the basement, as well as to interaction between fluids and phases from the four stages of alteration. This is suggestive of an increase in fluid-rock increased with time. Whole rock 87Sr/86Sr and d18O analyses of basalts from across the eastern flank of the JdFR reinforce petrological observations, with 87Sr/86Sr and d18O values slightly elevated above accepted pristine MORB values for this region. These results are consistent with an increase in the amount of fluid-rock interaction with time. Across the flank, enrichment in the 87Sr/86Sr and d18O relative to MORB, is influenced by a number of factors, including: local and regional variations in the crustal lithology and structure; the age of the crust; the extent of bulk rock alteration; and theoretically, the relative abundance of different isotopically-enriched secondary mineral phases in the crust.

Dam break wave of thixotropic fluid

Thixotropy is the characteristic of a fluid to form a gelled structure over time when it is not subjected to shearing, and to liquefy when agitated. Thixotropic fluids are commonly used in the construction industry (e.g., liquid concrete and drilling fluids), and related applications include some forms of mud flows and debris flows. This paper describes a basic study of dam break wave with thixotropic fluid. Theoretical considerations were developed based upon a kinematic wave approximation of the Saint-Venant equations down a prismatic sloping channel. A very simple thixotropic model, which predicts the basic theological trends of such fluids, was used. It describes the instantaneous state of fluid structure by a single parameter. The analytical solution of the basic flow motion and theology equations predicts three basic flow regimes depending upon the fluid properties and flow conditions, including the initial degree of jamming of the fluid (related to its time of restructuration at rest). These findings were successfully compared with systematic bentonite suspension experiments. The present work is the first theoretical analysis combining the basic principles of unsteady flow motion with a thixotropic fluid model and systematic laboratory experiments.

Seismic Response Analysis of a Full-Scale Base-Isolated Structure via Measurements and Modeling

The full-scale base-isolated structure studied in this dissertation is the only base-isolated building in South Island of New Zealand. It sustained hundreds of earthquake ground motions from September 2010 and well into 2012. Several large earthquake responses were recorded in December 2011 by NEES@UCLA and by GeoNet recording station nearby Christchurch Women's Hospital. The primary focus of this dissertation is to advance the state-of-the art of the methods to evaluate performance of seismic-isolated structures and the effects of soil-structure interaction by developing new data processing methodologies to overcome current limitations and by implementing advanced numerical modeling in OpenSees for direct analysis of soil-structure interaction.

This dissertation presents a novel method for recovering force-displacement relations within the isolators of building structures with unknown nonlinearities from sparse seismic-response measurements of floor accelerations. The method requires only direct matrix calculations (factorizations and multiplications); no iterative trial-and-error methods are required. The method requires a mass matrix, or at least an estimate of the floor masses. A stiffness matrix may be used, but is not necessary. Essentially, the method operates on a matrix of incomplete measurements of floor accelerations. In the special case of complete floor measurements of systems with linear dynamics, real modes, and equal floor masses, the principal components of this matrix are the modal responses. In the more general case of partial measurements and nonlinear dynamics, the method extracts a number of linearly-dependent components from Hankel matrices of measured horizontal response accelerations, assembles these components row-wise and extracts principal components from the singular value decomposition of this large matrix of linearly-dependent components. These principal components are then interpolated between floors in a way that minimizes the curvature energy of the interpolation. This interpolation step can make use of a reduced-order stiffness matrix, a backward difference matrix or a central difference matrix. The measured and interpolated floor acceleration components at all floors are then assembled and multiplied by a mass matrix. The recovered in-service force-displacement relations are then incorporated into the OpenSees soil structure interaction model.

Numerical simulations of soil-structure interaction involving non-uniform soil behavior are conducted following the development of the complete soil-structure interaction model of Christchurch Women's Hospital in OpenSees. In these 2D OpenSees models, the superstructure is modeled as two-dimensional frames in short span and long span respectively. The lead rubber bearings are modeled as elastomeric bearing (Bouc Wen) elements. The soil underlying the concrete raft foundation is modeled with linear elastic plane strain quadrilateral element. The non-uniformity of the soil profile is incorporated by extraction and interpolation of shear wave velocity profile from the Canterbury Geotechnical Database. The validity of the complete two-dimensional soil-structure interaction OpenSees model for the hospital is checked by comparing the results of peak floor responses and force-displacement relations within the isolation system achieved from OpenSees simulations to the recorded measurements. General explanations and implications, supported by displacement drifts, floor acceleration and displacement responses, force-displacement relations are described to address the effects of soil-structure interaction.

The Pico de Navas slump (Burgos, Spain): a large rocky landslide caused by underlying clayey sand fluidification

The Pico de Navas landslide was a large-magnitude rotational movement, affecting 50x106m3 of hard to soft rocks. The objectives of this study were: (1) to characterize the landslide in terms of geology, geomorphological features and geotechnical parameters; and (2) to obtain an adequate geomechanical model to comprehensively explain its rupture, considering topographic, hydro-geological and geomechanical conditions. The rupture surface crossed, from top to bottom: (a) more than 200 m of limestone and clay units of the Upper Cretaceous, affected by faults; and (b) the Albian unit of Utrillas facies composed of silty sand with clay (Kaolinite) of the Lower Cretaceous. This sand played an important role in the basal failure of the slide due to the influence of fine particles (silt and clay), which comprised on average more than 70% of the sand, and the high content presence of kaolinite (>40%) in some beds. Its geotechnical parameters are: unit weight (δ) = 19-23 KN/m3; friction angle (φ) = 13º-38º and cohesion (c) = 10-48 KN/m2. Its microstructure consists of accumulations of kaolinite crystals stuck to terrigenous grains, making clayey peds. We hypothesize that the presence of these aggregates was the internal cause of fluidification of this layer once wet. Besides the faulted structure of the massif, other conditioning factors of the movement were: the large load of the upper limestone layers; high water table levels; high water pore pressure; and the loss of strength due to wet conditions. The 3D simulation of the stability conditions concurs with our hypothesis. The landslide occurred in the Recent or Middle Holocene, certainly before at least 500 BC and possibly during a wet climate period. Today, it appears to be inactive. This study helps to understand the frequent slope instabilities all along the Iberian Range when facies Utrillas is present.

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

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

Induced Damping on Vibrating Circular Plates Submerged in still Fluid

When a structure vibrates immersed in a fluid it is known that the dynamic properties of the system are modified. The surrounding fluid will, in general, contribute to the inertia, the rigidity and the damping coefficient of the coupled fluid-structure system. For light structures, like spacecraft antennas, even when the fluid is air the contribution to the dynamic properties can be important. For not so light structures the ratio of the equivalent fluid/structure mass and rigidity can be very small and the fluid contribution could be neglected. For the ratio of equivalent fluid/structure damping both terms are of the same order and therefore the fluid contribution must be studied. The working life of the spacecraft structure would be on space and so without any surrounding fluid. The response of a spacecraft structure on its operational life would be attenuated by the structural damping alone but when the structure is dynamically tested on the earth the dynamic modal test is performed with the fluid surrounding it. The results thus are contaminated by the effects of the fluid. If the damping added by the fluid is of the same order as the structural damping the response of the structure in space can be quite different to the response predicted on earth. It is therefore desirable to have a method able to determine the amount of damping induced by the fluid and that should be subtracted of the total damping measured on the modal vibration test. In this work, a method for the determination of the effect of the surrounding fluid on the dynamic characteristics of a circular plate has been developed. The plate is assumed to vibrate harmonically with the vacuum modes and the generalized forces matrix due to the fluid is thus computed. For a compressible fluid this matrix is formed by complex numbers including terms of inertia, rigidity and damping. The matrix due to the fluid loading is determined by a boundary element method (BEM). The BEM used is of circular rings on the plate surface so the number of elements to obtain an accurate result is very low. The natural frequencies of the system are computed by an iteration procedure one by one and also the damping fluid contribution. Comparisons of the present method with various experimental data and other theories show the efficiency and accuracy of the method for any support condition of the plate.

Vortex Shedding Dynamics in Long Aspect-Ratio Aerodynamics Bodies

The focus of the current dissertation is to study qualitatively the underlying physics of vortex-shedding and wake dynamics in long aspect-ratio aerodynamics in incompressible viscous flow through the use of the KLE method. We carried out a long series of numerical experiments in the cases of flow around the cylinder at low Reynolds numbers. The study of flow at low Reynolds numbers provides an insight in the fluid physics and also plays a critical role when applying to stalled turbine rotors. Many of the conclusions about the qualitative nature of the physical mechanisms characterizing vortex formation, shedding and further interaction analyzed here at low Re could be extended to other Re regimes and help to understand the separation of the boundary layers in airfoils and other aerodynamic surfaces. In the long run, it aims to provide a better understanding of the complex multi-physics problems involving fluid-structure-control interaction through improved mathematical computational models of the multi-physics process. Besides the scientific conclusions produced, the research work on streamlined and bluff-body condition will also serve as a valuable guide for the future design of blade aerodynamics and the placement of wind turbines and hydrakinetic turbines, increasing the efficiency in the use of expensive workforce, supplies, and infrastructure. After the introductory section describing the main fields of application of wind power and hydrokinetic turbines, we describe the main features and theoretical background of the numerical method used here. Then, we present the analysis of the numerical experimentation results for the oscillatory regime right before the onset of vortex shedding for circular cylinders. We verified the wake length of the closed near-wake behind the cylinder and analysed the decay of the wake at the wake formation region, and then studied the St-Re relationship at the Reynolds numbers before the wake sheds compared to the experimental data. We found a theoretical model that describes the time evolution of the amplitude of fluctuations in the vorticity field on the twin vortex wake, which accurately matches the numerical results in terms of the frequency of the oscillation and rate of decay. We also proposed a model based on an analog circuit that is able to interpret the concerning flow by reducing the number of degrees of freedom. It follows the idea of the non-linear oscillator and resembles the dynamics mechanism of the closed near-wake with a common configured sine wave oscillator. This low-dimensional circuital model may also help to understand the underlying physical mechanisms, related to vorticity transport, that give origin to those oscillations.