ERROR : IOobject: volScalarField alpha.phase1 not stored. Use field.storePrevIter() at start of iteration.
In version v1806, the following bug has appeared in interIsoFoam. We have a test case corresponding to a single rising bubble.
The error is :
--> FOAM FATAL ERROR: previous iteration field IOobject: volScalarField alpha.phase1 readOpt: 0 writeOpt: 0 globalObject: 0 "/work/islin462_1/gametl/OpenFOAM/RUNS/HysingBenchmark_2D/RESULTS_interIsoFoam/160x320x001/0"
not stored. Use field.storePrevIter() at start of iteration.
From function const Foam::GeometricField<Type, PatchField, GeoMesh>& Foam::GeometricField<Type, PatchField, GeoMesh>::prevIter() const [with Type = double; PatchField = Foam::fvPatchField; GeoMesh = Foam::volMesh]
in file /home/issrvhome1/R15/gametl/work7/OpenFOAM/OpenFOAM-plus/src/OpenFOAM/lnInclude/GeometricField.C at line 900.
FOAM aborting
#0 Foam::error::printStack(Foam::Ostream&) at ??:? #1 (closed) Foam::error::abort() at ??:? #2 (closed) Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::prevIter() const at ??:? #3 (closed) ? at ??:? #4 (closed) __libc_start_main in /lib64/libc.so.6 #5 (closed) ? at ??:?
The problem might be related to $FOAM_APP/solvers/multiphase/interIsoFoam/alphaEqn.H line 14 where storePrevIter should be used instead of prevIter. Thanks in advance for checking this.
diff --git a/applications/solvers/multiphase/interIsoFoam/alphaEqn.H b/applications/solvers/multiphase/interIsoFoam/alphaEqn.H index c79afea..014f4aa 100644 --- a/applications/solvers/multiphase/interIsoFoam/alphaEqn.H +++ b/applications/solvers/multiphase/interIsoFoam/alphaEqn.H @@ -11,7 +11,7 @@ // To recalculate the alpha1 update in subsequent iterations, we // must store its current value before overwriting with the new // value
-
alpha1.prevIter();
-
alpha1.storePrevIter(); // Storing initial phi and U for use in later outer iterations. phi.storePrevIter(); U.storePrevIter();
Best regards
Lionel GAMET