Commit a0eeb3f0 authored by andy's avatar andy
Browse files

ENH: Removed unused entry from tutorial input dictionary

parent 63177c9d
......@@ -25,5 +25,7 @@
EaEqn.solve();
sources.correct(hea);
thermo.correct();
}
......@@ -2,7 +2,7 @@ if (ign.ignited())
{
volScalarField& heau = thermo.heu();
solve
fvScalarMatrix heauEqn
(
fvm::ddt(rho, heau) + mvConvection->fvmDiv(phi, heau)
+ (fvc::ddt(rho, K) + fvc::div(phi, K))*rho/thermo.rhou()
......@@ -23,5 +23,14 @@ if (ign.ignited())
// A possible solution would be to solve for ftu as well as ft.
//- fvm::div(muEff*fvc::grad(b)/(b + 0.001), heau)
//+ fvm::Sp(fvc::div(muEff*fvc::grad(b)/(b + 0.001)), heau)
==
sources(heau)
);
sources.constrain(heauEqn);
heauEqn.solve();
sources.correct(heau);
}
......@@ -15,5 +15,7 @@
if (pimple.momentumPredictor())
{
solve(UEqn == -fvc::grad(p));
sources.correct(U);
K = 0.5*magSqr(U);
}
......@@ -71,6 +71,8 @@ if (ign.ignited())
bEqn.solve();
sources.correct(b);
Info<< "min(b) = " << min(b).value() << endl;
......@@ -167,6 +169,8 @@ if (ign.ignited())
SuEqn.solve();
sources.correct(Su);
// Limit the maximum Su
// ~~~~~~~~~~~~~~~~~~~~
Su.min(SuMax);
......@@ -252,6 +256,8 @@ if (ign.ignited())
XiEqn.solve();
sources.correct(Xi);
// Correct boundedness of Xi
// ~~~~~~~~~~~~~~~~~~~~~~~~~
Xi.max(1.0);
......
......@@ -25,4 +25,6 @@ if (composition.contains("ft"))
sources.constrain(ftEqn);
ftEqn.solve();
sources.correct(ft);
}
......@@ -10,5 +10,7 @@
if (pimple.momentumPredictor())
{
solve(UEqn == -fvc::grad(p));
sources.correct(U);
K = 0.5*magSqr(U);
}
......@@ -27,5 +27,6 @@
)
);
sources.correct(U);
K = 0.5*magSqr(U);
}
......@@ -39,6 +39,8 @@ tmp<fv::convectionScheme<scalar> > mvConvection
YiEqn.solve(mesh.solver("Yi"));
sources.correct(Yi);
Yi.max(0.0);
Yt += Yi;
}
......@@ -82,6 +84,8 @@ tmp<fv::convectionScheme<scalar> > mvConvection
EEqn.solve();
sources.correct(he);
thermo.correct();
Info<< "min/max(T) = "
......
......@@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | Copyright (C) 2011 OpenFOAM Foundation
\\ / A nd | Copyright (C) 2011-2012 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
......@@ -43,6 +43,8 @@ Description
sources.constrain(rhoEqn);
rhoEqn.solve();
sources.correct(rho);
}
// ************************************************************************* //
......@@ -28,6 +28,8 @@
EEqn.solve();
sources.correct(he);
thermo.correct();
Info<< "min/max(T) = "
......
......@@ -15,5 +15,7 @@
if (pimple.momentumPredictor())
{
solve(UEqn == -fvc::grad(p));
sources.correct(U);
K = 0.5*magSqr(U);
}
......@@ -37,6 +37,8 @@ tmp<fv::convectionScheme<scalar> > mvConvection
YiEqn.solve(mesh.solver("Yi"));
sources.correct(Yi);
Yi.max(0.0);
Yt += Yi;
}
......
......@@ -27,7 +27,6 @@ momentumSource
{
fieldNames (U);
Ubar ( 0.1335 0 0 );
gradP 0;
}
}
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment