Nonlinear solve did not converge due to DIVERGED_FNORM_NAN iterations 1

[Emmanuel-324]

I get the error below from my simulation

I have read the debugging procedure on the moose website and even created a debug executable yet I am still confused how to approach this problem.

Below is the error:

etting Up
  Initializing
    Finished Initializing Equation Systems                                               [�[33m  1.04 s�[39m] [�[33m  397 MB�[39m]
  Finished Initializing                                                                  [�[33m  1.09 s�[39m] [�[33m  397 MB�[39m]
Using EXPERIMENTAL Stateful Material Property projection with Adaptivity!
Finished Setting Up                                                                      [�[33m  1.42 s�[39m] [�[33m  397 MB�[39m]
Framework Information:
MOOSE Version:           git commit d438806199 on 2025-01-17
LibMesh Version:         
PETSc Version:           3.22.1
SLEPc Version:           3.22.1
Current Time:            Fri Jan 31 16:02:28 2025
Executable Timestamp:    Fri Jan 17 02:55:05 2025

Input File(s):
  /home/emmanuel324/projects/belson324/Crystal_Plasticity/1-31-25/Test1/New_work5.i

Checkpoint:
  Wall Time Interval:      Every 3600 s
  User Checkpoint:         Disabled
  # Checkpoints Kept:      2
  Execute On:              TIMESTEP_END 

Parallelism:
  Num Processors:          1
  Num Threads:             1

Mesh: 
  Parallel Type:           replicated
  Mesh Dimension:          2
  Spatial Dimension:       2
  Nodes:                   56827
  Elems:                   52497
  Num Subdomains:          1

Nonlinear System:
  Num DOFs:                227308
  Num Local DOFs:          227308
  Variables:               { "disp_x" "disp_y" "eta0" "eta1" } 
  Finite Element Types:    "LAGRANGE" 
  Approximation Orders:    "FIRST" 

Auxiliary System:
  Num DOFs:                371809
  Num Local DOFs:          371809
  Variables:               "vonmises" "stress_xx" "temperature" { "eth_xx" "eth_yy" "fth_xx" "fth_yy" 
                             } 
  Finite Element Types:    "MONOMIAL" "MONOMIAL" "LAGRANGE" "MONOMIAL" 
  Approximation Orders:    "CONSTANT" "CONSTANT" "FIRST" "CONSTANT" 

Execution Information:
  Executioner:             Transient
  TimeStepper:             IterationAdaptiveDT
  TimeIntegrator(s):       ImplicitEuler
  Solver Mode:             Preconditioned JFNK
  PETSc Preconditioner:    lu mat_superlu_dist_fact: SamePattern  mat_superlu_dist_replacetinypivot: true 
  MOOSE Preconditioner:    SMP

Currently Executing
  Performing Initial Setup
    Finished Setting Up Initial Conditions                                               [�[33m  0.52 s�[39m] [�[33m  894 MB�[39m]
    Setting Up Materials
      Finished Computing Initial Material Values                                         [�[33m  1.25 s�[39m] [�[33m 1738 MB�[39m]
    Finished Setting Up Materials                                                        [�[33m  1.25 s�[39m] [�[33m 1738 MB�[39m]
    Currently Performing Initial Adaptivity
      Handling Mesh Changes..                                                            [�[33m 20.86 s�[39m] [�[33m 4301 MB�[39m]
    Finished Performing Initial Adaptivity                                               [�[33m 26.32 s�[39m] [�[33m 4351 MB�[39m]
  Finished Performing Initial Setup                                                      [�[33m 32.95 s�[39m] [�[33m 4386 MB�[39m]
  Outputting
    Outputting Step
      Handling Mesh Changes.                                                             [�[33m 14.66 s�[39m] [�[33m 4735 MB�[39m]
    Finished Outputting Step                                                             [�[33m 16.86 s�[39m] [�[33m 4735 MB�[39m]

Time Step 0, time = 0

Postprocessor Values:
+----------------+----------------+----------------+----------------+----------------+----------------+----------------+
| time           | eth_xx         | eth_yy         | fth_xx         | fth_yy         | stress_xx      | temperature    |
+----------------+----------------+----------------+----------------+----------------+----------------+----------------+
|   0.000000e+00 |   0.000000e+00 |   0.000000e+00 |   0.000000e+00 |   0.000000e+00 |   0.000000e+00 |   0.000000e+00 |
+----------------+----------------+----------------+----------------+----------------+----------------+----------------+

  Finished Outputting                                                                    [�[33m 16.87 s�[39m] [�[33m 4735 MB�[39m]

Time Step 1, time = 5e-08, dt = 5e-08

    Computing Residual....                                                               [�[33m 27.78 s�[39m] [�[33m 4735 MB�[39m]
    |residual|_2 of individual variables:
                     disp_x: 213504
                     disp_y: 106835
                     eta0:   0
                     eta1:   0
 0 Nonlinear |R| = �[32m2.387418e+05�[39m
    Computing Jacobian......                                                             [�[33m 35.94 s�[39m] [�[33m 4939 MB�[39m]
      0 Linear |R| = �[32m2.387418e+05�[39m

    Computing Residual.....                                                              [�[33m 32.23 s�[39m] [�[33m 8747 MB�[39m]
      1 Linear |R| = �[32m4.238450e+04�[39m
    Computing Residual.....                                                              [�[33m 32.20 s�[39m] [�[33m 8747 MB�[39m]
      2 Linear |R| = �[32m1.419121e+04�[39m
    Computing Residual.....                                                              [�[33m 31.77 s�[39m] [�[33m 8802 MB�[39m]
      3 Linear |R| = �[32m2.819610e+02�[39m
    Computing Residual.....                                                              [�[33m 30.16 s�[39m] [�[33m 8802 MB�[39m]
      4 Linear |R| = �[32m3.681539e-01�[39m
    Computing Residual.....                                                              [�[33m 30.19 s�[39m] [�[33m 8802 MB�[39m]

ComputeMultipleCrystalPlasticityStress_abs: Constitutive failure
To recover, the solution will fail and then be re-attempted with a reduced time step.

    |residual|_2 of individual variables:
                     disp_x: 213504
                     disp_y: 106835
                     eta0:   0
                     eta1:   0
 1 Nonlinear |R| = �[33m2.387418e+05�[39m
  Nonlinear solve did not converge due to DIVERGED_FNORM_NAN iterations 1
�[31m Solve Did NOT Converge!�[39m
  Finished Solving                                                                       [�[33m249.51 s�[39m] [�[33m 6646 MB�[39m]
Aborting as solve did not converge

Solve failed, cutting timestep.

Time Step 1, time = 2.5e-08, dt = 2.5e-08

    Computing Residual....                                                               [�[33m 26.06 s�[39m] [�[33m 6652 MB�[39m]
    |residual|_2 of individual variables:
                     disp_x: 213504
                     disp_y: 106835
                     eta0:   0
                     eta1:   0
 0 Nonlinear |R| = �[32m2.387418e+05�[39m
    Computing Jacobian......                                                             [�[33m 35.11 s�[39m] [�[33m 6653 MB�[39m]
      0 Linear |R| = �[32m2.387418e+05�[39m

    Computing Residual.....                                                              [�[33m 32.00 s�[39m] [�[33m 8862 MB�[39m]
      1 Linear |R| = �[32m4.238450e+04�[39m
    Computing Residual.....                                                              [�[33m 31.72 s�[39m] [�[33m 8862 MB�[39m]
      2 Linear |R| = �[32m1.419121e+04�[39m
    Computing Residual.....                                                              [�[33m 32.82 s�[39m] [�[33m 8917 MB�[39m]
      3 Linear |R| = �[32m2.819615e+02�[39m
    Computing Residual.....                                                              [�[33m 32.87 s�[39m] [�[33m 8917 MB�[39m]
      4 Linear |R| = �[32m3.679795e-01�[39m
    Computing Residual.....                                                              [�[33m 31.14 s�[39m] [�[33m 8917 MB�[39m]

ComputeMultipleCrystalPlasticityStress_abs: Constitutive failure
To recover, the solution will fail and then be re-attempted with a reduced time step.

    |residual|_2 of individual variables:
                     disp_x: 213504
                     disp_y: 106835
                     eta0:   0
                     eta1:   0
 1 Nonlinear |R| = �[33m2.387418e+05�[39m
  Nonlinear solve did not converge due to DIVERGED_FNORM_NAN iterations 1
�[31m Solve Did NOT Converge!�[39m
  Finished Solving                                                                       [�[33m247.58 s�[39m] [�[33m 7053 MB�[39m]
Aborting as solve did not converge

Solve failed, cutting timestep.

Time Step 1, time = 1.25e-08, dt = 1.25e-08

    Computing Residual....                                                               [�[33m 26.29 s�[39m] [�[33m 7054 MB�[39m]
    |residual|_2 of individual variables:
                     disp_x: 213504
                     disp_y: 106835
                     eta0:   0
                     eta1:   0
 0 Nonlinear |R| = �[32m2.387418e+05�[39m
    Computing Jacobian.....                                                              [�[33m 34.36 s�[39m] [�[33m 7054 MB�[39m]
      0 Linear |R| = �[32m2.387418e+05�[39m

    Computing Residual....                                                               [�[33m 29.51 s�[39m] [�[33m 8863 MB�[39m]
      1 Linear |R| = �[32m4.238450e+04�[39m
    Computing Residual.....                                                              [�[33m 31.05 s�[39m] [�[33m 8863 MB�[39m]
      2 Linear |R| = �[32m1.419121e+04�[39m
    Computing Residual.....                                                              [�[33m 32.76 s�[39m] [�[33m 8918 MB�[39m]
      3 Linear |R| = �[32m2.819612e+02�[39m
    Computing Residual.....                                                              [�[33m 33.19 s�[39m] [�[33m 8918 MB�[39m]
      4 Linear |R| = �[32m3.681913e-01�[39m
    Computing Residual.....                                                              [�[33m 32.13 s�[39m] [�[33m 8918 MB�[39m]

ComputeMultipleCrystalPlasticityStress_abs: Constitutive failure
To recover, the solution will fail and then be re-attempted with a reduced time step.

    |residual|_2 of individual variables:
                     disp_x: 213504
                     disp_y: 106835
                     eta0:   0
                     eta1:   0
 1 Nonlinear |R| = �[33m2.387418e+05�[39m
  Nonlinear solve did not converge due to DIVERGED_FNORM_NAN iterations 1
�[31m Solve Did NOT Converge!�[39m
  Finished Solving                                                                       [�[33m246.54 s�[39m] [�[33m 7346 MB�[39m]
Aborting as solve did not converge

Solve failed, cutting timestep.

Time Step 1, time = 6.25e-09, dt = 6.25e-09

    Computing Residual....                                                               [�[33m 27.54 s�[39m] [�[33m 7346 MB�[39m]
    |residual|_2 of individual variables:
                     disp_x: 213504
                     disp_y: 106835
                     eta0:   0
                     eta1:   0
 0 Nonlinear |R| = �[32m2.387418e+05�[39m
    Computing Jacobian......                                                             [�[33m 35.78 s�[39m] [�[33m 7346 MB�[39m]
      0 Linear |R| = �[32m2.387418e+05�[39m

    Computing Residual....                                                               [�[33m 29.61 s�[39m] [�[33m 8871 MB�[39m]
      1 Linear |R| = �[32m4.238450e+04�[39m
    Computing Residual.....                                                              [�[33m 30.83 s�[39m] [�[33m 8871 MB�[39m]
      2 Linear |R| = �[32m1.419121e+04�[39m
    Computing Residual.....                                                              [�[33m 30.68 s�[39m] [�[33m 8926 MB�[39m]
      3 Linear |R| = �[32m2.819610e+02�[39m
    Computing Residual.....                                                              [�[33m 30.20 s�[39m] [�[33m 8926 MB�[39m]
      4 Linear |R| = �[32m3.680848e-01�[39m
    Computing Residual.....                                                              [�[33m 31.54 s�[39m] [�[33m 8926 MB�[39m]

ComputeMultipleCrystalPlasticityStress_abs: Constitutive failure
To recover, the solution will fail and then be re-attempted with a reduced time step.

    |residual|_2 of individual variables:
                     disp_x: 213504
                     disp_y: 106835
                     eta0:   0
                     eta1:   0
 1 Nonlinear |R| = �[33m2.387418e+05�[39m
  Nonlinear solve did not converge due to DIVERGED_FNORM_NAN iterations 1
�[31m Solve Did NOT Converge!�[39m
  Finished Solving                                                                       [�[33m239.30 s�[39m] [�[33m 7556 MB�[39m]
Aborting as solve did not converge

Solve failed, cutting timestep.

Time Step 1, time = 3.125e-09, dt = 3.125e-09

    Computing Residual....                                                               [�[33m 27.87 s�[39m] [�[33m 7556 MB�[39m]
    |residual|_2 of individual variables:
                     disp_x: 213504
                     disp_y: 106835
                     eta0:   0
                     eta1:   0
 0 Nonlinear |R| = �[32m2.387418e+05�[39m
    Computing Jacobian......                                                             [�[33m 36.35 s�[39m] [�[33m 7556 MB�[39m]
      0 Linear |R| = �[32m2.387418e+05�[39m

    Computing Residual.....                                                              [�[33m 31.70 s�[39m] [�[33m 9040 MB�[39m]
      1 Linear |R| = �[32m4.238451e+04�[39m
    Computing Residual.....                                                              [�[33m 31.55 s�[39m] [�[33m 9040 MB�[39m]
      2 Linear |R| = �[32m1.419122e+04�[39m
    Computing Residual.....                                                              [�[33m 31.67 s�[39m] [�[33m 9040 MB�[39m]
      3 Linear |R| = �[32m2.819611e+02�[39m
    Computing Residual...

Below is my input script.

[Mesh]
  [file]
     type = FileMeshGenerator
     file = Two_phase_noload_out.e-s183
     use_for_exodus_restart = true
     
   []
 []
[Problem]
  allow_initial_conditions_with_restart = true
[]
 [Variables]
  [./disp_x]
  [../]
  [./disp_y]
  [../]
   
  
  # order parameter 0
  [./eta0]
    initial_from_file_var = eta1
  [../]
  # order parameter 1
  [./eta1]
    initial_from_file_var = eta3
  [../]
  []
  [ICs]
    [./eta1]
      variable = eta1
      type = RandomIC
      min = -0.1625
      max = 0.1625
      seed = 192
    [../]
  []
[AuxVariables]
  [./vonmises]
    order = CONSTANT
    family = MONOMIAL
  [../]
  [./stress_xx]
    order = CONSTANT
    family = MONOMIAL
    block = 0
  [../]
  [temperature]
    order = FIRST
    family = LAGRANGE
  []
  [eth_xx]
    order = CONSTANT
    family = MONOMIAL
  []
  [eth_yy]
    order = CONSTANT
    family = MONOMIAL
  []
  
  [fth_xx]
    order = CONSTANT
    family = MONOMIAL
  []
  [fth_yy]
    order = CONSTANT
    family = MONOMIAL
  []
  
[]


[AuxKernels]
  [./vonmises]
    type = RankTwoScalarAux
    rank_two_tensor = stress
    variable = vonmises
    scalar_type = VonMisesStress
    execute_on = timestep_end
#   block = 0
  [../]
  [./stress_xx]
    type = RankTwoAux
    rank_two_tensor = stress
    variable = stress_xx
    index_j = 0
    index_i = 0
    execute_on = timestep_end
    block = 0
  [../]
  [temperature]
    type = FunctionAux
    variable = temperature
    function = '300' # temperature increases at a constant rate
    execute_on = timestep_begin
  []
  [eth_xx]
    type = RankTwoAux
    variable = eth_xx
    rank_two_tensor =  phase0_thermal_eigenstrain
    index_j = 0
    index_i = 0
    execute_on = timestep_end
  []
  [eth_yy]
    type = RankTwoAux
    variable = eth_yy
    rank_two_tensor = phase0_thermal_eigenstrain
    index_j = 1
    index_i = 1
    execute_on = timestep_end
  []
  
  [fth_xx]
    type = RankTwoAux
    variable = fth_xx
    rank_two_tensor = phase0_thermal_deformation_gradient
    index_j = 0
    index_i = 0
    execute_on = timestep_end
  []
  [fth_yy]
    type = RankTwoAux
    variable = fth_yy
    rank_two_tensor = phase0_thermal_deformation_gradient
    index_j = 1
    index_i = 1
    execute_on = timestep_end
  []
  
[]

[BCs]
  [symmy]
    type = DirichletBC
    variable = disp_y
    boundary = bottom
    value = 0
  []
  [symmx]
    type = DirichletBC
    variable = disp_x
    boundary = left
    value = 0
  []
  
  [tdisp]
    type = FunctionDirichletBC
    variable = disp_x
    boundary = right
    function = '0.1*t'
  []
[]

[Materials]
  [elasticity_tensor_phase0]
    type = ComputeElasticityTensorCP
    C_ijkl = '2.906e5 1.87e5 1.87e5 2.906e5 1.87e5 2.906e5 1.142e5 1.142e5 1.142e5'
    fill_method = symmetric9
    base_name = phase0
  []
  [stress_phase0]
    type = ComputeMultipleCrystalPlasticityStress_abs
    crystal_plasticity_models = 'trial_xtalpl_phase0'
    eigenstrain_names = 'thermal_eigenstrain'
    tan_mod_type = exact
    rtol = 1e-08
    base_name = phase0
  []
  
  [trial_xtalpl_phase0]
    type = CrystalPlasticityKalidindiUpdate
    number_slip_systems = 12
    slip_sys_file_name = input_slip_sys.txt
    crystal_lattice_type = FCC
    resistance_tol = 0.01
    r = 1.0             
    h = 6000            
    t_sat = 598.5        
    gss_a = 1.5         
    ao = 0.001           
    xm = 0.017             
    gss_initial = 400
    base_name = phase0
  []
  [thermal_eigenstrain]
    type = ComputeCrystalPlasticityThermalEigenstrain
    eigenstrain_name = thermal_eigenstrain
    deformation_gradient_name = thermal_deformation_gradient
    temperature = temperature
    thermal_expansion_coefficients = '12.8e-06 12.8e-06 12.8e-06'
    base_name = phase0
  []
  
  [./strain_phase0]
    type = ComputeFiniteStrain
    displacements = 'disp_x disp_y'
    base_name = phase0
   # eigenstrain_names = 'eigenstrain_phase0'
  [../]

  
  [elasticity_tensor_phase1]
    type = ComputeElasticityTensorCP
    C_ijkl = '2.721e5 1.69e5 1.69e5 2.721e5 1.69e5 2.721e5 1.31e5 1.31e5 1.31e5'
    fill_method = symmetric9
    base_name = phase1
  []
  [stress_phase1]
    type = ComputeMultipleCrystalPlasticityStress_abs
    crystal_plasticity_models = 'trial_xtalpl_phase1'
    tan_mod_type = exact
    rtol = 1e-08
    base_name = phase1
  []
  [trial_xtalpl_phase1]
    type = CrystalPlasticityKalidindiUpdate
    number_slip_systems = 12
    slip_sys_file_name = input_slip_sys.txt
    crystal_lattice_type = FCC
    resistance_tol = 0.01
    r = 1.0             
    h = 6000            
    t_sat = 598.5        
    gss_a = 1.5         
    ao = 0.001           
    xm = 0.017             
    gss_initial = 316 
    base_name = phase1
  [] 
  [./strain_phase1]
    type = ComputeFiniteStrain
    displacements = 'disp_x disp_y'
    base_name = phase1
  [../]
 
   # Switching functions for each phase
   [./h0]
    type = SwitchingFunctionMultiPhaseMaterial
    phase_etas = eta0
    all_etas = 'eta0 eta1'
    h_name = h0
  [../]
  [./h1]
    type = SwitchingFunctionMultiPhaseMaterial
    phase_etas = eta1
    all_etas = 'eta0 eta1'
    h_name = h1
  [../]
 
  # Generate the global stress from the phase stresses
  [./global_stress]
    type = MultiPhaseStressMaterial
    phase_base = 'phase0 phase1'
    h          = 'h0     h1'
  [../]
[]

[Kernels]
  [./eta0_dt]
    type = TimeDerivative
    variable = eta0
  [../]
  [./eta1_dt]
    type = TimeDerivative
    variable = eta1
  [../]
  [./TensorMechanics]
    displacements = 'disp_x disp_y'
    strain = FINITE
    incremental = true
    add_variables = true
    generate_output = stress_xx
  [../]
[]

[Postprocessors]
  [stress_xx]
    type = ElementAverageValue
    variable = stress_xx
  []
  [eth_xx]
    type = ElementAverageValue
    variable = eth_xx
  []
  [eth_yy]
    type = ElementAverageValue
    variable = eth_yy
  []
  
  [fth_xx]
    type = ElementAverageValue
    variable = fth_xx
  []
  [fth_yy]
    type = ElementAverageValue
    variable = fth_yy
  []
  
  [temperature]
    type = ElementAverageValue
    variable = temperature
  []
[]

[Preconditioning]
  [smp]
    type = SMP
    full = true
  []
[]

[Executioner]
  type = Transient
  solve_type = 'PJFNK'

  petsc_options_iname = '-pc_type -pc_factor_mat_solver_package -snes_type'
  petsc_options_value = 'lu            superlu_dist          vinewtonrsls'  

#  petsc_options_iname = '-pc_type -pc_asm_overlap -sub_pc_type -ksp_type -ksp_gmres_restart'
#  petsc_options_value = ' asm      2              lu            gmres     200'
  l_max_its = 20
  nl_max_its = 10
  nl_rel_tol = 1.0e-8
  nl_abs_tol = 1.0e-9

  end_time = 100

  [./TimeStepper]
    type = IterationAdaptiveDT
    dt = 5e-8
    cutback_factor = 0.75
    growth_factor = 1.2
    optimal_iterations = 20
  [../]
  
    [./Adaptivity]
      initial_adaptivity = 1
      refine_fraction = 0.6
      coarsen_fraction = 0.1
      max_h_level = 1
    [../]
  

[]

[Outputs]
  csv = true
  exodus = true
  print_linear_residuals = true
  [console]
    type = Console
    max_rows = 5
  []
[]

[Debug]
  show_var_residual_norms = true
[]

[GiudGiud]

Hello

This is not a bug you dont need to use a debug version.
it s a physics model issue. see this failure in the log

ComputeMultipleCrystalPlasticityStress_abs: Constitutive failure
To recover, the solution will fail and then be re-attempted with a reduced time step.