|
[Sponsors] |
September 29, 2005, 05:47 |
Has anyone done run-time probi
|
#1 |
Senior Member
Håkan Nilsson
Join Date: Mar 2009
Location: Gothenburg, Sweden
Posts: 204
Rep Power: 18 |
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. |
|
September 29, 2005, 13:14 |
I would do something like:
-
|
#2 |
Senior Member
Mattijs Janssens
Join Date: Mar 2009
Posts: 1,419
Rep Power: 26 |
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. |
|
October 18, 2005, 11:42 |
Hi
How can I get the label f
|
#3 |
Member
olivier Petit
Join Date: Mar 2009
Location: Göteborg, Sweden
Posts: 67
Rep Power: 17 |
Hi
How can I get the label faces for volume cell ? thank you for helping |
|
October 18, 2005, 11:51 |
I think that faceCells (). it'
|
#4 |
Member
olivier Petit
Join Date: Mar 2009
Location: Göteborg, Sweden
Posts: 67
Rep Power: 17 |
I think that faceCells (). it's right?
thank u :-) |
|
October 18, 2005, 12:27 |
faceCells() gives you cells ar
|
#5 |
Senior Member
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,905
Rep Power: 33 |
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 |
|
October 18, 2005, 15:16 |
thank you Jasak,
Can you give
|
#6 |
Member
olivier Petit
Join Date: Mar 2009
Location: Göteborg, Sweden
Posts: 67
Rep Power: 17 |
thank you Jasak,
Can you give me the expression to use cellFaces()! |
|
October 19, 2005, 04:08 |
A 'cell' is a list of the face
|
#7 |
Senior Member
Mattijs Janssens
Join Date: Mar 2009
Posts: 1,419
Rep Power: 26 |
A 'cell' is a list of the face labels (section 6.1.1.3 in user guide)
const cell& cFaces = mesh.cells()[cellI]; |
|
July 7, 2009, 07:22 |
|
#8 | |
Senior Member
Sebastian Gatzka
Join Date: Mar 2009
Location: Frankfurt, Germany
Posts: 729
Rep Power: 20 |
Quote:
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!" |
||
July 7, 2009, 16:41 |
|
#9 |
Senior Member
Mattijs Janssens
Join Date: Mar 2009
Posts: 1,419
Rep Power: 26 |
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; |
|
July 7, 2009, 17:32 |
|
#10 |
Senior Member
Sebastian Gatzka
Join Date: Mar 2009
Location: Frankfurt, Germany
Posts: 729
Rep Power: 20 |
Sounds promising to me.
Furthermore, can you tell me the two following things:
__________________
Schrödingers wife: "What did you do to the cat? It's half dead!" |
|
|
|