
[Sponsors] 
December 26, 2020, 11:07 
Particles following the flow trajectory

#1 
New Member
Shabi
Join Date: Dec 2020
Posts: 16
Rep Power: 5 
Consider a simple 2D channel where the flow enters the domain from one side with the inlet velocity boundary condition with a constant value and it forms a boundary layer.
1) If i add many number of particles consistently to the inlet and neglect all inertial forces and let it to follow the flow trajectories, is it true that the volume fraction of the particle in the domain (C) is the solution of the advection equation, dC/dt + div (vC) = 0? If so, 2) If the inlet boundary condition for C is the constant value of 1, we know that in the steady state whole domain will be filled with the uniform value of C=1. But, it contradicts with what we expect with the particle case. Because if we add particles to the inlet consistently, because of the boundary layer, we expect a nonuniform value for the volume fraction in the steady state, because it is not possible to have a uniform distribution of particles in the domain in the steady state. Thank you in advance for sharing your ideas. 

December 26, 2020, 11:54 

#2 
Senior Member
Lucky
Join Date: Apr 2011
Location: Orlando, FL USA
Posts: 5,740
Rep Power: 66 
1) Yes
2) Recall conservation laws. If you add C=1 at the inlet, and there are no sources/sinks, then C=1 everywhere. 

January 4, 2021, 07:08 

#3 
New Member
Shabi
Join Date: Dec 2020
Posts: 16
Rep Power: 5 
Thanks. But how the particles can go to the area in the domain where the streamline is nearly zero? assume we have a solid circle in the middle of the domain and we have a moderate Reynolds number (no vortex behind the circle). So how the particles can fill the area behind the circle?


January 4, 2021, 10:26 

#4 
Senior Member
Join Date: Jul 2009
Posts: 357
Rep Power: 19 
If you have a boundary layer in the channel, then wouldn't there be gradients of the velocity near the walls? and wouldn't these gradients affect the distribution of C, so that your distribution is not uniform in the steady state?


January 4, 2021, 18:31 

#5  
New Member
q
Join Date: Sep 2020
Posts: 19
Rep Power: 6 
Quote:
How do you imagine particles while injecting just a single phase defining c=1 at the inlet? There is no distribution of c since there is only one phase. 

January 5, 2021, 04:29 

#6  
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,839
Rep Power: 73 
Quote:
You are wrong, the particle injection obeys to the lagrangian equation Dc/Dt=0 in a single phase, that is each particle follows a trajectori. In a divergence free velocity field, the eulerian equation counterpart is dc/dt+ div(vc)=0, if you integrate over the domain and apply Gauss, you see at the steady state that the inlet of c is balanced by the outlet of c. At the steady state (if exists), the trajectory and the streamlines are the same and c is constant along the streamline in the BL. 

January 5, 2021, 04:45 

#7  
New Member
q
Join Date: Sep 2020
Posts: 19
Rep Power: 6 
Quote:
I did not find the claim that he injects particles. Just assumes c=1 on the inlet. 

January 5, 2021, 04:47 

#8 
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,839
Rep Power: 73 

January 5, 2021, 04:50 

#9 
New Member
q
Join Date: Sep 2020
Posts: 19
Rep Power: 6 
Yes, I saw it. As can be seen I explicitly referred my reply to that statement: "2) If the inlet boundary condition for C is the constant value of 1, we know that in the steady state whole domain will be filled with the uniform value of C=1. But it contradicts with what we expect with the particle case... " Of course, it contradicts. So where am I wrong? 

January 5, 2021, 05:15 

#10  
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,839
Rep Power: 73 
Quote:
No contradiction... consider for example the analytic solution for the Poiseulle case. Immagine a uniform distribution of particles at the inlet and let them evolve until the particles reach the outlet. The particles close to a wall will have a slower velocity but they will finally reach the outlet so that all the domain is full of particles. The solution is steady and uniform. Of course, this model does not take into account a finite size of the particles. Now consider the Eulerian counterpart dc/dt+div (vc)=0. The uniform distribution of particles can be obtained by setting the uniform value c=1 at the inlet. Again if you wait until the transient is ended the final solution will show c=1 everywhere. There is no diffusion and no compressibility effects to change this value. 

