atgeirr
commented
2 years ago The convergence stalls, with a pattern like this:
Iter MB(O) MB(W) MB(G) MB(E) CNV(O) CNV(W) CNV(G) CNV(E)
0 0.000e+00 2.726e-19 0.000e+00 2.312e-18 0.000e+00 1.337e+01 0.000e+00 1.910e+02
1 0.000e+00 2.650e-09 0.000e+00 1.180e-01 0.000e+00 3.131e-02 0.000e+00 5.998e-01
2 0.000e+00 2.909e-09 0.000e+00 1.180e-01 0.000e+00 1.472e-03 0.000e+00 1.182e-01
3 0.000e+00 3.228e-09 0.000e+00 1.180e-01 0.000e+00 8.189e-04 0.000e+00 1.182e-01
4 0.000e+00 3.457e-09 0.000e+00 1.180e-01 0.000e+00 5.205e-04 0.000e+00 1.182e-01... no matter how much the timestep is cut. It could be that we have some error in scaling of the energy residual (have not looked at this in detail now), and that this is actually ok. Probably not, but you could test by increasing the CNV tolerance to something larger than 1.182e-1 and see if the solution seems ok? (That would be a far too coarse CNV tolerance for the fluids though. So only recommending this as a test.)
EdmundStephens
OPMUSER
Also the grid cells are non-orthogonal, which is never a good idea.
TEMP2 is more reasonable as there is a temperature gradient, no temperature discontinuity.
goncalvesmachadoc
[v2022.10] We are running some shallow geothermal evaluations, for which we have a layer subcropping to the sea bed and a geothermal gradient from 40 F (sea water) down to about 80 F at 600m depth. I am trying to use the TEMPI initial temperature to initialise the model, but this is triggering some sort of difficulty ...
Attached example runs:
When TEMPI is not a constant, the solver backup up in the first timestep:
Unclear what is the issue here ... might be something in the input, or something in Flow. I'd be very grateful if you could have a check in the debugger. Unfortunately in this case it is not possible to check in Eclipse as TEMPI is an E300-only keyword.
TEMPI1.zip TEMPI2.zip TEMPI3.zip