incompressibleDenseParticleFluid Class Reference

Solver module for transient flow of incompressible isothermal fluids coupled with particle clouds including the effect of the volume fraction of particles on the continuous phase, with optional mesh motion and change. More...

Inheritance diagram for incompressibleDenseParticleFluid:
Collaboration diagram for incompressibleDenseParticleFluid:

Public Member Functions

 TypeName ("incompressibleDenseParticleFluid")
 Runtime type information. More...
 
 incompressibleDenseParticleFluid (fvMesh &mesh)
 Construct from region mesh. More...
 
 incompressibleDenseParticleFluid (const incompressibleDenseParticleFluid &)=delete
 Disallow default bitwise copy construction. More...
 
virtual ~incompressibleDenseParticleFluid ()
 Destructor. More...
 
virtual void preSolve ()
 Called at the start of the time-step, before the PIMPLE loop. More...
 
virtual void moveMesh ()
 Called at the start of the PIMPLE loop to move the mesh. More...
 
virtual void prePredictor ()
 Called at the start of the PIMPLE loop. More...
 
virtual void momentumPredictor ()
 Construct and optionally solve the momentum equation. More...
 
virtual void thermophysicalPredictor ()
 Construct and solve the energy equation,. More...
 
virtual void pressureCorrector ()
 Construct and solve the pressure equation in the PISO loop. More...
 
virtual void postCorrector ()
 Correct the momentum and thermophysical transport modelling. More...
 
virtual void postSolve ()
 Called after the PIMPLE loop at the end of the time-step. More...
 
void operator= (const incompressibleDenseParticleFluid &)=delete
 Disallow default bitwise assignment. More...
 
- Public Member Functions inherited from fluidSolver
 TypeName ("fluidSolver")
 Runtime type information. More...
 
 fluidSolver (fvMesh &mesh)
 Construct from region mesh. More...
 
 fluidSolver (const fluidSolver &)=delete
 Disallow default bitwise copy construction. More...
 
virtual ~fluidSolver ()
 Destructor. More...
 
virtual scalar maxDeltaT () const
 Return the current maximum time-step for stable solution. More...
 
void operator= (const fluidSolver &)=delete
 Disallow default bitwise assignment. More...
 
- Public Member Functions inherited from solver
 TypeName ("solver")
 Runtime type information. More...
 
 declareRunTimeSelectionTable (autoPtr, solver, fvMesh,(fvMesh &mesh),(mesh))
 
 solver (fvMesh &mesh)
 Construct from region mesh. More...
 
 solver (const solver &)=delete
 Disallow default bitwise copy construction. More...
 
virtual ~solver ()
 Destructor. More...
 
bool transient () const
 
Foam::fvModelsfvModels () const
 Return the fvModels that are created on demand. More...
 
Foam::fvConstraintsfvConstraints () const
 Return the fvConstraints that are created on demand. More...
 
void operator= (const solver &)=delete
 Disallow default bitwise assignment. More...
 
- Public Member Functions inherited from regIOobject
 TypeName ("regIOobject")
 Runtime type information. More...
 
 regIOobject (const IOobject &, const bool isTime=false)
 Construct from IOobject. Optional flag for if IOobject is the. More...
 
 regIOobject (const regIOobject &)
 Copy constructor. More...
 
 regIOobject (const regIOobject &&)
 Move constructor. More...
 
 regIOobject (const regIOobject &, bool registerCopy)
 Copy constructor, transferring registry registration to copy. More...
 
 regIOobject (const word &newName, const regIOobject &, bool registerCopy)
 Copy constructor with new name, transferring registry registration. More...
 
 regIOobject (const IOobject &, const regIOobject &)
 Copy constructor with new IO parameters. More...
 
virtual ~regIOobject ()
 Destructor. More...
 
virtual bool global () const
 Return true if object is global, i.e. same for all processors. More...
 
virtual bool globalFile () const
 Return true if object is global, i.e. same for all processors. More...
 
const fileNamecaseName () const
 
fileName path () const
 Return complete path. More...
 
fileName objectPath () const
 Return complete path + object name. More...
 
fileName filePath () const
 Return complete path + object name if the file exists. More...
 
bool checkIn ()
 Add object to registry. More...
 
bool checkOut ()
 Remove object from registry. More...
 
