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

Start with DieselEngineFoam

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

Reply
 
LinkBack Thread Tools Display Modes
Old   April 10, 2006, 06:10
Default I tried to create a DieselEngi
  #1
Member
 
Tomislav Sencic
Join Date: Mar 2009
Posts: 42
Rep Power: 8
tsencic is on a distinguished road
I tried to create a DieselEngineFoam case, and it runs from CA=-180 to CA=-100 but then I obtain the usual: FOAM FATAL ERROR: attempt to use janafThermo <equation> out of temperature range 200->5000; T=5017

The mesh was created with k3prep (kiva preprocessor) and converted with kivaToFoam. It is a simple cylinder with polar coordinates (with central axis), 3359 cells, D=0.1m s=0.12m.

I tested it with checkMesh - everything looks o.k.

For the initial conditions I set up the data by combining the dieselFoam and engineFoam tutorials.

I have doubts about several imput data: in constant/injectionProperties, what are Cd and X?
How do I specify the start of injection?

I tried to play with deltaT, nCorectors, nNonOrtogonalCorectors, but it did not help.

I read some threads concerning the topic on this Message Board, but did not find the solution.

Where could I find the answers ( working tutorial, papers, commented code..) to start to use this solver?
tsencic is offline   Reply With Quote

Old   April 10, 2006, 13:51
Default are you using fixed temperatur
  #2
Member
 
stefan
Join Date: Mar 2009
Posts: 96
Rep Power: 8
stefanke is on a distinguished road
are you using fixed temperature bc ?

what is about the skewness and the non-orthg. ?
stefanke is offline   Reply With Quote

Old   April 11, 2006, 02:59
Default The temperature BC are fixedVa
  #3
Member
 
Tomislav Sencic
Join Date: Mar 2009
Posts: 42
Rep Power: 8
tsencic is on a distinguished road
The temperature BC are fixedValue uniform 373 (like for the kivaTest engineFoam tutorial).

The checkMesh gives me:
Non-ortogonality Max: 0.00127146 average: 0.00028. Non-ortogonality : OK
Max skewness = 72.9826 percent. Skewness OK
Mesh OK
tsencic is offline   Reply With Quote

Old   April 11, 2006, 03:04
Default try adabatic walls and tell me
  #4
Member
 
stefan
Join Date: Mar 2009
Posts: 96
Rep Power: 8
stefanke is on a distinguished road
try adabatic walls and tell me if the problem still exists.
stefanke is offline   Reply With Quote

Old   April 11, 2006, 05:11
Default It works with adiabatic walls
  #5
Member
 
Tomislav Sencic
Join Date: Mar 2009
Posts: 42
Rep Power: 8
tsencic is on a distinguished road
It works with adiabatic walls (zeroGradient)! Thank you for helping me.
I run it to ca=40, with injection and combustion. Then when I wanted to proceed from ca=40 it sais:
FOAM FATAL ERROR: attempt to use janafThermo <equation> out of temperature range 200->5000; T=5.5454e+212 (or such a value). What is the reason (and solution) for that?
I will try some other, more realistic meshes and injection setups with adiabatic walls.
tsencic is offline   Reply With Quote

Old   April 11, 2006, 07:42
Default I experienced the same problem
  #6
Member
 
stefan
Join Date: Mar 2009
Posts: 96
Rep Power: 8
stefanke is on a distinguished road
I experienced the same problem with fixed temperature bc (see bugs). I hope there is someone (perhaps Niklas) who can solve this problem. I think this is a bug of the dieselEngineFoam solver!
stefanke is offline   Reply With Quote

Old   April 28, 2006, 02:44
Default I did some progress but i stil
  #7
Member
 
Tomislav Sencic
Join Date: Mar 2009
Posts: 42
Rep Power: 8
tsencic is on a distinguished road
I did some progress but i still have many questions.

How can I set up a multiple holes injector or multiple injectors?

In constant/injectorProperties:
-what is X? I found something like molar fraction. Which moles in which?
-what is the unit of the massFlowRateProfile list entries?
tsencic is offline   Reply With Quote

Old   April 28, 2006, 04:53
Default >What you usually do is to use
  #8
Member
 
stefan
Join Date: Mar 2009
Posts: 96
Rep Power: 8
stefanke is on a distinguished road
>What you usually do is to use a sector mesh (i.e. 8 injetor holes -> 1/8 sector mesh) and apply cyclic boundary conditions.

> X defines the molecular concentration
> The massFlowRateProfile is a dimensionless profile. It only defines how the profile looks like not the total injected mass!
stefanke is offline   Reply With Quote

