<p>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.</p><p>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.</p><p> </p><p>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.</p> / Dissertation

Le domaine de l'interaction fluide-structure regroupe l'étude de tous les phénomènes présentant le couplage du mouvement d'une structure avec celui d'un fluide. La gamme des phénomènes étudiés est très étendue, allant de l'étude de cylindres vibrants dans des écoulements comme c'est le cas dans l'industrie nucléaire, à des structures vibrantes dans des écoulements turbulents, en passant par des phénomènes de surface libre dans des réservoirs. Cependant, la complexité des phénomènes étudiés se répercute par des coûts de calculs prohibitifs, ce qui nous amène à rechercher des modèles réduits dont le temps de calcul serait plus réaliste. Dans cette thèse, on va présenter les différents modèles d'interaction fluide-structure et on va mettre en avant le modèle adopté dans notre étude. La réduction du modèle ainsi que l'optimisation des structures vont être introduites dans un contexte de couplage. En introduisant les incertitudes, l'étude fiabiliste de même qu'une approche d'optimisation basée fiabilité vont être proposées. Les différentes méthodologies adoptées vont être validées numériquement et comparées expérimentalement / The domain of the fluid-structure interaction includes the study of all phenomena presenting the coupling of the motion of a structure with the one of a fluid. The range of the phenomena being studied is very extensive, going from the study of vibrating cylinders in the flow as is the case in the nuclear industry, to vibrating structures in turbulent flows, through the free surface phenomena in reservoirs. However, the complexity of the phenomena studied is reflected by the cost of the prohibitive calculations, which leads us to look for models with the computation time would be more realistic. In this thesis, we will present different models of fluid-structure interaction and we will put forward the model adopted in our study. Reducing the model as well as the optimization of the structures will be introduced into a coupling setting. By introducing uncertainties, the reliability study as well as an optimization based reliability approach will be proposed. The different methodologies adopted will be validated numerically and experimentally compared

Les canopées sont des assemblées compactes de plantes dont l'étude concerne de nombreuses problématiques environnementales. Des applications technologiques sont également envisageables. Les précédents travaux se sont principalement focalisés sur les écoulements au-dessus des canopées, considérées comme des rugosités de fond. La présence d'un point d'inflexion dans le profil de vitesses dans le fluide autorise le développement d'instabilités de type couche de mélange à l'interface. De plus, la prise en compte de la flexibilité des plantes complique le problème, car leur forme est modifiée par le courant pour réduire la traînée exercée sur elles --c'est le phénomène de reconfiguration-- mais elles ont également une dynamique propre qui peut éventuellement influencer l'écoulement. La démarche envisagée dans cette thèse est essentiellement expérimentale. Elle cherche à comprendre la réaction des tiges à différents types d'écoulements, afin d'étudier les effets collectifs inhérents à la canopée, et d'identifier les mécanismes communs qui en sont à l'origine. On utilise pour cela des tiges modèles très simples dans un canal étroit, ce qui assure une configuration quasi 2D et facilite les observations. Dans un premier temps, on étudie la réaction statique de la canopée à un écoulement établi. L'effet de la densité est très clair tant que les plantes sont assez proches, sinon elles se comportent comme si elles étaient seules. Ensuite, la canopée est soumise à un écoulement oscillant (houle), et, de la même façon, on étudie la différence de réaction entre une tige seule et une tige incluse dans une canopée. La troisième partie s'intéresse à la dynamique de la canopée soumise à un écoulement unidirectionnel, permettant le développement d'instabilités au sommet de la canopée. Le régime de grandes ondulations cohérentes de la canopée, apparenté au ‘monami' de la littérature, est caractérisé / Vegetation in rivers is often considered as a source of water resistance which slows down the water conveyance. It is also one of the main component for river equilibrium, insofar as it prevents body erosion by providing bed stabilization, it plays a vital role during floods. It is crucial for sediment transport, water quality and also shelter to provide the necessary habitat for the biodiversity of aquatic species. It is then useful to understand the mechanical behaviour of aquatic canopies resulting from the interaction between vegetation and a water flow. From land-use planning to river management, such a knowledge would also shed light upon plant biomechanics and improve bio-inspired engineering.Traditionally, studies on aquatic vegetation explored its influence on flow properties, like velocity distribution, wake dynamics, turbulence, water conveyance and sediment transport by considering it simply as a rigid or flexible roughness element.This thesis is an experimental work which aims at understanding how a model canopy reacts to a water flow depending on the canopy geometry and the flow conditions. Three different series of experiments are performed.First, the effect of density on the canopy reconfiguration and the corresponding drag reduction is investigated. The drag acting on the canopy, and also on individual sheets, is systematically measured. A strong sheltering effect exists as long as the spacing is smaller than a critical value depending on the sheet width.Then, the canopy is submitted to a wave flow to test its sensibility to a determined frequency. Each stem is found to act like a forced oscillator with a strong resonance at natural frequency (modified by canopy density).Finally, a parallel free flow allows mixing layer instabilities to develop above the canopy. Different behaviour are observed and characterized, especially the large coherent waving called 'monami'

