zhyang-dev/openfoam
1
0
Fork 0
mirror of https://develop.openfoam.com/Development/openfoam.git synced 2025-12-28 03:37:59 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

driftFluxFoam: Rationalisation of the handling of Udm and tauDm

Browse Source
This commit is contained in:
Henry
2014-03-05 16:04:10 +00:00
committed by Andrew Heather
parent 4acc7c1260
commit 09bb8b8f7f
13 changed files with 42 additions and 60 deletions
Download Patch File Download Diff File Expand all files Collapse all files

2
applications/solvers/multiphase/driftFluxFoam/UEqn.H
View File

@ -4,7 +4,7 @@
(
fvm::ddt(rho, U)
+ fvm::div(rhoPhi, U)
+ fvc::div(uRelModel.tau(), "div(phiUkm,Ukm)")
+ fvc::div(UdmModel.tauDm())
- fvm::laplacian(muEff, U)
- fvc::div(muEff*dev2(T(fvc::grad(U))))
==

2
applications/solvers/multiphase/driftFluxFoam/alphaEqnSubCycle.H
View File

@ -13,7 +13,7 @@
surfaceScalarField phir
(
mesh.Sf() & fvc::interpolate(uRelModel.Udm())
mesh.Sf() & fvc::interpolate(UdmModel.Udm())
);
if (nAlphaSubCycles > 1)

4
applications/solvers/multiphase/driftFluxFoam/createFields.H
View File

@ -86,7 +86,7 @@
// Relative Velocity
// ~~~~~~~~~~~~~~~~~
autoPtr<relativeVelocityModel> uRelModelPtr
autoPtr<relativeVelocityModel> UdmModelPtr
(
relativeVelocityModel::New
(
@ -95,7 +95,7 @@
)
);
relativeVelocityModel& uRelModel(uRelModelPtr());
relativeVelocityModel& UdmModel(UdmModelPtr());
// Turbulence

2
applications/solvers/multiphase/driftFluxFoam/driftFluxFoam.C
View File

@ -82,7 +82,7 @@ int main(int argc, char *argv[])
{
#include "alphaControls.H"
uRelModel.update();
UdmModel.correct();
#include "alphaEqnSubCycle.H"

8
applications/solvers/multiphase/driftFluxFoam/relativeVelocityModels/general/general.C
View File

@ -62,11 +62,11 @@ Foam::relativeVelocityModels::general::~general()
// * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * * //
Foam::tmp<Foam::volVectorField>
Foam::relativeVelocityModels::general::Ur() const
void Foam::relativeVelocityModels::general::correct()
{
return
V0_
Udm_ =
(rhoC_/rho())
*V0_
*(
exp(-a_*max(alphaD_ - residualAlpha_, scalar(0)))
- exp(-a1_*max(alphaD_ - residualAlpha_, scalar(0)))

4
applications/solvers/multiphase/driftFluxFoam/relativeVelocityModels/general/general.H
View File

@ -89,8 +89,8 @@ public:
// Member Functions
//- Relative velocity
virtual tmp<volVectorField> Ur() const;
//- Update the diffusion velocity
virtual void correct();
};

44
applications/solvers/multiphase/driftFluxFoam/relativeVelocityModels/relativeVelocityModel/relativeVelocityModel.C
View File

@ -85,23 +85,6 @@ Foam::relativeVelocityModel::relativeVelocityModel
alphaC_.mesh(),
dimensionedVector("Udm", dimVelocity, vector::zero),
mixture.U().boundaryField().types()
),
tau_
(
IOobject
(
"Udm",
alphaC_.time().timeName(),
alphaC_.mesh()
),
alphaC_.mesh(),
dimensionedSymmTensor
(
"Udm",
sqr(dimVelocity)*dimDensity,
symmTensor::zero
)
)
{}
@ -156,19 +139,28 @@ Foam::relativeVelocityModel::~relativeVelocityModel()
// * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * * //
void Foam::relativeVelocityModel::update()
tmp<volScalarField> Foam::relativeVelocityModel::rho() const
{
tmp<volVectorField> URel(Ur());
return alphaC_*rhoC_ + alphaD_*rhoD_;
}
tmp<volScalarField> betaC(alphaC_*rhoC_);
tmp<volScalarField> betaD(alphaD_*rhoD_);
tmp<volScalarField> rhoM(betaC() + betaD());
tmp<volVectorField> Udm = URel()*betaC()/rhoM;
tmp<volVectorField> Ucm = Udm() - URel;
tmp<volSymmTensorField> Foam::relativeVelocityModel::tauDm() const
{
volScalarField betaC(alphaC_*rhoC_);
volScalarField betaD(alphaD_*rhoD_);
Udm_ = Udm();
tau_ = betaD*sqr(Udm) + betaC*sqr(Ucm);
// Calculate the relative velocity of the continuous phase w.r.t the mean
volVectorField Ucm(betaD*Udm_/betaC);
return tmp<volSymmTensorField>
(
new volSymmTensorField
(
"tauDm",
betaD*sqr(Udm_) + betaC*sqr(Ucm)
)
);
}

18
applications/solvers/multiphase/driftFluxFoam/relativeVelocityModels/relativeVelocityModel/relativeVelocityModel.H
View File

@ -81,10 +81,7 @@ protected:
const dimensionedScalar& rhoD_;
//- Dispersed diffusion velocity
volVectorField Udm_;
//- Stress
volSymmTensorField tau_;
mutable volVectorField Udm_;
public:
@ -127,8 +124,8 @@ public:
// Member Functions
//- Calculate the relative velocity of the dispersed phase
virtual tmp<volVectorField> Ur() const = 0;
//- Return the mixture mean density
tmp<volScalarField> rho() const;
//- Return the diffusion velocity of the dispersed phase
const volVectorField& Udm() const
@ -137,13 +134,10 @@ public:
}
//- Return the stress tensor due to the phase transport
const volSymmTensorField& tau() const
{
return tau_;
}
tmp<volSymmTensorField> tauDm() const;
//- Update the stored diffusion velocity and stress
void update();
//- Update the diffusion velocity
virtual void correct() = 0;
};

8
applications/solvers/multiphase/driftFluxFoam/relativeVelocityModels/simple/simple.C
View File

@ -61,13 +61,9 @@ Foam::relativeVelocityModels::simple::~simple()
// * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * * //
Foam::tmp<Foam::volVectorField>
Foam::relativeVelocityModels::simple::Ur() const
void Foam::relativeVelocityModels::simple::correct()
{
return
V0_
*pow(scalar(10), -a_*max(alphaD_, scalar(0)))
/max(alphaC_, residualAlpha_);
Udm_ = (rhoC_/rho())*V0_*pow(scalar(10), -a_*max(alphaD_, scalar(0)));
}

4
applications/solvers/multiphase/driftFluxFoam/relativeVelocityModels/simple/simple.H
View File

@ -86,8 +86,8 @@ public:
// Member Functions
//- Relative velocity
virtual tmp<volVectorField> Ur() const;
//- Update the diffusion velocity
virtual void correct();
};