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.

Fluid-rock interaction during late stages of the Alpine exhumation

Résumé de la thèseLa fracturation des roches au cours de phases compressives ou extensives est un souvent évoquée pour expliquer la circulation de fluide au sein des roches cristallines. Dans le cadre de cette thèse, la circulation des fluides lors de l'exhumation tardive des Alpes a été étudiée en utilisant deux approches différentes: analyses structurales de la déformation fragile d'une part et analyses géochimiques des roches et des minéraux (isotopes stables, datations U/Pb, thermochronologie (U-Th)/He) d'autre part. Cette approche combinée a permis de mieux comprendre l'interaction existante entre les fluides métamorphiques et les fluides météoriques, ainsi que leur interaction avec les roches encaissantes. Le travail a été effectué dans la zone Pennique du Valais suisse.La première partie était focalisée sur la déformation fragile, le but étant de définir les différents types de déformations existantes et de déterminer l'âge relatif des différentes familles de failles. Dans la région d'étude, quatre domaines ont été distingués. Chacun d'eux comportent deux types de structures fragiles, certaines sont minéralisées alors que d'autre non. Au sein de chaque domaine, la direction principale des structures minéralisées correspond à l'orientation des accidents tectoniques majeurs de la région (Aosta- Ranzola Line au Sud, Rhône Line au Nord et Simplon Fault Zone à l'Est), alors que les structures non- minéralisées montrent des orientations plus variables. Ainsi, le premier type de structure est interprété comme résultant d'une dislocation tectonique alors que le deuxième type de structure résulterait d'une dislocation gravitaire locale. Il n'est néanmoins pas possible de classer chronologiquement la formation de ces deux types de structure ni d'attribuer un âge relatif aux changements d'orientation des contraintes majeures.La deuxième étude a été effectuée dans la région de la zone de faille du Simplon. Dans cette zone, la composition isotopique des minéraux ayant cristallisé à l'intérieur des fractures tardives permet de distinguer différents types de circulation de fluide. Les valeurs δ180 du quartz de la roche encaissante ainsi que ceux des veines tardives du bloque inférieur de la faille sont comparables. Ces valeurs indiquent un rééquilibrage et un tamponnage isotopique des fluides tardifs au contact de la roche encaissante lors de la fracturation de cette dernière et de la cristallisation des veines tardives. La même situation est observée dans la partie nord du bloque supérieur ainsi que dans sa partie sud. Ceci n'est néanmoins pas le cas pour la partie centrale du bloque supérieur où les valeurs isotopiques des minéraux dans les veines tardives sont approximativement 3 %o plus basses (avec des valeurs extrêmes négatifs), indiquant une contribution d'eau météorique aux fluides circulant dans les veines. Ces données suggèrent qu'une infiltration d'eau météorique a pu avoir lieu dans le bloque supérieur, où la fracturation des roches est plus intensive car le déplacement relatif le long de la faille y fut plus important, et la température maximale du métamorphisme plus basse. La troisième contribution traite de la géo-thermochronologie de la zone de contact entre la klippe de la Dent Blanche et la nappe de Tsaté. De petits zircons euhédraux ont été trouvés dans un plan de faille minéralisé (parallèle à la Faille du Rhône, voir première partie de l'étude), riche en hématite et quartz, de la zone d'étude. Les analyses U/Pb donnent des âges radiométriques autour de 270 - 280 Ma aux zircons extraits de la minéralisation ainsi que ceux extraits de la roche encaissante, ce qui correspond à l'âge de la nappe de la Dent Blanche et non celui de la nappe du Tsaté qui est elle-même classiquement interprétée comme une ophiolite Jurassique de l'Océan Liguro-Piémontais. Ces données suggèrent que les zircons contenus dans la veine ont été hérités de la roche encaissante. Les résultats (U-Th)/He indiquent un âge de refroidissement différent pour la roche encaissante (25.5 ± 2.0 Ma) que celui de la minéralisation (17.7 ±1.4 Ma). Le thermomètre isotopique quartz-hématite indique une température d'équilibre, et donc de mise en place de la minéralisation, d'environ 170 °C, température très proche de la température de -180 °C de fermeture du zircon pour le système (U-Th)/He. Ceci suggère que l'âge de refroidissement des zircons de la minéralisation correspond aussi à l'âge de formation de la faille.Thesis abstractFluid circulation in fractured rocks is a common process in geology, and it is generally the consequence of faulting and fracturing during both tectonic compression and extension. This thesis is focused on fluid circulation during late stages of the Alpine exhumation. After a structural analysis of the late brittle deformation of the studied samples, several analytical methods (stable isotope investigations, U/Pb radiometric dating, (U-Th)/He thermochronology) have been applied to understand the interaction of metamorphic and meteoric fluids with one another as well as with the host rock. This thesis is articulated around three study directions. All studies were conducted in the Penninic Zone of the Valais, Switzerland. The first study deals with late, brittle deformation and focuses on the different deformation styles and on the relative age of the different families of fractures. In order to do this, late brittle structures observed in four different domains have been subdivided as a function of the existence (or not) and type of mineralization. Comparisons between mineralized and non-mineralized strike directions for all four domains show that mineralized structures follow the strike orientation of major tectonic movements indicated in the Penninic Zone of the Valais (Aosta-Ranzola Line to the S, Rhône Line to the Ν and Simplon Fault Zone to the E), whereas non-mineralized fractures have a more variable strike orientation. This difference could be interpreted as indicative of tectonic-related faulting (mineralized structures) vs. local, collapse-related faulting (non-mineralized fractures), but it is not strong enough to indicate a relative age of the late brittle structures, and/or a change in the orientation of the strain field in post-Miocene times. The second studied area is focused on the Simplon Fault Zone (SFZ). Stable isotope analyses of minerals filling these late fractures indicate that there are two different fluid circulation systems in the footwall and hanging wall of the SFZ. In the footwall, δ180 values of quartz from both the host rock and the late veins range from +10 %o to +12 %o. This is consistent with buffering of circulating fluids by the host rock during fracturing and vein precipitation. In the hanging wall, δΙ80 values for quartz crystals from the host rock and the late veins are similar in both the northern and southern parts of the detachment that are both affected by the same degree of metamorphism (greenschist to the Ν and amphibolite to the S). This is not the case in the central part of the SFZ, where there is a jump from amphibolite facies in the footwall to greenschist facies in the hanging wall. δ,80 values for quartz from the hanging wall late veins are approximately 3.0 %o lower (down to negative values in some cases) than the values observed in the footwall These data suggest that infiltration of meteoric water may have occurred in the most fractured parts of the hanging wall, where relative displacement on the SFZ was the greatest and the peak temperature lower. In the less fractured footwall the δ180 values reflect a host rock-buffered system.The third study is focused on geo-thermochronology at the contact between the Dent Blanche nappe and the Tsaté nappe where small, euhedral zircons were found in a hematite- and quartz-rich mineralization on a late normal fault plane parallel to the Rhône Line (see first part of the study). U/Pb analysis indicates that the zircons - both in the late mineralization and in the host rock - have absolute radiometric ages clustering around 270 - 280 Ma, which is the accepted age for intrusive rocks from the Austroalpine Dent Blanche units but not for the Tsaté nappe. The latter is classically interpreted as an ophiolitic remnant of the Jurassic Liguro-Piemontais Ocean. U/Pb analyses suggest that zircons in late mineralization are all inherited from the host rock; however, results of (U-Th)/He analyses indicate that cooling ages for the host rocks are different to the cooling ages for the zircons in late mineralization. Indeed, the calculated cooling age for the Arolla gneiss is 25.5 ± 2.0 Ma, whilst the cooling age for the associated mineralized fault plane is 17.7 ±1.4 Ma. Oxygen stable isotope fractionation between quartz and hematite in the same late mineralization corresponds to temperatures of about 170 °C. The proximity of the calculated emplacement temperature for the mineralization and the lower accepted closure temperature for zircon in the (U-Th)/He system (-180 °C) imply that the age of 17.7 ± 1.4 Ma can also be interpreted as the formation age of this late brittle fault.Résumé grand publicLa circulation des fluides dans les roches fracturées est typique de nombreux processus géologiques, et très souvent est la conséquence de la fracturation des roches. Cette thèse aborde la question de la circulation des fluides pendant les dernières phases du soulèvement des Alpes. Après une analyse structurale de la fracturation directement sur le terrain, plusieurs méthodes géochimiques ont été appliquées pour comprendre l'interaction entre les différents fluides circulants, et avec leur propre roche mère. L'étude, concentrée sur trois directions principales, a été conduite dans la zone Pennique du Valais suisse. La première partie traite de la déformation cassante dans le secteur cité. L'analyse détaillée des fractures a permis de les subdiviser en structures minéralisées et non-minéralisées, sur quatre domaines différents. La comparaison entre les directions des structures minéralisées et non-minéralisées a permis de montrer que les premières suivent l'orientation des accidents tectoniques majeurs de la région, alors que les structures non- minéralisées ont une orientation plus variable. Cette différence pourrait être interprétée comme indication d'une dislocation tectonique (structures minéralisées) contre une dislocation gravitaire locale (structures non-minéralisées), mais elle n'est pas assez forte pour indiquer un âge relatif des structures tardives et/ou un changement de l'orientation des contraintes après -20 Ma vers le présent.A partir de ces observations, la deuxième étude est concentrée dans la région de la faille du Simplon. Les analyses géochimiques sur les minéraux remplissant les structures tardives indiquent qu'il y a deux différents systèmes de circulation des fluides dans les deux parties (toit et mur) de la faille. Dans le mur, les valeurs isotopiques des minéraux cristallisés à partir d'un fluide tardif sont les mêmes de ceux de la roche mère, donc il y a eu rééquilibration chimique entre fluide et roche pendant la fracturation de cette dernière et la précipitation des minéraux. Dans le toit, les valeurs isotopiques dans la roche mère et dans les minéraux des veines tardives sont comparables dans les parties Ν et S de la faille, où les roches du toit et du mur ont atteint une température maximale - pendant phase prograde de la formation des Alpes - comparable. Au contraire, dans la partie centrale, où le mur a atteint des températures maximales plus élevées par rapport au toit, les valeurs géochimiques des minéralisations tardives du toit sont parfois plus basses que les valeurs observées dans le mur. Ces données suggèrent que l'infiltration de l'eau de surface aurait pu se produire dans la partie plus fracturée du toit, où le déplacement relatif le long de la faille était majeur et les températures maximales mineures. Au contraire, les données géochimiques du mur de la partie centrale indiquent un système isotopique équilibré par la roche mère.La troisième partie de ce travail se base sur l'étude géochimique intégrée des isotopes stables d'Oxygène et radioactifs du Plomb, Uranium, Thorium et Hélium, auprès d'une faille normale minéralisée et des roches de la région à cheval entre deux nappes, la nappe de la Dent Blanche et la nappe de Tsaté. Ici, des petits zircons ont été trouvés dans la minéralisation citée, riche en hématite et quartz. L'analyse radiométrique Uranium/Plomb a montré que les zircons dans la minéralisation et dans les roches autour ont des âges comparables (autour 280 Ma). Cela signifie que les zircons dans la minéralisation tardive ont été hérités de la roche mère pendant la fracturation et la circulation des fluides tardives. De l'autre coté, les résultats des analyses Uranium-Thorium/Hélium indiquent que les âges de refroidissement pour les roches mères sont différents comparés aux âges de refroidissement pour les zircons dans la minéralisation tardive: ces derniers sont plus jeunes d'environ 8 Ma (autour 25 Ma et autour 17 Ma respectivement). Les analyses des isotopes de l'oxygène sur quartz et hématite dans la même minéralisation donnent une température de mise en place de cette dernière d'environ 170° C. La température de fermeture du système chimique des zircons dans le système (Uranium-Thorium)/Hélium est d'environ 180 °C: la proximité de ces deux températures implique que l'âge de refroidissement de la minéralisation tardive peut également être interprété comme âge de formation de la faille.

Finite element modeling of fluid-rock interaction problems in pore-fluid saturated hydrothermal/sedimentary basins

In order to use the finite element method for solving fluid-rock interaction problems in pore-fluid saturated hydrothermal/sedimentary basins effectively and efficiently, we have presented, in this paper, the new concept and numerical algorithms to deal with the fundamental issues associated with the fluid-rock interaction problems. These fundamental issues are often overlooked by some purely numerical modelers. (1) Since the fluid-rock interaction problem involves heterogeneous chemical reactions between reactive aqueous chemical species in the pore-fluid and solid minerals in the rock masses, it is necessary to develop the new concept of the generalized concentration of a solid mineral, so that two types of reactive mass transport equations, namely, the conventional mass transport equation for the aqueous chemical species in the pore-fluid and the degenerated mass transport equation for the solid minerals in the rock mass, can be solved simultaneously in computation. (2) Since the reaction area between the pore-fluid and mineral surfaces is basically a function of the generalized concentration of the solid mineral, there is a definite need to appropriately consider the dependence of the dissolution rate of a dissolving mineral on its generalized concentration in the numerical analysis. (3) Considering the direct consequence of the porosity evolution with time in the transient analysis of fluid-rock interaction problems; we have proposed the term splitting algorithm and the concept of the equivalent source/sink terms in mass transport equations so that the problem of variable mesh Peclet number and Courant number has been successfully converted into the problem of constant mesh Peclet and Courant numbers. The numerical results from an application example have demonstrated the usefulness of the proposed concepts and the robustness of the proposed numerical algorithms in dealing with fluid-rock interaction problems in pore-fluid saturated hydrothermal/sedimentary basins. (C) 2001 Elsevier Science B.V. All rights reserved.

Modeling the inelastic behavior of heat exchangers accounting for fluid-structure interactions

Magdeburg, Univ., Fak. für Maschinenbau, Diss., 2013

Stability and asymptotic analysis of a fluid-particle interaction model

We are interested in coupled microscopic/macroscopic models describing the evolution of particles dispersed in a fluid. The system consists in a Vlasov-Fokker-Planck equation to describe the microscopic motion of the particles coupled to the Euler equations for a compressible fluid. We investigate dissipative quantities, equilibria and their stability properties and the role of external forces. We also study some asymptotic problems, their equilibria and stability and the derivation of macroscopic two-phase models.

Numerical Investigations on Soil Structure Interaction of Multistorey Frames

Frames are the most widely used structural system for multistorey buildings. A building frame is a three dimensional discrete structure consisting of a number of high rise bays in two directions at right angles to each other in the vertical plane. Multistorey frames are a three dimensional lattice structure which are statically indeterminate. Frames sustain gravity loads and resist lateral forces acting on it. India lies at the north westem end of the Indo-Australian tectonic plate and is identified as an active tectonic area. Under horizontal shaking of the ground, horizontal inertial forces are generated at the floor levels of a multistorey frame. These lateral inertia forces are transferred by the floor slab to the beams, subsequently to the columns and finally to the soil through the foundation system. There are many parameters that affect the response of a structure to ground excitations such as, shape, size and geometry of the structure, type of foundation, soil characteristics etc. The Soil Structure Interaction (SS1) effects refer to the influence of the supporting soil medium on the behavior of the structure when it is subjected to different types of loads. Interaction between the structure and its supporting foundation and soil, which is a complete system, has been modeled with finite elements. Numerical investigations have been carried out on a four bay, twelve storeyed regular multistorey frame considering depth of fixity at ground level, at characteristic depth of pile and at full depth. Soil structure interaction effects have been studied by considering two models for soil viz., discrete and continuum. Linear static analysis has been conducted to study the interaction effects under static load. Free vibration analysis and further shock spectrum analysis has been conducted to study the interaction effects under time dependent loads. The study has been extended to four types of soil viz., laterite, sand, alluvium and layered.The structural responses evaluated in the finite element analysis are bending moment, shear force and axial force for columns, and bending moment and shear force for beams. These responses increase with increase in the founding depth; however these responses show minimal increase beyond the characteristic length of pile. When the soil structure interaction effects are incorporated in the analysis, the aforesaid responses of the frame increases upto the characteristic depth and decreases when the frame has been analysed for the full depth. It has been observed that shock spectrum analysis gives wide variation of responses in the frame compared to linear elastic analysis. Both increase and decrease in responses have been observed in the interior storeys. The good congruence shown by the two finite element models viz., discrete and continuum in linear static analysis has been absent in shock spectrum analysis.

Acoustic Structure Interaction For Floating Structure Under Moving Load

The present thesis concentrates largely on sound radiation from floating structure due to moving load

Structure-activity relationship of mastoparan analogs: effects of the number and positioning of Lys residues on secondary structure, interaction with membrane-mimetic systems and biological activity

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Advances in the study of soil-structure interaction effects on the dynamic response of piled structures

[EN]This Ph. D. thesis presents a simple and stable procedure for the estimation of periods and dampings of pile shear buildings taking soil-structure interaction into account. The coupled-system response is obtained by using a substructuring model. A boundary element-finite element coupling formulation is used to compute impedances and kinematic interaction factors of the pile group configurations under investigation. The proposed procedure is applied to perform parametric analyses to determine the influence of the main parameters of soil-structure interaction problems on the dynamic response of the superstructure. The scope of this thesis also encompasses the study of foundations including battered piles.

Dynamic structure-soil-structure interaction between nearby piled buildings under seismic excitation by BEM-FEM model

[EN]The dynamic throug-soil interaction between nearby pile supported structures in a viscoelastic half-space, under incident S and Rayleigh waves, is numerically studied. To this end, a three-dimensional viscoelastic BEM-FEM formulation for the dynamic analysis of piles and pile groups in the frequency domain is used, where soil is modelled by BEM and piles are simulated by one-dimensional finite elements as Bernouilli beams.

Structure-soil-structure interaction effects on the dynamic response of piled structures under obliquely-incident seismic shear waves.

[EN] This work studies the structure-soil-structure interaction (SSSI) effects on the dynamic response of nearby piled structures under obliquely-incident shear waves. For this purpose, a three-dimensional, frequency-domain, coupled boundary element-finite (BEM-FEM) model is used to analyse the response of configuration of three buildings aligned parallel to the horizontal component of the wave propagation direction.

Effects of soil-structure interaction on the dynamic properties and seismic response of piled structures.

[EN]This paper presents a simple and stable procedure for the estimation of periods and dampings of piled shear buildings taking soil-structure interaction into account. A substructuring methodology that incluedes the three-dimensional character of the foundations is used.

Aeroelastic stability of structures: flutter analysis using Computational Fluid Dynamics

Thanks to the increasing slenderness and lightness allowed by new construction techniques and materials, the effects of wind on structures became in the last decades a research field of great importance in Civil Engineering. Thanks to the advances in computers power, the numerical simulation of wind tunnel tests has became a valid complementary activity and an attractive alternative for the future. Due to its flexibility, during the last years, the computational approach gained importance with respect to the traditional experimental investigation. However, still today, the computational approach to fluid-structure interaction problems is not as widely adopted as it could be expected. The main reason for this lies in the difficulties encountered in the numerical simulation of the turbulent, unsteady flow conditions generally encountered around bluff bodies. This thesis aims at providing a guide to the numerical simulation of bridge deck aerodynamic and aeroelastic behaviour describing in detail the simulation strategies and setting guidelines useful for the interpretation of the results.

Analysis and mathematical modeling of wave-structure interaction

The aim of this thesis, included within the THESEUS project, is the development of a mathematical model 2DV two-phase, based on the existing code IH-2VOF developed by the University of Cantabria, able to represent together the overtopping phenomenon and the sediment transport. Several numerical simulations were carried out in order to analyze the flow characteristics on a dike crest. The results show that the seaward/landward slope does not affect the evolution of the flow depth and velocity over the dike crest whereas the most important parameter is the relative submergence. Wave heights decrease and flow velocities increase while waves travel over the crest. In particular, by increasing the submergence, the wave height decay and the increase of the velocity are less marked. Besides, an appropriate curve able to fit the variation of the wave height/velocity over the dike crest were found. Both for the wave height and for the wave velocity different fitting coefficients were determined on the basis of the submergence and of the significant wave height. An equation describing the trend of the dimensionless coefficient c_h for the wave height was derived. These conclusions could be taken into consideration for the design criteria and the upgrade of the structures. In the second part of the thesis, new equations for the representation of the sediment transport in the IH-2VOF model were introduced in order to represent beach erosion while waves run-up and overtop the sea banks during storms. The new model allows to calculate sediment fluxes in the water column together with the sediment concentration. Moreover it is possible to model the bed profile evolution. Different tests were performed under low-intensity regular waves with an homogeneous layer of sand on the bottom of a channel in order to analyze the erosion-deposition patterns and verify the model results.

Soil-structure interaction during tunnelling in urban area: observations and 3D numerical modelling

This work illustrates a soil-tunnel-structure interaction study performed by an integrated,geotechnical and structural,approach based on 3D finite element analyses and validated against experimental observations.The study aims at analysing the response of reinforced concrete framed buildings on discrete foundations in interaction with metro lines.It refers to the case of the twin tunnels of the Milan (Italy) metro line 5,recently built in coarse grained materials using EPB machines,for which subsidence measurements collected along ground and building sections during tunnelling were available.Settlements measured under freefield conditions are firstly back interpreted using Gaussian empirical predictions. Then,the in situ measurements’ analysis is extended to include the evolving response of a 9 storey reinforced concrete building while being undercrossed by the metro line.In the finite element study,the soil mechanical behaviour is described using an advanced constitutive model. This latter,when combined with a proper simulation of the excavation process, proves to realistically reproduce the subsidence profiles under free field conditions and to capture the interaction phenomena occurring between the twin tunnels during the excavation. Furthermore, when the numerical model is extended to include the building, schematised in a detailed manner, the results are in good agreement with the monitoring data for different stages of the twin tunnelling. Thus, they indirectly confirm the satisfactory performance of the adopted numerical approach which also allows a direct evaluation of the structural response as an outcome of the analysis. Further analyses are also carried out modelling the building with different levels of detail. The results highlight that, in this case, the simplified approach based on the equivalent plate schematisation is inadequate to capture the real tunnelling induced displacement field. The overall behaviour of the system proves to be mainly influenced by the buried portion of the building which plays an essential role in the interaction mechanism, due to its high stiffness.