CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM News & Announcements > OpenFOAM Announcements from Other Sources

A New Solver for Supersonic Combustion

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

Like Tree4Likes
  • 2 Post By nakul
  • 1 Post By nakul
  • 1 Post By hg2lf

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   December 7, 2010, 15:08
Default A New Solver for Supersonic Combustion
  #1
Senior Member
 
Nakul
Join Date: Apr 2010
Location: India
Posts: 147
Rep Power: 15
nakul is on a distinguished road
Hi,

I have modifed the rhoCentralFOAM solver for viscous supersonic flows to incorporate combustion and mixing. This solver is capable of simulating supersonic combustion. The combustion code has been taken from reactingFOAM. I have named it "rhoCentralSpeciesFOAM".

I have done a test case simulation with this and it works. But there are still few issues with the solver. The solver is not very stable. It is still quite sensitive to test case setup. Although the solver works but the results are not very accurate as the chemical reactions become more complex.

I request the OF community to have look at this solver. Any suggestions regarding any modifications are highly appreciated. Steps needed to integerate it with your OF have been specified in the following post.

-Nakul
Attached Files
File Type: gz rhoCentralSpeciesFoam.tar.gz (5.4 KB, 124 views)
ranjansm and 13722617 like this.
nakul is offline   Reply With Quote

Old   December 8, 2010, 09:42
Default
  #2
Senior Member
 
Holger Marschall
Join Date: Mar 2009
Location: Darmstadt, Germany
Posts: 125
Rep Power: 19
holger_marschall is on a distinguished road
Send a message via Skype™ to holger_marschall
Hi Nakul,

I encourage you to found a working group and a sub-project on this topic on the Extend-Project portal (www.extend-project.de). There you'll find all infrastructure to organize the development. E.g.
  • group discussion board (+group mailing)
  • group file sharing
  • project management (task management, messaging)
best regards,
Holger
__________________
Holger Marschall
web: http://www.holger-marschall.info
mail: holgermarschall@yahoo.de
holger_marschall is offline   Reply With Quote

Old   June 4, 2012, 17:07
Default Any progress on the supersonc combustion sover
  #3
New Member
 
Ranjan S. Mehta
Join Date: Apr 2010
Location: Huntsville, Alabama
Posts: 6
Rep Power: 15
ranjansm is on a distinguished road
Hi Nakul

I was wondering if you have developed the supersonic combustion solver any further, or put it in the openfoam-extend or anything like that ?

Thanks
Ranjan
ranjansm is offline   Reply With Quote

Old   June 4, 2012, 17:24
Default
  #4
New Member
 
Ranjan S. Mehta
Join Date: Apr 2010
Location: Huntsville, Alabama
Posts: 6
Rep Power: 15
ranjansm is on a distinguished road
Quote:
Originally Posted by nakul View Post
Hi,

I have modifed the rhoCentralFOAM solver for viscous supersonic flows to incorporate combustion and mixing. This solver is capable of simulating supersonic combustion. The combustion code has been taken from reactingFOAM. I have named it "rhoCentralSpeciesFOAM".

I have done a test case simulation with this and it works. But there are still few issues with the solver.

-Nakul
Nakul,

Do you happen to have the test case that you ran for this solver at all ?
Thanks
Ranjan
ranjansm is offline   Reply With Quote

Old   June 5, 2012, 02:42
Default
  #5
Senior Member
 
Nakul
Join Date: Apr 2010
Location: India
Posts: 147
Rep Power: 15
nakul is on a distinguished road
Hi,

No I didn't find time to work on it any further but I am interested in its further development.

I have run a test case using this solver. You may tell me what you expect to do with the solver and I would see if I could help with that.

-Nakul
nakul is offline   Reply With Quote

Old   June 6, 2012, 12:07
Default test case for rhoCentralSpeciesFoam
  #6
New Member
 
Ranjan S. Mehta
Join Date: Apr 2010
Location: Huntsville, Alabama
Posts: 6
Rep Power: 15
ranjansm is on a distinguished road
Hi Nakul

I may have some dedicated time to work with you and get this rhoCentralSpeciesFoam in order.

