CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM

DieselFoam spray model

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

Like Tree2Likes
  • 2 Post By qndfekjb

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   February 6, 2011, 20:41
Default DieselFoam spray model
  #1
New Member
 
Viliya Shen
Join Date: Feb 2011
Posts: 1
Rep Power: 0
qndfekjb is on a distinguished road
Dear All,

I'd like to simulate diesel spray and ignition for
my diploma thesis . SInce im a beginner i want to run some test case like Aachenbomb.
Where can I find these test cases ?
or hat anyone some spray cases, and can you so kindly share with me?
qndfekjb@gmail.com

Thank you.
Sincerely,

Viliya
kebi112358 and saeedshayae like this.
qndfekjb is offline   Reply With Quote

Old   February 7, 2011, 01:38
Default
  #2
Senior Member
 
Join Date: Jun 2009
Posts: 100
Rep Power: 16
alastormoody11 is on a distinguished road
Hi,

there is a very good presentation detailing spray breakup using spray breakup models ETAB, TAB, CAB using ANSYS CFX, the presentation contains images of the spray development at specific instants as well as certain quantitative data. So you might check that out, it will be better to ensure that your models as correctly simulating the spray breakup and development before going for a combustion simulation.
alastormoody11 is offline   Reply With Quote

Old   February 7, 2011, 10:05
Default
  #3
Senior Member
 
Elvis
Join Date: Mar 2009
Location: Sindelfingen, Germany
Posts: 620
Blog Entries: 6
Rep Power: 24
elvis will become famous soon enough
You will find aachenBomb
OpenFOAM-1.5-dev:
https://openfoam-extend.svn.sourcefo...am/aachenBomb/
OpenFOAM-1.6-ext:
http://openfoam-extend.git.sourcefor...39082b;hb=HEAD

Last edited by elvis; February 7, 2011 at 10:23.
elvis is offline   Reply With Quote

Old   February 29, 2012, 08:01
Default
  #4
Senior Member
 
shoeb khan
Join Date: Nov 2011
Posts: 179
Rep Power: 14
shk12345 is on a distinguished road
Quote:
Originally Posted by alastormoody11 View Post
Hi,

there is a very good presentation detailing spray breakup using spray breakup models ETAB, TAB, CAB using ANSYS CFX, the presentation contains images of the spray development at specific instants as well as certain quantitative data. So you might check that out, it will be better to ensure that your models as correctly simulating the spray breakup and development before going for a combustion simulation.
can you give some links about those presentation.I shall be highly greatfull.

thanks
shk
shk12345 is offline   Reply With Quote

Old   February 29, 2012, 14:09
Default
  #5
Senior Member
 
mturcios777's Avatar
 
Marco A. Turcios
Join Date: Mar 2009
Location: Vancouver, BC, Canada
Posts: 740
Rep Power: 28
mturcios777 will become famous soon enough
If you want to run diesel spray simulations in OF 2.1.x, be aware that the dieselSpray class has been deprecated and is not longer included. The new solver sprayFoam is now used, you can find the AachenBomb tutorial in $FOAM_TUTORIALS/lagrangian/sprayFoam/aachenBomb.
mturcios777 is offline   Reply With Quote

Old   March 5, 2012, 15:09
Default
  #6
New Member
 
SamboReal
Join Date: Feb 2011
Posts: 9
Rep Power: 15
sambo is on a distinguished road
Be careful using sprayFoam.
Due that boiling is not included, the lagrangian spray model is only for low temperature effects.
Check this link: http://www.openfoam.org/mantisbt/view.php?id=346
sambo is offline   Reply With Quote

Old   March 5, 2012, 17:00
Default
  #7
Super Moderator
 
niklas's Avatar
 
Niklas Nordin
Join Date: Mar 2009
Location: Stockholm, Sweden
Posts: 693
Rep Power: 29
niklas will become famous soon enoughniklas will become famous soon enough
not if you're using git, its been added

https://github.com/OpenFOAM/OpenFOAM...26abd04b29c583
niklas is offline   Reply With Quote

