[shinjanghosh]

Hi FOamers,

I am currently trying to write a boundary condition which reads like this :
q = −ρ(w · v) − ρ wnvn − μ(n · ∇wn) − ω1ρvnvn
(q is the variable which needs a patch update.)

I'm having trouble including the effective viscosity (which contains the molecular and turbulent viscosity) in the patch update. These are the expression I have used for the member functions in the boundary condition code :

const fvsPatchField<scalar>& phip =
patch().lookupPatchField<surfaceScalarField, scalar>("phi");

const fvsPatchField<scalar>& phiap =
patch().lookupPatchField<surfaceScalarField, scalar>("phia");

const fvPatchField<vector>& Up =
patch().lookupPatchField<volVectorField, vector>("U");

const fvPatchField<vector>& Uap =
patch().lookupPatchField<volVectorField, vector>("Ua");

const fvPatchField<scalar>& Tp =
patch().lookupPatchField<volScalarField, scalar>("T");


scalarField nueff = rasModel.nuEff()().boundaryField()

The last line is where I'm having the error as to how to access the nu_effective value from the solution itself like I did for the other quantities. Please let me know what you think.

Thanks,
Shinjan

[simrego]

Hi!


Just a guess, sorry if I'm wrong, but this should be fine for any turbulence model not just RAS:
turbulence->nuEff();
But maybe you will have to look up for "turbulence" somehow in the BC.