ftEqn.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 if (composition.contains("ft"))
13 {
14  volScalarField& ft = composition.Y("ft");
15 
16  fvScalarMatrix ftEqn
17  (
18  fvm::ddt(rho, ft)
19  + mvConvection->fvmDiv(phi, ft)
20  - fvm::laplacian(turbulence->alphaEff(), ft)
21  ==
22  fvOptions(rho, ft)
23  );
24 
25  fvOptions.constrain(ftEqn);
26 
27  ftEqn.solve();
28 
29  fvOptions.correct(ft);
30 }
Info<< "Creating field dpdt\n"<< endl;volScalarField dpdt(IOobject("dpdt", runTime.timeName(), mesh), mesh, dimensionedScalar("dpdt", p.dimensions()/dimTime, 0));Info<< "Creating field kinetic energy K\n"<< endl;volScalarField K("K", 0.5 *magSqr(U));volScalarField p_rgh(IOobject("p_rgh", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE), mesh);p_rgh=p-rho *gh;mesh.setFluxRequired(p_rgh.name());multivariateSurfaceInterpolationScheme< scalar >::fieldTable fields
Definition: createFields.H:127
tmp< GeometricField< Type, fvPatchField, volMesh > > laplacian(const GeometricField< Type, fvPatchField, volMesh > &vf, const word &name)
Definition: fvcLaplacian.C:45
autoPtr< compressible::turbulenceModel > turbulence
Definition: createFields.H:23
tmp< GeometricField< Type, fvPatchField, volMesh > > ddt(const dimensioned< Type > dt, const fvMesh &mesh)
Definition: fvcDdt.C:45
dynamicFvMesh & mesh
autoPtr< BasicCompressibleTurbulenceModel > New(const volScalarField &rho, const volVectorField &U, const surfaceScalarField &phi, const typename BasicCompressibleTurbulenceModel::transportModel &transport, const word &propertiesName)
fv::IOoptionList & fvOptions
fvMatrix< scalar > fvScalarMatrix
Definition: fvMatricesFwd.H:42
basicMultiComponentMixture & composition
Definition: createFields.H:35
GeometricField< scalar, fvPatchField, volMesh > volScalarField
Definition: volFieldsFwd.H:52
surfaceScalarField & phi
tmp< fv::convectionScheme< scalar > > mvConvection(fv::convectionScheme< scalar >::New( mesh, fields, phi, mesh.divScheme("div(phi,ft_b_ha_hau)") ))