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

driftFluxFoam::mixtureViscosityModels: Added the basic Quemada viscosity model for for colloidal dispersions

Browse Source 
currently without strain-rate dependency.

Class
    Foam::mixtureViscosityModels::Quemada

Description
     Quemada viscosity model for for colloidal dispersions.

     References:
     \verbatim
         Quemada, D. (1998).
         Rheological modelling of complex fluids. I.
         The concept of effective volume fraction revisited.
         The European Physical Journal-Applied Physics, 1(1), 119-127.
    \endverbatim

Usage
    Example usage:
    \verbatim
    viscosityModel  Quemada;

    alphaMax    0.6; // Maximum dispersed phase-fraction (packing fraction)
    q           2;   // Exponent, defaults to 2

    rho         1996;
    \endverbatim
This commit is contained in:
Henry Weller
2022-01-26 15:25:30 +00:00
parent 4623ece721
commit f4b497bd81
6 changed files with 251 additions and 12 deletions
Download Patch File Download Diff File Expand all files Collapse all files

8
applications/solvers/multiphase/driftFluxFoam/mixtureViscosityModels/BinghamPlastic/BinghamPlastic.H
View File

@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Copyright (C) 2014-2021 OpenFOAM Foundation
\\ / A nd | Copyright (C) 2014-2022 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@ -79,7 +79,7 @@ public:
//- Destructor
~BinghamPlastic()
virtual ~BinghamPlastic()
{}
@ -87,14 +87,14 @@ public:
//- Return the mixture viscosity
// given the viscosity of the continuous phase
tmp<volScalarField> mu
virtual tmp<volScalarField> mu
(
const volScalarField& muc,
const volVectorField& U
) const;
//- Read phaseProperties dictionary
bool read();
virtual bool read();
};

1
applications/solvers/multiphase/driftFluxFoam/mixtureViscosityModels/Make/files
View File

@ -3,5 +3,6 @@ mixtureViscosityModel/mixtureViscosityModelNew.C
plastic/plastic.C
BinghamPlastic/BinghamPlastic.C
slurry/slurry.C
Quemada/Quemada.C
LIB = $(FOAM_LIBBIN)/libdriftFluxTransportModels

109
applications/solvers/multiphase/driftFluxFoam/mixtureViscosityModels/Quemada/Quemada.C Normal file
View File

@ -0,0 +1,109 @@
/*---------------------------------------------------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Copyright (C) 2022 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
This file is part of OpenFOAM.
OpenFOAM is free software: you can redistribute it and/or modify it
under the terms of the GNU General Public License as published by
the Free Software Foundation, either version 3 of the License, or
(at your option) any later version.
OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
for more details.
You should have received a copy of the GNU General Public License
along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
\*---------------------------------------------------------------------------*/
#include "Quemada.H"
#include "volFields.H"
#include "addToRunTimeSelectionTable.H"
// * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
namespace Foam
{
namespace mixtureViscosityModels
{
defineTypeNameAndDebug(Quemada, 0);
addToRunTimeSelectionTable
(
mixtureViscosityModel,
Quemada,
dictionary
);
}
}
// * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
Foam::mixtureViscosityModels::Quemada::Quemada
(
const fvMesh& mesh,
const word& group
)
:
mixtureViscosityModel(mesh, group),
alphaMax_
(
"alphaMax",
dimless,
optionalSubDict(typeName + "Coeffs").lookup("alphaMax")
),
q_(optionalSubDict(typeName + "Coeffs").lookupOrDefault("q", scalar(2))),
alpha_
(
mesh.lookupObject<volScalarField>
(
IOobject::groupName
(
lookupOrDefault<word>("alpha", "alpha"),
group
)
)
)
{}
// * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * * //
Foam::tmp<Foam::volScalarField>
Foam::mixtureViscosityModels::Quemada::mu
(
const volScalarField& muc,
const volVectorField& U
) const
{
return muc*pow(1.0 - alpha_/alphaMax_, -q_);
}
bool Foam::mixtureViscosityModels::Quemada::read()
{
if (mixtureViscosityModel::read())
{
const dictionary& dict = optionalSubDict(typeName + "Coeffs");
dict.lookup("alphaMax") >> alphaMax_;
dict.lookup("q") >> q_;
return true;
}
else
{
return false;
}
}
// ************************************************************************* //

