|
[Sponsors] |
September 14, 2012, 13:05 |
Problems with surface sampling
|
#1 |
New Member
Matteo Carpentieri
Join Date: Mar 2009
Posts: 28
Rep Power: 17 |
Hi all,
I'm trying to sample the pressure values on a wall surface. I can easily see the values in paraView, but when I execute the sample utility I don't get any output files. The folder "surface/10000/" is created but it's empty. I'm using OpenFOAM 2.1.x and this is my sampleDict file: Code:
surfaceFormat raw; interpolationScheme cellPoint; fields ( p ); surfaces ( wallPressure { type patch; patches (buildingWall); interpolate true; triangulate false; } ); Using the sample utility with sets, instead of surface, on a set of sample points, works well. Can anyone help me? I'm sure I'm doing a very silly mistake, but I can't see it... Thank you |
|
September 17, 2012, 03:52 |
|
#2 | |
Member
|
Quote:
Then surface is probably called something different. You can check by using "patchSummary". This will list all boundary conditions including the full names of all patches/walls. Could be something like "buildingWall_patch0" perhaps ...
__________________
~~~_/)~~~ |
||
September 17, 2012, 07:53 |
|
#3 | |
New Member
Matteo Carpentieri
Join Date: Mar 2009
Posts: 28
Rep Power: 17 |
Quote:
However this is not the case. I've tried several different patches and always got the same result. This is the output from patchSummary, for example: Code:
[...] wall: buildingWall scalar nut generic scalar k generic scalar p zeroGradient scalar epsilon generic vector U fixedValue symmTensor R generic [...] |
||
September 17, 2012, 17:46 |
|
#4 |
Member
|
hmm ... maybe try this:
Code:
setFormat raw; surfaceFormat vtk; interpolationScheme cellPoint; fields ( p ); surfaces ( wallPressure { type patch; patches (buildingWall); interpolate true; triangulate false; } );
__________________
~~~_/)~~~ |
|
September 18, 2012, 06:56 |
|
#5 |
New Member
Matteo Carpentieri
Join Date: Mar 2009
Posts: 28
Rep Power: 17 |
Thanks again SirWombat.
Tried that, and several other combinations, but unfortunately no results at all... |
|
September 18, 2012, 08:08 |
|
#6 |
Member
|
Hi Matteo, just too bad!
You may wanna try a different approach: use a basePoint and normalVector instead, i.e. like this: Code:
setFormat raw; surfaceFormat vtk; interpolationScheme cellPoint; fields ( p ); surfaces ( wallplane{ type plane; basePoint (0 0.5 0); normalVector (0 1 0); interpolate true; } ); Hope this gets you somewhere! Greetings, Jan
__________________
~~~_/)~~~ |
|
September 19, 2012, 05:18 |
|
#7 | |
New Member
Matteo Carpentieri
Join Date: Mar 2009
Posts: 28
Rep Power: 17 |
Quote:
no luck. It seems that surface sampling is not working at all, regardless the surface chosen... I'm going to try with runtime function objects, but, of course, I have to re-run the simulation. Could it this be a bug? Is anyone else having problems with surface sampling? Thanks again Matteo |
||
September 19, 2012, 05:27 |
|
#8 |
Member
|
Hi Matteo, could you post the output of the sample-utility? Have you tried other fields (i.e. U)? What solver did you use?
As you said you can see the data in paraview, i guess your setup and output should be ok. the sample utility really should work then. I never had problems using 'sample', but i have only used it with simpleFoam (1 phase, steady), LTSInterFoam (2 phase, steady), interFoam (2phase, transient). Hope someone else has an idea here. G Jan
__________________
~~~_/)~~~ |
|
September 19, 2012, 06:00 |
|
#9 |
New Member
Matteo Carpentieri
Join Date: Mar 2009
Posts: 28
Rep Power: 17 |
Here is the output of the sample command:
Code:
/*---------------------------------------------------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 2.1.x | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ Build : 2.1.x-c1ffc2074f88 Exec : sample -latestTime Date : Sep 19 2012 Time : 09:56:17 Host : "EnFlo-MC01" PID : 11479 Case : /home/mcarpentieri/OpenFOAM/mcarpentieri-2.1.x/run/HRModUrb/simpleC_0 nProcs : 1 sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE). fileModificationChecking : Monitoring run-time modified files using timeStampMaster allowSystemOperations : Disallowing user-supplied system call operations // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // Create time Create mesh for time = 1811 Reading surface description: wallPressure Time = 1811 End Thanks Matteo |
|
September 19, 2012, 06:05 |
|
#10 |
Member
|
hmmm ... if its not the surface description i am really lost on this one ... no idea whats going wrong. sorry! G jan
__________________
~~~_/)~~~ |
|
September 19, 2012, 09:13 |
|
#11 | |
New Member
Matteo Carpentieri
Join Date: Mar 2009
Posts: 28
Rep Power: 17 |
Quote:
By the way, the function object sampledPatch is working perfectly. I just have to re-run all the simulations, now. |
||
November 9, 2012, 12:00 |
|
#12 |
Member
Albert Tong
Join Date: Dec 2010
Location: Perth, WA, Australia
Posts: 76
Blog Entries: 1
Rep Power: 16 |
Hi Matteo,What's sampledPatch look like? Can you please give an example? I searched on-line and the tutorials of 0F2.1.1 but cannot find. Thanks.
__________________
Kind regards, Albert |
|
Tags |
patch, sample surface, sampledict |
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
[snappyHexMesh] snappyHexMesh aborting | Tobi | OpenFOAM Meshing & Mesh Conversion | 0 | November 10, 2010 04:23 |
[snappyHexMesh] Boundary layer generation problems | ivan_cozza | OpenFOAM Meshing & Mesh Conversion | 0 | October 6, 2010 14:47 |
time step in free surface problems | KtoTo | Siemens | 4 | June 26, 2007 08:03 |
icem surface problems | joe | Main CFD Forum | 0 | February 11, 2004 18:20 |
free convection heat transfer from a heated horizontal surface through a liquid to a thin cooled fin | Kaushik | FLUENT | 1 | May 8, 2000 07:47 |