Old   March 5, 2012, 17:13
Default
  #8
New Member
 
SamboReal
Join Date: Feb 2011
Posts: 9
Rep Power: 15
sambo is on a distinguished road
Thank you very much niklas.
I did not know that
sambo is offline   Reply With Quote

Old   March 5, 2012, 18:12
Default multi-hole nozzle
  #9
New Member
 
SamboReal
Join Date: Feb 2011
Posts: 9
Rep Power: 15
sambo is on a distinguished road
By the way.
Is there a multi-hole nozzle model (5 to 7).
So far I just used the coneNozzleInjection and did not found any multi-hole nozzle model (/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel)
sambo is offline   Reply With Quote

Old   March 6, 2012, 09:03
Default
  #10
Super Moderator
 
niklas's Avatar
 
Niklas Nordin
Join Date: Mar 2009
Location: Stockholm, Sweden
Posts: 693
Rep Power: 29
niklas will become famous soon enoughniklas will become famous soon enough
No there is no multihole injector available at the moment.

I have written my own injector though.
It works almost the same as the one available in the old dieselSpray library.
The difference is that it doesnt inject a parcel from each hole every time,
but it randomly picks one of the available holes to inject the parcel from.
That means that you should increase the injected parcels to get better
statistics and that you can have small differences in the injected mass
between the holes.

If you want to test it, you can unpack it under
src/lagrangian/intermediate/submodels/Kinematic/InjectionModel

and then, in the file
src/lagrangian/intermediate/parcels/include/makeReactingParcelInjectionModels.H

you add
#include "MultiHoleInjector.H"
to the include-statements

and a bit lower you add the MultiHoleInjection-line to the makeInjectionModel templates.

Code:
     makeInjectionModelType(PatchInjection, CloudType);                        \
    makeInjectionModelType(MultiHoleInjection, CloudType);                    \
    makeInjectionModelType(ReactingLookupTableInjection, CloudType);
then its simplest to just do a wclean, wmake libso, in the intermediate directory, to recompile
the library with the new injector type.
You might have to recompile sprayFoam as well, but you will see that yourself
if the MultiHoleInjector doesnt appear among the available injector types)
Attached Files
File Type: gz MultiHoleInjection.tar.gz (3.5 KB, 98 views)
niklas is offline   Reply With Quote

Old   March 8, 2012, 17:02
Default
  #11
New Member
 
SamboReal
Join Date: Feb 2011
Posts: 9
Rep Power: 15
sambo is on a distinguished road
Thank you very much niklas. I wiil have a look.
Actually do you use the "ORourke" stochasticCollisionModel.
If you use use, which coalescence is sensible. I do not really understand this value even though I had a look at the ORourke equations.
sambo is offline   Reply With Quote

Old   March 9, 2012, 02:26
Default
  #12
Super Moderator
 
niklas's Avatar
 
Niklas Nordin
Join Date: Mar 2009
Location: Stockholm, Sweden
Posts: 693
Rep Power: 29
niklas will become famous soon enoughniklas will become famous soon enough
coalescence is just a boolean
so if you set it to on, parcels can coalesce,
if you set it to off, they will not exchange mass, only momentum.
niklas is offline   Reply With Quote

Old   August 20, 2012, 07:43
Default Droplet coalescence when collision model is off
  #13
New Member
 
fairus
Join Date: Apr 2010
Posts: 20
Rep Power: 15
mfmohdyasin is on a distinguished road
Dear Niklas,

I simulated a non-evaporating spray and turn off the collision model,
evaporation model, and heat transfer. I'm using TAB breakup model.
However, when I sample the number of droplets in axial direction over
a period of steady spray, I have a decrease in droplet number as shown
below (The number of droplet has been normalized, it is in the order of 1000).Is the decrease in droplet number due to coalescence?
Please comment.

Best,
Fairus
Attached Images
File Type: png axialDropletcompare.png (13.6 KB, 88 views)
mfmohdyasin is offline   Reply With Quote

Old   August 20, 2012, 07:53
Default
  #14
Super Moderator
 
