[landau]

Hello!

I'm a very beginner in OpenFoam programming, and I would like to ask the community some help with an issue I have not been able to find any information of in this forum, or even in google.

First of all I am modelling the wind over complex terrain. I have introduced some source terms in the equations that depend on the vertical coordinate. The problem is that I need to fix the reference for the z coordinate over the ground and not over the plane z=0. Therefore I need to subtract from every cell z coordinate in the volume mesh the corresponding height of the coincident (x,y) point of the ground.

I have think of the function FoamatchWave:atchWave().distance()

For this reason I've used (I don't think that properly) this function as follows:

  
label patchID = mesh.boundaryMesh().findPatchID("terrain");
scalarField d;//distances are stored in this variable
patchWave w(mesh, patchID, true);
d = w.distance();

The only distance it has computed is a negative one:

204901{-1e+15}

and obviously it gives an error but I do not know why.

Thank you in advance.

Héctor

[clapointe]

Have you looked into wallDist (https://github.com/OpenFOAM/OpenFOAM...ist/wallDist.H)? It should be what you need.

Caelan

[landau]

Thank you so much for your suggestion. The fact is that I saw this function but I didn't know how to extract a new volScalarField from the function wallDist. Could you guide me a little?

I know that the functions that wallDist has return a newmesh with updateMesh. But my goal is obtain a volscalarField with all the cell distances to the required patch.

Thank you so much for your patience.

Héctor

[clapointe]

Fair warning -- I'm not sure if I've used it before, but I've seen it in e.g. turbulence models. I pulled this snippet from the kOmegaSST base... it might help :

y_(wallDist::New(this->mesh_).y())

.

Caelan

[landau]

Thank you so much Caelan. So much appreciated! At the end it was much easier than I thought.