zhyang-dev/openfoam
1
0
Fork 0
mirror of https://develop.openfoam.com/Development/openfoam.git synced 2025-11-28 03:28:01 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

BUG: First occurence of alpha1.prevIter() causes crash (closes #911)

Browse Source 
- in interIsoFoam/alphaEqn.H changed to alpha1.prevIter() to
  alpha1.storePrevIter().
This commit is contained in:
Johan Roenby
2018-07-19 23:36:42 +02:00
committed by Mark Olesen
parent a6e735e5ce
commit 40618a3a86
1 changed files with 9 additions and 8 deletions
Download Patch File Download Diff File Expand all files Collapse all files

3
applications/solvers/multiphase/interIsoFoam/alphaEqn.H
View File

@ -11,7 +11,7 @@
// To recalculate the alpha1 update in subsequent iterations, we // To recalculate the alpha1 update in subsequent iterations, we
// must store its current value before overwriting with the new // must store its current value before overwriting with the new
// value // value
alpha1.prevIter(); alpha1.storePrevIter();
// Storing initial phi and U for use in later outer iterations. // Storing initial phi and U for use in later outer iterations.
phi.storePrevIter(); phi.storePrevIter();
U.storePrevIter(); U.storePrevIter();
@ -21,6 +21,7 @@
// Resetting alpha1 to value before advection in first PIMPLE // Resetting alpha1 to value before advection in first PIMPLE
// iteration. // iteration.
alpha1 = alpha1.prevIter(); alpha1 = alpha1.prevIter();
// Setting U and phi with which to advect interface. // Setting U and phi with which to advect interface.
U = 0.5*U.prevIter() + 0.5*U; U = 0.5*U.prevIter() + 0.5*U;
phi = 0.5*phi.prevIter() + 0.5*phi; phi = 0.5*phi.prevIter() + 0.5*phi;