New example cases: Lid-driven cavity (compressible), double Mach, and more

[sbryngelson]

Simulate and add several new examples/test cases to the examples/ directory:

I will accept PRs that add these individually (they do not have to be added all at once). Please add small (in terms of MB) visualization images of the results to the example/ directory.

I will also accept additions to this list, particularly problems in 3D or multiphase 2D problems.

Update: More problems!

Can find refs in doi:10.1016/j.cpc.2024.109236

• 1D: Entropy wave advection problem
• 1D: Shu--Osher, repeat some others from Florian preprint. Sod and Lax shock tubes.
• 2D: 2D Riemann problem' sec 2.5.1 of paper above
• 2D: Mixing layer
• 3D: Inviscid and viscous isotropic Taylor--Green vortex
• 3D: HIT: Forced and decaying
• 3D: Shear turbulence

Original:

• Shu-Osher problem (1D):

  • https://arxiv.org/pdf/2110.05461.pdf

• Lax shock tube problem (1D):

  • https://arxiv.org/pdf/2106.01738.pdf (eqn 59)

• Titarev-Toro problem (1D):

  • https://arxiv.org/pdf/2110.05461.pdf

• Riemann test problem (2D):

  • https://epubs.siam.org/doi/10.1137/0914082
  • https://arxiv.org/pdf/2110.05461.pdf

• Rayleigh-Taylor instability (2D):

  • https://arxiv.org/pdf/2110.05461.pdf (example 4.7)
  • https://www.sciencedirect.com/science/article/abs/pii/S0021999101968925?via%3Dihub
  • https://arxiv.org/pdf/2106.01738.pdf (example 4.7)

• Lid-driven compressible cavity (2D), finite Reynolds number:

  • https://doi.org/10.1016/j.cpc.2022.108527
  • https://arxiv.org/pdf/2106.01738.pdf (example 4.12)

• Double-Mach reflection (2D):

  • https://doi.org/10.1016/j.camwa.2022.12.003
  • https://arxiv.org/pdf/2110.05461.pdf (example 4.8)
  • https://arxiv.org/pdf/2106.01738.pdf (example 4.9)

• Kelvin-Helmholtz instability (2D):

  • https://arxiv.org/pdf/2106.01738.pdf (example 4.8)

• Viscous shock tube (2D):

  • https://arxiv.org/pdf/2110.05461.pdf (example 4.9)
  • https://doi.org/10.1016/j.compfluid.2008.06.008
  • https://arxiv.org/pdf/2204.00393.pdf

• Double periodic shear layer (2D):

  • https://arxiv.org/pdf/2206.14891.pdf (sec 4.2)

• 2D single-material triple point problem here

  • arxiv.org/abs/2312.17042

• Shock-cylinder interaction (2D):

  • https://arxiv.org/pdf/2106.01738.pdf (example 4.17)

• Viscous Richtmeyer-Meshkov instability (2D):

  • https://arxiv.org/pdf/2106.01738.pdf (example 4.18)

• Underwater explosion (2D):

  • https://doi.org/10.1016/j.compfluid.2015.07.010 (fig 12 and 15, sec 4.4)

• Compressible decaying isotropic turbulence (2D and 3D!):

  • https://doi.org/10.1016/j.cpc.2022.108527

• Spherical underwater explosion in a rigid cylindrical container (3D)

  • https://doi.org/10.1016/j.compfluid.2015.07.010 (sec 4.5)

• Liquid jet in crossflow (3D)

  • https://doi.org/10.1016/j.jcp.2014.09.005 (fig 36)
  • https://doi.org/10.1016/j.ijmultiphaseflow.2022.104303
  • https://doi.org/10.1016/j.ast.2019.05.017

• Jet-spray [SprayAB] (3D)

  • https://doi.org/10.1016/j.compfluid.2024.106186
  • https://ecn.sandia.gov/diesel-spray-combustion/target-condition/spray-ab/

[wilfonba]

@sbryngelson How do you want to handle the cases that will require hardcoding analytic initial conditions?

[sbryngelson]

@wilfonba for cases with analytic functions, I think we should add some functionality to make this possible without commenting/uncommenting source code. I'm open to suggestions. Some example possibilities are

• Allowing analytic functions in the Python dictionary
• Adding new functions as new patch types
• etc.

[wilfonba]

What do you think of having the user supply an additional input file with fypp functions that can be included in m_patches.f90?

[sbryngelson]

I'm not sure I understand, I would have to see a small example I think.

[wilfonba]

I'll put something together

[henryleberre]

@wilfonba I’d be happy to discuss possible toolchain implementations.