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

Cantera

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

Reply
 
LinkBack Thread Tools Display Modes
Old   August 7, 2013, 10:58
Default
  #101
TBO
Member
 
Join Date: May 2013
Location: Netherlands
Posts: 30
Rep Power: 4
TBO is on a distinguished road
Quote:
Originally Posted by gschaider View Post
Which SVN is not available? I tried the cantera-link and the two in section 2 ("Install ...") work too ... just not in the browser. Remove the svn:// and they will too

Thank you for the answer. Some people mention that is is also possible to modify the solver of later openfoam versions to run steady state chemistry. For me this is the preferred solution, since we are already running two different OpenFOAM versions (1.6 and 2.1.1), so we like to avoid using a third one. Can anyone tell me which steps are required to change reactingFoam into steadyreactingFoam (or even better, does anyone already have this solver).

I already had a quite extensive look on different fora, but I didn't find a clear answer yet ...
TBO is offline   Reply With Quote

Old   June 6, 2014, 06:53
Default reactingFoam to steadyState
  #102
Member
 
James
Join Date: Jul 2013
Posts: 37
Rep Power: 4
ni-openfoam-user is on a distinguished road
Hi,

Did you ever get an answer to your question: "Can anyone tell me which steps are required to change reactingFoam into steadyreactingFoam (or even better, does anyone already have this solver)?. "

I am also looking to run reactingFOAM in steady state. I am concerned at the moment with a high pressure release, I have switched off combustion and chemistry. What else do I need to modify to run in steady state.

Many thanks,

James
ni-openfoam-user is offline   Reply With Quote

Old   June 6, 2014, 10:50
Default
  #103
TBO
Member
 
Join Date: May 2013
Location: Netherlands
Posts: 30
Rep Power: 4
TBO is on a distinguished road
Dear James,

I never got a reply to this question, I did some tries with changing to pseudo-steady state solvers by changing the ddt schemes, however I was not succesfull with that (see also steady state combustion posts). Later I moved to the FGM based solver libOpenSmoke. This one works very well as steady state combustion solver.

Regards
TBO is offline   Reply With Quote

Old   March 18, 2015, 07:10
Default
  #104
New Member
 
Join Date: Feb 2015
Posts: 22
Rep Power: 2
Stefano Puggelli is on a distinguished road
Hi TBO,
I am working on the FGM combustion modeling in OpenFOAM and I find the only one already available by Kroger. I am testing it and I having some problems with the flamelet generator based on Cantera.
You say that you moved to the FGM model based on the libOpenSmoke and I also had this idea, but is there an FGM model already available with libOpenSmoke (because I see that it can be used only for non premixed combustion with only the equations for mixture fraction) or did you build on your own?
Stefano Puggelli is offline   Reply With Quote

Old   March 18, 2015, 10:11
Default
  #105
TBO
Member
 
Join Date: May 2013
Location: Netherlands
Posts: 30
Rep Power: 4
TBO is on a distinguished road
Hi Stefano,
I used this model for non premixed combustion, so the libOpenSmoke based solver was suitable for my simulations. Maybe repost the question on the libOpenSmoke forum to see if someone else has experience with using this solver for premixed combustion.

Regards,
TBO
TBO is offline   Reply With Quote

Old   March 18, 2015, 10:36
Default
  #106
New Member
 
Join Date: Feb 2015
Posts: 22
Rep Power: 2
Stefano Puggelli is on a distinguished road
Ok,
maybe I will ask on the libOpenSmoke forum.
Thank you,
Stefano
Stefano Puggelli is offline   Reply With Quote

Old   April 15, 2015, 11:51
Default
  #107
New Member
 
Ali Kadar
Join Date: Oct 2014
Location: Delft
Posts: 24
Rep Power: 2
flowAlways is on a distinguished road
Hello everyone,

has there been any new development on alternateSteadyReactingFoam which is now compatible only with old versions of OpenFOAM.
The transient solver reactingFoam is good but is not very appropriate considering the computational costs.
Is there any new development for a steady state combustion solver in OpenFOAM. ?

These are the related thread
modify reactingFoam in order to make it a steady state solver
Steady state chemistry
https://github.com/Unofficial-Extend...rnateChemistry
__________________
A good solution is one which does justice to the inner nature of the problem- Cornelius Lanczos in a letter to Albert Einstein on March 9, 1947
flowAlways is offline   Reply With Quote

Old   April 16, 2015, 02:29
Default
  #108
TBO
Member
 
Join Date: May 2013
Location: Netherlands
Posts: 30
Rep Power: 4
TBO is on a distinguished road
Quote:
Originally Posted by flowAlways View Post
Is there any new development for a steady state combustion solver in OpenFOAM. ?
Hi Ali,

Yes, there are some developments. By default OpenFoam has now the solver LTSReactingFoam, which is a steady state version of reactingFoam (available in OF 2.3).

Besides that, also e.g. the solvers based on libOpenSmoke might be worth considering

Regards,
TBO
TBO is offline   Reply With Quote

Old   April 17, 2015, 02:47
Default
  #109
