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

Description and Implementation of SprayFoam

Register Blogs Community New Posts Updated Threads Search

Like Tree2Likes
  • 1 Post By mm66
  • 1 Post By mm66

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   February 19, 2015, 12:46
Cool Description and Implementation of SprayFoam
  #1
New Member
 
Juan David Rodriguez P
Join Date: Jan 2015
Location: Milano
Posts: 20
Rep Power: 11
JuanRodriguez is on a distinguished road
Hi.

I have read already Anton Persson's report describing dieselFoam. I have noticed this solver and sprayFoam have many differences.

Is there a description/report of sprayFoam>

I ask this becuase I dont know what does Constant Properties in sprayCloudProperties refer to. Also, I wan to know how to set the spray injection temperature and pressure. I know in 0/T or 0/p one can set boundary conditions for these scalars but for the general duct where the spray takes place (the mass of air, for example)

Thank you
JuanRodriguez is offline   Reply With Quote

Old   February 19, 2015, 14:00
Default
  #2
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
You can set the spray temperature (stagnation temperature) in T0. If you want a pressure driven flow in the spray you need to select a coneNozzleInjector as the type, and flowType as pressureDrivenVelocity. You will need a keyword called Pinj which is the injection pressure.
mturcios777 is offline   Reply With Quote

Old   February 20, 2015, 12:54
Lightbulb
  #3
New Member
 
Juan David Rodriguez P
Join Date: Jan 2015
Location: Milano
Posts: 20
Rep Power: 11
JuanRodriguez is on a distinguished road
Hi Marcos,

Do you know what parameter describes liquid evaporation? I used Y(component)l in paraview but its scale was from 1 to 2, which seems odd, does not it? If it is not odd at all, then what does it mean?

By the way, I am having trouble setting the BC to get an inner flow of air. Could you please give me some guide with respect to this?

Thank you.
Attached Files
File Type: gz 0.tar.gz (1.6 KB, 45 views)
JuanRodriguez is offline   Reply With Quote

Old   February 23, 2015, 17:23
Default
  #4
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
What do you mean by the liquid evaporation? The species fractions, parcel mass, etc? Ycomponent is the index of the species in the parcel, my guess is you have 3 species?
mturcios777 is offline   Reply With Quote

Old   February 23, 2015, 18:24
Default
  #5
New Member
 
Juan David Rodriguez P
Join Date: Jan 2015
Location: Milano
Posts: 20
Rep Power: 11
JuanRodriguez is on a distinguished road
Hi Marco,
Oh, I see. Yes indeed, I am working with three species (N2, O2 and Water).
My question regarding liquid evaporation was if there exists a variable in OF that represents a vapor fraction of the sprayed species, or that contains information about phase change.

Thank you, Marco.
JuanRodriguez is offline   Reply With Quote

Old   February 23, 2015, 18:50
Default
  #6
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
For lagrangian spray models, diameter would give you some indication of evaporation. The vapour fraction of sprayed species would be the species fraction in the continuous phase I would think. When you plot H2O it should show it increasing where the droplets are evaporating.
mturcios777 is offline   Reply With Quote

Old   February 25, 2015, 09:24
Default
  #7
New Member
 
Juan David Rodriguez P
Join Date: Jan 2015
Location: Milano
Posts: 20
Rep Power: 11
JuanRodriguez is on a distinguished road
Thank you, Marco.
I will follow your recommendation.
JuanRodriguez is offline   Reply With Quote

Old   February 26, 2015, 09:56
Default
  #8
New Member
 
Juan David Rodriguez P
Join Date: Jan 2015
Location: Milano
Posts: 20
Rep Power: 11
JuanRodriguez is on a distinguished road
I am investigating how sprayFoam works. And I have seen that inside the solver folder there are only pEqn.H, rhoEqn.h and UEqn.H, but sprayFoam.C uses also YEqn.H and EEqn.H. My question is where does sprayFoam.C load them from?

Thanks.
JuanRodriguez is offline   Reply With Quote

