Sphere Drop + VOF (#530)

Co-authored-by: Wyatt James Horne <[email protected]>

reg_tests/CTestList.cmake

@@ -138,6 +138,7 @@ add_test_r(variableDensNonUniform 2)
 add_test_r(femPassiveScalar 4)
 add_test_r(femFluidsVortex 4)
 add_test_r(re1kChannel 4)
+add_test_r(vofSphereDrop 4)
 add_test_r_rst(kEpsPipe 8)
 if(ENABLE_TIOGA)
   add_test_r(oversetSphereTIOGA 8)

reg_tests/test_files/vofSphereDrop/vofSphereDrop.i

@@ -0,0 +1,184 @@
+Simulation:
+  name: NaluSim
+
+linear_solvers:
+
+  - name: solve_scalar
+    type: tpetra
+    method: gmres
+    preconditioner: sgs 
+    tolerance: 1e-5
+    max_iterations: 50
+    kspace: 50
+    output_level: 0
+
+  - name: solve_cont
+    type: tpetra
+    method: gmres 
+    preconditioner: muelu
+    tolerance: 1e-5
+    max_iterations: 150
+    kspace: 75
+    output_level: 0
+    muelu_xml_file_name: ../../xml/milestone.xml
+
+realms:
+
+  - name: realm_1
+    mesh: ../../mesh/sphere_drop.g
+    use_edges: no       
+    automatic_decomposition_type: rcb
+
+    time_step_control:
+     target_courant: 1.0
+     time_step_change_factor: 2.5
+   
+    equation_systems:
+      name: theEqSys
+      max_iterations: 2 
+   
+      solver_system_specification:
+        velocity: solve_scalar
+        volume_of_fluid: solve_scalar
+        pressure: solve_cont
+
+      systems:
+        - LowMachEOM:
+            name: myLowMach
+            max_iterations: 1 
+            convergence_tolerance: 1e-5
+        - VolumeOfFluid:
+            name: myV
+            max_iterations: 1
+            convergence_tolerance: 1.e-2
+
+    initial_conditions:
+      - constant: ic_1
+        target_name: [ovr_mesh, bg_mesh]
+        value:
+          pressure: 0
+          velocity: [0.0,0.0]
+
+      - user_function: ic_2
+        target_name: [ovr_mesh, bg_mesh]
+        user_function_name:
+          volume_of_fluid: RayleighTaylor
+        user_function_parameters:
+          volume_of_fluid: [0.001, 1.0, 2.0, 0.0001]
+
+    material_properties:
+      target_name: [ovr_mesh, bg_mesh]
+      specifications:
+        - name: density
+          type: volume_of_fluid
+          phase_one: 1.0
+          phase_two: 1.5
+
+        - name: viscosity
+          type: volume_of_fluid
+          phase_one: 0.001
+          phase_two: 0.001
+
+    boundary_conditions:
+
+    - inflow_boundary_condition: top_bot
+      target_name: surface_1
+      inflow_user_data:
+        velocity: [0,0,0]
+
+    - open_boundary_condition: top_bot2
+      target_name: surface_3
+      open_user_data:
+        velocity: [0,0,0]
+        pressure: 0
+
+    - wall_boundary_condition: sides1
+      target_name: surface_4
+      wall_user_data:
+        velocity: [0,0,0]
+
+    - wall_boundary_condition: sides2
+      target_name: surface_2
+      wall_user_data:
+        velocity: [0,0,0]
+
+    - overset_boundary_condition: sphere_edge
+      overset_user_data:
+        percent_overlap: 20.0
+        percent_overlap_inner: 30.0
+        background_block: bg_mesh
+        overset_block: ovr_mesh
+        overset_surface: surface_5
+        background_cut_block: block_10
+        background_cut_surface: surface_20
+
+    - wall_boundary_condition: bc_sphere
+      target_name: surface_6
+      wall_user_data:
+        user_function_name:
+         velocity: mesh_motion
+        user_function_string_parameters:
+         velocity: [mmSphere_ss5]
+
+    solution_options:
+      name: myOptions
+
+
+      mesh_motion:
+
+        - name: mmBackground
+          target_name: [bg_mesh]
+          omega: 0.0
+
+        - name: mmSphere_ss5
+          target_name: [ovr_mesh]
+          include_six_dof: yes
+          body_omega: [0.0,0.0,0.0]
+          body_mass: 0.981747
+          body_velocity: [0,0,0]
+          body_cc_disp: [0,3,0]
+          principal_moments_inertia: [0.0,0.0,1.0]
+          forcing_surface: [surface_6]
+          applied_force: [0.0,-9.81,0.0]
+          compute_centroid: yes
+
+      use_consolidated_solver_algorithm: yes
+      use_consolidated_face_elem_bc_algorithm: yes
+
+      options:
+        - element_source_terms:
+            momentum: [lumped_momentum_time_derivative, advection_diffusion, buoyancy]
+            continuity: [advection]
+            volume_of_fluid: [vof, sucv_nso]
+
+        - user_constants:
+            gravity: [0.0, -10.0]
+            reference_density: 0.0
+
+    output:
+      output_data_base_name: sphere_drop.e
+      output_frequency: 5 
+      output_node_set: no 
+      output_variables:
+       - velocity
+       - pressure
+       - mesh_displacement
+       - fringe_node
+       - dpdx
+       - dudx
+       - intersected_element
+       - mesh_velocity
+       - density
+
+Time_Integrators:
+  - StandardTimeIntegrator:
+      name: ti_1
+      start_time: 0
+      termination_step_count: 20
+      time_step: 0.005
+      time_stepping_type: fixed
+      time_step_count: 0
+      second_order_accuracy: no
+
+      realms:
+        - realm_1

reg_tests/test_files/vofSphereDrop/vofSphereDrop.norm.gold

@@ -0,0 +1,20 @@
+0.02432116329066863 1 0.005
+0.001966663692513982 2 0.01
+0.0003478067418379743 3 0.015
+0.000128807125321143 4 0.02
+0.0001186665578123386 5 0.025
+6.159589527807734e-05 6 0.03
+2.85084248670952e-05 7 0.035
+1.804536444197644e-05 8 0.04
+1.696814237987382e-05 9 0.045
+1.82415952817005e-05 10 0.05
+2.013158697076528e-05 11 0.055
+2.227717639712913e-05 12 0.06
+2.459260313863444e-05 13 0.065
+2.704331805762574e-05 14 0.07
+2.961771509211273e-05 15 0.075
+3.230282985030123e-05 16 0.08
+3.50955110186631e-05 17 0.085
+3.796971046972843e-05 18 0.09
+4.093875506528796e-05 19 0.095
+4.399335630296404e-05 20 0.1