kc_jacobian_tester clean in 1D with vel_x, disp_x, p, T

idaholab · Nov 30, 2018 · d4f7cfe · d4f7cfe
1 parent c031719
commit d4f7cfe
Showing 1 changed file with 97 additions and 98 deletions.
diff --git a/test/tests/kc/kc_jacobian_tester.i b/test/tests/kc/kc_jacobian_tester.i
@@ -17,7 +17,7 @@
   ymax = 1
   zmin = -1
   zmax = 0
-  nx = 1
+  nx = 2
   ny = 1
   nz = 1
   # displacements = 'disp_x disp_y disp_z'
@@ -42,8 +42,8 @@
   # [./vel_z]
   # [../]
 
-  # [./T]
-  # [../]
+  [./T]
+  [../]
 
   # [./disp_y]
   # [../]
@@ -64,13 +64,14 @@
   #   type = Diffusion
   #   variable = disp_z
   # [../]
-  # [./mesh_x]
-  #   type = INSConvectedMesh
-  #   variable = vel_x
-  #   disp_x = disp_x
-  #   disp_y = disp_y
-  #   # disp_z = disp_z
-  # [../]
+  [./mesh_x]
+    type = INSConvectedMesh
+    variable = vel_x
+    disp_x = disp_x
+    # disp_y = disp_y
+    # disp_z = disp_z
+    use_displaced_mesh = true
+  [../]
   # [./mesh_y]
   #   type = INSConvectedMesh
   #   variable = vel_x
@@ -99,12 +100,12 @@
     use_displaced_mesh = true
   [../]
 
-  # # x-momentum, time
-  # [./x_momentum_time]
-  #   type = INSADMomentumTimeDerivative
-  #   variable = vel_x
-  #   use_displaced_mesh = true
-  # [../]
+  # x-momentum, time
+  [./x_momentum_time]
+    type = INSADMomentumTimeDerivative
+    variable = vel_x
+    use_displaced_mesh = true
+  [../]
 
   # x-momentum, space
   [./x_momentum_space]
@@ -156,31 +157,31 @@
   #   use_displaced_mesh = true
   # [../]
 
-#  # temperature
-#  [./temperature_time]
-#    type = INSADTemperatureTimeDerivative
-#    variable = T
-#    use_displaced_mesh = true
-#  [../]
-
-#  [./temperature_space]
-#    type = INSADTemperature
-#    variable = T
-#    u = vel_x
-#    v = vel_y
-#    # w = vel_z
-#    p = p
-#    use_displaced_mesh = true
-#  [../]
+ # temperature
+ [./temperature_time]
+   type = INSADTemperatureTimeDerivative
+   variable = T
+   use_displaced_mesh = true
+ [../]
+
+ [./temperature_space]
+   type = INSADTemperature
+   variable = T
+   u = vel_x
+   # v = vel_y
+   # w = vel_z
+   p = p
+   use_displaced_mesh = true
+ [../]
 []
 
-# [BCs]
-#   [./x_no_disp]
-#     type = DirichletBC
-#     variable = disp_x
-#     boundary = 'bottom right left top back'
-#     value = 0
-#   [../]
+[BCs]
+  [./x_no_disp]
+    type = DirichletBC
+    variable = disp_x
+    boundary = 'right'
+    value = 0
+  [../]
 #   [./y_no_disp]
 #     type = DirichletBC
 #     variable = disp_y
@@ -194,12 +195,12 @@
 #     value = 0
 #   [../]
 
-#   [./x_no_slip]
-#     type = DirichletBC
-#     variable = vel_x
-#     boundary = 'bottom right left top back'
-#     value = 0.0
-#   [../]
+  [./x_no_slip]
+    type = DirichletBC
+    variable = vel_x
+    boundary = 'right'
+    value = 0.0
+  [../]
 
 #   [./y_no_slip]
 #     type = DirichletBC
@@ -215,44 +216,42 @@
 #     value = 0.0
 #   [../]
 
