Merge branch 'master' of https://github.com/OpenFOAM/OpenFOAM-dev

applications/utilities/postProcessing/noise/noise.C

@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2011-2016 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2011-2017 OpenFOAM Foundation
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -133,7 +133,7 @@ int main(int argc, char *argv[])
     #include "createFields.H"
 
     Info<< "Reading data file" << endl;
-    Function1Types::CSV<scalar> pData("pressure", dict, "Data");
+    Function1Types::CSV<scalar> pData("pressure", dict.subDict("pressureData"));
 
     // time history data
     const scalarField t(pData.x());

etc/caseDicts/postProcessing/lagrangian/dsmcFields

@@ -0,0 +1,27 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     |
+    \\  /    A nd           | Web:      www.OpenFOAM.org
+     \\/     M anipulation  |
+-------------------------------------------------------------------------------
+Description
+    Calculate intensive fields:
+    - UMean
+    - translationalT
+    - internalT
+    - overallT
+    from averaged extensive fields from a DSMC calculation.
+
+\*---------------------------------------------------------------------------*/
+
+type            dsmcFields;
+libs            ("liblagrangianFunctionObjects.so");
+
+fields          (rhoNMean rhoMMean momentumMean linearKEMean internalEMean
+                 iDofMean fDMean);
+
+executeControl  writeTime;
+writeControl    writeTime;
+
+// ************************************************************************* //

src/OpenFOAM/Make/files

@@ -76,8 +76,14 @@ primitives/quaternion/quaternion.C
 primitives/septernion/septernion.C
 primitives/triad/triad.C
 
-/* functions, data entries */
+/* Run-time selectable functions */
 primitives/functions/Function1/makeDataEntries.C
+primitives/functions/Function1/ramp/ramp.C
+primitives/functions/Function1/linear/linear.C
+primitives/functions/Function1/quadratic/quadratic.C
+primitives/functions/Function1/quarterSine/quarterSine.C
+primitives/functions/Function1/quarterCosine/quarterCosine.C
+primitives/functions/Function1/halfCosine/halfCosine.C
 primitives/functions/Polynomial/polynomialFunction.C
 
 primitives/subModelBase/subModelBase.C

src/OpenFOAM/db/dictionary/dictionary.H

@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2011-2016 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2011-2017 OpenFOAM Foundation
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -334,7 +334,19 @@ public:
             ) const;
 
             //- Find and return a T,
-            //  if not found return the given default value
+            //  if not found throw a fatal error.
+            //  If recursive, search parent dictionaries.
+            //  If patternMatch, use regular expressions.
+            template<class T>
+            T lookupType
+            (
+                const word&,
+                bool recursive=false,
+                bool patternMatch=true
+            ) const;
+
+            //- Find and return a T,
+            //  if not found return the given default value.
             //  If recursive, search parent dictionaries.
             //  If patternMatch, use regular expressions.
             template<class T>

src/OpenFOAM/db/dictionary/dictionaryTemplates.C

@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2011-2015 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2011-2017 OpenFOAM Foundation
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -28,6 +28,30 @@ License
 
 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
 
+template<class T>
+T Foam::dictionary::lookupType
+(
+    const word& keyword,
+    bool recursive,
+    bool patternMatch
+) const
+{
+    const entry* entryPtr = lookupEntryPtr(keyword, recursive, patternMatch);
+
+    if (entryPtr == nullptr)
+    {
+        FatalIOErrorInFunction
+        (
+            *this
+        )   << "keyword " << keyword << " is undefined in dictionary "
+            << name()
+            << exit(FatalIOError);
+    }
+
+    return pTraits<T>(entryPtr->stream());
+}
+
+
 template<class T>
 T Foam::dictionary::lookupOrDefault
 (

src/OpenFOAM/memory/autoPtr/autoPtr.H

@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2011-2016 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2011-2017 OpenFOAM Foundation
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -125,6 +125,9 @@ public:
             //- Return const object pointer
             inline const T* operator->() const;
 
+            //- Take over the object pointer from parameter
+            inline void operator=(T*);
+
             //- Take over the object pointer from parameter
             inline void operator=(const autoPtr<T>&);
 };

src/OpenFOAM/memory/autoPtr/autoPtrI.H

@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2011-2016 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2011-2017 OpenFOAM Foundation
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -190,6 +190,13 @@ inline const T* Foam::autoPtr<T>::operator->() const
 }
 
 