But if you can give me the test case that you used, it will help me greatly.

I compiled your solver with OF-1.7.1 and ran a case, which I have, but I am sure my BCs are a big problem,

Quickly go temperature out of bounds etc.

We can work together, and as you have been suggested already, that it may be a good idea to put it in a open repository so multiple users can suggest changes etc.

Let me know!!
Thanks
Ranjan
ranjansm is offline   Reply With Quote

Old   June 7, 2012, 02:45
Default
  #7
Senior Member
 
Nakul
Join Date: Apr 2010
Location: India
Posts: 147
Rep Power: 15
nakul is on a distinguished road
Hi Ranjan

Please tell me your mail-id and I will mail you the test case by Saturday. We can discuss its setup after that. But I would advice that if you have dedicated time to work on supersonic combustion, then create a test case on your own.

Besides I don't have all the validation data for my test case so it would be a good idea if you prepare a case with large amount of validation data readily available. We can then discuss how to solve that problem. Also you may tell me what exactly are you planning to study in supersonic combustion because there are quite a few problems that one can solve in this area.

-Nakul
qinfei likes this.
nakul is offline   Reply With Quote

Old   August 3, 2012, 12:45
Default
  #8
New Member
 
qf
Join Date: Sep 2011
Posts: 3
Rep Power: 14
qinfei is on a distinguished road
Hi Nakul and ranjansm

I interest in simulation of combustion in high speed flow, I think rhoCentralSpeciesFoam will be very helpful to me. Could you give me the test case that you used, it will help me greatly. Then I will post the results of the simulation, and may develop it further. Thanks in advance.
qinfei is offline   Reply With Quote

Old   August 4, 2012, 17:06
Default solver
  #9
New Member
 
Ranjan S. Mehta
Join Date: Apr 2010
Location: Huntsville, Alabama
Posts: 6
Rep Power: 15
ranjansm is on a distinguished road
Hi qinfei,
I am working iwth Nakul and one more engineer here.
I think we have ironed out most of the issues. There are some good validation cases available too.
As soon as things look good enough to getting out, I will let you know
Thanks!
Ranjan

Quote:
Originally Posted by qinfei View Post
Hi Nakul and ranjansm

I interest in simulation of combustion in high speed flow, I think rhoCentralSpeciesFoam will be very helpful to me. Could you give me the test case that you used, it will help me greatly. Then I will post the results of the simulation, and may develop it further. Thanks in advance.
ranjansm is offline   Reply With Quote

Old   August 5, 2012, 07:26
Default
  #10
New Member
 
qf
Join Date: Sep 2011
Posts: 3
Rep Power: 14
qinfei is on a distinguished road
Hi ranjansm,

Thank your reply, I am waiting for your results, good luck!
qinfei is offline   Reply With Quote

Old   August 30, 2012, 19:53
Default
  #11
Member
 
Chris
Join Date: Aug 2012
Location: Calgary, Alberta, Canada
Posts: 77
Rep Power: 13
ChrisA is on a distinguished road
Hi Nakul, Ranjan,

Your work sounds like the exact thing I'm working on over here, implementation of combustion/reactions from reactionFoam in our modified rhocentralfoam solver. Where are you guys at with the solver? Would it be possible to get one of your test cases?
ChrisA is offline   Reply With Quote

Old   December 27, 2012, 03:30
Smile sound speed calculated not properly
  #12
New Member
 
何刚
Join Date: Jan 2011
Posts: 25
Rep Power: 15
hg2lf is on a distinguished road
When you calculate the sound speed c,you do like this:
c=sqrt(thermo.Cp()/(hermo.Cp()-R)*rPsi),R is the commom specific heat : 8.314j/mol.K,You should do like this:
c=sqrt(thermo.Cp()/(thermo.Cp()-rPsi/T)*rPsi).
Best regards.
tatu likes this.

Last edited by hg2lf; December 27, 2012 at 03:30. Reason: wrong
hg2lf is offline   Reply With Quote

Old   December 29, 2012, 01:14
Default something wrong
  #13
New Member
 
