Development and validation of Fluid Structure Interaction Methods for an Oscillating Wave Surge Converter

The term fatigue loads on the Oyster Oscillating Wave Surge Converter (OWSC) is used to describe hydrostatic loads due to water surface elevation with quasi-static changes of state. Therefore a procedure to implement hydrostatic pressure distributions into finite element analysis of the structure is desired. Currently available experimental methods enable one to measure time variant water surface elevation at discrete locations either on or around the body of the scale model during tank tests. This paper discusses the development of a finite element analysis procedure to implement time variant, spatially distributed hydrostatic pressure derived from discretely measured water surface elevation. The developed method can process differently resolved (temporal and spatial) input data and approximate the elevation over the flap faces with user defined properties. The structural loads, namely the forces and moments on the body can then be investigated by post processing the numerical results. This method offers the possibility to process surface elevation or hydrostatic pressure data from computational fluid dynamics simulations and can thus be seen as a first step to a fluid-structure interaction model.

Development of a Vehicle-Bridge-Soil Dynamic Interaction Model for Scour Damage Modelling

Damage detection in bridges using vibration-based methods is an area of growing research interest. Improved assessment
methodologies combined with state-of-the-art sensor technology are rapidly making these approaches applicable for real-world
structures. Applying these techniques to the detection and monitoring of scour around bridge foundations has remained
challenging; however this area has gained attraction in recent years. Several authors have investigated a range of methods but
there is still significant work required to achieve a rounded and widely applicable methodology to detect and monitor scour.This
paper presents a novel Vehicle-Bridge-Soil Dynamic Interaction (VBSDI) model which can be used to simulate the effect of scour
on an integral bridge. The model outputs dynamic signals which can be analysed to determine modal parameters and the variation
of these parameters with respect to scour can be examined.The key novelty of this model is that it is the first numerical model for
simulating scour that combines a realistic vehicle loadingmodel with a robust foundation soil responsemodel.This paper provides a
description of the model development and explains the mathematical theory underlying themodel. Finally a case study application
of the model using typical bridge, soil, and vehicle properties is provided.

Acoustic Structure Interaction For Floating Structure Under Moving Load

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

The analysis of ranked observations of soil structure using indicator geostatistics

Structure is an important physical feature of the soil that is associated with water movement, the soil atmosphere, microorganism activity and nutrient uptake. A soil without any obvious organisation of its components is known as apedal and this state can have marked effects on several soil processes. Accurate maps of topsoil and subsoil structure are desirable for a wide range of models that aim to predict erosion, solute transport, or flow of water through the soil. Also such maps would be useful to precision farmers when deciding how to apply nutrients and pesticides in a site-specific way, and to target subsoiling and soil structure stabilization procedures. Typically, soil structure is inferred from bulk density or penetrometer resistance measurements and more recently from soil resistivity and conductivity surveys. To measure the former is both time-consuming and costly, whereas observations made by the latter methods can be made automatically and swiftly using a vehicle-mounted penetrometer or resistivity and conductivity sensors. The results of each of these methods, however, are affected by other soil properties, in particular moisture content at the time of sampling, texture, and the presence of stones. Traditional methods of observing soil structure identify the type of ped and its degree of development. Methods of ranking such observations from good to poor for different soil textures have been developed. Indicator variograms can be computed for each category or rank of structure and these can be summed to give the sum of indicator variograms (SIV). Observations of the topsoil and subsoil structure were made at four field sites where the soil had developed on different parent materials. The observations were ranked by four methods and indicator and the sum of indicator variograms were computed and modelled for each method of ranking. The individual indicators were then kriged with the parameters of the appropriate indicator variogram model to map the probability of encountering soil with the structure represented by that indicator. The model parameters of the SIVs for each ranking system were used with the data to krige the soil structure classes, and the results are compared with those for the individual indicators. The relations between maps of soil structure and selected wavebands from aerial photographs are examined as basis for planning surveys of soil structure. (C) 2007 Elsevier B.V. All rights reserved.

