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

driftFluxFoam::relativeVelocityModel: Updated dispersed-phase flow acceleration estimate

Browse Source 
// Dispersed phase velocity
    // const volVectorField Ud(mixture_.U() + Udm_);

    // Use the mixture rather than the dispersed-phase velocity to approximate
    // the dispersed-phase acceleration to improve stability as only the mixture
    // momentum equation is coupled to continuity and pressure
    //
    // This approximation is valid only in the limit of small drift-velocity.
    // For large drift-velocity an Euler-Euler approach should be used in
    // which both the continuous and dispersed-phase momentum equations are
    // solved and coupled to the pressure.
    const volVectorField& Ud = mixture_.U();

    return g_ - (Ud & fvc::grad(Ud));
This commit is contained in:
Henry Weller
2022-05-17 14:01:59 +01:00
parent 0ead4cbc81
commit 583fc4aec0
1 changed files with 12 additions and 1 deletions
Download Patch File Download Diff File Expand all files Collapse all files

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

@ -144,7 +144,18 @@ Foam::relativeVelocityModel::~relativeVelocityModel()
Foam::tmp<Foam::volVectorField>
Foam::relativeVelocityModel::acceleration() const
{
const volVectorField Ud(mixture_.U() + Udm_);
// Dispersed phase velocity
// const volVectorField Ud(mixture_.U() + Udm_);
// Use the mixture rather than the dispersed-phase velocity to approximate
// the dispersed-phase acceleration to improve stability as only the mixture
// momentum equation is coupled to continuity and pressure
//
// This approximation is valid only in the limit of small drift-velocity.
// For large drift-velocity an Euler-Euler approach should be used in
// which both the continuous and dispersed-phase momentum equations are
// solved and coupled to the pressure.
const volVectorField& Ud = mixture_.U();
return g_ - (Ud & fvc::grad(Ud));
}