alphaDiffusionEqn.H
Go to the documentation of this file.
1 {
2  fvScalarMatrix alpha1Eqn
3  (
7  (
8  volScalarField("Dab", Dab + alphatab*turbulence->nut()),
9  alpha1
10  )
11  );
12 
13  alpha1Eqn.solve();
14 
15  alpha2 = 1.0 - alpha1;
16  rhoPhi += alpha1Eqn.flux()*(rho1 - rho2);
17 }
18 
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
fvMatrix< scalar > fvScalarMatrix
Definition: fvMatricesFwd.H:42
GeometricField< scalar, fvPatchField, volMesh > volScalarField
Definition: volFieldsFwd.H:52
rho1
Definition: pEqn.H:124
rho2
Definition: pEqn.H:125
volScalarField & alpha1
Definition: createFields.H:15