[CGuaq]

Hello,

I try to modify stochastic dispersion model on Open FOAM to improved my results concerning particles deposition. I found a tutorial from Chalmers University:

http://www.tfd.chalmers.se/~hani/kur..._Xu_Report.pdf

I changed the stochasticDispersionRAS.C and the stochasticDispersionRAS.H like in the tutorial but it used Open FOAM 4 and I use Open FOAM 9. My problem is on the following part of the code (from line 96 to 103):


const word turbName= IOobject::groupName
(
turbulenceModel:: propertiesName,
this->owner().U().group()
);
const turbulenceModel& turbModel = obr.lookupObject<turbulenceModel>(turbName);


I changed all "turbulenceModel" by "momentumTransportModel" concidering the modifications between Open FOAM 4 et 9 but "propertiesName" is no longer defined in momentumTransportModel.C and momentumTransportModel.H like it was on Open FOAM 4. I don't know how to modify this part of the code to make it works on Open FOAM 9.
I don't have enough experience developing code on OpenFoam yet, so please can someone help me ?

Kind regards,

[francois]

Any idea where to find some helpful pointers on this particular point ?

[SmileMax]

Originally Posted by CGuaq

const word turbName= IOobject::groupName
(
turbulenceModel:: propertiesName,
this->owner().U().group()
);
const turbulenceModel& turbModel = obr.lookupObject<turbulenceModel>(turbName);

In fact, turbulenceModel:: propertiesName corresponds to the generic turbulence dictionary name, such as turbulenceProperties in OpenFOAM(.com) or momentumTransport in OpenFOAM(version 8+, .org) , and this->owner().U().group() corresponds to your phase name.

so you can use the following simplified form：

const word turbName= IOobject::groupName
(
    "momentumTransport.<phase name>"
);

[CGuaq]

Thanks for your reply,

Unfortunatly it does not work, are you sure it is momentumTransport.<phase name> and not momentumTransportModel.<phase name> ? Because OpenFOAM suggested me to use "momentumTransportModel_H" instead of "momentumTransport" when I compile the code.
Also, how can I find my phase name to replace in your simplified form because "phase name" is not declared.
Thank you in advance for your answer

[SmileMax]

Originally Posted by CGuaq

Thanks for your reply,

Unfortunatly it does not work, are you sure it is momentumTransport.<phase name> and not momentumTransportModel.<phase name> ? Because OpenFOAM suggested me to use "momentumTransportModel_H" instead of "momentumTransport" when I compile the code.
Also, how can I find my phase name to replace in your simplified form because "phase name" is not declared.
Thank you in advance for your answer

You'd better be able to provide the error message output by the compiler.

[CGuaq]

Sorry it is my first time on CFD Online.

I think I found the solution by compararing files like shearStress.C or yPlus.C, so I changed my previous code into:

const word turbName= IOobject::groupName
                (
                    momentumTransportModel::typeName,
                    phaseName_
                );

The error message has disappeared since this modification. But now I have a new error message concerning the following part of the code (line 43 to 55):

yP_
    	(
    	IOobject
    	(
    	"yP",
    	this->owner_.mesh().time().typeName(),
    	this->owner_.mesh(),
    	IOobject::NO_READ,
    	IOobject::AUTO_WRITE
    	),
    	this->owner_.mesh(),
    	dimensionedScalar("yP", dimless, 0.0)
    	)

This is the error message:
error_message.jpg

i'm not sure what that means.

[SmileMax]

Originally Posted by CGuaq

Sorry it is my first time on CFD Online.

I think I found the solution by compararing files like shearStress.C or yPlus.C, so I changed my previous code into:

const word turbName= IOobject::groupName
                (
                    momentumTransportModel::typeName,
                    phaseName_
                );

The error message has disappeared since this modification. But now I have a new error message concerning the following part of the code (line 43 to 55):

yP_
        (
        IOobject
        (
        "yP",
        this->owner_.mesh().time().typeName(),
        this->owner_.mesh(),
        IOobject::NO_READ,
        IOobject::AUTO_WRITE
        ),
        this->owner_.mesh(),
        dimensionedScalar("yP", dimless, 0.0)
        )

This is the error message:
Attachment 91684

i'm not sure what that means.

This seems to be a basic syntax error in the C++ language.

[CGuaq]

Do you have any idea what i can modify to make it works on OpenFOAM 9 ? I tried modifications but nothing works.
thank you in advance