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

ReactingFoam

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

Reply
 
LinkBack Thread Tools Display Modes
Old   August 24, 2007, 07:59
Default Its because the chemistry mech
  #61
Super Moderator
 
niklas's Avatar
 
Niklas Nordin
Join Date: Mar 2009
Location: Stockholm, Sweden
Posts: 693
Rep Power: 21
niklas will become famous soon enoughniklas will become famous soon enough
Its because the chemistry mechanism is constructed at the standard temperature.

You can safely set the temperature range to 200K.
niklas is offline   Reply With Quote

Old   August 24, 2007, 09:32
Default Hallo Niklas, I have set the
  #62
Member
 
Marco Moscaritolo
Join Date: Mar 2009
Location: Bergamo, Italy
Posts: 33
Rep Power: 9
mavimo is on a distinguished road
Hallo Niklas,
I have set the temperature range to 200K and it work fine!

Thanks
Marco
mavimo is offline   Reply With Quote

Old   September 13, 2007, 18:39
Default Hi All, I'm running reactin
  #63
Member
 
victor
Join Date: Mar 2009
Location: mexico city, MX
Posts: 50
Rep Power: 9
torvic is on a distinguished road
Hi All,

I'm running reactingFoam in OF-1.4 64 bits and after a short time I get the following.

foam@ws-linux:~/OpenFOAM/foam-1.4/run/combustor/cpaperext/pruebas/prueba49b> #0 Foam::error::printStack(Foam:stream&)
#1 Foam::sigFpe::sigFpeHandler(int)
#2 ??
#3 Foam::sqrt(Foam::Field<double>&, Foam::UList<double> const&)
#4 Foam::tmp<foam::geometricfield<double,> > Foam::sqrt<foam::fvpatchfield,>(Foam::tmp<foam::ge ometricfield<double,> > const&)
#5 main
#6 __libc_start_main
#7 __gxx_personality_v0 at ../sysdeps/x86_64/elf/start.S:116

[1]+ Floating point exceptionreactingFoam $FOAM_RUN/combustor/cpaperext/pruebas prueba49b >reac1.log

But, if i run the same case in another machine with OF-1.4 32 bits the simulation goes on without displaying that message.

Just for clarity the 64 bit machine has SLED-10 SP1 and the 32 bit one OpenSUSE 10.2

Any comment is appreciated

thanks

V
torvic is offline   Reply With Quote

Old   September 13, 2007, 19:23
Default 1) possibly floating point exc
  #64
Super Moderator
 
Mattijs Janssens
Join Date: Mar 2009
Posts: 1,416
Rep Power: 17
mattijs is on a distinguished road
1) possibly floating point exception (sigfpe) behaviour is different on both machines. Do you have $FOAM_SIGFPE set on both machines.

2) There might be slightly different behaviour of floating point calculations, causing one architecture to have a number just within acceptable range whereas the other doesn't.

Add some printing of min and max of your field to your application and see how large it becomes.
mattijs is offline   Reply With Quote

Old   September 13, 2007, 19:40
Default Hi Mattijs Thanks so much f
  #65
Member
 
victor
Join Date: Mar 2009
Location: mexico city, MX
Posts: 50
Rep Power: 9
torvic is on a distinguished road
Hi Mattijs

Thanks so much for the quick response. As for $FOAM_SIGFPE, it's not set.

However, the 32 bit already blowed up, but it was in the 0.088 time step, whereas with the 64 bit it was in the 0.009. Moreover, the message is slighlty different, pointing to the chemistry model. I add it just for clarity:
Solving chemistry
#0 Foam::error::printStack(Foam:stream&)
#1 Foam::sigFpe::sigFpeHandler(int)
#2 Uninterpreted: [0xb7ef2420]
#3 Foam::chemistryModel::solve(double, double)
#4 main
#5 __libc_start_main
#6 __gxx_personality_v0 at /usr/src/packages/BUILD/glibc-2.3/csu/../sysdeps/i386/elf/start.S:122
Floating point exception


I will try your suggestions of printing min/max.

thanks so much again

V
torvic is offline   Reply With Quote

