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

Les and diesel

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

Reply
 
LinkBack Thread Tools Display Modes
Old   February 12, 2009, 09:55
Default I would like to use a LES appr
  #1
New Member
 
Antonio Attili
Join Date: Mar 2009
Location: Rome, Italy
Posts: 21
Rep Power: 8
wizard1301 is on a distinguished road
I would like to use a LES approach instead of the RANS one to simulate the diesel spray


is there someone who knows where I could start?
wizard1301 is offline   Reply With Quote

Old   February 12, 2009, 11:03
Default You could start from the Xoodl
  #2
Member
 
ville vuorinen
Join Date: Mar 2009
Posts: 63
Rep Power: 8
ville is on a distinguished road
You could start from the Xoodles solver and add the dieselSpray stuff from that RANS solver to Xoodles.
All you need to do in principle is:

1) add relevant include files to Xoodles (that's basically the spray class)
2) add spray momentum source to UEqn.H
3) add the lines where you create the dieselSpray
4) take the input files from the dieselSpray (is that the name of the RANS solver!!? Don't remember,
something like that) tutorial where the inputfiles injectorProperties and sprayProperties are
5) I think you need to change the thermodynamic approach to hcombustionThermo (instead of hhu) for some reason.. BTW, COULD SOMEONE TELL WHY THIS NEEDED to be done?? Perhaps it was only an instruction I got some years ago since I didn't want to study combustion..
6) I think if you can figure out those things you
are nearly done - it shouldn't take long since it is mostly copy-pasting to Xoodles from an existing solver ..

7) And remember to modify the files in Make : take example from the dieselSpray tutorials.

8) the aachenbomb case gives you the input files.

Good luck!

-Ville
ville is offline   Reply With Quote

Old   February 12, 2009, 11:08
Default Hi, By the way, I have a ques
  #3
Member
 
ville vuorinen
Join Date: Mar 2009
Posts: 63
Rep Power: 8
ville is on a distinguished road
Hi,
By the way, I have a question about the thermodynamic approach in the main top level solver and then e.g. spray class: could someone
please explain about the thermodynamic approach the spray class expects? So does the spray class expect the actual solver to use the hcombustionthermo approach or is the class 'blind' to the thermodynamic approach otherwise than the information that goes to the constructor?

Best regards,
Ville
ville is offline   Reply With Quote

Old   February 13, 2009, 04:14
Default I guess Im the one to answer t
  #4
Super Moderator
 
niklas's Avatar
 
Niklas Nordin
Join Date: Mar 2009
Location: Stockholm, Sweden
Posts: 693
Rep Power: 19
niklas will become famous soon enough
I guess Im the one to answer this...

but your question is abit tricky and I think what you really is after is something else.

the spray class depends on that the gaseous properties are the same as those defined
in sprayThermoTypes.H.

The liquids have real multicomponent treatment, this also have some implications on the gaseous treatment. It needs to be a multicomponent mixture.
It should be pretty easy to use these to sum it up and collect them in a fuel source term though.
but you will still need the individual component concentrations to evaluate the evaporation rate.

There are a few places where it also assumes concentrations to be the same in the entire cell, so if you want to split the cells into burnt/unburned fractions, you need to rewrite some of the code. Mainly how it estimates evaporation times and saturation, so its nothing serious.

so I would say there is a slight thermo-dependency.
niklas is offline   Reply With Quote

Old   February 13, 2009, 04:50
Default Thank you very much, Would
  #5
Member
 
ville vuorinen
Join Date: Mar 2009
Posts: 63
Rep Power: 8
ville is on a distinguished road
Thank you very much,
Would it please be possible to still clarify with a few words a little bit more on my confusion - especially question 2) ..

1) a thermodynamic approach is chosen in the top-level code to define the gaseous phase thermodynamics -> I would now understand that depending on what kind of a spray I'm simulating, the thermo-dependency comes to more/less into play. Thus, if I just study a non-evaporating cloud, 'solid-particle' cold-injection type situation with no reactions etc, then the spray class will only need the knowledge of gas density and velocity which go to the constructor (I guess that's all the class needs to know about the flow field to update parcel positions).

2) I still don't quite get why the original hhuCombustionThermo from Xoodles needs to be changed to hCombustionThermo if Xoodles is modified to be a spray solver? Is there some restriction on using
e.g. the basicThermo approach (as in coodles) to
create a cold spray injection LES solver starting
from the coodles?

Best regards,
Ville
ville is offline   Reply With Quote

Old   February 13, 2009, 06:05
Default 1. yup. 2. hhu means that t
  #6
Super Moderator
 
niklas's Avatar
 
Niklas Nordin
Join Date: Mar 2009
Location: Stockholm, Sweden
Posts: 693
Rep Power: 19
niklas will become famous soon enough
1. yup.

2. hhu means that the thermodynamics is split into 2 fractions of burned and unburned gas, you also have 2 temperatures. There is no implementation of multicomponent mixtures for hhu.

I dont know what you mean with restriction, but if you look in the
src/thermophysicalModels/combustion/hCombustionThermo/hCombustionThermos.C
you will see all the available combinations and it should be possible, in theory, to modify the speciesProperties typedef in sprayThermoTypes.H to the type of thermo you would like to use.
niklas is offline   Reply With Quote

Old   February 13, 2009, 06:59
Default Right, think I got it, thanks
  #7
Member
 
ville vuorinen
Join Date: Mar 2009
Posts: 63
Rep Power: 8
ville is on a distinguished road
Right, think I got it, thanks once again!
-V
ville is offline   Reply With Quote

