Creating a case for uncoupledKinematicParcelFoam (OpenFOAM 1.6)

andrewryan · November 23, 2009, 12:58

Dear Foamers

I'm trying to understand how the uncoupledKinematicParcelFoam solver works. Unfortunately I could not find a tutorial case for it. Therefore I'm trying to create my own case.

From the description:

Code:

Transient solver for the passive transport of a single kinematic particle could.

Uses a pre- calculated velocity field to evolve the cloud.

So far I figured out that it's looking for the initial positions in 0/lagrangian/kinematicCloud/positions .

The file has the format:

Code:

n
(
( x y z) cellID
..
)

Where n is the number of particles and each row contains the position of the individual particle and the cellID in which it is. (Is this correct?)

This is exactly the problem: How can I figure out in which cell the particle is, given the position. I saw that in some other tutorial (for example the evaporationTest) cases this file is written in all timestep directories (except 0) by the solver.

Whats the way to generate the initial position file? (for now I just want to have a few particles somewhere in my simple mesh)

How can I figure out the id of a cell in paraview? When using foamToVTK I could select to colorize the wireframe by cellIDs [1], this however only works when the mesh is really small. Also the cell IDs did not show up the parallel openfoam reader. Is there another way to see the cell id (for example by clicking on a cell)?

[1] http://emerging-patterns.com/tmp/CellIDs.png

maksen · March 7, 2010, 04:29

Hi Andrew,

did you succeed with the uncoupledKinematicParcelFoam? I was looking for a tutorial but found your post...

could you make your case work?

cheers
markus

andrewryan · March 7, 2010, 04:35

No so far I just used icoLagrangianFoam..

maksen · March 8, 2010, 02:36

alright..
i played with the rhoPisoTwinParcelFoam yesterday which does particle tracking, however i haven't found out how to visualize them properly ...

cheers
markus

andrewryan · March 8, 2010, 02:46

> however i haven't found out how to visualize them properly ...

http://openfoamwiki.net/index.php/Ma...gian_particles

foamToVTK seems to be the only option here..

maksen · March 8, 2010, 04:42

yes, i agree... i "glyphed" the p-tracks yesterday and it seemed to work...

how are your experiences with icoLagrangianFoam. found a tutorial on the net which says that its "only a demo"....

later
markus

andrewryan · March 8, 2010, 12:08

Well I only looked at it to see how LPT works..

This here could be interesting for you: http://www.tfd.chalmers.se/~hani/kur...angianFoam.pdf

maksen · March 8, 2010, 12:52

thats the tutorial i found as well...thanks

cheers
markus

randolph · July 19, 2017, 11:09

Quote:

Originally Posted by andrewryan

No so far I just used icoLagrangianFoam..

isn't icoLagrangianFoam quite limited? why not use icoUncoupledKinematicParcelFoam?

November 23, 2009, 12:58	Creating a case for uncoupledKinematicParcelFoam (OpenFOAM 1.6)	#1
andrewryan Member Andrew Ryan Join Date: Mar 2009 Posts: 47 Rep Power: 17	Dear Foamers I'm trying to understand how the uncoupledKinematicParcelFoam solver works. Unfortunately I could not find a tutorial case for it. Therefore I'm trying to create my own case. From the description: Code: Transient solver for the passive transport of a single kinematic particle could. Uses a pre- calculated velocity field to evolve the cloud. So far I figured out that it's looking for the initial positions in 0/lagrangian/kinematicCloud/positions . The file has the format: Code: n ( ( x y z) cellID .. ) Where n is the number of particles and each row contains the position of the individual particle and the cellID in which it is. (Is this correct?) This is exactly the problem: How can I figure out in which cell the particle is, given the position. I saw that in some other tutorial (for example the evaporationTest) cases this file is written in all timestep directories (except 0) by the solver. Whats the way to generate the initial position file? (for now I just want to have a few particles somewhere in my simple mesh) How can I figure out the id of a cell in paraview? When using foamToVTK I could select to colorize the wireframe by cellIDs [1], this however only works when the mesh is really small. Also the cell IDs did not show up the parallel openfoam reader. Is there another way to see the cell id (for example by clicking on a cell)? [1] http://emerging-patterns.com/tmp/CellIDs.png maysmech likes this.

March 7, 2010, 04:29	uncoupledKinematicParcelFoam	#2
maksen New Member Markus Trenker Join Date: Feb 2010 Location: Vienna, Austria Posts: 11 Rep Power: 16	Hi Andrew, did you succeed with the uncoupledKinematicParcelFoam? I was looking for a tutorial but found your post... could you make your case work? cheers markus

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
OpenFOAM 1.6 GGI and mergeMeshes	hansel	OpenFOAM Bugs	4	April 14, 2010 05:16
OpenFoam 1.6 Permissions in Ubuntu 9.04	AlanR	OpenFOAM Installation	14	December 28, 2009 19:45
paraFoam reader for OpenFOAM 1.6	smart	OpenFOAM Installation	13	November 16, 2009 21:41
FoamX for OpenFOAM 1.6	vitor	OpenFOAM	4	October 29, 2009 06:24

March 7, 2010, 04:35		#3
andrewryan Member Andrew Ryan Join Date: Mar 2009 Posts: 47 Rep Power: 17	No so far I just used icoLagrangianFoam..

March 8, 2010, 02:36		#4
maksen New Member Markus Trenker Join Date: Feb 2010 Location: Vienna, Austria Posts: 11 Rep Power: 16	alright.. i played with the rhoPisoTwinParcelFoam yesterday which does particle tracking, however i haven't found out how to visualize them properly ... cheers markus

March 8, 2010, 02:46		#5
andrewryan Member Andrew Ryan Join Date: Mar 2009 Posts: 47 Rep Power: 17	> however i haven't found out how to visualize them properly ... http://openfoamwiki.net/index.php/Ma...gian_particles foamToVTK seems to be the only option here..

March 8, 2010, 04:42		#6
maksen New Member Markus Trenker Join Date: Feb 2010 Location: Vienna, Austria Posts: 11 Rep Power: 16	yes, i agree... i "glyphed" the p-tracks yesterday and it seemed to work... how are your experiences with icoLagrangianFoam. found a tutorial on the net which says that its "only a demo".... later markus

March 8, 2010, 12:08		#7
andrewryan Member Andrew Ryan Join Date: Mar 2009 Posts: 47 Rep Power: 17	Well I only looked at it to see how LPT works.. This here could be interesting for you: http://www.tfd.chalmers.se/~hani/kur...angianFoam.pdf

March 8, 2010, 12:52		#8
maksen New Member Markus Trenker Join Date: Feb 2010 Location: Vienna, Austria Posts: 11 Rep Power: 16	thats the tutorial i found as well...thanks cheers markus