Shock Droplet (2D)

Reference: Panchal et. al., A Seven-Equation Diffused Interface Method for Resolved Multiphase Flows, JCP, 475 (2023)

Initial Condition

Result

2D Riemann Test (2D)

Reference: Chamarthi, A., & Hoffmann, N., & Nishikawa, H., & Frankel S. (2023). Implicit gradients based conservative numerical scheme for compressible flows. arXiv:2110.05461

Density Initial Condition

Density

Density Final Condition

Density Norms

3D Weak Scaling

The 3D_weak_scaling case depends on two parameters:

The number of MPI ranks (procs): As procs increases, the problem size per rank remains constant. procs is determined using information provided to the case file by mfc.sh run.
GPU memory usage per rank (gbpp): As gbpp increases, the problem size per rank increases and the number of timesteps decreases so that wall times consistent. gbpp is a user-defined optional argument to the case.py file. It can be specified right after the case filepath when invoking mfc.sh run.

Weak scaling benchmarks can be produced by keeping gbpp constant and varying procs.

For example, to run a weak scaling test that uses ~4GB of GPU memory per rank on 8 2-rank nodes with case optimization, one could:

./mfc.sh run examples/3D_weak_scaling/case.py 4 -t pre_process simulation          \
             -e batch -p mypartition -N 8 -n 2 -w "01:00:00" -# "MFC Weak Scaling" \
             --case-optimization -j 32

Shu-Osher problem (1D)

Reference: C. W. Shu, S. Osher, Efficient implementation of essentially non-oscillatory shock-capturing schemes, Journal of Computational Physics 77 (2) (1988) 439–471. doi:10.1016/0021-9991(88)90177-5.

Initial Condition

Result

Titarev-Toro problem (1D)

Reference: V. A. Titarev, E. F. Toro, Finite-volume WENO schemes for three-dimensional conservation laws, Journal of Computational Physics 201 (1) (2004) 238–260.

Initial Condition

Result

2D Hardcodied IC Example

Initial Condition

Result

Lax shock tube problem (1D)

Reference: P. D. Lax, Weak solutions of nonlinear hyperbolic equations and their numerical computation, Communications on pure and applied mathematics 7 (1) (1954) 159–193.

Initial Condition

Result

Isentropic vortex problem (2D)

Reference: Coralic, V., & Colonius, T. (2014). Finite-volume Weno scheme for viscous compressible multicomponent flows. Journal of Computational Physics, 274, 95–121. https://doi.org/10.1016/j.jcp.2014.06.003

Density

Density Norms

Lid-Driven Cavity Problem (2D)

Reference: Bezgin, D. A., & Buhendwa A. B., & Adams N. A. (2022). JAX-FLUIDS: A fully-differentiable high-order computational fluid dynamics solver for compressible two-phase flows. arXiv:2203.13760

Reference: Ghia, U., & Ghia, K. N., & Shin, C. T. (1982). High-re solutions for incompressible flow using the Navier-Stokes equations and a multigrid method. Journal of Computational Physics, 48, 387-411

Video: https://youtube.com/shorts/JEP28scZrBM?feature=share

Final Condition

Centerline Velocities