A New Method for Modeling Free Surface Flows and Fluid-structure Interaction with Ocean Applications

The computational modeling of ocean waves and ocean-faring devices poses numerous challenges. Among these are the need to stably and accurately represent both the fluid-fluid interface between water and air as well as the fluid-structure interfaces arising between solid devices and one or more fluids. As techniques are developed to stably and accurately balance the interactions between fluid and structural solvers at these boundaries, a similarly pressing challenge is the development of algorithms that are massively scalable and capable of performing large-scale three-dimensional simulations on reasonable time scales. This dissertation introduces two separate methods for approaching this problem, with the first focusing on the development of sophisticated fluid-fluid interface representations and the second focusing primarily on scalability and extensibility to higher-order methods.

We begin by introducing the narrow-band gradient-augmented level set method (GALSM) for incompressible multiphase Navier-Stokes flow. This is the first use of the high-order GALSM for a fluid flow application, and its reliability and accuracy in modeling ocean environments is tested extensively. The method demonstrates numerous advantages over the traditional level set method, among these a heightened conservation of fluid volume and the representation of subgrid structures.

Next, we present a finite-volume algorithm for solving the incompressible Euler equations in two and three dimensions in the presence of a flow-driven free surface and a dynamic rigid body. In this development, the chief concerns are efficiency, scalability, and extensibility (to higher-order and truly conservative methods). These priorities informed a number of important choices: The air phase is substituted by a pressure boundary condition in order to greatly reduce the size of the computational domain, a cut-cell finite-volume approach is chosen in order to minimize fluid volume loss and open the door to higher-order methods, and adaptive mesh refinement (AMR) is employed to focus computational effort and make large-scale 3D simulations possible. This algorithm is shown to produce robust and accurate results that are well-suited for the study of ocean waves and the development of wave energy conversion (WEC) devices.

A finite volume approach to dynamic fluid structure interaction

Abstract not available

Effect of soil structure on temporal and spatial dynamics of bacteria

Soil is a complex heterogeneous system comprising of highly variable and dynamic micro-habitats that have significant impacts on the growth and activity of resident microbiota. A question addressed in this research is how soil structure affects the temporal dynamics and spatial distribution of bacteria. Using repacked microcosms, the effect of bulk-density, aggregate sizes and water content on growth and distribution of introduced Pseudomonas fluorescens and Bacillus subtilis bacteria was determined. Soil bulk-density and aggregate sizes were altered to manipulate the characteristics of the pore volume where bacteria reside and through which distribution of solutes and nutrients is controlled. X-ray CT was used to characterise the pore geometry of repacked soil microcosms. Soil porosity, connectivity and soil-pore interface area declined with increasing bulk-density. In samples that differ in pore geometry, its effect on growth and extent of spread of introduced bacteria was investigated. The growth rate of bacteria reduced with increasing bulk-density, consistent with a significant difference in pore geometry. To measure the ability of bacteria to spread thorough soil, placement experiments were developed. Bacteria were capable of spreading several cm’s through soil. The extent of spread of bacteria was faster and further in soil with larger and better connected pore volumes. To study the spatial distribution in detail, a methodology was developed where a combination of X-ray microtopography, to characterize the soil structure, and fluorescence microscopy, to visualize and quantify bacteria in soil sections was used. The influence of pore characteristics on distribution of bacteria was analysed at macro- and microscales. Soil porosity, connectivity and soil-pore interface influenced bacterial distribution only at the macroscale. The method developed was applied to investigate the effect of soil pore characteristics on the extent of spread of bacteria introduced locally towards a C source in soil. Soil-pore interface influenced spread of bacteria and colonization, therefore higher bacterial densities were found in soil with higher pore volumes. Therefore the results in this showed that pore geometry affects the growth and spread of bacteria in soil. The method developed showed showed how thin sectioning technique can be combined with 3D X-ray CT to visualize bacterial colonization of a 3D pore volume. This novel combination of methods is a significant step towards a full mechanistic understanding of microbial dynamics in structured soils.

