sampling when mesh has moving points

Phicau · March 4, 2013, 03:50

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, 04: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, 04: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, 11: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, 12: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, 12:24

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

Phicau · March 4, 2013, 12: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.

March 4, 2013, 03:50	sampling when mesh has moving points	#1
Phicau Senior Member Pablo Higuera Join Date: Jan 2011 Location: Auckland Posts: 627 Rep Power: 19	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

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
[snappyHexMesh] SnappyHexMesh for internal Flow	vishwa	OpenFOAM Meshing & Mesh Conversion	24	June 27, 2016 09:54
Question on moving mesh, mesh velocity is really small!	ripperjack	Main CFD Forum	2	April 28, 2014 14:37
3D Hybrid Mesh Errors	DarrenC	ANSYS Meshing & Geometry	11	August 5, 2013 07:42
Dynamic Mesh moving interface help	akash.iitb	FLUENT	0	August 24, 2010 00:53
[snappyHexMesh] external flow with snappyHexMesh	chelvistero	OpenFOAM Meshing & Mesh Conversion	11	January 15, 2010 20:43

March 4, 2013, 04:34		#2
ngj Senior Member Niels Gjoel Jacobsen Join Date: Mar 2009 Location: Copenhagen, Denmark Posts: 1,903 Rep Power: 37	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

March 4, 2013, 04:41		#3
Phicau Senior Member Pablo Higuera Join Date: Jan 2011 Location: Auckland Posts: 627 Rep Power: 19	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!

March 4, 2013, 11:42		#4
fredo490 Senior Member HECKMANN Frédéric Join Date: Jul 2010 Posts: 249 Rep Power: 17	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?

March 4, 2013, 12:19		#5
Phicau Senior Member Pablo Higuera Join Date: Jan 2011 Location: Auckland Posts: 627 Rep Power: 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

March 4, 2013, 12:24		#6
fredo490 Senior Member HECKMANN Frédéric Join Date: Jul 2010 Posts: 249 Rep Power: 17	Can you use a post treatment software like EnSight or ParaView? It would make the thing easier (but longer)

March 4, 2013, 12:36		#7
Phicau Senior Member Pablo Higuera Join Date: Jan 2011 Location: Auckland Posts: 627 Rep Power: 19	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.