New Member
 
Ali Kadar
Join Date: Oct 2014
Location: Delft
Posts: 24
Rep Power: 2
flowAlways is on a distinguished road
Thank you for pointing that. I did test runs using the LTSreactingFoam, but unfortunately it fails(after 27 time iterations) with the following error(attached).
It runs perfectly ok with reactingFoam and I get good results.
I could not understand the parameters in the PIMPLE Algorithm too.
What are alphaTemp, rDeltaTSmoothingCoeff, rDeltaTDampingCoeff.
and why there is no relaxation ?

I think I will have to look at the code to understand what they are doing and what is the idea of local time stepping(LTS). I know that if you want to go fast in time you have to solve more pressure-velocity couplings per iteration i.e. nOuterCorrectors > 1 with under-relaxation(idea behind SIMPLE). But they are not doing that!!. Whats the idea behind LTS ?

Any ideas would be really helpful !!

Code:
PIMPLE
{
    momentumPredictor no;
    nOuterCorrectors  1;
    nCorrectors     1;
    nNonOrthogonalCorrectors 0;

    maxDeltaT       1e-2;
    maxCo           1;
    alphaTemp       0.05;
    rDeltaTSmoothingCoeff 1;
    rDeltaTDampingCoeff 1;
}

relaxationFactors
{
    fields
    {
    }
    equations
    {
        ".*" 1;
    }
}
Code:
Time = 26

Time scales min/max:
    Flow        = 1.28808e-05, 0.01
    Temperature = 6.78274e-06, 1e+300
    Overall     = 6.78274e-06, 0.01
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG:  Solving for O2, Initial residual = 0.075075, Final residual = 0.00549963, No Iterations 4
DILUPBiCG:  Solving for H2O, Initial residual = 0.963064, Final residual = 0.075786, No Iterations 3
DILUPBiCG:  Solving for CH4, Initial residual = 0.978202, Final residual = 0.0737773, No Iterations 3
DILUPBiCG:  Solving for CO2, Initial residual = 0.962777, Final residual = 0.0756163, No Iterations 3
DILUPBiCG:  Solving for h, Initial residual = 0.996642, Final residual = 0.0833314, No Iterations 2
min/max(T) = 297.502, 2074.17
DICPCG:  Solving for p, Initial residual = 0.846321, Final residual = 9.13219e-07, No Iterations 226
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.14244e-05, global = -5.85091e-08, cumulative = 1.33362e-08
DILUPBiCG:  Solving for epsilon, Initial residual = 0.000435871, Final residual = 2.96746e-07, No Iterations 1
bounding epsilon, min: -45699 max: 4.85026e+08 average: 332153
DILUPBiCG:  Solving for k, Initial residual = 8.21421e-09, Final residual = 8.21421e-09, No Iterations 0
ExecutionTime = 32.79 s  ClockTime = 34 s

cellSource volumeTemperature output:
    average(sampledSurface) for T = 1285.89

faceSource faceTemperature output:
    average(outlet) for T = 1113.5

Time = 27

Time scales min/max:
    Flow        = 6.88187e-06, 0.01
    Temperature = -6.73481e+07, 1e+300
    Overall     = 6.88187e-06, 0.01
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
#0  Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam230/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1  Foam::sigFpe::sigHandler(int) in "/opt/openfoam230/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2   in "/lib/x86_64-linux-gnu/libc.so.6"
#3  Foam::EulerImplicit<Foam::chemistryModel<Foam::psiChemistryModel, Foam::sutherlandTransport<Foam::species::thermo<Foam::janafThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleEnthalpy> > > >::solve(Foam::Field<double>&, double&, double&, double&, double&) const in "/opt/openfoam230/platforms/linux64GccDPOpt/lib/libchemistryModel.so"
#4  double Foam::chemistryModel<Foam::psiChemistryModel, Foam::sutherlandTransport<Foam::species::thermo<Foam::janafThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleEnthalpy> > >::solve<Foam::Field<double> >(Foam::Field<double> const&) in "/opt/openfoam230/platforms/linux64GccDPOpt/lib/libchemistryModel.so"
#5  Foam::combustionModels::laminar<Foam::combustionModels::psiChemistryCombustion>::correct() in "/opt/openfoam230/platforms/linux64GccDPOpt/lib/libcombustionModels.so"
#6  Foam::combustionModels::PaSR<Foam::combustionModels::psiChemistryCombustion>::correct() in "/opt/openfoam230/platforms/linux64GccDPOpt/lib/libcombustionModels.so"
#7  
 at /home/opencfd/OpenFOAM/OpenFOAM-2.3.0/src/OpenFOAM/lnInclude/autoPtrI.H:174
#8  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#9  
 in "/opt/openfoam230/platforms/linux64GccDPOpt/bin/LTSReactingFoam"
Floating point exception (core dumped)
__________________
A good solution is one which does justice to the inner nature of the problem- Cornelius Lanczos in a letter to Albert Einstein on March 9, 1947
flowAlways 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



All times are GMT -4. The time now is 00:08.