Oil and gas exploitation in deep waters has become more than just a profit business to be a daily necessity, since the world energy matrix is based on fossil components. Risers are offshore structures that are intimately linked with oil and gas exploitation and those are subjected to a great variety of effects in field, e.g., marine currents, Vortex Induced Vibration (VIV), heave motion caused by gravitational waves, non-linear contact with the sea floor, and many others. Riser dynamics is essentially non-linear and experimental tests in real scale are almost impossible due to a great variety of control parameters acting concomitantly. Small-scale models are a better experimental approach. Nevertheless, there are many structural and hydrodynamical parameters to be evaluated. Considering only vertical risers in the present work, Galerkin\'s modal decomposition is used in order to reduce the dynamics of a vertical flexible cylinder to a few linear modes in which the majority of energy and information are contained. From the modal analysis, added mass and structural parameters damping of a vertical flexible cylinder using data obtained from free-decay tests performed both in water and in air are evaluated. Finally, a modal Mathieu-Hill oscillator with non-linear damping is constructed and, based on aStrutt diagram, modal stability under parametric resonance is discussed. / Exploração de óleo e gás em bacias de águas profundas tem-se tornado mais do que apenas uma economia lucrativa, para ser uma necessidade diária, já que a matriz energética mundial está baseada em componentes fósseis. Risers são estruturas offshore ligadas intimamente com a exploração de óleo e gás e essas estão sujeitas a uma grande variedade de efeitos na operação, e.g., correntes marítimas, Vibrações Induzidas por Vórtices (VIV), movimento de heave causado por ondas gravitacionais, contato não-linear com o solo marinho, entre outros. Dinâmica de risers é essencialmente não-linear e testes experimentais em escala real são praticamente impossíveis devido a uma enorme variedade de parâmetros de controle agindo concomitantemente. Modelos em escala reduzida são uma abordagem experimental mais conveniente. Não obstante, há muitos parâmetros estruturais e hidrodinâmicos a serem determinados. Considerando apenas risers verticais no trabalho presente, a decomposição modal de Galerkin é usada a fim de reduzir a dinâmica de um cilindro fléxivel vertical a alguns modos lineares em que a maior parte da energia e informação estão contidos. A partir da análise modal, parâmetros de massa adicional e amortecimento estrutural de um cilindro flexível vertical são obtidos usando testes de decaimento livre conduzidos na água e no ar. Finalmente, um oscilador modal de Mathieu-Hill com amortecimento não-linear é proposto e, baseado em um diagrama de Strutt, a estabilidade modal sob excitação de ressonânica paramétrica é discutida.

