CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Programming & Development

Calculate distance from vol to patch cells

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   June 10, 2019, 07:31
Default Calculate distance from vol to patch cells
  #1
New Member
 
HBR
Join Date: Jun 2018
Posts: 6
Rep Power: 4
landau is on a distinguished road
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:
Code:
  
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:
Code:
204901{-1e+15}
and obviously it gives an error but I do not know why.

Thank you in advance.

Héctor
landau is offline   Reply With Quote

Old   June 10, 2019, 08:11
Default
  #2
Senior Member
 
Join Date: Aug 2015
Posts: 440
Rep Power: 11
clapointe is on a distinguished road
Have you looked into wallDist (https://github.com/OpenFOAM/OpenFOAM...ist/wallDist.H)? It should be what you need.

Caelan
clapointe is offline   Reply With Quote

Old   June 10, 2019, 08:21
Default
  #3
New Member
 
HBR
Join Date: Jun 2018
Posts: 6
Rep Power: 4
landau is on a distinguished road
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
landau is offline   Reply With Quote

Old   June 10, 2019, 08:43
Default
  #4
Senior Member
 
Join Date: Aug 2015
Posts: 440
Rep Power: 11
clapointe is on a distinguished road
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 :

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

Caelan
clapointe is offline   Reply With Quote

Old   June 11, 2019, 03:46
Default
  #5
New Member
 
HBR
Join Date: Jun 2018
Posts: 6
Rep Power: 4
landau is on a distinguished road
Thank you so much Caelan. So much appreciated! At the end it was much easier than I thought.
landau is offline   Reply With Quote

Reply

Tags
distance, patch, vertical distance

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
[snappyHexMesh] SnappyHexMesh running killed! Mark JIN OpenFOAM Meshing & Mesh Conversion 3 March 12, 2020 18:50
y+ and u+ values with low-Re RANS turbulence models: utility + testcase florian_krause OpenFOAM 113 June 25, 2019 01:58
[snappyHexMesh] sHM layer process keeps getting killed MBttR OpenFOAM Meshing & Mesh Conversion 4 August 15, 2016 04:21
[snappyHexMesh] snappyHexMesh won't work - zeros everywhere! sc298 OpenFOAM Meshing & Mesh Conversion 2 March 27, 2011 22:11
[snappyHexMesh] snappyHexMesh aborting Tobi OpenFOAM Meshing & Mesh Conversion 0 November 10, 2010 04:23


All times are GMT -4. The time now is 13:21.