CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

Using faceSets for probing

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

Reply
 
LinkBack Thread Tools Display Modes
Old   September 29, 2005, 05:47
Default Has anyone done run-time probi
  #1
Senior Member
 
Håkan Nilsson
Join Date: Mar 2009
Location: Gothenburg, Sweden
Posts: 185
Rep Power: 7
hani is on a distinguished road
Has anyone done run-time probing using faceSets? (c.f. probing in oodles) Or can someone give me a hint on how to do it?

I want to use faceSets to output the results at those 'internal patches' during the computation. The results on the faceSets should be output similar to how the external patches are output.

The main reasons why I want to do this are that I am mostly interested in the results at particular cuts of the domain, and that the disk won't allow the whole field to be written for all the needed timesteps. It would also be a nice tool to see how the solution propagates while making the computations. It would also save lots of time when post-processing in for instance Ensight, where I would then be able to read only the most interesting parts of the domain.

I got the following proposition By Mattijs Janssens on Thursday, September 22, 2005 - 10:35 am:

- you can make them (faceSets) into external patches using the splitMesh utility (discussed before)

- create a cellSet from the cells on one/both sides of the faceSet using the cellSet utility. Visualize the results on the cells with foamToVTK

As I understand it, these propositions require that the whole filed is saved each timestep, which is impossible for large computations that run for many timesteps.

Håkan.
hani is offline   Reply With Quote

Old   September 29, 2005, 13:14
Default I would do something like: -
  #2
Super Moderator
 
Mattijs Janssens
Join Date: Mar 2009
Posts: 1,416
Rep Power: 15
mattijs is on a distinguished road
I would do something like:
- get the sorted face labels from the faceSet
- get the values on the faces and dump them like probes, so preserving time directory.
- extend e.g. foamToVTK to read and convert these.
mattijs is offline   Reply With Quote

Old   October 18, 2005, 11:42
Default Hi How can I get the label f
  #3
Member
 
olivier Petit
Join Date: Mar 2009
Location: Göteborg, Sweden
Posts: 67
Rep Power: 7
olivier is on a distinguished road
Hi
How can I get the label faces for volume cell ?
thank you for helping
olivier is offline   Reply With Quote

Old   October 18, 2005, 11:51
Default I think that faceCells (). it'
  #4
Member
 
olivier Petit
Join Date: Mar 2009
Location: Göteborg, Sweden
Posts: 67
Rep Power: 7
olivier is on a distinguished road
I think that faceCells (). it's right?
thank u :-)
olivier is offline   Reply With Quote

Old   October 18, 2005, 12:27
Default faceCells() gives you cells ar
  #5
Senior Member
 
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,741
Rep Power: 20
hjasak will become famous soon enough
faceCells() gives you cells around the face. cellFaces() gives you faces of the cell.

Enjoy,

Hrv
__________________
Hrvoje Jasak
Providing commercial FOAM/OpenFOAM and CFD Consulting: http://wikki.co.uk
hjasak is offline   Reply With Quote

Old   October 18, 2005, 15:16
Default thank you Jasak, Can you give
  #6
Member
 
olivier Petit
Join Date: Mar 2009
Location: Göteborg, Sweden
Posts: 67
Rep Power: 7
olivier is on a distinguished road
thank you Jasak,
Can you give me the expression to use cellFaces()!
olivier is offline   Reply With Quote

Old   October 19, 2005, 04:08
Default A 'cell' is a list of the face
  #7
Super Moderator
 
Mattijs Janssens
Join Date: Mar 2009
Posts: 1,416
Rep Power: 15
mattijs is on a distinguished road
A 'cell' is a list of the face labels (section 6.1.1.3 in user guide)

const cell& cFaces = mesh.cells()[cellI];
mattijs is offline   Reply With Quote

Old   July 7, 2009, 07:22
Default
  #8
Senior Member
 
sega's Avatar
 
Sebastian Gatzka
Join Date: Mar 2009
Location: Frankfurt, Germany
Posts: 729
Rep Power: 10
sega is on a distinguished road
Quote:
Originally Posted by mattijs View Post
I would do something like:
- get the sorted face labels from the faceSet
- get the values on the faces and dump them like probes, so preserving time directory.
- extend e.g. foamToVTK to read and convert these.
Dear mattuijs.

Can you tell me how I can get the values at the faces?

I suppose I have to run faceSet first.
But what do I have to do with the list of face numbers to get the values at their face centers?!
__________________
Schrödingers wife: "What did you do to the cat? It's half dead!"
sega is offline   Reply With Quote

Old   July 7, 2009, 16:41
Default
  #9
Super Moderator
 
Mattijs Janssens
Join Date: Mar 2009
Posts: 1,416
Rep Power: 15
mattijs is on a distinguished road
Either you already have a field on faces (e.g. a surfaceScalarField) or you'll need to interpolate a volField (fvc::interpolate) to create one. Then you can just index the field using the face label (if it is an internal face). E.g.

surfaceScalarField s = fvc::interpolate(p);
Pout<< s[10] << endl;
mattijs is offline   Reply With Quote

Old   July 7, 2009, 17:32
Default
  #10
Senior Member
 
sega's Avatar
 
Sebastian Gatzka
Join Date: Mar 2009
Location: Frankfurt, Germany
Posts: 729
Rep Power: 10
sega is on a distinguished road
Sounds promising to me.
Furthermore, can you tell me the two following things:
  • How can I read the face index out of the faceSet output file?
  • Is the output of your command s[index] the value at the face center?
  • How can I write the face center value (and maybe its coordinates) into a file?
Good night.
__________________
Schrödingers wife: "What did you do to the cat? It's half dead!"
sega 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



All times are GMT -4. The time now is 05:31.