Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Excutioner settings to solve elastodynamic explicitly #157

Open
lyyc199586 opened this issue Jan 26, 2024 · 2 comments
Open

Excutioner settings to solve elastodynamic explicitly #157

lyyc199586 opened this issue Jan 26, 2024 · 2 comments

Comments

@lyyc199586
Copy link
Contributor

lyyc199586 commented Jan 26, 2024

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.

@hugary1995
Copy link
Owner

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

@lyyc199586
Copy link
Contributor Author

Now at idaholab/moose#26653

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants