37 namespace regionModels
39 namespace areaSurfaceFilmModels
51 const word& modelType,
114 for (
int corr=1; corr<=
nCorr_; corr++)
SolverPerformance< Type > solve(faMatrix< Type > &, const dictionary &)
Solve returning the solution statistics given convergence tolerance.
uniformDimensionedVectorField g_
Gravity.
void constrain(faMatrix< Type > &eqn)
Apply constraints to equation.
A list of keyword definitions, which are a keyword followed by a number of values (eg...
Foam::fa::options & faOptions()
Return faOptions.
areaScalarField gn_
Normal gravity field.
autoPtr< filmTurbulenceModel > turbulence_
Turbulence model.
void relax(const scalar alpha)
Relax matrix (for steady-state solution).
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 tmp< GeometricField< Type, faePatchField, edgeMesh > > interpolate(const GeometricField< Type, faPatchField, areaMesh > &tvf, const edgeScalarField &faceFlux, Istream &schemeData)
Interpolate field onto faces using scheme given by Istream.
edgeScalarField phif_
Film momentum flux.
constexpr char nl
The newline '\n' character (0x0a)
dimensionedScalar h0_
Smallest numerical thickness.
areaVectorField USp_
Momentum source.
areaScalarField h_
Film hight.
Ostream & endl(Ostream &os)
Add newline and flush stream.
areaVectorField Uf_
Film velocity.
edgeScalarField phi2s_
Film height flux.
tmp< areaScalarField > pg() const
Map primary static pressure.
defineTypeNameAndDebug(kinematicThinFilm, 0)
label nFilmCorr_
Number of film thickness correctors.
areaScalarField rho_
Density [kg/m3].
const edgeScalarField & magLe() const
Return edge length magnitudes.
Macros for easy insertion into run-time selection tables.
void correct(GeometricField< Type, faPatchField, areaMesh > &field)
Apply correction to field.
areaScalarField sigma_
Surface tension [m/s2].
label nOuterCorr_
Number of outer correctors.
addToRunTimeSelectionTable(liquidFilmBase, kinematicThinFilm, dictionary)
A class for handling words, derived from Foam::string.
virtual void evolveRegion()
Evolve the film.
const edgeVectorField & Le() const
Return edge length vectors.
tmp< GeometricField< Type, faPatchField, areaMesh > > laplacian(const GeometricField< Type, faPatchField, areaMesh > &vf, const word &name)
tmp< faVectorMatrix > correct(areaVectorField &U)
Return (net) force system.
virtual void postEvolveRegion()
Post-evolve film.
kinematicThinFilm(const word &modelType, const fvMesh &mesh, const dictionary &dict)
Construct from components and dict.
const faMesh & regionMesh() const
Return the region mesh database.
MinMax< Type > gMinMax(const FieldField< Field, Type > &f)
forceList forces_
Transfer with the continuous phase.
tmp< faMatrix< Type > > ddt(const GeometricField< Type, faPatchField, areaMesh > &vf)
tmp< faMatrix< Type > > div(const edgeScalarField &flux, const GeometricField< Type, faPatchField, areaMesh > &vf, const word &name)
int debug
Static debugging option.
scalarMinMax gMinMaxMag(const FieldField< Field, Type > &f)
areaScalarField pf_
Film pressure.
tmp< GeometricField< Type, faePatchField, edgeMesh > > lnGrad(const GeometricField< Type, faPatchField, areaMesh > &vf, const word &name)
virtual void preEvolveRegion()
Pre-evolve film.
areaScalarField ppf_
Primary region pressure.
virtual void postEvolveRegion()
Post-evolve film.
areaScalarField pnSp_
Normal pressure by particles.
void relax(const scalar alpha)
Relax field (for steady-state solution).
Mesh data needed to do the Finite Volume discretisation.
virtual void preEvolveRegion()
Pre-evolve film.
A special matrix type and solver, designed for finite area solutions of scalar equations. Face addressing is used to make all matrix assembly and solution loops vectorise.
Template specialisation for scalar faMatrix.
void correctThermoFields()
Correct thermo.
void correctBoundaryConditions()
Correct boundary field.
messageStream Info
Information stream (stdout output on master, null elsewhere)
void storePrevIter() const
Store the field as the previous iteration value.
label nCorr_
Number of PISO-like correctors.
const areaVectorField & faceAreaNormals() const
Return face area normals.
tmp< GeometricField< typename outerProduct< vector, Type >::type, faPatchField, areaMesh >> grad(const GeometricField< Type, faePatchField, edgeMesh > &ssf)
areaScalarField rhoSp_
Mass source.
#define InfoInFunction
Report an information message using Foam::Info.
const dimensionSet dimVelocity
Switch momentumPredictor_
Momentum predictor.