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

eddy dissipation concept + infinitely fast chemistry

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

Reply
 
LinkBack Thread Tools Display Modes
Old   July 20, 2010, 10:00
Default eddy dissipation concept + infinitely fast chemistry
  #1
Senior Member
 
Aram Amouzandeh
Join Date: Mar 2009
Location: Vienna, Vienna, Austria
Posts: 184
Rep Power: 8
mabinty is on a distinguished road
dear all!

i implemented the eddy dissipation concept combustion model with infinitely fast chemistry according to the paper of Magnussen: "The Eddy Dissipation Concept A Bridge Between Science And Technology", Lisbon, 2005, based on fireFoam, using a k-eps turbulence model. the code runs and the flow field looks plausible at first sight (in 3D). but when integrating the reaction rate (or the heat release due to burnt fuel) over the whole computational domain, I can see that the reaction rate is overestimated, i.e. more fuel burns as is injected. I tried to increase spatial resolution as well as to change the value for the turbulent dissipation rate at fuel inlet but always observe the same behavior. Interestingly the reaction rate matches quiet well with the injected amount of fuel for the 2D tutorial case applying k-eps turbulence model.

Anybody experienced similar problems??

I would greatly appreciate your comments!!

Regards,
Aram
mabinty is offline   Reply With Quote

Old   September 23, 2010, 05:27
Default
  #2
Senior Member
 
Aram Amouzandeh
Join Date: Mar 2009
Location: Vienna, Vienna, Austria
Posts: 184
Rep Power: 8
mabinty is on a distinguished road
dear all!

i figured out that the problem is not the combustion model. additional to the convective input of fuel (prescribed by mass flow), fuel enters the domain by diffusion due to the sharp gradient of fuel mass fraction fu at inlet. first i tried to add a long fuel inlet pipe but as a result of low inlet velocities the combustion process "travels" back to the inlet causing the same problem as before. Hence, i decided to set the diffusion coefficient alphaEff in the transport equation of fu (fuhsEqn.H) to zero at inlet, what solved the problem:

Code:
    alphaEffZeroDiff = turbulence->alphaEff(); 
    label inletPatchID = mesh.boundaryMesh().findPatchID("inlet");
    alphaEffZeroDiff.boundaryField()[inletPatchID] -= 
        alphaEffZeroDiff.boundaryField()[inletPatchID];
    fvScalarMatrix R = combustion->R(fu);

    {
        fvScalarMatrix fuEqn
        (
            fvm::ddt(rho, fu)
          + mvConvection->fvmDiv(phi, fu)
          - fvm::laplacian(alphaEffZeroDiff, fu)
         ==
            R
        );

        fuEqn.relax();
        fuEqn.solve();
    }
I m sure there is a nicer way to do that, e.g. through a boundary condition where the diffusion coeff. at a certain patch for the considered variable is set to zero. I tried to create a "fixedValueZeroDiffusivity" BC as follows; in fixedValueZeroDiffusivityFvPatchScalarField::updat eCoeffs():

Code:
    const label patchI = patch().index();

    const RASModel& rasModel = db().lookupObject<RASModel>("RASProperties");

    scalarField alphaEffp = rasModel.alphaEff()().boundaryField()[patchI];
    
    rasModel.alphaEff()().boundaryField()[patchI] -= alphaEffp;

    fixedValueFvPatchScalarField::updateCoeffs();
To allow alphaEff to be overwritten at the boundary i added a non-constant function of alphaEff() in kEpsilon.H. When analysing the results i can see that alphaEff is not changed. Any ideas?? I greatly appreciate your comments.

Thx in advance,
Aram
mabinty is offline   Reply With Quote

Old   September 12, 2012, 14:01
Default
  #3
New Member
 
Michael
Join Date: Dec 2011
Location: Geneva
Posts: 28
Rep Power: 5
lithos is on a distinguished road
Dear mabinty,

I am working with the solver since a few month and I discovered the problem very lately. So I will see to confirm it as soon as possible. In fact, I still haven't found a possibility to compute the diffusive flux.

It looks like this is still unsolved: http://lists.cfd-online.com/pipermai...ne/004327.html

also:

http://www.extend-project.de/user-fo...library-update

Cheers,

Michael

Last edited by lithos; September 13, 2012 at 10:01.
lithos is offline   Reply With Quote

Old   September 13, 2012, 04:39
Default
  #4
New Member
 
Michael
Join Date: Dec 2011
Location: Geneva
Posts: 28
Rep Power: 5
lithos is on a distinguished road
Sorry for the double post due to a missing reload, questions already raised above (Apart from deactivating completely the diffusion, is there any way to compute the diffusive mass flow through a patch?).
lithos is offline   Reply With Quote

Old   September 17, 2012, 04:55
Default
  #5
Senior Member
 
Aram Amouzandeh
Join Date: Mar 2009
Location: Vienna, Vienna, Austria
Posts: 184
Rep Power: 8
mabinty is on a distinguished road
hey Michael,

my problem was not to compute diffusive fluxes. for a variable phi and a diffusion coefficient alphaEff i d try it like that

Code:
volScalarField dDiff = 
        - fvc::laplacian(turbulence->alphaEff(), phi);
Instead of setting diffusion to zero a boundary condition which balances convective and diffusive fluxes in order to prescribe a total in/outflow of a specific variable at a certain patch solved my problem (see attachment).

cheers,
aram
Attached Files
File Type: gz turbulentTotalFlowRateAdvDiff.gz (2.6 KB, 36 views)
mabinty 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
Eddy dissipation concept thanz OpenFOAM Running, Solving & CFD 1 March 26, 2010 04:32
Eddy Dissipation Concept adrian CFX 0 October 18, 2007 03:57
Eddy Dissipation Concept adrian FLUENT 0 October 18, 2007 03:55
Eddy Dissipation Concept adrian Main CFD Forum 0 October 18, 2007 03:54
eddy dissipation concept sisilxp FLUENT 5 May 13, 2004 05:13


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