exaFOAM/foam-extend4.1-coherent-io
Archived
1
1
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity

Reverting to old (correct) consistency faceU calculation. Vuko Vukcevic.

Browse source
This commit is contained in:
Hrvoje Jasak 2019-07-18 15:01:46 +01:00
commit 1f5b9fa132
1 changed files with 7 additions and 7 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

14
src/finiteVolume/cfdTools/general/solutionControl/solutionControl/solutionControl.C
View file

@ -1019,19 +1019,12 @@ void Foam::solutionControl::reconstructTransientVelocity
const volScalarField& p const volScalarField& p
) const ) const
{ {
// If the mesh is moving, flux needs to be relative before boundary
// conditions for velocity are corrected. VV and IG, 4/Jan/2016.
fvc::makeRelative(phi, U);
// Reconstruct the velocity using all the components from original equation // Reconstruct the velocity using all the components from original equation
U = 1.0/(1.0/rAU + ddtUEqn.A())* U = 1.0/(1.0/rAU + ddtUEqn.A())*
( (
U/rAU + ddtUEqn.H() - fvc::grad(p) U/rAU + ddtUEqn.H() - fvc::grad(p)
); );
// Correct boundary conditions with relative flux
U.correctBoundaryConditions();
// Get name and the corresponding index // Get name and the corresponding index
const word& UName = U.name(); const word& UName = U.name();
const label i = indices_[UName]; const label i = indices_[UName];
@ -1074,6 +1067,13 @@ void Foam::solutionControl::reconstructTransientVelocity
// Now that the normal component is zero, add the normal component from // Now that the normal component is zero, add the normal component from
// conservative flux // conservative flux
faceU += phi*rSf; faceU += phi*rSf;
// If the mesh is moving, flux needs to be relative before boundary
// conditions for velocity are corrected. VV and IG, 4/Jan/2016.
fvc::makeRelative(phi, U);
// Correct boundary conditions with relative flux
U.correctBoundaryConditions();
} }