+template<class T>
+inline void Foam::autoPtr<T>::operator=(T* p)
+{
+    reset(p);
+}
+
+
 template<class T>
 inline void Foam::autoPtr<T>::operator=(const autoPtr<T>& ap)
 {

src/OpenFOAM/primitives/functions/Function1/CSV/CSV.C

@@ -203,18 +203,16 @@ template<class Type>
 Foam::Function1Types::CSV<Type>::CSV
 (
     const word& entryName,
-    const dictionary& dict,
-    const word& ext
+    const dictionary& dict
 )
 :
-    TableBase<Type>(entryName, dict.subDict(entryName + ext)),
-    coeffs_(dict.subDict(entryName + ext)),
-    nHeaderLine_(readLabel(coeffs_.lookup("nHeaderLine"))),
-    refColumn_(readLabel(coeffs_.lookup("refColumn"))),
-    componentColumns_(coeffs_.lookup("componentColumns")),
-    separator_(coeffs_.lookupOrDefault<string>("separator", string(","))[0]),
-    mergeSeparators_(readBool(coeffs_.lookup("mergeSeparators"))),
-    fName_(coeffs_.lookup("file"))
+    TableBase<Type>(entryName, dict),
+    nHeaderLine_(readLabel(dict.lookup("nHeaderLine"))),
+    refColumn_(readLabel(dict.lookup("refColumn"))),
+    componentColumns_(dict.lookup("componentColumns")),
+    separator_(dict.lookupOrDefault<string>("separator", string(","))[0]),
+    mergeSeparators_(readBool(dict.lookup("mergeSeparators"))),
+    fName_(dict.lookup("file"))
 {
     if (componentColumns_.size() != pTraits<Type>::nComponents)
     {

src/OpenFOAM/primitives/functions/Function1/CSV/CSV.H

@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2011-2016 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2011-2017 OpenFOAM Foundation
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -25,9 +25,11 @@ Class
     Foam::Function1Types::CSV
 
 Description
-    Templated CSV container data entry.  Reference column is always a scalar,
-    e.g. time
+    Templated CSV function.
 
+    Reference column is always a scalar, e.g. time.
+
+    Usage:
     \verbatim
         <entryName> csvFile;
         <entryName>Coeffs
@@ -75,9 +77,6 @@ class CSV
 {
     // Private data
 
-        //- Coefficients dictionary (for convenience on reading)
-        dictionary coeffs_;
-
         //- Number header lines
         label nHeaderLine_;
 
@@ -121,8 +120,7 @@ public:
         CSV
         (
             const word& entryName,
-            const dictionary& dict,
-            const word& ext = "Coeffs"
+            const dictionary& dict
         );
 
         //- Copy constructor

src/OpenFOAM/primitives/functions/Function1/Constant/Constant.H

@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2011-2016 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2011-2017 OpenFOAM Foundation
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -25,9 +25,9 @@ Class
     Foam::Function1Types::Constant
 
 Description
-    Templated basic entry that holds a constant value.
+    Templated function that returns a constant value.
 
-    Usage - for entry \<entryName\> having the value <value>:
+    Usage - for entry \<entryName\> returning the value <value>:
     \verbatim
         <entryName>    constant  <value>
     \endverbatim

src/OpenFOAM/primitives/functions/Function1/Function1/Function1.H

@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2011-2016 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2011-2017 OpenFOAM Foundation
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -195,6 +195,14 @@ public:
         add##SS##Type##ConstructorToTable_;
 
 
+#define makeScalarFunction1(SS)                                                \
+                                                                               \
+    defineTypeNameAndDebug(SS, 0);                                             \
+                                                                               \
+    Function1<scalar>::adddictionaryConstructorToTable<SS>                     \
+        add##SS##ConstructorToTable_;
+
+
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
 #ifdef NoRepository

src/OpenFOAM/primitives/functions/Function1/Function1/Function1New.C

@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2011-2016 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2011-2017 OpenFOAM Foundation
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -34,39 +34,70 @@ Foam::autoPtr<Foam::Function1<Type>> Foam::Function1<Type>::New
     const dictionary& dict
 )
 {
-    Istream& is(dict.lookup(entryName, false));
-
-    token firstToken(is);
-    word Function1Type;
-
-    if (!firstToken.isWord())
+    if (dict.isDict(entryName))
     {
-        is.putBack(firstToken);
-        return autoPtr<Function1<Type>>
-        (
-            new Function1Types::Constant<Type>(entryName, is)
-        );
+        const dictionary& coeffsDict(dict.subDict(entryName));
+
+        const word Function1Type(coeffsDict.lookup("type"));
+
+        typename dictionaryConstructorTable::iterator cstrIter =
+            dictionaryConstructorTablePtr_->find(Function1Type);
+
+        if (cstrIter == dictionaryConstructorTablePtr_->end())
+        {
+            FatalErrorInFunction
+                << "Unknown Function1 type "
+                << Function1Type << " for Function1 "
+                << entryName << nl << nl
+                << "Valid Function1 types are:" << nl
+                << dictionaryConstructorTablePtr_->sortedToc() << nl
+                << exit(FatalError);
+        }
+
+        return cstrIter()(entryName, coeffsDict);
     }
     else
     {
-        Function1Type = firstToken.wordToken();
+        Istream& is(dict.lookup(entryName, false));
+
+        token firstToken(is);
+        word Function1Type;
+
+        if (!firstToken.isWord())
+        {
+            is.putBack(firstToken);
+            return autoPtr<Function1<Type>>
+            (
+                new Function1Types::Constant<Type>(entryName, is)
+            );
+        }
+        else
+        {
+            Function1Type = firstToken.wordToken();
+        }
+
+        typename dictionaryConstructorTable::iterator cstrIter =
+            dictionaryConstructorTablePtr_->find(Function1Type);
+
+        if (cstrIter == dictionaryConstructorTablePtr_->end())
+        {
+            FatalErrorInFunction
+                << "Unknown Function1 type "
+                << Function1Type << " for Function1 "
+                << entryName << nl << nl
+                << "Valid Function1 types are:" << nl
+                << dictionaryConstructorTablePtr_->sortedToc() << nl
+                << exit(FatalError);
+        }
+
+        return cstrIter()
+        (
+            entryName,
+            dict.found(entryName + "Coeffs")
+          ? dict.subDict(entryName + "Coeffs")
+          : dict
+        );
     }
-
-    typename dictionaryConstructorTable::iterator cstrIter =
-        dictionaryConstructorTablePtr_->find(Function1Type);
-
-    if (cstrIter == dictionaryConstructorTablePtr_->end())
-    {
-        FatalErrorInFunction
-            << "Unknown Function1 type "
-            << Function1Type << " for Function1 "
-            << entryName << nl << nl
-            << "Valid Function1 types are:" << nl
-            << dictionaryConstructorTablePtr_->sortedToc() << nl
-            << exit(FatalError);
-    }
-
-    return cstrIter()(entryName, dict);
 }
 
 

src/OpenFOAM/primitives/functions/Function1/One/OneConstant.C

@@ -0,0 +1,84 @@
+/*---------------------------------------------------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     |
+    \\  /    A nd           | Copyright (C) 2017 OpenFOAM Foundation
+     \\/     M anipulation  |
+-------------------------------------------------------------------------------
+License
+    This file is part of OpenFOAM.
+
+    OpenFOAM is free software: you can redistribute it and/or modify it
+    under the terms of the GNU General Public License as published by
+    the Free Software Foundation, either version 3 of the License, or
+    (at your option) any later version.
+
+    OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
+    ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
+    FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
+    for more details.
+
+    You should have received a copy of the GNU General Public License
+    along with OpenFOAM.  If not, see <http://www.gnu.org/licenses/>.
+
+\*---------------------------------------------------------------------------*/
+
+#include "OneConstant.H"
+
+// * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
+
+template<class Type>
+Foam::Function1Types::OneConstant<Type>::OneConstant(const word& entryName)
+:
+    Function1<Type>(entryName)
+{}
+
+
+template<class Type>
+Foam::Function1Types::OneConstant<Type>::OneConstant
+(
+    const word& entryName,
+    const dictionary& dict
+)
+:
+    Function1<Type>(entryName)
+{}
+
+
+// * * * * * * * * * * * * * * * * Destructor  * * * * * * * * * * * * * * * //
+
+template<class Type>
+Foam::Function1Types::OneConstant<Type>::~OneConstant()
+{}
+
+
+// * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
+
+template<class Type>
+Type Foam::Function1Types::OneConstant<Type>::value(const scalar x) const
+{
+    return pTraits<Type>::one;
+}
+
+
+template<class Type>
+Type Foam::Function1Types::OneConstant<Type>::integrate
+(
+    const scalar x1,
+    const scalar x2
+) const
+{
+    return (x2 - x1)*pTraits<Type>::one;
+}
+
+
+template<class Type>
+void Foam::Function1Types::OneConstant<Type>::writeData(Ostream& os) const
+{
+    Function1<Type>::writeData(os);
+
+    os  << token::END_STATEMENT << nl;
+}
+
+
+// ************************************************************************* //

src/OpenFOAM/primitives/functions/Function1/One/OneConstant.H

@@ -0,0 +1,120 @@
+/*---------------------------------------------------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     |
+    \\  /    A nd           | Copyright (C) 2017 OpenFOAM Foundation
+     \\/     M anipulation  |
+-------------------------------------------------------------------------------
+License
+    This file is part of OpenFOAM.
+
+    OpenFOAM is free software: you can redistribute it and/or modify it
+    under the terms of the GNU General Public License as published by
+    the Free Software Foundation, either version 3 of the License, or
+    (at your option) any later version.
+
+    OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
+    ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
+    FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
+    for more details.
+
+    You should have received a copy of the GNU General Public License
+    along with OpenFOAM.  If not, see <http://www.gnu.org/licenses/>.
+
+Class
+    Foam::Function1Types::OneConstant
+
+Description
+    Templated function that returns the corresponding 1 (one).
+
+    Usage:
+    \verbatim
+        <entryName> one;
+    \endverbatim
+
+SourceFiles
+    OneConstant.C
+
+\*---------------------------------------------------------------------------*/
+
+#ifndef OneConstant_H
+#define OneConstant_H
+
+#include "Function1.H"
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+namespace Foam
+{
+namespace Function1Types
+{
+
+/*---------------------------------------------------------------------------*\
+                           Class OneConstant Declaration
+\*---------------------------------------------------------------------------*/
+
+template<class Type>
+class OneConstant
+:
+    public Function1<Type>
+{
+    // Private Member Functions
+
+        //- Disallow default bitwise assignment
+        void operator=(const OneConstant<Type>&);
+
+
+public:
+
+    // Runtime type information
+    TypeName("one");
+
+
+    // Constructors
+
+        //- Construct from entry name
+        OneConstant(const word& entryName);
+
+        //- Construct from entry name and dictionary
+        OneConstant(const word& entryName, const dictionary& dict);
+
+        //- Construct and return a clone
+        virtual tmp<Function1<Type>> clone() const
+        {
+            return tmp<Function1<Type>>(new OneConstant<Type>(*this));
+        }
+
+
+    //- Destructor
+    virtual ~OneConstant();
+
+
+    // Member Functions
+
+        //- Return constant value
+        Type value(const scalar) const;
+
+        //- Integrate between two values
+        Type integrate(const scalar x1, const scalar x2) const;
+
+        //- Write in dictionary format
+        virtual void writeData(Ostream& os) const;
+};
+
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+} // End namespace Function1Types
+} // End namespace Foam
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+#ifdef NoRepository
+    #include "OneConstant.C"
+#endif
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+#endif
+
+// ************************************************************************* //

src/OpenFOAM/primitives/functions/Function1/Sine/Sine.C

@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2016 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2016-2017 OpenFOAM Foundation
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -43,13 +43,12 @@ template<class Type>
 Foam::Function1Types::Sine<Type>::Sine
 (
     const word& entryName,
-    const dictionary& dict,
-    const word& ext
+    const dictionary& dict
 )
 :
     Function1<Type>(entryName)
 {
-    read(dict.subDict(entryName + ext));
+    read(dict);
 }
 
 

src/OpenFOAM/primitives/functions/Function1/Sine/Sine.H

@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2016 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2016-2017 OpenFOAM Foundation
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -132,8 +132,7 @@ public:
         Sine
         (
             const word& entryName,
-            const dictionary& dict,
-            const word& ext = "Coeffs"
+            const dictionary& dict
         );
 
         //- Copy constructor

src/OpenFOAM/primitives/functions/Function1/Square/Square.C

@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2016 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2016-2017 OpenFOAM Foundation
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -43,13 +43,12 @@ template<class Type>
 Foam::Function1Types::Square<Type>::Square
 (
     const word& entryName,
-    const dictionary& dict,
-    const word& ext
+    const dictionary& dict
 )
 :
     Function1<Type>(entryName)
 {
-    read(dict.subDict(entryName + ext));
+    read(dict);
 }
 
 

src/OpenFOAM/primitives/functions/Function1/Square/Square.H

@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2016 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2016-2017 OpenFOAM Foundation
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -139,8 +139,7 @@ public:
         Square
         (
             const word& entryName,
-            const dictionary& dict,
-            const word& ext = "Coeffs"
+            const dictionary& dict
         );
 
         //- Copy constructor

src/OpenFOAM/primitives/functions/Function1/Table/Table.H

@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2011-2016 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2011-2017 OpenFOAM Foundation
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -25,10 +25,12 @@ Class
     Foam::Function1Types::Table
 
 Description
-    Templated table container data entry. Items are stored in a list of
-    Tuple2's. First column is always stored as scalar entries. Data is read
-    in Tuple2 form, e.g. for an entry \<entryName\> that is (scalar, vector):
+    Templated table container function.
 
+    Items are stored in a list of Tuple2's. First column is always stored as
+    scalar entries. Data is read in Tuple2 form.
+
+    Usage:
     \verbatim
         <entryName>   table
         (

src/OpenFOAM/primitives/functions/Function1/TableFile/TableFile.C

@@ -34,11 +34,10 @@ Foam::Function1Types::TableFile<Type>::TableFile
     const dictionary& dict
 )
 :
-    TableBase<Type>(entryName, dict.subDict(entryName + "Coeffs")),
+    TableBase<Type>(entryName, dict),
     fName_("none")
 {
-    const dictionary coeffs(dict.subDict(entryName + "Coeffs"));
-    coeffs.lookup("file") >> fName_;
+    dict.lookup("file") >> fName_;
 
     fileName expandedFile(fName_);
     IFstream is(expandedFile.expand());

src/OpenFOAM/primitives/functions/Function1/TableFile/TableFile.H

@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2011-2016 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2011-2017 OpenFOAM Foundation
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -25,8 +25,9 @@ Class
     Foam::Function1Types::TableFile
 
 Description
-    Templated table container data entry where data is read from file.
+    Templated table container function where data is read from file.
 
+    Usage:
     \verbatim
         <entryName> tableFile;
         <entryName>Coeffs
@@ -37,7 +38,7 @@ Description
         }
     \endverbatim
 
-    Items are stored in a list of Tuple2's. First column is always stored as
+    Data is stored as a list of Tuple2's. First column is always stored as
     scalar entries.  Data is read in the form, e.g. for an entry \<entryName\>
     that is (scalar, vector):
     \verbatim
@@ -47,7 +48,6 @@ Description
         );
     \endverbatim
 
-
 SourceFiles
     TableFile.C
 

src/OpenFOAM/primitives/functions/Function1/Zero/ZeroConstant.C

@@ -0,0 +1,77 @@
+/*---------------------------------------------------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     |
+    \\  /    A nd           | Copyright (C) 2017 OpenFOAM Foundation
+     \\/     M anipulation  |
+-------------------------------------------------------------------------------
+License
+    This file is part of OpenFOAM.
+
+    OpenFOAM is free software: you can redistribute it and/or modify it
+    under the terms of the GNU General Public License as published by
+    the Free Software Foundation, either version 3 of the License, or
+    (at your option) any later version.
+
+    OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
+    ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
+    FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
+    for more details.
+
+    You should have received a copy of the GNU General Public License
+    along with OpenFOAM.  If not, see <http://www.gnu.org/licenses/>.
+
+\*---------------------------------------------------------------------------*/
+
+#include "ZeroConstant.H"
+
+// * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
+
+template<class Type>
+Foam::Function1Types::ZeroConstant<Type>::ZeroConstant
+(
+    const word& entryName,
+    const dictionary& dict
+)
+:
+    Function1<Type>(entryName)
+{}
+
+
+// * * * * * * * * * * * * * * * * Destructor  * * * * * * * * * * * * * * * //
+
+template<class Type>
+Foam::Function1Types::ZeroConstant<Type>::~ZeroConstant()
+{}
+
+
+// * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
+
+template<class Type>
+Type Foam::Function1Types::ZeroConstant<Type>::value(const scalar x) const
+{
+    return pTraits<Type>::zero;
+}
+
+
+template<class Type>
+Type Foam::Function1Types::ZeroConstant<Type>::integrate
+(
+    const scalar x1,
+    const scalar x2
+) const
+{
+    return pTraits<Type>::zero;
+}
+
+
+template<class Type>
+void Foam::Function1Types::ZeroConstant<Type>::writeData(Ostream& os) const
+{
+    Function1<Type>::writeData(os);
+
+    os  << token::END_STATEMENT << nl;
+}
+
+
+// ************************************************************************* //

src/OpenFOAM/primitives/functions/Function1/Zero/ZeroConstant.H

@@ -0,0 +1,117 @@
+/*---------------------------------------------------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     |
+    \\  /    A nd           | Copyright (C) 2017 OpenFOAM Foundation
+     \\/     M anipulation  |
+-------------------------------------------------------------------------------
+License
+    This file is part of OpenFOAM.
+
+    OpenFOAM is free software: you can redistribute it and/or modify it
+    under the terms of the GNU General Public License as published by
+    the Free Software Foundation, either version 3 of the License, or
+    (at your option) any later version.
+
+    OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
+    ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
+    FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
+    for more details.
+
+    You should have received a copy of the GNU General Public License
+    along with OpenFOAM.  If not, see <http://www.gnu.org/licenses/>.
+
+Class
+    Foam::Function1Types::ZeroConstant
+
+Description
+    Templated function that returns the corresponding 0 (zero).
+
+    Usage:
+    \verbatim
+        <entryName> zero;
+    \endverbatim
+
+SourceFiles
+    ZeroConstant.C
+
+\*---------------------------------------------------------------------------*/
+
+#ifndef ZeroConstant_H
+#define ZeroConstant_H
+
+#include "Function1.H"
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+namespace Foam
+{
+namespace Function1Types
+{
+
+/*---------------------------------------------------------------------------*\
+                           Class ZeroConstant Declaration
+\*---------------------------------------------------------------------------*/
+
+template<class Type>
+class ZeroConstant
+:
+    public Function1<Type>
+{
+    // Private Member Functions
+
+        //- Disallow default bitwise assignment
+        void operator=(const ZeroConstant<Type>&);
+
+
+public:
+
+    // Runtime type information
+    TypeName("zero");
+
+
+    // Constructors
+
+        //- Construct from entry name and dictionary
+        ZeroConstant(const word& entryName, const dictionary& dict);
+
+        //- Construct and return a clzero
+        virtual tmp<Function1<Type>> clzero() const
+        {
+            return tmp<Function1<Type>>(new ZeroConstant<Type>(*this));
+        }
+
+
+    //- Destructor
+    virtual ~ZeroConstant();
+
+
+    // Member Functions
+
+        //- Return constant value
+        Type value(const scalar) const;
+
+        //- Integrate between two values
+        Type integrate(const scalar x1, const scalar x2) const;
+
+        //- Write in dictionary format
+        virtual void writeData(Ostream& os) const;
+};
+
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+} // End namespace Function1Types
+} // End namespace Foam
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+#ifdef NoRepository
+    #include "ZeroConstant.C"
+#endif
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+#endif
+
+// ************************************************************************* //

src/OpenFOAM/primitives/functions/Function1/halfCosine/halfCosine.C

@@ -0,0 +1,66 @@
+/*---------------------------------------------------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     |
+    \\  /    A nd           | Copyright (C) 2017 OpenFOAM Foundation
+     \\/     M anipulation  |
+-------------------------------------------------------------------------------
+License
+    This file is part of OpenFOAM.
+
+    OpenFOAM is free software: you can redistribute it and/or modify it
+    under the terms of the GNU General Public License as published by
+    the Free Software Foundation, either version 3 of the License, or
+    (at your option) any later version.
+
+    OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
+    ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
+    FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
+    for more details.
+
+    You should have received a copy of the GNU General Public License
+    along with OpenFOAM.  If not, see <http://www.gnu.org/licenses/>.
+
+\*---------------------------------------------------------------------------*/
+
+#include "halfCosine.H"
+#include "mathematicalConstants.H"
+
+// * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
+
+namespace Foam
+{
+namespace Function1Types
+{
+    makeScalarFunction1(halfCosine);
+}
+}
+
+
+// * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
+
+Foam::Function1Types::halfCosine::halfCosine
+(
+    const word& entryName,
+    const dictionary& dict
+)
+:
+    ramp(entryName, dict)
+{}
+
+
+// * * * * * * * * * * * * * * * * Destructor  * * * * * * * * * * * * * * * //
+
+Foam::Function1Types::halfCosine::~halfCosine()
+{}
+
+
+// * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
+
+Foam::scalar Foam::Function1Types::halfCosine::value(const scalar t) const
+{
+    return 0.5*(1 - cos(constant::mathematical::pi*linearRamp(t)));
+}
+
+
+// ************************************************************************* //

src/OpenFOAM/primitives/functions/Function1/halfCosine/halfCosine.H

@@ -0,0 +1,101 @@
+/*---------------------------------------------------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     |
+    \\  /    A nd           | Copyright (C) 2017 OpenFOAM Foundation
+     \\/     M anipulation  |
+-------------------------------------------------------------------------------
+License
+    This file is part of OpenFOAM.
+
+    OpenFOAM is free software: you can redistribute it and/or modify it
+    under the terms of the GNU General Public License as published by
+    the Free Software Foundation, either version 3 of the License, or
+    (at your option) any later version.
+
+    OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
+    ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
+    FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
+    for more details.
+
+    You should have received a copy of the GNU General Public License
+    along with OpenFOAM.  If not, see <http://www.gnu.org/licenses/>.
+
+Class
+    Foam::Function1Types::halfCosine
+
+Description
+    Half-cosine ramp function starting from 0 and increasing to 1 from \c start
+    over the \c duration and remaining at 1 thereafter.
+
+See also
+    Foam::Function1Types::ramp
+
+SourceFiles
+    halfCosine.C
+
+\*---------------------------------------------------------------------------*/
+
+#ifndef halfCosine_H
+#define halfCosine_H
+
+#include "ramp.H"
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+namespace Foam
+{
+namespace Function1Types
+{
+
+/*---------------------------------------------------------------------------*\
+                           Class halfCosine Declaration
+\*---------------------------------------------------------------------------*/
+
+class halfCosine
+:
+    public ramp
+{
+    // Private Member Functions
+
+        //- Disallow default bitwise assignment
+        void operator=(const halfCosine&);
+
+
+public:
+
+    // Runtime type information
+    TypeName("halfCosine");
+
+
+    // Constructors
+
+        //- Construct from entry name and dictionary
+        halfCosine
+        (
+            const word& entryName,
+            const dictionary& dict
+        );
+
+
+    //- Destructor
+    virtual ~halfCosine();
+
+
+    // Member Functions
+
+        //- Return value for time t
+        scalar value(const scalar t) const;
+};
+
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+} // End namespace Function1Types
+} // End namespace Foam
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+#endif
+
+// ************************************************************************* //

src/OpenFOAM/primitives/functions/Function1/linear/linear.C

@@ -0,0 +1,65 @@
+/*---------------------------------------------------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     |
+    \\  /    A nd           | Copyright (C) 2017 OpenFOAM Foundation
+     \\/     M anipulation  |
+-------------------------------------------------------------------------------
+License
+    This file is part of OpenFOAM.
+
+    OpenFOAM is free software: you can redistribute it and/or modify it
+    under the terms of the GNU General Public License as published by
+    the Free Software Foundation, either version 3 of the License, or
+    (at your option) any later version.
+
+    OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
+    ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
+    FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
+    for more details.
+
+    You should have received a copy of the GNU General Public License
+    along with OpenFOAM.  If not, see <http://www.gnu.org/licenses/>.
+
+\*---------------------------------------------------------------------------*/
+
+#include "linear.H"
+
+// * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
+
+namespace Foam
+{
+namespace Function1Types
+{
+    makeScalarFunction1(linear);
+}
+}
+
+
+// * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
+
+Foam::Function1Types::linear::linear
+(
+    const word& entryName,
+    const dictionary& dict
+)
+:
+    ramp(entryName, dict)
+{}
+
+
+// * * * * * * * * * * * * * * * * Destructor  * * * * * * * * * * * * * * * //
+
+Foam::Function1Types::linear::~linear()
+{}
+
+
+// * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
+
+Foam::scalar Foam::Function1Types::linear::value(const scalar t) const
+{
+    return linearRamp(t);
+}
+
+
+// ************************************************************************* //

src/OpenFOAM/primitives/functions/Function1/linear/linear.H

@@ -0,0 +1,101 @@
+/*---------------------------------------------------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     |
+    \\  /    A nd           | Copyright (C) 2017 OpenFOAM Foundation
+     \\/     M anipulation  |
+-------------------------------------------------------------------------------
+License
+    This file is part of OpenFOAM.
+
+    OpenFOAM is free software: you can redistribute it and/or modify it
+    under the terms of the GNU General Public License as published by
+    the Free Software Foundation, either version 3 of the License, or
+    (at your option) any later version.
+
+    OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
+    ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
+    FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
+    for more details.
+
+    You should have received a copy of the GNU General Public License
+    along with OpenFOAM.  If not, see <http://www.gnu.org/licenses/>.
+
+Class
+    Foam::Function1Types::linear
+
+Description
+    Linear ramp function starting from 0 and increasing linearly to 1 from \c
+    start over the \c duration and remaining at 1 thereafter.
+
+See also
+    Foam::Function1Types::ramp
+
+SourceFiles
+    linear.C
+
+\*---------------------------------------------------------------------------*/
+
+#ifndef linear_H
+#define linear_H
+
+#include "ramp.H"
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+namespace Foam
+{
+namespace Function1Types
+{
+
+/*---------------------------------------------------------------------------*\
+                           Class linear Declaration
+\*---------------------------------------------------------------------------*/
+
+class linear
+:
+    public ramp
+{
+    // Private Member Functions
+
+        //- Disallow default bitwise assignment
+        void operator=(const linear&);
+
+
+public:
+
+    // Runtime type information
+    TypeName("linear");
+
+
+    // Constructors
+
+        //- Construct from entry name and dictionary
+        linear
+        (
+            const word& entryName,
+            const dictionary& dict
+        );
+
+
+    //- Destructor
+    virtual ~linear();
+
+
+    // Member Functions
+
+        //- Return value for time t
+        scalar value(const scalar t) const;
+};
+
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+} // End namespace Function1Types
+} // End namespace Foam
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+#endif
+
+// ************************************************************************* //

src/OpenFOAM/primitives/functions/Function1/makeDataEntries.C

@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2011-2016 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2011-2017 OpenFOAM Foundation
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -24,6 +24,8 @@ License
 \*---------------------------------------------------------------------------*/
 
 #include "Constant.H"
+#include "ZeroConstant.H"
+#include "OneConstant.H"
 #include "PolynomialEntry.H"
 #include "Sine.H"
 #include "Square.H"
@@ -38,6 +40,8 @@ License
 #define makeFunction1s(Type)                                                   \
     makeFunction1(Type);                                                       \
     makeFunction1Type(Constant, Type);                                         \
+    makeFunction1Type(ZeroConstant, Type);                                     \
+    makeFunction1Type(OneConstant, Type);                                      \
     makeFunction1Type(Polynomial, Type);                                       \
     makeFunction1Type(Sine, Type);                                             \
     makeFunction1Type(Square, Type);                                           \

src/OpenFOAM/primitives/functions/Function1/quadratic/quadratic.C

@@ -0,0 +1,65 @@
+/*---------------------------------------------------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     |
+    \\  /    A nd           | Copyright (C) 2017 OpenFOAM Foundation
+     \\/     M anipulation  |
+-------------------------------------------------------------------------------
+License
+    This file is part of OpenFOAM.
+
+    OpenFOAM is free software: you can redistribute it and/or modify it
+    under the terms of the GNU General Public License as published by
+    the Free Software Foundation, either version 3 of the License, or
+    (at your option) any later version.
+
+    OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
+    ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
+    FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
+    for more details.
+
+    You should have received a copy of the GNU General Public License
+    along with OpenFOAM.  If not, see <http://www.gnu.org/licenses/>.
+
+\*---------------------------------------------------------------------------*/
+
+#include "quadratic.H"
+
+// * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
+
+namespace Foam
+{
+namespace Function1Types
+{
+    makeScalarFunction1(quadratic);
+}
+}
+
+
+// * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
+
+Foam::Function1Types::quadratic::quadratic
+(
+    const word& entryName,
+    const dictionary& dict
+)
+:
+    ramp(entryName, dict)
+{}
+
+
+// * * * * * * * * * * * * * * * * Destructor  * * * * * * * * * * * * * * * //
+
+Foam::Function1Types::quadratic::~quadratic()
+{}
+
+
+// * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
+
+Foam::scalar Foam::Function1Types::quadratic::value(const scalar t) const
+{
+    return sqr(linearRamp(t));
+}
+
+
+// ************************************************************************* //

src/OpenFOAM/primitives/functions/Function1/quadratic/quadratic.H

@@ -0,0 +1,101 @@
+/*---------------------------------------------------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     |
+    \\  /    A nd           | Copyright (C) 2017 OpenFOAM Foundation
+     \\/     M anipulation  |
+-------------------------------------------------------------------------------
+License
+    This file is part of OpenFOAM.
+
+    OpenFOAM is free software: you can redistribute it and/or modify it
+    under the terms of the GNU General Public License as published by
+    the Free Software Foundation, either version 3 of the License, or
+    (at your option) any later version.
+
+    OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
+    ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
+    FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
+    for more details.
+
+    You should have received a copy of the GNU General Public License
+    along with OpenFOAM.  If not, see <http://www.gnu.org/licenses/>.
+
+Class
+    Foam::Function1Types::quadratic
+
+Description
+    Quadratic ramp function starting from 0 and increasing quadratically to 1
+    from \c t_0 over the \c duration and remaining at 1 thereafter.
+
+See also
+    Foam::Function1Types::ramp
+
+SourceFiles
+    quadratic.C
+
+\*---------------------------------------------------------------------------*/
+
+#ifndef quadratic_H
+#define quadratic_H
+
+#include "ramp.H"
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+namespace Foam
+{
+namespace Function1Types
+{
+
+/*---------------------------------------------------------------------------*\
+                           Class quadratic Declaration
+\*---------------------------------------------------------------------------*/
+
+class quadratic
+:
+    public ramp
+{
+    // Private Member Functions
+
+        //- Disallow default bitwise assignment
+        void operator=(const quadratic&);
+
+
+public:
+
+    // Runtime type information
+    TypeName("quadratic");
+
+
+    // Constructors
+
+        //- Construct from entry name and dictionary
+        quadratic
+        (
+            const word& entryName,
+            const dictionary& dict
+        );
+
+
+    //- Destructor
+    virtual ~quadratic();
+
+
+    // Member Functions
+
+        //- Return value for time t
+        scalar value(const scalar t) const;
+};
+
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+} // End namespace Function1Types
+} // End namespace Foam
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+#endif
+
+// ************************************************************************* //

src/OpenFOAM/primitives/functions/Function1/quarterCosine/quarterCosine.C

@@ -0,0 +1,66 @@
+/*---------------------------------------------------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     |
+    \\  /    A nd           | Copyright (C) 2017 OpenFOAM Foundation
+     \\/     M anipulation  |
+-------------------------------------------------------------------------------
+License
+    This file is part of OpenFOAM.
+
+    OpenFOAM is free software: you can redistribute it and/or modify it
+    under the terms of the GNU General Public License as published by
+    the Free Software Foundation, either version 3 of the License, or
+    (at your option) any later version.
+
+    OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
+    ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
+    FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
+    for more details.
+
+    You should have received a copy of the GNU General Public License
+    along with OpenFOAM.  If not, see <http://www.gnu.org/licenses/>.
+
+\*---------------------------------------------------------------------------*/
+
+#include "quarterCosine.H"
+#include "mathematicalConstants.H"
+
+// * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
+
+namespace Foam
+{
+namespace Function1Types
+{
+    makeScalarFunction1(quarterCosine);
+}
+}
+
+
+// * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
+
+Foam::Function1Types::quarterCosine::quarterCosine
+(
+    const word& entryName,
+    const dictionary& dict
+)
+:
+    ramp(entryName, dict)
+{}
+
+
+// * * * * * * * * * * * * * * * * Destructor  * * * * * * * * * * * * * * * //
+
+Foam::Function1Types::quarterCosine::~quarterCosine()
+{}
+
+
+// * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
+
+Foam::scalar Foam::Function1Types::quarterCosine::value(const scalar t) const
+{
+    return 1 - cos(0.5*constant::mathematical::pi*linearRamp(t));
+}
+
+
+// ************************************************************************* //

src/OpenFOAM/primitives/functions/Function1/quarterCosine/quarterCosine.H

@@ -0,0 +1,101 @@
+/*---------------------------------------------------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     |
+    \\  /    A nd           | Copyright (C) 2017 OpenFOAM Foundation
+     \\/     M anipulation  |
+-------------------------------------------------------------------------------
+License
+    This file is part of OpenFOAM.
+
+    OpenFOAM is free software: you can redistribute it and/or modify it
+    under the terms of the GNU General Public License as published by
+    the Free Software Foundation, either version 3 of the License, or
+    (at your option) any later version.
+
+    OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
+    ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
+    FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
+    for more details.
+
+    You should have received a copy of the GNU General Public License
+    along with OpenFOAM.  If not, see <http://www.gnu.org/licenses/>.
+
+Class
+    Foam::Function1Types::quarterCosine
+
+Description
+    Quarter-cosine ramp function starting from 0 and increasing to 1 from \c
+    start over the \c duration and remaining at 1 thereafter.
+
+See also
+    Foam::Function1Types::ramp
+
+SourceFiles
+    quarterCosine.C
+
+\*---------------------------------------------------------------------------*/
+
+#ifndef quarterCosine_H
+#define quarterCosine_H
+
+#include "ramp.H"
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+namespace Foam
+{
+namespace Function1Types
+{
+
+/*---------------------------------------------------------------------------*\
+                           Class quarterCosine Declaration
+\*---------------------------------------------------------------------------*/
+
+class quarterCosine
+:
+    public ramp
+{
+    // Private Member Functions
+
+        //- Disallow default bitwise assignment
+        void operator=(const quarterCosine&);
+
+
+public:
+
+    // Runtime type information
+    TypeName("quarterCosine");
+
+
+    // Constructors
+
+        //- Construct from entry name and dictionary
+        quarterCosine
+        (
+            const word& entryName,
+            const dictionary& dict
+        );
+
+
+    //- Destructor
+    virtual ~quarterCosine();
+
+
+    // Member Functions
+
+        //- Return value for time t
+        scalar value(const scalar t) const;
+};
+
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+} // End namespace Function1Types
+} // End namespace Foam
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+#endif
+
+// ************************************************************************* //

src/OpenFOAM/primitives/functions/Function1/quarterSine/quarterSine.C

@@ -0,0 +1,66 @@
+/*---------------------------------------------------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     |
+    \\  /    A nd           | Copyright (C) 2017 OpenFOAM Foundation
+     \\/     M anipulation  |
+-------------------------------------------------------------------------------
+License
+    This file is part of OpenFOAM.
+
+    OpenFOAM is free software: you can redistribute it and/or modify it
+    under the terms of the GNU General Public License as published by
+    the Free Software Foundation, either version 3 of the License, or
+    (at your option) any later version.
+
+    OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
+    ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
+    FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
+    for more details.
+
+    You should have received a copy of the GNU General Public License
+    along with OpenFOAM.  If not, see <http://www.gnu.org/licenses/>.
+
+\*---------------------------------------------------------------------------*/
+
+#include "quarterSine.H"
+#include "mathematicalConstants.H"
+
+// * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
+
+namespace Foam
+{
+namespace Function1Types
+{
+    makeScalarFunction1(quarterSine);
+}
+}
+
+
+// * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
+
+Foam::Function1Types::quarterSine::quarterSine
+(
+    const word& entryName,
+    const dictionary& dict
+)
+:
+    ramp(entryName, dict)
+{}
+
+
+// * * * * * * * * * * * * * * * * Destructor  * * * * * * * * * * * * * * * //
+
+Foam::Function1Types::quarterSine::~quarterSine()
+{}
+
+
+// * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
+
+Foam::scalar Foam::Function1Types::quarterSine::value(const scalar t) const
+{
+    return sin(0.5*constant::mathematical::pi*linearRamp(t));
+}
+
+
+// ************************************************************************* //

src/OpenFOAM/primitives/functions/Function1/quarterSine/quarterSine.H

@@ -0,0 +1,101 @@
+/*---------------------------------------------------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     |
+    \\  /    A nd           | Copyright (C) 2017 OpenFOAM Foundation
+     \\/     M anipulation  |
+-------------------------------------------------------------------------------
+License
+    This file is part of OpenFOAM.
+
+    OpenFOAM is free software: you can redistribute it and/or modify it
+    under the terms of the GNU General Public License as published by
+    the Free Software Foundation, either version 3 of the License, or
+    (at your option) any later version.
+
+    OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
+    ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
+    FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
+    for more details.
+
+    You should have received a copy of the GNU General Public License
+    along with OpenFOAM.  If not, see <http://www.gnu.org/licenses/>.
+
+Class
+    Foam::Function1Types::quarterSine
+
+Description
+    Quarter-sine ramp function starting from 0 and increasing to 1 from \c start
+    over the \c duration and remaining at 1 thereafter.
+
+See also
+    Foam::Function1Types::ramp
+
+SourceFiles
+    quarterSine.C
+
+\*---------------------------------------------------------------------------*/
+
+#ifndef quarterSine_H
+#define quarterSine_H
+
+#include "ramp.H"
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+namespace Foam
+{
+namespace Function1Types
+{
+
+/*---------------------------------------------------------------------------*\
+                           Class quarterSine Declaration
+\*---------------------------------------------------------------------------*/
+
+class quarterSine
+:
+    public ramp
+{
+    // Private Member Functions
+
+        //- Disallow default bitwise assignment
+        void operator=(const quarterSine&);
+
+
+public:
+
+    // Runtime type information
+    TypeName("quarterSine");
+
+
+    // Constructors
+
+        //- Construct from entry name and dictionary
+        quarterSine
+        (
+            const word& entryName,
+            const dictionary& dict
+        );
+
+
+    //- Destructor
+    virtual ~quarterSine();
+
+
+    // Member Functions
+
+        //- Return value for time t
+        scalar value(const scalar t) const;
+};
+
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+} // End namespace Function1Types
+} // End namespace Foam
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+#endif
+
+// ************************************************************************* //

src/OpenFOAM/primitives/functions/Function1/ramp/ramp.C

@@ -0,0 +1,69 @@
+/*---------------------------------------------------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     |
+    \\  /    A nd           | Copyright (C) 2017 OpenFOAM Foundation
+     \\/     M anipulation  |
+-------------------------------------------------------------------------------
+License
+    This file is part of OpenFOAM.
+
+    OpenFOAM is free software: you can redistribute it and/or modify it
+    under the terms of the GNU General Public License as published by
+    the Free Software Foundation, either version 3 of the License, or
+    (at your option) any later version.
+
+    OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
+    ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
+    FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
+    for more details.
+
+    You should have received a copy of the GNU General Public License
+    along with OpenFOAM.  If not, see <http://www.gnu.org/licenses/>.
+
+\*---------------------------------------------------------------------------*/
+
+#include "ramp.H"
+
+// * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
+
+void Foam::Function1Types::ramp::read(const dictionary& coeffs)
+{
+    start_ = coeffs.lookupOrDefault<scalar>("start", 0);
+    duration_ = coeffs.lookupType<scalar>("duration");
+}
+
+
+Foam::Function1Types::ramp::ramp
+(
+    const word& entryName,
+    const dictionary& dict
+)
+:
+    Function1<scalar>(entryName)
+{
+    read(dict);
+}
+
+
+// * * * * * * * * * * * * * * * * Destructor  * * * * * * * * * * * * * * * //
+
+Foam::Function1Types::ramp::~ramp()
+{}
+
+
+// * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
+
+void Foam::Function1Types::ramp::writeData(Ostream& os) const
+{
+    Function1<scalar>::writeData(os);
+    os  << token::END_STATEMENT << nl;
+    os  << indent << word(this->name() + "Coeffs") << nl;
+    os  << indent << token::BEGIN_BLOCK << incrIndent << nl;
+    os.writeKeyword("start") << start_ << token::END_STATEMENT << nl;
+    os.writeKeyword("duration") << duration_ << token::END_STATEMENT << nl;
+    os  << decrIndent << indent << token::END_BLOCK << endl;
+}
+
+
+// ************************************************************************* //

src/OpenFOAM/primitives/functions/Function1/ramp/ramp.H

@@ -0,0 +1,150 @@
+/*---------------------------------------------------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     |
+    \\  /    A nd           | Copyright (C) 2017 OpenFOAM Foundation
+     \\/     M anipulation  |
+-------------------------------------------------------------------------------
+License
+    This file is part of OpenFOAM.
+
+    OpenFOAM is free software: you can redistribute it and/or modify it
+    under the terms of the GNU General Public License as published by
+    the Free Software Foundation, either version 3 of the License, or
+    (at your option) any later version.
+
+    OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
+    ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
+    FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
+    for more details.
+
+    You should have received a copy of the GNU General Public License
+    along with OpenFOAM.  If not, see <http://www.gnu.org/licenses/>.
+
+Class
+    Foam::Function1Types::ramp
+
+Description
+    Ramp function base class for the set of scalar functions starting from 0 and
+    increasing monotonically to 1 from \c start over the \c duration and
+    remaining at 1 thereafter.
+
+    Usage:
+    \verbatim
+        <entryName> <rampFunction>;
+        <entryName>Coeffs
+        {
+            start     10;
+            duration  20;
+        }
+    \endverbatim
+    or
+    \verbatim
+        <entryName>
+        {
+            type      <rampFunction>;
+            start     10;
+            duration  20;
+        }
+    \endverbatim
+
+    Where:
+    \table
+        Property | Description  | Required | Default value
+        start    | Start time   | no       | 0
+        duration | Duration     | yes      |
+    \endtable
+
+See also
+    Foam::Function1
+
+SourceFiles
+    ramp.C
+
+\*---------------------------------------------------------------------------*/
+
+#ifndef ramp_H
+#define ramp_H
+
+#include "Function1.H"
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+namespace Foam
+{
+namespace Function1Types
+{
+
+/*---------------------------------------------------------------------------*\
+                           Class ramp Declaration
+\*---------------------------------------------------------------------------*/
+
+class ramp
+:
+    public Function1<scalar>
+{
+protected:
+
+    // Protected data
+
+        //- Start-time of the ramp function
+        scalar start_;
+
+        //- Duration of the ramp function
+        scalar duration_;
+
+        //- Simple linear ramp function
+        //  which form the basis of many more complex ramp functions
+        inline scalar linearRamp(const scalar t) const
+        {
+            return max(min((t - start_)/duration_, 1), 0);
+        }
+
+
+private:
+
+    // Private Member Functions
+
+        //- Read the coefficients from the given dictionary
+        void read(const dictionary& coeffs);
+
+        //- Disallow default bitwise assignment
+        void operator=(const ramp&);
+
+
+public:
+
+    // Constructors
+
+        //- Construct from entry name and dictionary
+        ramp
+        (
+            const word& entryName,
+            const dictionary& dict
+        );
+
+
+    //- Destructor
+    virtual ~ramp();
+
+
+    // Member Functions
+
+        //- Return value for time t
+        scalar value(const scalar t) const = 0;
+
+        //- Write in dictionary format
+        virtual void writeData(Ostream& os) const;
+};
+
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+} // End namespace Function1Types
+} // End namespace Foam
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+#endif
+
+// ************************************************************************* //

src/combustionModels/EDC/EDC.C

@@ -0,0 +1,248 @@
+/*---------------------------------------------------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     |
+    \\  /    A nd           | Copyright (C) 2017 OpenFOAM Foundation
+     \\/     M anipulation  |
+-------------------------------------------------------------------------------
+License
+    This file is part of OpenFOAM.
+
+    OpenFOAM is free software: you can redistribute it and/or modify it
+    under the terms of the GNU General Public License as published by
+    the Free Software Foundation, either version 3 of the License, or
+    (at your option) any later version.
+
+    OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
+    ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
+    FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
+    for more details.
+
+    You should have received a copy of the GNU General Public License
+    along with OpenFOAM.  If not, see <http://www.gnu.org/licenses/>.
+
+\*---------------------------------------------------------------------------*/
+
+#include "EDC.H"
+
+// * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
+
+template<class Type>
+Foam::combustionModels::EDC<Type>::EDC
+(
+    const word& modelType,
+    const fvMesh& mesh,
+    const word& combustionProperties,
+    const word& phaseName
+)
+:
+    laminar<Type>(modelType, mesh, combustionProperties, phaseName),
+    version_
+    (
+        EDCversionNames
+        [
+            this->coeffs().lookupOrDefault
+            (
+                "version",
+                word(EDCversionNames[EDCdefaultVersion])
+            )
+        ]
+    ),
+    C1_(this->coeffs().lookupOrDefault("C1", 0.05774)),
+    C2_(this->coeffs().lookupOrDefault("C2", 0.5)),
+    Cgamma_(this->coeffs().lookupOrDefault("Cgamma", 2.1377)),
+    Ctau_(this->coeffs().lookupOrDefault("Ctau", 0.4083)),
+    exp1_(this->coeffs().lookupOrDefault("exp1", EDCexp1[int(version_)])),
+    exp2_(this->coeffs().lookupOrDefault("exp2", EDCexp2[int(version_)])),
+    kappa_
+    (
+        IOobject
+        (
+            IOobject::groupName(typeName + ":kappa", phaseName),
+            mesh.time().timeName(),
+            mesh,
+            IOobject::NO_READ,
+            IOobject::AUTO_WRITE
+        ),
+        mesh,
+        dimensionedScalar("kappa", dimless, 0)
+    )
+{}
+
+
+// * * * * * * * * * * * * * * * * Destructor  * * * * * * * * * * * * * * * //
+
+template<class Type>
+Foam::combustionModels::EDC<Type>::~EDC()
+{}
+
+
+// * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * * //
+
+template<class Type>
+void Foam::combustionModels::EDC<Type>::correct()
+{
+    if (this->active())
+    {
+        tmp<volScalarField> tepsilon(this->turbulence().epsilon());
+        const volScalarField& epsilon = tepsilon();
+
+        tmp<volScalarField> tmu(this->turbulence().mu());
+        const volScalarField& mu = tmu();
+
+        tmp<volScalarField> tk(this->turbulence().k());
+        const volScalarField& k = tk();
+
+        tmp<volScalarField> trho(this->rho());
+        const volScalarField& rho = trho();
+
+        scalarField tauStar(epsilon.size(), 0);
+
+        if (version_ == EDCversions::v2016)
+        {
+            tmp<volScalarField> ttc(this->chemistryPtr_->tc());
+            const volScalarField& tc = ttc();
+
+            forAll(tauStar, i)
+            {
+                const scalar nu = mu[i]/(rho[i] + SMALL);
+
+                const scalar Da =
+                    max(min(sqrt(nu/(epsilon[i] + SMALL))/tc[i], 10), 1e-10);
+
+                const scalar ReT = sqr(k[i])/(nu*epsilon[i] + SMALL);
+                const scalar CtauI = min(C1_/(Da*sqrt(ReT + 1)), 2.1377);
+
+                const scalar CgammaI =
+                    max(min(C2_*sqrt(Da*(ReT + 1)), 5), 0.4082);
+
+                const scalar gammaL =
+                    CgammaI*pow025(nu*epsilon[i]/(sqr(k[i]) + SMALL));
+
+                tauStar[i] = CtauI*sqrt(nu/(epsilon[i] + SMALL));
+
+                if (gammaL >= 1)
+                {
+                    kappa_[i] = 1;
+                }
+                else
+                {
+                    kappa_[i] =
+                        max
+                        (
+                            min
+                            (
+                                pow(gammaL, exp1_)/(1 - pow(gammaL, exp2_)),
+                                1
+                            ),
+                            0
+                        );
+                }
+            }
+        }
+        else
+        {
+            forAll(tauStar, i)
+            {
+                const scalar nu = mu[i]/(rho[i] + SMALL);
+                const scalar gammaL =
+                    Cgamma_*pow025(nu*epsilon[i]/(sqr(k[i]) + SMALL));
+
+                tauStar[i] = Ctau_*sqrt(nu/(epsilon[i] + SMALL));
+                if (gammaL >= 1)
+                {
+                    kappa_[i] = 1;
+                }
+                else
+                {
+                    kappa_[i] =
+                        max
+                        (
+                            min
+                            (
+                                pow(gammaL, exp1_)/(1 - pow(gammaL, exp2_)),
+                                1
+                            ),
+                            0
+                        );
+                }
+            }
+        }
+
+        this->chemistryPtr_->solve(tauStar);
+    }
+}
+
+
+template<class Type>
+Foam::tmp<Foam::fvScalarMatrix>
+Foam::combustionModels::EDC<Type>::R(volScalarField& Y) const
+{
+    return kappa_*laminar<Type>::R(Y);
+}
+
+
+template<class Type>
+Foam::tmp<Foam::volScalarField>
+Foam::combustionModels::EDC<Type>::Qdot() const
+{
+    tmp<volScalarField> tQdot
+    (
+        new volScalarField
+        (
+            IOobject
+            (
+                IOobject::groupName(typeName + ":Qdot", this->phaseName_),
+                this->mesh().time().timeName(),
+                this->mesh(),
+                IOobject::NO_READ,
+                IOobject::NO_WRITE,
+                false
+            ),
+            this->mesh(),
+            dimensionedScalar("Qdot", dimEnergy/dimVolume/dimTime, 0)
+        )
+    );
+
+    if (this->active())
+    {
+        tQdot.ref() = kappa_*this->chemistryPtr_->Qdot();
+    }
+
+    return tQdot;
+}
+
+
+template<class Type>
+bool Foam::combustionModels::EDC<Type>::read()
+{
+    if (Type::read())
+    {
+        version_ =
+        (
+            EDCversionNames
+            [
+                this->coeffs().lookupOrDefault
+                (
+                    "version",
+                    word(EDCversionNames[EDCdefaultVersion])
+                )
+            ]
+        );
+        C1_ = this->coeffs().lookupOrDefault("C1", 0.05774);
+        C2_ = this->coeffs().lookupOrDefault("C2", 0.5);
+        Cgamma_ = this->coeffs().lookupOrDefault("Cgamma", 2.1377);
+        Ctau_ = this->coeffs().lookupOrDefault("Ctau", 0.4083);
+        exp1_ = this->coeffs().lookupOrDefault("exp1", EDCexp1[int(version_)]);
+        exp2_ = this->coeffs().lookupOrDefault("exp2", EDCexp2[int(version_)]);
+
+        return true;
+    }
+    else
+    {
+        return false;
+    }
+}
+
+
+// ************************************************************************* //

src/combustionModels/EDC/EDC.H

@@ -0,0 +1,214 @@
+/*---------------------------------------------------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     |
+    \\  /    A nd           | Copyright (C) 2017 OpenFOAM Foundation
+     \\/     M anipulation  |
+-------------------------------------------------------------------------------
+License
+    This file is part of OpenFOAM.
+
+    OpenFOAM is free software: you can redistribute it and/or modify it
+    under the terms of the GNU General Public License as published by
+    the Free Software Foundation, either version 3 of the License, or
+    (at your option) any later version.
+
+    OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
+    ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
+    FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
+    for more details.
+
+    You should have received a copy of the GNU General Public License
+    along with OpenFOAM.  If not, see <http://www.gnu.org/licenses/>.
+
+Class
+    Foam::combustionModels::EDC
+
+Description
+    Eddy Dissipation Concept (EDC) turbulent combustion model.
+
+    This model considers that the reaction occurs in the regions of the flow
+    where the dissipation of turbulence kinetic energy takes place (fine
+    structures). The mass fraction of the fine structures and the mean residence
+    time are provided by an energy cascade model.
+
+    There are many versions and developments of the EDC model, 4 of which are
+    currently supported in this implementation: v1981, v1996, v2005 and
+    v2016.  The model variant is selected using the optional \c version entry in
+    the \c EDCCoeffs dictionary, \eg
+
+    \verbatim
+        EDCCoeffs
+        {
+            version v2016;
+        }
+    \endverbatim
+
+    The default version is \c v2015 if the \c version entry is not specified.
+
+    Model versions and references:
+    \verbatim
+        Version v2005:
+
+            Cgamma = 2.1377
+            Ctau = 0.4083
+            kappa = gammaL^exp1 / (1 - gammaL^exp2),
+
+            where exp1 = 2, and exp2 = 2.
+
+            Magnussen, B. F. (2005, June).
+            The Eddy Dissipation Concept -
+            A Bridge Between Science and Technology.
+            In ECCOMAS thematic conference on computational combustion
+            (pp. 21-24).
+
+        Version v1981:
+
+            Changes coefficients exp1 = 3 and exp2 = 3
+
+            Magnussen, B. (1981, January).
+            On the structure of turbulence and a generalized
+            eddy dissipation concept for chemical reaction in turbulent flow.
+            In 19th Aerospace Sciences Meeting (p. 42).
+
+        Version v1996:
+
+            Changes coefficients exp1 = 2 and exp2 = 3
+
+            Gran, I. R., & Magnussen, B. F. (1996).
+            A numerical study of a bluff-body stabilized diffusion flame.
+            Part 2. Influence of combustion modeling and finite-rate chemistry.
+            Combustion Science and Technology, 119(1-6), 191-217.
+
+        Version v2016:
+
+            Use local constants computed from the turbulent Da and Re numbers.
+
+            Parente, A., Malik, M. R., Contino, F., Cuoci, A., & Dally, B. B.
+            (2016).
+            Extension of the Eddy Dissipation Concept for
+            turbulence/chemistry interactions to MILD combustion.
+            Fuel, 163, 98-111.
+    \endverbatim
+
+SourceFiles
+    EDC.C
+
+\*---------------------------------------------------------------------------*/
+
+#ifndef EDC_H
+#define EDC_H
+
+#include "../laminar/laminar.H"
+#include "NamedEnum.H"
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+namespace Foam
+{
+namespace combustionModels
+{
+
+//- EDC model versions
+enum class EDCversions
+{
+    v1981,
+    v1996,
+    v2005,
+    v2016
+};
+
+extern const NamedEnum<EDCversions, 4> EDCversionNames;
+extern const EDCversions EDCdefaultVersion;
+
+const scalar EDCexp1[] = {3, 2, 2, 2};
+const scalar EDCexp2[] = {3, 3, 2, 2};
+
+/*---------------------------------------------------------------------------*\
+                            Class EDC Declaration
+\*---------------------------------------------------------------------------*/
+
+template<class Type>
+class EDC
+:
+    public laminar<Type>
+{
+    // Private data
+
+        //- The selected model version
+        EDCversions version_;
+
+        scalar C1_;
+        scalar C2_;
+        scalar Cgamma_;
+        scalar Ctau_;
+        scalar exp1_;
+        scalar exp2_;
+
+        //- Mixing parameter
+        volScalarField kappa_;
+
+
+    // Private Member Functions
+
+        //- Disallow copy construct
+        EDC(const EDC&);
+
+        //- Disallow default bitwise assignment
+        void operator=(const EDC&);
+
+
+public:
+
+    //- Runtime type information
+    TypeName("EDC");
+
+
+    // Constructors
+
+        //- Construct from components
+        EDC
+        (
+            const word& modelType,
+            const fvMesh& mesh,
+            const word& combustionProperties,
+            const word& phaseName
+        );
+
+
+    //- Destructor
+    virtual ~EDC();
+
+
+    // Member Functions
+
+        //- Correct combustion rate
+        virtual void correct();
+
+        //- Fuel consumption rate matrix.
+        virtual tmp<fvScalarMatrix> R(volScalarField& Y) const;
+
+        //- Heat release rate [kg/m/s3]
+        virtual tmp<volScalarField> Qdot() const;
+
+        //- Update properties from given dictionary
+        virtual bool read();
+};
+
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+} // End namespace combustionModels
+} // End namespace Foam
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+#ifdef NoRepository
+    #include "EDC.C"
+#endif
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+#endif
+
+// ************************************************************************* //

src/combustionModels/EDC/EDCs.C

@@ -0,0 +1,61 @@
+/*---------------------------------------------------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     |
+    \\  /    A nd           | Copyright (C) 2017 OpenFOAM Foundation
+     \\/     M anipulation  |
+-------------------------------------------------------------------------------
+License
+    This file is part of OpenFOAM.
+
+    OpenFOAM is free software: you can redistribute it and/or modify it
+    under the terms of the GNU General Public License as published by
+    the Free Software Foundation, either version 3 of the License, or
+    (at your option) any later version.
+
+    OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
+    ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
+    FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
+    for more details.
+
+    You should have received a copy of the GNU General Public License
+    along with OpenFOAM.  If not, see <http://www.gnu.org/licenses/>.
+
+\*---------------------------------------------------------------------------*/
+
+#include "makeCombustionTypes.H"
+
+#include "psiChemistryCombustion.H"
+#include "rhoChemistryCombustion.H"
+#include "EDC.H"
+
+// * * * * * * * * * * * * * Static Member Data  * * * * * * * * * * * * * * //
+
+template<>
+const char* Foam::NamedEnum
+<
+    Foam::combustionModels::EDCversions,
+    4
+>::names[] =
+{
+    "v1981",
+    "v1996",
+    "v2005",
+    "v2016"
+};
+
+const Foam::NamedEnum<Foam::combustionModels::EDCversions, 4>
+    Foam::combustionModels::EDCversionNames;
+
+const Foam::combustionModels::EDCversions
+Foam::combustionModels::EDCdefaultVersion
+(
+    Foam::combustionModels::EDCversions::v2005
+);
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+makeCombustionTypes(EDC, psiChemistryCombustion, psiCombustionModel);
+makeCombustionTypes(EDC, rhoChemistryCombustion, rhoCombustionModel);
+
+// ************************************************************************* //

src/combustionModels/Make/files

@@ -14,6 +14,7 @@ diffusion/diffusions.C
 infinitelyFastChemistry/infinitelyFastChemistrys.C
 
 PaSR/PaSRs.C
+EDC/EDCs.C
 
 laminar/laminars.C
 

src/functionObjects/utilities/systemCall/Make/files

@@ -1,4 +0,0 @@
-systemCall.C
-systemCallFunctionObject.C
-
-LIB = $(FOAM_LIBBIN)/libsystemCall

src/rigidBodyMeshMotion/rigidBodyMeshMotion/rigidBodyMeshMotion.C

@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2016 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2016-2017 OpenFOAM Foundation
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -30,6 +30,7 @@ License
 #include "pointConstraints.H"
 #include "uniformDimensionedFields.H"
 #include "forces.H"
+#include "OneConstant.H"
 #include "mathematicalConstants.H"
 
 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
@@ -115,6 +116,7 @@ Foam::rigidBodyMeshMotion::rigidBodyMeshMotion
     test_(coeffDict().lookupOrDefault<Switch>("test", false)),
     rhoInf_(1.0),
     rhoName_(coeffDict().lookupOrDefault<word>("rho", "rho")),
+    ramp_(nullptr),
     curTimeIndex_(-1)
 {
     if (rhoName_ == "rhoInf")
@@ -122,6 +124,15 @@ Foam::rigidBodyMeshMotion::rigidBodyMeshMotion
         rhoInf_ = readScalar(coeffDict().lookup("rhoInf"));
     }
 
+    if (coeffDict().found("ramp"))
+    {
+        ramp_ = Function1<scalar>::New("ramp", coeffDict());
+    }
+    else
+    {
+        ramp_ = new Function1Types::OneConstant<scalar>("ramp");
+    }
+
     const dictionary& bodiesDict = coeffDict().subDict("bodies");
 
     forAllConstIter(IDLList<entry>, bodiesDict, iter)
@@ -232,10 +243,12 @@ void Foam::rigidBodyMeshMotion::solve()
         curTimeIndex_ = this->db().time().timeIndex();
     }
 
+    const scalar ramp = ramp_->value(t.value());
+
     if (db().foundObject<uniformDimensionedVectorField>("g"))
     {
         model_.g() =
-            db().lookupObject<uniformDimensionedVectorField>("g").value();
+            ramp*db().lookupObject<uniformDimensionedVectorField>("g").value();
     }
 
     if (test_)
@@ -270,7 +283,7 @@ void Foam::rigidBodyMeshMotion::solve()
             functionObjects::forces f("forces", db(), forcesDict);
             f.calcForcesMoment();
 
-            fx[bodyID] = spatialVector(f.momentEff(), f.forceEff());
+            fx[bodyID] = ramp*spatialVector(f.momentEff(), f.forceEff());
         }
 
         model_.solve