Estimation and evaluation of dynamic properties as indicators of changes on soil structure in sugarcane fields of São Paulo State - Brazil

The indiscriminate management and use of soils without moisture control has changed the structure of it due to the increment of the traffic by agricultural machines through the years, causing in consequence, a soil compaction and yield reduction in the areas of intensive traffic. The purpose of this work was to estimate and to evaluate the performance of preconsolidation pressure of the soil and shear stress as indicators of changes on soil structure in fields cropped with sugarcane, as well as the impact of management processes in an Eutrorthox soil structure located in São Paulo State. The experimental field was located in Piracicaba's rural area (São Paulo State, Brazil) and has been cropped with sugarcane, in the second harvest cycle. The soil was classified by Empresa Brasileira de Pesquisa Agropecuária (EMBRAPA) [Empresa Brasileira de Pesquisa Agropecuária (EMBRAPA), 1999. Centro Nacional de Pesquisa de Solos. Sistema Brasileiro de Classificao de Solos, Empresa Brasileira de Pesquisa Agropecuária (EMBRAPA), Brasilia, 412 pp.] as an Eutrorthox. Undisturbed samples were collected and georeferenced in a grid of 60 m x 60 m from two depths: 0-0.10 m (superficial layer - SL) and in the layer of greatest mechanical resistance (LGMR), previously identified by cone index (CI). The investigated variables were pressure preconsolidation (sigma(p)), apparent cohesion (c) and internal friction angle (phi). The conclusions from the results were that the SLSC was predicted satisfactorily from up as a function of soil moisture; thus, decisions about machinery size and loading (contact pressures) can be taken. Apparent cohesion (c), internal friction angle (phi) and the Coulomb equation were significantly altered by traffic intensity. The sigma(p), c and phi maps were shown to be important tools to localize and visualize soil compaction and mechanical resistance zones. They constitute a valuable resource to evaluate the traffic impact in areas cropped with sugarcane in State of São Paulo, Brazil. (C) 2008 Elsevier B.V. All rights reserved.

Geosynthetic-Reinforced Soil Structure Built to Restore Original Geometry of a 12 m High Slope After Failure Using Local Soil

This paper presents a study case in which a geosynthetic-reinforced soil (GRS) structure was used to rebuild a 12 m high slope after its failure. The failed slope is located between the parking lot of a private company and a public school. Due to surrounding structures restrictions, this project required a solution with rapidity in execution. In addition, as a requirement established by its owner, this structure should recover the original geometry of the slope. Besides the importance regarding surrounding constructions, an interesting aspect of this study case relies on the versatility of geosynthetic materials. A woven geotextile was used as reinforcement. Five other geosynthetic materials were used in this study case. Facing comprised a geocell filled with local soil cover and grass mats, resulting in a green facing. A geonet was used to hold the grass mats in place before grass roots development. Regarding the drainage system, geocomposite drains and geopipes were installed to drain subsurface water. A nonwoven geotextile was used as filter in drainage trenches, which were placed near the structure toe. Additionally to the GRS structure, the lower portion of the slope was reinforced with soil nailing technique. The face of the nailed soil portion was covered with sandbags and shotcrete. It emphasizes the flexibility of GRS structures regarding their application with other technical options in Geotechnical Engineering. The economic aspect of this study case also deserves attention. It did not require soil transportation and other design and construction steps, e.g. concrete structures design and construction.

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)

A 2D BEM-FEM approach for time harmonic fluid-structure interaction analysis of thin elastic bodies.

[EN]This paper deals with two-dimensional time harmonic fluid-structure interaction problems when the fluid is at rest, and the elastic bodies have small thicknesses. A BEM-FEM numerical approach is used, where the BEM is applied to the fluid, and the structural FEM is applied to the thin elastic bodies.

Fluid-structure interaction by a coupled lattice Boltzmann-finite element approach

