[zizi63]

Hi Kevin,
Maybe like the coal example we should use two separate kinematic cloud properties, if it works for you will you please inform me?

[Scram_1]

Hello!
Has anyone managed to inject particles of different densities through two different patches?

My case is a Y - junction with 2 inlets and 1 outlet. I want to analyze the mixing pattern of the particles. For this, I want my particles to be colored differently, say white particles injected from the top inlet and black from the bottom inlet. If this were the case then I can easily see the mixing pattern. The problem is, since the particles are generated from the same kinematic cloud, they are all colored the same. I don't think it's possible to generate 2 kinematic clouds using icoUncoupledKinematicParcelFoam because in the .C file, the cloud is generated using kinematicCloud.evolve(). So if I have 'two' clouds, the second one will overwrite the first one and effectively there will only be one cloud. I was thinking if particles of different densities are injected, then I can color them differently which will serve my purpose. The easiest way would be to find a way to color the particles based on the patch from which they are injected in ParaView itself. But I couldn't find an option to that effect.

Also, I tried using the coalCloud feature in coalChemistryFoam but it requires a lot of auxiliary information like radiation models, chemistry properties etc. and it's not allowing me to turn off radiation.

Best,
Ram

[a.weber]

Hi,
I'm also planning to do two clouds of totally different properties.... e.g. gas bubbles and catalysator particles.
Differences in calculation are important because:
bubbles should collide and coalesce with each other, cat particles should not.

What i did so far is to just add a second cloud in the createClouds.C:

Info<< "\nConstructing bubble cloud" << endl;
basicSprayCloud bubbles
(
    "bubbles",
    rho,
    U,
    g,
    slgThermo
);

Info<< "\nConstructing catalysator cloud" << endl;
basicSprayCloud cat
(
    "cat",
    rho,
    U,
    g,
    slgThermo
);

Then I just solve them both in my solver.C:

Info<< "Time = " << runTime.timeName() << nl << endl;
        bubbles.evolve();
        cat.evolve();
...

This works indeed, I get Info about both Clouds in the Output:

Solving 3-D cloud bubbles

--> Cloud: bubbles injector: model1
    Added 2 new parcels

Cloud: bubbles
    Current number of parcels       = 3
    Current mass in system          = 8.55318e-08
    Linear momentum                 = (7.09132e-11 1.15702e-11 2.83297e-08)
   |Linear momentum|                = 2.83298e-08
    Linear kinetic energy           = 4.69921e-09
    model1:
        number of parcels added     = 3
        mass introduced             = 8.55318e-08
    Parcels absorbed into film      = 0
    New film detached parcels       = 0
    Parcel fate (number, mass)
      - escape                      = 0, 0
      - stick                       = 0, 0
    Temperature min/max             = 300, 300
    Mass transfer phase change      = 0
    D10, D32, Dmax (mu)             = 3555.08, 3555.08, 3555.11
    Liquid penetration 95% mass (m) = 0.00136594
    hold-up of gasphase             = 0.000230061 %


Solving 3-D cloud cat

--> Cloud: cat injector: model1
    Added 1 new parcels

Cloud: cat
    Current number of parcels       = 2
    Current mass in system          = 1.04145e-09
    Linear momentum                 = (4.5255e-15 3.02526e-15 4.09251e-11)
   |Linear momentum|                = 4.09251e-11
    Linear kinetic energy           = 9.02176e-13
    model1:
        number of parcels added     = 2
        mass introduced             = 4.18879e-09
    Parcels absorbed into film      = 0
    New film detached parcels       = 0
    Parcel fate (number, mass)
      - escape                      = 0, 0
      - stick                       = 0, 0
    Temperature min/max             = 300, 300
    Mass transfer phase change      = 0
    D10, D32, Dmax (mu)             = 100, 100, 100
    Liquid penetration 95% mass (m) = 0.000996569
    hold-up of gasphase             = 0.000230061 %

BUT: Only my first cloud "bubbles" is saved to my lagrangian directory.... the second one "cat" is missing....

Any ideas?

Edit: OK, weird but I just renamed from "cat" to "catalyst" and it is working now.

Best regards
Andreas

[Danath]

i 've done the same

the spheres from the first cloud pass through the spheres from the second cloud


any ideas ?

Originally Posted by a.weber

Hi,
I'm also planning to do two clouds of totally different properties.... e.g. gas bubbles and catalysator particles.
Differences in calculation are important because:
bubbles should collide and coalesce with each other, cat particles should not.

What i did so far is to just add a second cloud in the createClouds.C:

Info<< "\nConstructing bubble cloud" << endl;
basicSprayCloud bubbles
(
    "bubbles",
    rho,
    U,
    g,
    slgThermo
);

Info<< "\nConstructing catalysator cloud" << endl;
basicSprayCloud cat
(
    "cat",
    rho,
    U,
    g,
    slgThermo
);

Then I just solve them both in my solver.C:

Info<< "Time = " << runTime.timeName() << nl << endl;
        bubbles.evolve();
        cat.evolve();
...

This works indeed, I get Info about both Clouds in the Output:

Solving 3-D cloud bubbles

--> Cloud: bubbles injector: model1
    Added 2 new parcels

Cloud: bubbles
    Current number of parcels       = 3
    Current mass in system          = 8.55318e-08
    Linear momentum                 = (7.09132e-11 1.15702e-11 2.83297e-08)
   |Linear momentum|                = 2.83298e-08
    Linear kinetic energy           = 4.69921e-09
    model1:
        number of parcels added     = 3
        mass introduced             = 8.55318e-08
    Parcels absorbed into film      = 0
    New film detached parcels       = 0
    Parcel fate (number, mass)
      - escape                      = 0, 0
      - stick                       = 0, 0
    Temperature min/max             = 300, 300
    Mass transfer phase change      = 0
    D10, D32, Dmax (mu)             = 3555.08, 3555.08, 3555.11
    Liquid penetration 95% mass (m) = 0.00136594
    hold-up of gasphase             = 0.000230061 %


Solving 3-D cloud cat

--> Cloud: cat injector: model1
    Added 1 new parcels

Cloud: cat
    Current number of parcels       = 2
    Current mass in system          = 1.04145e-09
    Linear momentum                 = (4.5255e-15 3.02526e-15 4.09251e-11)
   |Linear momentum|                = 4.09251e-11
    Linear kinetic energy           = 9.02176e-13
    model1:
        number of parcels added     = 2
        mass introduced             = 4.18879e-09
    Parcels absorbed into film      = 0
    New film detached parcels       = 0
    Parcel fate (number, mass)
      - escape                      = 0, 0
      - stick                       = 0, 0
    Temperature min/max             = 300, 300
    Mass transfer phase change      = 0
    D10, D32, Dmax (mu)             = 100, 100, 100
    Liquid penetration 95% mass (m) = 0.000996569
    hold-up of gasphase             = 0.000230061 %

BUT: Only my first cloud "bubbles" is saved to my lagrangian directory.... the second one "cat" is missing....

Any ideas?

Edit: OK, weird but I just renamed from "cat" to "catalyst" and it is working now.

Best regards
Andreas

[a.weber]

In my case, im ok with two seperate clouds, it was intended to not interfrere with each other.

I don't think that there is an easy way to make two clouds collide.. at least not in the original code.

You could try coding a force that is based on theta() (phase fraction) to push away some particles in crowded areas, but it would not stop spheres from moving over each other...

best regards,
Andy