CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

Injection of Lagrangian particles

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   November 23, 2010, 03:28
Default Injection of Lagrangian particles
  #1
Senior Member
 
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,894
Rep Power: 26
alberto will become famous soon enoughalberto will become famous soon enough
Hello,

I have a doubt about the Lagrangian classes. I would like to simulate the injection of two jets of particles using the O'Rourke collision model.

My case is extremely simple: particles are injected from two patches at constant velocity and constant diameter, and they can collide.

It seems collision models is only available using the spray sub-models, which adds the complication of using injectors (there does not seem to be an injector which uses a patch, at least if my understanding is correct).

Could you give some hint to set this case up?

Thanks
__________________
Alberto Passalacqua

GeekoCFD - A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as live DVD/USB, hard drive image and virtual image.
OpenQBMM - An open-source implementation of quadrature-based moment methods
alberto is offline   Reply With Quote

Old   November 23, 2010, 19:33
Default
  #2
Member
 
sradl's Avatar
 
Stefan Radl
Join Date: Mar 2009
Location: Graz, Austria
Posts: 82
Rep Power: 9
sradl is on a distinguished road
Hey Alberto!

I did some LPT using injections in a region (e.g., a sphere), but not in a domain. I imagine one could write a routine do inject particles near an inlet.

I'd be interested to work on an implementation of the MP-PIC approach of O'Rourke into twoPhaseEulerFoam. We have started a project on that here recently, and would love to share code of LPT (coupled to LES, using an eddy interaction model)

cheers
Stefan
sradl is offline   Reply With Quote

Old   November 23, 2010, 22:38
Default
  #3
Senior Member
 
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,894
Rep Power: 26
alberto will become famous soon enoughalberto will become famous soon enough
Hi Stefan,

thanks for your answer. I guess definedInjector can do that. My case is a pure Lagrangian simulation without fluid, so I do not really need an inlet, just an area next to a wall from where particles come out. I was trying to find the easiest way to do this, since injectors want a whole bunch of information I do not really need! :-)

About MP-PIC, you probably need a more robust solver for dense flows than twoPhaseEulerFoam, since the key in some of the recent implementations of MP-PIC is enforcing the particle packing by means of the Eulerian equations, to evaluate the acceleration term responsible of preventing packing implicitly.
I saw you joined my project on extend, so you know I am working on that ;-)

What kind of applications would you like to do with MP-PIC (I saw you are in Prof Sundaresan's group, so we might have a chance to meet sometime).

Best,
__________________
Alberto Passalacqua

GeekoCFD - A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as live DVD/USB, hard drive image and virtual image.
OpenQBMM - An open-source implementation of quadrature-based moment methods
alberto is offline   Reply With Quote

Old   November 24, 2010, 05:04
Default
  #4
New Member
 
Astrid Mahrla
Join Date: May 2010
Location: Munich, Germany
Posts: 22
Rep Power: 8
AMahrla is on a distinguished road
Hi Alberto!

There is a PatchInjectionModel in OpenFOAM..
It's used in the solver rhoPisoTwinParcelFoam with KinematicParcel class and can be found in src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/PatchInjection..

Hope this helps!

Best,

Astrid
AMahrla is offline   Reply With Quote

Old   November 24, 2010, 09:57
Default
  #5
Member
 
sradl's Avatar
 
Stefan Radl
Join Date: Mar 2009
Location: Graz, Austria
Posts: 82
Rep Power: 9
sradl is on a distinguished road
Hey Alberto,

thanks for the useful input. Yeah, I'm aware of the stability issues associated with MP-PIC - Snider is working on that for more than 10 years...

I think the best way is to have a stable and two-fluid solver, and knowing how to stabilize the particle pressure calculation near the packing limit. I'd be glad to contribute to this development in the twoPhaseEulerPimpleFoam code.

Are there already implementations of MP-PIC in OF available to have a quick look at that?

My goal is compare MP-PIC with E-E based on filtered simulations (similar to the work of Igci et al.), and then go to CFB.

cheers
Stefan
sradl is offline   Reply With Quote

Old   November 24, 2010, 11:16
Default
  #6
Senior Member
 
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,894
Rep Power: 26
alberto will become famous soon enoughalberto will become famous soon enough
Quote:
Originally Posted by sradl View Post
Hey Alberto,

thanks for the useful input. Yeah, I'm aware of the stability issues associated with MP-PIC - Snider is working on that for more than 10 years...
Yes, I do not have a hands-on view (never coded anything related to MP-PIC), but I think their issues are related to the exchange between Euler-Euler and Euler-Lagrange framework, which forces them to use quite strong relaxation. A similar problem was discussed at ICMF last year by a "famous commercial code vendor", which is implementing MP-PIC under another name

Quote:
I think the best way is to have a stable and two-fluid solver, and knowing how to stabilize the particle pressure calculation near the packing limit. I'd be glad to contribute to this development in the twoPhaseEulerPimpleFoam code.
Thanks

Quote:
Are there already implementations of MP-PIC in OF available to have a quick look at that?
Not that I know. I remember seeing some thread talking about MP-PIC on this forum, but it was about intentions of implementing it.

Quote:
My goal is compare MP-PIC with E-E based on filtered simulations (similar to the work of Igci et al.), and then go to CFB.
That's interesting!

Best,
__________________
Alberto Passalacqua

GeekoCFD - A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as live DVD/USB, hard drive image and virtual image.
OpenQBMM - An open-source implementation of quadrature-based moment methods
alberto is offline   Reply With Quote

Old   November 24, 2010, 11:17
Default
  #7
Senior Member
 
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,894
Rep Power: 26
alberto will become famous soon enoughalberto will become famous soon enough
Quote:
Originally Posted by AMahrla View Post
Hi Alberto!

There is a PatchInjectionModel in OpenFOAM..
It's used in the solver rhoPisoTwinParcelFoam with KinematicParcel class and can be found in src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/PatchInjection..
Thanks Mahrla
__________________
Alberto Passalacqua

GeekoCFD - A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as live DVD/USB, hard drive image and virtual image.
OpenQBMM - An open-source implementation of quadrature-based moment methods
alberto is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Injection of specific number of particles in DPM DH FLUENT 11 October 20, 2011 07:50
Sampling tool for Lagrangian particles sankarv OpenFOAM 1 July 1, 2010 11:34
How to set correct mass flow rate for lagrangian particles ? sankarv OpenFOAM 0 April 19, 2010 11:40
Mapping Lagrangian particles gschaider OpenFOAM Bugs 1 April 1, 2009 05:05
injection particles zenith FLUENT 0 April 3, 2007 00:56


All times are GMT -4. The time now is 04:49.