-#   [./T_cold]
-#     type = DirichletBC
-#     variable = T
-#     boundary = 'back'
-#     value = 300
-#   [../]
-
-#   [./weld_flux]
-#     type = GaussianWeldEnergyFluxBC
-#     variable = T
-#     boundary = 'front'
-#     reff = 1
-#     F0 = 1
-#     R = 1
-#     x_beam_coord = 0
-#     y_beam_coord = 0
-#     z_beam_coord = 0
-#     use_displaced_mesh = true
-#   [../]
-# []
+  [./T_cold]
+    type = DirichletBC
+    variable = T
+    boundary = 'right'
+    value = 1
+  [../]
+[]
 
-# [ADBCs]
-#   [./radiation_flux]
-#     type = RadiationEnergyFluxBC
-#     variable = T
-#     boundary = 'front'
-#     ff_temp = 1
-#     use_displaced_mesh = true
-#   [../]
+[ADBCs]
+  [./radiation_flux]
+    type = RadiationEnergyFluxBC
+    variable = T
+    boundary = 'left'
+    ff_temp = 1
+    use_displaced_mesh = true
+  [../]
+  [./weld_flux]
+    type = GaussianWeldEnergyFluxBC
+    variable = T
+    boundary = 'left'
+    reff = 1
+    F0 = 1
+    R = 1
+    x_beam_coord = 0
+    y_beam_coord = 0
+    z_beam_coord = 0
+    use_displaced_mesh = true
+  [../]
 
-#   [./vapor_recoil_x]
-#     type = VaporRecoilPressureMomentumFluxBC
-#     temperature = T
-#     variable = vel_x
-#     boundary = 'front'
-#     component = 0
-#     use_displaced_mesh = true
-#   [../]
+  [./vapor_recoil_x]
+    type = VaporRecoilPressureMomentumFluxBC
+    variable = vel_x
+    boundary = 'left'
+    component = 0
+    use_displaced_mesh = true
+  [../]
 
 #   [./vapor_recoil_y]
 #     type = VaporRecoilPressureMomentumFluxBC
@@ -272,12 +271,12 @@
 #     use_displaced_mesh = true
 #   [../]
 
-#   [./displace_x_top]
-#     type = DisplaceBoundaryBC
-#     boundary = 'front'
-#     variable = 'disp_x'
-#     velocity = 'vel_x'
-#   [../]
+  [./displace_x_top]
+    type = DisplaceBoundaryBC
+    boundary = 'left'
+    variable = 'disp_x'
+    velocity = 'vel_x'
+  [../]
 #   [./displace_y_top]
 #     type = DisplaceBoundaryBC
 #     boundary = 'front'
@@ -290,12 +289,12 @@
 #     variable = 'disp_z'
 #     velocity = 'vel_z'
 #   [../]
-# []
+[]
 
 [ADMaterials]
   [./kc_fits]
     type = CrazyKCPlantFits
-    # temperature = T
+    temperature = T
     c_mu1 = 1
     c_mu2 = 1
     c_mu3 = 1
@@ -313,7 +312,7 @@
     type = CrazyKCPlantFitsBoundary
     # boundary = 'front'
     boundary = 'left'
-    # temperature = T
+    temperature = T
     c_mu0 = 1
     ap0 = 1
     ap1 = 1
@@ -350,13 +349,13 @@
   end_time = 10000
   dt = 1
   dtmin = 1
-  num_steps = 1
+  num_steps = 2
   petsc_options = '-snes_converged_reason -ksp_converged_reason -options_left -snes_test_jacobian -snes_test_jacobian_view'
   petsc_options_iname = '-pc_type -pc_factor_shift_type -pc_factor_shift_amount -snes_linesearch_minlambda -pc_factor_mat_solver_type -ksp_gmres_restart'
   petsc_options_value = 'lu       NONZERO               1e-15                   1e-3                       superlu_dist               100'
 
   line_search = 'none'
-  nl_max_its = 1
+  nl_max_its = 10
   l_max_its = 100
 []
 
@@ -428,10 +427,10 @@
     max = 3.9
     variable = p
   [../]
-  # [./T]
-  #   type = RandomIC
-  #   min = 0.1
-  #   max = 3.9
-  #   variable = T
-  # [../]
+  [./T]
+    type = RandomIC
+    min = 0.1
+    max = 3.9
+    variable = T
+  [../]
 []