66 int main(
int argc,
char *argv[])
70 "Passive scalar transport finiteArea equation solver." 77 #include "createFaFields.H" 78 #include "createVolFields.H" 80 Info<<
"Total mass of surfactant: " 81 <<
sum(
Cs.internalField()*aMesh.S()) <<
endl;
83 Info<<
"\nStarting time loop\n" <<
endl;
100 vsm.mapToVolume(
Cs,
Cvf.boundaryFieldRef());
105 Info<<
"Total mass of surfactant: " 106 <<
sum(
Cs.internalField()*aMesh.S()) <<
endl;
dimensioned< Type > sum(const DimensionedField< Type, GeoMesh > &f1)
const volSurfaceMapping vsm(aMesh)
tmp< GeometricField< Type, faPatchField, areaMesh > > div(const GeometricField< Type, faePatchField, edgeMesh > &ssf)
Ostream & endl(Ostream &os)
Add newline and flush stream.
tmp< GeometricField< Type, faPatchField, areaMesh > > ddt(const dimensioned< Type > dt, const faMesh &mesh)
volScalarField Cvf(IOobject("Cvf", runTime.timeName(), mesh.thisDb(), IOobject::NO_READ, IOobject::AUTO_WRITE), mesh, dimensionedScalar(dimless/dimLength, Zero))
tmp< GeometricField< Type, faPatchField, areaMesh > > laplacian(const GeometricField< Type, faPatchField, areaMesh > &vf, const word &name)
edgeScalarField phis(IOobject("phis", runTime.timeName(), aMesh.thisDb(), IOobject::NO_READ, IOobject::NO_WRITE), linearEdgeInterpolate(Us) &aMesh.Le())
faMatrix< scalar > faScalarMatrix
messageStream Info
Information stream (stdout output on master, null elsewhere)