[ivan_cozza]

Hi Foamers,
I need to write a b.c. for a scalar field for which the boundary field must be an accurate interpolation of the internal field.
Is there any ready to use interpolation tool?
I'm managing to use linearInterpolation(), like in phi construction, but does it use only linear scheme? What I have to do to use higher order schemes (cubic)?

Thanks, Ivan

[ivan_cozza]

Ok, yesterday I wrote down some errors (the stuff about linearInterpolate).

Now my question is: is possible using fvc::interpolate() to extrapolate the boundary field of a volScalarField starting from the internal field values?

To add some infos, I need to do it as in a piece of my codes there's a scalar field that's not a solution of a PDE, but is just evaluated algebrically, and I need some good extimation of its boundary values in order to do some other calculations (zeroGradient b.c. it's not enough and the field is time-varying, so I cannot impose at the beginnig the value).

Thank for the help!!!
Ivan

[aliqasemi]

This is an old post but:

for any field X, the function fvc::interpolate(X) should consider the boundary conditions specified for X to interpolate it to boundary faces. i.e. if you assign a correct boundary condition it will evaluate the values at boundary faces correctly. I am not so good at the boundary conditions to help you in choosing the correct BC. By default, it uses zeroGradient BC, I think.

[BSengupta]

Hi! I am sorry to bump into such old post. But I am getting similar problem.


I am having a velocity quantity named URel and b.c. at wall is zeroGradient. I need to calculate a quanity name betta, which is given by:

betta = rho*URel.n/(U_Inf*rho_Inf)

where n is unit normal from wall.



I am defining betta_surf as:

surfaceScalarField betta_surf
     (
         IOobject
         (
             "betta_surf",
             runTime.timeName(),
             mesh,
             IOobject::NO_READ,
             IOobject::NO_WRITE
         ),
         (fvc::interpolate(URel))&(mesh.Sf()/mesh.magSf())
     );

I am trying to print its value at wall by:

Info<<"betta surf : "<<betta_surf.boundaryField()[wallID]<<endl;

But I am always getting it to be zero

betta surf : type            calculated;
 value           uniform 0;

although the cell center value of all cells adjacent to wall is non zero.


What mistake am I making, can anyone kindly let me know if possible. If any other information is required I shall be happy to share.


Thank you.
p { margin-bottom: 0.1in; line-height: 115% } a:link { so-language: zxx }
p { margin-bottom: 0.1in; line-height: 115% } a:link { so-language: zxx }
p { margin-bottom: 0.1in; line-height: 115% } a:link { so-language: zxx }