PAVEL: Decorative Patterns with Packed Volumetric Elements

Many real-world hand-crafted objects are decorated with elements that are packed onto the object's surface and deformed to cover it as much as possible. Examples are artisanal ceramics and metal jewelry. Inspired by these objects, we present a method to enrich surfaces with packed volumetric decorations. Our algorithm works by first determining the locations in which to add the decorative elements and then removing the non-physical overlap between them while preserving the decoration volume. For the placement, we support several strategies depending on the desired overall motif. To remove the overlap, we use an approach based on implicit deformable models creating the qualitative effect of plastic warping while avoiding expensive and hard-to-control physical simulations. Our decorative elements can be used to enhance virtual surfaces, as well as 3D-printed pieces, by assembling the decorations onto real-surfaces to obtain tangible reproductions.Comment: 11 page

A Moving Boundary Flux Stabilization Method for Cartesian Cut-Cell Grids using Directional Operator Splitting

An explicit moving boundary method for the numerical solution of time-dependent hyperbolic conservation laws on grids produced by the intersection of complex geometries with a regular Cartesian grid is presented. As it employs directional operator splitting, implementation of the scheme is rather straightforward. Extending the method for static walls from Klein et al., Phil. Trans. Roy. Soc., A367, no. 1907, 4559-4575 (2009), the scheme calculates fluxes needed for a conservative update of the near-wall cut-cells as linear combinations of standard fluxes from a one-dimensional extended stencil. Here the standard fluxes are those obtained without regard to the small sub-cell problem, and the linear combination weights involve detailed information regarding the cut-cell geometry. This linear combination of standard fluxes stabilizes the updates such that the time-step yielding marginal stability for arbitrarily small cut-cells is of the same order as that for regular cells. Moreover, it renders the approach compatible with a wide range of existing numerical flux-approximation methods. The scheme is extended here to time dependent rigid boundaries by reformulating the linear combination weights of the stabilizing flux stencil to account for the time dependence of cut-cell volume and interface area fractions. The two-dimensional tests discussed include advection in a channel oriented at an oblique angle to the Cartesian computational mesh, cylinders with circular and triangular cross-section passing through a stationary shock wave, a piston moving through an open-ended shock tube, and the flow around an oscillating NACA 0012 aerofoil profile.Comment: 30 pages, 27 figures, 3 table

A stable FSI algorithm for light rigid bodies in compressible flow

In this article we describe a stable partitioned algorithm that overcomes the added mass instability arising in fluid-structure interactions of light rigid bodies and inviscid compressible flow. The new algorithm is stable even for bodies with zero mass and zero moments of inertia. The approach is based on a local characteristic projection of the force on the rigid body and is a natural extension of the recently developed algorithm for coupling compressible flow and deformable bodies. Normal mode analysis is used to prove the stability of the approximation for a one-dimensional model problem and numerical computations confirm these results. In multiple space dimensions the approach naturally reveals the form of the added mass tensors in the equations governing the motion of the rigid body. These tensors, which depend on certain surface integrals of the fluid impedance, couple the translational and angular velocities of the body. Numerical results in two space dimensions, based on the use of moving overlapping grids and adaptive mesh refinement, demonstrate the behavior and efficacy of the new scheme. These results include the simulation of the difficult problem of a shock impacting an ellipse of zero mass.Comment: 32 pages, 20 figure

A multi-physics methodology for the simulation of reactive flow and elastoplastic structural response

