[k.farnagh]

Hi dears,

I want to lookup or invoke Ub (inflow air velocity at inlet) to drag model (.c file).

Ub is valued in "0" directory and is volVectorField type but its not important for me to keep type class, in other word scalar output is better for me.

I apply "twoPhaseEulerFoam" solver of 2.1.1 OF and drag model is "GidaspowErgunWenYu".

Would you please guide me.

[k.farnagh]

no body no suggestion?

[epi_c]

Hi Kasra,

Maybe you can try

const label patchI = U.mesh().boundaryMesh().findPatchID("inlet");
const vectorField Uinlet = U.boundaryField()[patchI];

[k.farnagh]

Hi Dear Jianzhi,

Thanks for reply, As I'm newcomer in coding would you please explain where should I put them. .c file of drag model --> member function or somewhere else?
furthermore should I include some .H file, change Constructor or ...?

Sincerely Yours
Kasra

[k.farnagh]

Hi,

I've added below code to end of Member Function of "GidaspowErgunWenYu.c" and nothing else.

Info<< "Looking up field Ub\n" << endl;
const volVectorField& Ub = Ub.mesh().lookupObject<volVectorField>("Ub");

Info<< "Reading inlet velocity uInlet\n" << endl;
scalar UInlet = 0.0;
const label inletI = Ub.mesh().boundaryMesh().findPatchID("inlet");
const fvPatchVectorField& fvp = Ub.boundaryField()[inletI];

if (fvp.size())
{
UInlet = fvp[0].y();
}
reduce(UInlet, maxOp<scalar>());

I've compiled new drag model and encountered just below warning

dragModels/GidaspowErgunWenYulogvelo/GidaspowErgunWenYulogvelo.C:98:36: warning: ‘Ub’ may be used uninitialized in this function [-Wmaybe-uninitialized] const volVectorField& Ub = Ub.mesh().lookupObject<volVectorField>("Ub");

However I've ran a case with new drag model and it has stopped calculation when looking for Ub field.

would you please guide me to solve it. thanks again Jianzhi.

Sincerely Yours,
Kasra