January 5, 2021, 05:26 

#11 
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,839
Rep Power: 73 

January 21, 2021, 17:57 

#12  
New Member
Shabi
Join Date: Dec 2020
Posts: 16
Rep Power: 5 
Quote:
Or lets formulate like this. Assume that we don't have a constant inlet velocity (we have e.g. a parabolic one) and then if we inject the particles uniformly, we can clearly see the difference between the volume fraction of the particles near the inlet. However, if we solve the advection equation, we see the uniform constant value of C in the steady state in the whole domain... 

January 22, 2021, 04:53 

#13  
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,839
Rep Power: 73 
Quote:


January 22, 2021, 06:41 

#14  
New Member
Shabi
Join Date: Dec 2020
Posts: 16
Rep Power: 5 
Quote:
But if we solve the advection equation exactly with the same configuration, we will get the same value of "C" in the whole domain. 

January 22, 2021, 11:24 

#15  
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,839
Rep Power: 73 
Quote:
But you are considering from a side a discrete particle representation and on the other a continuous representation of the convection discretized on a grid. What you see in the latter case is the same phenomenology, only differently represented. 

January 23, 2021, 17:21 

#16  
New Member
Shabi
Join Date: Dec 2020
Posts: 16
Rep Power: 5 
Quote:
What we know so far: 1) The number of particles in some area per unit volume is more than other area. 2) the scalar field 'C' is constant everywhere in the steady state. Recalling the start of this topic, i wanted to understand the relation between the volume fraction of the particles and the scalar field 'C'. I think if 1 and 2 is correct, the only explanation for this conflict is that the number of particles per unit volume does not represent the volume fraction. I hope i could clarify my question. you say that the convection discretized on a grid is the same phenomenology. But, the number of particles per unit volume is not constant, whereas 'C' is constant everywhere. Last edited by Shabi; January 24, 2021 at 06:18. 

January 24, 2021, 07:32 

#17  
Senior Member
Lefteris
Join Date: Oct 2011
Location: UK
Posts: 341
Rep Power: 15 
Quote:
I can't make sense of what you're asking. Would you like to clarify this please?
__________________
Lefteris 

January 24, 2021, 08:19 

#18  
New Member
Shabi
Join Date: Dec 2020
Posts: 16
Rep Power: 5 
Quote:
What is clear is how a advection equation is solved. Assuming that the inlet boundary condition of C is 1, we know that in the steady state C is 1 in the whole domain. So, i would like to resemble this process with very small particles, i.e., transport very small particles and make them to follow the flow trajectories. Those particles are not going to have any influence on the flow. They only follow the flow trajectories. Similar to what the continuous scalar field 'C' does. if the number of my particles reaches infinity, i expect to converge to the continuous field 'C'. 

January 24, 2021, 08:28 

#19  
Senior Member
Lefteris
Join Date: Oct 2011
Location: UK
Posts: 341
Rep Power: 15 
Quote:
__________________
Lefteris 

January 24, 2021, 08:34 

#20 
New Member
Shabi
Join Date: Dec 2020
Posts: 16
Rep Power: 5 
I am not sure, actually i am not doing any real simulation using any method. It is more related to understanding the basics. I saw the conflict here, and i cannot find the answer for that. We are in an incompressible flow. 'C' fills the domain gradually and in the steady state we get a uniform distribution of 'C'. And it is what we expect from the advection equation and the incompressibility. So, why can't we repeat this process with infinite number of tiny 'massless' particles?


Tags 
advection equation, particle concentration, transport equation, volume fraction. 
Thread Tools  Search this Thread 
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Will the results of steady state solver and transient solver be same?  carye  OpenFOAM Running, Solving & CFD  9  December 28, 2019 05:21 
How to simulate particles in a gas flow  sara  OpenFOAM Running, Solving & CFD  13  October 8, 2019 05:12 
Choosing a solver to model particles at rest in a flow field  Nigel Bruce Khan  OpenFOAM Running, Solving & CFD  0  March 22, 2019 08:35 
Particles in the flow  nissen21  Fluent Multiphase  5  May 21, 2014 19:21 
Interpolate particles together with flow  asal  FLUENT  0  March 30, 2014 14:22 