
[Sponsors] 
December 7, 2010, 14:08 
A New Solver for Supersonic Combustion

#1 
Senior Member
Nakul
Join Date: Apr 2010
Location: India
Posts: 147
Rep Power: 16 
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 

December 8, 2010, 08:42 

#2 
Senior Member

Hi Nakul,
I encourage you to found a working group and a subproject on this topic on the ExtendProject portal (www.extendproject.de). There you'll find all infrastructure to organize the development. E.g.
Holger
__________________
Holger Marschall web: http://www.holgermarschall.info mail: holgermarschall@yahoo.de 

June 4, 2012, 16:07 
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 
Hi Nakul
I was wondering if you have developed the supersonic combustion solver any further, or put it in the openfoamextend or anything like that ? Thanks Ranjan 

June 4, 2012, 16:24 

#4  
New Member
Ranjan S. Mehta
Join Date: Apr 2010
Location: Huntsville, Alabama
Posts: 6
Rep Power: 15 
Quote:
Do you happen to have the test case that you ran for this solver at all ? Thanks Ranjan 

June 5, 2012, 01:42 

#5 
Senior Member
Nakul
Join Date: Apr 2010
Location: India
Posts: 147
Rep Power: 16 
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 

June 6, 2012, 11:07 
test case for rhoCentralSpeciesFoam

#6 
New Member
Ranjan S. Mehta
Join Date: Apr 2010
Location: Huntsville, Alabama
Posts: 6
Rep Power: 15 
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 OF1.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 

June 7, 2012, 01:45 

#7 
Senior Member
Nakul
Join Date: Apr 2010
Location: India
Posts: 147
Rep Power: 16 
Hi Ranjan
Please tell me your mailid 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 

August 3, 2012, 11:45 

#8 
New Member
qf
Join Date: Sep 2011
Posts: 3
Rep Power: 14 
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. 

August 4, 2012, 16:06 
solver

#9  
New Member
Ranjan S. Mehta
Join Date: Apr 2010
Location: Huntsville, Alabama
Posts: 6
Rep Power: 15 
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:


August 5, 2012, 06:26 

#10 
New Member
qf
Join Date: Sep 2011
Posts: 3
Rep Power: 14 
Hi ranjansm,
Thank your reply, I am waiting for your results, good luck! 

August 30, 2012, 18:53 

#11 
Member
Chris
Join Date: Aug 2012
Location: Calgary, Alberta, Canada
Posts: 77
Rep Power: 13 
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? 

December 27, 2012, 02:30 
sound speed calculated not properly

#12 
New Member
何刚
Join Date: Jan 2011
Posts: 25
Rep Power: 15 
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. Last edited by hg2lf; December 27, 2012 at 02:30. Reason: wrong 

December 29, 2012, 00:14 
something wrong

#13 
New Member
何刚
Join Date: Jan 2011
Posts: 25
Rep Power: 15 
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?


February 13, 2013, 09:59 

#14  
New Member
Tatu Pinomaa
Join Date: Oct 2012
Location: Finland
Posts: 16
Rep Power: 13 
Quote:
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 ); 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 ); 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 Tatu Last edited by tatu; February 19, 2013 at 07:37. 

February 13, 2013, 12:12 

#15 
Member
Chris
Join Date: Aug 2012
Location: Calgary, Alberta, Canada
Posts: 77
Rep Power: 13 
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 12:58. 

February 19, 2013, 06:44 
The matter with p

#16 
New Member
何刚
Join Date: Jan 2011
Posts: 25
Rep Power: 15 
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? 

February 19, 2013, 07:34 

#17  
New Member
Tatu Pinomaa
Join Date: Oct 2012
Location: Finland
Posts: 16
Rep Power: 13 
Quote:
Tatu 

February 19, 2013, 07:47 

#18 
New Member
何刚
Join Date: Jan 2011
Posts: 25
Rep Power: 15 
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 ?


February 19, 2013, 07:48 

#19 
New Member
Tatu Pinomaa
Join Date: Oct 2012
Location: Finland
Posts: 16
Rep Power: 13 

February 27, 2024, 09:44 

#20 
Senior Member
Dongyue Li
Join Date: Jun 2012
Location: Beijing, China
Posts: 838
Rep Power: 17 
`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+ + phip) 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 CFDbased forum: http://www.cfdchina.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 

Thread Tools  Search this Thread 
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Working directory via command line  Luiz  CFX  4  March 6, 2011 20:02 
premixed combustion  parallel solver error  Craig  FLUENT  0  October 14, 2008 15:07 
why the solver reject it? Anyone with experience?  bearcat  CFX  6  April 28, 2008 14:08 
Hydrogen Air combustion in a combustion chamber  popi  CFX  7  July 11, 2007 18:40 
Need Similar Combustion Tutorial as my model  Arnab  Siemens  0  May 24, 2005 18:23 