Skip to content
Snippets Groups Projects
Commit 7322f73a authored by Andrew Heather's avatar Andrew Heather Committed by Andrew Heather
Browse files

Merge branch 'feature-isoAdvectorWithMorphingMeshes' into 'develop'

Introduced changes required to make isoAdvector and interIsoFoam work with...

See merge request OpenFOAM-plus!257
parents afb9c1d9 61c87b55
Branches
Tags
No related merge requests found
{ {
// Temporarily making U relative to mesh motion
if (mesh.moving())
{
U -= fvc::reconstruct(mesh.phi());
}
// Updating alpha1 // Updating alpha1
advector.advect(); advector.advect();
// Making U absolute again after advection step
if (mesh.moving())
{
U += fvc::reconstruct(mesh.phi());
}
#include "rhofs.H" #include "rhofs.H"
rhoPhi = advector.getRhoPhi(rho1f, rho2f); rhoPhi = advector.getRhoPhi(rho1f, rho2f);
......
...@@ -6,7 +6,7 @@ ...@@ -6,7 +6,7 @@
\\/ M anipulation | \\/ M anipulation |
------------------------------------------------------------------------------- -------------------------------------------------------------------------------
isoAdvector | Copyright (C) 2016-2017 DHI isoAdvector | Copyright (C) 2016-2017 DHI
Modified work | Copyright (C) 2018 Johan Roenby Modified work | Copyright (C) 2019 Johan Roenby
------------------------------------------------------------------------------- -------------------------------------------------------------------------------
License License
...@@ -852,6 +852,12 @@ void Foam::isoAdvection::advect() ...@@ -852,6 +852,12 @@ void Foam::isoAdvection::advect()
// Do the isoAdvection on surface cells // Do the isoAdvection on surface cells
timeIntegratedFlux(); timeIntegratedFlux();
// Adjust alpha for mesh motion
if (mesh_.moving())
{
alpha1In_ *= (mesh_.Vsc0()/mesh_.Vsc());
}
// Adjust dVf for unbounded cells // Adjust dVf for unbounded cells
limitFluxes(); limitFluxes();
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment