From ff3747fbd3c0d30e6fddc75119dc10b116ec96a4 Mon Sep 17 00:00:00 2001 From: Michael Kuhn Date: Thu, 10 Oct 2024 16:20:11 -0600 Subject: [PATCH] add regression test with init_wave_profile --- test/CMakeLists.txt | 1 + .../ow_linear_init_waves.inp | 70 +++++++++++++++++++ 2 files changed, 71 insertions(+) create mode 100644 test/test_files/ow_linear_init_waves/ow_linear_init_waves.inp diff --git a/test/CMakeLists.txt b/test/CMakeLists.txt index 64f787e225..cda1c2129c 100644 --- a/test/CMakeLists.txt +++ b/test/CMakeLists.txt @@ -246,6 +246,7 @@ add_test_re(zalesak_disk_scalar_vel) add_test_re(rain_drop) add_test_re(inertial_drop) add_test_re(ow_linear) +add_test_re(ow_linear_init_waves) add_test_re(ow_stokes) add_test_re(scalar_advection_uniform) add_test_re(scalar_advection_refined) diff --git a/test/test_files/ow_linear_init_waves/ow_linear_init_waves.inp b/test/test_files/ow_linear_init_waves/ow_linear_init_waves.inp new file mode 100644 index 0000000000..63e8efb925 --- /dev/null +++ b/test/test_files/ow_linear_init_waves/ow_linear_init_waves.inp @@ -0,0 +1,70 @@ +#¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨# +# SIMULATION STOP # +#.......................................# +time.stop_time = 1 # Max (simulated) time to evolve +time.max_step = 20 # Max number of time steps + +#¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨# +# TIME STEP COMPUTATION # +#.......................................# +time.fixed_dt = 0.01 # Use this constant dt if > 0 +time.cfl = 0.95 # CFL factor +time.use_force_cfl= false +#¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨# +# INPUT AND OUTPUT # +#.......................................# +time.plot_interval = 10 # Steps between plot files +time.checkpoint_interval = -1 # Steps between checkpoint files + +#¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨# +# PHYSICS # +#.......................................# +transport.model = TwoPhaseTransport +transport.viscosity_fluid1=0.0 +transport.viscosity_fluid2=0.0 +transport.turbulent_prandtl = 0.3333 +turbulence.model = Laminar + +incflo.physics = MultiPhase OceanWaves +OceanWaves.label = Wave1 +OceanWaves.Wave1.type = LinearWaves +OceanWaves.Wave1.initialize_wave_field = true +OceanWaves.Wave1.wave_height=0.03 +OceanWaves.Wave1.wave_length=2.0 +OceanWaves.Wave1.water_depth=1 +OceanWaves.Wave1.relax_zone_gen_length=1.5 +OceanWaves.Wave1.numerical_beach_length=1.5 +OceanWaves.Wave1.numerical_beach_length_factor=2.0 +MultiPhase.density_fluid1=1000. +MultiPhase.density_fluid2=1. +ICNS.source_terms = GravityForcing +ICNS.use_perturb_pressure = true +ICNS.reconstruct_true_pressure = true +#¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨# +# ADAPTIVE MESH REFINEMENT # +#.......................................# +amr.n_cell = 192 16 32 # Grid cells at coarsest AMRlevel +amr.max_level = 1 +#¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨# +# GEOMETRY # +#.......................................# +geometry.prob_lo = 0.0 -0.3 -1 # Lo corner coordinates +geometry.prob_hi = 7.2 0.3 0.2 # Hi corner coordinates +geometry.is_periodic = 0 1 0 # Periodicity x y z (0/1) + +xlo.type = "wave_generation" +xhi.type = "slip_wall" + +zlo.type = "slip_wall" +zhi.type = "slip_wall" + +tagging.labels = refine0 +tagging.refine0.type = GeometryRefinement +tagging.refine0.shapes = b0 +tagging.refine0.level = 0 + +tagging.refine0.b0.type = box +tagging.refine0.b0.origin = 0.0 -0.3 -0.1 +tagging.refine0.b0.xaxis = 7.2 0.0 0.0 +tagging.refine0.b0.yaxis = 0.0 0.6 0.0 +tagging.refine0.b0.zaxis = 0.0 0.0 0.2 \ No newline at end of file