[usergk]

Hello FOAMers,

I am trying to modify the following equation:


solve
(
fvm::ddt(rho, Yi)
+ mvConvection->fvmDiv(phi, Yi)
- fvm::laplacian
(
(TF*E*(turbulence->mu()+ turbulence->mut())),
Yi, "laplacian(muEff,Yi)"
)
==
parcels.Srho(i)
+ E*chemistry.RR(i)().dimensionedInternalField()/TF
+ massSource.Su(i)
);


Here TF is a scalar (constant value in each cell of the domain); however, I want to replace it by a volScalarField (a dynamic value calculated at each time step and varies across the cells).


Can someone provide the syntax of how to go about it? I've tried a few ways but they don't seem to be working.

Many thanks,
gk

[santiagomarquezd]

Hi, see in createFields.H file and add a new TF field.

Regards.