Old   June 23, 2009, 10:01
Default
  #8
Senior Member
 
xinguang cui
Join Date: Mar 2009
Posts: 116
Rep Power: 8
flying is on a distinguished road
I am also interested in the using LES approach instead of the RANS one to simulate the diesel spray.

But it is very sad that I am not sure about how to solve "Y" because I don't know the equation when I choose LES. I used the method as below to solve the Y, but I am not sure about it. could anyone make sure my method?

solve
(
fvm::ddt(rho, Yi)
+ mvConvection->fvmDiv(phi, Yi)
- fvm::laplacian(turbulence->muEff(), Yi)
==
dieselSpray.evaporationSource(i),
// + kappa*chemistry.RR(i),
mesh.solver("Yi")
);

In fact, There is no combustion in my case, so I didn't care about the chemistry and I delete "kappa*chemistry.RR(i)".If I care about the Chemistry, how do treat the term of "kappa*chemistry.RR(i)", which is generated in the RAS?
flying is offline   Reply With Quote

Old   June 23, 2009, 10:01
Default
  #9
Senior Member
 
xinguang cui
Join Date: Mar 2009
Posts: 116
Rep Power: 8
flying is on a distinguished road
Thanks and best wishes!
flying is offline   Reply With Quote

Old   July 9, 2009, 10:41
Question
  #10
Member
 
amin
Join Date: May 2009
Posts: 59
Rep Power: 8
az1362f is on a distinguished road
Hello Dear Ville;

I have some question about LES Spray simulation,please describe these points by more details:

first you said :"2) add spray momentum source to UEqn.H"

1-where is UEqn.h file? because there is no UEqn.h file in Xoodles directory. is it in Xifoam solver folder?

second you said:"4) take the input files from the dieselSpray (is that the name of the RANS solver!!? Don't remember,
something like that) tutorial where the inputfiles injectorProperties and sprayProperties are"

2-what do you mean? do you mean just copy and paste aachenbomb case in Xoodles tutorial folder?

third you said:"5) I think you need to change the thermodynamic approach to hcombustionThermo (instead of hhu) for some reason.. BTW, COULD SOMEONE TELL WHY THIS NEEDED to be done?? Perhaps it was only an instruction I got some years ago since I didn't want to study combustion.."

how can I change hhucobustionthermo to hcombustionthermo?
you mean simply change it in thermophysicalproperties file of acchenbomb case?
does it need to use chemkin folder in aachenbomb case for Xoodles solver?

fourth you said:"3) add the lines where you create the dieselSpray"


you mean just add these lines to Xoodles solver:

Info << "Evolving Spray" << endl;

dieselSpray.evolve();

or something else.
I try to make this solver but i face more problem.
Please help me.

Thanks
regard
az1362f is offline   Reply With Quote

Old   July 23, 2009, 16:46
Default
  #11
Senior Member
 
Sandy Lee
Join Date: Mar 2009
Posts: 207
Rep Power: 9
sandy is on a distinguished road
How to set BC about LES??
sandy is offline   Reply With Quote

Old   July 25, 2009, 06:06
Default
  #12
Senior Member
 
xinguang cui
Join Date: Mar 2009
Posts: 116
Rep Power: 8
flying is on a distinguished road
There is example in tutorial!
flying is offline   Reply With Quote

Old   July 25, 2009, 06:07
Default
  #13
Senior Member
 
xinguang cui
Join Date: Mar 2009
Posts: 116
Rep Power: 8
flying is on a distinguished road
Hello Dear Ville and az1362f:

Do you care about the question below?


I am also interested in the using LES approach instead of the RANS one to simulate the diesel spray.

But it is very sad that I am not sure about how to solve "Y" because I don't know the equation when I choose LES. I used the method as below to solve the Y, but I am not sure about it. could anyone make sure my method?

solve
(
fvm::ddt(rho, Yi)
+ mvConvection->fvmDiv(phi, Yi)
- fvm::laplacian(turbulence->muEff(), Yi)
==
dieselSpray.evaporationSource(i),
// + kappa*chemistry.RR(i),
mesh.solver("Yi")
);

In fact, There is no combustion in my case, so I didn't care about the chemistry and I delete "kappa*chemistry.RR(i)".If I care about the Chemistry, how do treat the term of "kappa*chemistry.RR(i)", which is generated in the RAS?

Thanks and best wishes!

Xinguang Cui
flying is offline   Reply With Quote

Old   July 31, 2009, 09:13
Default
  #14
Senior Member
 
Sandy Lee
Join Date: Mar 2009
Posts: 207
Rep Power: 9
sandy is on a distinguished road
Hi Cui, do you know how to specify the inlet value of k to oneEqEddy model?
sandy is offline   Reply With Quote

Old   July 31, 2009, 11:43
Default
  #15
Senior Member
 
xinguang cui
Join Date: Mar 2009
Posts: 116
Rep Power: 8
flying is on a distinguished road
I am really sorry. I don't know how to give it.
flying 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
DIESEL COMBUSTION Benzamia FLUENT 8 January 22, 2012 01:33
DPM for diesel sprays Cemil FLUENT 0 November 26, 2007 05:13
diesel spary adi FLUENT 0 March 29, 2006 04:42
diesel combustion Marco FLUENT 0 September 14, 2005 08:27
RNG k-e & Diesel simulation Ossi CD-adapco 0 August 10, 2001 12:25


All times are GMT -4. The time now is 15:27.