bEqn.H
Go to the documentation of this file.
1 tmp<fv::convectionScheme<scalar>> mvConvection
2 (
4  (
5  mesh,
6  fields,
7  phi,
8  mesh.divScheme("div(phi,ft_b_ha_hau)")
9  )
10 );
11 
12 volScalarField Db("Db", turbulence->muEff());
13 
14 if (ign.ignited())
15 {
16  // Calculate the unstrained laminar flame speed
17  // ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
19 
20  const volScalarField& Xi = flameWrinkling->Xi();
21 
22  // progress variable
23  // ~~~~~~~~~~~~~~~~~
24  volScalarField c("c", 1.0 - b);
25 
26  // Unburnt gas density
27  // ~~~~~~~~~~~~~~~~~~~
28  volScalarField rhou(thermo.rhou());
29 
30  // Calculate flame normal etc.
31  // ~~~~~~~~~~~~~~~~~~~~~~~~~~~
32 
33  // volVectorField n(fvc::grad(b));
35 
36  volScalarField mgb("mgb", mag(n));
37 
38  dimensionedScalar dMgb("dMgb", mgb.dimensions(), small);
39 
40  {
41  volScalarField bc(b*c);
42 
43  dMgb += 1.0e-3*
44  (bc*mgb)().weightedAverage(mesh.V())
45  /(bc.weightedAverage(mesh.V()) + small);
46  }
47 
48  mgb += dMgb;
49 
50  surfaceVectorField Sfhat(mesh.Sf()/mesh.magSf());
52  nfVec += Sfhat*(fvc::snGrad(b) - (Sfhat & nfVec));
53  nfVec /= (mag(nfVec) + dMgb);
54  surfaceScalarField nf("nf", mesh.Sf() & nfVec);
55  n /= mgb;
56 
57  #include "StCorr.H"
58 
59  // Calculate turbulent flame speed flux
60  // ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
61  surfaceScalarField phiSt("phiSt", fvc::interpolate(rhou*StCorr*St)*nf);
62 
63  #include "StCourantNo.H"
64 
65  Db = flameWrinkling->Db();
66 
67  // Create b equation
68  // ~~~~~~~~~~~~~~~~~
69  fvScalarMatrix bEqn
70  (
71  betav*fvm::ddt(rho, b)
72  + mvConvection->fvmDiv(phi, b)
73  + fvm::div(phiSt, b)
74  - fvm::Sp(fvc::div(phiSt), b)
75  - fvm::laplacian(Db, b)
76  ==
77  betav*fvModels.source(rho, b)
78  );
79 
80 
81  // Add ignition cell contribution to b-equation
82  // ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
83  #include "ignite.H"
84 
85  // Solve for b
86  // ~~~~~~~~~~~
87  bEqn.relax();
88 
90 
91  bEqn.solve();
92 
94 
95  Info<< "min(b) = " << min(b).value() << endl;
96 
97  if (composition.contains("ft"))
98  {
99  volScalarField& ft = composition.Y("ft");
100 
101  Info<< "Combustion progress = "
102  << 100*(1.0 - b)().weightedAverage(mesh.V()*ft).value() << "%"
103  << endl;
104  }
105  else
106  {
107  Info<< "Combustion progress = "
108  << 100*(1.0 - b)().weightedAverage(mesh.V()).value() << "%"
109  << endl;
110  }
111 
112  // Correct the flame-wrinkling
113  flameWrinkling->correct(mvConvection);
114 
115  St = Xi*Su;
116 }
fvMatrix< scalar > fvScalarMatrix
Definition: fvMatricesFwd.H:42
fluidReactionThermo & thermo
Definition: createFields.H:28
basicSpecieMixture & composition
tmp< GeometricField< Type, fvPatchField, volMesh > > div(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
Definition: fvcDiv.C:47
GeometricField< vector, fvsPatchField, surfaceMesh > surfaceVectorField
Calculates and outputs the mean and maximum Courant Numbers.
const dimensionedScalar b
Wien displacement law constant: default SI units: [m K].
Definition: createFields.H:27
Ostream & endl(Ostream &os)
Add newline and flush stream.
Definition: Ostream.H:251
Info<< "Calculating turbulent flame speed field St\"<< endl;volScalarField St(IOobject("St", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::AUTO_WRITE), flameWrinkling->Xi() *Su);multivariateSurfaceInterpolationScheme< scalar >::fieldTable fields
Definition: createFields.H:228
dimensionedScalar StCorr("StCorr", dimless, 1.0)
GeometricField< vector, fvPatchField, volMesh > volVectorField
Definition: volFieldsFwd.H:58
tmp< GeometricField< Type, fvPatchField, volMesh > > laplacian(const GeometricField< Type, fvPatchField, volMesh > &vf, const word &name)
Definition: fvcLaplacian.C:45
const dimensionedScalar c
Speed of light in a vacuum.
tmp< GeometricField< Type, fvPatchField, volMesh > > ddt(const dimensioned< Type > dt, const fvMesh &mesh)
Definition: fvcDdt.C:45
tmp< fv::convectionScheme< scalar > > mvConvection(fv::convectionScheme< scalar >::New(mesh, fields, phi, mesh.divScheme("div(phi,ft_b_ha_hau)")))
GeometricField< scalar, fvPatchField, volMesh > volScalarField
Definition: volFieldsFwd.H:57
Info<< "Creating field dpdt\"<< endl;volScalarField dpdt(IOobject("dpdt", runTime.timeName(), mesh), mesh, dimensionedScalar(p.dimensions()/dimTime, 0));Info<< "Creating field kinetic energy K\"<< endl;volScalarField K("K", 0.5 *magSqr(U));Info<< "Creating the unstrained laminar flame speed\"<< endl;autoPtr< laminarFlameSpeed > unstrainedLaminarFlameSpeed(laminarFlameSpeed::New(thermo))
Definition: createFields.H:90
dynamicFvMesh & mesh
Foam::fvConstraints & fvConstraints
autoPtr< BasicCompressibleMomentumTransportModel > New(const volScalarField &rho, const volVectorField &U, const surfaceScalarField &phi, const typename BasicCompressibleMomentumTransportModel::transportModel &transport)
phi
Definition: correctPhi.H:3
Info<< "Reading field U\"<< endl;volVectorField U(IOobject("U", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE), mesh);volScalarField rho(IOobject("rho", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::AUTO_WRITE), thermo.rho());volVectorField rhoU(IOobject("rhoU", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::NO_WRITE), rho *U);volScalarField rhoE(IOobject("rhoE", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::NO_WRITE), rho *(e+0.5 *magSqr(U)));surfaceScalarField pos(IOobject("pos", runTime.timeName(), mesh), mesh, dimensionedScalar(dimless, 1.0));surfaceScalarField neg(IOobject("neg", runTime.timeName(), mesh), mesh, dimensionedScalar(dimless, -1.0));surfaceScalarField phi("phi", fvc::flux(rhoU));Info<< "Creating turbulence model\"<< endl;autoPtr< compressible::momentumTransportModel > turbulence(compressible::momentumTransportModel::New(rho, U, phi, thermo))
Definition: createFields.H:94
dimensioned< Type > min(const dimensioned< Type > &, const dimensioned< Type > &)
static tmp< GeometricField< Type, fvsPatchField, surfaceMesh > > interpolate(const GeometricField< Type, fvPatchField, volMesh > &tvf, const surfaceScalarField &faceFlux, Istream &schemeData)
Interpolate field onto faces using scheme given by Istream.
const tmp< volScalarField::Internal > & Sp
Definition: alphaSuSp.H:7
bool constrain(fvMatrix< Type > &eqn) const
Apply constraints to an equation.
Foam::fvModels & fvModels
volScalarField Db("Db", turbulence->muEff())
tmp< GeometricField< typename outerProduct< vector, Type >::type, fvPatchField, volMesh >> reconstruct(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
messageStream Info
dimensioned< scalar > mag(const dimensioned< Type > &)
label n
const tmp< volScalarField::Internal > & Su
Definition: alphaSuSp.H:6
GeometricField< scalar, fvsPatchField, surfaceMesh > surfaceScalarField
tmp< GeometricField< Type, fvsPatchField, surfaceMesh > > snGrad(const GeometricField< Type, fvPatchField, volMesh > &vf, const word &name)
Definition: fvcSnGrad.C:45
autoPtr< XiModel > flameWrinkling
Create the flame-wrinkling model.
Definition: createFields.H:202