30 #include "twoPhaseSystem.H" 39 namespace diameterModels
54 const dictionary&
dict 57 auto* ctorPtr = dictionaryConstructorTable(
type);
66 *dictionaryConstructorTablePtr_
70 return autoPtr<IATEsource>(ctorPtr(iate,
dict));
86 *(otherPhase().
rho() - phase().
rho())
89 *
pow(
max(1 - phase(), scalar(0)), 1.75);
101 return max(Ur()*phase().d()/otherPhase().
nu(), scalar(1.0
e-3));
115 (16/
Re)*(1 + 0.15*
pow(
Re, 0.687)),
130 *(otherPhase().rho() - phase().rho())
142 *(otherPhase().rho() - phase().rho())
149 return otherPhase().rho()*
sqr(Ur())*phase().d()/
fluid().sigma();
errorManipArg< error, int > exit(error &err, const int errNo=1)
dimensioned< typename typeOfMag< Type >::type > mag(const dimensioned< Type > &dt)
UniformDimensionedField< vector > uniformDimensionedVectorField
label max(const labelHashSet &set, label maxValue=labelMin)
Find the max value in labelHashSet, optionally limited by second argument.
dimensionedSymmTensor sqr(const dimensionedVector &dv)
static autoPtr< IATEsource > New(const word &type, const IATE &iate, const dictionary &dict)
dimensionedScalar sqrt(const dimensionedScalar &ds)
dimensionedScalar pow025(const dimensionedScalar &ds)
tmp< volScalarField > Ur() const
Return the bubble relative velocity.
defineTypeNameAndDebug(constant, 0)
label k
Boltzmann constant.
tmp< volScalarField > We() const
Return the bubble Webber number.
dimensionedScalar sigma("sigma", dimMass/sqr(dimTime), transportProperties)
GeometricField< scalar, fvPatchField, volMesh > volScalarField
fileName::Type type(const fileName &name, const bool followLink=true)
Return the file type: DIRECTORY or FILE, normally following symbolic links.
const dimensionedScalar e
Elementary charge.
tmp< volScalarField > CD() const
Return the bubble drag coefficient.
label min(const labelHashSet &set, label minValue=labelMax)
Find the min value in labelHashSet, optionally limited by second argument.
tmp< volScalarField > Ut() const
Return the bubble turbulent velocity.
const uniformDimensionedVectorField & g
Info<< "Reading field U\"<< endl;volVectorField U(IOobject("U", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE), mesh);volScalarField rho(IOobject("rho", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::AUTO_WRITE), thermo.rho());volVectorField rhoU(IOobject("rhoU", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::NO_WRITE), rho *U);volScalarField rhoE(IOobject("rhoE", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::NO_WRITE), rho *(e+0.5 *magSqr(U)));surfaceScalarField pos(IOobject("pos", runTime.timeName(), mesh), mesh, dimensionedScalar("pos", dimless, 1.0));surfaceScalarField neg(IOobject("neg", runTime.timeName(), mesh), mesh, dimensionedScalar("neg", dimless, -1.0));surfaceScalarField phi("phi", fvc::flux(rhoU));Info<< "Creating turbulence model\"<< endl;autoPtr< compressible::turbulenceModel > turbulence(compressible::turbulenceModel::New(rho, U, phi, thermo))
static const gravity & New(const word &name, const Time &runTime)
Return named gravity field cached or construct on Time.
scalarField Re(const UList< complex > &cmplx)
Extract real component.
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
dimensionedScalar pow3(const dimensionedScalar &ds)
dimensionedScalar pow4(const dimensionedScalar &ds)
tmp< volScalarField > Eo() const
Return the bubble Eotvos number.
Pointer management similar to std::unique_ptr, with some additional methods and type checking...
tmp< volScalarField > Mo() const
Return the bubble Morton number.
A class for managing temporary objects.
#define FatalIOErrorInLookup(ios, lookupTag, lookupName, lookupTable)
Report an error message using Foam::FatalIOError.
defineRunTimeSelectionTable(binaryBreakupModel, dictionary)
tmp< volScalarField > Re() const
Return the bubble Reynolds number.
IOerror FatalIOError
Error stream (stdout output on all processes), with additional 'FOAM FATAL IO ERROR' header text and ...