void addWatch ()
 Add file watch on object (if registered and READ_IF_MODIFIED) More...
 
bool registered () const
 Is this object registered with the registry? More...
 
bool ownedByRegistry () const
 Is this object owned by the registry? More...
 
void store ()
 Transfer ownership of this object to its registry. More...
 
void release ()
 Release ownership of this object from its registry. More...
 
label eventNo () const
 Event number at last update. More...
 
labeleventNo ()
 Event number at last update. More...
 
bool upToDate (const regIOobject &) const
 Return true if up-to-date with respect to given object. More...
 
bool upToDate (const regIOobject &, const regIOobject &) const
 Return true if up-to-date with respect to given objects. More...
 
bool upToDate (const regIOobject &, const regIOobject &, const regIOobject &) const
 Return true if up-to-date with respect to given objects. More...
 
bool upToDate (const regIOobject &, const regIOobject &, const regIOobject &, const regIOobject &) const
 Return true if up-to-date with respect to given objects. More...
 
void setUpToDate ()
 Set up to date (obviously) More...
 
virtual void rename (const word &newName)
 Rename. More...
 
bool headerOk ()
 Read and check header info. More...
 
IstreamreadStream (const word &, const bool read=true)
 Return Istream and check object type against that given. More...
 
void close ()
 Close Istream. More...
 
virtual bool readData (Istream &)
 Virtual readData function. More...
 
virtual bool read ()
 Read object. More...
 
const labelListwatchIndices () const
 Return file-monitoring handles. More...
 
labelListwatchIndices ()
 Return file-monitoring handles. More...
 
virtual bool modified () const
 Return true if the object's file (or files for objectRegistry) More...
 
virtual bool readIfModified ()
 Read object if modified (as set by call to modified) More...
 
virtual bool writeObject (IOstream::streamFormat, IOstream::versionNumber, IOstream::compressionType, const bool write) const
 Write using given format, version and compression. More...
 
virtual bool write (const bool write=true) const
 Write using setting from DB. More...
 
void operator= (const IOobject &)
 
const fileNamecaseName (const bool global) const
 
fileName path (const bool global) const
 Return complete path including the processor sub-directory. More...
 
fileName objectPath (const bool global) const
 Return complete path + object name including the processor. More...
 
fileName filePath (const word &typeName, const bool global) const
 Return complete path + object name if the file exists. More...
 
- Public Member Functions inherited from IOobject
 TypeName ("IOobject")
 Runtime type information. More...
 
 IOobject (const word &name, const fileName &instance, const objectRegistry &registry, readOption r=NO_READ, writeOption w=NO_WRITE, bool registerObject=true)
 Construct from name, instance, registry, io options. More...
 
 IOobject (const word &name, const fileName &instance, const fileName &local, const objectRegistry &registry, readOption r=NO_READ, writeOption w=NO_WRITE, bool registerObject=true)
 Construct from name, instance, local, registry, io options. More...
 
 IOobject (const fileName &path, const objectRegistry &registry, readOption r=NO_READ, writeOption w=NO_WRITE, bool registerObject=true)
 Construct from path, registry, io options. More...
 
 IOobject (const IOobject &io, const objectRegistry &registry)
 Construct from copy resetting registry. More...
 
 IOobject (const IOobject &io, const word &name)
 Construct from copy resetting name. More...
 
 IOobject (const IOobject &io)=default
 Copy constructor. More...
 
autoPtr< IOobjectclone () const
 Clone. More...
 
autoPtr< IOobjectclone (const objectRegistry &registry) const
 Clone resetting registry. More...
 
virtual ~IOobject ()
 Destructor. More...
 
const Timetime () const
 Return time. More...
 
const objectRegistrydb () const
 Return the local objectRegistry. More...
 
const wordname () const
 Return name. More...
 
const wordheaderClassName () const
 Return name of the class name read from header. More...
 
wordheaderClassName ()
 Return name of the class name read from header. More...
 
stringnote ()
 Return non-constant access to the optional note. More...
 
const stringnote () const
 Return the optional note. More...
 
bool & registerObject ()
 Register object created from this IOobject with registry if true. More...
 
bool registerObject () const
 Register object created from this IOobject with registry if true. More...
 
readOption readOpt () const
 
readOptionreadOpt ()
 
writeOption writeOpt () const
 
writeOptionwriteOpt ()
 
word group () const
 Return group (extension part of name) More...
 
word member () const
 Return member (name without the extension) More...
 
const fileNamerootPath () const
 
const fileNamecaseName (const bool global) const
 
fileNameinstance () const
 Return the instance directory, constant, system, <time> etc. More...
 
void updateInstance () const
 If the instance is a time directory update to the current time. More...
 
const fileNamelocal () const
 
fileName path (const bool global) const
 Return complete path including the processor sub-directory. More...
 
fileName objectPath (const bool global) const
 Return complete path + object name including the processor. More...
 
fileName relativePath () const
 Return the path relative to the case directory. More...
 
fileName relativeObjectPath () const
 Return complete relativePath + object name. More...
 
fileName filePath (const word &typeName, const bool global) const
 Return complete path + object name if the file exists. More...
 
bool readHeader (Istream &)
 Read header. More...
 
bool headerOk ()
 Read header of local object without type-checking. More...
 
template<class Type >
void warnNoRereading () const
 Helper: warn that type does not support re-reading. More...
 
bool writeHeader (Ostream &) const
 Write header. More...
 
bool writeHeader (Ostream &, const word &objectType) const
 Write header. Allow override of type. More...
 
bool good () const
 
bool bad () const
 
InfoProxy< IOobjectinfo () const
 Return info proxy. More...
 
void operator= (const IOobject &)
 
template<class Name >
Foam::word groupName (Name name, const word &group)
 

Public Attributes

const volScalarFieldp
 Reference to the pressure field. More...
 
const volVectorFieldUc
 Reference to the continuous phase velocity field. More...
 
const surfaceScalarFieldphic
 Reference to the continuous phase volumetric-flux field. More...
 
const volScalarFieldalphac
 Reference continuous phase-fraction. More...
 
- Public Attributes inherited from solver
const fvMeshmesh
 Region mesh. More...
 
const TimerunTime
 Time. More...
 
pimpleNoLoopControl pimple
 PIMPLE inner-loop controls. More...
 

Protected Member Functions

void correctCoNum ()
 Correct the cached Courant numbers. More...
 
void continuityErrors ()
 Calculate and print the continuity errors. More...
 
virtual void correctPressure ()
 Construct the pressure equation. More...
 
- Protected Member Functions inherited from fluidSolver
void readControls ()
 Read controls. More...
 
void meshCourantNo () const
 Check mesh Courant numbers for moving mesh cases. More...
 
void correctCoNum (const surfaceScalarField &phi)
 Correct the cached Courant numbers. More...
 
void correctCoNum (const volScalarField &rho, const surfaceScalarField &phi)
 Correct the cached Courant numbers. More...
 
void continuityErrors (const surfaceScalarField &phi)
 Calculate and print the continuity errors. More...
 
void continuityErrors (const volScalarField &rho, const volScalarField &thermoRho, const surfaceScalarField &phi)
 Calculate and print the continuity errors. More...
 
- Protected Member Functions inherited from regIOobject
bool readHeaderOk (const IOstream::streamFormat defaultFormat, const word &typeName)
 Read header, check readOpt flags and read data if necessary. More...
 
- Protected Member Functions inherited from IOobject
void setBad (const string &)
 Set the object state to bad. More...
 
template<class Type >
bool typeHeaderOk (const bool checkType)
 Read header using typeGlobalFile to find file. More...
 

Protected Attributes

const word continuousPhaseName
 Name of the continuous phase. More...
 
volScalarField p_
 Pressure field. More...
 
Foam::pressureReference pressureReference
 Pressure reference. More...
 
uniformDimensionedVectorField g
 Acceleration due to gravity. More...
 
volVectorField Uc_
 Continuous phase velocity field. More...
 
surfaceScalarField phic_
 Continuous phase flux field. More...
 
autoPtr< viscosityModelviscosity
 Continuous phase kinematic viscosity model. More...
 
volScalarField rhoc
 Continuous phase density. More...
 
volScalarField muc
 Continuous phase viscosity. More...
 
volScalarField alphac_
 Continuous phase-fraction. More...
 
scalar alphacMin
 Minimum continuous phase-fraction. More...
 
