CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Post-Processing (https://www.cfd-online.com/Forums/openfoam-post-processing/)
-   -   sampling when mesh has moving points (https://www.cfd-online.com/Forums/openfoam-post-processing/114025-sampling-when-mesh-has-moving-points.html)

Phicau March 4, 2013 02:50

sampling when mesh has moving points
 
Dear all,

I am trying to sample along a line for a mesh where points move, hence, when the initial cells where the line is located move away I get an error and the sampling process is aborted:

Code:

Time = 0

Time = 0.05

Time = 0.1

Time = 0.15


--> FOAM FATAL ERROR:
None of the neighbours of face 887211 contains point (44.9001 0.000100006 1.105)

    From function sampledSet::getCell(const label, const point&)
    in file sampledSet/sampledSet/sampledSet.C at line 105.

FOAM aborting

Any clue on how to tell OpenFOAM to read the "points" file from the time folder instead from the polymesh one?

Just to clarify, I want to measure always at the same location, it does not matter at which cell.

Best,

Pablo

ngj March 4, 2013 03:34

Hi Pablo,

What version are you using? I have had similar problems in the official release, but never in 1.6-ext. I never came across a solution, though.

Good luck,

Niels

Phicau March 4, 2013 03:41

Hi Niels,

we are currently using 2.1.1.

Good to know, I will then try to dig my way out of here programming...

Thanks!

fredo490 March 4, 2013 10:42

Does this problem happen when only a part of your line I'd out of the domain or only when all the line gets out?
By the way, what can of deformation do you use? Translating?

Phicau March 4, 2013 11:19

Hi,

actually my line is always within the domain. The problem arises when any location of the line starts to belong to other cell than the initial one.

I am using a modified version of dynamicMotionSolverFvMesh.

Best,

Pablo

fredo490 March 4, 2013 11:24

Can you use a post treatment software like EnSight or ParaView? It would make the thing easier (but longer)

Phicau March 4, 2013 11:36

Definitely.

There is another option which I currently use: sample for free surface as an isosurface, and then by using a python script I obtain the time series of free surface elevation anywhere. This is kind of brute force, though.


All times are GMT -4. The time now is 10:46.