何刚
Join Date: Jan 2011
Posts: 25
Rep Power: 15
hg2lf is on a distinguished road
when I modified the calculation of sound speed c, I use rhoCentralSpeciesFoam to simulate a cold flow in a laval nozzle,but the temperature field seems to be wrong. H =hs+1/2*U^2 does not keep constant. Somebody who research more deeply into the solver?
hg2lf is offline   Reply With Quote

Old   February 13, 2013, 10:59
Default
  #14
New Member
 
Tatu Pinomaa
Join Date: Oct 2012
Location: Finland
Posts: 16
Rep Power: 13
tatu is on a distinguished road
Quote:
Originally Posted by hg2lf View Post
when I modified the calculation of sound speed c, I use rhoCentralSpeciesFoam to simulate a cold flow in a laval nozzle,but the temperature field seems to be wrong. H =hs+1/2*U^2 does not keep constant. Somebody who research more deeply into the solver?
Nakul's phiHp is equal to the one with internal energy, i.e. it is
Code:
    surfaceScalarField phiEp
        (
            aphiv_pos*(rho_pos*(hs_pos + 0.5*magSqr(U_pos)) + p_pos)
          + aphiv_neg*(rho_neg*(hs_neg + 0.5*magSqr(U_neg)) + p_neg)
          + aSf*p_pos - aSf*p_neg
        );
Shouldn't this equation look different if we replace the pressure term with internal energy (rho*e = rho*hs - p). Could removing the first two pressure terms work? In other words:
Code:
    surfaceScalarField phiHp
        (
            aphiv_pos*(rho_pos*(hs_pos + 0.5*magSqr(U_pos)) )
          + aphiv_neg*(rho_neg*(hs_neg + 0.5*magSqr(U_neg)) )
          + aSf*p_pos - aSf*p_neg
        );
This seems to bring my underestimated temperature field closer to the value it should be in my case (according to reactingFoam and experiments).

There's also something fishy going on with the diffusive energy corrector equation. I get very promising temperature values if I comment out the last two terms:
Code:
fvm::ddt(rho, hs) - fvc::ddt(rho, hs)
 - fvm::laplacian(turbulence->alphaEff(), hs) // alphaEff = alpha + alphat 
//+ fvc::laplacian(turbulence->alpha(), hs)  
//- fvc::laplacian(k, T)   // turbulent + laminar energy diffusion via T
I guess the last two terms are there to drive energy and temperature consistent, i.e. hs --> CpT + p/rho + magSqr(U). I still don't fully understand the two last terms, e.g. why k is laminar + turbulent, but alpha is laminar?



Tatu

Last edited by tatu; February 19, 2013 at 08:37.
tatu is offline   Reply With Quote

Old   February 13, 2013, 13:12
Default
  #15
Member
 
Chris
Join Date: Aug 2012
Location: Calgary, Alberta, Canada
Posts: 77
Rep Power: 13
ChrisA is on a distinguished road
It's worth noting that the thermal conductivity (k) here is calculated using a prandlt number equal to 1 (or constant, depending on if they fixed the bug where it was hard coded equal to 1) this could be a source of discrepancy, especially in the boundary layer. Anyway, sutherland is able to calculate k (kappa in sutherland), giving it proper variability with respect to temperature. If you're using a different viscosity model, it should have an equation to calculate k. There just need to be functions implemented to calculate it for the domain and pass it up to the solver (you can copy the way Cp is written for the relevant libraries).

Also, with regards to your question, I don't believe k has a turbulent value as that term governs the energy transfer due to conduction, which shouldn't change as a result of turbulent flow (directly anyway).

The alpha bits I'm a bit more foggy on, my focus case is laminar so I haven't checked out the turbulence much.

It's been awhile since I've checked this stuff in the libraries, feel free to correct anything you think might be wrong.

Last edited by ChrisA; February 13, 2013 at 13:58.
ChrisA is offline   Reply With Quote

Old   February 19, 2013, 07:44
Default The matter with p
  #16
New Member
 
何刚
Join Date: Jan 2011
Posts: 25
Rep Power: 15
hg2lf is on a distinguished road
Hi,Tatu,

Thanks for your reply. But don't you think the phiHp has nothing to do with p. I think the phiHp should be like this:

surfaceScalarField phiHp
(
aphiv_pos*(rho_pos*(hs_pos + 0.5*magSqr(U_pos)) )
+ aphiv_neg*(rho_neg*(hs_neg + 0.5*magSqr(U_neg)) )
);
so ,what do you think of this?
hg2lf is offline   Reply With Quote

Old   February 19, 2013, 08:34
Default
  #17
New Member
 
Tatu Pinomaa
Join Date: Oct 2012
Location: Finland
Posts: 16
Rep Power: 13
tatu is on a distinguished road
Quote:
Originally Posted by hg2lf View Post
Hi,Tatu,

Thanks for your reply. But don't you think the phiHp has nothing to do with p. I think the phiHp should be like this:

surfaceScalarField phiHp
(
aphiv_pos*(rho_pos*(hs_pos + 0.5*magSqr(U_pos)) )
+ aphiv_neg*(rho_neg*(hs_neg + 0.5*magSqr(U_neg)) )
);
so ,what do you think of this?
I'm not completely sure what is the meaning of the term aSf*p_pos - aSf*p_neg (anyone?). It could be there to capture "propagation of the stress" through cells due to compressibility effects. Anyway, I think it should be left in the equation. Can anyone confirm this?

Tatu
tatu is offline   Reply With Quote

Old   February 19, 2013, 08:47
Default
  #18
New Member
 
何刚
Join Date: Jan 2011
Posts: 25
Rep Power: 15
hg2lf is on a distinguished road
When you compare the equation in the solver and the governing equation,you could find that the phiHp should not inclue p. Maybe the term aSf*p_pos - aSf*p_neg has something to do with the scheme of rho central upwind ?
hg2lf is offline   Reply With Quote

Old   February 19, 2013, 08:48
Default
  #19
New Member
 
Tatu Pinomaa
Join Date: Oct 2012
Location: Finland
Posts: 16
Rep Power: 13
tatu is on a distinguished road
Quote:
Originally Posted by hg2lf View Post
When you compare the equation in the solver and the governing equation,you could find that the phiHp should not inclue p. Maybe the term aSf*p_pos - aSf*p_neg has something to do with the scheme of rho central upwind ?
Yep, that's what I meant.

Tatu
tatu is offline   Reply With Quote

Old   February 27, 2024, 10:44
Default
  #20
Senior Member
 
Dongyue Li
Join Date: Jun 2012
Location: Beijing, China
Posts: 838
Rep Power: 17
sharonyue is on a distinguished road
`aSf*p_pos - aSf*p_neg` is not related to the scheme. In the following code, p is simply advected by 0.5*(phi+p+ + phi-p-) instead of 0.5*(phi(KT scheme)p+ + phi(KT scheme)p-). Written in this way, they offset with each.

Oh you can think it in this way, if you neglect `aSf*p_pos - aSf*p_neg`, KT scheme would be used to discretize \nabla\cdot(Up), but it should not.

See Eqn 27 in http://dyfluid.com/rhoCentralFoam.html

```
surfaceScalarField phiEp
(
aphiv_pos*(rho_pos*(hs_pos + 0.5*magSqr(U_pos)) + p_pos)
+ aphiv_neg*(rho_neg*(hs_neg + 0.5*magSqr(U_neg)) + p_neg)
+ aSf*p_pos - aSf*p_neg
);
```
__________________
My OpenFOAM algorithm website: http://dyfluid.com
By far the largest Chinese CFD-based forum: http://www.cfd-china.com/category/6/openfoam
We provide lots of clusters to Chinese customers, and we are considering to do business overseas: http://dyfluid.com/DMCmodel.html
sharonyue 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
Working directory via command line Luiz CFX 4 March 6, 2011 21:02
premixed combustion - parallel solver error Craig FLUENT 0 October 14, 2008 16:07
why the solver reject it? Anyone with experience? bearcat CFX 6 April 28, 2008 15:08
Hydrogen Air combustion in a combustion chamber popi CFX 7 July 11, 2007 19:40
Need Similar Combustion Tutorial as my model Arnab Siemens 0 May 24, 2005 19:23


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