[som]

Hello. I have some scalar field and need to find its value for all cells near some boundary. I need to have array or something like that with values on every timestep. How can i do it? I use some modification of pimpleDyMFoam solver.

[nimasam]

your question is vague, however maybe you can use scalarField type to save your data in it

[som]

Thanks for reply.

I'm asking not about type where i can save my data. Actualy I need to printf value of scalar field near the specified boundary.

[nimasam]

you can show it in terminal with Info
but if you want to save it in file you show use
IOStream

[som]

yep, i understand that. But the question is how to find cells which are near the specific boundary.

What I have
1) Boundary
2) Field in space

I need to know (i.e. print it to console) field value in the boundary layer.

[doubtsincfd]

const fvPatchList& patches = mesh.boundary();
forAll(patches,patchi)
{
const fvPatch& currPatch = patches[patchi];
if(currPatch.name()=="nameOfThePatch")
{
Info<<U.boundaryField()[patchi].patchInternalField()<<endl;
}
}

Let me know if this works
patchInternalField() gives you the values at cell centers of cells near the boundary

[som]

doubtsincfd

Thanks for your response. It seems you've shown me right direction.

[som]

It works. One more stupid question. How can i take the maximum the value of a patch?

[doubtsincfd]

there is an in-built function max(fieldName). I havent used it but it should be pretty straightforward.

http://www.cfd-online.com/Forums/ope...alarfield.html

link might help

[som]

It works.Thank you!