Old   March 1, 2015, 20:13
Default
  #9
Member
 
Neilson Whit
Join Date: Aug 2011
Posts: 74
Rep Power: 15
wolfindark is on a distinguished road
Quote:
Originally Posted by JuanRodriguez View Post
I am investigating how sprayFoam works. And I have seen that inside the solver folder there are only pEqn.H, rhoEqn.h and UEqn.H, but sprayFoam.C uses also YEqn.H and EEqn.H. My question is where does sprayFoam.C load them from?

Thanks.
I think It uses the ones reactingParcelFoam. Please check it...
wolfindark is offline   Reply With Quote

Old   March 9, 2015, 18:02
Default
  #10
New Member
 
Juan David Rodriguez P
Join Date: Jan 2015
Location: Milano
Posts: 20
Rep Power: 11
JuanRodriguez is on a distinguished road
I also have that intuition, Neilson. But I don't know how to confirm it.
Anyway, sprayFoam seems to calculate at each iteration these three measures, D32, D10, and Dmax. Other solvers, like reactingParcelFoam don't. How can I register these values with that solver? Is it possible?

Thanks/
JuanRodriguez is offline   Reply With Quote

Old   March 10, 2015, 03:47
Default
  #11
Member
 
Likun
Join Date: Feb 2013
Posts: 52
Rep Power: 13
Likun is on a distinguished road
Send a message via Skype™ to Likun
Neilson is correct. If you check ''sprayFoam.dep", you will find
Quote:
sprayFoam.dep: ../reactingParcelFoam/YEqn.H
sprayFoam.dep: ../reactingParcelFoam/EEqn.H
which show that sprayFoam uses YEqn.H and EEqn.H from reactingParcelFoam.

Best,
Likun
Likun is offline   Reply With Quote

Old   June 13, 2018, 10:00
Default sprayFoam, penetration, evaporation, breakup
  #12
New Member
 
heruitian
Join Date: Jun 2018
Posts: 6
Rep Power: 8
heruitian is on a distinguished road
Hi,
I am trying to simulate the n-dodecane evaporation based on sprayFoam in OF. With the C12H26 file in each time documents, I define the vapor penetration just like the definition in other commercial software. But the vapor penetration is shorter than the liquid penetration in the initial time. It's unacceptable in the experiments. I wonder if there is mistakes in vapor penetration definition or the setting values, such as KHRT parameters in OF,which can let the droplet evaporate faster in the initial time?
Thank you!
heruitian is offline   Reply With Quote

Old   October 29, 2020, 02:53
Default pressure bondary condition in sprayFoam
  #13
Member
 
Riddhideep Biswas
Join Date: May 2020
Posts: 30
Rep Power: 6
Rid@foam is on a distinguished road
Hello everyone.
I am simulating a case in sprayFoam and I want to keep an atmospheric pressure boundary condition to the outlet of the domain.
What value should I give for this?
Should I give 0 Pa(i.e. the gauge pressure) or should I give 1e5 Pa(i.e. the absolute pressure)?
Anyone please help me as I am new to sprayFoam solver and do not understand it well.
Rid@foam is offline   Reply With Quote

Old   November 10, 2020, 12:11
Question Pressure driven spray
  #14
Member
 
Join Date: Mar 2019
Posts: 81
Rep Power: 7
mm66 is on a distinguished road
Dear Foamers,

