Particle position each timestep (during simulation)

lindner · August 27, 2010, 07:32

Hi, i have the following problem:
I need to know how many particles are there in the zone 1 and zone 2 of a room.

An example:
10 meter long room (x direction), while the first 5 meters are the zone 1 and the other 5 meters are zone 2.
The particle position file outputs the coordinate (x y z) of each particle, so it will be like:
if x < 5, sum to zone 1,
else, sum to zone 2

The actual output while solving the uncoupledKinematicParcelFoam is

Cloud: kinematicCloud
Total number of parcels added = 1210
Total mass introduced = 1.01368723e-17
Current number of parcels = 1210
Current mass in system = 1.01368723e-17
ExecutionTime = 0.61 s ClockTime = 1 s

and my objective is to add the following
Current number of parcels in Zone 1 = xxx
Current number of parcels in Zone 2 = yyy

Is it possible to make it in the solver?
Or it's better to make a separate code to extract this information?
By the way, i don't know c++ programming =(

Thank you, and if anyone needs I can send the room from Annex 20 with the uncoupledKinematicParcelFoam working.

PS: Using OpenFOAM 1.7.x

Edison_Ge · September 10, 2010, 01:58

rather than modify the code, you can acturally find the particle position under time/lagrangian/postion. sort them in the dirction of interest in Excel or matlab. Then your answer would be obvious.

lindner · September 10, 2010, 04:22

Thank you for the answer Edison, this really works for simple cases if there is not many zones.
But the problem is to have more 3D zones where the coordinates are complicated.
I did a simple fortran code and its working by the moment.

Thank you for the tips

August 27, 2010, 07:32	Particle position each timestep (during simulation)	#1
lindner Member Join Date: Mar 2010 Posts: 42 Rep Power: 16	Hi, i have the following problem: I need to know how many particles are there in the zone 1 and zone 2 of a room. An example: 10 meter long room (x direction), while the first 5 meters are the zone 1 and the other 5 meters are zone 2. The particle position file outputs the coordinate (x y z) of each particle, so it will be like: if x < 5, sum to zone 1, else, sum to zone 2 The actual output while solving the uncoupledKinematicParcelFoam is Cloud: kinematicCloud Total number of parcels added = 1210 Total mass introduced = 1.01368723e-17 Current number of parcels = 1210 Current mass in system = 1.01368723e-17 ExecutionTime = 0.61 s ClockTime = 1 s and my objective is to add the following Current number of parcels in Zone 1 = xxx Current number of parcels in Zone 2 = yyy Is it possible to make it in the solver? Or it's better to make a separate code to extract this information? By the way, i don't know c++ programming =( Thank you, and if anyone needs I can send the room from Annex 20 with the uncoupledKinematicParcelFoam working. PS: Using OpenFOAM 1.7.x

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Particle Trajectory calculation	yueroo	FLUENT	2	July 13, 2009 14:32
Display Particle Position on a Plane	Jana	FLUENT	2	July 13, 2007 03:46
Timestep for transient simulation	houman	CFX	2	April 28, 2005 04:44
timestep for combustion simulation	yangqing	FLUENT	0	August 18, 2002 09:46
Particle Trajectory calculation when particle hit	yueroo	FLUENT	2	February 28, 2001 05:26

September 10, 2010, 01:58		#2
Edison_Ge Member edison Join Date: May 2009 Location: Australia Posts: 35 Rep Power: 17	rather than modify the code, you can acturally find the particle position under time/lagrangian/postion. sort them in the dirction of interest in Excel or matlab. Then your answer would be obvious.

September 10, 2010, 04:22		#3
lindner Member Join Date: Mar 2010 Posts: 42 Rep Power: 16	Thank you for the answer Edison, this really works for simple cases if there is not many zones. But the problem is to have more 3D zones where the coordinates are complicated. I did a simple fortran code and its working by the moment. Thank you for the tips