raccoon icon indicating copy to clipboard operation
raccoon copied to clipboard

Excutioner settings to solve elastodynamic explicitly

Open lyyc199586 opened this issue 5 months ago • 2 comments

Bug Description

I am trying to solve elastodynamic part explicitly with Central difference time integrator, but getting DIVERGED_PC_FAILED errors:

Time Step 1, time = 2e-05, dt = 2e-05
Linear solve did not converge due to DIVERGED_PC_FAILED iterations 0
               PC failed due to FACTOR_NUMERIC_ZEROPIVOT 
 Solve Did NOT Converge!
Aborting as solve did not converge

Steps to Reproduce

For example, in tutorials for mode1_brittle_fracture, I change the Excutioner block to

[Executioner]
  type = Transient
  start_time = 0
  end_time = 3.5e-3
  dt = 2e-5
  [TimeIntegrator]
    type = CentralDifference
  []
[]

and change the stress divergence and add inertiaforce, with some density

[Kernels]
  [solid_x]
    type = ADDynamicStressDivergenceTensors
    variable = disp_x
    component = 0
  []
  [solid_y]
    type = ADDynamicStressDivergenceTensors
    variable = disp_y
    component = 1
    save_in = fy
  []
  [inertia_x]
    type = ADInertialForce
    variable = disp_x
  []
  [inertia_y]
    type = ADInertialForce
    variable = disp_y
  []
[]

Impact

I actually did the same thing for the elastodynamic solve in nano-pulsed simulation before and it works... I don't know why it cannot run now.

lyyc199586 avatar Jan 26 '24 18:01 lyyc199586

Can you reproduce this using a MOOSE test input file and open a bug report there?

hugary1995 avatar Jan 26 '24 18:01 hugary1995

Now at https://github.com/idaholab/moose/issues/26653

lyyc199586 avatar Jan 26 '24 19:01 lyyc199586