Old   April 28, 2006, 07:30
Default >> How can I set up a multiple
  #9
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
>> How can I set up a multiple holes injector or multiple injectors?
injectorProperties is just a list of injectors so you can add as many as you want.
(
{
injector1....
}
{
injector2....
}
)

X is the liquid volume fraction of the
species defined in liquidFuelComponents in
thermophysicalProperties.

the unit is dimensionless, it is scaled so that the integrated profile corresponds to the injected mass.

N
niklas is offline   Reply With Quote

Old   July 13, 2006, 06:10
Default I did some simulations with di
  #10
Member
 
Tomislav Sencic
Join Date: Mar 2009
Posts: 42
Rep Power: 8
tsencic is on a distinguished road
I did some simulations with dieselEngineFoam. The results (total pressure and temperature) are not bad, but I still have many doubts and things to improove.

1. I tried to compare total dQ. I calculated it by dQ.weighedAverage(mesh.V()).value()
The shape is nice, but the value (I supose the unit is J/s) is about 50 x higher than the one obtained with VIBE function in qD simulations. Is it another unit (which?) or I am doing something wrong?

2. What is the difference between atomisation (atomisationModel) and breakup (braeakupModel)? I thaught it was a different word for the same thing, and that there was the difference between primary and secondary breakup/atomisatin.

3. What is the difference between injectorTypes unitInjector and commonRailInjector?

4. The sprayProperties. Which injectorModel is appropriate to simulate the injection (unit injector or common rail) in a big diesel engine (bore x stroke = 0.23 x 0.25 m). And the other sprayProperties (atomisationModel, breakupModel,...)? It would be easiest if someone would show me his sprayProperties file of a succesful simulation.
tsencic is offline   Reply With Quote

Old   July 13, 2006, 06:37
Default 1. I don't know what you are d
  #11
Member
 
chris book
Join Date: Mar 2009
Posts: 85
Rep Power: 8
chris1980 is on a distinguished road
1. I don't know what you are doing in detail. Without providing a detailed information about your calculation of the RoHR we are not able to help you.

To something like this:

loop over all species and cells and use

dQ[celli] -= chemistry.specieThermo()[i].H(T[celli])*RRi[celli]*mesh.V()[celli];

to calculte the rate of heate release due to the chemical reactions.

2.
atomisation -> primary break up
break up -> secondary break up

3. Have look at the different properties in the injector directory

4. The choice of spay sub models depends on the application. In your case I would prefer a commonRailInjector, instead of using a atomisationModel define an initial doplet size distribution and KHRT breakupModel.

hth
chris1980 is offline   Reply With Quote

Old   August 30, 2006, 04:11
Default I expected the same "strange"
  #12
Member
 
chris book
Join Date: Mar 2009
Posts: 85
Rep Power: 8
chris1980 is on a distinguished road
I expected the same "strange" behaviour but I cannot find a solution to get rid of it.

This problem breaks my head
chris1980 is offline   Reply With Quote

Old   August 30, 2006, 04:23
Default Yupp, I saw the same problem s
  #13
Member
 
stefan
Join Date: Mar 2009
Posts: 96
Rep Power: 8
stefanke is on a distinguished road
Yupp, I saw the same problem some time ago but I don't find the time to investigate the problem further.

In my case the aspect ratio of the squish was very poor because the cells are squeezed from BTC to TDC. This may impact the spray so this strange behaviour shows up.

Fell free to test some different mesh setups (lower the layers in the squish => better aspect ratio at TDC) to see if my guess becomes true.
stefanke is offline   Reply With Quote

Old   September 8, 2006, 06:10
Default Any news releated to this stra
  #14
Member
 
chris book
Join Date: Mar 2009
Posts: 85
Rep Power: 8
chris1980 is on a distinguished road
Any news releated to this strange spray behaviour?
chris1980 is offline   Reply With Quote

Old   September 8, 2006, 08:35
Default I have also experienced this,
  #15
New Member
 
Rikard Mikalsen
Join Date: Mar 2009
Posts: 10
Rep Power: 8
kukikano is on a distinguished road
I have also experienced this, and I got better behaviour when I moved the injector further into the cylinder (i.e. below the highly compressed mesh in the liner region and away from the cylinderHead wall). I figured it was either the poor aspect ratio that caused the problems, like Stefan says, or the fact that the fuel spray was too close to the wall.

I didn't investigate it further, though.
kukikano is offline   Reply With Quote

Old   November 15, 2006, 12:17
Default Hi, How could I define an inj
  #16
Member
 
