GiudGiud
commented
4 months ago only converges with SVD for me. Weird that your log reports SVD too. For LU it seems we have got to find a way around the saddle point. Field split is a good idea.
Framework Information:
MOOSE Version: git commit f155eb1de0 on 2024-02-07
LibMesh Version:
PETSc Version: 3.20.3
SLEPc Version: 3.20.1
Current Time: Fri Feb 23 08:57:56 2024
Executable Timestamp: Fri Feb 23 08:57:33 2024
Parallelism:
Num Processors: 1
Num Threads: 1
Mesh:
Parallel Type: replicated
Mesh Dimension: 2
Spatial Dimension: 2
Nodes: 231
Elems: 50
Num Subdomains: 1
Nonlinear System:
Num DOFs: 528
Num Local DOFs: 528
Variables: "velocity" "p"
Finite Element Types: "LAGRANGE_VEC" "LAGRANGE"
Approximation Orders: "SECOND" "FIRST"
Auxiliary System:
Num DOFs: 132
Num Local DOFs: 132
Variables: { "vel_x" "vel_y" }
Finite Element Types: "LAGRANGE"
Approximation Orders: "FIRST"
Execution Information:
Executioner: Transient
TimeStepper: ConstantDT
TimeIntegrator: ImplicitEuler
Solver Mode: NEWTON
PETSc Preconditioner: svd
MOOSE Preconditioner: SMP
Time Step 0, time = 0
Time Step 1, time = 3.1536e+06, dt = 3.1536e+06
0 Nonlinear |R| = 1.379140e+00
SVD: condition number 2.606499026776e+17, 0 of 528 singular values are (nearly) zero
SVD: smallest singular values: 2.195040813607e-12 2.530633109456e-12 2.615343334311e-12 2.894437182073e-12 3.350991855258e-12
SVD: largest singular values : 5.683357729119e+05 5.707121183379e+05 5.707121185051e+05 5.721371742527e+05 5.721371744402e+05
0 Linear |R| = 1.379140e+00
1 Linear |R| = 1.561705e-04
2 Linear |R| = 2.384225e-10
1 Nonlinear |R| = 2.903115e-01
SVD: condition number 2.604980373976e+17, 0 of 528 singular values are (nearly) zero
SVD: smallest singular values: 2.195036259922e-12 2.531284041482e-12 2.756510373859e-12 3.258149672264e-12 3.510944603970e-12
SVD: largest singular values : 5.653752781411e+05 5.693788673245e+05 5.693788674770e+05 5.718026375429e+05 5.718026377264e+05
0 Linear |R| = 2.903115e-01
1 Linear |R| = 8.757300e-06
2 Linear |R| = 3.709858e-12
2 Nonlinear |R| = 1.875023e-01
SVD: condition number 2.603920693564e+17, 0 of 528 singular values are (nearly) zero
SVD: smallest singular values: 2.195041227021e-12 2.531440905243e-12 2.801797642962e-12 3.359261010733e-12 3.585186601310e-12
SVD: largest singular values : 5.633758897990e+05 5.684621465368e+05 5.684621466793e+05 5.715713272458e+05 5.715713274265e+05
0 Linear |R| = 1.875023e-01
1 Linear |R| = 1.640546e-06
3 Nonlinear |R| = 1.314256e-01
SVD: condition number 2.603280438807e+17, 0 of 528 singular values are (nearly) zero
SVD: smallest singular values: 2.195029882786e-12 2.531549754907e-12 2.830355619618e-12 3.397564421603e-12 3.655003677263e-12
SVD: largest singular values : 5.620622803166e+05 5.678814582849e+05 5.678814584172e+05 5.714278354673e+05 5.714278356453e+05
0 Linear |R| = 1.314256e-01
1 Linear |R| = 3.598852e-07
4 Nonlinear |R| = 9.242926e-02
SVD: condition number 2.602838352382e+17, 0 of 528 singular values are (nearly) zero
SVD: smallest singular values: 2.195036426968e-12 2.531640226087e-12 2.852728682074e-12 3.416459039081e-12 3.712867607931e-12
SVD: largest singular values : 5.614043428864e+05 5.675302304444e+05 5.675302305751e+05 5.713324995213e+05 5.713324996987e+05
0 Linear |R| = 9.242926e-02
1 Linear |R| = 1.668183e-07
5 Nonlinear |R| = 6.394064e-02
SVD: condition number 2.602634396391e+17, 0 of 528 singular values are (nearly) zero
SVD: smallest singular values: 2.195044276215e-12 2.531668181939e-12 2.864802025031e-12 3.424864010039e-12 3.748604785635e-12
SVD: largest singular values : 5.610564135954e+05 5.673623024893e+05 5.673623026190e+05 5.712897733107e+05 5.712897734878e+05
0 Linear |R| = 6.394064e-02
1 Linear |R| = 1.727536e-07
6 Nonlinear |R| = 4.366106e-02
SVD: condition number 2.602311157190e+17, 0 of 528 singular values are (nearly) zero
SVD: smallest singular values: 2.195036768602e-12 2.531726406256e-12 2.872062553037e-12 3.429377209786e-12 3.771178954683e-12
SVD: largest singular values : 5.604155129635e+05 5.670687703345e+05 5.670687704622e+05 5.712168671611e+05 5.712168673376e+05
0 Linear |R| = 4.366106e-02
1 Linear |R| = 9.937027e-08
7 Nonlinear |R| = 2.968398e-02
SVD: condition number 2.602127231219e+17, 0 of 528 singular values are (nearly) zero
SVD: smallest singular values: 2.195028511643e-12 2.531737663486e-12 2.876602300732e-12 3.432026758850e-12 3.785735810138e-12
SVD: largest singular values : 5.599882131494e+05 5.668898899128e+05 5.668898900380e+05 5.711743461690e+05 5.711743463447e+05
0 Linear |R| = 2.968398e-02
1 Linear |R| = 5.105166e-08
8 Nonlinear |R| = 2.007425e-02
SVD: condition number 2.601990935315e+17, 0 of 528 singular values are (nearly) zero
SVD: smallest singular values: 2.195038493058e-12 2.531723589503e-12 2.879441852861e-12 3.433614541935e-12 3.795055674349e-12
SVD: largest singular values : 5.597047358694e+05 5.667738265226e+05 5.667738266455e+05 5.711470259852e+05 5.711470261604e+05
0 Linear |R| = 2.007425e-02
1 Linear |R| = 3.289593e-08
9 Nonlinear |R| = 1.351553e-02
SVD: condition number 2.601923778405e+17, 0 of 528 singular values are (nearly) zero
SVD: smallest singular values: 2.195038300463e-12 2.531759713096e-12 2.881226100512e-12 3.434613915513e-12 3.801038459157e-12
SVD: largest singular values : 5.595510634927e+05 5.667109365518e+05 5.667109366736e+05 5.711322346737e+05 5.711322348486e+05
0 Linear |R| = 1.351553e-02
1 Linear |R| = 1.300519e-08
10 Nonlinear |R| = 9.081960e-03
AdrienWehrle
./diuca-opt -i problems/iceslab_ad_2d_SI.igives
and the following till minimum timestep.
yet @GiudGiud mentioned it was converging on their machine in https://github.com/AdrienWehrle/diuca/pull/4 (and I'm using this version of the input file). Were you converging on a different version of
iceslab_ad_2d_SI.i? If not, it may be because of a different version of PETSC/MOOSE...