surfaceScalarField alphacf
 Interpolated continuous phase-fraction. More...
 
surfaceScalarField alphaPhic
 Continuous phase volumetric-flux field. More...
 
autoPtr< phaseIncompressible::momentumTransportModelmomentumTransport
 Pointer to the momentum transport model. More...
 
parcelCloudsclouds
 
autoPtr< volVectorFieldFd
 Dispersed phase drag force. More...
 
autoPtr< volScalarFieldDc
 
autoPtr< surfaceScalarFieldDcf
 Face continuous-dispersed phase drag coefficient. More...
 
autoPtr< surfaceScalarFieldphid
 Effective volumetric flux of the dispersed phase. More...
 
autoPtr< surfaceVectorFieldUcf
 Pointer to the surface momentum field. More...
 
tmp< fvVectorMatrixtUcEqn
 Cached momentum matrix. More...
 
- Protected Attributes inherited from fluidSolver
bool correctPhi
 Switch to correct the flux after mesh change. More...
 
scalar CoNum
 Current maximum Courant number for time-step control. More...
 
- Protected Attributes inherited from solver
fvMeshmesh_
 Region mesh. More...
 
bool steady
 Switch for steady-state. More...
 
bool LTS
 Switch for local time step transient operation. More...
 

Additional Inherited Members

- Public Types inherited from IOobject
enum  objectState { GOOD , BAD }
 Enumeration defining the valid states of an IOobject. More...
 
enum  readOption { MUST_READ , MUST_READ_IF_MODIFIED , READ_IF_PRESENT , NO_READ }
 Enumeration defining the read options. More...
 
enum  writeOption { AUTO_WRITE = 0 , NO_WRITE = 1 }
 Enumeration defining the write options. More...
 
enum  fileCheckTypes { timeStamp , timeStampMaster , inotify , inotifyMaster }
 Enumeration defining the file checking options. More...
 
- Static Public Member Functions inherited from solver
static autoPtr< solverNew (const word &solverName, fvMesh &mesh)
 Select, construct and return the solver. More...
 
static void load (const word &solverName)
 Load the specified solver library. More...
 
- Static Public Member Functions inherited from regIOobject
template<class Type >
static Type & store (Type *)
 Transfer ownership of the given object pointer to its registry. More...
 
template<class Type >
static Type & store (autoPtr< Type > &)
 Transfer ownership of the given object pointer to its registry. More...
 
- Static Public Member Functions inherited from IOobject
static bool fileNameComponents (const fileName &path, fileName &instance, fileName &local, word &name)
 Split path into instance, local, name components. More...
 
template<class Name >
static word groupName (Name name, const word &group)
 
static word group (const word &name)
 Return group (extension part of name) More...
 
static word member (const word &name)
 Return member (name without the extension) More...
 
template<class Stream >
static Stream & writeBanner (Stream &os, bool noHint=false)
 Write the standard OpenFOAM file/dictionary banner. More...
 
template<class Stream >
static Stream & writeDivider (Stream &os)
 Write the standard file section divider. More...
 
template<class Stream >
static Stream & writeEndDivider (Stream &os)
 Write the standard end file divider. More...
 
- Static Public Attributes inherited from solver
static scalar deltaTFactor = 1.2
 deltaT increase factor More...
 
- Static Public Attributes inherited from regIOobject
static float fileModificationSkew
 
- Static Public Attributes inherited from IOobject
static constexpr const char * foamFile = "FoamFile"
 Keyword for the FoamFile header sub-dictionary. More...
 
static const NamedEnum< fileCheckTypes, 4 > fileCheckTypesNames
 
static fileCheckTypes fileModificationChecking
 Type of file modification checking. More...
 

Detailed Description

Solver module for transient flow of incompressible isothermal fluids coupled with particle clouds including the effect of the volume fraction of particles on the continuous phase, with optional mesh motion and change.

Uses the flexible PIMPLE (PISO-SIMPLE) solution for time-resolved and pseudo-transient and steady simulations.

Optional fvModels and fvConstraints are provided to enhance the simulation in many ways including adding various sources, constraining or limiting the solution.

Reference:

    Greenshields, C. J., & Weller, H. G. (2022).
    Notes on Computational Fluid Dynamics: General Principles.
    CFD Direct Ltd.: Reading, UK.