src/rigidBodyMeshMotion/rigidBodyMeshMotion/rigidBodyMeshMotion.H

@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2016 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2016-2017 OpenFOAM Foundation
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -40,6 +40,7 @@ SourceFiles
 
 #include "displacementMotionSolver.H"
 #include "rigidBodyMotion.H"
+#include "ramp.H"
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
@@ -116,6 +117,9 @@ class rigidBodyMeshMotion
         //  as rhoInf
         word rhoName_;
 
+        //- Ramp the forces according to the specified function and period
+        autoPtr<Function1<scalar>> ramp_;
+
         //- Current time index (used for updating)
         label curTimeIndex_;
 

src/thermophysicalModels/chemistryModel/chemistryModel/TDACChemistryModel/TDACChemistryModel.C

@@ -616,6 +616,8 @@ Foam::scalar Foam::TDACChemistryModel<CompType, ThermoType>::solve
     scalar solveChemistryCpuTime_ = 0;
     scalar searchISATCpuTime_ = 0;
 
+    this->resetTabulationResults();
+
     // Average number of active species
     scalar nActiveSpecies = 0;
     scalar nAvg = 0;

src/thermophysicalModels/chemistryModel/chemistryModel/TDACChemistryModel/TDACChemistryModel.H

@@ -283,6 +283,8 @@ public:
             void setTabulationResultsGrow(const label celli);
 
             void setTabulationResultsRetrieve(const label celli);
