CFD Online URL
[Sponsors]
Home > Forums > OpenFOAM Bugs

wallDist error in OF-1.5.x

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

Reply
 
LinkBack Thread Tools Display Modes
Old   August 19, 2010, 12:27
Default wallDist error in OF-1.5.x
  #1
Member
 
Elisabet Mas de les Valls
Join Date: Mar 2009
Location: Barcelona, Spain
Posts: 60
Rep Power: 7
elisabet is on a distinguished road
Hi everybody,

I've experienced some problems with the wallDist function for the case shown below. It consists on two banana shaped channels (inlet at left and outlet at right) and the U-bend that connects both channels. I need to evaluate the wall distance from the patch market in black (at the left of the figure)

banana_shape.png

The wallDist function seems to return the good wall distance for the inlet channel but for the outlet channel it somehow follows what would be a fluid trajectory:

banana_wallDist.png

In order to obtain the right wall distance I've had to create a new geometry without the separation of both channels and, over it, apply the wallDist function. Afterwards, interpolate the result to the real geometry:

unifiedBanana_wallDist.png

In my case, as I needed the wall distance as a pre-processor and not during the simulation, the above mentioned solution was good enough (despite not elegant at all...).

I guess the solution is not straightforward but for those who like maths challenges... here you are!

Regards,

Elisabet
elisabet is offline   Reply With Quote

Old   August 19, 2010, 13:06
Default
  #2
Super Moderator
 
Mattijs Janssens
Join Date: Mar 2009
Posts: 1,416
Rep Power: 15
mattijs is on a distinguished road
wallDist calculates the distance to any 'wall' patches. Where are your walls?

Can you post/send the mesh?

Mattijs
mattijs is offline   Reply With Quote

Old   August 19, 2010, 13:41
Default
  #3
Member
 
Elisabet Mas de les Valls
Join Date: Mar 2009
Location: Barcelona, Spain
Posts: 60
Rep Power: 7
elisabet is on a distinguished road
Dear Matiijs,

For the pre-processor I modified the boundary file in order to have one unique wall type boundary condition: the one that I want as a reference for the wallDist function (shown in black in the first figure).

I'll send you the mesh in a private mail.


Elisabet
elisabet is offline   Reply With Quote

Old   August 20, 2010, 07:11
Default
  #4
Super Moderator
 
Mattijs Janssens
Join Date: Mar 2009
Posts: 1,416
Rep Power: 15
mattijs is on a distinguished road
The wall distance is used in turbulence models to determine the effect of the wall so the current behaviour seems to be correct - it is the shortest path through fluid to the wall - it does not 'jump' gaps. If you want to change the behaviour you'll have to change finiteVolume/fvMesh/wallDist to do the calculation yourself.
mattijs is offline   Reply With Quote

Old   August 20, 2010, 17:45
Default
  #5
Member
 
Elisabet Mas de les Valls
Join Date: Mar 2009
Location: Barcelona, Spain
Posts: 60
Rep Power: 7
elisabet is on a distinguished road
Yes, I first thought this, but when I tried to define the inlet as the 'wall' boundary and find the wall distance, wallDist function jumped the gap.

See for instance the U-shaped channel defined by blockMeshDict.txt. When the bottom is defined as the unique 'wall', the output of wallDist is: bottom_wall.png

However, when the inlet (left and bottom boundary) is the unique 'wall', the output is:
inlet_wall.png

hence, in this case, it jumps the gap!
elisabet is offline   Reply With Quote

Old   August 23, 2010, 17:14
Default
  #6
Super Moderator
 
Mattijs Janssens
Join Date: Mar 2009
Posts: 1,416
Rep Power: 15
mattijs is on a distinguished road
Hi Elisabet,

you are indeed right - it does jump the gap. The information travels on the mesh but then the distance is calculated 'as the crow flies'.

Attached an interesting hack in the form of replacements for wallDist.C, patchWave.[CH] and wallPoint*[CH]. It calculates the cumulative distance through the mesh which I think is what you want for your case. It also reads the wallDistance if it is already there so you can calculate your own as a preprocessing step.