I am trying to simulate a pure pressure driven flow through a nozzle using speayFoam in OpenFoam1812. I would like to obtain the mass flow rate through the nozzle. Here is the relevant info from sprayCloudProperties file:
Code:
    injectionModels
    {
        model1
        {
            type            coneNozzleInjection;
            SOI             0;
            massTotal       1e-5;
            parcelBasisType mass;
            injectionMethod disc;
            flowType        pressureDrivenVelocity;
            Pinj	    constant 1e6;
            outerDiameter   2e-4;
            innerDiameter   0;
            duration        50e-3;
            position        (0 0.5 0);
            direction       (0 -1 0);
            parcelsPerSecond 20000000;
            flowRateProfile constant 1;
The problem is that (to my knowledge) both "massTotal" and "duration" should be defined. Consequently, wouldn't the mass flow rate be imposed to the domain as mass flow = "massTotal" / "duration"?

I really appreciate your thoughts on this.

Regards,
MJ
raj kumar saini likes this.
mm66 is offline   Reply With Quote

Old   November 10, 2020, 12:15
Default
  #15
Member
 
Join Date: Mar 2019
Posts: 81
Rep Power: 7
mm66 is on a distinguished road
Quote:
Originally Posted by Rid@foam View Post
Hello everyone.
I am simulating a case in sprayFoam and I want to keep an atmospheric pressure boundary condition to the outlet of the domain.
What value should I give for this?
Should I give 0 Pa(i.e. the gauge pressure) or should I give 1e5 Pa(i.e. the absolute pressure)?
Anyone please help me as I am new to sprayFoam solver and do not understand it well.
Dear Riddhideep,

I believe you should use the absolute pressure (something like the following):

Code:
dimensions          [1 -1 -2 0 0 0 0];

internalField       uniform 1e06;

boundaryField
{
    walls
    {
        type            zeroGradient;
    }
    outlet
    {
        type            fixedValue;
        value           uniform 1e06;
    }
}
Hope this helps

Cheers,
MJ
raj kumar saini likes this.
mm66 is offline   Reply With Quote

Old   March 9, 2022, 09:26
Default
  #16
Member
 
Join Date: Apr 2014
Location: N/A
Posts: 50
Rep Power: 12
FluentStarter is on a distinguished road
Quote:
Originally Posted by mm66 View Post
Dear Foamers,

I am trying to simulate a pure pressure driven flow through a nozzle using speayFoam in OpenFoam1812. I would like to obtain the mass flow rate through the nozzle. Here is the relevant info from sprayCloudProperties file:
Code:
    injectionModels
    {
        model1
        {
            type            coneNozzleInjection;
            SOI             0;
            massTotal       1e-5;
            parcelBasisType mass;
            injectionMethod disc;
            flowType        pressureDrivenVelocity;
            Pinj	    constant 1e6;
            outerDiameter   2e-4;
            innerDiameter   0;
            duration        50e-3;
            position        (0 0.5 0);
            direction       (0 -1 0);
            parcelsPerSecond 20000000;
            flowRateProfile constant 1;
The problem is that (to my knowledge) both "massTotal" and "duration" should be defined. Consequently, wouldn't the mass flow rate be imposed to the domain as mass flow = "massTotal" / "duration"?

I really appreciate your thoughts on this.

Regards,
MJ

Did you manage to answer this question?
FluentStarter is offline   Reply With Quote

Old   July 27, 2023, 03:44
Default How to change the working medium in sprayfoam from C7H16 to liquid NH3
  #17
New Member
 
mingyang
Join Date: Jul 2023
Posts: 1
Rep Power: 0
Mingyang is on a distinguished road
Dear Foamers,

In the calculation example of sprayFoam delivered with OpenFoam, the working medium is hydrocarbon C7H16. I want to use this calculation example to adapt the working medium to liquid NH3 for numerical simulation, observation and study of its spray effect. At present, the first step is to know how to convert the spray working medium in the tutorials.I uploaded my case files to the git hub,https://github.com/Fancy2121/sprayfo...m/aachenBomb,I would appreciate some advice.
Mingyang is offline   Reply With Quote

Reply

Tags
dieselfoam, injection temperature, sprayfoam


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
Mistype in Darcy-Forchheimer porosity model description makaveli_lcf OpenFOAM Programming & Development 6 December 11, 2014 17:07
Implementation: Break-up model of Coulaloglou and Tavlarides a.weber OpenFOAM Programming & Development 7 September 8, 2014 06:43
Implementation of the Low-Reynolds-number Chien model om3ro Main CFD Forum 0 March 24, 2010 19:16


All times are GMT -4. The time now is 18:12.