Simulating water-entry/exit problems using Eulerian-Lagrangian and fully-Eulerian fictitious domain methods within the open-source IBAMR library

In this paper we employ two implementations of the fictitious domain (FD) method to simulate water-entry and water-exit problems and demonstrate their ability to simulate practical marine engineering problems. In FD methods, the fluid momentum equation is extended within the solid domain using an additional body force that constrains the structure velocity to be that of a rigid body. Using this formulation, a single set of equations is solved over the entire computational domain. The constraint force is calculated in two distinct ways: one using an Eulerian-Lagrangian framework of the immersed boundary (IB) method and another using a fully-Eulerian approach of the Brinkman penalization (BP) method. Both FSI strategies use the same multiphase flow algorithm that solves the discrete incompressible Navier-Stokes system in conservative form. A consistent transport scheme is employed to advect mass and momentum in the domain, which ensures numerical stability of high density ratio multiphase flows involved in practical marine engineering applications. Example cases of a free falling wedge (straight and inclined) and cylinder are simulated, and the numerical results are compared against benchmark cases in literature.Comment: The current paper builds on arXiv:1901.07892 and re-explains some parts of it for the reader's convenienc

Simulations of propelling and energy harvesting articulated bodies via vortex particle-mesh methods

The emergence and understanding of new design paradigms that exploit flow induced mechanical instabilities for propulsion or energy harvesting demands robust and accurate flow structure interaction numerical models. In this context, we develop a novel two dimensional algorithm that combines a Vortex Particle-Mesh (VPM) method and a Multi-Body System (MBS) solver for the simulation of passive and actuated structures in fluids. The hydrodynamic forces and torques are recovered through an innovative approach which crucially complements and extends the projection and penalization approach of Coquerelle et al. and Gazzola et al. The resulting method avoids time consuming computation of the stresses at the wall to recover the force distribution on the surface of complex deforming shapes. This feature distinguishes the proposed approach from other VPM formulations. The methodology was verified against a number of benchmark results ranging from the sedimentation of a 2D cylinder to a passive three segmented structure in the wake of a cylinder. We then showcase the capabilities of this method through the study of an energy harvesting structure where the stocking process is modeled by the use of damping elements

A parallel interaction potential approach coupled with the immersed boundary method for fully resolved simulations of deformable interfaces and membranes

In this paper we show and discuss the use of a versatile interaction potential approach coupled with an immersed boundary method to simulate a variety of flows involving deformable bodies. In particular, we focus on two kinds of problems, namely (i) deformation of liquid-liquid interfaces and (ii) flow in the left ventricle of the heart with either a mechanical or a natural valve. Both examples have in common the two-way interaction of the flow with a deformable interface or a membrane. The interaction potential approach (de Tullio & Pascazio, Jou. Comp. Phys., 2016; Tanaka, Wada and Nakamura, Computational Biomechanics, 2016) with minor modifications can be used to capture the deformation dynamics in both classes of problems. We show that the approach can be used to replicate the deformation dynamics of liquid-liquid interfaces through the use of ad-hoc elastic constants. The results from our simulations agree very well with previous studies on the deformation of drops in standard flow configurations such as deforming drop in a shear flow or a cross flow. We show that the same potential approach can also be used to study the flow in the left ventricle of the heart. The flow imposed into the ventricle interacts dynamically with the mitral valve (mechanical or natural) and the ventricle which are simulated using the same model. Results from these simulations are compared with ad- hoc in-house experimental measurements. Finally, a parallelisation scheme is presented, as parallelisation is unavoidable when studying large scale problems involving several thousands of simultaneously deforming bodies on hundreds of distributed memory computing processors

Index to NASA Tech Briefs, 1975

This index contains abstracts and four indexes--subject, personal author, originating Center, and Tech Brief number--for 1975 Tech Briefs

Generalized thick strip modelling for vortex-induced vibration of long flexible cylinders

We propose a generalized strip modelling method that is computationally efficient for the VIV prediction of long flexible cylinders in three-dimensional incompressible flow. In order to overcome the shortcomings of conventional strip-theory-based 2D models, the fluid domain is divided into “thick” strips, which are sufficiently thick to locally resolve the small scale turbulence effects and three dimensionality of the flow around the cylinder. An attractive feature of the model is that we independently construct a three-dimensional scale resolving model for individual strips, which have local spanwise scale along the cylinder's axial direction and are only coupled through the structural model of the cylinder. Therefore, this approach is able to cover the full spectrum for fully resolved 3D modelling to 2D strip theory. The connection between these strips is achieved through the calculation of a tensioned beam equation, which is used to represent the dynamics of the flexible body. In the limit, however, a single “thick” strip would fill the full 3D domain. A parallel Fourier spectral/hp element method is employed to solve the 3D flow dynamics in the strip-domain, and then the VIV response prediction is achieved through the strip-structure interactions. Numerical tests on both laminar and turbulent flows as well as the comparison against the fully resolved DNS are presented to demonstrate the applicability of this approach

Real-time lattice boltzmann shallow waters method for breaking wave simulations

We present a new approach for the simulation of surfacebased fluids based in a hybrid formulation of Lattice Boltzmann Method for Shallow Waters and particle systems. The modified LBM can handle arbitrary underlying terrain conditions and arbitrary fluid depth. It also introduces a novel method for tracking dry-wet regions and moving boundaries. Dynamic rigid bodies are also included in our simulations using a two-way coupling. Certain features of the simulation that the LBM can not handle because of its heightfield nature, as breaking waves, are detected and automatically turned into splash particles. Here we use a ballistic particle system, but our hybrid method can handle more complex systems as SPH. Both the LBM and particle systems are implemented in CUDA, although dynamic rigid bodies are simulated in CPU. We show the effectiveness of our method with various examples which achieve real-time on consumer-level hardware.Peer ReviewedPostprint (author's final draft