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

Progress in particle-based multiscale and hybrid methods for flow applications

This work focuses on the review of particle-based multiscale and hybrid methods that have surfaced in the field of fluid mechanics over the last 20 years. We consider five established particle methods: molecular dynamics, direct simulation Monte Carlo, lattice Boltzmann method, dissipative particle dynamics and smoothed-particle hydrodynamics. A general description is given on each particle method in conjunction with multiscale and hybrid applications. An analysis on the length scale separation revealed that current multiscale methods only bridge across scales which are of the order of O(102)−O(103) and that further work on complex geometries and parallel code optimisation is needed to increase the separation. Similarities between methods are highlighted and combinations discussed. Advantages, disadvantages and applications of each particle method have been tabulated as a reference

EFFICIENT PARTICLE-BASED VISCOUS FLUID SIMULATION WITH VIDEO-GUIDED REAL-TO-VIRTUAL PARAMETER TRANSFER

Viscous fluids, such as honey and molten chocolate, are common materials frequently seen in our daily life. These viscous fluids exhibit characteristic behaviors. Capturing and understanding such dynamics have been required for various applications. Although recent research made advances in simulating the viscous fluid dynamics, still many challenges are left to be addressed. In this dissertation, I present novel techniques to more efficiently and accurately simulate viscous fluid dynamics and propose a parameter identification framework to facilitate the tedious parameter tuning steps for viscous materials. In fluid simulation, enforcing the incompressibility robustly and efficiently is essential. One known challenge is how to set appropriate boundary conditions for free surfaces and solid boundaries. I propose a new boundary handling approach for an incompressible particle-based solver based on the connectivity analysis for simulation particles. Another challenge is that previously proposed techniques do not scale well. To address this, I propose a new multilevel particle-based solver which constructs the hierarchy of simulation particles. These techniques improve the robustness and efficiency achieving the nearly linear scaling unlike previous approaches. To simulate characteristic behaviors of viscous fluids, such as coiling and buckling phenomena and adhesion to other materials, it is necessary to develop a specialized solver. I propose a stable and efficient particle-based solver for simulating highly viscous fluids by using implicit integration with the full form of viscosity. To simulate more accurate interactions with solid objects, I propose a new two-way fluid-solid coupling method for viscous fluids via the unified minimization. These approaches also improve the robustness and efficiency while generating rotational and sticky behaviors of viscous fluids. One important challenge for the physically-based simulation is that it is not obvious how to choose appropriate material parameters to generate our desirable behaviors of simulated materials. I propose a parameter identification framework that helps to tune material parameters for viscous fluids with example video data captured from real world fluid phenomena. This framework identifies viscosity parameters for the real viscous fluids while estimating the hidden variables for the fluids, and enables the parameter transfer from the real world to virtual environment.Doctor of Philosoph

Multiscale Universal Interface: A Concurrent Framework for Coupling Heterogeneous Solvers

Concurrently coupled numerical simulations using heterogeneous solvers are powerful tools for modeling multiscale phenomena. However, major modifications to existing codes are often required to enable such simulations, posing significant difficulties in practice. In this paper we present a C++ library, i.e. the Multiscale Universal Interface (MUI), which is capable of facilitating the coupling effort for a wide range of multiscale simulations. The library adopts a header-only form with minimal external dependency and hence can be easily dropped into existing codes. A data sampler concept is introduced, combined with a hybrid dynamic/static typing mechanism, to create an easily customizable framework for solver-independent data interpretation. The library integrates MPI MPMD support and an asynchronous communication protocol to handle inter-solver information exchange irrespective of the solvers' own MPI awareness. Template metaprogramming is heavily employed to simultaneously improve runtime performance and code flexibility. We validated the library by solving three different multiscale problems, which also serve to demonstrate the flexibility of the framework in handling heterogeneous models and solvers. In the first example, a Couette flow was simulated using two concurrently coupled Smoothed Particle Hydrodynamics (SPH) simulations of different spatial resolutions. In the second example, we coupled the deterministic SPH method with the stochastic Dissipative Particle Dynamics (DPD) method to study the effect of surface grafting on the hydrodynamics properties on the surface. In the third example, we consider conjugate heat transfer between a solid domain and a fluid domain by coupling the particle-based energy-conserving DPD (eDPD) method with the Finite Element Method (FEM).Comment: The library source code is freely available under the GPLv3 license at http://www.cfm.brown.edu/repo/release/MUI

ACM Transactions on Graphics

This paper presents a liquid simulation technique that enforces the incompressibility condition using a stream function solve instead of a pressure projection. Previous methods have used stream function techniques for the simulation of detailed single-phase flows, but a formulation for liquid simulation has proved elusive in part due to the free surface boundary conditions. In this paper, we introduce a stream function approach to liquid simulations with novel boundary conditions for free surfaces, solid obstacles, and solid-fluid coupling. Although our approach increases the dimension of the linear system necessary to enforce incompressibility, it provides interesting and surprising benefits. First, the resulting flow is guaranteed to be divergence-free regardless of the accuracy of the solve. Second, our free-surface boundary conditions guarantee divergence-free motion even in the un-simulated air phase, which enables two-phase flow simulation by only computing a single phase. We implemented this method using a variant of FLIP simulation which only samples particles within a narrow band of the liquid surface, and we illustrate the effectiveness of our method for detailed two-phase flow simulations with complex boundaries, detailed bubble interactions, and two-way solid-fluid coupling

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