[msman]

I have used the same geometry for other solver and its works.

Files are attached herewith where i have made modifications according to my case

[msman]

Hi again,

can anyone explain the meaning of some terms used in KimematicCloudProperties dictionary

type patchInjection; //Type of injection
parcelBasisType fixed; // ?
patchName upperInlet; // start patch for injection
U0 (18.7939 -6.8404 0); // ?
nParticle 1; // 1 is for steady state condition? right?
parcelsPerSecond 1390885; // number of particles entering per second
sizeDistribution // ?
{
type normal; // ?
normalDistribution
{
expectation 650e-6; // ?
variance 25e-6;
minValue 500e-6;
maxValue 800e-6;
}
}
flowRateProfile constant 1; // ?
massTotal 0; // ?
SOI 0; // start of injection
duration 60; // ?

Kindly rectify my mistake

Your kind help is highly appreciated

Thank you

[ansubru]

Hi!!

I will try and answer the terms you have marked with a question mark -

Hi again,

can anyone explain the meaning of some terms used in KinematicCloudProperties dictionary

type patchInjection; //--> Type of injection (you can have either patch injection where you define a surface in the geometry from where particles are injected or other types, have look at injection types.
parcelBasisType fixed; // --> Don't really concern yourself about this, leave it as FIXED
patchName upperInlet; // start patch for injection (the name of the patch where injection is initiated)
U0 (18.7939 -6.8404 0); // Particle velocity at injecton
nParticle 1; // This is a condition for parcels, if you change this to say 100, it means one parcel represensts 100 particles, if this value is 1 then you are treating the system as particle based and not parcel based
parcelsPerSecond 1390885; // number of particles entering per second (correct ))
sizeDistribution // If you want a specific particle size distribution you specify this condition (refer size distribution types in OF guides online)
{
type normal; // Type of distribution can be normal, rossin-ramler etc. (read more in OF guides)
normalDistribution // (read about normal distribution to understand what these terms) mean
{
expectation 650e-6; // ?
variance 25e-6;
minValue 500e-6;
maxValue 800e-6;
}
}
flowRateProfile constant 1; // --> type of flow rate profile again for most laminar flow cases this would be 1, else you will have to change them
massTotal 0; //--> This is used to fix the number of particles that should be injected into the system (or you can use it to control the rate of particles injected per time-step).. play around with it to better understand
SOI 0; // start of injection (correct)
duration 60; // --) Total duration of the injection, in this case the particles are injected until 60 secs of simulation time



I would advice you to thoroughly understand the tutorial case, before proceeding. Everything you need is there, and you should also dig into the source code for the 'lagrangian' function in openfoam (Openfoam x./src/lagrnagian..)

Regards

ansubru

[ansubru]

Hi again!!

Looking at your geometry, I can see some potential regions where you might have numerical inconsistencies. Your geometry is too open-ended. I don't know if this is correct, but you could try and close the pores (they seem to open in the sides).. ensure that the pores have just an opening in the center and the remaining portions are fully sealed (mainly around the circumference)

Regards

ansubru

[msman]

Hi again,

thanks ansubru for you kind help. I am about to successfully simulate my first part of the project.

These open pores shown in last figure are completed sealed.

Now I wanna know about the boundry condition for particle located in kinematicCloudproperties dictionary. Usally for velocity and pressure the Bcs for 2D case are like so

walls
{
type empty;
}


but i am not sure what should be the bc for particles and also I couldn't find from langrangian tutorials too. currently I am using

walls
{
type stick;

}

Is it correct or should I also need to use it "empty" for particles?

[ansubru]

Hi again!!

The properties of the Lagrangian phase are specified in the kinematiccloudproperties file. Generally you can specify, stick or rebound or escape depending on your system. Generally particles would rebound from walls, but if you specifically want to capture them ( for some special reasons) then you can use stick. The conditions are specie based on what the system requires.

Regards

ansubru

[panachristos]

Hallo,

patchInjection injects the parcels randomly over the patch, which is a randomness to my results that i do not want. So, is there any way to modify patchInjection so that it distributes the parcels over the patch under a certain pattern??

If not possible is there any other model to be used instead? manualInjection is not an option for me since i have thousands of particles.

I am looking forward to give me any suggestions on how to solve this.

Regards,
Christos

[dussa]

Hi,

I see that this thread has not been used in a while, but hopefully someone can help me out.

I am trying to simulate a solid particle and fluid system, where the fluid is more dense than the particles causing them to float. I have just now tried taking the DPMFoam/Goldschmidt case and changing only the fluid and particle density (I also reduced the domain size and number of particles injected to speed up the simulation). I also removed the inlet velocity for further simplicity.

When the simulation is with air the particles fall to the bottom of the geometry as expected and quickly lose most of their kinetic energy forming a layer with very little movement.

When I try this with the case altered by using water (rho = 1000 kg/m3) and reducing the particle density (rho0 = 750) I find that the particles move in a very erratic manner, and do not simply float from their injection location to the top of the fluid domain. They seem to be attracted to the walls, and away from the centre of the fluid domain as you can see in the attached images after 1 second of physical time.

Does anyone know why this could be? I have seen some posts talk about additional forces that need to be included if a dense fluid like water is used, but I am not sure which ones would be required here. If someone coudld point me in the direction of where I can look at what contact models, fluid-particle forces models etc are available that would be much appreciated.

Kind regards,
Robert