November 20, 2023 at 5:55 pm
Rocco Bombardieri
Subscriber

Interestingly enough, reducing the time step by keyword *CONTROL_TIMESTEP  help achieve a good solution. Same soluction could be obtained using the implicit solver with the following settings. 

*CONTROL_ACCURACY
$#     osu       inn    pidosu      iacc    
         1         4         0         0
*CONTROL_IMPLICIT_AUTO
$#   iauto    iteopt    itewin     dtmin     dtmax     dtexp     kfail    kcycle
         1        30        10       0.0    1.0E-4       0.0         0         0
*CONTROL_IMPLICIT_DYNAMICS
$#   imass     gamma      beta    tdybir    tdydth    tdybur     irate     alpha
         1       0.6      0.38       0.01.00000E281.00000E28         0       0.0
*CONTROL_IMPLICIT_GENERAL
$#  imflag       dt0    imform      nsbs       igs     cnstn      form    zero_v
         1    1.0E-4         2         1         1         0         0         0
*CONTROL_IMPLICIT_SOLUTION
$#  nsolvr    ilimit    maxref     dctol     ectol     rctol     lstol    abstol
        12        15        15     0.001      0.011.00000E10       0.91.0000E-20
$#   dnorm    diverg     istif   nlprint    nlnorm   d3itctl     cpchk     
         1         1         1         1         0         1         0
$#  arcctl    arcdir    arclen    arcmth    arcdmp    arcpsi    arcalf    arctim
         0         0       0.0         1         2         0         0         0
$#   lsmtd     lsdir      irad      srad      awgt      sred    
         5         2       0.0       0.0       0.0       0.0