Lock-Exchange Implementation

[surbhidhiman]

Motivation
A potential example of the lock-exchange problem involves two fluids of different densities separated by a vertical barrier on a slope, and the system evolves into an underwater density current. My goal is to use this case to deepen my understanding of computational fluid dynamics, boundary conditions, and density-driven mixing. My ideal outcome is to produce a fully documented example script that meets the standards of the existing examples.

Proposed Implementation:

• Set up a RectilinearGrid
  • ImmersedBoundaryGrid to help define the slope in the grid
  • Topology (Bounded, Flat, Bounded)
• Set up a NonhydrostaticModel with a buoyancy tracer
  • Boussinesq approximation
• Set up the initial conditions
  • Two fluids of different densities are separated by a vertical barrier that is removed at t=0 (lock exchange)
• Run simulation

This outline is still a work in progress. I plan to continue refining this implementation and including some interesting new feature, and I will keep this issue updated with my progress. If there’s anything I should consider or any feedback from the community, I’d be grateful to hear it. Thank you!

[simone-silvestri]

We already have a lock-exchange setup in the validation case at https://github.com/CliMA/Oceananigans.jl/blob/main/validation/z_star_coordinate/lock_release.jl

Would be quite simple to spruce it up by adding some comments and a nice visualization, and convert it into an example.

[surbhidhiman]

We already have a lock-exchange setup in the validation case at https://github.com/CliMA/Oceananigans.jl/blob/main/validation/z_star_coordinate/lock_release.jl

Would be quite simple to spruce it up by adding some comments and a nice visualization, and convert it into an example.

@simone-silvestri Thank you for sharing this!