We propose a numerical methodology for the numerical simulation of distinct, interacting physical processes described by a combination of compressible, inert and reactive forms of the Euler equations, multiphase equations and elastoplastic equations. These systems of equations are usually solved by coupling finite element and CFD models. Here we solve them simultaneously, by recasting all the equations in the same, hyperbolic form and solving them on the same grid with the same finite-volume numerical schemes. The proposed compressible, multiphase, hydrodynamic formulation can employ a hierarchy of five reactive and non-reactive flow models, which allows simple to more involved applications to be directly described by the appropriate selection. The communication between the hydrodynamic and elastoplastic systems is facilitated by means of mixed-material Riemann solvers at the boundaries of the systems, which represent physical material boundaries. To this end we derive approximate mixed Riemann solvers for each pair of the above models based on characteristic equations. The components for reactive flow and elastoplastic solid modelling are validated separately before presenting validation for the full, coupled systems. Multi-dimensional use cases demonstrate the suitability of the reactive flow-solid interaction methodology in the context of impact-driven initiation of reactive flow and structural response due to violent reaction in automotive (e.g. car crash) or defence (e.g. explosive reactive armour) applications. Several types of explosives (C4, Detasheet, nitromethane, gaseous fuel) in gaseous, liquid and solid state are considered.This work was supported by Jaguar Land Rover and the UK-EPSRC grant EP/K014188/1 as part of the jointly funded Programme for Simulation Innovation

A Multi-physics Methodology for Four States of Matter

Abstract: We propose a numerical methodology for the simultaneous numerical simulation of four states of matter: gas, liquid, elastoplastic solids, and plasma. The distinct, interacting physical processes are described by a combination of compressible, inert, and reactive forms of the Euler equations, multi-phase equations, elastoplastic equations, and resistive MHD equations. Combinations of systems of equations are usually solved by coupling finite element for solid modelling and CFD models for fluid modelling or including material effects through boundary conditions rather than full material discretisation. Our simultaneous solution methodology lies on the recasting of all the equations in the same, hyperbolic form allowing their solution on the same grid with the same finite volume numerical schemes. We use a combination of sharp- and diffuse-interface methods to track or capture material interfaces, depending on the application. The communication between the distinct systems of equations (i.e., materials separated by sharp interfaces) is facilitated by means of mixed-material Riemann solvers at the boundaries of the systems, which represent physical material boundaries. To this end, we derive approximate mixed-material Riemann solvers for each pair of the above models based on characteristic equations. To demonstrate the applicability of the new methodology, we consider a case study, where we investigate the possibility of ignition of a combustible gas that lies over a liquid in a metal container that is struck by a plasma arc akin to a lightning strike. We study the effect of the metal container material and its conductivity on the ignition of the combustible gas, as well as the effects of an additional dielectric coating, the sensitivity of the gas, and differences between scenarios with sealed and pre-damaged metal surfaces

A Multi-physics Methodology for Four States of Matter

Abstract: We propose a numerical methodology for the simultaneous numerical simulation of four states of matter: gas, liquid, elastoplastic solids, and plasma. The distinct, interacting physical processes are described by a combination of compressible, inert, and reactive forms of the Euler equations, multi-phase equations, elastoplastic equations, and resistive MHD equations. Combinations of systems of equations are usually solved by coupling finite element for solid modelling and CFD models for fluid modelling or including material effects through boundary conditions rather than full material discretisation. Our simultaneous solution methodology lies on the recasting of all the equations in the same, hyperbolic form allowing their solution on the same grid with the same finite volume numerical schemes. We use a combination of sharp- and diffuse-interface methods to track or capture material interfaces, depending on the application. The communication between the distinct systems of equations (i.e., materials separated by sharp interfaces) is facilitated by means of mixed-material Riemann solvers at the boundaries of the systems, which represent physical material boundaries. To this end, we derive approximate mixed-material Riemann solvers for each pair of the above models based on characteristic equations. To demonstrate the applicability of the new methodology, we consider a case study, where we investigate the possibility of ignition of a combustible gas that lies over a liquid in a metal container that is struck by a plasma arc akin to a lightning strike. We study the effect of the metal container material and its conductivity on the ignition of the combustible gas, as well as the effects of an additional dielectric coating, the sensitivity of the gas, and differences between scenarios with sealed and pre-damaged metal surfaces