+
+            inline void resetTabulationResults();
 };
 
 

src/thermophysicalModels/chemistryModel/chemistryModel/TDACChemistryModel/TDACChemistryModelI.H

@@ -165,4 +165,14 @@ Foam::TDACChemistryModel<CompType, ThermoType>::specieComp()
 }
 
 
+template<class CompType, class ThermoType>
+void Foam::TDACChemistryModel<CompType, ThermoType>::resetTabulationResults()
+{
+    forAll(tabulationResults_, tabi)
+    {
+        tabulationResults_[tabi] = 2;
+    }
+}
+
+
 // ************************************************************************* //

tutorials/combustion/reactingFoam/RAS/DLR_A_LTS/0/CH4

@@ -0,0 +1,67 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  dev                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    location    "0";
+    object      CH4;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 0 0 0 0 0 0];
+
+internalField   uniform 0;
+
+boundaryField
+{
+    inletfuel
+    {
+        type            fixedValue;
+        value           uniform 0.245642;
+    }
+    inletair
+    {
+        type            fixedValue;
+        value           uniform 0;
+    }
+    outlet
+    {
+        type            inletOutlet;
+        inletValue      uniform 0;
+        value           uniform 0;
+    }
+    axis
+    {
+        type            empty;
+    }
+    leftside
+    {
+        type            zeroGradient;
+    }
+    burnerwall
+    {
+        type            zeroGradient;
+    }
+    burnertip
+    {
+        type            zeroGradient;
+    }
+    front
+    {
+        type            wedge;
+    }
+    back
+    {
+        type            wedge;
+    }
+}
+
+
+// ************************************************************************* //

tutorials/combustion/reactingFoam/RAS/DLR_A_LTS/0/G

