51 int main(
int argc,
char *argv[])
59 #include "createFields.H"
66 Info<<
"\nStarting time loop\n" <<
endl;
68 while (runTime.loop())
70 Info<<
"Time = " << runTime.userTimeName() <<
nl <<
endl;
77 divR.source() =
flowMask & divR.source();
99 U += (
Ubar - UbarStar);
100 gradP += (
Ubar - UbarStar)/(1.0/
UEqn.
A())().weightedAverage(mesh.V());
111 if (runTime.writeTime())
116 Info<<
"ExecutionTime = " << runTime.elapsedCpuTime() <<
" s"
117 <<
" ClockTime = " << runTime.elapsedClockTime() <<
" s"
static void noParallel()
Remove the parallel options.
const Type & value() const
Return const reference to value.
Finite volume constraints.
bool constrain(fvMatrix< Type > &eqn) const
Apply constraints to an equation.
A special matrix type and solver, designed for finite volume solutions of scalar equations....
void relax(const scalar alpha)
Relax matrix (for steady-state solution).
SolverPerformance< Type > solve(const dictionary &)
Solve segregated or coupled returning the solution statistics.
tmp< volScalarField > A() const
Return the central coefficient.
tmp< fvMatrix< Type > > source(const VolField< Type > &field) const
Return source for an equation.
virtual void correct()
Correct the fvModels.
Abstract base class for all fluid physical properties.
int main(int argc, char *argv[])
Ostream & endl(Ostream &os)
Add newline and flush stream.
dimensionedVector gradP("gradP", dimensionSet(0, 1, -2, 0, 0), Zero)
autoPtr< incompressible::momentumTransportModel > turbulence(incompressible::momentumTransportModel::New(U, phi, viscosity))
dimensionedVector Ubar("Ubar", dimVelocity, viscosity)