zhyang-dev/OpenFOAM-12
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

interFoam,compressibleInterFoam,multiphaseEulerFoam: Ensure divU is registered for mapping

Browse Source 
divU is cached before mesh-motion and mapped post-motion to drive the optional pcorr flux
update to ensure the fluxes are conservative following mesh motion and change.
This commit is contained in:
Henry Weller
2021-04-15 10:44:04 +01:00
parent e5dd6117b2
commit 3de28f8736
3 changed files with 6 additions and 3 deletions
Download Patch File Download Diff File Expand all files Collapse all files

3
applications/solvers/multiphase/compressibleInterFoam/compressibleInterFoam.C
View File

@ -106,7 +106,8 @@ int main(int argc, char *argv[])
if (correctPhi)
{
divU = volScalarField::New
// Construct and register divU for mapping
divU = new volScalarField
(
"divU0",
fvc::div(fvc::absolute(phi, U))

3
applications/solvers/multiphase/interFoam/interFoam.C
View File

@ -111,7 +111,8 @@ int main(int argc, char *argv[])
&& !isType<twoPhaseChangeModels::noPhaseChange>(phaseChange)
)
{
divU = volScalarField::New
// Construct and register divU for mapping
divU = new volScalarField
(
"divU0",
fvc::div(fvc::absolute(phi, U))

3
applications/solvers/multiphase/multiphaseEulerFoam/multiphaseEulerFoam/multiphaseEulerFoam.C
View File

@ -136,7 +136,8 @@ int main(int argc, char *argv[])
correctPhi
)
{
divU = volScalarField::New
// Construct and register divU for mapping
divU = new volScalarField
(
"divU0",
fvc::div