@@ -0,0 +1,42 @@
+/*--------------------------------*- C++ -*----------------------------------* \
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  dev                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    object      G;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [1 0 -3 0 0 0 0];
+
+internalField   uniform 0;
+
+boundaryField
+{
+    ".*"
+    {
+        type            MarshakRadiation;
+        T               T;
+        emissivityMode  lookup;
+        emissivity      uniform 1.0;
+        value           uniform 0;
+    }
+
+    front
+    {
+        type            wedge;
+    }
+    back
+    {
+        type            wedge;
+    }
+}
+
+// ************************************************************************* //

tutorials/combustion/reactingFoam/RAS/DLR_A_LTS/0/H2

@@ -0,0 +1,67 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  dev                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    location    "0";
+    object      H2;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 0 0 0 0 0 0];
+
+internalField   uniform 0;
+
+boundaryField
+{
+    inletfuel
+    {
+        type            fixedValue;
+        value           uniform 0.046496;
+    }
+    inletair
+    {
+        type            fixedValue;
+        value           uniform 0;
+    }
+    outlet
+    {
+        type            inletOutlet;
+        inletValue      uniform 0;
+        value           uniform 0;
+    }
+    axis
+    {
+        type            empty;
+    }
+    leftside
+    {
+        type            zeroGradient;
+    }
+    burnerwall
+    {
+        type            zeroGradient;
+    }
+    burnertip
+    {
+        type            zeroGradient;
+    }
+    front
+    {
+        type            wedge;
+    }
+    back
+    {
+        type            wedge;
+    }
+}
+
+
+// ************************************************************************* //

tutorials/combustion/reactingFoam/RAS/DLR_A_LTS/0/H2O

@@ -0,0 +1,67 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  dev                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    location    "0";
+    object      H2;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 0 0 0 0 0 0];
+
+internalField   uniform 0.005008;
+
+boundaryField
+{
+    inletfuel
+    {
+        type            fixedValue;
+        value           uniform 0;
+    }
+    inletair
+    {
+        type            fixedValue;
+        value           uniform 0.005008;
+    }
+    outlet
+    {
+        type            inletOutlet;
+        inletValue      uniform 0.005008;
+        value           uniform 0.005008;
+    }
+    axis
+    {
+        type            empty;
+    }
+    leftside
+    {
+        type            zeroGradient;
+    }
+    burnerwall
+    {
+        type            zeroGradient;
+    }
+    burnertip
+    {
+        type            zeroGradient;
+    }
+    front
+    {
+        type            wedge;
+    }
+    back
+    {
+        type            wedge;
+    }
+}
+
+
+// ************************************************************************* //

tutorials/combustion/reactingFoam/RAS/DLR_A_LTS/0/N2

@@ -0,0 +1,67 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  dev                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    location    "0";
+    object      N2;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 0 0 0 0 0 0];
+
+internalField   uniform 0.763149;
+
+boundaryField
+{
+    inletfuel
+    {
+        type            fixedValue;
+        value           uniform 0.707861;
+    }
+    inletair
+    {
+        type            fixedValue;
+        value           uniform 0.763149;
+    }
+    outlet
+    {
+        type            inletOutlet;
+        inletValue      uniform 0.763149;
+        value           uniform 0.763149;
+    }
+    axis
+    {
+        type            empty;
+    }
+    leftside
+    {
+        type            zeroGradient;
+    }
+    burnerwall
+    {
+        type            zeroGradient;
+    }
+    burnertip
+    {
+        type            zeroGradient;
+    }
+    front
+    {
+        type            wedge;
+    }
+    back
+    {
+        type            wedge;
+    }
+}
+
+
+// ************************************************************************* //

tutorials/combustion/reactingFoam/RAS/DLR_A_LTS/0/O2

@@ -0,0 +1,67 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  dev                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    location    "0";
+    object      O2;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 0 0 0 0 0 0];
+
+internalField  uniform 0.231843;
+
+boundaryField
+{
+    inletfuel
+    {
+        type            fixedValue;
+        value           uniform 0;
+    }
+    inletair
+    {
+        type            fixedValue;
+        value           uniform 0.231843;
+    }
+        outlet
+    {
+        type            inletOutlet;
+        inletValue      uniform 0.231843;
+        value           uniform 0.231843;
+    }
+    axis
+    {
+        type            empty;
+    }
+    leftside
+    {
+        type            zeroGradient;
+    }
+    burnerwall
+    {
+        type            zeroGradient;
+    }
+    burnertip
+    {
+        type            zeroGradient;
+    }
+    front
+    {
+        type            wedge;
+    }
+    back
+    {
+        type            wedge;
+    }
+}
+
+
+// ************************************************************************* //

tutorials/combustion/reactingFoam/RAS/DLR_A_LTS/0/T.orig

@@ -0,0 +1,67 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  dev                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    location    "0";
+    object      T;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 0 0 1 0 0 0];
+
+internalField   uniform 292;
+
+boundaryField
+{
+    inletfuel
+    {
+        type            fixedValue;
+        value           uniform 292;
+    }
+    inletair
+    {
+        type            fixedValue;
+        value           uniform 292;
+    }
+    outlet
+    {
+        type            inletOutlet;
+        inletValue      uniform 292;
+        value           uniform 292;
+    }
+    axis
+    {
+        type            empty;
+    }
+    leftside
+    {
+        type            zeroGradient;
+    }
+    burnerwall
+    {
+        type            zeroGradient;
+    }
+    burnertip
+    {
+        type            zeroGradient;
+    }
+    front
+    {
+        type            wedge;
+    }
+    back
+    {
+        type            wedge;
+    }
+}
+
+
+// ************************************************************************* //

tutorials/combustion/reactingFoam/RAS/DLR_A_LTS/0/U

@@ -0,0 +1,69 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  dev                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volVectorField;
+    location    "0";
+    object      U;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 1 -1 0 0 0 0];
+
+internalField   uniform (0 0 0.3);
+
+boundaryField
+{
+    inletfuel
+    {
+        type            fixedValue;
+        value           uniform (0 0 42.2);
+    }
+    inletair
+    {
+        type            fixedValue;
+        value           uniform (0 0 0.3);
+    }
+    outlet
+    {
+        type            pressureInletOutletVelocity;
+        value           uniform (0 0 0.3);
+    }
+    axis
+    {
+        type            empty;
+    }
+    leftside
+    {
+        type            fixedValue;
+        value           uniform (0 0 0.3);
+    }
+    burnerwall
+    {
+        type            fixedValue;
+        value           uniform (0 0 0);
+    }
+    burnertip
+    {
+        type            fixedValue;
+        value           uniform (0 0 0);
+    }
+    front
+    {
+        type            wedge;
+    }
+    back
+    {
+        type            wedge;
+    }
+}
+
+
+// ************************************************************************* //

tutorials/combustion/reactingFoam/RAS/DLR_A_LTS/0/Ydefault

@@ -0,0 +1,67 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  dev                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    location    "0";
+    object      Ydefault;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 0 0 0 0 0 0];
+
+internalField   uniform 0;
+
+boundaryField
+{
+    inletfuel
+    {
+        type            fixedValue;
+        value           uniform 0;
+    }
+    inletair
+    {
+        type            fixedValue;
+        value           uniform 0;
+    }
+    outlet
+    {
+        type            inletOutlet;
+        inletValue      uniform 0;
+        value           uniform 0;
+    }
+    axis
+    {
+        type            empty;
+    }
+    leftside
+    {
+        type            zeroGradient;
+    }
+    burnerwall
+    {
+        type            zeroGradient;
+    }
+    burnertip
+    {
+        type            zeroGradient;
+    }
+    front
+    {
+        type            wedge;
+    }
+    back
+    {
+        type            wedge;
+    }
+}
+
+
+// ************************************************************************* //

tutorials/combustion/reactingFoam/RAS/DLR_A_LTS/0/alphat

@@ -0,0 +1,72 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  dev                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    location    "0";
+    object      alphat;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [1 -1 -1 0 0 0 0];
+
+internalField   uniform 0;
+
+boundaryField
+{
+    inletfuel
+    {
+        type            calculated;
+        value           uniform 0;
+    }
+    inletair
+    {
+        type            calculated;
+        value           uniform 0;
+    }
+    outlet
+    {
+        type            calculated;
+        value           uniform 0;
+    }
+    axis
+    {
+        type            empty;
+    }
+    leftside
+    {
+        type            compressible::alphatWallFunction;
+        Prt             0.85;
+        value           uniform 0;
+    }
+    burnerwall
+    {
+        type            compressible::alphatWallFunction;
+        Prt             0.85;
+        value           uniform 0;
+    }
+    burnertip
+    {
+        type            compressible::alphatWallFunction;
+        Prt             0.85;
+        value           uniform 0;
+    }
+    front
+    {
+        type            wedge;
+    }
+    back
+    {
+        type            wedge;
+    }
+}
+
+
+// ************************************************************************* //

tutorials/combustion/reactingFoam/RAS/DLR_A_LTS/0/epsilon

@@ -0,0 +1,85 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  dev                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    location    "0";
+    object      epsilon;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 2 -3 0 0 0 0];
+
+internalField   uniform 200;
+
+boundaryField
+{
+    inletfuel
+    {
+        type            turbulentMixingLengthDissipationRateInlet;
+        mixingLength    0.005;
+        phi             phi;
+        k               k;
+        value           uniform 200;
+    }
+    inletair
+    {
+        type            turbulentMixingLengthDissipationRateInlet;
+        mixingLength    0.005;
+        phi             phi;
+        k               k;
+        value           uniform 200;
+    }
+    outlet
+    {
+        type            inletOutlet;
+        inletValue      uniform 200;
+        value           uniform 200;
+    }
+    axis
+    {
+        type            empty;
+    }
+    leftside
+    {
+        type            epsilonWallFunction;
+        Cmu             0.09;
+        kappa           0.41;
+        E               9.8;
+        value           uniform 200;
+    }
+    burnerwall
+    {
+        type            epsilonWallFunction;
+        Cmu             0.09;
+        kappa           0.41;
+        E               9.8;
+        value           uniform 200;
+    }
+    burnertip
+    {
+        type            epsilonWallFunction;
+        Cmu             0.09;
+        kappa           0.41;
+        E               9.8;
+        value           uniform 200;
+    }
+    front
+    {
+        type            wedge;
+    }
+    back
+    {
+        type            wedge;
+    }
+}
+
+
+// ************************************************************************* //

tutorials/combustion/reactingFoam/RAS/DLR_A_LTS/0/k

@@ -0,0 +1,72 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  dev                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    location    "0";
+    object      k;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 2 -2 0 0 0 0];
+
+internalField   uniform 1;
+
+boundaryField
+{
+    inletfuel
+    {
+        type            turbulentIntensityKineticEnergyInlet;
+        intensity       0.08;
+        value           uniform 1;
+    }
+    inletair
+    {
+        type            turbulentIntensityKineticEnergyInlet;
+        intensity       0.02;
+        value           uniform 1;
+    }
+    outlet
+    {
+        type            inletOutlet;
+        inletValue      uniform 1;
+        value           uniform 1;
+    }
+    axis
+    {
+        type            empty;
+    }
+    leftside
+    {
+        type            kqRWallFunction;
+        value           uniform 1;
+    }
+    burnerwall
+    {
+        type            kqRWallFunction;
+        value           uniform 1;
+    }
+    burnertip
+    {
+        type            kqRWallFunction;
+        value           uniform 1;
+    }
+    front
+    {
+        type            wedge;
+    }
+    back
+    {
+        type            wedge;
+    }
+}
+
+
+// ************************************************************************* //

tutorials/combustion/reactingFoam/RAS/DLR_A_LTS/0/nut

@@ -0,0 +1,78 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  dev                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    location    "0";
+    object      mut;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 2 -1 0 0 0 0];
+
+internalField   uniform 0;
+
+boundaryField
+{
+    inletfuel
+    {
+        type            calculated;
+        value           uniform 0;
+    }
+    inletair
+    {
+        type            calculated;
+        value           uniform 0;
+    }
+    outlet
+    {
+        type            calculated;
+        value           uniform 0;
+    }
+    axis
+    {
+        type            empty;
+    }
+    leftside
+    {
+        type            nutkWallFunction;
+        Cmu             0.09;
+        kappa           0.41;
+        E               9.8;
+        value           uniform 0;
+    }
+    burnerwall
+    {
+        type            nutkWallFunction;
+        Cmu             0.09;
+        kappa           0.41;
+        E               9.8;
+        value           uniform 0;
+    }
+    burnertip
+    {
+        type            nutkWallFunction;
+        Cmu             0.09;
+        kappa           0.41;
+        E               9.8;
+        value           uniform 0;
+    }
+    front
+    {
+        type            wedge;
+    }
+    back
+    {
+        type            wedge;
+    }
+}
+
+
+// ************************************************************************* //

tutorials/combustion/reactingFoam/RAS/DLR_A_LTS/0/p

@@ -0,0 +1,64 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  dev                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    location    "0";
+    object      p;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [1 -1 -2 0 0 0 0];
+
+internalField   uniform 101325;
+
+boundaryField
+{
+    inletfuel
+    {
+        type            zeroGradient;
+    }
+    inletair
+    {
+        type            zeroGradient;
+    }
+    outlet
+    {
+        type            totalPressure;
+        p0              $internalField;
+    }
+    axis
+    {
+        type            empty;
+    }
+    leftside
+    {
+        type            zeroGradient;
+    }
+    burnerwall
+    {
+        type            zeroGradient;
+    }
+    burnertip
+    {
+        type            zeroGradient;
+    }
+    front
+    {
+        type            wedge;
+    }
+    back
+    {
+        type            wedge;
+    }
+}
+
+
+// ************************************************************************* //

tutorials/combustion/reactingFoam/RAS/DLR_A_LTS/Allclean

@@ -0,0 +1,11 @@
+#!/bin/sh
+cd ${0%/*} || exit 1    # Run from this directory
+
+# Source tutorial clean functions
+. $WM_PROJECT_DIR/bin/tools/CleanFunctions
+
+cleanCase
+
+rm -f 0/T
+
+#------------------------------------------------------------------------------

tutorials/combustion/reactingFoam/RAS/DLR_A_LTS/Allrun

@@ -0,0 +1,23 @@
+#!/bin/sh
+cd ${0%/*} || exit 1    # Run from this directory
+
+# Source tutorial run functions
+. $WM_PROJECT_DIR/bin/tools/RunFunctions
+
+# Set application name
+application=`getApplication`
+
+rm -f 0/T
+cp 0/T.orig 0/T
+
+runApplication chemkinToFoam \
+               chemkin/grimech30.dat chemkin/thermo30.dat chemkin/transportProperties \
+               constant/reactionsGRI constant/thermo.compressibleGasGRI
+
+runApplication blockMesh
+runApplication setFields
+runApplication decomposePar -force
+runParallel $application
+runApplication reconstructPar
+
+#------------------------------------------------------------------------------

tutorials/combustion/reactingFoam/RAS/DLR_A_LTS/chemkin/grimech30.dat

@@ -0,0 +1,328 @@
+! GRI-Mech Version 3.0 7/30/99  CHEMKIN-II format
+! See README30 file at anonymous FTP site unix.sri.com, directory gri;
+! WorldWideWeb home page http://www.me.berkeley.edu/gri_mech/ or
+! through http://www.gri.org , under 'Basic  Research',
+! for additional information, contacts, and disclaimer
+ELEMENTS
+O  H  C  N  AR
+END
+SPECIES
+H2      H       O       O2      OH      H2O     HO2     H2O2
+C       CH      CH2     CH2(S)  CH3     CH4     CO      CO2
+HCO     CH2O    CH2OH   CH3O    CH3OH   C2H     C2H2    C2H3
+C2H4    C2H5    C2H6    HCCO    CH2CO   HCCOH   N2      AR
+C3H7    C3H8    CH2CHO  CH3CHO
+END
+!THERMO
+! Insert GRI-Mech thermodynamics here or use in default file
+!END
+REACTIONS
+2O+M<=>O2+M                              1.200E+17   -1.000        .00
+H2/ 2.40/ H2O/15.40/ CH4/ 2.00/ CO/ 1.75/ CO2/ 3.60/ C2H6/ 3.00/ AR/  .83/
+O+H+M<=>OH+M                             5.000E+17   -1.000        .00
+H2/2.00/ H2O/6.00/ CH4/2.00/ CO/1.50/ CO2/2.00/ C2H6/3.00/ AR/ .70/
+O+H2<=>H+OH                              3.870E+04    2.700    6260.00
+O+HO2<=>OH+O2                            2.000E+13     .000        .00
+O+H2O2<=>OH+HO2                          9.630E+06    2.000    4000.00
+O+CH<=>H+CO                              5.700E+13     .000        .00
+O+CH2<=>H+HCO                            8.000E+13     .000        .00
+O+CH2(S)<=>H2+CO                         1.500E+13     .000        .00
+O+CH2(S)<=>H+HCO                         1.500E+13     .000        .00
+O+CH3<=>H+CH2O                           5.060E+13     .000        .00
+O+CH4<=>OH+CH3                           1.020E+09    1.500    8600.00
+O+CO(+M)<=>CO2(+M)                       1.800E+10     .000    2385.00
+   LOW/ 6.020E+14     .000    3000.00/
+H2/2.00/ O2/6.00/ H2O/6.00/ CH4/2.00/ CO/1.50/ CO2/3.50/ C2H6/3.00/ AR/ .50/
+O+HCO<=>OH+CO                            3.000E+13     .000        .00
+O+HCO<=>H+CO2                            3.000E+13     .000        .00
+O+CH2O<=>OH+HCO                          3.900E+13     .000    3540.00
+O+CH2OH<=>OH+CH2O                        1.000E+13     .000        .00
+O+CH3O<=>OH+CH2O                         1.000E+13     .000        .00
+O+CH3OH<=>OH+CH2OH                       3.880E+05    2.500    3100.00
+O+CH3OH<=>OH+CH3O                        1.300E+05    2.500    5000.00
+O+C2H<=>CH+CO                            5.000E+13     .000        .00
+O+C2H2<=>H+HCCO                          1.350E+07    2.000    1900.00
+O+C2H2<=>OH+C2H                          4.600E+19   -1.410   28950.00
+O+C2H2<=>CO+CH2                          6.940E+06    2.000    1900.00
+O+C2H3<=>H+CH2CO                         3.000E+13     .000        .00
+O+C2H4<=>CH3+HCO                         1.250E+07    1.830     220.00
+O+C2H5<=>CH3+CH2O                        2.240E+13     .000        .00
+O+C2H6<=>OH+C2H5                         8.980E+07    1.920    5690.00
+O+HCCO<=>H+2CO                           1.000E+14     .000        .00
+O+CH2CO<=>OH+HCCO                        1.000E+13     .000    8000.00
+O+CH2CO<=>CH2+CO2                        1.750E+12     .000    1350.00
+O2+CO<=>O+CO2                            2.500E+12     .000   47800.00
+O2+CH2O<=>HO2+HCO                        1.000E+14     .000   40000.00
+H+O2+M<=>HO2+M                           2.800E+18    -.860        .00
+O2/ .00/ H2O/ .00/ CO/ .75/ CO2/1.50/ C2H6/1.50/ N2/ .00/ AR/ .00/
+H+2O2<=>HO2+O2                           2.080E+19   -1.240        .00
+H+O2+H2O<=>HO2+H2O                       11.26E+18    -.760        .00
+H+O2+N2<=>HO2+N2                         2.600E+19   -1.240        .00
+H+O2+AR<=>HO2+AR                         7.000E+17    -.800        .00
+H+O2<=>O+OH                              2.650E+16    -.6707  17041.00
+2H+M<=>H2+M                              1.000E+18   -1.000        .00
+H2/ .00/ H2O/ .00/ CH4/2.00/ CO2/ .00/ C2H6/3.00/ AR/ .63/
+2H+H2<=>2H2                              9.000E+16    -.600        .00
+2H+H2O<=>H2+H2O                          6.000E+19   -1.250        .00
+2H+CO2<=>H2+CO2                          5.500E+20   -2.000        .00
+H+OH+M<=>H2O+M                           2.200E+22   -2.000        .00
+H2/ .73/ H2O/3.65/ CH4/2.00/ C2H6/3.00/ AR/ .38/
+H+HO2<=>O+H2O                            3.970E+12     .000     671.00
+H+HO2<=>O2+H2                            4.480E+13     .000    1068.00
+H+HO2<=>2OH                              0.840E+14     .000     635.00
+H+H2O2<=>HO2+H2                          1.210E+07    2.000    5200.00
+H+H2O2<=>OH+H2O                          1.000E+13     .000    3600.00
+H+CH<=>C+H2                              1.650E+14     .000        .00
+H+CH2(+M)<=>CH3(+M)                      6.000E+14     .000        .00
+     LOW  /  1.040E+26   -2.760   1600.00/
+     TROE/   .5620  91.00  5836.00  8552.00/
+H2/2.00/ H2O/6.00/ CH4/2.00/ CO/1.50/ CO2/2.00/ C2H6/3.00/ AR/ .70/
+H+CH2(S)<=>CH+H2                         3.000E+13     .000        .00
+H+CH3(+M)<=>CH4(+M)                      13.90E+15    -.534     536.00
+     LOW  /  2.620E+33   -4.760   2440.00/
+     TROE/   .7830   74.00  2941.00  6964.00 /
+H2/2.00/ H2O/6.00/ CH4/3.00/ CO/1.50/ CO2/2.00/ C2H6/3.00/ AR/ .70/
+H+CH4<=>CH3+H2                           6.600E+08    1.620   10840.00
+H+HCO(+M)<=>CH2O(+M)                     1.090E+12     .480    -260.00
+     LOW  /  2.470E+24   -2.570    425.00/
+     TROE/   .7824  271.00  2755.00  6570.00 /
+H2/2.00/ H2O/6.00/ CH4/2.00/ CO/1.50/ CO2/2.00/ C2H6/3.00/ AR/ .70/
+H+HCO<=>H2+CO                            7.340E+13     .000        .00
+H+CH2O(+M)<=>CH2OH(+M)                   5.400E+11     .454    3600.00
+     LOW  /  1.270E+32   -4.820   6530.00/
+     TROE/   .7187  103.00  1291.00  4160.00 /
+H2/2.00/ H2O/6.00/ CH4/2.00/ CO/1.50/ CO2/2.00/ C2H6/3.00/
+H+CH2O(+M)<=>CH3O(+M)                    5.400E+11     .454    2600.00
+     LOW  /  2.200E+30   -4.800   5560.00/
+     TROE/   .7580   94.00  1555.00  4200.00 /
+H2/2.00/ H2O/6.00/ CH4/2.00/ CO/1.50/ CO2/2.00/ C2H6/3.00/
+H+CH2O<=>HCO+H2                          5.740E+07    1.900    2742.00
+H+CH2OH(+M)<=>CH3OH(+M)                  1.055E+12     .500      86.00
+     LOW  /  4.360E+31   -4.650   5080.00/
+     TROE/   .600  100.00  90000.0  10000.0 /
+H2/2.00/ H2O/6.00/ CH4/2.00/ CO/1.50/ CO2/2.00/ C2H6/3.00/
+H+CH2OH<=>H2+CH2O                        2.000E+13     .000        .00
+H+CH2OH<=>OH+CH3                         1.650E+11     .650    -284.00
+H+CH2OH<=>CH2(S)+H2O                     3.280E+13    -.090     610.00
+H+CH3O(+M)<=>CH3OH(+M)                   2.430E+12     .515      50.00
+     LOW  /  4.660E+41   -7.440   14080.0/
+     TROE/   .700  100.00  90000.0 10000.00 /
+H2/2.00/ H2O/6.00/ CH4/2.00/ CO/1.50/ CO2/2.00/ C2H6/3.00/
+H+CH3O<=>H+CH2OH                         4.150E+07    1.630    1924.00
+H+CH3O<=>H2+CH2O                         2.000E+13     .000        .00
+H+CH3O<=>OH+CH3                          1.500E+12     .500    -110.00
+H+CH3O<=>CH2(S)+H2O                      2.620E+14    -.230    1070.00
+H+CH3OH<=>CH2OH+H2                       1.700E+07    2.100    4870.00
+H+CH3OH<=>CH3O+H2                        4.200E+06    2.100    4870.00
+H+C2H(+M)<=>C2H2(+M)                     1.000E+17   -1.000        .00
+     LOW  /  3.750E+33   -4.800   1900.00/
+     TROE/   .6464  132.00  1315.00  5566.00 /
+H2/2.00/ H2O/6.00/ CH4/2.00/ CO/1.50/ CO2/2.00/ C2H6/3.00/ AR/ .70/
+H+C2H2(+M)<=>C2H3(+M)                    5.600E+12     .000    2400.00
+     LOW  /  3.800E+40   -7.270   7220.00/
+     TROE/   .7507   98.50  1302.00  4167.00 /
+H2/2.00/ H2O/6.00/ CH4/2.00/ CO/1.50/ CO2/2.00/ C2H6/3.00/ AR/ .70/
+H+C2H3(+M)<=>C2H4(+M)                    6.080E+12     .270     280.00
+     LOW  /  1.400E+30   -3.860   3320.00/
+     TROE/   .7820  207.50  2663.00  6095.00 /
+H2/2.00/ H2O/6.00/ CH4/2.00/ CO/1.50/ CO2/2.00/ C2H6/3.00/ AR/ .70/
+H+C2H3<=>H2+C2H2                         3.000E+13     .000        .00
+H+C2H4(+M)<=>C2H5(+M)                    0.540E+12     .454    1820.00
+     LOW  /  0.600E+42   -7.620   6970.00/
+     TROE/   .9753  210.00   984.00  4374.00 /
+H2/2.00/ H2O/6.00/ CH4/2.00/ CO/1.50/ CO2/2.00/ C2H6/3.00/ AR/ .70/
+H+C2H4<=>C2H3+H2                         1.325E+06    2.530   12240.00
+H+C2H5(+M)<=>C2H6(+M)                    5.210E+17    -.990    1580.00
+     LOW  /  1.990E+41   -7.080   6685.00/
+     TROE/   .8422  125.00  2219.00  6882.00 /
+H2/2.00/ H2O/6.00/ CH4/2.00/ CO/1.50/ CO2/2.00/ C2H6/3.00/ AR/ .70/
+H+C2H5<=>H2+C2H4                         2.000E+12     .000        .00
+H+C2H6<=>C2H5+H2                         1.150E+08    1.900    7530.00
+H+HCCO<=>CH2(S)+CO                       1.000E+14     .000        .00
+H+CH2CO<=>HCCO+H2                        5.000E+13     .000    8000.00
+H+CH2CO<=>CH3+CO                         1.130E+13     .000    3428.00
+H+HCCOH<=>H+CH2CO                        1.000E+13     .000        .00
+H2+CO(+M)<=>CH2O(+M)                     4.300E+07    1.500   79600.00
+     LOW  /  5.070E+27   -3.420  84350.00/
+     TROE/   .9320  197.00  1540.00 10300.00 /
+H2/2.00/ H2O/6.00/ CH4/2.00/ CO/1.50/ CO2/2.00/ C2H6/3.00/ AR/ .70/
+OH+H2<=>H+H2O                            2.160E+08    1.510    3430.00
+2OH(+M)<=>H2O2(+M)                       7.400E+13    -.370        .00
+     LOW  /  2.300E+18    -.900  -1700.00/
+     TROE/   .7346   94.00  1756.00  5182.00 /
+H2/2.00/ H2O/6.00/ CH4/2.00/ CO/1.50/ CO2/2.00/ C2H6/3.00/ AR/ .70/
+2OH<=>O+H2O                              3.570E+04    2.400   -2110.00
+OH+HO2<=>O2+H2O                          1.450E+13     .000    -500.00
+ DUPLICATE
+OH+H2O2<=>HO2+H2O                        2.000E+12     .000     427.00
+ DUPLICATE
+OH+H2O2<=>HO2+H2O                        1.700E+18     .000   29410.00
+ DUPLICATE
+OH+C<=>H+CO                              5.000E+13     .000        .00
+OH+CH<=>H+HCO                            3.000E+13     .000        .00
+OH+CH2<=>H+CH2O                          2.000E+13     .000        .00
+OH+CH2<=>CH+H2O                          1.130E+07    2.000    3000.00
+OH+CH2(S)<=>H+CH2O                       3.000E+13     .000        .00
+OH+CH3(+M)<=>CH3OH(+M)                   2.790E+18   -1.430    1330.00
+     LOW  /  4.000E+36   -5.920   3140.00/
+     TROE/   .4120  195.0  5900.00  6394.00/
+H2/2.00/ H2O/6.00/ CH4/2.00/ CO/1.50/ CO2/2.00/ C2H6/3.00/
+OH+CH3<=>CH2+H2O                         5.600E+07    1.600    5420.00
+OH+CH3<=>CH2(S)+H2O                      6.440E+17   -1.340    1417.00
+OH+CH4<=>CH3+H2O                         1.000E+08    1.600    3120.00
+OH+CO<=>H+CO2                            4.760E+07    1.228      70.00
+OH+HCO<=>H2O+CO                          5.000E+13     .000        .00
+OH+CH2O<=>HCO+H2O                        3.430E+09    1.180    -447.00
+OH+CH2OH<=>H2O+CH2O                      5.000E+12     .000        .00
+OH+CH3O<=>H2O+CH2O                       5.000E+12     .000        .00
+OH+CH3OH<=>CH2OH+H2O                     1.440E+06    2.000    -840.00
+OH+CH3OH<=>CH3O+H2O                      6.300E+06    2.000    1500.00
+OH+C2H<=>H+HCCO                          2.000E+13     .000        .00
+OH+C2H2<=>H+CH2CO                        2.180E-04    4.500   -1000.00
+OH+C2H2<=>H+HCCOH                        5.040E+05    2.300   13500.00
+OH+C2H2<=>C2H+H2O                        3.370E+07    2.000   14000.00
+OH+C2H2<=>CH3+CO                         4.830E-04    4.000   -2000.00
+OH+C2H3<=>H2O+C2H2                       5.000E+12     .000        .00
+OH+C2H4<=>C2H3+H2O                       3.600E+06    2.000    2500.00
+OH+C2H6<=>C2H5+H2O                       3.540E+06    2.120     870.00
+OH+CH2CO<=>HCCO+H2O                      7.500E+12     .000    2000.00
+2HO2<=>O2+H2O2                           1.300E+11     .000   -1630.00
+ DUPLICATE
+2HO2<=>O2+H2O2                           4.200E+14     .000   12000.00
+ DUPLICATE
+HO2+CH2<=>OH+CH2O                        2.000E+13     .000        .00
+HO2+CH3<=>O2+CH4                         1.000E+12     .000        .00
+HO2+CH3<=>OH+CH3O                        3.780E+13     .000        .00
+HO2+CO<=>OH+CO2                          1.500E+14     .000   23600.00
+HO2+CH2O<=>HCO+H2O2                      5.600E+06    2.000   12000.00
+C+O2<=>O+CO                              5.800E+13     .000     576.00
+C+CH2<=>H+C2H                            5.000E+13     .000        .00
+C+CH3<=>H+C2H2                           5.000E+13     .000        .00
+CH+O2<=>O+HCO                            6.710E+13     .000        .00
+CH+H2<=>H+CH2                            1.080E+14     .000    3110.00
+CH+H2O<=>H+CH2O                          5.710E+12     .000    -755.00
+CH+CH2<=>H+C2H2                          4.000E+13     .000        .00
+CH+CH3<=>H+C2H3                          3.000E+13     .000        .00
+CH+CH4<=>H+C2H4                          6.000E+13     .000        .00
+CH+CO(+M)<=>HCCO(+M)                     5.000E+13     .000        .00
+     LOW  /  2.690E+28   -3.740   1936.00/
+     TROE/   .5757  237.00  1652.00  5069.00 /
+H2/2.00/ H2O/6.00/ CH4/2.00/ CO/1.50/ CO2/2.00/ C2H6/3.00/ AR/ .70/
+CH+CO2<=>HCO+CO                          1.900E+14     .000   15792.00
+CH+CH2O<=>H+CH2CO                        9.460E+13     .000    -515.00
+CH+HCCO<=>CO+C2H2                        5.000E+13     .000        .00
+CH2+O2=>OH+H+CO                          5.000E+12     .000    1500.00
+CH2+H2<=>H+CH3                           5.000E+05    2.000    7230.00
+2CH2<=>H2+C2H2                           1.600E+15     .000   11944.00
+CH2+CH3<=>H+C2H4                         4.000E+13     .000        .00
+CH2+CH4<=>2CH3                           2.460E+06    2.000    8270.00
+CH2+CO(+M)<=>CH2CO(+M)                   8.100E+11     .500    4510.00
+     LOW  /  2.690E+33   -5.110   7095.00/
+     TROE/   .5907  275.00  1226.00  5185.00 /
+H2/2.00/ H2O/6.00/ CH4/2.00/ CO/1.50/ CO2/2.00/ C2H6/3.00/ AR/ .70/
+CH2+HCCO<=>C2H3+CO                       3.000E+13     .000        .00
+CH2(S)+N2<=>CH2+N2                       1.500E+13     .000     600.00
+CH2(S)+AR<=>CH2+AR                       9.000E+12     .000     600.00
+CH2(S)+O2<=>H+OH+CO                      2.800E+13     .000        .00
+CH2(S)+O2<=>CO+H2O                       1.200E+13     .000        .00
+CH2(S)+H2<=>CH3+H                        7.000E+13     .000        .00
+CH2(S)+H2O(+M)<=>CH3OH(+M)               4.820E+17   -1.160    1145.00
+     LOW  /  1.880E+38   -6.360   5040.00/
+     TROE/   .6027  208.00  3922.00  10180.0 /
+H2/2.00/ H2O/6.00/ CH4/2.00/ CO/1.50/ CO2/2.00/ C2H6/3.00/
+CH2(S)+H2O<=>CH2+H2O                     3.000E+13     .000        .00
+CH2(S)+CH3<=>H+C2H4                      1.200E+13     .000    -570.00
+CH2(S)+CH4<=>2CH3                        1.600E+13     .000    -570.00
+CH2(S)+CO<=>CH2+CO                       9.000E+12     .000        .00
+CH2(S)+CO2<=>CH2+CO2                     7.000E+12     .000        .00
+CH2(S)+CO2<=>CO+CH2O                     1.400E+13     .000        .00
+CH2(S)+C2H6<=>CH3+C2H5                   4.000E+13     .000    -550.00
+CH3+O2<=>O+CH3O                          3.560E+13     .000   30480.00
+CH3+O2<=>OH+CH2O                         2.310E+12     .000   20315.00
+CH3+H2O2<=>HO2+CH4                       2.450E+04    2.470    5180.00
+2CH3(+M)<=>C2H6(+M)                      6.770E+16   -1.180     654.00
+     LOW  /  3.400E+41   -7.030   2762.00/
+     TROE/   .6190  73.20  1180.00  9999.00 /
+H2/2.00/ H2O/6.00/ CH4/2.00/ CO/1.50/ CO2/2.00/ C2H6/3.00/ AR/ .70/
+2CH3<=>H+C2H5                            6.840E+12     .100   10600.00
+CH3+HCO<=>CH4+CO                         2.648E+13     .000        .00
+CH3+CH2O<=>HCO+CH4                       3.320E+03    2.810    5860.00
+CH3+CH3OH<=>CH2OH+CH4                    3.000E+07    1.500    9940.00
+CH3+CH3OH<=>CH3O+CH4                     1.000E+07    1.500    9940.00
+CH3+C2H4<=>C2H3+CH4                      2.270E+05    2.000    9200.00
+CH3+C2H6<=>C2H5+CH4                      6.140E+06    1.740   10450.00
+HCO+H2O<=>H+CO+H2O                       1.500E+18   -1.000   17000.00
+HCO+M<=>H+CO+M                           1.870E+17   -1.000   17000.00
+H2/2.00/ H2O/ .00/ CH4/2.00/ CO/1.50/ CO2/2.00/ C2H6/3.00/
+HCO+O2<=>HO2+CO                          13.45E+12     .000     400.00
+CH2OH+O2<=>HO2+CH2O                      1.800E+13     .000     900.00
+CH3O+O2<=>HO2+CH2O                       4.280E-13    7.600   -3530.00
+C2H+O2<=>HCO+CO                          1.000E+13     .000    -755.00
+C2H+H2<=>H+C2H2                          5.680E+10    0.900    1993.00
+C2H3+O2<=>HCO+CH2O                       4.580E+16   -1.390    1015.00
+C2H4(+M)<=>H2+C2H2(+M)                   8.000E+12     .440   86770.00
+     LOW  /  1.580E+51   -9.300  97800.00/
+     TROE/   .7345  180.00  1035.00  5417.00 /
+H2/2.00/ H2O/6.00/ CH4/2.00/ CO/1.50/ CO2/2.00/ C2H6/3.00/ AR/ .70/
+C2H5+O2<=>HO2+C2H4                       8.400E+11     .000    3875.00
+HCCO+O2<=>OH+2CO                         3.200E+12     .000     854.00
+2HCCO<=>2CO+C2H2                         1.000E+13     .000        .00
+O+CH3=>H+H2+CO                           3.370E+13     .000        .00
+O+C2H4<=>H+CH2CHO                        6.700E+06    1.830     220.00
+O+C2H5<=>H+CH3CHO                        1.096E+14     .000        .00
+OH+HO2<=>O2+H2O                          0.500E+16     .000   17330.00
+  DUPLICATE
+OH+CH3=>H2+CH2O                          8.000E+09     .500   -1755.00
+CH+H2(+M)<=>CH3(+M)                      1.970E+12     .430    -370.00
+   LOW/ 4.820E+25  -2.80  590.0 /
+   TROE/ .578  122.0  2535.0  9365.0 /
+H2/2.00/ H2O/6.00/ CH4/2.00/ CO/1.50/ CO2/2.00/ C2H6/3.00/ AR/ .70/
+CH2+O2=>2H+CO2                           5.800E+12     .000    1500.00
+CH2+O2<=>O+CH2O                          2.400E+12     .000    1500.00
+CH2+CH2=>2H+C2H2                         2.000E+14     .000   10989.00
+CH2(S)+H2O=>H2+CH2O                      6.820E+10     .250    -935.00
+C2H3+O2<=>O+CH2CHO                       3.030E+11     .290      11.00
+C2H3+O2<=>HO2+C2H2                       1.337E+06    1.610    -384.00
+O+CH3CHO<=>OH+CH2CHO                     2.920E+12     .000    1808.00
+O+CH3CHO=>OH+CH3+CO                      2.920E+12     .000    1808.00
+O2+CH3CHO=>HO2+CH3+CO                    3.010E+13     .000   39150.00
+H+CH3CHO<=>CH2CHO+H2                     2.050E+09    1.160    2405.00
+H+CH3CHO=>CH3+H2+CO                      2.050E+09    1.160    2405.00
+OH+CH3CHO=>CH3+H2O+CO                    2.343E+10    0.730   -1113.00
+HO2+CH3CHO=>CH3+H2O2+CO                  3.010E+12     .000   11923.00
+CH3+CH3CHO=>CH3+CH4+CO                   2.720E+06    1.770    5920.00
+H+CH2CO(+M)<=>CH2CHO(+M)                 4.865E+11    0.422   -1755.00
+    LOW/ 1.012E+42  -7.63  3854.0/
+    TROE/ 0.465  201.0  1773.0  5333.0 /
+H2/2.00/ H2O/6.00/ CH4/2.00/ CO/1.50/ CO2/2.00/ C2H6/3.00/ AR/ .70/
+O+CH2CHO=>H+CH2+CO2                      1.500E+14     .000       .00
+O2+CH2CHO=>OH+CO+CH2O                    1.810E+10     .000       .00
+O2+CH2CHO=>OH+2HCO                       2.350E+10     .000       .00
+H+CH2CHO<=>CH3+HCO                       2.200E+13     .000       .00
+H+CH2CHO<=>CH2CO+H2                      1.100E+13     .000       .00
+OH+CH2CHO<=>H2O+CH2CO                    1.200E+13     .000       .00
+OH+CH2CHO<=>HCO+CH2OH                    3.010E+13     .000       .00
+CH3+C2H5(+M)<=>C3H8(+M)                  .9430E+13     .000       .00
+     LOW/ 2.710E+74  -16.82  13065.0 /
+     TROE/ .1527  291.0  2742.0  7748.0 /
+H2/2.00/ H2O/6.00/ CH4/2.00/ CO/1.50/ CO2/2.00/ C2H6/3.00/ AR/ .70/
+O+C3H8<=>OH+C3H7                         1.930E+05    2.680   3716.00
+H+C3H8<=>C3H7+H2                         1.320E+06    2.540   6756.00
+OH+C3H8<=>C3H7+H2O                       3.160E+07    1.800    934.00
+C3H7+H2O2<=>HO2+C3H8                     3.780E+02    2.720   1500.00
+CH3+C3H8<=>C3H7+CH4                      0.903E+00    3.650   7154.00
+CH3+C2H4(+M)<=>C3H7(+M)                  2.550E+06    1.600   5700.00
+      LOW/ 3.00E+63  -14.6  18170./
+      TROE/ .1894  277.0  8748.0  7891.0 /
+H2/2.00/ H2O/6.00/ CH4/2.00/ CO/1.50/ CO2/2.00/ C2H6/3.00/ AR/ .70/
+O+C3H7<=>C2H5+CH2O                       9.640E+13     .000       .00
+H+C3H7(+M)<=>C3H8(+M)                    3.613E+13     .000       .00
+      LOW/ 4.420E+61  -13.545  11357.0/
+      TROE/ .315  369.0  3285.0  6667.0 /
+H2/2.00/ H2O/6.00/ CH4/2.00/ CO/1.50/ CO2/2.00/ C2H6/3.00/ AR/ .70/
+H+C3H7<=>CH3+C2H5                        4.060E+06    2.190    890.00
+OH+C3H7<=>C2H5+CH2OH                     2.410E+13     .000       .00
+HO2+C3H7<=>O2+C3H8                       2.550E+10    0.255   -943.00
+HO2+C3H7=>OH+C2H5+CH2O                   2.410E+13     .000       .00
+CH3+C3H7<=>2C2H5                         1.927E+13   -0.320       .00
+END

tutorials/combustion/reactingFoam/RAS/DLR_A_LTS/chemkin/thermo30.dat

@@ -0,0 +1,218 @@
+THERMO ALL
+   250.000  1000.000  5000.000
+! GRI-Mech Version 3.0 Thermodynamics released 7/30/99
+! NASA Polynomial format for CHEMKIN-II
+! see README file for disclaimer
+O                 L 1/90O   1               G   200.000  3500.000  1000.000    1
+ 2.56942078E+00-8.59741137E-05 4.19484589E-08-1.00177799E-11 1.22833691E-15    2
+ 2.92175791E+04 4.78433864E+00 3.16826710E+00-3.27931884E-03 6.64306396E-06    3
+-6.12806624E-09 2.11265971E-12 2.91222592E+04 2.05193346E+00                   4
+O2                TPIS89O   2               G   200.000  3500.000  1000.000    1
+ 3.28253784E+00 1.48308754E-03-7.57966669E-07 2.09470555E-10-2.16717794E-14    2
+-1.08845772E+03 5.45323129E+00 3.78245636E+00-2.99673416E-03 9.84730201E-06    3
+-9.68129509E-09 3.24372837E-12-1.06394356E+03 3.65767573E+00                   4
+H                 L 7/88H   1               G   200.000  3500.000  1000.000    1
+ 2.50000001E+00-2.30842973E-11 1.61561948E-14-4.73515235E-18 4.98197357E-22    2
+ 2.54736599E+04-4.46682914E-01 2.50000000E+00 7.05332819E-13-1.99591964E-15    3
+ 2.30081632E-18-9.27732332E-22 2.54736599E+04-4.46682853E-01                   4
+H2                TPIS78H   2               G   200.000  3500.000  1000.000    1
+ 3.33727920E+00-4.94024731E-05 4.99456778E-07-1.79566394E-10 2.00255376E-14    2
+-9.50158922E+02-3.20502331E+00 2.34433112E+00 7.98052075E-03-1.94781510E-05    3
+ 2.01572094E-08-7.37611761E-12-9.17935173E+02 6.83010238E-01                   4
+OH                RUS 78O   1H   1          G   200.000  3500.000  1000.000    1
+ 3.09288767E+00 5.48429716E-04 1.26505228E-07-8.79461556E-11 1.17412376E-14    2
+ 3.85865700E+03 4.47669610E+00 3.99201543E+00-2.40131752E-03 4.61793841E-06    3
+-3.88113333E-09 1.36411470E-12 3.61508056E+03-1.03925458E-01                   4
+H2O               L 8/89H   2O   1          G   200.000  3500.000  1000.000    1
+ 3.03399249E+00 2.17691804E-03-1.64072518E-07-9.70419870E-11 1.68200992E-14    2
+-3.00042971E+04 4.96677010E+00 4.19864056E+00-2.03643410E-03 6.52040211E-06    3
+-5.48797062E-09 1.77197817E-12-3.02937267E+04-8.49032208E-01                   4
+HO2               L 5/89H   1O   2          G   200.000  3500.000  1000.000    1
+ 4.01721090E+00 2.23982013E-03-6.33658150E-07 1.14246370E-10-1.07908535E-14    2
+ 1.11856713E+02 3.78510215E+00 4.30179801E+00-4.74912051E-03 2.11582891E-05    3
+-2.42763894E-08 9.29225124E-12 2.94808040E+02 3.71666245E+00                   4
+H2O2              L 7/88H   2O   2          G   200.000  3500.000  1000.000    1
+ 4.16500285E+00 4.90831694E-03-1.90139225E-06 3.71185986E-10-2.87908305E-14    2
+-1.78617877E+04 2.91615662E+00 4.27611269E+00-5.42822417E-04 1.67335701E-05    3
+-2.15770813E-08 8.62454363E-12-1.77025821E+04 3.43505074E+00                   4
+C                 L11/88C   1               G   200.000  3500.000  1000.000    1
+ 2.49266888E+00 4.79889284E-05-7.24335020E-08 3.74291029E-11-4.87277893E-15    2
+ 8.54512953E+04 4.80150373E+00 2.55423955E+00-3.21537724E-04 7.33792245E-07    3
+-7.32234889E-10 2.66521446E-13 8.54438832E+04 4.53130848E+00                   4
+CH                TPIS79C   1H   1          G   200.000  3500.000  1000.000    1
+ 2.87846473E+00 9.70913681E-04 1.44445655E-07-1.30687849E-10 1.76079383E-14    2
+ 7.10124364E+04 5.48497999E+00 3.48981665E+00 3.23835541E-04-1.68899065E-06    3
+ 3.16217327E-09-1.40609067E-12 7.07972934E+04 2.08401108E+00                   4
+CH2               L S/93C   1H   2          G   200.000  3500.000  1000.000    1
+ 2.87410113E+00 3.65639292E-03-1.40894597E-06 2.60179549E-10-1.87727567E-14    2
+ 4.62636040E+04 6.17119324E+00 3.76267867E+00 9.68872143E-04 2.79489841E-06    3
+-3.85091153E-09 1.68741719E-12 4.60040401E+04 1.56253185E+00                   4
+CH2(S)            L S/93C   1H   2          G   200.000  3500.000  1000.000    1
+ 2.29203842E+00 4.65588637E-03-2.01191947E-06 4.17906000E-10-3.39716365E-14    2
+ 5.09259997E+04 8.62650169E+00 4.19860411E+00-2.36661419E-03 8.23296220E-06    3
+-6.68815981E-09 1.94314737E-12 5.04968163E+04-7.69118967E-01                   4
+CH3               L11/89C   1H   3          G   200.000  3500.000  1000.000    1
+ 2.28571772E+00 7.23990037E-03-2.98714348E-06 5.95684644E-10-4.67154394E-14    2
+ 1.67755843E+04 8.48007179E+00 3.67359040E+00 2.01095175E-03 5.73021856E-06    3
+-6.87117425E-09 2.54385734E-12 1.64449988E+04 1.60456433E+00                   4
+CH4               L 8/88C   1H   4          G   200.000  3500.000  1000.000    1
+ 7.48514950E-02 1.33909467E-02-5.73285809E-06 1.22292535E-09-1.01815230E-13    2
+-9.46834459E+03 1.84373180E+01 5.14987613E+00-1.36709788E-02 4.91800599E-05    3
+-4.84743026E-08 1.66693956E-11-1.02466476E+04-4.64130376E+00                   4
+CO                TPIS79C   1O   1          G   200.000  3500.000  1000.000    1
+ 2.71518561E+00 2.06252743E-03-9.98825771E-07 2.30053008E-10-2.03647716E-14    2
+-1.41518724E+04 7.81868772E+00 3.57953347E+00-6.10353680E-04 1.01681433E-06    3
+ 9.07005884E-10-9.04424499E-13-1.43440860E+04 3.50840928E+00                   4
+CO2               L 7/88C   1O   2          G   200.000  3500.000  1000.000    1
+ 3.85746029E+00 4.41437026E-03-2.21481404E-06 5.23490188E-10-4.72084164E-14    2
+-4.87591660E+04 2.27163806E+00 2.35677352E+00 8.98459677E-03-7.12356269E-06    3
+ 2.45919022E-09-1.43699548E-13-4.83719697E+04 9.90105222E+00                   4
+HCO               L12/89H   1C   1O   1     G   200.000  3500.000  1000.000    1
+ 2.77217438E+00 4.95695526E-03-2.48445613E-06 5.89161778E-10-5.33508711E-14    2
+ 4.01191815E+03 9.79834492E+00 4.22118584E+00-3.24392532E-03 1.37799446E-05    3
+-1.33144093E-08 4.33768865E-12 3.83956496E+03 3.39437243E+00                   4
+CH2O              L 8/88H   2C   1O   1     G   200.000  3500.000  1000.000    1
+ 1.76069008E+00 9.20000082E-03-4.42258813E-06 1.00641212E-09-8.83855640E-14    2
+-1.39958323E+04 1.36563230E+01 4.79372315E+00-9.90833369E-03 3.73220008E-05    3
+-3.79285261E-08 1.31772652E-11-1.43089567E+04 6.02812900E-01                   4
+CH2OH             GUNL93C   1H   3O   1     G   200.000  3500.000  1000.000    1
+ 3.69266569E+00 8.64576797E-03-3.75101120E-06 7.87234636E-10-6.48554201E-14    2
+-3.24250627E+03 5.81043215E+00 3.86388918E+00 5.59672304E-03 5.93271791E-06    3
+-1.04532012E-08 4.36967278E-12-3.19391367E+03 5.47302243E+00                   4
+CH3O              121686C   1H   3O   1     G   250.00   3000.00   1000.000    1
+ 0.03770799E+02 0.07871497E-01-0.02656384E-04 0.03944431E-08-0.02112616E-12    2
+ 0.12783252E+03 0.02929575E+02 0.02106204E+02 0.07216595E-01 0.05338472E-04    3
+-0.07377636E-07 0.02075610E-10 0.09786011E+04 0.13152177E+02                   4
+CH3OH             L 8/88C   1H   4O   1     G   200.000  3500.000  1000.000    1
+ 1.78970791E+00 1.40938292E-02-6.36500835E-06 1.38171085E-09-1.17060220E-13    2
+-2.53748747E+04 1.45023623E+01 5.71539582E+00-1.52309129E-02 6.52441155E-05    3
+-7.10806889E-08 2.61352698E-11-2.56427656E+04-1.50409823E+00                   4
+C2H               L 1/91C   2H   1          G   200.000  3500.000  1000.000    1
+ 3.16780652E+00 4.75221902E-03-1.83787077E-06 3.04190252E-10-1.77232770E-14    2
+ 6.71210650E+04 6.63589475E+00 2.88965733E+00 1.34099611E-02-2.84769501E-05    3
+ 2.94791045E-08-1.09331511E-11 6.68393932E+04 6.22296438E+00                   4
+C2H2              L 1/91C   2H   2          G   200.000  3500.000  1000.000    1
+ 4.14756964E+00 5.96166664E-03-2.37294852E-06 4.67412171E-10-3.61235213E-14    2
+ 2.59359992E+04-1.23028121E+00 8.08681094E-01 2.33615629E-02-3.55171815E-05    3
+ 2.80152437E-08-8.50072974E-12 2.64289807E+04 1.39397051E+01                   4
+C2H3              L 2/92C   2H   3          G   200.000  3500.000  1000.000    1
+ 3.01672400E+00 1.03302292E-02-4.68082349E-06 1.01763288E-09-8.62607041E-14    2
+ 3.46128739E+04 7.78732378E+00 3.21246645E+00 1.51479162E-03 2.59209412E-05    3
+-3.57657847E-08 1.47150873E-11 3.48598468E+04 8.51054025E+00                   4
+C2H4              L 1/91C   2H   4          G   200.000  3500.000  1000.000    1
+ 2.03611116E+00 1.46454151E-02-6.71077915E-06 1.47222923E-09-1.25706061E-13    2
+ 4.93988614E+03 1.03053693E+01 3.95920148E+00-7.57052247E-03 5.70990292E-05    3
+-6.91588753E-08 2.69884373E-11 5.08977593E+03 4.09733096E+00                   4
+C2H5              L12/92C   2H   5          G   200.000  3500.000  1000.000    1
+ 1.95465642E+00 1.73972722E-02-7.98206668E-06 1.75217689E-09-1.49641576E-13    2
+ 1.28575200E+04 1.34624343E+01 4.30646568E+00-4.18658892E-03 4.97142807E-05    3
+-5.99126606E-08 2.30509004E-11 1.28416265E+04 4.70720924E+00                   4
+C2H6              L 8/88C   2H   6          G   200.000  3500.000  1000.000    1
+ 1.07188150E+00 2.16852677E-02-1.00256067E-05 2.21412001E-09-1.90002890E-13    2
+-1.14263932E+04 1.51156107E+01 4.29142492E+00-5.50154270E-03 5.99438288E-05    3
+-7.08466285E-08 2.68685771E-11-1.15222055E+04 2.66682316E+00                   4
+CH2CO             L 5/90C   2H   2O   1     G   200.000  3500.000  1000.000    1
+ 4.51129732E+00 9.00359745E-03-4.16939635E-06 9.23345882E-10-7.94838201E-14    2
+-7.55105311E+03 6.32247205E-01 2.13583630E+00 1.81188721E-02-1.73947474E-05    3
+ 9.34397568E-09-2.01457615E-12-7.04291804E+03 1.22156480E+01                   4
+HCCO              SRIC91H   1C   2O   1     G   250.00   4000.00   1000.000    1
+ 0.56282058E+01 0.40853401E-02-0.15934547E-05 0.28626052E-09-0.19407832E-13    2
+ 0.19327215E+05-0.39302595E+01 0.22517214E+01 0.17655021E-01-0.23729101E-04    3
+ 0.17275759E-07-0.50664811E-11 0.20059449E+05 0.12490417E+02                   4
+HCCOH              SRI91C   2O   1H   2     G   250.000  5000.000  1000.000    1
+ 0.59238291E+01 0.67923600E-02-0.25658564E-05 0.44987841E-09-0.29940101E-13    2
+ 0.72646260E+04-0.76017742E+01 0.12423733E+01 0.31072201E-01-0.50866864E-04    3
+ 0.43137131E-07-0.14014594E-10 0.80316143E+04 0.13874319E+02                   4
+H2CN               41687H   2C   1N   1     G   250.00   4000.000  1000.000    1
+ 0.52097030E+01 0.29692911E-02-0.28555891E-06-0.16355500E-09 0.30432589E-13    2
+ 0.27677109E+05-0.44444780E+01 0.28516610E+01 0.56952331E-02 0.10711400E-05    3
+-0.16226120E-08-0.23511081E-12 0.28637820E+05 0.89927511E+01                   4
+HCN               GRI/98H   1C   1N   1     G   200.000  6000.000  1000.000    1
+ 0.38022392E+01 0.31464228E-02-0.10632185E-05 0.16619757E-09-0.97997570E-14    2
+ 0.14407292E+05 0.15754601E+01 0.22589886E+01 0.10051170E-01-0.13351763E-04    3
+ 0.10092349E-07-0.30089028E-11 0.14712633E+05 0.89164419E+01                   4
+HNO               And93 H   1N   1O   1     G   200.000  6000.000  1000.000    1
+ 0.29792509E+01 0.34944059E-02-0.78549778E-06 0.57479594E-10-0.19335916E-15    2
+ 0.11750582E+05 0.86063728E+01 0.45334916E+01-0.56696171E-02 0.18473207E-04    3
+-0.17137094E-07 0.55454573E-11 0.11548297E+05 0.17498417E+01                   4
+N                 L 6/88N   1               G   200.000  6000.000  1000.000    1
+ 0.24159429E+01 0.17489065E-03-0.11902369E-06 0.30226245E-10-0.20360982E-14    2
+ 0.56133773E+05 0.46496096E+01 0.25000000E+01 0.00000000E+00 0.00000000E+00    3
+ 0.00000000E+00 0.00000000E+00 0.56104637E+05 0.41939087E+01                   4
+NNH               T07/93N   2H   1          G   200.000  6000.000  1000.000    1
+ 0.37667544E+01 0.28915082E-02-0.10416620E-05 0.16842594E-09-0.10091896E-13    2
+ 0.28650697E+05 0.44705067E+01 0.43446927E+01-0.48497072E-02 0.20059459E-04    3
+-0.21726464E-07 0.79469539E-11 0.28791973E+05 0.29779410E+01                   4
+N2O               L 7/88N   2O   1          G   200.000  6000.000  1000.000    1
+ 0.48230729E+01 0.26270251E-02-0.95850874E-06 0.16000712E-09-0.97752303E-14    2
+ 0.80734048E+04-0.22017207E+01 0.22571502E+01 0.11304728E-01-0.13671319E-04    3
+ 0.96819806E-08-0.29307182E-11 0.87417744E+04 0.10757992E+02                   4
+NH                And94 N   1H   1          G   200.000  6000.000  1000.000    1
+ 0.27836928E+01 0.13298430E-02-0.42478047E-06 0.78348501E-10-0.55044470E-14    2
+ 0.42120848E+05 0.57407799E+01 0.34929085E+01 0.31179198E-03-0.14890484E-05    3
+ 0.24816442E-08-0.10356967E-11 0.41880629E+05 0.18483278E+01                   4
+NH2               And89 N   1H   2          G   200.000  6000.000  1000.000    1
+ 0.28347421E+01 0.32073082E-02-0.93390804E-06 0.13702953E-09-0.79206144E-14    2
+ 0.22171957E+05 0.65204163E+01 0.42040029E+01-0.21061385E-02 0.71068348E-05    3
+-0.56115197E-08 0.16440717E-11 0.21885910E+05-0.14184248E+00                   4
+NH3               J 6/77N   1H   3          G   200.000  6000.000  1000.000    1
+ 0.26344521E+01 0.56662560E-02-0.17278676E-05 0.23867161E-09-0.12578786E-13    2
+-0.65446958E+04 0.65662928E+01 0.42860274E+01-0.46605230E-02 0.21718513E-04    3
+-0.22808887E-07 0.82638046E-11-0.67417285E+04-0.62537277E+00                   4
+NO                RUS 78N   1O   1          G   200.000  6000.000  1000.000    1
+ 0.32606056E+01 0.11911043E-02-0.42917048E-06 0.69457669E-10-0.40336099E-14    2
+ 0.99209746E+04 0.63693027E+01 0.42184763E+01-0.46389760E-02 0.11041022E-04    3
+-0.93361354E-08 0.28035770E-11 0.98446230E+04 0.22808464E+01                   4
+NO2               L 7/88N   1O   2          G   200.000  6000.000  1000.000    1
+ 0.48847542E+01 0.21723956E-02-0.82806906E-06 0.15747510E-09-0.10510895E-13    2
+ 0.23164983E+04-0.11741695E+00 0.39440312E+01-0.15854290E-02 0.16657812E-04    3
+-0.20475426E-07 0.78350564E-11 0.28966179E+04 0.63119917E+01                   4
+HCNO              BDEA94H   1N   1C   1O   1G   250.000  5000.000  1382.000    1
+ 6.59860456E+00 3.02778626E-03-1.07704346E-06 1.71666528E-10-1.01439391E-14    2
+ 1.79661339E+04-1.03306599E+01 2.64727989E+00 1.27505342E-02-1.04794236E-05    3
+ 4.41432836E-09-7.57521466E-13 1.92990252E+04 1.07332972E+01                   4
+HOCN              BDEA94H   1N   1C   1O   1G   250.000  5000.000  1368.000    1
+ 5.89784885E+00 3.16789393E-03-1.11801064E-06 1.77243144E-10-1.04339177E-14    2
+-3.70653331E+03-6.18167825E+00 3.78604952E+00 6.88667922E-03-3.21487864E-06    3
+ 5.17195767E-10 1.19360788E-14-2.82698400E+03 5.63292162E+00                   4
+HNCO              BDEA94H   1N   1C   1O   1G   250.000  5000.000  1478.000    1
+ 6.22395134E+00 3.17864004E-03-1.09378755E-06 1.70735163E-10-9.95021955E-15    2
+-1.66599344E+04-8.38224741E+00 3.63096317E+00 7.30282357E-03-2.28050003E-06    3
+-6.61271298E-10 3.62235752E-13-1.55873636E+04 6.19457727E+00                   4
+NCO               EA 93 N   1C   1O   1     G   200.000  6000.000  1000.000    1
+ 0.51521845E+01 0.23051761E-02-0.88033153E-06 0.14789098E-09-0.90977996E-14    2
+ 0.14004123E+05-0.25442660E+01 0.28269308E+01 0.88051688E-02-0.83866134E-05    3
+ 0.48016964E-08-0.13313595E-11 0.14682477E+05 0.95504646E+01                   4
+CN                HBH92 C   1N   1          G   200.000  6000.000  1000.000    1
+ 0.37459805E+01 0.43450775E-04 0.29705984E-06-0.68651806E-10 0.44134173E-14    2
+ 0.51536188E+05 0.27867601E+01 0.36129351E+01-0.95551327E-03 0.21442977E-05    3
+-0.31516323E-09-0.46430356E-12 0.51708340E+05 0.39804995E+01                   4
+HCNN              SRI/94C   1N   2H   1     G   250.000  5000.000  1000.000    1
+ 0.58946362E+01 0.39895959E-02-0.15982380E-05 0.29249395E-09-0.20094686E-13    2
+ 0.53452941E+05-0.51030502E+01 0.25243194E+01 0.15960619E-01-0.18816354E-04    3
+ 0.12125540E-07-0.32357378E-11 0.54261984E+05 0.11675870E+02                   4
+N2                121286N   2               G   250.000  5000.000  1000.000    1
+ 0.02926640E+02 0.14879768E-02-0.05684760E-05 0.10097038E-09-0.06753351E-13    2
+-0.09227977E+04 0.05980528E+02 0.03298677E+02 0.14082404E-02-0.03963222E-04    3
+ 0.05641515E-07-0.02444854E-10-0.10208999E+04 0.03950372E+02                   4
+AR                120186AR  1               G   250.000  5000.000  1000.000    1
+ 0.02500000E+02 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00    2
+-0.07453750E+04 0.04366000E+02 0.02500000E+02 0.00000000E+00 0.00000000E+00    3
+ 0.00000000E+00 0.00000000E+00-0.07453750E+04 0.04366000E+02                   4
+C3H8              L 4/85C   3H   8          G   250.000  5000.000  1000.000    1
+ 0.75341368E+01 0.18872239E-01-0.62718491E-05 0.91475649E-09-0.47838069E-13    2
+-0.16467516E+05-0.17892349E+02 0.93355381E+00 0.26424579E-01 0.61059727E-05    3
+-0.21977499E-07 0.95149253E-11-0.13958520E+05 0.19201691E+02                   4
+C3H7              L 9/84C   3H   7          G   250.000  5000.000  1000.000    1
+ 0.77026987E+01 0.16044203E-01-0.52833220E-05 0.76298590E-09-0.39392284E-13    2
+ 0.82984336E+04-0.15480180E+02 0.10515518E+01 0.25991980E-01 0.23800540E-05    3
+-0.19609569E-07 0.93732470E-11 0.10631863E+05 0.21122559E+02                   4
+CH3CHO            L 8/88C   2H   4O   1     G   200.000  6000.000  1000.000    1
+ 0.54041108E+01 0.11723059E-01-0.42263137E-05 0.68372451E-09-0.40984863E-13    2
+-0.22593122E+05-0.34807917E+01 0.47294595E+01-0.31932858E-02 0.47534921E-04    3
+-0.57458611E-07 0.21931112E-10-0.21572878E+05 0.41030159E+01                   4
+CH2CHO            SAND86O   1H   3C   2     G   250.000  5000.000  1000.000    1
+ 0.05975670E+02 0.08130591E-01-0.02743624E-04 0.04070304E-08-0.02176017E-12    2
+ 0.04903218E+04-0.05045251E+02 0.03409062E+02 0.10738574E-01 0.01891492E-04    3
+-0.07158583E-07 0.02867385E-10 0.15214766E+04 0.09558290E+02                   4
+END

tutorials/combustion/reactingFoam/RAS/DLR_A_LTS/chemkin/transportProperties

@@ -0,0 +1,45 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  dev                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    location    "chemkin";
+    object      transportProperties;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+".*"
+{
+    transport
+    {
+        As 1.512e-06;
+        Ts 120.;
+    }
+}
+
+"H2"
+{
+    transport
+    {
+        As 6.362e-07;
+        Ts 72.;
+    }
+}
+
+"CO2"
+{
+    transport
+    {
+        As 1.572e-06;
+        Ts 240.;
+    }
+}
+
+// ************************************************************************* //

tutorials/combustion/reactingFoam/RAS/DLR_A_LTS/constant/chemistryProperties

@@ -0,0 +1,146 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  dev                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    location    "constant";
+    object      chemistryProperties;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+chemistryType
+{
+    chemistrySolver   ode;
+    chemistryThermo   psi;
+    TDAC              on;
+}
+
+importantSpecies
+{
+    CH4;
+    H2O;
+    O2;
+    CO2;
+}
+
+chemistry           on;
+
+initialChemicalTimeStep 1e-07;
+
+odeCoeffs
+{
+    solver          Rosenbrock34;  // Rosenbrock34, seulex or rodas23
+    absTol          1e-12;
+    relTol          1e-7;
+}
+
+reduction
+{
+    // Activate reduction
+    active  on;
+
+    // Switch logging of the reduction statistics and performance
+    log         on;
+
+    // Tolerance depends on the reduction method, see details for each method
+    tolerance   1e-4;
+
+    // Available methods: DRG, DAC, DRGEP, PFA, EFA
+    method DAC;
+
+    // Search initiating set (SIS) of species, needed for most methods
+    initialSet
+    {
+        CO;
+        CH4;
+        HO2;
+    }
+
+    // For DAC, option to automatically change the SIS switch from HO2 to H2O
+    // and CO to CO2, + disable fuel
+    automaticSIS    off;
+
+    // When automaticSIS, the method needs to know the fuel
+    fuelSpecies
+    {
+        CH4 1;
+    }
+}
+
+tabulation
+{
+    // Activate tabulation
+    active      on;
+
+    // Switch logging of the tabulation statistics and performance
+    log         on;
+
+    printProportion    off;
+
+    printNumRetrieve   off;
+
+    // Tolerance used for retrieve and grow
+    tolerance   1e-4;
+
+    // ISAT is the only method currently available
+    method    ISAT;
+
+    // Scale factors used in the definition of the ellipsoid of accuracy
+    scaleFactor
+    {
+        otherSpecies 1;
+        Temperature  1000;
+        Pressure     1e15;
+        deltaT       0.5;
+    }
+
+    // Maximum number of leafs stored in the binary tree
+    maxNLeafs  2000;
+
+    // Maximum life time of the leafs (in time steps) used in unsteady
+    // simulations to force renewal of the stored chemPoints and keep the tree
+    // small
+    chPMaxLifeTime 100;
+
+    // Maximum number of growth allowed on a chemPoint to avoid distorted
+    // chemPoints
+    maxGrowth  10;
+
+    // Number of time steps between analysis of the tree to remove old
+    // chemPoints or try to balance it
+    checkEntireTreeInterval  5;
+
+    // Parameters used to decide whether to balance or not if the tree's depth
+    // is larger than maxDepthFactor*log2(nLeafs) then balance the tree
+    maxDepthFactor   2;
+
+    // Try to balance the tree only if the size of the tree is greater
+    minBalanceThreshold 30;
+
+
+    // Activate the use of a MRU (most recently used) list
+    MRURetrieve false;
+
+    // Maximum size of the MRU list
+    maxMRUSize 0;
+
+    // Allow to grow points
+    growPoints  true;
+
+    // When mechanism reduction is used, new dimensions might be added
+    // maxNumNewDim set the maximum number of new dimensions added during a
+    // growth
+    maxNumNewDim 10;
+
+    variableTimeStep true;
+}
+
+
+// ************************************************************************* //

tutorials/combustion/reactingFoam/RAS/DLR_A_LTS/constant/chemistryProperties.test

@@ -0,0 +1,143 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  dev                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    location    "constant";
+    object      chemistryProperties;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+chemistryType
+{
+    chemistrySolver   ode;
+    chemistryThermo   psi;
+    TDAC              on;
+}
+
+chemistry           on;
+
+importantSpecies
+{
+    CH4;
+    H2O;
+    O2;
+    CO2;
+}
+
+initialChemicalTimeStep 1e-07;
+
+odeCoeffs
+{
+    solver          Rosenbrock34;  // Rosenbrock34, seulex or rodas23
+    absTol          1e-12;
+    relTol          1e-7;
+}
+
+reduction
+{
+    // Activate reduction
+    active  on;
+
+    // Switch logging of the reduction statistics and performance
+    log         on;
+
+    // Tolerance depends on the reduction method, see details for each method
+    tolerance   1e-4;
+
+    // Available methods: DRG, DAC, DRGEP, PFA, EFA
+    method DAC;
+
+    // Search initiating set (SIS) of species, needed for most methods
+    initialSet
+    {
+        CO;
+        CH4;
+        HO2;
+    }
+
+    // For DAC, option to automatically change the SIS switch from HO2 to H2O
+    // and CO to CO2, + disable fuel
+    automaticSIS    off;
+
+    // When automaticSIS, the method needs to know the fuel
+    fuelSpecies
+    {
+        CH4 1;
+    }
+}
+
+tabulation
+{
+    // Activate tabulation
+    active      on;
+
+    // Switch logging of the tabulation statistics and performance
+    log         on;
+
+    printProportion    off;
+
+    printNumRetrieve   off;
+
+    // Tolerance used for retrieve and grow
+    tolerance   0.003;
+
+    // ISAT is the only method currently available
+    method    ISAT;
+
+    // Scale factors used in the definition of the ellipsoid of accuracy
+    scaleFactor
+    {
+        otherSpecies 1;
+        Temperature  10000;
+        Pressure     1e15;
+        deltaT       1;
+    }
+
+    // Maximum number of leafs stored in the binary tree
+    maxNLeafs  5000;
+
+    // Maximum life time of the leafs (in time steps) used in unsteady
+    // simulations to force renewal of the stored chemPoints and keep the tree
+    // small
+    chPMaxLifeTime 1000;
+
+    // Maximum number of growth allowed on a chemPoint to avoid distorted
+    // chemPoints
+    maxGrowth  100;
+
+    // Number of time steps between analysis of the tree to remove old
+    // chemPoints or try to balance it
+    checkEntireTreeInterval  500;
+
+    // Parameters used to decide whether to balance or not if the tree's depth
+    // is larger than maxDepthFactor*log2(nLeafs) then balance the tree
+    maxDepthFactor   2;
+
+    // Try to balance the tree only if the size of the tree is greater
+    minBalanceThreshold 30;
+
+    // Activate the use of a MRU (most recently used) list
+    MRURetrieve false;
+
+    // Maximum size of the MRU list
+    maxMRUSize 0;
+
+    // Allow to grow points
+    growPoints  true;
+
+    // When mechanism reduction is used, new dimensions might be added
+    // maxNumNewDim set the maximum number of new dimensions added during a
+    // growth
+    maxNumNewDim 10;
+}
+
+
+// ************************************************************************* //

tutorials/combustion/reactingFoam/RAS/DLR_A_LTS/constant/combustionProperties

@@ -0,0 +1,27 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  dev                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    location    "constant";
+    object      combustionProperties;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+combustionModel  EDC<psiChemistryCombustion>;
+
+active  true;
+
+EDCCoeffs
+{
+    version v2005;
+}
+
+// ************************************************************************* //

tutorials/combustion/reactingFoam/RAS/DLR_A_LTS/constant/g

@@ -0,0 +1,21 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  dev                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       uniformDimensionedVectorField;
+    location    "constant";
+    object      g;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 1 -2 0 0 0 0];
+value           (0 0 -9.81);
+
+// ************************************************************************* //

tutorials/combustion/reactingFoam/RAS/DLR_A_LTS/constant/thermophysicalProperties

@@ -0,0 +1,36 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  dev                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    location    "constant";
+    object      thermophysicalProperties;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+thermoType
+{
+
+    type            hePsiThermo;
+    mixture         reactingMixture;
+    transport       sutherland;
+    thermo          janaf;
+    energy          sensibleEnthalpy;
+    equationOfState perfectGas;
+    specie          specie;
+}
+
+inertSpecie N2;
+
+chemistryReader foamChemistryReader;
+foamChemistryFile "$FOAM_CASE/constant/reactionsGRI";
+foamChemistryThermoFile "$FOAM_CASE/constant/thermo.compressibleGasGRI";
+
+// ************************************************************************* //

tutorials/combustion/reactingFoam/RAS/DLR_A_LTS/constant/turbulenceProperties

@@ -0,0 +1,33 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  dev                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    location    "constant";
+    object      turbulenceProperties;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+simulationType   RAS;
+
+RAS
+{
+    RASModel        kEpsilon;
+
+    kEpsilonCoeffs
+    {
+        Prt 0.85;
+    }
+
+    turbulence      on;
+    printCoeffs     on;
+}
+
+// ************************************************************************* //

tutorials/combustion/reactingFoam/RAS/DLR_A_LTS/system/blockMeshDict

@@ -0,0 +1,209 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  dev                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    object      blockMeshDict;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+convertToMeters 0.001;
+
+R1X 3.99619288633;
+R2X 4.49571699712;
+R3X 179.82867988473;
+R4X 189.81916210055;
+
+R1Y 0.17447754946;
+R2Y 0.19628724314;
+R3Y 7.85148972576;
+R4Y 8.28768359941;
+
+R1Ym -0.17447754946;
+R2Ym -0.19628724314;
+R3Ym -7.85148972576;
+R4Ym -8.28768359941;
+
+L  1000;
+Lm -20;
+
+vertices
+(
+    (0       0        0)    // 0
+    ($R1X    $R1Y     0)    // 1
+    ($R1X    $R1Y     $L)   // 2
+
+    (0       0        $L)   // 3
+    ($R1X    $R1Ym    0)    // 4
+    ($R1X    $R1Ym    $L)   // 5
+
+    (0       0        $Lm)  // 6
+    ($R1X    $R1Y     $Lm)  // 7
+    ($R1X    $R1Ym    $Lm)  // 8
+
+    ($R2X    $R2Y     0)    // 9
+    ($R2X    $R2Ym    0)    // 10
+    ($R2X    $R2Y     $L)   // 11
+    ($R2X    $R2Ym    $L)   // 12
+
+    ($R3X    $R3Y     0)    // 13
+    ($R3X    $R3Ym    0)    // 14
+    ($R3X    $R3Y     $L)   // 15
+    ($R3X    $R3Ym    $L)   // 16
+
+    ($R2X    $R2Y     $Lm)  // 17
+    ($R3X    $R3Y     $Lm)  // 18
+    ($R3X    $R3Ym    $Lm)  // 19
+    ($R2X    $R2Ym    $Lm)  // 20
+
+    ($R4X    $R4Y     0)    // 21
+    ($R4X    $R4Ym    0)    // 22
+    ($R4X    $R4Y     $L)   // 23
+    ($R4X    $R4Ym    $L)   // 24
+);
+
+nFuel          4;
+nBurner        1;
+nCoflow        30;
+nExternal      2;
+nLength        90;
+nLengthReverse 4;
+
+gradingFuel          1;
+gradingCoflow        6;
+gradingLength        12.;
+gradingLengthInverse 0.5;
+
+blocks
+(
+    // Fuel
+    hex (6 8 7 6 0 4 1 0) ($nFuel 1 $nLengthReverse)
+    simpleGrading ($gradingFuel  1 $gradingLengthInverse )
+
+    hex (0 4 1 0 3 5 2 3) ($nFuel  1 $nLength)
+    simpleGrading ($gradingFuel 1 $gradingLength)
+
+    // Wall
+    hex (4 10 9 1 5 12 11 2) ($nBurner 1 $nLength)
+    simpleGrading (1 1 $gradingLength)
+
+    // Coflow
+    hex (20 19 18 17 10 14 13 9) ($nCoflow 1 $nLengthReverse)
+    simpleGrading ($gradingCoflow  1 $gradingLengthInverse)
+    hex (10 14 13 9 12 16 15 11) ($nCoflow 1 $nLength)
+    simpleGrading ($gradingCoflow  1 $gradingLength)
+
+    // External wall
+    hex (14 22 21 13 16 24 23 15) ($nExternal 1 $nLength)
+    simpleGrading (1  1 $gradingLength)
+);
+
+boundary
+(
+    inletfuel
+    {
+        type patch;
+        faces
+        (
+            (6 8 7 6)
+        );
+    }
+
+    inletair
+    {
+        type patch;
+        faces
+        (
+            (19 20 17 18)
+        );
+    }
+
+    outlet
+    {
+        type patch;
+        faces
+        (
+            (5 3 3 2)
+            (12 11 2 5)
+            (16 15 11 12)
+            (15 16 24 23)
+        );
+    }
+
+    axis
+    {
+        type empty;
+        faces
+        (
+            (3 0 0 3)
+            (0 6 6 0)
+        );
+    }
+
+    leftside
+    {
+        type wall;
+        faces
+        (
+            (14 13 21 22)
+            (19 18 13 14)
+            (23 24 22 21)
+        );
+    }
+
+    burnerwall
+    {
+        type wall;
+        faces
+        (
+            (8 7 1 4)
+            (10 9 17 20)
+        );
+    }
+
+    burnertip
+    {
+        type wall;
+        faces
+        (
+            (4 1 9 10)
+        );
+    }
+
+    front
+    {
+        type wedge;
+        faces
+        (
+            (1 0 3 2)
+            (7 6 0 1)
+            (9 1 2 11)
+            (13 9 11 15)
+            (18 17 9 13)
+            (15 23 21 13)
+        );
+    }
+
+    back
+    {
+        type wedge;
+        faces
+        (
+            (5 3 0 4)
+            (4 0 6 8)
+            (12 5 4 10)
+            (16 12 10 14)
+            (19 14 10 20)
+            (16 14 22 24)
+        );
+    }
+);
+
+// ************************************************************************* //

tutorials/combustion/reactingFoam/RAS/DLR_A_LTS/system/controlDict

@@ -0,0 +1,48 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  dev                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    location    "system";
+    object      controlDict;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+application     reactingFoam;
+
+startFrom       latestTime;
+
+startTime       0;
+
+stopAt          endTime;
+
+endTime         70000;
+
+deltaT          1;
+
+writeControl    runTime;
+
+writeInterval   1000;
+
+purgeWrite      0;
+
+writeFormat     binary;
+
+writePrecision  10;
+
+writeCompression no;
+
+timeFormat      general;
+
+timePrecision   6;
+
+runTimeModifiable true;
+
+// ************************************************************************* //

tutorials/combustion/reactingFoam/RAS/DLR_A_LTS/system/decomposeParDict

@@ -0,0 +1,45 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  dev                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    location    "system";
+    object      decomposeParDict;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+numberOfSubdomains 6;
+
+method          simple;
+
+simpleCoeffs
+{
+    n               (1 1 6);
+    delta           0.001;
+}
+
+hierarchicalCoeffs
+{
+    n               ( 1 1 1 );
+    delta           0.001;
+    order           xyz;
+}
+
+manualCoeffs
+{
+    dataFile        "";
+}
+
+distributed     no;
+
+roots           ( );
+
+
+// ************************************************************************* //

tutorials/combustion/reactingFoam/RAS/DLR_A_LTS/system/fvSchemes

@@ -0,0 +1,59 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  dev                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    location    "system";
+    object      fvSchemes;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+ddtSchemes
+{
+    default         localEuler;
+}
+
+gradSchemes
+{
+    default         Gauss linear;
+}
+
+divSchemes
+{
+    default             none;
+
+    div(phi,U)          Gauss limitedLinearV 1;
+    div(phi,Yi)         Gauss limitedLinear01 1;
+    div(phi,h)          Gauss limitedLinear 1;
+    div(phi,K)          Gauss limitedLinear 1;
+    div(phid,p)         Gauss limitedLinear 1;
+    div(phi,epsilon)    Gauss limitedLinear 1;
+    div(phi,Yi_h)       Gauss limitedLinear01 1;
+    div(phi,k)          Gauss limitedLinear 1;
+    div(((rho*nuEff)*dev2(T(grad(U)))))     Gauss linear;
+}
+
+laplacianSchemes
+{
+    default         Gauss linear orthogonal;
+}
+
+interpolationSchemes
+{
+    default         linear;
+}
+
+snGradSchemes
+{
+    default         orthogonal;
+}
+
+
+// ************************************************************************* //

tutorials/combustion/reactingFoam/RAS/DLR_A_LTS/system/fvSolution

@@ -0,0 +1,89 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  dev                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    location    "system";
+    object      fvSolution;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+solvers
+{
+    "rho.*"
+    {
+        solver          diagonal;
+    }
+
+    p
+    {
+        solver           PCG;
+        preconditioner   DIC;
+        tolerance        1e-6;
+        relTol           0.01;
+    }
+
+    pFinal
+    {
+        $p;
+        relTol           0;
+    }
+
+    "(U|h|k|epsilon)"
+    {
+        solver          PBiCGStab;
+        preconditioner  DILU;
+        tolerance       1e-6;
+        relTol          0.1;
+    }
+
+    "(U|h|k|epsilon)Final"
+    {
+        $U;
+    }
+
+    Yi
+    {
+        solver          PBiCG;
+        preconditioner  DILU;
+        tolerance       1e-8;
+        relTol          0.1;
+    }
+}
+
+PIMPLE
+{
+    momentumPredictor   yes;
+    nOuterCorrectors    1;
+    nCorrectors         2;
+    nNonOrthogonalCorrectors 0;
+
+    maxDeltaT           1e-4;
+    maxCo               0.25;
+    alphaTemp           0.05;
+    alphaY              0.05;
+    Yref
+    {
+        O2                  0.1;
+        CH4                 0.1;
+    }
+    rDeltaTSmoothingCoeff   0.05;
+    rDeltaTDampingCoeff     1;
+}
+
+relaxationFactors
+{
+    equations
+    {
+        ".*"    1;
+    }
+}
+
+// ************************************************************************* //

tutorials/combustion/reactingFoam/RAS/DLR_A_LTS/system/setFieldsDict

@@ -0,0 +1,36 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  dev                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    location    "system";
+    object      setFieldsDict;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+defaultFieldValues
+(
+    volScalarFieldValue T 292
+);
+
+regions
+(
+    boxToCell
+    {
+        box (0.002 -0.01 0.005) (0.02 0.01 0.055);
+        fieldValues
+        (
+            volScalarFieldValue T 2200
+        );
+    }
+);
+
+
+// ************************************************************************* //

tutorials/combustion/reactingFoam/RAS/SandiaD_LTS/0.orig/CH4

@@ -0,0 +1,69 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  dev                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    location    "0";
+    object      CH4;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 0 0 0 0 0 0];
+
+internalField   uniform 0;
+
+boundaryField
+{
+    inletCH4
+    {
+        type            fixedValue;
+        value           uniform 0.1561;
+    }
+
+    wallOutside
+    {
+        type            zeroGradient;
+    }
+
+    wallTube
+    {
+        type            zeroGradient;
+    }
+
+    inletPilot
+    {
+        type            fixedValue;
+        value           uniform 0;
+    }
+
+    inletAir
+    {
+        type            fixedValue;
+        value           uniform 0;
+    }
+
+    outlet
+    {
+        type            zeroGradient;
+    }
+
+    frontAndBack_pos
+    {
+        type            wedge;
+    }
+
+    frontAndBack_neg
+    {
+        type            wedge;
+    }
+}
+
+
+// ************************************************************************* //

tutorials/combustion/reactingFoam/RAS/SandiaD_LTS/0.orig/CO

@@ -0,0 +1,69 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  dev                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    location    "0";
+    object      CO;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 0 0 0 0 0 0];
+
+internalField   uniform 0;
+
+boundaryField
+{
+    wallTube
+    {
+        type            zeroGradient;
+    }
+
+    outlet
+    {
+        type            zeroGradient;
+    }
+
+    inletPilot
+    {
+        type            fixedValue;
+        value           uniform 4.07e-3;
+    }
+
+    inletAir
+    {
+        type            fixedValue;
+        value           uniform 0;
+    }
+
+    wallOutside
+    {
+        type            zeroGradient;
+    }
+
+    inletCH4
+    {
+        type            fixedValue;
+        value           uniform 0;
+    }
+
+    frontAndBack_pos
+    {
+        type            wedge;
+    }
+
+    frontAndBack_neg
+    {
+        type            wedge;
+    }
+}
+
+
+// ************************************************************************* //

tutorials/combustion/reactingFoam/RAS/SandiaD_LTS/0.orig/CO2

@@ -0,0 +1,69 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  dev                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    location    "0";
+    object      CO2;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 0 0 0 0 0 0];
+
+internalField   uniform 0;
+
+boundaryField
+{
+    wallTube
+    {
+        type            zeroGradient;
+    }
+
+    outlet
+    {
+        type            zeroGradient;
+    }
+
+    inletPilot
+    {
+        type            fixedValue;
+        value           uniform 0.1098;
+    }
+
+    inletAir
+    {
+        type            fixedValue;
+        value           uniform 0;
+    }
+
+    wallOutside
+    {
+        type            zeroGradient;
+    }
+
+    inletCH4
+    {
+        type            fixedValue;
+        value           uniform 0;
+    }
+
+    frontAndBack_pos
+    {
+        type            wedge;
+    }
+
+    frontAndBack_neg
+    {
+        type            wedge;
+    }
+}
+
+
+// ************************************************************************* //

tutorials/combustion/reactingFoam/RAS/SandiaD_LTS/0.orig/H

@@ -0,0 +1,69 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  dev                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    location    "0";
+    object      H;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 0 0 0 0 0 0];
+
+internalField   uniform 0;
+
+boundaryField
+{
+    wallTube
+    {
+        type            zeroGradient;
+    }
+
+    outlet
+    {
+        type            zeroGradient;
+    }
+
+    inletPilot
+    {
+        type            fixedValue;
+        value           uniform 2.48e-5;
+    }
+
+    inletAir
+    {
+        type            fixedValue;
+        value           uniform 0;
+    }
+
+    wallOutside
+    {
+        type            zeroGradient;
+    }
+
+    inletCH4
+    {
+        type            fixedValue;
+        value           uniform 0;
+    }
+
+    frontAndBack_pos
+    {
+        type            wedge;
+    }
+
+    frontAndBack_neg
+    {
+        type            wedge;
+    }
+}
+
+
+// ************************************************************************* //

tutorials/combustion/reactingFoam/RAS/SandiaD_LTS/0.orig/H2

@@ -0,0 +1,69 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  dev                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    location    "0";
+    object      H2;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 0 0 0 0 0 0];
+
+internalField   uniform 0;
+
+boundaryField
+{
+    wallTube
+    {
+        type            zeroGradient;
+    }
+
+    outlet
+    {
+        type            zeroGradient;
+    }
+
+    inletPilot
+    {
+        type            fixedValue;
+        value           uniform 1.29e-4;
+    }
+
+    inletAir
+    {
+        type            fixedValue;
+        value           uniform 0;
+    }
+
+    wallOutside
+    {
+        type            zeroGradient;
+    }
+
+    inletCH4
+    {
+        type            fixedValue;
+        value           uniform 0;
+    }
+
+    frontAndBack_pos
+    {
+        type            wedge;
+    }
+
+    frontAndBack_neg
+    {
+        type            wedge;
+    }
+}
+
+
+// ************************************************************************* //

tutorials/combustion/reactingFoam/RAS/SandiaD_LTS/0.orig/H2O

@@ -0,0 +1,69 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  dev                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    location    "0";
+    object      H2O;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 0 0 0 0 0 0];
+
+internalField   uniform 0;
+
+boundaryField
+{
+    wallTube
+    {
+        type            zeroGradient;
+    }
+
+    outlet
+    {
+        type            zeroGradient;
+    }
+
+    inletPilot
+    {
+        type            fixedValue;
+        value           uniform 0.0942;
+    }
+
+    inletAir
+    {
+        type            fixedValue;
+        value           uniform 0;
+    }
+
+    wallOutside
+    {
+        type            zeroGradient;
+    }
+
+    inletCH4
+    {
+        type            fixedValue;
+        value           uniform 0;
+    }
+
+    frontAndBack_pos
+    {
+        type            wedge;
+    }
+
+    frontAndBack_neg
+    {
+        type            wedge;
+    }
+}
+
+
+// ************************************************************************* //

tutorials/combustion/reactingFoam/RAS/SandiaD_LTS/0.orig/N2

@@ -0,0 +1,69 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  dev                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    location    "0";
+    object      N2;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 0 0 0 0 0 0];
+
+internalField   uniform 0.77;
+
+boundaryField
+{
+    inletCH4
+    {
+        type            fixedValue;
+        value           uniform 0.6473;
+    }
+
+    wallOutside
+    {
+        type            zeroGradient;
+    }
+
+    wallTube
+    {
+        type            zeroGradient;
+    }
+
+    inletPilot
+    {
+        type            fixedValue;
+        value           uniform 0.7342;
+    }
+
+    inletAir
+    {
+        type            fixedValue;
+        value           uniform 0.77;
+    }
+
+    outlet
+    {
+        type            zeroGradient;
+    }
+
+    frontAndBack_pos
+    {
+        type            wedge;
+    }
+
+    frontAndBack_neg
+    {
+        type            wedge;
+    }
+}
+
+
+// ************************************************************************* //

tutorials/combustion/reactingFoam/RAS/SandiaD_LTS/0.orig/O

@@ -0,0 +1,69 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  dev                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    location    "0";
+    object      O;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 0 0 0 0 0 0];
+
+internalField   uniform 0;
+
+boundaryField
+{
+    wallTube
+    {
+        type            zeroGradient;
+    }
+
+    outlet
+    {
+        type            zeroGradient;
+    }
+
+    inletPilot
+    {
+        type            fixedValue;
+        value           uniform 7.47e-4;
+    }
+
+    inletAir
+    {
+        type            fixedValue;
+        value           uniform 0;
+    }
+
+    wallOutside
+    {
+        type            zeroGradient;
+    }
+
+    inletCH4
+    {
+        type            fixedValue;
+        value           uniform 0;
+    }
+
+    frontAndBack_pos
+    {
+        type            wedge;
+    }
+
+    frontAndBack_neg
+    {
+        type            wedge;
+    }
+}
+
+
+// ************************************************************************* //

tutorials/combustion/reactingFoam/RAS/SandiaD_LTS/0.orig/O2

@@ -0,0 +1,69 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  dev                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    location    "0";
+    object      O2;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 0 0 0 0 0 0];
+
+internalField   uniform 0.23;
+
+boundaryField
+{
+    inletCH4
+    {
+        type            fixedValue;
+        value           uniform 0.1966;
+    }
+
+    wallOutside
+    {
+        type            zeroGradient;
+    }
+
+    wallTube
+    {
+        type            zeroGradient;
+    }
+
+    inletPilot
+    {
+        type            fixedValue;
+        value           uniform 0.054;
+    }
+
+    inletAir
+    {
+        type            fixedValue;
+        value           uniform 0.23;
+    }
+
+    outlet
+    {
+        type            zeroGradient;
+    }
+
+    frontAndBack_pos
+    {
+        type            wedge;
+    }
+
+    frontAndBack_neg
+    {
+        type            wedge;
+    }
+}
+
+
+// ************************************************************************* //

tutorials/combustion/reactingFoam/RAS/SandiaD_LTS/0.orig/OH

@@ -0,0 +1,69 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  dev                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    location    "0";
+    object      OH;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 0 0 0 0 0 0];
+
+internalField   uniform 0;
+
+boundaryField
+{
+    wallTube
+    {
+        type            zeroGradient;
+    }
+
+    outlet
+    {
+        type            zeroGradient;
+    }
+
+    inletPilot
+    {
+        type            fixedValue;
+        value           uniform 0.0028;
+    }
+
+    inletAir
+    {
+        type            fixedValue;
+        value           uniform 0;
+    }
+
+    wallOutside
+    {
+        type            zeroGradient;
+    }
+
+    inletCH4
+    {
+        type            fixedValue;
+        value           uniform 0;
+    }
+
+    frontAndBack_pos
+    {
+        type            wedge;
+    }
+
+    frontAndBack_neg
+    {
+        type            wedge;
+    }
+}
+
+
+// ************************************************************************* //

tutorials/combustion/reactingFoam/RAS/SandiaD_LTS/0.orig/T

@@ -0,0 +1,69 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  dev                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    location    "0";
+    object      T;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 0 0 1 0 0 0];
+
+internalField   uniform 300;
+
+boundaryField
+{
+    inletCH4
+    {
+        type            fixedValue;
+        value           uniform 294;
+    }
+
+    wallOutside
+    {
+        type            zeroGradient;
+    }
+
+    wallTube
+    {
+        type            zeroGradient;
+    }
+
+    inletPilot
+    {
+        type            fixedValue;
+        value           uniform 1880;
+    }
+
+    inletAir
+    {
+        type            fixedValue;
+        value           uniform 291;
+    }
+
+    outlet
+    {
+        type            zeroGradient;
+    }
+
+    frontAndBack_pos
+    {
+        type            wedge;
+    }
+
+    frontAndBack_neg
+    {
+        type            wedge;
+    }
+}
+
+
+// ************************************************************************* //

tutorials/combustion/reactingFoam/RAS/SandiaD_LTS/0.orig/U

@@ -0,0 +1,70 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  dev                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volVectorField;
+    location    "0";
+    object      U;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 1 -1 0 0 0 0];
+
+internalField   uniform (0 0 0.9);
+
+boundaryField
+{
+    wallTube
+    {
+        type            noSlip;
+    }
+
+    outlet
+    {
+        type            pressureInletOutletVelocity;
+        value           $internalField;
+    }
+
+    inletPilot
+    {
+        type            fixedValue;
+        value           uniform (0 0 11.4);
+    }
+
+    inletAir
+    {
+        type            fixedValue;
+        value           uniform (0 0 0.9);
+    }
+
+    wallOutside
+    {
+        type            zeroGradient;
+    }
+
+    inletCH4
+    {
+        type            fixedValue;
+        value           uniform (0 0 49.6);
+    }
+
+    frontAndBack_pos
+    {
+        type            wedge;
+    }
+
+    frontAndBack_neg
+    {
+        type            wedge;
+    }
+}
+
+
+// ************************************************************************* //

tutorials/combustion/reactingFoam/RAS/SandiaD_LTS/0.orig/Ydefault

@@ -0,0 +1,69 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  dev                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    location    "0";
+    object      Ydefault;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 0 0 0 0 0 0];
+
+internalField   uniform 0;
+
+boundaryField
+{
+    wallTube
+    {
+        type            zeroGradient;
+    }
+
+    outlet
+    {
+        type            zeroGradient;
+    }
+
+    inletPilot
+    {
+        type            fixedValue;
+        value           uniform 0;
+    }
+
+    inletAir
+    {
+        type            fixedValue;
+        value           uniform 0;
+    }
+
+    wallOutside
+    {
+        type            zeroGradient;
+    }
+
+    inletCH4
+    {
+        type            fixedValue;
+        value           uniform 0;
+    }
+
+    frontAndBack_pos
+    {
+        type            wedge;
+    }
+
+    frontAndBack_neg
+    {
+        type            wedge;
+    }
+}
+
+
+// ************************************************************************* //

tutorials/combustion/reactingFoam/RAS/SandiaD_LTS/0.orig/alphat

@@ -0,0 +1,74 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  dev                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    location    "0";
+    object      alphat;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [1 -1 -1 0 0 0 0];
+
+internalField   uniform 0;
+
+boundaryField
+{
+    wallTube
+    {
+        type            compressible::alphatWallFunction;
+        Prt             0.85;
+        value           $internalField;
+    }
+
+    outlet
+    {
+        type            calculated;
+        value           $internalField;
+    }
+
+    inletPilot
+    {
+        type            calculated;
+        value           $internalField;
+    }
+
+    inletAir
+    {
+        type            calculated;
+        value           $internalField;
+    }
+
+    wallOutside
+    {
+        type            compressible::alphatWallFunction;
+        Prt             0.85;
+        value           $internalField;
+    }
+
+    inletCH4
+    {
+        type            calculated;
+        value           $internalField;
+    }
+
+    frontAndBack_pos
+    {
+        type            wedge;
+    }
+
+    frontAndBack_neg
+    {
+        type            wedge;
+    }
+}
+
+
+// ************************************************************************* //

tutorials/combustion/reactingFoam/RAS/SandiaD_LTS/0.orig/epsilon

@@ -0,0 +1,84 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  dev                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    location    "0";
+    object      epsilon;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 2 -3 0 0 0 0];
+
+internalField   uniform 30000;
+
+boundaryField
+{
+    wallTube
+    {
+        type            epsilonWallFunction;
+        Cmu             0.09;
+        kappa           0.41;
+        E               9.8;
+        value           uniform 30000;
+    }
+
+    outlet
+    {
+        type            zeroGradient;
+    }
+
+    inletPilot
+    {
+        type            turbulentMixingLengthDissipationRateInlet;
+        mixingLength    0.000735;
+        phi             phi;
+        k               k;
+        value           uniform 1;
+    }
+
+    inletAir
+    {
+        type            turbulentMixingLengthDissipationRateInlet;
+        mixingLength    0.019677;
+        value           uniform 1;
+    }
+
+    wallOutside
+    {
+        type            epsilonWallFunction;
+        Cmu             0.09;
+        kappa           0.41;
+        E               9.8;
+        value           uniform 30000;
+    }
+
+    inletCH4
+    {
+        type            turbulentMixingLengthDissipationRateInlet;
+        mixingLength    0.000504;
+        phi             phi;
+        k               k;
+        value           uniform 1;
+    }
+
+    frontAndBack_pos
+    {
+        type            wedge;
+    }
+
+    frontAndBack_neg
+    {
+        type            wedge;
+    }
+}
+
+
+// ************************************************************************* //

tutorials/combustion/reactingFoam/RAS/SandiaD_LTS/0.orig/k

@@ -0,0 +1,74 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  dev                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    location    "0";
+    object      k;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 2 -2 0 0 0 0];
+
+internalField   uniform 30;
+
+boundaryField
+{
+    wallTube
+    {
+        type            kqRWallFunction;
+        value           uniform 30;
+    }
+
+    outlet
+    {
+        type            zeroGradient;
+    }
+
+    inletPilot
+    {
+        type            turbulentIntensityKineticEnergyInlet;
+        intensity       0.0628;
+        value           uniform 1;
+    }
+
+    inletAir
+    {
+        type            turbulentIntensityKineticEnergyInlet;
+        intensity       0.0471;
+        value           uniform 1;
+    }
+
+    wallOutside
+    {
+        type            kqRWallFunction;
+        value           uniform 30;
+    }
+
+    inletCH4
+    {
+        type            turbulentIntensityKineticEnergyInlet;
+        intensity       0.0458;
+        value           uniform 1;
+    }
+
+    frontAndBack_pos
+    {
+        type            wedge;
+    }
+
+    frontAndBack_neg
+    {
+        type            wedge;
+    }
+}
+
+
+// ************************************************************************* //

tutorials/combustion/reactingFoam/RAS/SandiaD_LTS/0.orig/nut

@@ -0,0 +1,78 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  dev                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    location    "0";
+    object      nut;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 2 -1 0 0 0 0];
+
+internalField   uniform 0;
+
+boundaryField
+{
+    wallTube
+    {
+        type            nutkWallFunction;
+        Cmu             0.09;
+        kappa           0.41;
+        E               9.8;
+        value           $internalField;
+    }
+
+    outlet
+    {
+        type            calculated;
+        value           $internalField;
+    }
+
+    inletPilot
+    {
+        type            calculated;
+        value           $internalField;
+    }
+
+    inletAir
+    {
+        type            calculated;
+        value           $internalField;
+    }
+
+    wallOutside
+    {
+        type            nutkWallFunction;
+        Cmu             0.09;
+        kappa           0.41;
+        E               9.8;
+        value           $internalField;
+    }
+
+    inletCH4
+    {
+        type            calculated;
+        value           $internalField;
+    }
+
+    frontAndBack_pos
+    {
+        type            wedge;
+    }
+
+    frontAndBack_neg
+    {
+        type            wedge;
+    }
+}
+
+
+// ************************************************************************* //

tutorials/combustion/reactingFoam/RAS/SandiaD_LTS/0.orig/p

@@ -0,0 +1,67 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  dev                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    location    "0";
+    object      p;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [1 -1 -2 0 0 0 0];
+
+internalField   uniform 100000;
+
+boundaryField
+{
+    wallTube
+    {
+        type            zeroGradient;
+    }
+
+    outlet
+    {
+        type            totalPressure;
+        p0              $internalField;
+    }
+
+    inletPilot
+    {
+        type            zeroGradient;
+    }
+
+    inletAir
+    {
+        type            zeroGradient;
+    }
+
+    wallOutside
+    {
+        type            zeroGradient;
+    }
+
+    inletCH4
+    {
+        type            zeroGradient;
+    }
+
+    frontAndBack_pos
+    {
+        type            wedge;
+    }
+
+    frontAndBack_neg
+    {
+        type            wedge;
+    }
+}
+
+
+// ************************************************************************* //

tutorials/combustion/reactingFoam/RAS/SandiaD_LTS/Allclean

@@ -0,0 +1,11 @@
+#!/bin/sh
+cd ${0%/*} || exit 1    # Run from this directory
+
+# Source tutorial clean functions
+. $WM_PROJECT_DIR/bin/tools/CleanFunctions
+
+cleanCase
+
+rm -rf 0
+
+#------------------------------------------------------------------------------

tutorials/combustion/reactingFoam/RAS/SandiaD_LTS/Allrun

@@ -0,0 +1,36 @@
+#!/bin/sh
+cd ${0%/*} || exit 1    # Run from this directory
+
+# Source tutorial run functions
+. $WM_PROJECT_DIR/bin/tools/RunFunctions
+
+# Set application name
+application=`getApplication`
+
+rm -f 0
+cp -r 0.orig 0
+
+runApplication chemkinToFoam \
+               chemkin/grimech30.dat chemkin/thermo30.dat chemkin/transportProperties \
+               constant/reactionsGRI constant/thermo.compressibleGasGRI
+
+runApplication blockMesh
+runApplication setFields
+
+
+# Run the application without chemistry until 1500 to let the flow field develop
+foamDictionary -entry "writeInterval" -set "1500" system/controlDict
+foamDictionary -entry "endTime" -set "1500" system/controlDict
+foamDictionary -entry "chemistry" -set "off" constant/chemistryProperties
+
+runApplication $application
+
+
+# Run with chemistry until flame reach its full size
+foamDictionary -entry "writeInterval" -set "100" system/controlDict
+foamDictionary -entry "endTime" -set "5000" system/controlDict
+foamDictionary -entry "chemistry" -set "on" constant/chemistryProperties
+
+runApplication -o $application
+
+#------------------------------------------------------------------------------