Source files
See also
Foam::solvers::fluidSolver Foam::solvers::incompressibleFluid

Definition at line 76 of file incompressibleDenseParticleFluid.H.

Constructor & Destructor Documentation

◆ incompressibleDenseParticleFluid() [1/2]

◆ incompressibleDenseParticleFluid() [2/2]

Disallow default bitwise copy construction.

◆ ~incompressibleDenseParticleFluid()

Destructor.

Definition at line 239 of file incompressibleDenseParticleFluid.C.

Member Function Documentation

◆ correctCoNum()

void correctCoNum ( )
protected

Correct the cached Courant numbers.

Definition at line 52 of file incompressibleDenseParticleFluid.C.

References incompressibleDenseParticleFluid::phic.

Referenced by incompressibleDenseParticleFluid::incompressibleDenseParticleFluid().

Here is the caller graph for this function:

◆ continuityErrors()

void continuityErrors ( )
protected

Calculate and print the continuity errors.

Definition at line 58 of file incompressibleDenseParticleFluid.C.

References fluidSolver::continuityErrors().

Here is the call graph for this function:

◆ correctPressure()

void correctPressure ( )
protectedvirtual

Construct the pressure equation.

and correct the pressure and velocity

Definition at line 39 of file correctPressure.C.

References fvMatrix< Type >::A(), Foam::adjustPhi(), incompressibleDenseParticleFluid::alphac, incompressibleDenseParticleFluid::alphacf, fvConstraints::constrain(), Foam::constrainHbyA(), Foam::constrainPressure(), GeometricField< Type, PatchField, GeoMesh >::correctBoundaryConditions(), pimpleNoLoopControl::correctNonOrthogonal(), Foam::fvc::correctUf(), incompressibleDenseParticleFluid::Dc, incompressibleDenseParticleFluid::Dcf, Foam::fvc::ddt(), Foam::fvc::ddtCorr(), Foam::fvc::div(), incompressibleDenseParticleFluid::Fd, nonOrthogonalSolutionControl::finalNonOrthogonalIter(), fvMatrix< Type >::flux(), Foam::fvc::flux(), solver::fvConstraints(), incompressibleDenseParticleFluid::g, fvMatrix< Type >::H(), HbyA, Foam::fvc::interpolate(), Foam::fvm::laplacian(), Foam::fvc::makeAbsolute(), Foam::fvc::makeRelative(), solver::mesh, GeometricField< Type, PatchField, GeoMesh >::needReference(), incompressibleDenseParticleFluid::p, incompressibleDenseParticleFluid::p_, incompressibleDenseParticleFluid::phic, incompressibleDenseParticleFluid::phic_, incompressibleDenseParticleFluid::phid, solver::pimple, Foam::fvc::reconstruct(), pressureReference::refCell(), pressureReference::refValue(), GeometricField< Type, PatchField, GeoMesh >::relax(), fvMatrix< Type >::setReference(), fvMesh::Sf(), fvMatrix< Type >::solve(), incompressibleDenseParticleFluid::tUcEqn, incompressibleDenseParticleFluid::Uc, incompressibleDenseParticleFluid::Uc_, and incompressibleDenseParticleFluid::Ucf.

Here is the call graph for this function:

◆ TypeName()

TypeName ( "incompressibleDenseParticleFluid"  )

Runtime type information.

◆ preSolve()

void preSolve ( )
virtual

Called at the start of the time-step, before the PIMPLE loop.

Implements solver.

Definition at line 246 of file incompressibleDenseParticleFluid.C.

References IOobject::AUTO_WRITE, Foam::endl(), fvModels(), Foam::Info, Foam::fvc::interpolate(), fvModels::preUpdateMesh(), and IOobject::READ_IF_PRESENT.

Here is the call graph for this function:

◆ moveMesh()

◆ prePredictor()

void prePredictor ( )
virtual

Called at the start of the PIMPLE loop.

Implements solver.

Definition at line 287 of file incompressibleDenseParticleFluid.C.

References lduMatrix::diag(), Foam::dimAcceleration, Foam::dimless, Foam::dimTime, Foam::fvc::flux(), Foam::fvc::interpolate(), Foam::max(), viscosity::nu(), pimple(), fvMatrix< Type >::source(), and Foam::Zero.

