Merge branch 'master' of ssh://git.code.sf.net/p/foam-extend/foam-extend-3.2
[foam-extend-3.2.git] / src / thermophysicalModels / chemistryModel / chemistrySolver / odeSolver / odeChemistrySolver.C
blob71e914038f440c974dc5356061ec7fafbdc975b7
1 /*---------------------------------------------------------------------------*\
2   =========                 |
3   \\      /  F ield         | foam-extend: Open Source CFD
4    \\    /   O peration     | Version:     3.2
5     \\  /    A nd           | Web:         http://www.foam-extend.org
6      \\/     M anipulation  | For copyright notice see file Copyright
7 -------------------------------------------------------------------------------
8 License
9     This file is part of foam-extend.
11     foam-extend is free software: you can redistribute it and/or modify it
12     under the terms of the GNU General Public License as published by the
13     Free Software Foundation, either version 3 of the License, or (at your
14     option) any later version.
16     foam-extend is distributed in the hope that it will be useful, but
17     WITHOUT ANY WARRANTY; without even the implied warranty of
18     MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
19     General Public License for more details.
21     You should have received a copy of the GNU General Public License
22     along with foam-extend.  If not, see <http://www.gnu.org/licenses/>.
24 \*---------------------------------------------------------------------------*/
26 #include "odeChemistrySolver.H"
27 #include "ODEChemistryModel.H"
29 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
31 template<class CompType, class ThermoType>
32 Foam::ode<CompType, ThermoType>::ode
34     ODEChemistryModel<CompType, ThermoType>& model,
35     const word& modelName
38     chemistrySolver<CompType, ThermoType>(model, modelName),
39     coeffsDict_(model.subDict(modelName + "Coeffs")),
40     solverName_(coeffsDict_.lookup("ODESolver")),
41     odeSolver_(ODESolver::New(solverName_, model)),
42     eps_(readScalar(coeffsDict_.lookup("eps"))),
43     scale_(readScalar(coeffsDict_.lookup("scale")))
47 // * * * * * * * * * * * * * * * * Destructor  * * * * * * * * * * * * * * * //
49 template<class CompType, class ThermoType>
50 Foam::ode<CompType, ThermoType>::~ode()
54 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
56 template<class CompType, class ThermoType>
57 Foam::scalar Foam::ode<CompType, ThermoType>::solve
59     scalarField& c,
60     const scalar T,
61     const scalar p,
62     const scalar t0,
63     const scalar dt
64 ) const
66     label nSpecie = this->model_.nSpecie();
67     scalarField& c1 = this->model_.coeffs();
69     // copy the concentration, T and P to the total solve-vector
70     for (label i=0; i<nSpecie; i++)
71     {
72         c1[i] = c[i];
73     }
74     c1[nSpecie] = T;
75     c1[nSpecie+1] = p;
77     scalar dtEst = dt;
79     odeSolver_->solve
80     (
81         t0,
82         t0 + dt,
83         eps_,
84         dtEst
85     );
87     for (label i=0; i<c.size(); i++)
88     {
89         c[i] = max(0.0, c1[i]);
90     }
92     return dtEst;
96 // ************************************************************************* //