5 U -= fvc::reconstruct(
mesh.phi());
15 U += fvc::reconstruct(
mesh.phi());
28 const volScalarField& porosity =
tporosity.cref();
29 rhoPhi *= scalar(1)/fvc::interpolate(porosity);
37Info<<
"Phase-1 volume fraction = "
const volScalarField & alpha1
surfaceScalarField rho2f(fvc::interpolate(rho2))
surfaceScalarField rho1f(fvc::interpolate(rho1))
const bool porosityEnabled(porosityProperties.getOrDefault< bool >("porosityEnabled", false))
tmp< volScalarField > tporosity
messageStream Info
Information stream (stdout output on master, null elsewhere)
Ostream & endl(Ostream &os)
Add newline and flush stream.
label min(const labelHashSet &set, label minValue=labelMax)
Find the min value in labelHashSet, optionally limited by second argument.
Info<< "Creating temperaturePhaseChangeTwoPhaseMixture\n"<< endl;autoPtr< temperaturePhaseChangeTwoPhaseMixture > mixture
isoAdvection advector(alpha1, phi, U)