ville vuorinen
Join Date: Mar 2009
Posts: 63
Rep Power: 8
ville is on a distinguished road
Hi,
How could I define an injector with multiple holes?
As mentioned:
"injectorProperties is just a list of injectors so you can add as many as you want.
(
{
injector1....
}
{
injector2....
}
) "

But where do I manifest that I want to 'repeat'
the single hole injector N times and put these
injectors at given positions pointing at given
directions?

-Ville
ville is offline   Reply With Quote

Old   April 27, 2007, 09:31
Default Hi all, I'm trying to defin
  #17
New Member
 
Daniele Ettorre
Join Date: Mar 2009
Location: Milan, Italy
Posts: 1
Rep Power: 0
daniele is on a distinguished road
Hi all,

I'm trying to define an injector with multiple holes
but I have thi s ERROR:

/*---------------------------------------------------------------------------*\
| ========= | |
| \ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \ / O peration | Version: 1.3.6 |
| \ / A nd | Web: http://www.openfoam.org |
| \/ M anipulation | |
\*---------------------------------------------------------------------------*/

Exec : dieselFoam . sprayTot
Date : Apr 26 2007
Time : 12:03:33
Host : klunk
PID : 13683
Root : /locals/klunkLoc/sprayCNR25
Case : sprayTot
Nprocs : 1
Create time

Create mesh for time = 0.0019835


Reading thermophysicalProperties
Selecting thermodynamics package hMixtureThermo<reactingmixture>
Selecting chemistryReader chemkinReader
Reading field U

Reading/calculating face flux field phi

Creating turbulence model.

Selecting turbulence model kEpsilon
Creating field DpDt

Constructing chemical mechanism
Selecting ODE solver SIBS
chemistryModel::chemistryModel: Number of species = 5 and reactions = 1
Selecting sootModel noSootModel

Reading environmentalProperties
Reading combustion properties

Constructing Spray
Selecting injectorType commonRailInjector
injectionPressureProfile_.size() = 2, massFlowRateProfile_.size() = 12
end constructor. in commonRail
Selecting atomizationModel off
Selecting dragModel standardDragModel
Selecting evaporationModel saturateEvaporationModel
Selecting heatTransferModel RanzMarshall
Selecting wallModel remove
Selecting breakupModel ReitzKHRT
Selecting collisionModel off
Selecting dispersionModel off
Selecting injectorModel hollowConeInjector
Selecting pdfType RosinRammler



--> FOAM FATAL ERROR : object is not allocated

From function
autoPtr<t>::operator->()
in file /gamma/dettorre/OpenFOAM/OpenFOAM-1.3.6/src/OpenFOAM/lnInclude/autoPtrI.H at line 165.

FOAM aborting


I put some info and I found the error in the file:

/gamma/dettorre/OpenFOAM/OpenFOAM-1.3.6/src/lagrangian/dieselSpray/spraySubModel s/injectorModel/hollowCone/hollowCone.C

at line

vector dir = sm.injectors()[i].properties()->direction();

but I can't understand what is the problem.
Can somebody help me?
daniele is offline   Reply With Quote

Old   April 27, 2007, 09:36
Default probably there is something wr
  #18
Member
 
stefan
Join Date: Mar 2009
Posts: 96
Rep Power: 8
stefanke is on a distinguished road
probably there is something wrong in your injector dictionary. Please post the latter.
stefanke is offline   Reply With Quote

Old   June 28, 2007, 01:12
Default Hi. I run dieselEngineFoam of
  #19
nishio
Guest
 
Posts: n/a
Hi.
I run dieselEngineFoam of tutorials.
Automatically,injection start 0deg. but I want to change injection timing .
Which file have this parameter??
Please tell me.
Thanks.
  Reply With Quote

Old   June 28, 2007, 03:16
Default You have to modify the massFlo
  #20
Member
 
Tomislav Sencic
Join Date: Mar 2009
Posts: 42
Rep Power: 8
tsencic is on a distinguished road
You have to modify the massFlowRateProfile in constat/injectorProperties
tsencic 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
DieselEngineFoam stefanke OpenFOAM Pre-Processing 42 December 3, 2008 23:53
Disable injection in dieselEngineFoam francesco OpenFOAM Running, Solving & CFD 0 November 26, 2008 02:23
Strange pressure with dieselEngineFoam tsencic OpenFOAM Bugs 1 December 12, 2007 05:39
Needing some advise about dieselEngineFoam adorean OpenFOAM Running, Solving & CFD 33 November 20, 2007 00:27
Cannot start using CFX-5.7.1 Atit CFX 8 February 16, 2005 07:40


All times are GMT -4. The time now is 03:45.