Development of a Coupling Model for Fluid-Structure Interaction using the Mesh-free Finite Element Method and the Lattice Boltzmann Method

In the presented thesis work, the meshfree method with distance fields was coupled with the lattice Boltzmann method to obtain solutions of fluid-structure interaction problems. The thesis work involved development and implementation of numerical algorithms, data structure, and software. Numerical and computational properties of the coupling algorithm combining the meshfree method with distance fields and the lattice Boltzmann method were investigated. Convergence and accuracy of the methodology was validated by analytical solutions. The research was focused on fluid-structure interaction solutions in complex, mesh-resistant domains as both the lattice Boltzmann method and the meshfree method with distance fields are particularly adept in these situations. Furthermore, the fluid solution provided by the lattice Boltzmann method is massively scalable, allowing extensive use of cutting edge parallel computing resources to accelerate this phase of the solution process. The meshfree method with distance fields allows for exact satisfaction of boundary conditions making it possible to exactly capture the effects of the fluid field on the solid structure.

Determination of stresses around a cylindrical single pile caused by horizontal movements of soils with the finitive element method

This study was carried out with the aim of modeling in 2D, in plain strain, the movement of a soft cohesive soil around a pile, in order to enable the determination of stresses resulting along the pile, per unit length. The problem in study fits into the large deformations problem and can be due to landslide, be close of depth excavations, to be near of zones where big loads are applied in the soil, etc. In this study is used an constitutive Elasto-Plastic model with the failure criterion of Mohr-Coulomb to model the soil behavior. The analysis is developed considering the soil in undrained conditions. To the modeling is used the finite element program PLAXIS, which use the Updated Lagrangian - Finite Element Method (UL-FEM). In this work, special attention is given to the soil-pile interaction, where is presented with some detail the formulation of the interface elements and some studies for a better understand of his behavior. It is developed a 2-D model that simulates the effect of depth allowing the study of his influence in the stress distribution around the pile. The results obtained give an important base about how behaves the movement of the soil around a pile, about how work the finite element program PLAXIS and how is the stress distribution around the pile. The analysis demonstrate that the soil-structure interaction modeled with the UL-FEM and interface elements is more appropriate to small deformations problems.

Pontes integrais: estado da arte

Dissertação para obtenção do grau de mestre em Engenharia Civil

Dynamic stiffness of the piled foundation of a onshore wind turbine

One of today's biggest concerns is the increase of energetic needs, especially in the developed countries. Among various clean energies, wind energy is one of the technologies that assume greater importance on the sustainable development of humanity. Despite wind turbines had been developed and studied over the years, there are phenomena that haven't been yet fully understood. This work studies the soil-structure interaction that occurs on a wind turbine's foundation composed by a group of piles that is under dynamic loads caused by wind. This problem assumes special importance when the foundation is implemented on locations where safety criteria are very demanding, like the case of a foundation mounted on a dike. To the phenomenon of interaction between two piles and the soil between them it's given the name of pile-soil-pile interaction. It is known that such behavior is frequency dependent, and therefore, on this work evaluation of relevant frequencies for the intended analysis is held. During the development of this thesis, two methods were selected in order to assess pile-soil-pile interaction, being one of analytical nature and the other of numerical origin. The analytical solution was recently developed and its called Generalized pile-soil-pile theory, while for the numerical method the commercial nite element software PLAXIS 3D was used. A study of applicability of the numerical method is also done comparing the given solution by the nite element methods with a rigorous solution widely accepted by the majority of the authors.

Design Methodology for Corrugated Metal Pipe Tiedowns, HR-362, Phase II, 1995

