Files
openfoam/tutorials/multiphase/reactingTwoPhaseEulerFoam/laminar/bubbleColumnEvaporating/system/continuityFunctions
Henry Weller c983670c91 functionObjects: Changed options 'outputControl' -> 'writeControl' and 'outputInterval' -> 'writeInterval'
for consistency with the time controls in controlDict and to avoid
unnecessary confusion.  All code and tutorials have been updated.

The old names 'outputControl' and 'outputInterval' are but supported for
backward compatibility but deprecated.
2016-05-12 11:38:11 +01:00

90 lines
2.5 KiB
C++

/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: dev |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
inletMassFlowRate
{
type faceSource;
functionObjectLibs ("libfieldFunctionObjects.so");
fields
(
alphaRhoPhi.gas
alphaRhoPhi.liquid
);
valueOutput false;
log true;
surfaceFormat null;
source patch;
sourceName inlet;
operation sum;
writeControl timeStep;
}
outletMassFlowRate
{
type faceSource;
functionObjectLibs ("libfieldFunctionObjects.so");
fields
(
alphaRhoPhi.gas
alphaRhoPhi.liquid
);
valueOutput false;
log true;
surfaceFormat null;
source patch;
sourceName outlet;
operation sum;
writeControl timeStep;
}
totalMass
{
type coded;
functionObjectLibs ("libutilityFunctionObjects.so");
redirectType error;
code
#{
const volScalarField& alphaGas =
mesh().lookupObject<volScalarField>("alpha.gas");
const volScalarField& alphaLiquid =
mesh().lookupObject<volScalarField>("alpha.liquid");
const volScalarField& rhoGas =
mesh().lookupObject<volScalarField>("thermo:rho.gas");
const volScalarField& rhoLiquid =
mesh().lookupObject<volScalarField>("thermo:rho.liquid");
const volScalarField& dmdt =
mesh().lookupObject<volScalarField>("dmdt.gasAndLiquid");
const scalarField& v = mesh().V();
Info<< "coded totalMass output:" << endl
<< " volIntegrate(all) for alpha.gas*rho.gas = "
<< gSum(alphaGas*rhoGas*v) << endl
<< " volIntegrate(all) for alpha.liquid*rho.liquid = "
<< gSum(alphaLiquid*rhoLiquid*v) << endl
<< " volIntegrate(all) for dmdt = "
<< gSum(dmdt*v) << endl
<< endl;
#};
}
// ************************************************************************* //