[holzmichel]

Hello all,

can someone explain me what the number after the lagrangian position (x y z) n means, please. I have tried to find out the meaning but I come to no solution .

Thanks

Michél

[niklas]

It is the normal of the plane that the parcel moves in for 2D wedge simlations.

[holzmichel]

thank you for the answer

[niklas]

actually, i was wrong.

its the cell number of the cell it is currently in.
the vectors are all written last.

[holzmichel]

well. are you sure? because i played with the number and the particles start every time at the same position but has different movements properties.

[holzmichel]

ok. you are right.
but where can i found the function for this?

[niklas]

Originally Posted by holzmichel

ok. you are right.
but where can i found the function for this?

src/lagrangian/basic/Particle/Particle.H

[holzmichel]

thank you again

[wikstrom]

Hi,

I want to trace passive particles through a steady state RANS-solution, and have written a simple code using kinemaitcCloud. Seems to work well at first, but as the time steps march on, the clock time for each time step increases immensely. What might be the reason to this behaviour? Sure, as the particles enter further into the solution the flow speed increases, but not tremendously so. See the attached graph showing how clock time increases with time steps. Number of parcels are ~1e4 and the case size is about 2M cells from snappyHexMesh.

Cheers
Niklas

[gschaider]

Hi Niklas!

Originally Posted by wikstrom

I want to trace passive particles through a steady state RANS-solution, and have written a simple code using kinemaitcCloud. Seems to work well at first, but as the time steps march on, the clock time for each time step increases immensely. What might be the reason to this behaviour? Sure, as the particles enter further into the solution the flow speed increases, but not tremendously so. See the attached graph showing how clock time increases with time steps. Number of parcels are ~1e4 and the case size is about 2M cells from snappyHexMesh.

First question: does the number of particles stay constant during the simulation? (I guess you already checked for that). If for instance injected particles get trapped in some cavity and won't leave the system that would explain the increase

Does the time spent in the lagrangian solver increase or in the other parts? To find an answer to things like that I use this extension:
http://www.cfd-online.com/Forums/ope...profiling.html
If you prefer to use valgrind (like Henry suggested in that thread) I'd suggest that you restart from one of the later timesteps (where the calculation-time already exploded) because callgrind multiplies the time a run needs(at least for me, but maybe I'm using it wrong).

If the number of particles in the system doesn't increase and the lagrangian part is the one that eats up the CPU-time then I guess the reason is that the time-step during the lagrangian calculation becomes very small (either because particles have to travel through very small cells or because particles often pass zones were the velocity relative to the gas is high) but that is just a guess

Bernhard

[wikstrom]

Thank you Bernhard for the answer!

The answer probably lies in small local time steps as the particles enter the 'interesting region'. To be realistic, sadly, I will not have the time to profile the solver; It was years since I used Valgrind I'm afraid...

again, Thank you.

N.

[wikstrom]

Oh, it does not look that hard to use :-) I'll try it.

N.