Mattijs
Attached Files
File Type: gz wallDist.tar.gz (5.8 KB, 6 views)
mattijs is offline   Reply With Quote

Old   August 27, 2010, 17:48
Default
  #7
Member
 
Elisabet Mas de les Valls
Join Date: Mar 2009
Location: Barcelona, Spain
Posts: 60
Rep Power: 7
elisabet is on a distinguished road
Dear Mattijs,

Sorry for the delay, but I've been really busy these days....

I've saved your new files in the src directory and run Allwmake. There are some problems related with the usage of 'distSqr' in:

1. wallPointDataI.H: there are no arguments to ‘distSqr’ that depend on a template parameter, so a declaration of ‘distSqr’ must be available

2. wallPointYPlusI.H: ‘distSqr’ was not declared in this scope

3. inverseFaceDistanceDiffusivity.C: ‘const class Foam::wallPoint’ has no member named ‘distSqr’

4. patchDataWave.C: ‘const class Foam::wallPointData<Foam::Vector<double> >’ has no member named ‘distSqr’

5. wallPointDataI.H: ‘distSqr’ was not declared in this scope

Any suggestions?

Elisabet
elisabet is offline   Reply With Quote

Old   August 30, 2010, 12:59
Default
  #8
Super Moderator
 
Mattijs Janssens
Join Date: Mar 2009
Posts: 1,416
Rep Power: 15
mattijs is on a distinguished road
My fault - tested it just with a copy of a solver - they will not be drop-in replacements. Just copy your favourite solver and add the files (and add wallPoint.C, patchWave.C, wallDist.C to the Make/files).
mattijs is offline   Reply With Quote

Old   September 7, 2010, 11:56
Default
  #9
Member
 
Elisabet Mas de les Valls
Join Date: Mar 2009
Location: Barcelona, Spain
Posts: 60
Rep Power: 7
elisabet is on a distinguished road
Hi Mattijs,

I've arrived today from real holidays (those without any available PC), so probably I'm not reaaly concentrated today and I missed something.

I've saved all your files (.C and .H) in my application folder. In order to make it compile, I've changed the line 74 in wallDist.C:

old:
Code:
labelHashSet wallPatchIDs(getPatchIDs<wallPolyPatch>());
new:
Code:
labelHashSet wallPatchIDs(getPatchIDs(wallPolyPatch::typeName));
but it still provides the wrong result when I call the function
Code:
wallDist(mesh).y()
what am I missing?

elisabet
elisabet is offline   Reply With Quote

Old   September 7, 2010, 12:20
Default
  #10
Super Moderator
 
Mattijs Janssens
Join Date: Mar 2009
Posts: 1,416
Rep Power: 15
mattijs is on a distinguished road
Attached my little test program (a stripped down yPlusRas) which just calculates and writes wallDist. It compiles fine under 1.7.x, probably needs some small changes for 1.5.x.
Attached Files
File Type: gz dumpWallDist.tar.gz (6.8 KB, 5 views)
mattijs is offline   Reply With Quote

Old   September 8, 2010, 05:13
Default
  #11
Member
 
Elisabet Mas de les Valls
Join Date: Mar 2009
Location: Barcelona, Spain
Posts: 60
Rep Power: 7
elisabet is on a distinguished road
Ok, again the same results but I've realized that, based on the geometry I attached some days before:

1. If 'inlet' is the unique patch: your dumpWallDist works perfectly

2. If the 'refWall' is the unique patch: keeps giving the same result as before

so.... almost you get it!

elisabet
elisabet is offline   Reply With Quote

Reply

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Help needed: building 1.5.x on Lucid sripplinger OpenFOAM Installation 3 June 10, 2010 23:01
calculating wall distance wallDist volker OpenFOAM 8 May 12, 2010 15:29
Differences between 1.5.x and 1.5-dev carsten OpenFOAM 20 October 12, 2009 10:18
OpenFOAM 1.5.x package - CentOS 5.3 x86_64 linnemann OpenFOAM Installation 7 July 30, 2009 04:14
OpenFoam 1.5.x installation gareth__it_power OpenFOAM Installation 5 April 1, 2009 04:27


All times are GMT -4. The time now is 23:47.