In this thesis, a strategy to model the behavior of fluids and their interaction with deformable bodies is proposed. The fluid domain is modeled by using the lattice Boltzmann method, thus analyzing the fluid dynamics by a mesoscopic point of view. It has been proved that the solution provided by this method is equivalent to solve the Navier-Stokes equations for an incompressible flow with a second-order accuracy. Slender elastic structures idealized through beam finite elements are used. Large displacements are accounted for by using the corotational formulation. Structural dynamics is computed by using the Time Discontinuous Galerkin method. Therefore, two different solution procedures are used, one for the fluid domain and the other for the structural part, respectively. These two solvers need to communicate and to transfer each other several information, i.e. stresses, velocities, displacements. In order to guarantee a continuous, effective, and mutual exchange of information, a coupling strategy, consisting of three different algorithms, has been developed and numerically tested. In particular, the effectiveness of the three algorithms is shown in terms of interface energy artificially produced by the approximate fulfilling of compatibility and equilibrium conditions at the fluid-structure interface. The proposed coupled approach is used in order to solve different fluid-structure interaction problems, i.e. cantilever beams immersed in a viscous fluid, the impact of the hull of the ship on the marine free-surface, blood flow in a deformable vessels, and even flapping wings simulating the take-off of a butterfly. The good results achieved in each application highlight the effectiveness of the proposed methodology and of the C++ developed software to successfully approach several two-dimensional fluid-structure interaction problems.

Multilevel Domain Decomposition Algorithms for Monolithic Fluid-Structure Interaction Problems with Application to Haemodynamics

Finite element techniques for solving the problem of fluid-structure interaction of an elastic solid material in a laminar incompressible viscous flow are described. The mathematical problem consists of the Navier-Stokes equations in the Arbitrary Lagrangian-Eulerian formulation coupled with a non-linear structure model, considering the problem as one continuum. The coupling between the structure and the fluid is enforced inside a monolithic framework which computes simultaneously for the fluid and the structure unknowns within a unique solver. We used the well-known Crouzeix-Raviart finite element pair for discretization in space and the method of lines for discretization in time. A stability result using the Backward-Euler time-stepping scheme for both fluid and solid part and the finite element method for the space discretization has been proved. The resulting linear system has been solved by multilevel domain decomposition techniques. Our strategy is to solve several local subproblems over subdomain patches using the Schur-complement or GMRES smoother within a multigrid iterative solver. For validation and evaluation of the accuracy of the proposed methodology, we present corresponding results for a set of two FSI benchmark configurations which describe the self-induced elastic deformation of a beam attached to a cylinder in a laminar channel flow, allowing stationary as well as periodically oscillating deformations, and for a benchmark proposed by COMSOL multiphysics where a narrow vertical structure attached to the bottom wall of a channel bends under the force due to both viscous drag and pressure. Then, as an example of fluid-structure interaction in biomedical problems, we considered the academic numerical test which consists in simulating the pressure wave propagation through a straight compliant vessel. All the tests show the applicability and the numerical efficiency of our approach to both two-dimensional and three-dimensional problems.

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.

Numerical modelling of fluid-structure interaction with application in hemodynamics

The thesis deals with numerical algorithms for fluid-structure interaction problems with application in blood flow modelling. It starts with a short introduction on the mathematical description of incompressible viscous flow with non-Newtonian viscosity and a moving linear viscoelastic structure. The mathematical model consists of the generalized Navier-Stokes equation used for the description of fluid flow and the generalized string model for structure movement. The arbitrary Lagrangian-Eulerian approach is used in order to take into account moving computational domain. A part of the thesis is devoted to the discussion on the non-Newtonian behaviour of shear-thinning fluids, which is in our case blood, and derivation of two non-Newtonian models frequently used in the blood flow modelling. Further we give a brief overview on recent fluid-structure interaction schemes with discussion about the difficulties arising in numerical modelling of blood flow. Our main contribution lies in numerical and experimental study of a new loosely-coupled partitioned scheme called the kinematic splitting fluid-structure interaction algorithm. We present stability analysis for a coupled problem of non-Newtonian shear-dependent fluids in moving domains with viscoelastic boundaries. Here, we assume both, the nonlinearity in convective as well is diffusive term. We analyse the convergence of proposed numerical scheme for a simplified fluid model of the Oseen type. Moreover, we present series of experiments including numerical error analysis, comparison of hemodynamic parameters for the Newtonian and non-Newtonian fluids and comparison of several physiologically relevant computational geometries in terms of wall displacement and wall shear stress. Numerical analysis and extensive experimental study for several standard geometries confirm reliability and accuracy of the proposed kinematic splitting scheme in order to approximate fluid-structure interaction problems.

Neural network modelling of the wave-structure interaction processes

This thesis presents a new Artificial Neural Network (ANN) able to predict at once the main parameters representative of the wave-structure interaction processes, i.e. the wave overtopping discharge, the wave transmission coefficient and the wave reflection coefficient. The new ANN has been specifically developed in order to provide managers and scientists with a tool that can be efficiently used for design purposes. The development of this ANN started with the preparation of a new extended and homogeneous database that collects all the available tests reporting at least one of the three parameters, for a total amount of 16’165 data. The variety of structure types and wave attack conditions in the database includes smooth, rock and armour unit slopes, berm breakwaters, vertical walls, low crested structures, oblique wave attacks. Some of the existing ANNs were compared and improved, leading to the selection of a final ANN, whose architecture was optimized through an in-depth sensitivity analysis to the training parameters of the ANN. Each of the selected 15 input parameters represents a physical aspect of the wave-structure interaction process, describing the wave attack (wave steepness and obliquity, breaking and shoaling factors), the structure geometry (submergence, straight or non-straight slope, with or without berm or toe, presence or not of a crown wall), or the structure type (smooth or covered by an armour layer, with permeable or impermeable core). The advanced ANN here proposed provides accurate predictions for all the three parameters, and demonstrates to overcome the limits imposed by the traditional formulae and approach adopted so far by some of the existing ANNs. The possibility to adopt just one model to obtain a handy and accurate evaluation of the overall performance of a coastal or harbor structure represents the most important and exportable result of the work.

Assessment Of Dynamic Properties Of A Crowd Model For Human-Structure Interaction Modeling

As lightweight and slender structural elements are more frequently used in the design, large scale structures become more flexible and susceptible to excessive vibrations. To ensure the functionality of the structure, dynamic properties of the occupied structure need to be estimated during the design phase. Traditional analysis method models occupants simply as an additional mass; however, research has shown that human occupants could be better modeled as an additional degree-of- freedom. In the United Kingdom, active and passive crowd models are proposed by the Joint Working Group as a result of a series of analytical and experimental research. It is expected that the crowd models would yield a more accurate estimation to the dynamic response of the occupied structure. However, experimental testing recently conducted through a graduate student project at Bucknell University indicated that the proposed passive crowd model might be inaccurate in representing the impact on the structure from the occupants. The objective of this study is to provide an assessment of the validity of the crowd models proposed by JWG through comparing the dynamic properties obtained from experimental testing data and analytical modeling results. The experimental data used in this study was collected by Firman in 2010. The analytical results were obtained by performing a time-history analysis on a finite element model of the occupied structure. The crowd models were created based on the recommendations from the JWG combined with the physical properties of the occupants during the experimental study. During this study, SAP2000 was used to create the finite element models and to implement the analysis; Matlab and ME¿scope were used to obtain the dynamic properties of the structure through processing the time-history analysis results from SAP2000. The result of this study indicates that the active crowd model could quite accurately represent the impact on the structure from occupants standing with bent knees while the passive crowd model could not properly simulate the dynamic response of the structure when occupants were standing straight or sitting on the structure. Future work related to this study involves improving the passive crowd model and evaluating the crowd models with full-scale structure models and operating data.