If you would like to contribute to our gallery, please contact us.

Atomization of a Liquid Jet

All-Mach, Compressible Multiphase Flow

Simulation of primary atomization of a liquid (water) jet in a Mach 1.94 supersonic crossflow (air). An all-Mach, compressible multiphase flow solver is developed within the framework of AMReX. The Interface Reconstruction Library (IRL) was used to implement an efficient, discretely conservative, unsplit, geometric volume-of-fluid (VOF) transport scheme. The liquid-gas interface in each mixed cell is represented using piecewise linear interface calculation (PLIC). The movie shows streamwise velocity contours with the 4 level mesh and an isocontour of VOF with the finest level of refinement.

A Robust All-Mach Multiphase Flow Algorithm for High-Fidelity Simulations of Compressible Atomization, M.B. Kuhn and O. Desjardins, ILASS-Americas 30th Annual Conference on Liquid Atomization and Spray Systems, Tempe, AZ, May 2019.

An all-Mach multiphase flow solver using block-structured AMR, M. Natarajan, R. Chiodi, M. Kuhn and O. Desjardins, ILASS-Americas 30th Annual Conference on Liquid Atomization and Spray Systems, Tempe, AZ, May 2019.

Image courtesy M. Natarajan, R. Chiodi, M. B. Kuhn, O. Desjardins, Computational Thermo-Fluids Laboratory, Cornell University.

Microstructure Evolution Using Solid Mechanics

Simulation of microstructure evolution in a polycrystalline solid using Alamo

In this simulation, the multiphase field method is used and the evolution equation is integrated explicitly. Microstructure evolution is driven by boundary curvature (as in high temperature annealing) which causes coarsening. The microstructure is initialized using a Voronoi tesselation with 40 initial grains. The simulation has three levels of mesh refinement, and was run on the Texas Advanced Computing Center Stampede2 computer with 512 MPI processes for 10 hours.

Image courtesy of Brandon Runnels, Multiscale Materials Modeling Group, University of Colorado, Colorado Springs, CO.

Detonation Propagation and Failure

Modeling Compressible Reactive Gas Dynamics

Simulation of detonation propagation and failure by diffraction with HyBurn. The detonation is initiated in a reactant layer bounded by high-temperature products that have a very low acoustic impedance. The simulation used with 6 levels of refinement with a refinement ratio of 256 between the finest and coarsest levels. The detonation is initiated by a series of high-pressure, high-temperature spots. The detonation propagates steadily until it encounters a step change in the height of the reactant layer. The leading shock of the detonation weakens as it diffracts around the step, resulting in a decoupled shock and flame. The movie shows the temperature field and follows the detonation.

Image courtesy Brayden Roque, Hsiao-Chi Li, and Ryan Houim

Three-dimensional Hydrogen Jet

Detailed adaptive simulation of a burning Hydrogen jet with RNS

Three-dimensional premixed Hydrogen/Air flame computed with the RNS code. RNS is a block-structured AMR code that solves the compressible reactive Navier-Stokes equations with detailed models for the chemistry, and is based on high-order numerical methods (AMLSDC and WENO) that achieve fourth-order accuracy in both time and space.

  • A Fourth-Order Adaptive Mesh Refinement Algorithm for the Multicomponent, Reacting Compressible Navier-Stokes Equations , M. Emmett, E. Motheau, W. Zhang, M. Minion and J. B. Bell, submitted for publication, 2018.[arxiv]
Image courtesy Emmanuel Motheau

Yield-stress fluids

Capturing transient behaviour of strain-rate-dependent rheological models

Lid-driven cavity problem for a Papanastasiou-regularised Bingham fluid with Reynold's number 1000, Bingham number 1 and regularisation parameter 0.0025. The heatmap shows the effective viscosity distribution, while the black contour lines illustrate the location of the surface where the stress magnitude equals a characteristic threshold value. In order to accurately resolve this yield surface, we utilise adaptive mesh refinement with stress-triggered cell tagging in three layers.

  • Time-dependent viscoplastic fluid flow simulations in two and three dimensions, K. Sverdrup, N. Nikiforakis and A. Almgren, in preparation, 2018, arXiv:1803.00417
Image courtesy of Knut Sverdrup, Laboratory for Scientific Computing, University of Cambridge.

Dimethyl Ether Jet

Detailed adaptive simulation of a burning Dimethyl Ether jet with RNS

An instantaneous snapshot of the temperature field of a Dimethyl Ether flame computed with the RNS code. RNS is a block-structured AMR code that solves the compressible reactive Navier-Stokes equations with detailed models for the chemistry, and is based on high-order numerical methods (AMLSDC and WENO) that achieve fourth-order accuracy in both time and space.

  • A Fourth-Order Adaptive Mesh Refinement Algorithm for the Multicomponent, Reacting Compressible Navier-Stokes Equations , M. Emmett, E. Motheau, W. Zhang, M. Minion and J. B. Bell, submitted for publication, 2018.[arxiv]
Image courtesy Emmanuel Motheau

Laser Wakefield Acceleration

Modeling the interaction between plasma and electrons

Simulation of laser wakefield acceleration performed with WarpX. The laser pulse propagates from left to right in a uniform plasma. A moving window is used, i.e. the simulation box travels at the speed of light to follow the laser pulse. The central slice of plasma electrons is shown as transparent white dots. A cavity free of plasma electrons forms in the laser wake, where an electron bunch (solid white dots) is accelerated. The colormap shows the longitudinal electric field in the wake. The white box in the center shows the mesh-refined area.

Image courtesy Maxence Thévenet & the WarpX team.

Shock Reflection

Compressible Gas Dynamics with AMR Embedded Boundaries

A compressible gas dynamics shock reflection using an embedded boundary representation of the ramp. The colors represent the density field. There are 3 total levels of refinement. The code for this simulation is available in the AMReX tutorial, amrex/Tutorials/EB/CNS/Exec/ShockRef/.

Black hole Collisions

Code Generation and Simulation of the Spacetime Evolution of Black Hole Mergers

The video to the right shows the spacetime evolution of two equal mass spinning (Kerr) black holes merging into a single black hole with outward propagating gravitational waves. The evolution was simulated using the Z4c formulation of the Einstein equations wtih 6 levels of AMR in AMReX. The complex equations of motion were generated using the STvAR package, designed for converting symbolic/tensorial forms to executable code for AMReX.

The STvAR package was inspired by Professor Zach Etienne and the NRPy+ project at West Virginia University.

Videos and Images courtesy of Adam Peterson and Don Wilcox.

Black Hole Advertising

Simulation of Large Systems of Black Holes with Heavily Designed Initial Conditions

The video to the left shows an amusing application of the STvAR code generators and AMReX to spacetime evolution. Click on the image to see what happens and then read the description below.

The video shows the simulation of a very large system of black holes, starting from heavily designed initial conditions using the STvAR package. The entire video is run in reverse.

Videos and Images courtesy of Adam Peterson and Don Wilcox.