Bugfix: removed SuSp stabilisation of the momentum equation

This commit is contained in:
Hrvoje Jasak 2017-02-06 11:45:17 +00:00
parent 13f27d9376
commit e98722db65
2 changed files with 2 additions and 6 deletions

View file

@ -52,9 +52,7 @@
trAU().rename("rAU");
}
// Insert the additional components. Note this will destroy the H and A
UEqn += fvm::SuSp(-divPhi, U) + divPhi*U;
// Under-relax momentum. Note this will destroy the H and A
UEqn.relax();
// Insert momentum equation

View file

@ -14,9 +14,7 @@
rAU = 1.0/UEqn.A();
// Insert the additional components. Note this will destroy the H and A
UEqn += fvm::SuSp(-divPhi, U) + divPhi*U;
// Under-relax momentum. Note this will destroy the H and A
UEqn.relax();
UpEqn.insertEquation(0, UEqn);