niklas's Avatar
 
Niklas Nordin
Join Date: Mar 2009
Location: Stockholm, Sweden
Posts: 693
Rep Power: 29
niklas will become famous soon enoughniklas will become famous soon enough
How can it be if it is turned off?
how do you do the sampling?
niklas is offline   Reply With Quote

Old   August 20, 2012, 08:57
Default
  #15
New Member
 
fairus
Join Date: Apr 2010
Posts: 20
Rep Power: 15
mfmohdyasin is on a distinguished road
Once the number of parcel reached steady state, I start outputting the result over a period of time.
Than, I sample the droplets from each time directory according to their axial positions and the number of droplet is the same as the total number of parcel in respective axial position.

I'm not sure how the droplet number can be reduced. Obviously, we expect the droplet number to increase in axial direction since there is no means for it to decrease since evaporation and coalescense are off. No droplets are removed through the domain boundaries between the spray BC location and
the outlet.

Fairus
mfmohdyasin is offline   Reply With Quote

Old   August 20, 2012, 08:58
Default
  #16
Super Moderator
 
niklas's Avatar
 
Niklas Nordin
Join Date: Mar 2009
Location: Stockholm, Sweden
Posts: 693
Rep Power: 29
niklas will become famous soon enoughniklas will become famous soon enough
how do you do the sampling?
code please
niklas is offline   Reply With Quote

Old   August 20, 2012, 11:29
Default
  #17
New Member
 
fairus
Join Date: Apr 2010
Posts: 20
Rep Power: 15
mfmohdyasin is on a distinguished road
Niklas,

I'm using the following shell scripts to do the sampling. There are 2 main scripts that need to be run in the following sequence:

1) from the case directory:

./runppc2Atom (sample and put droplets in different axial positions in "ppd" folder)

2) from the ppd folder:
./breakupDropletSample (output the number of droplet vs axial position)

Sorry if it appears a bit untidy. I've to separate them because the script need to serve multiple purposes, so I need to have them separated.

Fairus
Attached Files
File Type: txt runppc2Atom.txt (676 Bytes, 39 views)
File Type: txt postLagPrep.txt (1.6 KB, 19 views)
File Type: txt postLag.txt (1.2 KB, 21 views)
File Type: txt breakupDropletSample.txt (1.5 KB, 35 views)
mfmohdyasin is offline   Reply With Quote

Old   August 27, 2012, 17:03
Default
  #18
New Member
 
Join Date: Jun 2012
Posts: 14
Rep Power: 13
bigeddy is on a distinguished road
Hey mfmohdyasin, Usually for cfd codes involving sprays and specifically breakup models, its very costly to simulate a droplet with a diameter with order say for e.g 1E-6, and its rather pointless. With brekaup models, hundreds of these small droplets would be produced especially if your spray is supposed to be steady state. So if a droplet get to this size, the code usually adds its mass and momentum to the continuous phase and subsequently removes it from the computation.

This might explain what you are seeing. Look for some sort of parameter in openfoam that specifies a minimum droplet diameter before removal from domain. Myt be better to look into D10 and D32 droplet diameters sampled across your domain space.

Also are you counting droplet parcels (the lagrangian objects) or the number of real droplets contained within each parcel?? Breakup models modify the latter parameter depending on the new droplet diameter computed via the breakup model procedure in order to conserve mass in the parcel.

Last edited by bigeddy; August 27, 2012 at 17:31.
bigeddy is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Superlinear speedup in OpenFOAM 13 msrinath80 OpenFOAM Running, Solving & CFD 18 March 3, 2015 06:36
DieselFoam spray thumthae OpenFOAM Running, Solving & CFD 98 December 24, 2014 16:55
Problem: hollow cone spray with TAB model. chenhu OpenFOAM 1 July 3, 2009 12:25
axi-symmetric mesh of dieselFoam az1362f OpenFOAM Running, Solving & CFD 2 June 27, 2009 12:32
the role of the atomization model Jay Siemens 2 March 7, 2008 21:16


All times are GMT -4. The time now is 13:05.