Problems about the `eddyDissipation` model

[TimoLin]

Description

• Compiling problem
  If one want to use the eddyDisspation modek, the following thermo class should be modified:

From:

deepflame-dev/applications/solvers/dfLowMachFoam/createFields.H

Lines 7 to 9 in 1be82b6

	// fluidThermo* pThermo = new hePsiThermo<psiThermo, CanteraMixture>(mesh, word::null);
	fluidThermo* pThermo = new heRhoThermo<rhoThermo, CanteraMixture>(mesh, word::null);
	fluidThermo& thermo = *pThermo;

To:

 fluidThermo* pThermo = new heRhoThermo<rhoThermo, dfSingleStepReactingMixture>(mesh, word::null); 

• Simulation problem
  When the eddyDissipation model is set in combustionProperties, the simulation throws the following error:

[7] --> FOAM FATAL ERROR:
[7]
    lookup of turbulenceProperties from objectRegistry region0 successful
    but it is not a LES, it is a kEpsilon
[7]
[7]     From function const Type& Foam::objectRegistry::lookupObject(const Foam::word&) const [with Type = Foam::LESModel<Foam::EddyDiffusivity<Foam::ThermalDiffusivity<Foam::CompressibleTurbulenceModel<Foam::fluidThermo> > > >]
[7]     in file /public/home/chenzhi/dftorch/OpenFOAM-7/src/OpenFOAM/lnInclude/objectRegistryTemplates.C at line 178.
[7]
FOAM parallel run aborting

Releated code:

deepflame-dev/src/dfCombustionModels/eddyDissipationModel/eddyDissipationModel.C

Lines 93 to 104 in 1be82b6

	eddyDissipationModel<ReactionThermo>::rtDiff() const
	{
	const CanteraMixture& mixture_ = dynamic_cast<const CanteraMixture&>(this->thermo());
	const volScalarField& YO2 = mixture_.Y("O2");
	const compressible::LESModel& lesModel =
	YO2.db().lookupObject<compressible::LESModel>
	(
	turbulenceModel::propertiesName
	);
	return Cd_*this->thermo().alpha()/this->rho()/sqr(lesModel.delta());
	}

[TimoLin]

Possible fix patch

Use mesh.V()^{1/3} to calculate the delta directly?