Esta tese apresenta o desenvolvimento e aplicação de modelos de turbulência, transição laminar-turbulenta e de interações fluido-estrutura ao escoamento externo em cilindro rígido estacionário e em vibrações induzidas por vórtices. Tais desenvolvimentos foram realizados no código ReFRESCO, baseado em técnicas de dinâmica de fluidos computacional (CFD). Realizou-se um estudo quanto ao desempenho do modelo k- SST em extensa faixa de números de Reynolds, segundo o qual se identificaram as deficiências de modelagem para este escoamento. A modelagem adaptativa das escalas (SAS) e o modelo de transição por correlações locais (LCTM), ambos combinados ao SST, melhoraram a aderência aos resultados experimentais para este escoamento, em uma contribuição original deste trabalho. A aplicação de técnicas de verificação e validação possibilitou a estimação de incertezas e erros para os modelos e números de Reynolds e também de identificada como outra contribuição deste trabalho. A combinação da modelagem em SST, SAS e LCTM com movimentos impostos de realizada para números de Reynolds moderados, diferentes frequências e amplitudes de vibração, algo que poucas publicações abordam em detalhes. Com relação aos movimentos livres, este trabalho traz contribuições com a aplicação dos modelos SST e SAS ao estudo de vibrações induzidas por vórtices em dois graus de liberdade, baixa razão de massa e números de Reynolds moderados, mais altos do que normalmente observados na literatura. Por fim, a investigação da importância relativa de efeitos da turbulência aos casos de movimentos livres e impostos, com relação ao caso de cilindro estacionário, comprovou a conjetura formulada na parte inicial deste trabalho, no que tange à escolha do modelo de turbulência em determinadas aplicações. Tal escolha mostrou-se menos decisiva no caso do cilindro em movimento imposto e ainda menos nos movimentos livres, em comparação ao caso estacionário, uma vez que a resposta em movimentos do corpo filtra grande parte dos efeitos turbulentos de ordem superior. Esta observação mostra-se relevante, uma vez que pode permitir simplificações na modelagem e aplicação de ferramentas de CFD em uma classe importante de projetos de engenharia. / This thesis presents the development, implementation and application of turbulence and laminar-turbulent transition models and fuid-structure capabilities to address the vortexshedding and vortex-induced vibrations of a rigid cylinder. These numerical developments have been carried out in the computational fuid dynamics (CFD) code ReFRESCO. In the current work, an investigation of the performance of the turbulence modeling with k- SST in a broad range of Reynolds numbers is carried out identifying its modeling deficiencies for this fow. The implementation and systematic application of the scale adaptive simulations (SAS) and the local correlation transition model (LCTM), both combined with the SST, have improved the agreement with experimental results for the cylinder ow, in a novel contribution of this work. The application of verification and validation technique has allowed the estimation of numerical errors and uncertainties for the diferent models. That is also identified as a contribution of this thesis. The combination of SST modeling with imposed motions is carried out as well as with the SAS and LCTM for moderate Reynolds numbers, diferent vibration frequencies and amplitudes, which is considered novel, as few publications address this issue in extent. Regarding the free-moving cylinder capabilities, the present work brings contributions with the application of SST and SASSST with free-moving cylinder for the study of VIV of two degrees of-freedom, low mass ratio and moderate Reynolds numbers, higher than commonly seen in the literature. Finally, the investigation of the relative importance of turbulence effects on the freemoving cylinder and the imposed-motions case, with respect to the fixed case is carried out. A natural conjecture that has been raised early on this work and proved correct is that, for engineering applications, the choice of turbulence modeling strategy is less decisive when the cylinder is moving with prescribed motion and even less stringent, for free motions as the body response filters most of the higher order turbulence effects. That is a relevant observation as it might allow modeling simplifications and the application of CFD tools to a range of engineering problems.

Les animaux volants et flottants ont développé des façons efficaces de produire l'écoulement de fluide qui génère les forces désirées pour leur locomotion. Cette thèse est placée dans ce contexte interdisciplinaire et utilise des simulations numériques pour étudier ces problèmes d'interaction fluides-structure, et les applique au vol des insectes et à la nage des poissons. Basée sur les travaux existants sur les obstacles mobiles rigides, une méthode numérique a été développée, permettant également la simulation des obstacles déformables et fournissant une polyvalence et précision accrues dans le cas des obstacles rigides. Nous appliquons cette méthode d'abord aux insectes avec des ailes rigides, où le corps et d'autres détails, tels que les pattes et les antennes, peuvent être inclus. Après la présentation de tests de validation détaillée, nous procédons à l'étude d'un modèle de bourdon dans un écoulement turbulent pleinement développé. Nos simulations montrent que les perturbations turbulentes affectent les insectes volants d'une manière différente de celle des avions aux ailes fixées et conçues par l'humain. Dans le cas de ces derniers, des perturbations en amont peuvent déclencher des transitions dans la couche limite, tandis que les premiers ne présentent pas de changements systématiques dans les forces aérodynamiques. Nous concluons que les insectes se trouvent plutôt confrontés à des problèmes de contrôle dans un environnement turbulent qu'à une détérioration de la production de force. Lors de l‘étape suivante, nous concevons un modèle solide, basé sur une équation de barre monodimensionnelle, et nous passons à la simulation des systèmes couplés fluide–structure. / Flying and swimming animals have developed efficient ways to produce the fluid flow that generates the desired forces for their locomotion. These bio-inspired problems couple fluid dynamics and solid mechanics with complex geometries and kinematics. The present thesis is placed in this interdisciplinary context and uses numerical simulations to study these fluid--structure interaction problems with applications in insect flight and swimming fish. Based on existing work on rigid moving obstacles, using an efficient Fourier discretization, a numerical method has been developed, which allows the simulation of flexible, deforming obstacles as well, and provides enhanced versatility and accuracy in the case of rigid obstacles. The method relies on the volume penalization method and the fluid discretization is still based on a Fourier discretization. We first apply this method to insects with rigid wings, where the body and other details, such as the legs and antennae, can be included. After presenting detailed validation tests, we proceed to studying a bumblebee model in fully developed turbulent flow. Our simulations show that turbulent perturbations affect flapping insects in a different way than human-designed fixed-wing aircrafts. While in the latter, upstream perturbations can cause transitions in the boundary layer, the former do not present systematical changes in aerodynamic forces. We conclude that insects rather face control problems in a turbulent environment than a deterioration in force production. In the next step, we design a solid model, based on a one--dimensional beam equation, and simulate coupled fluid--solid systems.

