[alimea]

Hi all

I want to define a volScalarField in one of the viscous model codes. I saw this function that is set 0 for initial value of a dimensionedScalar:

virtual const dimensionedScalar etaP()
{
return (dimensionedScalar ("zeroU", dimensionSet(...), 0) )
}

but it does't work for a volScalarField :

virtual const volScalarField etaP()
{
return (volScalarField ("zeroU", dimensionSet(...), 0) )
}

this problem is in createFields file in solvers' code. When I want to define a volScalarField in createFields file, compiler doesn't allow me to set it in special value like 10 for all of the field like:

volScalarField AA = 0;

If I write

virtual const volScalarField etaP()
{
volScalarField A = dimensionedScalar ("zeroU", dimensionSet(...), 0);
return A;
}

compiler give me this error:

error: conversion from 'Foam::dimensionedScalar {aka Foam::dimensioned<double>}' to non-scalar type 'Foam::volScalarField {aka Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>}' requested

Could you plz help me to do that?

Thanks

[Cyp]

See slides 94 to 102

Cheers,
Cyprien

[alimea]

Originally Posted by Cyp

See slides 94 to 102

Cheers,
Cyprien

Dear Cyprien

Thanks for your reply.

Those are very useful, but at the middle of a viscousity model code, I can't use non of them!
If I want to define as

volScalarField etaP
 (
                IOobject
                (
                    "etaP",
                    runTime.timeName(),
                    mesh,
                    IOobject::NO_READ,
                    IOobject::NO_WRITE
	   ),
                mesh
);

the runTime and mesh objects are not defined!

If I want to use

volScalrField etaP ("etaP", 0.*p);

p or any other volScalarField isn't defined here!

So, How to do that?

Regards,

[jpeter3]

You need to get mesh and time objects, they do not exist in every class under this names.

Maybe try to Access them through velocity field, something like:

U_.mesh() instead of mesh

and

U_.time() instead of runTime

[alimea]

Originally Posted by jpeter3

You need to get mesh and time objects, they do not exist in every class under this names.

Maybe try to Access them through velocity field, something like:

U_.mesh() instead of mesh

and

U_.time() instead of runTime

Thanks for your reply.

could you plz explain it more?
I have written it as:

volScalarField etaP
 (
                IOobject
                (
                    "etaP",
                    runTime.timeName(),
                    mesh,
                    IOobject::NO_READ,
                    IOobject::NO_WRITE
	   ),
                mesh
);

that runTime and mesh are not defined!
so as you mentioned I have changed them to:

volScalarField etaP
 (
                IOobject
                (
                    "etaP",
                    U_.time(),
                    U_.mesh(),
                    IOobject::NO_READ,
                    IOobject::NO_WRITE
	   ),
                U_.mesh()
);

but it gives me this error:

error: no matching function for call to 'Foam::IOobject::IOobject(const char [4], const Foam::Time&, const Mesh&, Foam::IOobject::readOption, Foam::IOobject::writeOption)'
),
^

Regards,

[jpeter3]

Yes, that is exactly what I meant. Only one thing should be changed:

in original, you have written: runTime.timeName()

That is a call to function timeName() of the runTime.

So, when you replace runTime by U_.time(), you still need to call the function timeName():

...

IOobject
(
"etaP",
U_.time().timeName(),
U_.mesh(),
IOobject::NO_READ,
IOobject::NO_WRITE
),

...

[alimea]

Originally Posted by jpeter3

Yes, that is exactly what I meant. Only one thing should be changed:

in original, you have written: runTime.timeName()

That is a call to function timeName() of the runTime.

So, when you replace runTime by U_.time(), you still need to call the function timeName():

...

IOobject
(
"etaP",
U_.time().timeName(),
U_.mesh(),
IOobject::NO_READ,
IOobject::NO_WRITE
),

...

Thenk you,
but it gives me again this error:

error: no matching function for call to 'Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::GeometricField(Foam::IOobject, <unresolved overloaded function type>)'
);

what does it say?

Regards,

[alimea]

Originally Posted by jpeter3

Yes, that is exactly what I meant. Only one thing should be changed:

in original, you have written: runTime.timeName()

That is a call to function timeName() of the runTime.

So, when you replace runTime by U_.time(), you still need to call the function timeName():

...

IOobject
(
"etaP",
U_.time().timeName(),
U_.mesh(),
IOobject::NO_READ,
IOobject::NO_WRITE
),

...

Dear JPeternel

That is solved!
I was wrong in last "mesh" component.

Thank you