added reg. test for the transition model

reg_tests/CTestList.cmake

@@ -215,6 +215,7 @@ if(NOT ENABLE_CUDA AND NOT ENABLE_ROCM)
     add_test_r(ablNeutralNGPHypreSegregated 2)
     add_test_r(airfoilRANSEdgeNGPHypre.rst 2)
     add_test_r(fsiTurbineSurrogate 4)
+    add_test_r(airfoilSST_Gamma_Trans 4)
 
     if (ENABLE_TRILINOS_SOLVERS)
       add_test_r(aslNeutralEdgeSST 4)

reg_tests/test_files/airfoilSST_Gamma_Trans/airfoilSST_Gamma_Trans.norm.gold

@@ -0,0 +1,5 @@
+9886408.858641412 1 0.0146628
+3938.732212071296 2 0.0293255
+43.78119071147142 3 0.0439883
+26.35564999344441 4 0.058651
+18.24183547831864 5 0.0733138

reg_tests/test_files/airfoilSST_Gamma_Trans/airfoilSST_Gamma_Trans.yaml

@@ -0,0 +1,316 @@
+
+Simulations:
+- name: sim1
+  time_integrator: ti_1
+  optimizer: opt1
+
+# Hypre memory and other configurations
+hypre_config:
+    memory_pool_mbs: 3000
+    use_vendor_sgemm: no
+
+linear_solvers:
+  # solver for the momentum
+  - name: solve_mom
+    type: hypre
+    method: hypre_gmres
+    preconditioner: boomerAMG
+    tolerance: 1e-5
+    max_iterations: 100
+    kspace: 20
+    segregated_solver: yes
+    write_matrix_files: no
+    #write_preassembly_matrix_files: yes
+    recompute_preconditioner_frequency: 1
+    reuse_linear_system: yes
+    simple_hypre_matrix_assemble: yes
+    dump_hypre_matrix_stats: no
+
+    bamg_max_levels: 1
+    bamg_relax_type: 12
+    bamg_num_sweeps: 2
+    bamg_cycle_type: 1
+    bamg_relax_order: 0
+
+    bamg_output_level: 0
+    bamg_debug: 0
+    output_level: 0
+
+  # solver for the scalars TKE, SDR
+  - name: solve_scalar
+    type: hypre
+    method: hypre_gmres
+    preconditioner: boomerAMG
+    tolerance: 1e-5
+    max_iterations: 100
+    kspace: 20
+    write_matrix_files: no
+    #write_preassembly_matrix_files: yes
+    recompute_preconditioner_frequency: 1
+    reuse_linear_system: yes
+    simple_hypre_matrix_assemble: yes
+    dump_hypre_matrix_stats: no
+
+    bamg_max_levels: 1
+    bamg_relax_type: 12
+    bamg_num_sweeps: 2
+    bamg_cycle_type: 1
+    bamg_relax_order: 0
+
+    bamg_output_level: 0
+    bamg_debug: 0
+    output_level: 0
+
+  # solver for the pressure
+  - name: solve_elliptic
+    type: hypre
+    method: hypre_gmres
+    preconditioner: boomerAMG
+    tolerance: 1e-5
+    max_iterations: 100
+    kspace: 40
+    write_matrix_files:  no
+    #write_preassembly_matrix_files: yes
+    recompute_preconditioner_frequency: 1
+    reuse_linear_system: yes
+    simple_hypre_matrix_assemble: yes
+    dump_hypre_matrix_stats: no
+
+    #bamg_max_levels: 7
+    bamg_coarsen_type: 8
+    bamg_interp_type: 6
+    bamg_relax_type: 11
+    bamg_cycle_type: 1
+    bamg_relax_order: 0
+    bamg_keep_transpose: 1
+
+    bamg_num_sweeps: 2
+    #bamg_num_down_sweeps: 2
+    #bamg_num_up_sweeps: 2
+    #bamg_num_coarse_sweeps: 1
+
+    bamg_strong_threshold: 0.3243198904981685
+    bamg_trunc_factor: 0.4227143073461007
+    bamg_pmax_elmts: 3
+
+    bamg_agg_num_levels: 1
+    bamg_agg_interp_type: 7
+    bamg_agg_pmax_elmts: 3
+
+    bamg_output_level: 0
+    bamg_debug: 0
+    output_level: 0
+
+  # solver for the wall
+  - name: solve_wall
+    type: hypre
+    method: hypre_gmres
+    preconditioner: boomerAMG
+    tolerance: 1e-5
+    max_iterations: 100
+    kspace: 40
+    write_matrix_files:  no
+    #write_preassembly_matrix_files: yes
+    recompute_preconditioner_frequency: 1
+    reuse_linear_system: yes
+    simple_hypre_matrix_assemble: yes
+    dump_hypre_matrix_stats: no
+
+    bamg_max_levels: 7
+    bamg_coarsen_type: 8
+    bamg_interp_type: 6
+    bamg_relax_type: 11
+    bamg_num_sweeps: 2
+    bamg_cycle_type: 1
+    bamg_relax_order: 0
+    bamg_trunc_factor: 0.1
+    bamg_agg_num_levels: 2
+    bamg_agg_interp_type: 5
+    bamg_agg_pmax_elmts: 3
+    bamg_pmax_elmts: 3
+    bamg_strong_threshold: 0.25
+
+    bamg_output_level: 0
+    bamg_debug: 0
+    output_level: 0
+
+realms:
+
+  - name: realm_1
+    mesh: ../../mesh/du91w2_airfoil.exo
+    automatic_decomposition_type: rcb
+    use_edges: yes
+
+    time_step_control:
+     target_courant: 1000.0
+     time_step_change_factor: 1.05
+
+    equation_systems:
+      name: theEqSys
+      max_iterations: 4
+
+      solver_system_specification:
+        velocity: solve_mom
+        turbulent_ke: solve_scalar
+        specific_dissipation_rate: solve_scalar
+        pressure: solve_elliptic
+        ndtw: solve_wall
+        gamma_transition: solve_scalar
+
+      systems:
+
+        - WallDistance:
+            name: myNDTW
+            max_iterations: 1
+            convergence_tolerance: 1e-8
+
+        - LowMachEOM:
+            name: myLowMach
+            max_iterations: 1
+            convergence_tolerance: 1e-8
+
+        - ShearStressTransport:
+            name: mySST
+            max_iterations: 1
+            convergence_tolerance: 1e-8
+
+    initial_conditions:
+      - constant: ic_1
+        target_name: [Flow-QUAD]
+        value:
+          pressure: 0
+          velocity: [34.1, 0.0]
+          turbulent_ke: 0.00392448375
+          specific_dissipation_rate: 460.34999999999997
+          gamma_transition: 1.0
+
+    material_properties:
+      target_name: [Flow-QUAD]
+      specifications:
+        - name: density
+          type: constant
+          value: 1.225
+        - name: viscosity
+          type: constant
+          value: 1.0443125000000002e-05
+
+    boundary_conditions:
+
+    - wall_boundary_condition: bc_wall
+      target_name: airfoil
+      wall_user_data:
+        velocity: [0,0]
+        use_wall_function: no
+        turbulent_ke: 0.0
+
+
+    - inflow_boundary_condition: bc_inflow
+      target_name: inlet
+      inflow_user_data:
+        velocity: [34.1, 0.0]
+        turbulent_ke: 0.00392448375
+        specific_dissipation_rate: 460.34999999999997
+
+    - open_boundary_condition: bc_open
+      target_name: outlet
+      open_user_data:
+        velocity: [0,0]
+        pressure: 0.0
+        turbulent_ke: 0.00392448375
+        specific_dissipation_rate: 460.34999999999997
+
+    solution_options:
+      name: myOptions
+      turbulence_model: sst
+      transition_model: true
+      projected_timescale_type: momentum_diag_inv #### Use 1/diagA formulation
+
+      options:
+        - hybrid_factor:
+            velocity: 1.0
+            turbulent_ke: 1.0
+            specific_dissipation_rate: 1.0
+            gamma_transition: 1.0
+
+        - upw_factor:
+            velocity: 1.0
+            turbulent_ke: 0.0
+            specific_dissipation_rate: 0.0
+            gamma_transition: 0.0
+
+        - alpha_upw:
+            velocity: 1.0
+            turbulent_ke: 1.0
+            specific_dissipation_rate: 1.0
+            gamma_transition: 1.0
+
+        - noc_correction:
+            pressure: yes
+
+        - limiter:
+            pressure: no
+            velocity: yes
+            turbulent_ke: yes
+            specific_dissipation_rate: yes
+
+        - projected_nodal_gradient:
+            velocity: element
+            pressure: element
+            turbulent_ke: element
+            specific_dissipation_rate: element
+            gamma_transition: element
+
+        - relaxation_factor:
+            velocity: 0.7
+            pressure: 0.3
+            turbulent_ke: 0.7
+            specific_dissipation_rate: 0.7
+            gamma_transition: 0.7
+
+        - turbulence_model_constants:
+            fsti: 0.15
+
+    #Use this to set the mean angle of attack
+    mesh_transformation:
+      - name: interior
+        mesh_parts: [ Flow-QUAD ]
+        motion:
+         - type: rotation
+           angle: 0.0
+           axis: [0.0, 0.0, -1.0]
+           centroid: [0.25, 0.0, 0.0]
+
+
+    post_processing:
+    - type: surface
+      physics: surface_force_and_moment
+      output_file_name: results/forces.dat
+      frequency: 1
+      parameters: [0.25, 0]
+      target_name: airfoil
+
+    output:
+      output_data_base_name: results/du91w2.e
+      output_frequency: 5000
+      output_node_set: no
+      output_variables:
+       - velocity
+       - pressure
+       - turbulent_ke
+       - specific_dissipation_rate
+       - turbulent_viscosity
+       - minimum_distance_to_wall
+       - gamma_transition
+
+Time_Integrators:
+  - StandardTimeIntegrator:
+      name: ti_1
+      start_time: 0
+      time_step: 0.014662756598240468
+      termination_step_count: 5
+      time_stepping_type: fixed
+      time_step_count: 0
+      second_order_accuracy: yes
+
+      realms:
+        - realm_1