Multidisciplinary design optimization (MDO) refers to the process of designing systems characterized by the interaction of multiple interconnected disciplines. High-fidelity MDO usually requires large computational resources due to the computational cost of achieving multidisciplinary consistent solutions by coupling high-fidelity physics-based solvers. Gradient-based minimization algorithms are generally applied to find local minima, due to their efficiency in solving problems with a large number of design variables. This represents a limitation to performing global MDO and integrating black-box type analysis tools, usually not providing gradient information. The latter issues generally inhibit a wide use of MDO in complex industrial applications. An architecture named multi-criterion adaptive sampling MDO (MCAS-MDO) is presented in the current research for complex simulation-based applications. This research aims at building a global derivative-free optimization tool able to employ high-fidelity/expensive black-box solvers for the analysis of the disciplines. MCAS-MDO is a surrogate-based architecture featuring a variable level of coupling among the disciplines and is driven by a multi-criterion adaptive sampling (MCAS) assessing coupling and sampling uncertainties. MCAS uses the dynamic radial basis function surrogate model to identify the optimal solution and explore the design space through parallel infill of new solutions. The MCAS-MDO is tested versus a global derivative-free multidisciplinary feasible (MDF) approach, which solves fully-coupled multidisciplinary analyses, for two analytical test problems. Evaluation metrics include number of function evaluations required to achieve the optimal solution and sample distribution. The MCAS-MDO outperforms the MDF showing a faster convergence by clustering refined function evaluations in the optimum region. The architecture is applied to a steady fluid-structure interaction (FSI) problem, namely the design of a tapered three-dimensional carbon fiber-reinforced plastic hydrofoil for minimum drag. The objective is the design of shape and composite material layout subject to hydrodynamic, structural, and geometrical constraints. Experimental data are available for the original configuration of the hydrofoil and allow validating the FSI analysis, which is performed coupling computational fluid dynamics, solving the Reynolds averaged Navier-Stokes equations, and finite elements, solving the structural equation of elastic motion. Hydrofoil forces, tip displacement, and tip twist are evaluated for several materials providing qualitative agreement with the experiments and confirming the need for the two-way versus one-way coupling approach in case of significantly compliant structures. The free-form deformation method is applied to generate shape modifications of the hydrofoil geometry. To reduce the global computational expense of the optimization, a design space assessment and dimensionality reduction based on the Karhunen–Loève expansion (KLE) is performed off-line, i.e. without the need for high-fidelity simulations. It provides with a selection of design variables for the problem at hand through basis rotation and re-parametrization. By using the KLE, an efficient design space is identified for the current problem and the number of design variables is reduced by 92%. A sensitivity analysis is performed prior to the optimization to assess the variability associated with the shape design variables and the composite material design variable, i.e. the fiber orientation. These simulations are used to initialize the surrogate model for the optimization, which is carried out for two models: one in aluminum and one in composite material. The optimized designs are assessed by comparison with the original models through evaluation of the flow field, pressure distribution on the body, and deformation under the hydrodynamic load. The drag of the aluminum and composite material hydrofoils is reduced by 4 and 11%, respectively, increasing the hydrodynamic efficiency by 4 and 7%. The optimized designs are obtained by evaluating approximately 100 designs. The quality of the results indicates that global derivative-free MDO of complex engineering applications using expensive black-box solvers can be achieved at a feasible computational cost by minimizing the design space dimensionality and performing an intelligent sampling to train the surrogate-based optimization.

Le développement des systèmes houlomoteurs ainsi que la gestion du littoral reposent sur une bonne compréhension des mécanismes liés aux interactions houle-structure. Dans cette thèse, nous nous intéressons à l'étude d'un champ de structures flexibles soumises à des ondes de surface, en vue de développer un système qui puisse à la fois atténuer les vagues et absorber l'énergie qui leur est associée de manière efficace. Les résultats présentés se basent autour d'expériences réalisées dans des installations de petite échelle, dans lesquelles la disposition spatiale des objets flexibles est le principal paramètre étudié. Dans un premier temps, nous caractérisons notre champ modèle afin d'évaluer l'influence de divers paramètres (configuration, flexibilité, fréquences des vagues) sur la distribution de l'énergie dans le système. Sur la base de ces résultats, nous développons ensuite un modèle d'interférences permettant de décrire les observations globales du système à partir de paramètres locaux connus, associés à une portion unitaire du champ. Ce modèle nous sert ensuite d'outil pour l'exploration d'une multitude de configurations spatiales, afin de déterminer le choix optimal vis-à-vis de l'atténuation et de l'absorption des vagues incidentes. Enfin, une campagne de mesures supplémentaire est utilisée afin d'expliquer les résultats obtenus avec le modèle et d'identifier les principes sous-jacents à cette optimisation / Understanding the mechanisms involved in wave-structure interactions is of high interest for the development of efficient wave energy harvesters as well as for coastal management. In this thesis, we study the interactions of surface waves with a model array of slender flexible structures, in view of developing an efficient system for both attenuating and harvesting wave energy. The presented results are based around experimental investigations, by means of small scale facilities, in which the spatial arrangement of the flexible objects is the key parameter of study. The model array is first characterised by evaluating the role played by various parameters (configuration, flexibility, wave frequency) on the energy distribution in our system. Following these first observations, an interference model is then developed in order to describe the observed global effects of the array on both the wave field and the blade dynamics, based on known local parameters of a unit item of the array. This model then serves as a tool for exploring many possible array configurations, in order to determine the optimal choice regarding both the attenuation and the absorption of the imposed waves. A final experimental study is presented, in which the key results from the interference model are evaluated and the underlying principles of array optimisation are identified

