Add 2D 2nd-order jacobian test

test/tests/kc/2nd-order-2d-jacobian-tester.i

@@ -0,0 +1,386 @@
+[GlobalParams]
+  gravity = '0 0 0'
+  pspg = true
+  supg = true
+  laplace = true
+  integrate_p_by_parts = true
+  convective_term = true
+  transient_term = true
+  order = SECOND
+[]
+
+[Mesh]
+  type = GeneratedMesh
+  dim = 2
+  xmin = -1
+  xmax = 1
+  ymin = -1
+  ymax = 0
+  nx = 3
+  ny = 3
+  displacements = 'disp_x disp_y'
+  elem_type = QUAD9
+[]
+
+[Problem]
+  kernel_coverage_check = false
+[]
+
+[Variables]
+  [./disp_x]
+  [../]
+  [./disp_y]
+  [../]
+  [./p]
+  [../]
+  [./vel_x]
+  [../]
+  [./vel_y]
+  [../]
+  [./T]
+  [../]
+[]
+
+[Kernels]
+  [./disp_x]
+    type = Diffusion
+    variable = disp_x
+  [../]
+  [./disp_y]
+    type = Diffusion
+    variable = disp_y
+  [../]
+[]
+
+[ADKernels]
+  [./mesh_x]
+    type = INSConvectedMesh
+    variable = vel_x
+    disp_x = disp_x
+    disp_y = disp_y
+    temperature = T
+    use_displaced_mesh = true
+  [../]
+  [./mesh_y]
+    type = INSConvectedMesh
+    variable = vel_y
+    disp_x = disp_x
+    disp_y = disp_y
+    temperature = T
+    use_displaced_mesh = true
+  [../]
+  [./mesh_T]
+    type = INSTemperatureConvectedMesh
+    variable = T
+    disp_x = disp_x
+    disp_y = disp_y
+    use_displaced_mesh = true
+  [../]
+
+# mass
+  [./mass]
+    type = INSADMass
+    variable = p
+    u = vel_x
+    v = vel_y
+    p = p
+    use_displaced_mesh = true
+    temperature = T
+  [../]
+
+  # x-momentum, time
+  [./x_momentum_time]
+    type = INSADMomentumTimeDerivative
+    variable = vel_x
+    use_displaced_mesh = true
+    temperature = T
+  [../]
+
+  # x-momentum, space
+  [./x_momentum_space]
+    type = INSADMomentumBase
+    variable = vel_x
+    u = vel_x
+    v = vel_y
+    p = p
+    temperature = T
+    component = 0
+    use_displaced_mesh = true
+  [../]
+
+  # y-momentum, time
+  [./y_momentum_time]
+    type = INSADMomentumTimeDerivative
+    variable = vel_y
+    use_displaced_mesh = true
+    temperature = T
+  [../]
+
+  # y-momentum, space
+  [./y_momentum_space]
+    type = INSADMomentumBase
+    variable = vel_y
+    u = vel_x
+    v = vel_y
+    p = p
+    temperature = T
+    component = 1
+    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
+   p = p
+   use_displaced_mesh = true
+ [../]
+[]
+
+[BCs]
+  [./x_no_disp]
+    type = DirichletBC
+    variable = disp_x
+    boundary = 'bottom'
+    value = 0
+  [../]
+  [./y_no_disp]
+    type = DirichletBC
+    variable = disp_y
+    boundary = 'bottom'
+    value = 0
+  [../]
+
+  [./x_no_slip]
+    type = DirichletBC
+    variable = vel_x
+    boundary = 'left right bottom'
+    value = 0.0
+  [../]
+
+  [./y_no_slip]
+    type = DirichletBC
+    variable = vel_y
+    boundary = 'left right bottom'
+    value = 0.0
+  [../]
+
+  [./T_cold]
+    type = DirichletBC
+    variable = T
+    boundary = 'bottom'
+    value = 1
+  [../]
+[]
+
+[ADBCs]
+  [./radiation_flux]
+    type = RadiationEnergyFluxBC
+    variable = T
+    boundary = 'top'
+    ff_temp = 1
+    use_displaced_mesh = true
+  [../]
+  [./weld_flux]
+    type = GaussianWeldEnergyFluxBC
+    variable = T
+    boundary = 'top'
+    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
+    variable = vel_x
+    boundary = 'top'
+    component = 0
+    use_displaced_mesh = true
+    temperature = T
+  [../]
+
+  [./vapor_recoil_y]
+    type = VaporRecoilPressureMomentumFluxBC
+    variable = vel_y
+    boundary = 'top'
+    component = 1
+    use_displaced_mesh = true
+    temperature = T
+  [../]
+
+  [./surface_x]
+    type = SurfaceTensionBC
+    variable = vel_x
+    boundary = 'top'
+    component = 0
+    use_displaced_mesh = true
+    temperature = T
+  [../]
+
+  [./surface_y]
+    type = SurfaceTensionBC
+    variable = vel_y
+    boundary = 'top'
+    component = 1
+    use_displaced_mesh = true
+    temperature = T
+  [../]
+
+[./displace_x_top]
+    type = DisplaceBoundaryBC
+    boundary = 'top'
+    variable = 'disp_x'
+    velocity = 'vel_x'
+  [../]
+  [./displace_y_top]
+    type = DisplaceBoundaryBC
+    boundary = 'top'
+    variable = 'disp_y'
+    velocity = 'vel_y'
+  [../]
+[]
+
+[ADMaterials]
+  [./kc_fits]
+    type = CrazyKCPlantFits
+    temperature = T
+    c_mu1 = 1
+    c_mu2 = 1
+    c_mu3 = 1
+    Tmax = 3
+    Tl = 2
+    T90 = 1
+    beta = 1
+    c_k0 = 1
+    c_k1 = 1
+    c_cp0 = 1
+    c_cp1 = 1
+    c_rho0 = 1
+  [../]
+  [./boundary]
+    type = CrazyKCPlantFitsBoundary
+    use_displaced_mesh = true
+    boundary = 'top'
+    temperature = T
+    c_mu0 = 1
+    ap0 = 1
+    ap1 = 1
+    ap2 = 1
+    ap3 = 1
+    bp0 = 1
+    bp1 = 1
+    bp2 = 1
+    bp3 = 1
+    Tb = 1
+    Tbound1 = 2
+    Tbound2 = 3
+    alpha = -1
+    T0 = 2
+  [../]
+[]
+
+[Materials]
+  [./const]
+    type = GenericConstantMaterial
+    prop_names = 'abs sb_constant'
+    prop_values = '1 1'
+  [../]
+[]
+
+[Preconditioning]
+  [./SMP]
+    type = SMP
+    full = true
+    solve_type = 'NEWTON'
+  [../]
+[]
+
+[Executioner]
+  type = Transient
+  end_time = 10000
+  dt = 1
+  dtmin = 1
+  num_steps = 3
+  petsc_options = '-snes_converged_reason -ksp_converged_reason -options_left -snes_test_jacobian'
+  petsc_options_iname = '-pc_type -pc_factor_shift_type -pc_factor_shift_amount -snes_linesearch_minlambda -pc_factor_mat_solver_type -ksp_gmres_restart -snes_test_err'
+  petsc_options_value = 'lu       NONZERO               1e-15                   1e-3                       superlu_dist               100		 1e-9'
+
+  line_search = 'none'
+  nl_max_its = 10
+  l_max_its = 100
+[]
+
+[Outputs]
+  [./exodus]
+    type = Exodus
+    output_material_properties = true
+    show_material_properties = 'mu rc_pressure surface_term_curvature surface_term_gradient1 surface_term_gradient2'
+  [../]
+  [./dofmap]
+    type = DOFMap
+    execute_on = 'initial'
+  [../]
+  checkpoint = true
+[]
+
+[Debug]
+  show_var_residual_norms = true
+[]
+
+[Postprocessors]
+  [./num_dofs]
+    type = NumDOFs
+    system = 'NL'
+  [../]
+[]
+
+[ICs]
+  [./vel_x]
+    type = RandomIC
+    min = 0.1
+    max = 0.2
+    variable = vel_x
+  [../]
+  [./vel_y]
+    type = RandomIC
+    min = 0.1
+    max = 0.2
+    variable = vel_y
+  [../]
+  [./disp_x]
+    type = RandomIC
+    min = 0.1
+    max = 0.2
+    variable = disp_x
+  [../]
+  [./disp_y]
+    type = RandomIC
+    min = 0.1
+    max = 0.2
+    variable = disp_y
+  [../]
+  [./p]
+    type = RandomIC
+    min = 0.1
+    max = 3.9
+    variable = p
+  [../]
+  [./T]
+    type = RandomIC
+    min = 0.1
+    max = 3.9
+    variable = T
+  [../]
+[]

test/tests/kc/tests

@@ -14,4 +14,11 @@
     run_sim = True
     min_parallel = 6
   [../]
+  [./2d-jacobian-2nd-order]
+    type = PetscJacobianTester
+    input = 2nd-order-2d-jacobian-tester.i
+    ratio_tol = 1e-7
+    difference_tol = 1e-4
+    run_sim = True
+  [../]
 []