[CGuaq]

Hello everyone,

I'm working on a turbulent particulate flow in a horizontal 2D channel, I use v2f model for turbulence and I would like to use v2 field values (returned by v2f model) in my dispersion model. For now I use correlations from DNS to calculate velocity fluctuations on fluid but I want to be more general with v2 values.

I tried to add v2 values in file "momentumTransportModel.H" like this (same way as for k, epsilon, ...):

virtual tmp<volScalarField> v2() const = 0;

And I also modified DispersioRASModel.C by adding :

template<class CloudType>
Foam::tmp<Foam::volScalarField>
Foam::DispersionRASModel<CloudType>::v2Model() const
{
    const objectRegistry& obr = this->owner().mesh();
    const word turbName =
        IOobject::groupName
        (
            momentumTransportModel::typeName,
            this->owner().U().group()
        );

    if (obr.foundObject<momentumTransportModel>(turbName))
    {
        const momentumTransportModel& model =
            obr.lookupObject<momentumTransportModel>(turbName);
        return model.v2();
    }
    else
    {
        FatalErrorInFunction
            << "Turbulence model not found in mesh database" << nl
            << "Database objects include: " << obr.sortedToc()
            << abort(FatalError);

        return tmp<volScalarField>(nullptr);
    }
}

To finish, in my dispersion model code I use following line to call v2 values in each cell (again like k and espilon in stochasticDispersionRAS model)

const scalar v2_ = this->v2Ptr_->primitiveField()[celli];

I can compile my code and launch my simulation but when I display (on parafoam) what returns me this line I have a very strange scalar field ...

So can someone help me to find an other solution to access to my v2 field values and use them in my dispersion code ?? Or tell me if I made a mistake in my code.

Thanks in advance for your help !!

Best regards

[CGuaq]

Problem solved

[Tobermory]

Can you consider adding your solution, instead of just saying that you solved it? It will help others in the future if they have the same problem, and therefore helps increase the utility of the site? Thanks in advance.

[CGuaq]

Yes of course, Sorry.
To solve my problem and add v2 field values, I just added following lines in my turbulent dispersion file:

const objectRegistry& obr = this->owner().mesh();
const volScalarField& v2_ = obr.lookupObject<volScalarField>("v2");    
const scalar v2 = v2_[celli];

I guess it's possible to get others volScalarFields "saved" in the mesh with this method.