[slimak]

Dear Foamers,

I am trying to implement my own solver in OpenFOAM 10. The solver is a version of scalarTransportFoam, where I want to introduce turbulent contribution to the dispersion term (for heat transfer this would be something like turbulent Prandtl, Prt, number or Schmidt, Sct, number for species transfer). This would make the laplacian term look like fvm::laplacian(DT+nut/Prt, T)
with the extra term being turbulent diffusivity (nut) divided by Prt. For simplicity I assume Prt is constant. My struggle is to introduce the nut to the equation via createFields.H.
In previous version of OF (I guess 5 or 6, OF org branch by CFDdirect), I was able to compile code with
singlePhaseTransportModel laminarTransport(U, phi);


autoPtr<incompressible::turbulenceModel> turbulence //turbulentTransportModel
(
incompressible::turbulenceModel::New(U, phi, laminarTransport)
);
volScalarField nut = turbulence->nut();

which created nut. But this is not the case in the current version. My guess is, that it has something to do with
#include"fvModels.H"

replacing
#include "singlePhaseTransportModel.H"
#include "turbulentTransportModel.H"
but I was not able to solve the issue. There are some threads on this topic (How to access nu and nut values for a Turbulent Solver), but they seem to be from 2014-2017. My source code for the modified solver are included. Do you have any hints?

Thank you in advance.
p { margin-bottom: 0.1in; line-height: 115% }a:link { so-language: zxx }