The ultimate goal of this research is to develop the capability of direct numerical simulation of a flow containing numerous rigid finite size particles. In order to reach this goal, we have implemented the smoothed profile method (SPM) in the University of Iowa in-house solver, pELAFINT3D and overcame several challenges related to the method. This includes a proposed formulation for the interface thickness and many validations and comparisons with experimental data as well as with a second-order accurate sharp interface method. As one of the issues related to low-density particles is the instability, SPM has been improved by developing to a fully implicit scheme. Moreover, use of higher order integration formulation and implementation of Euler parameters have been shown to be helpful in stabilization of the calculations. To preserve the efficiency when the number of the particles increases, local mesh refinement is shown to be a very effective tool. A revised version of SPM that has only one projection step is proposed to improve the efficiency of the method. A comprehensive efficiency study is performed and it has been shown that the new method is less expensive in problems with high added mass effect when strongly coupled fluid-structure interaction schemes are used. Moreover, the code is massively parallelized using MPI and PETSc libraries. The parallelization includes I/O, operations leading to construction of the linear solver as well as the solver itself. Simulation of a particle laden flow involves particles collisions. Two novel collision models are suggested which are able to avoid particle overlapping for arbitrary shape particles. The methods are efficient as they are not involved with extra grid refinement related to implementing lubrication forces. The issue of handling continuously changing number of particles in a particle laden flow is solved by implementation of a linked list data structure for the particles. By studying a flow over a constricted region we showed that the platelets’ activation is more likely to happen for the particles that pass from the middle of the upper bump region because those particles will have longer exposure time to the high shear flow behind the bump. PDF contour of particles’ presence show the more concentrated presence of the particles near the bump. Moreover, the interaction of RBCs and platelets pushes the platelets toward the wall after the bottom wall.

Patient-specific multi-physics simulations have the potential to improve the diagnosis, treatment, and scientific inquiry of heart valve dynamics. It has been shown that the flow characteristics within the left ventricle are important to correctly capture the aortic and mitral valve motion and corresponding fluid dynamics, motivating the use of patient-specific imaging to describe the aortic and mitral valve geometries as well as the motion of the left ventricle (LV). The LV position can be captured at several time points in the cardiac cycle, such that its motion can be prescribed a priori as a Dirichlet boundary condition during a simulation. Valve leaflet motion, however, should be computed from soft-tissue models and incorporated using fully-coupled Fluid Structure Interaction (FSI) algorithms. While FSI simulations have in part or wholly been achieved by multiple groups, to date, no high-throughput models have been developed, which are needed for use in a clinical environment. This project seeks to enable patient-derived moving LV boundary conditions, and has been developed for use with a previously developed immersed boundary, fixed Cartesian grid FSI framework. One challenge in specifying LV motion from medical images stems from the low temporal resolution available. Typical imaging modalities contain only tens of images during the cardiac cycle to describe the change in position of the left ventricle. This temporal resolution is significantly lower than the time resolution needed to capture fluid dynamics of a highly deforming heart valve, and thus an approach to describe intermediate positions of the LV is necessary. Here, we propose a primarily Eulerian means of representing LV displacement. This is a natural extension, since an Eulerian framework is employed in the CFD model to describe the large displacement of the heart valve leaflets. This approach to using Eulerian interface representation is accomplished by applying “morphing” techniques commonly used in the field of computer graphics. For the approach developed in the current work, morphing is adapted to the unique characteristics of a Cartesian grid flow solver which presents challenges of adaptive mesh refinement, narrow band approach, parallel domain decomposition, and the need to supply a local surface velocity to the flow solver that describes both normal and tangential motion. This is accomplished by first generating a skeleton from the Eulerian interface representation, and deforming the skeleton between image frames to determine bulk displacement. After supplying bulk displacement, local displacement is determined using the Eulerian fields. The skeletons are also utilized to automate the simulation setup to track the locations upstream and downstream where the system inflow/outflow boundary conditions are to be applied, which in the current approach, are not limited to Cartesian domain boundaries.