Here is the call graph for this function:

◆ momentumPredictor()

◆ thermophysicalPredictor()

void thermophysicalPredictor ( )
virtual

Construct and solve the energy equation,.

convert to temperature and update thermophysical and transport properties

Implements solver.

Definition at line 354 of file incompressibleDenseParticleFluid.C.

◆ pressureCorrector()

void pressureCorrector ( )
virtual

Construct and solve the pressure equation in the PISO loop.

Implements solver.

Definition at line 358 of file incompressibleDenseParticleFluid.C.

References pimple().

Here is the call graph for this function:

◆ postCorrector()

void postCorrector ( )
virtual

Correct the momentum and thermophysical transport modelling.

Implements solver.

Definition at line 369 of file incompressibleDenseParticleFluid.C.

References pimple().

Here is the call graph for this function:

◆ postSolve()

void postSolve ( )
virtual

Called after the PIMPLE loop at the end of the time-step.

Implements solver.

Definition at line 379 of file incompressibleDenseParticleFluid.C.

◆ operator=()

void operator= ( const incompressibleDenseParticleFluid )
delete

Disallow default bitwise assignment.

Member Data Documentation

◆ continuousPhaseName

const word continuousPhaseName
protected

Name of the continuous phase.

Definition at line 84 of file incompressibleDenseParticleFluid.H.

◆ p_

volScalarField p_
protected

Pressure field.

Definition at line 90 of file incompressibleDenseParticleFluid.H.

Referenced by incompressibleDenseParticleFluid::correctPressure().

◆ pressureReference

Pressure reference.

Definition at line 93 of file incompressibleDenseParticleFluid.H.

◆ g

◆ Uc_

◆ phic_

◆ viscosity

Continuous phase kinematic viscosity model.

Definition at line 111 of file incompressibleDenseParticleFluid.H.

◆ rhoc

volScalarField rhoc
protected

Continuous phase density.

Definition at line 117 of file incompressibleDenseParticleFluid.H.

◆ muc

volScalarField muc
protected

Continuous phase viscosity.

Definition at line 120 of file incompressibleDenseParticleFluid.H.

◆ alphac_

volScalarField alphac_
protected

Continuous phase-fraction.

Definition at line 123 of file incompressibleDenseParticleFluid.H.

Referenced by incompressibleDenseParticleFluid::incompressibleDenseParticleFluid().

◆ alphacMin

scalar alphacMin
protected

Minimum continuous phase-fraction.

Definition at line 126 of file incompressibleDenseParticleFluid.H.

Referenced by incompressibleDenseParticleFluid::incompressibleDenseParticleFluid().

◆ alphacf

◆ alphaPhic

◆ momentumTransport

◆ clouds

◆ Fd

◆ Dc

◆ Dcf

autoPtr<surfaceScalarField> Dcf
protected

Face continuous-dispersed phase drag coefficient.

Definition at line 155 of file incompressibleDenseParticleFluid.H.

Referenced by incompressibleDenseParticleFluid::correctPressure(), and incompressibleDenseParticleFluid::momentumPredictor().

◆ phid

autoPtr<surfaceScalarField> phid
protected

Effective volumetric flux of the dispersed phase.

Definition at line 158 of file incompressibleDenseParticleFluid.H.

Referenced by incompressibleDenseParticleFluid::correctPressure(), and incompressibleDenseParticleFluid::momentumPredictor().

◆ Ucf

autoPtr<surfaceVectorField> Ucf
protected

Pointer to the surface momentum field.

used to recreate the flux after mesh-change

Definition at line 162 of file incompressibleDenseParticleFluid.H.

Referenced by incompressibleDenseParticleFluid::correctPressure(), and incompressibleDenseParticleFluid::moveMesh().

◆ tUcEqn

tmp<fvVectorMatrix> tUcEqn
protected

Cached momentum matrix.

shared between the momentum predictor and pressure corrector

Definition at line 166 of file incompressibleDenseParticleFluid.H.

Referenced by incompressibleDenseParticleFluid::correctPressure(), and incompressibleDenseParticleFluid::momentumPredictor().

◆ p

◆ Uc

◆ phic

◆ alphac


The documentation for this class was generated from the following files: