A unified Eulerian framework for multimaterial continuum mechanics

A framework for simulating the interactions between multiple different continua is presented. Each constituent material is governed by the same set of equations, differing only in terms of their equations of state and strain dissipation functions. The interfaces between any combination of fluids, solids, and vacuum are handled by a new Riemann Ghost Fluid Method, which is agnostic to the type of material on either side (depending only on the desired boundary conditions). The Godunov-Peshkov-Romenski (GPR) model is used for modelling the continua (having recently been used to solve a range of problems involving Newtonian and non-Newtonian fluids, and elastic and elastoplastic solids), and this study represents a novel approach for handling multimaterial problems under this model. The resulting framework is simple, yet capable of accurately reproducing a wide range of different physical scenarios. It is demonstrated here to accurately reproduce analytical results for known Riemann problems, and to produce expected results in other cases, including some featuring heat conduction across interfaces, and impact-induced deformation and detonation of combustible materials. The framework thus has the potential to streamline development of simulation software for scenarios involving multiple materials and phases of matter, by reducing the number of different systems of equations that require solvers, and cutting down on the amount of theoretical work required to deal with the interfaces between materials.EPSRC Centre for Doctoral Training in Computational Methods for Materials Science under grant EP/L015552/

Surface-Only Simulation of Fluids

Surface-only simulation methods for fluid dynamics are those that perform computation only on a surface representation, without relying on any volumetric discretization. Such methods have superior asymptotic complexity in time and memory than the traditional volumetric discretization approaches, and thus are more tractable for simulation of complex fluid phenomena. Although for most computer graphics applications and many engineering applications, the interior flow inside the fluid phases is typically not of interest, the vast majority of existing numerical techniques still rely on discretization of the volumetric domain. My research first tackles the mesh-based surface tracking problem in the multimaterial setting, and then proposes surface-only simulation solutions for two scenarios: the soap-films and bubbles, and the general 3D liquids. Throughout these simulation approaches, all computation takes place on the surface, and volumetric discretization is entirely eliminated

Accurate Sharp Interface Scheme for Multimaterials

We present a method to capture the evolution of a contact discontinuity separating two different material. A locally non-conservative scheme allows an accurate and stable simulation while the interface is kept sharp. Numerical illustrations include problems involving fluid and elastic problems

MFC: An open-source high-order multi-component, multi-phase, and multi-scale compressible flow solver

MFC is an open-source tool for solving multi-component, multi-phase, and bubbly compressible flows. It is capable of efficiently solving a wide range of flows, including droplet atomization, shock–bubble interaction, and bubble dynamics. We present the 5- and 6-equation thermodynamically-consistent diffuse-interface models we use to handle such flows, which are coupled to high-order interface-capturing methods, HLL-type Riemann solvers, and TVD time-integration schemes that are capable of simulating unsteady flows with strong shocks. The numerical methods are implemented in a flexible, modular framework that is amenable to future development. The methods we employ are validated via comparisons to experimental results for shock–bubble, shock–droplet, and shock–water-cylinder interaction problems and verified to be free of spurious oscillations for material-interface advection and gas–liquid Riemann problems. For smooth solutions, such as the advection of an isentropic vortex, the methods are verified to be high-order accurate. Illustrative examples involving shock–bubble-vessel-wall and acoustic–bubble-net interactions are used to demonstrate the full capabilities of MFC

LES of an Inclined Jet into a Supersonic Turbulent Crossflow

This short article describes flow parameters, numerical method, and animations of the fluid dynamics video "LES of an Inclined Jet into a Supersonic Turbulent Crossflow" (http://ecommons.library.cornell.edu/bitstream/1813/14073/3/GFM-2009.mpg [high-resolution] and http://ecommons.library.cornell.edu/bitstream/1813/14073/2/GFM-2009-web.m1v [low-resolution] video). We performed large-eddy simulation with the sub-grid scale (LES-SGS) stretched-vortex model of momentum and scalar transport to study the gas-dynamics interactions of a helium inclined round jet into a supersonic ($M=3.6$) turbulent (\Reth$=13\times10^3$) air flow over a flat surface. The video shows the temporal development of Mach-number and magnitude of density-gradient in the mid-span plane, and isosurface of helium mass-fraction and \lam_2 (vortical structures). The identified vortical structures are sheets, tilted tubes, and discontinuous rings. The vortical structures are shown to be well correlated in space and time with helium mass-fraction isosurface ($Y_{\rm He}=0.25$).Comment: 7 pages, 1 figure, 1 table, article describing fluid dynamics video submitted to Gallery of Fluid Motion, APS-DFD 200

Linear stability analysis of capillary instabilities for concentric cylindrical shells

Motivated by complex multi-fluid geometries currently being explored in fibre-device manufacturing, we study capillary instabilities in concentric cylindrical flows of $N$ fluids with arbitrary viscosities, thicknesses, densities, and surface tensions in both the Stokes regime and for the full Navier--Stokes problem. Generalizing previous work by Tomotika (N=2), Stone & Brenner (N=3, equal viscosities) and others, we present a full linear stability analysis of the growth modes and rates, reducing the system to a linear generalized eigenproblem in the Stokes case. Furthermore, we demonstrate by Plateau-style geometrical arguments that only axisymmetric instabilities need be considered. We show that the N=3 case is already sufficient to obtain several interesting phenomena: limiting cases of thin shells or low shell viscosity that reduce to N=2 problems, and a system with competing breakup processes at very different length scales. The latter is demonstrated with full 3-dimensional Stokes-flow simulations. Many $N > 3$ cases remain to be explored, and as a first step we discuss two illustrative $N \to \infty$ cases, an alternating-layer structure and a geometry with a continuously varying viscosity

Doctor of Philosophy

dissertationVisualizing surfaces is a fundamental technique in computer science and is frequently used across a wide range of fields such as computer graphics, biology, engineering, and scientific visualization. In many cases, visualizing an interface between boundaries can provide meaningful analysis or simplification of complex data. Some examples include physical simulation for animation, multimaterial mesh extraction in biophysiology, flow on airfoils in aeronautics, and integral surfaces. However, the quest for high-quality visualization, coupled with increasingly complex data, comes with a high computational cost. Therefore, new techniques are needed to solve surface visualization problems within a reasonable amount of time while also providing sophisticated visuals that are meaningful to scientists and engineers. In this dissertation, novel techniques are presented to facilitate surface visualization. First, a particle system for mesh extraction is parallelized on the graphics processing unit (GPU) with a red-black update scheme to achieve an order of magnitude speed-up over a central processing unit (CPU) implementation. Next, extending the red-black technique to multiple materials showed inefficiencies on the GPU. Therefore, we borrow the underlying data structure from the closest point method, the closest point embedding, and the particle system solver is switched to hierarchical octree-based approach on the GPU. Third, to demonstrate that the closest point embedding is a fast, flexible data structure for surface particles, it is adapted to unsteady surface flow visualization at near-interactive speeds. Finally, the closest point embedding is a three-dimensional dense structure that does not scale well. Therefore, we introduce a closest point sparse octree that allows the closest point embedding to scale to higher resolution. Further, we demonstrate unsteady line integral convolution using the closest point method