129
applications/solvers/multiphase/driftFluxFoam/mixtureViscosityModels/Quemada/Quemada.H Normal file
View File

@ -0,0 +1,129 @@
/*---------------------------------------------------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Copyright (C) 2022 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
This file is part of OpenFOAM.
OpenFOAM is free software: you can redistribute it and/or modify it
under the terms of the GNU General Public License as published by
the Free Software Foundation, either version 3 of the License, or
(at your option) any later version.
OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
for more details.
You should have received a copy of the GNU General Public License
along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
Class
Foam::mixtureViscosityModels::Quemada
Description
Quemada viscosity model for for colloidal dispersions.
References:
\verbatim
Quemada, D. (1998).
Rheological modelling of complex fluids. I.
The concept of effective volume fraction revisited.
The European Physical Journal-Applied Physics, 1(1), 119-127.
\endverbatim
Usage
Example usage:
\verbatim
viscosityModel Quemada;
alphaMax 0.6; // Maximum dispersed phase-fraction (packing fraction)
q 2; // Exponent, defaults to 2
rho 1996;
\endverbatim
SourceFiles
Quemada.C
\*---------------------------------------------------------------------------*/
#ifndef Quemada_H
#define Quemada_H
#include "mixtureViscosityModel.H"
#include "dimensionedScalar.H"
#include "volFieldsFwd.H"
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
namespace Foam
{
namespace mixtureViscosityModels
{
/*---------------------------------------------------------------------------*\
Class Quemada Declaration
\*---------------------------------------------------------------------------*/
class Quemada
:
public mixtureViscosityModel
{
// Private data
//- Maximum dispersed phase-fraction (packing fraction)
dimensionedScalar alphaMax_;
//- Exponent (defaults to 2)
scalar q_;
//- Dispersed phase fraction
const volScalarField& alpha_;
public:
//- Runtime type information
TypeName("Quemada");
// Constructors
//- Construct from components
Quemada(const fvMesh& mesh, const word& group);
//- Destructor
virtual ~Quemada()
{}
// Member Functions
//- Return the mixture viscosity
// given the viscosity of the continuous phase
virtual tmp<volScalarField> mu
(
const volScalarField& muc,
const volVectorField& U
) const;
//- Read phaseProperties dictionary
virtual bool read();
};
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
} // End namespace mixtureViscosityModels
} // End namespace Foam
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
#endif
// ************************************************************************* //

8
applications/solvers/multiphase/driftFluxFoam/mixtureViscosityModels/plastic/plastic.H
View File

@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Copyright (C) 2014-2021 OpenFOAM Foundation
\\ / A nd | Copyright (C) 2014-2022 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@ -90,7 +90,7 @@ public:
//- Destructor
~plastic()
virtual ~plastic()
{}
@ -98,14 +98,14 @@ public:
//- Return the mixture viscosity
// given the viscosity of the continuous phase
tmp<volScalarField> mu
virtual tmp<volScalarField> mu
(
const volScalarField& muc,
const volVectorField& U
) const;
//- Read phaseProperties dictionary
bool read();
virtual bool read();
};

8
applications/solvers/multiphase/driftFluxFoam/mixtureViscosityModels/slurry/slurry.H
View File

@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Copyright (C) 2014-2021 OpenFOAM Foundation
\\ / A nd | Copyright (C) 2014-2022 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@ -85,7 +85,7 @@ public:
//- Destructor
~slurry()
virtual ~slurry()
{}
@ -93,14 +93,14 @@ public:
//- Return the mixture viscosity
// given the viscosity of the continuous phase
tmp<volScalarField> mu
virtual tmp<volScalarField> mu
(
const volScalarField& muc,
const volVectorField& U
) const;
//- Read phaseProperties dictionary
bool read();
virtual bool read();
};