|
[Sponsors] |
September 11, 2012, 13:51 |
mapFields problem with 3D
|
#1 |
Senior Member
Robert Castilla
Join Date: Apr 2009
Location: Spain
Posts: 109
Rep Power: 17 |
Hi,
when I try to map fields from a case to another with a 3D geometry made with SHM, I found a lot of warnings like these: --> FOAM Warning : From function Foam::List<Foam::tetIndices> Foam:olyMeshTetDecomposition::faceTetIndices(con st polyMesh&, label, label) in file meshes/polyMesh/polyMeshTetDecomposition/polyMeshTetDecomposition.C at line 565 No base point for face 514084, 4(350741 350091 350739 348992), produces a valid tet decomposition. May be they are not important, but after the mapping the Courant number begins to growth and finally it totally diverges. The solver is the pimpleFoam. Do you think that this warning with the mapping can be related to the divergence of the solution? Or it is unrelevant and the cause of the divergence is in another place? I have to add that when I did it with a similar (not exactly the same) case with 2D, I got no warning with the mapping, and the simulation ran fine. Thank you Robert |
|
September 11, 2012, 18:14 |
|
#2 |
Senior Member
Marco A. Turcios
Join Date: Mar 2009
Location: Vancouver, BC, Canada
Posts: 740
Rep Power: 28 |
Are you mapping consistently or with mapped/cutting patches? And more importantly, what do your mapped fields look like?
|
|
September 12, 2012, 11:43 |
|
#3 |
Senior Member
Robert Castilla
Join Date: Apr 2009
Location: Spain
Posts: 109
Rep Power: 17 |
They are consistent (option -consistent in mapFields)
I don't understand completely what do you mean with what the fields look like... The flow is 3D and it is difficult to show velocity and pressure in a graphic. But, anyway, I have the impression that this warning is more related to the mesh than to physical fields. Why do you think that this is important? I would like also to know how to mark and show these problematic faces in paraview. Regards Robert |
|
September 12, 2012, 11:59 |
|
#4 |
Senior Member
Marco A. Turcios
Join Date: Mar 2009
Location: Vancouver, BC, Canada
Posts: 740
Rep Power: 28 |
Hi Robert,
The reason I ask about the fields is that when mapping to a snappyHexMesh case, since the boundaries will not be exactly the same you probably won't be able to successfully map all cells, particularly if the target domain is slightly larger than the source domain (look at the user guide entry for mapping fields to see what I mean; basically any cells in your target domain that are outside the source domain geometry don't have a meaningful value given to them). To fix something like this, have a look at the following thread: http://www.cfd-online.com/Forums/ope...mapfields.html Good luck! |
|
September 13, 2012, 10:37 |
|
#5 |
Senior Member
Robert Castilla
Join Date: Apr 2009
Location: Spain
Posts: 109
Rep Power: 17 |
Hi, Marco,
finally I got the point. I was mapping to a dynamic mesh which was deformed. I didn't realize that in the deformation process some faces became wrong oriented, and that was the cause of the warning. Anyway, it seems that it is not the cause of the divergence. Thank you very much for your help! Robert |
|
October 17, 2012, 02:56 |
|
#6 |
Senior Member
Robert Castilla
Join Date: Apr 2009
Location: Spain
Posts: 109
Rep Power: 17 |
Related to this same issue, I have a doubt.
If I have a dynamic mesh, the points coordinates are changing in time, and they are stored in the file polyMesh/points in each time directory. When I interpolate FROM a time directory in a case, does the mapFields take the points coordinates from the time directory or from the constant/polyMesh directory? I have the impression that it should take it from the time directory, but I have doubts from the results that I obtain. Regards Robert |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
UDF compiling problem | Wouter | Fluent UDF and Scheme Programming | 6 | June 6, 2012 04:43 |
Problem in implementing cht | tilek | CFX | 3 | May 8, 2011 08:39 |
natural convection problem for a CHT problem | Se-Hee | CFX | 2 | June 10, 2007 06:29 |
Adiabatic and Rotating wall (Convection problem) | ParodDav | CFX | 5 | April 29, 2007 19:13 |
Is this problem well posed? | Thomas P. Abraham | Main CFD Forum | 5 | September 8, 1999 14:52 |