This investigation is the final phase of a three part study whose overall objectives were to determine if a restraining force is required to prevent inlet uplift failures in corrugated metal pipe (CMP) installations, and to develop a procedure for calculating the required force when restraint is required. In the initial phase of the study (HR-306), the extent of the uplift problem in Iowa was determined and the forces acting on a CMP were quantified. In the second phase of the study (HR- 332), laboratory and field tests were conducted. Laboratory tests measured the longitudinal stiffness ofCMP and a full scale field test on a 3.05 m (10 ft) diameter CMP with 0.612 m (2 ft) of cover determined the soil-structure interaction in response to uplift forces. Reported herein are the tasks that were completed in the final phase of the study. In this phase, a buried 2.44 m (8 ft) CMP was tested with and without end-restraint and with various configurations of soil at the inlet end of the pipe. A total of four different soil configurations were tested; in all tests the soil cover was constant at 0.61 m (2 ft). Data from these tests were used to verify the finite element analysis model (FEA) that was developed in this phase of the research. Both experiments and analyses indicate that the primary soil contribution to uplift resistance occurs in the foreslope and that depth of soil cover does not affect the required tiedown force. Using the FEA, design charts were developed with which engineers can determine for a given situation if restraint force is required to prevent an uplift failure. If an engineer determines restraint is needed, the design charts provide the magnitude of the required force. The design charts are applicable to six gages of CMP for four flow conditions and two types of soil.

Investigation of High Density Polyethylene Pipe for Highway Applications, HR-373, Phase 1, 1996

In the past, culvert pipes were made only of corrugated metal or reinforced concrete. In recent years, several manufacturers have made pipe of lightweight plastic - for example, high density polyethylene (HDPE) - which is considered to be viscoelastic in its structural behavior. It appears that there are several highway applications in which HDPE pipe would be an economically favorable alternative. However, the newness of plastic pipe requires the evaluation of its performance, integrity, and durability; A review of the Iowa Department of Transportation Standard Specifications for Highway and Bridge Construction reveals limited information on the use of plastic pipe for state projects. The objective of this study was to review and evaluate the use of HDPE pipe in roadway applications. Structural performance, soil-structure interaction, and the sensitivity of the pipe to installation was investigated. Comprehensive computerized literature searches were undertaken to define the state-of-the-art in the design and use of HDPE pipe in highway applications. A questionnaire was developed and sent to all Iowa county engineers to learn of their use of HDPE pipe. Responses indicated that the majority of county engineers were aware of the product but were not confident in its ability to perform as well as conventional materials. Counties currently using HDPE pipe in general only use it in driveway crossings. Originally, we intended to survey states as to their usage of HDPE pipe. However, a few weeks after initiation of the project, it was learned that the Tennessee DOT was in the process of making a similar survey of state DOT's. Results of the Tennessee survey of states have been obtained and included in this report. In an effort to develop more confidence in the pipe's performance parameters, this research included laboratory tests to determine the ring and flexural stiffness of HDPE pipe provided by various manufacturers. Parallel plate tests verified all specimens were in compliance with ASTM specifications. Flexural testing revealed that pipe profile had a significant effect on the longitudinal stiffness and that strength could not be accurately predicted on the basis of diameter alone. Realizing that the soil around a buried HDPE pipe contributes to the pipe stiffness, the research team completed a limited series of tests on buried 3 ft-diameter HDPE pipe. The tests simulated the effects of truck wheel loads above the pipe and were conducted with two feet of cover. These tests indicated that the type and quality of backfill significantly influences the performance of HDPE pipe. The tests revealed that the soil envelope does significantly affect the performance of HDPE pipe in situ, and after a certain point, no additional strength is realized by increasing the quality of the backfill.

Be-Fe coupled method for the analysis of 3D elastodynamic problems in the time domain

By coupling the Boundary Element Method (BEM) and the Finite Element Method (FEM) an algorithm that combines the advantages of both numerical processes is developed. The main aim of the work concerns the time domain analysis of general three-dimensional wave propagation problems in elastic media. In addition, mathematical and numerical aspects of the related BE-, FE- and BE/FE-formulations are discussed. The coupling algorithm allows investigations of elastodynamic problems with a BE- and a FE-subdomain. In order to observe the performance of the coupling algorithm two problems are solved and their results compared to other numerical solutions.

Análise do comportamento de estacas do tipo broca escavada com trado mecânico, solicitadas por esforços transversais