Old   September 19, 2008, 05:34
Default Dear All I am working with re
  #66
New Member
 
Vijayaratnam Piradeepan
Join Date: Mar 2009
Posts: 6
Rep Power: 9
piradeepan is on a distinguished road
Dear All
I am working with reacting Foam solver.
Please tell me, how I can set ignition for only 1 second.
Vijay
piradeepan is offline   Reply With Quote

Old   January 13, 2009, 08:54
Default I want to model H2 dispersion
  #67
Member
 
Hamed Aghajani
Join Date: Mar 2009
Location: London, UK
Posts: 77
Rep Power: 9
haghajani is on a distinguished road
I want to model H2 dispersion through a hole in a closed box, containing O2 in it.

To do this, I disabled chemistry and turbulentReacrion (off) in chemistryProperties of reactingFoam. The problem is, the solution is very slow due to reactingFoam timing.

1- What are the equations which reactingfoam solve?
2 Is reactingfoam the proper solver? to model such a problem!

Is there any tutorial for dispersion of one Gas in to 2nd Gas?

Many thanks,
Hamed.
haghajani is offline   Reply With Quote

Old   January 13, 2009, 11:01
Default I want to model H2 dispersion
  #68
Member
 
Hamed Aghajani
Join Date: Mar 2009
Location: London, UK
Posts: 77
Rep Power: 9
haghajani is on a distinguished road
I want to model H2 dispersion through a hole in a closed box, containing O2 in it.

To do this, I disabled chemistry and turbulentReacrion (off) in chemistryProperties of reactingFoam. The problem is, the solution is very slow due to reactingFoam timing.

1- What are the equations which reactingfoam solve?
2 Is reactingfoam the proper solver? to model such a problem!

Is there any tutorial for dispersion of one Gas in to 2nd Gas?

Many thanks,
Hamed.
haghajani is offline   Reply With Quote

Old   January 16, 2009, 09:49
Default Hi Hamed, source code for "
  #69
Member
 
Ville Tossavainen
Join Date: Mar 2009
Location: Helsinki, Finland
Posts: 60
Rep Power: 9
villet is on a distinguished road
Hi Hamed,

source code for "reactingFoam" can be found at "$FOAM_APP/solvers/combustion/reactingFoam". Some files are used from "XiFoam" solver (../XiFoam directory).

You can solve your problem using "reactingFoam" in the way you explained. But if your problem is purely mixing/dispersion problem of two species, you could start wit somewhat simpler solver. I was thinking adding a passive scalar (representing H2/O2) to compressible/incompressible turbulent/laminar solver. You can find examples on this forum and OF wiki (adding temperature equation).
villet is offline   Reply With Quote

Old   January 16, 2009, 11:07
Default Hi Ville, Thanks for your k
  #70
Member
 
Hamed Aghajani
Join Date: Mar 2009
Location: London, UK
Posts: 77
Rep Power: 9
haghajani is on a distinguished road
Hi Ville,

Thanks for your kind reply,

The problem is purely Mixing/dispersion.
To do this, first I run the PitzDaily case, H2 from upper inlet and O2 from lower inlet. The mixing behaviour is acceptable. But, by changing the Geometry and the problem a bit (evacuation of H2 stream in a closed container of O2), apart from the slowness of solution (due to time-step control nature of reactingFoam), and considering effect of Gravity in "enviromentalProperties", the solution goes crazy!!! It seems mass equation does not satisfies and some mass misses and the flow field is really far from thoughts.
Anyway, I'll try on the compressible/incompressible solvers.

Reg.
Hamed
haghajani is offline   Reply With Quote

Old   February 18, 2009, 01:47
Default hi how and were to define bo
  #71
New Member
 
pavan
Join Date: Mar 2009
Location: BANGALORE, KARNATAKA, INDIA
Posts: 7
Rep Power: 9
viji is on a distinguished road
hi
how and were to define boundary condition of separate inlets for air and fuel in a combustion problem.. i am working in XIfoam
viji is offline   Reply With Quote

Old   September 28, 2009, 07:44
Default reactingFoam - Ignition in Cell Zone - Ignition Temperature
  #72
New Member
 
Join Date: Mar 2009
Posts: 1
Rep Power: 0
Javier is on a distinguished road
Hello,

Im working with reactingFoam case, and I would like to know how to define the start of combustion (ignition) only in a group of cells (cell zone), not in the whole system.

Also I would like to know where the autoignition temperature must be introduced.

Thank you very much in advance!
Javier is offline   Reply With Quote

Old   September 28, 2009, 11:26
Default
  #73
Member
 
Christof Benz
Join Date: Mar 2009
Posts: 52
Rep Power: 9
chbenz is on a distinguished road
Hi,

please have a look at http://openfoamwiki.net/index.php/Contrib_reactingFoam
Another way is to set T high enough to ignite your composition. You could use setField. An example is here : $FOAM_APP/utilities/setFields

Christof
chbenz is offline   Reply With Quote

Old   December 16, 2009, 10:05
Default concerning the reactingFoam test case
  #74
New Member
 
Join Date: Sep 2009
Posts: 2
Rep Power: 0
aortwein is on a distinguished road
Quote:
Originally Posted by finch View Post
No, I did not change anything in the case, except the paths to CHEMKINFile and CHEMKINThermoFile. To make sure, I unpacked the archive again and ran the case. I got exactly the same error. Did Hannes, or anyone else, run the case all the way to completion?
It is possible to run the case. There is a minor error in the test case making trouble: in the "T"-file in the "0"-directory, the outlet is set to "type fixedValue" with a value of 800 K. By changing that to "type zeroGradient", the test case should run.

Greets, Andreas
aortwein is offline   Reply With Quote

Old   October 28, 2010, 13:53
Default Cmix
  #75
New Member
 
Silvano
Join Date: Aug 2010
Location: Chicago /Torino Us/Italy
Posts: 11
Rep Power: 8
SilPaut is on a distinguished road
Hi everybody!

i'm trying to simulate a dump combustor (turbulent/premixed case) with reactingFoam (OF-1.7).
The mixture is propane/air.

I have some troubles to equilibrate the the k-epsilon model. Actually it runs pretty well but I get a laminarization in correspondence of dQ "line". Which mean that mut and alphat drop down and there is not enough transport and diffusion of heat. Thus the trend of temperature along the radius (it is an axis-symmetric geometry) is very sharp comparing with the one from Fluent and experimental datas .

so I played a whit the Prt number and Cmix and i get good results with Prt=0.9 and Cmix=5.

My question is: does it make sense use a such as big value for Cmix?


so the coeffs i'm using are:

Code:
kEpsilonCoeffs
    {
        Cmu         0.09;
        C1          1.44;
        C2          1.92;
        C3          -0.33
        sigmak      1.0;
        sigmaEps    1.3;
        Prt         0.9;   
    }
and
Code:
Cmix =5
Thank you!!
Silvano
SilPaut is offline   Reply With Quote

Old   September 11, 2013, 15:22
Default regarding reactingFoam solver
  #76
Member
 
sonu
Join Date: Jul 2013
Location: delhi
Posts: 91
Rep Power: 5
yash.aesi is on a distinguished road
greetings oll ,

i am using OF 2.2 and trying to solve my case with the reacting Foam solver
the doubt in my mind is that do i need to use chemkin folder along with 0, constant and system in OF 2.2 as it already have a reaction and thermo. file in constant folder as given in counterflowflame2D.



thanks in advance
Regards
sonu .
yash.aesi 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
mvConvection in reactingFoam smehdi609 OpenFOAM Running, Solving & CFD 5 August 27, 2011 05:00
DieselFoam and ReactingFoam matteo_rosa_sentinella OpenFOAM Pre-Processing 4 September 28, 2009 10:35
ReactingFoam solver muthukaalai OpenFOAM Running, Solving & CFD 1 June 16, 2008 13:36
ReactingFoam without reactions lasb OpenFOAM Running, Solving & CFD 5 June 10, 2008 08:50
ReactingFoam error prashant24983 OpenFOAM Running, Solving & CFD 3 October 4, 2007 04:54


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