This source is designed to model the effect of solidification and melting processes, e.g. windhield defrosting, within a specified region. The phase change occurs at the melting temperature, Tmelt
.
More...
Public Types | |
enum | thermoMode { mdThermo, mdLookup } |
Options for the thermo mode specification. More... | |
Public Types inherited from cellSetOption | |
enum | selectionModeType { smAll, smCellSet, smCellZone, smPoints, smGeometric, smCellType } |
Enumeration for selection mode types. More... | |
Public Member Functions | |
TypeName ("solidificationMeltingSource") | |
Runtime type information. More... | |
solidificationMeltingSource (const word &sourceName, const word &modelType, const dictionary &dict, const fvMesh &mesh) | |
Construct from explicit source name and mesh. More... | |
solidificationMeltingSource (const solidificationMeltingSource &)=delete | |
No copy construct. More... | |
void | operator= (const solidificationMeltingSource &)=delete |
No copy assignment. More... | |
~solidificationMeltingSource ()=default | |
Destructor. More... | |
virtual void | addSup (fvMatrix< scalar > &eqn, const label fieldi) |
Add explicit contribution to enthalpy equation. More... | |
virtual void | addSup (fvMatrix< vector > &eqn, const label fieldi) |
Add implicit contribution to momentum equation. More... | |
virtual void | addSup (const volScalarField &rho, fvMatrix< scalar > &eqn, const label fieldi) |
Add explicit contribution to compressible enthalpy equation. More... | |
virtual void | addSup (const volScalarField &rho, fvMatrix< vector > &eqn, const label fieldi) |
Add implicit contribution to compressible momentum equation. More... | |
virtual bool | read (const dictionary &dict) |
Read source dictionary. More... | |
Public Member Functions inherited from cellSetOption | |
TypeName ("cellSetOption") | |
Runtime type information. More... | |
cellSetOption (const word &name, const word &modelType, const dictionary &dict, const fvMesh &mesh) | |
Construct from components. More... | |
virtual | ~cellSetOption ()=default |
Destructor. More... | |
scalar | timeStart () const noexcept |
Return const access to the time start. More... | |
scalar | duration () const noexcept |
Return const access to the duration. More... | |
bool | inTimeLimits (const scalar timeValue) const |
True if within time limits. More... | |
bool | useSubMesh () const noexcept |
True if sub-selection should be used. More... | |
selectionModeType | selectionMode () const noexcept |
Return the cell selection mode. More... | |
const wordRes & | selectionNames () const noexcept |
Return const access to the selection names (set or zone selection) More... | |
const wordRe & | zoneName () const |
Return const access to the first set/zone name. More... | |
scalar | V () const noexcept |
Return const access to the total cell volume. More... | |
const labelList & | cells () const noexcept |
Return const access to the cell selection. More... | |
scalar | timeStart (scalar val) noexcept |
Adjust the time start, return the old value. More... | |
scalar | duration (scalar val) noexcept |
Adjust the duration, return the old value. More... | |
virtual bool | isActive () |
Is the source active? More... | |
const word & | cellSetName () const |
The name of the cell set/zone [as a word] for "cellSet" / "cellZone" selection modes) More... | |
Public Member Functions inherited from option | |
TypeName ("option") | |
Runtime type information. More... | |
declareRunTimeSelectionTable (autoPtr, option, dictionary,(const word &name, const word &modelType, const dictionary &dict, const fvMesh &mesh),(name, modelType, dict, mesh)) | |
option (const word &name, const word &modelType, const dictionary &dict, const fvMesh &mesh) | |
Construct from components. More... | |
autoPtr< option > | clone () const |
Return clone. More... | |
virtual | ~option ()=default |
Destructor. More... | |
const word & | name () const noexcept |
Return const access to the source name. More... | |
const fvMesh & | mesh () const noexcept |
Return const access to the mesh database. More... | |
const dictionary & | coeffs () const noexcept |
Return dictionary. More... | |
bool | active () const noexcept |
True if source is active. More... | |
void | setApplied (const label fieldi) |
Set the applied flag to true for field index fieldi. More... | |
bool | active (const bool on) noexcept |
Change source active flag, return previous value. More... | |
virtual label | applyToField (const word &fieldName) const |
Return index of field name if found in fieldNames list. More... | |
virtual void | checkApplied () const |
Check that the source has been applied. More... | |
virtual void | addSup (fvMatrix< symmTensor > &eqn, const label fieldi) |
virtual void | addSup (fvMatrix< sphericalTensor > &eqn, const label fieldi) |
virtual void | addSup (fvMatrix< tensor > &eqn, const label fieldi) |
virtual void | addSup (const volScalarField &rho, fvMatrix< symmTensor > &eqn, const label fieldi) |
virtual void | addSup (const volScalarField &rho, fvMatrix< sphericalTensor > &eqn, const label fieldi) |
virtual void | addSup (const volScalarField &rho, fvMatrix< tensor > &eqn, const label fieldi) |
virtual void | addSup (const volScalarField &alpha, const volScalarField &rho, fvMatrix< scalar > &eqn, const label fieldi) |
virtual void | addSup (const volScalarField &alpha, const volScalarField &rho, fvMatrix< vector > &eqn, const label fieldi) |
virtual void | addSup (const volScalarField &alpha, const volScalarField &rho, fvMatrix< symmTensor > &eqn, const label fieldi) |
virtual void | addSup (const volScalarField &alpha, const volScalarField &rho, fvMatrix< sphericalTensor > &eqn, const label fieldi) |
virtual void | addSup (const volScalarField &alpha, const volScalarField &rho, fvMatrix< tensor > &eqn, const label fieldi) |
virtual void | constrain (fvMatrix< scalar > &eqn, const label fieldi) |
virtual void | constrain (fvMatrix< vector > &eqn, const label fieldi) |
virtual void | constrain (fvMatrix< sphericalTensor > &eqn, const label fieldi) |
virtual void | constrain (fvMatrix< symmTensor > &eqn, const label fieldi) |
virtual void | constrain (fvMatrix< tensor > &eqn, const label fieldi) |
virtual void | correct (volScalarField &field) |
virtual void | correct (volVectorField &field) |
virtual void | correct (volSphericalTensorField &field) |
virtual void | correct (volSymmTensorField &field) |
virtual void | correct (volTensorField &field) |
virtual void | postProcessSens (scalarField &sensField, const word &fieldName=word::null, const word &designVariablesName=word::null) |
virtual void | postProcessSens (vectorField &sensField, const word &fieldName=word::null, const word &designVariablesName=word::null) |
virtual void | postProcessSens (tensorField &sensField, const word &fieldName=word::null, const word &designVariablesName=word::null) |
virtual void | writeHeader (Ostream &) const |
Write the source header information. More... | |
virtual void | writeFooter (Ostream &) const |
Write the source footer information. More... | |
virtual void | writeData (Ostream &) const |
Write the source properties. More... | |
Static Public Attributes | |
static const Enum< thermoMode > | thermoModeTypeNames_ |
Names for thermoMode. More... | |
Static Public Attributes inherited from cellSetOption | |
static const Enum< selectionModeType > | selectionModeTypeNames_ |
List of selection mode type names. More... | |
Additional Inherited Members | |
Static Public Member Functions inherited from option | |
static autoPtr< option > | New (const word &name, const dictionary &dict, const fvMesh &mesh) |
Return a reference to the selected fvOption model. More... | |
Public Attributes inherited from option | |
bool | log |
Switch write log to Info. More... | |
Protected Member Functions inherited from cellSetOption | |
void | setSelection (const dictionary &dict) |
Set cell selection name or points selection from dictionary input. More... | |
void | setCellSelection () |
Set the cell selection based on user input selection mode. More... | |
void | setVol () |
Recalculate the volume. More... | |
Protected Member Functions inherited from option | |
void | resetApplied () |
Resize/reset applied flag list for all fieldNames_ entries. More... | |
Protected Attributes inherited from cellSetOption | |
scalar | timeStart_ |
Start time of fvOption. More... | |
scalar | duration_ |
Duration of fvOption execution starting from timeStart. More... | |
selectionModeType | selectionMode_ |
Cell selection mode. More... | |
wordRes | selectionNames_ |
Face selection names (for set or zone selections) More... | |
List< point > | points_ |
List of points for "points" selectionMode. More... | |
dictionary | geometricSelection_ |
Dictionary entries for "geometric" (topoSetCellSource) selection. More... | |
labelList | cells_ |
Set of cells to apply source to. More... | |
scalar | V_ |
Sum of cell volumes. More... | |
Protected Attributes inherited from option | |
const word | name_ |
Source name. More... | |
const word | modelType_ |
Model type. More... | |
const fvMesh & | mesh_ |
Reference to the mesh database. More... | |
dictionary | dict_ |
Top level source dictionary. More... | |
dictionary | coeffs_ |
Dictionary containing source coefficients. More... | |
wordList | fieldNames_ |
Field names to apply source to - populated by derived models. More... | |
List< bool > | applied_ |
Applied flag list - corresponds to each fieldNames_ entry. More... | |
bool | active_ |
Source active flag. More... | |
This source is designed to model the effect of solidification and melting processes, e.g. windhield defrosting, within a specified region. The phase change occurs at the melting temperature, Tmelt
.
The presence of the solid phase in the flow field is incorporated into the model as a momentum porosity contribution; the energy associated with the phase change is added as an enthalpy contribution.
References:
Voller, V. R., & Prakash, C. (1987). A fixed grid numerical modelling methodology for convection-diffusion mushy region phase-change problems. International Journal of Heat and Mass Transfer, 30(8), 1709-1719. DOI:10.1016/0017-9310(87)90317-6 Swaminathan, C. R., & Voller, V. R. (1992). A general enthalpy method for modeling solidification processes. Metallurgical transactions B, 23(5), 651-664. DOI:10.1007/BF02649725
The model generates a field <name>
:alpha1 which can be visualised to to show the melt distribution as a fraction [0-1].
constant/fvOptions
: solidificationMeltingSource1 { // Mandatory entries (unmodifiable) type solidificationMeltingSource; // Mandatory entries (runtime modifiable) Tmelt 273; L 334000; thermoMode <thermoModeName>; rhoRef 800; beta 5e-6; // Optional entries (runtime modifiable) relax 0.9; T <Tname>; rho <rhoName>; U <Uname>; phi <phiName>; Cu 1e5; q 1e-2; // Conditional optional entries (runtime modifiable) // when thermoMode=lookup Cp Cp; // Conditional mandatory entries (runtime modifiable) // when Cp=CpRef CpRef 1000; // Mandatory/Optional (inherited) entries ... }
where the entries mean:
Property | Description | Type | Reqd | Dflt |
---|---|---|---|---|
type | Type name: solidificationMeltingSource | word | yes | - |
Tmelt | Melting temperature [K] | scalar | yes | - |
L | Latent heat of fusion [J/kg] | scalar | yes | - |
thermoMode | Thermo mode | word | yes | - |
rhoRef | Reference (solid) density | scalar | yes | - |
beta | Thermal expansion coefficient [1/K] | scalar | yes | - |
relax | Relaxation factor [0-1] | scakar | no | 0.9 |
T | Name of operand temperature field | word | no | T |
rho | Name of operand density field | word | no | rho |
U | Name of operand velocity field | word | no | U |
phi | Name of operand flux field | word | no | phi |
Cu | Mushy region momentum sink coefficient [1/s] | scalar | no | 1e5 |
q | Coefficient used in porosity calc | scalar | no | 1e-2 |
Cp | Name of specific heat capacity field | word | cndtnl | Cp |
CpRef | Specific heat capacity value | scalar | cndtnl | - |
The inherited entries are elaborated in:
Options for the thermoMode
entry:
thermo | Access Cp information from database lookup | Access Cp information by looking up from dictionary
Definition at line 247 of file solidificationMeltingSource.H.
enum thermoMode |
Options for the thermo mode specification.
Enumerator | |
---|---|
mdThermo | |
mdLookup |
Definition at line 258 of file solidificationMeltingSource.H.
solidificationMeltingSource | ( | const word & | sourceName, |
const word & | modelType, | ||
const dictionary & | dict, | ||
const fvMesh & | mesh | ||
) |
Construct from explicit source name and mesh.
Definition at line 156 of file solidificationMeltingSource.C.
References Foam::abort(), dictionary::dictName(), Foam::FatalError, FatalErrorInFunction, and option::resetApplied().
|
delete |
No copy construct.
|
default |
Destructor.
TypeName | ( | "solidificationMeltingSource" | ) |
Runtime type information.
|
delete |
No copy assignment.
|
virtual |
Add explicit contribution to enthalpy equation.
Reimplemented from option.
Definition at line 226 of file solidificationMeltingSource.C.
References Foam::apply().
Add implicit contribution to momentum equation.
Reimplemented from option.
Definition at line 247 of file solidificationMeltingSource.C.
References Cp, Foam::ensightOutput::debug, lduMatrix::diag(), Foam::endl(), forAll, g, Foam::Info, IOobject::name(), gravity::New(), Foam::pow3(), fvMatrix< Type >::psi(), fvMatrix< Type >::source(), Sp, Foam::sqr(), Su, Foam::type(), and update().
|
virtual |
Add explicit contribution to compressible enthalpy equation.
Reimplemented from option.
Definition at line 236 of file solidificationMeltingSource.C.
References Foam::apply(), and rho.
|
virtual |
Add implicit contribution to compressible momentum equation.
Reimplemented from option.
Definition at line 284 of file solidificationMeltingSource.C.
|
virtual |
Read source dictionary.
Reimplemented from cellSetOption.
Definition at line 295 of file solidificationMeltingSource.C.
References dict, and cellSetOption::read().
|
static |
Names for thermoMode.
Definition at line 267 of file solidificationMeltingSource.H.
Copyright © 2011-2018 OpenFOAM | OPENFOAM® is a registered trademark of OpenCFD Ltd.