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.

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.

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.

Investigation of Human-Structure Interaction Through Experimental and Analytical Studies

Vibration serviceability is a widely recognized design criterion for assembly-type structures, such as stadiums, that are likely subjected to rhythmic human-induced excitation. Human-induced excitation of a structure occurs from the movement of the occupants such as walking, running, jumping, or dancing. Vibration serviceability is based on the level of comfort that people have with the vibrations of a structure. Current design guidance uses the natural frequency of the structure to assess vibration serviceability. However, a phenomenon known as human-structure interaction suggests that there is a dynamic interaction between the structure and passive occupants, altering the natural frequency of the system. Human-structure interaction is dependent on many factors, including the dynamic properties of the structure, posture of the occupants, and relative size of the crowd. It is unknown if the shift in natural frequency due to humanstructure interaction is significant enough to warrant consideration in the design process. This study explores the interface of both structural and crowd characteristics through experimental testing to determine if human-structure interaction should be considered because of its potential impact on serviceability assessment. An experimental test structure that represents the dynamic properties of a cantilevered stadium structure was designed and constructed. Experimental modal analysis was implemented to determine the dynamic properties of the empty test structure and when occupied with up to seven people arranged in different locations and postures. Comparisons of the dynamic properties were made between the empty and occupied testing configurations and analytical results from the use of a dynamic crowd model recommended from the Joint Working Group of Europe. Data trends lead to the development of a refined dynamic crowd model. This dynamic model can be used in conjunction with a finite element model of the test structure to estimate the dynamic influence due to human-structure interaction due to occupants standing with straight knees. In the future, the crowd model will be refined and can aid in assessing the dynamic properties of in-service stadium structures.

Experimental Evaluation of the Influence of Human-Structure Interaction for Vibration Serviceability

The effects of human-structure interaction on the dynamic performance of occupied structures have long been observed. The inclusion of the effects of human-structure interaction is important to ensure that the dynamic response of a structure is not overestimated. Previous observations, both in service and in the laboratory, have yielded results indicating that the effects are dependent on the natural frequency of the structure, the posture of the occupants, and the mass ratio of the occupants to the structure. These results are noteworthy, but are limited in their application,because the data are sparse and are only pertinent to a specific set of characteristics identified in a given study. To examine these characteristics simultaneously and consistently, an experimental test structure was designed with variable properties to replicate a variety of configurations within a controlled setting focusing on the effects of passive occupants. Experimental modal analysis techniques were employed to both the empty and occupied conditions of the structure and the dynamic properties associated with each condition were compared. Results similar to previous investigations were observed, including both an increase and a decrease in natural frequency of the occupied structure with respect to the empty structure, as well as the identification of a second mode of vibration. The damping of the combined system was higher for all configurations. Overall, this study provides a broad data set representing a wide array of configurations. The experimental results of this study were used to assess current recommendations for the dynamic properties of a crowd to analytically predict the effects of human-structure interaction. The experimental results were used to select a set of properties for passive, standing occupants and develop a new model that can more accurately represent the behavior of the human-structure system as experimentally measured in this study.

Solution structure and interaction of cupiennin 1a, a spider venom peptide, with phospholipid bilayers

The solution structure of cupiennin 1a, a 35 residue, basic antibacterial peptide isolated from the venom of the spider Cupiennius salei, has been determined by nuclear magnetic resonance (NMR) spectroscopy. The peptide was found to adopt a helix−hinge−helix structure in a membrane mimicking solvent. The hinge may play a role in allowing the amphipathic N-terminal helix and polar C-terminal helix to orient independently upon membrane binding, in order to achieve maximal antibacterial efficacy. Solid-state 31P and 2H NMR was used to further study the effects of cupiennin 1a on the dynamic properties of lipid membranes, using zwitterionic chain deuterated dimyristoylphosphatidylcholine (d54-DMPC) and anionic dimyristoylphosphatidylglycerol (DMPG) multilamellar vesicles. In d54-DMPC alone, cupiennin 1a caused a decrease in the 31P chemical shift anisotropy, indicating some interaction with the lipid head groups, and a decrease in order over the entire acyl chain. In contrast, for the mixed (d54-DMPC/DMPG) lipid system cupiennin 1a appeared to induce lateral separation of the two lipids as evidenced by the 31P spectra, in which the peptide preferentially interacted with DMPG. Little effect was observed on the deuterated acyl chain order parameters in the d54-DMPC/DMPG model membranes. Furthermore, 31P NMR relaxation measurements confirmed a differential effect on the lipid motions depending upon the membrane composition. Therefore, subtle differences are likely in the mechanism by which cupiennin 1a causes membrane lysis in either prokaryotic or eukaryotic cells, and may explain the specific spectrum of activity.

Structural response and dynamics of fluid-structure-control interaction in wind turbine blades

Reducing the uncertainties related to blade dynamics by the improvement of the quality of numerical simulations of the fluid structure interaction process is a key for a breakthrough in wind-turbine technology. A fundamental step in that direction is the implementation of aeroelastic models capable of capturing the complex features of innovative prototype blades, so they can be tested at realistic full-scale conditions with a reasonable computational cost. We make use of a code based on a combination of two advanced numerical models implemented in a parallel HPC supercomputer platform: First, a model of the structural response of heterogeneous composite blades, based on a variation of the dimensional reduction technique proposed by Hodges and Yu. This technique has the capacity of reducing the geometrical complexity of the blade section into a stiffness matrix for an equivalent beam. The reduced 1-D strain energy is equivalent to the actual 3-D strain energy in an asymptotic sense, allowing accurate modeling of the blade structure as a 1-D finite-element problem. This substantially reduces the computational effort required to model the structural dynamics at each time step. Second, a novel aerodynamic model based on an advanced implementation of the BEM(Blade ElementMomentum) Theory; where all velocities and forces are re-projected through orthogonal matrices into the instantaneous deformed configuration to fully include the effects of large displacements and rotation of the airfoil sections into the computation of aerodynamic forces. This allows the aerodynamic model to take into account the effects of the complex flexo-torsional deformation that can be captured by the more sophisticated structural model mentioned above. In this thesis we have successfully developed a powerful computational tool for the aeroelastic analysis of wind-turbine blades. Due to the particular features mentioned above in terms of a full representation of the combined modes of deformation of the blade as a complex structural part and their effects on the aerodynamic loads, it constitutes a substantial advancement ahead the state-of-the-art aeroelastic models currently available, like the FAST-Aerodyn suite. In this thesis, we also include the results of several experiments on the NREL-5MW blade, which is widely accepted today as a benchmark blade, together with some modifications intended to explore the capacities of the new code in terms of capturing features on blade-dynamic behavior, which are normally overlooked by the existing aeroelastic models.

The interaction of generic structure and interpersonal relations in two-party e-chat discourse

On the basis of a corpus of e-chat IRC exchanges (approximately 10,000 words in total) between Greek- and English-speaking speakers, the paper establishes a typical generic structure for two-party IRC exchanges, by focusing on how participants are oriented towards an ideal schema of phases and acts, as well as on how their interpersonal concerns contribute to the shaping of this schema. It is found that IRC interlocutors are primarily concerned with establishing contact with each other, while the (ideational) development of topic seems to be a less pressing need. The signaling of interpersonal relations is pervasive throughout e-chat discourse, as seen both in the range of devices developed and the two free elements of the generic schema, that is conversation play and channel check. It is also found that the accomplishment of the generic schema in each IRC exchange crucially depends on the acts of negotiation performed by the initiator and the responder.

NMR structure of a complex between the VirB9/VirB7 interaction domains of the pKM101 type IV secretion system.

Type IV secretion (T4S) systems translocate DNA and protein effectors through the double membrane of Gram-negative bacteria. The paradigmatic T4S system in Agrobacterium tumefaciens is assembled from 11 VirB subunits and VirD4. Two subunits, VirB9 and VirB7, form an important stabilizing complex in the outer membrane. We describe here the NMR structure of a complex between the C-terminal domain of the VirB9 homolog TraO (TraO(CT)), bound to VirB7-like TraN from plasmid pKM101. TraO(CT) forms a beta-sandwich around which TraN winds. Structure-based mutations in VirB7 and VirB9 of A. tumefaciens show that the heterodimer interface is conserved. Opposite this interface, the TraO structure shows a protruding three-stranded beta-appendage, and here, we supply evidence that the corresponding region of VirB9 of A. tumefaciens inserts in the membrane and protrudes extracellularly. This complex structure elucidates the molecular basis for the interaction between two essential components of a T4S system.

Interaction of Plectin with Keratins 5 and 14: Dependence on Several Plectin Domains and Keratin Quaternary Structure.

Plectin, a cytolinker of the plakin family, anchors the intermediate filament (IF) network formed by keratins 5 and 14 (K5/K14) to hemidesmosomes, junctional adhesion complexes in basal keratinocytes. Genetic alterations of these proteins cause epidermolysis bullosa simplex (EBS) characterized by disturbed cytoarchitecture and cell fragility. The mechanisms through which mutations located after the documented plectin IF-binding site, composed of the plakin-repeat domain (PRD) B5 and the linker, as well as mutations in K5 or K14, lead to EBS remain unclear. We investigated the interaction of plectin C terminus, encompassing four domains, the PRD B5, the linker, the PRD C, and the C extremity, with K5/K14 using different approaches, including a rapid and sensitive fluorescent protein-binding assay, based on enhanced green fluorescent protein-tagged proteins (FluoBACE). Our results demonstrate that all four plectin C-terminal domains contribute to its association with K5/K14 and act synergistically to ensure efficient IF binding. The plectin C terminus predominantly interacted with the K5/K14 coil 1 domain and bound more extensively to K5/K14 filaments compared with monomeric keratins or IF assembly intermediates. These findings indicate a multimodular association of plectin with K5/K14 filaments and give insights into the molecular basis of EBS associated with pathogenic mutations in plectin, K5, or K14 genes.Journal of Investigative Dermatology advance online publication, 10 July 2014; doi:10.1038/jid.2014.255.