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

compressibleMultiphaseInterFoam: Added support for fvModels and fvConstraints

Browse Source
This commit is contained in:
Henry Weller
2022-02-08 16:25:10 +00:00
parent 83f0335151
commit 13b78294ff
8 changed files with 81 additions and 61 deletions
Download Patch File Download Diff File Expand all files Collapse all files

57
applications/solvers/multiphase/compressibleInterFoam/UEqn.H
View File

@ -1,34 +1,33 @@
fvVectorMatrix UEqn
fvVectorMatrix UEqn
(
fvm::ddt(rho, U) + fvm::div(rhoPhi, U) - fvm::Sp(contErr, U)
+ MRF.DDt(rho, U)
+ turbulence.divDevTau(U)
==
fvModels.source(rho, U)
);
UEqn.relax();
fvConstraints.constrain(UEqn);
if (pimple.momentumPredictor())
{
solve
(
fvm::ddt(rho, U) + fvm::div(rhoPhi, U)
- fvm::Sp(contErr, U)
+ MRF.DDt(rho, U)
+ turbulence.divDevTau(U)
UEqn
==
fvModels.source(rho, U)
fvc::reconstruct
(
(
mixture.surfaceTensionForce()
- ghf*fvc::snGrad(rho)
- fvc::snGrad(p_rgh)
) * mesh.magSf()
)
);
UEqn.relax();
fvConstraints.constrain(U);
fvConstraints.constrain(UEqn);
if (pimple.momentumPredictor())
{
solve
(
UEqn
==
fvc::reconstruct
(
(
mixture.surfaceTensionForce()
- ghf*fvc::snGrad(rho)
- fvc::snGrad(p_rgh)
) * mesh.magSf()
)
);
fvConstraints.constrain(U);
K = 0.5*magSqr(U);
}
K = 0.5*magSqr(U);
}

4
applications/solvers/multiphase/compressibleInterFoam/compressibleAlphaEqnSubCycle.H
View File

@ -61,7 +61,7 @@ volScalarField::Internal contErr
(
(
fvc::ddt(rho) + fvc::div(rhoPhi)
- (fvModels.source(alpha1, mixture.thermo1().rho())&rho1)
- (fvModels.source(alpha2, mixture.thermo2().rho())&rho2)
- (fvModels.source(alpha1, rho1)&rho1)
- (fvModels.source(alpha2, rho2)&rho2)
)()
);

2
applications/solvers/multiphase/compressibleMultiphaseInterFoam/Make/options
View File

@ -19,4 +19,6 @@ EXE_LIBS = \
-lmomentumTransportModels \
-lcompressibleMomentumTransportModels \
-lfiniteVolume \
-lfvModels \
-lfvConstraints \
-lmeshTools

17
applications/solvers/multiphase/compressibleMultiphaseInterFoam/TEqn.H
View File

@ -1,18 +1,25 @@
{
fvScalarMatrix TEqn
(
fvm::ddt(rho, T)
+ fvm::div(mixture.rhoPhi(), T)
fvm::ddt(rho, T) + fvm::div(mixture.rhoPhi(), T) - fvm::Sp(contErr, T)
- fvm::laplacian(mixture.alphaEff(rho*turbulence->nut()), T)
+ (
fvc::div(fvc::absolute(phi, U), p)
+ fvc::ddt(rho, K) + fvc::div(mixture.rhoPhi(), K)
)*mixture.rCv()
fvc::div(fvc::absolute(phi, U), p)()() // - contErr/rho*p
+ (fvc::ddt(rho, K) + fvc::div(mixture.rhoPhi(), K))()()
- (U()&(fvModels.source(rho, U)&U)()) - contErr*K
)*mixture.rCv()()
==
fvModels.source(rho, T)
);
TEqn.relax();
fvConstraints.constrain(TEqn);
TEqn.solve();
fvConstraints.constrain(T);
mixture.correctThermo();
mixture.correct();
}

51
applications/solvers/multiphase/compressibleMultiphaseInterFoam/UEqn.H
View File

@ -1,27 +1,32 @@
fvVectorMatrix UEqn
fvVectorMatrix UEqn
(
fvm::ddt(rho, U) + fvm::div(mixture.rhoPhi(), U) - fvm::Sp(contErr, U)
+ turbulence->divDevTau(U)
==
fvModels.source(rho, U)
);
UEqn.relax();
fvConstraints.constrain(UEqn);
if (pimple.momentumPredictor())
{
solve
(
fvm::ddt(rho, U)
+ fvm::div(mixture.rhoPhi(), U)
+ turbulence->divDevTau(U)
UEqn
==
fvc::reconstruct
(
(
mixture.surfaceTensionForce()
- ghf*fvc::snGrad(rho)
- fvc::snGrad(p_rgh)
) * mesh.magSf()
)
);
UEqn.relax();
fvConstraints.constrain(U);
if (pimple.momentumPredictor())
{
solve
(
UEqn
==
fvc::reconstruct
(
(
mixture.surfaceTensionForce()
- ghf*fvc::snGrad(rho)
- fvc::snGrad(p_rgh)
) * mesh.magSf()
)
);
K = 0.5*magSqr(U);
}
K = 0.5*magSqr(U);
}

6
applications/solvers/multiphase/compressibleMultiphaseInterFoam/compressibleMultiphaseInterFoam.C
View File

@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Copyright (C) 2013-2021 OpenFOAM Foundation
\\ / A nd | Copyright (C) 2013-2022 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@ -39,6 +39,8 @@ Description
#include "compressibleMultiphaseMixture.H"
#include "compressibleMomentumTransportModels.H"
#include "pimpleControl.H"
#include "fvModels.H"
#include "fvConstraints.H"
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
@ -80,7 +82,7 @@ int main(int argc, char *argv[])
{
mixture.solve();
solve(fvm::ddt(rho) + fvc::div(mixture.rhoPhi()));
#include "contErr.H"
#include "UEqn.H"
#include "TEqn.H"

3
applications/solvers/multiphase/compressibleMultiphaseInterFoam/createFields.H
View File

@ -57,3 +57,6 @@ autoPtr<compressible::momentumTransportModel> turbulence
Info<< "Creating field kinetic energy K\n" << endl;
volScalarField K("K", 0.5*magSqr(U));
#include "createFvModels.H"
#include "createFvConstraints.H"

2
applications/solvers/multiphase/compressibleMultiphaseInterFoam/pEqn.H
View File

@ -41,6 +41,7 @@
(
fvc::ddt(rho) + thermo.psi()*correction(fvm::ddt(p_rgh))
+ fvc::div(phi, rho) - fvc::Sp(fvc::div(phi), rho)
- (fvModels.source(phase(), rho)&rho)
).ptr()
);
@ -113,6 +114,7 @@
U = HbyA
+ rAU*fvc::reconstruct((phig + p_rghEqnIncomp.flux())/rAUf);
U.correctBoundaryConditions();
fvConstraints.constrain(U);
}
}