Pós-graduação em Engenharia Civil - FEIS

Numerical combination for nonlinear analysis of structures coupled to layered soils

This paper presents an alternative coupling strategy between the Boundary Element Method (BEM) and the Finite Element Method (FEM) in order to create a computational code for the analysis of geometrical nonlinear 2D frames coupled to layered soils. The soil is modeled via BEM, considering multiple inclusions and internal load lines, through an alternative formulation to eliminate traction variables on subregions interfaces. A total Lagrangean formulation based on positions is adopted for the consideration of the geometric nonlinear behavior of frame structures with exact kinematics. The numerical coupling is performed by an algebraic strategy that extracts and condenses the equivalent soil's stiffness matrix and contact forces to be introduced into the frame structures hessian matrix and internal force vector, respectively. The formulation covers the analysis of shallow foundation structures and piles in any direction. Furthermore, the piles can pass through different layers. Numerical examples are shown in order to illustrate and confirm the accuracy and applicability of the proposed technique.

Numerical Model for the Dynamic Analysis of Pile Foundations

Programa de doctorado: Sistemas Inteligentes y Aplicaciones Numéricas en Ingeniería Instituto Universitario (SIANI)

Influence of pile inclination angle on the dynamic properties and seismic response of piled structures.

[EN] This paper aims to contribute to clarify whether the use of battered piles has a positive or negative influence on the dynamic response of deep foundations and superstructures. For this purpose, the dynamic response of slender and non-slender structures supported on several configurations of 2X2 and 3X3 pile groups including battered elements is obtained through a procedure based on a substructuring model whick takes soil-structure interaction into account.

Tendencias actuales en la construcción sismorresistente para edificios en hormigón armado

Este estudio aborda la recopilación de nuevas tendencias del diseño sismorresistente, enfocándose en la técnica del aislamiento de base, por ser la más efectiva, difundida y utilizada; y el análisis de las ventajas que puede tener una edificación que aplica dicha técnica, desde el punto de vista estructural y económico. Se elige la tipología más frecuente o común de edificios de hormigón armado propensos a ser aislados, que en este caso es un hospital, cuyo modelo empotrado se somete a varias normas sismorresistentes comparando principalmente fuerzas de cortante basal, y considerando la interacción suelo-estructura; para asistir a este cálculo se desarrolla un programa de elementos viga de 6 gdl por nodo en código Matlab. El modelo aislado incluye el análisis de tres combinaciones de tipos de aisladores HDR, LPR y FPS, alternando modelos lineales simplificados de 1 y 3 gdl por piso, evaluando diferencias de respuestas de la estructura, y procediendo a la elección de la combinación que de resultados más convenientes; para la modelación no lineal de cada sistema de aislamiento se utiliza el método explícito de diferencias centrales. Finalmente, se realiza un análisis comparativo de daños esperados en el caso de la ocurrencia del sismo de diseño, utilizando el método rápido y tomando como referencia el desplazamiento espectral del último piso; llegando a dar conclusiones y recomendaciones para el uso de sistemas de aislamiento. This study addresses the collection of new seismic design trends, focusing on base isolation technique, as the most effective and widely used, and the analysis of the advantages in buildings that apply this technique, from the structurally and economically point of view. Choosing the most common types of concrete buildings likely to be isolated, which in this case is a hospital, the fix model is subjected to various seismic codes mainly comparing base shear forces, and considering the soil-structure interaction; for this calculation attend a program of bars 6 dof per node is made in Matlab code. The isolated model includes analysis of three types of isolators combinations HDR, LPR and FPS, alternating simplified linear model of 1 and 3 dof per floor, evaluating differences in the response of the structure, and proceeding to the choice of the combination of results more convenient; for modeling nonlinear each insulation system, the explicit central difference method is used. Finally, a comparative analysis of expected damage in the case of the design earthquake, using a fast combined method and by reference to the spectral displacement of the top floor; reaching conclusions and give recommendations for the use of insulation systems.