CFD solver for laminar reacting flows with detailed kinetic mechanisms based on OpenFOAM and OpenSMOKE++ Library
- Eigen (http://eigen.tuxfamily.org/index.php?title=Main_Page)
- RapidXML (http://rapidxml.sourceforge.net/)
- Boost C++ (http://www.boost.org/)
- OpenSMOKE++ ([email protected])
- Intel MKL (https://software.intel.com/en-us/intel-mkl)
- ISATLib ([email protected])
- ODEPACK (http://computation.llnl.gov/casc/odepack/odepack_home.html)
- DVODE (http://computation.llnl.gov/casc/odepack/odepack_home.html)
- DASPK (http://www.engineering.ucsb.edu/~cse/software.html)
- Sundials (http://computation.llnl.gov/casc/sundials/main.html)
- MEBDF (http://wwwf.imperial.ac.uk/~jcash/IVP_software/readme.html)
- RADAU (http://www.unige.ch/~hairer/software.html)
Three different options are available to compile the code, according to the level of support for the solution of ODE systems. The ISATLib is needed only if you want to apply the In Situ Adaptive Tabulation (ISAT) technique.
-
Minimalist: no external, optional libraries are required. Only the native OpenSMOKE++ ODE solver can be used.
-
Minimalist + Intel MKL: only the native OpenSMOKE++ ODE solver can be used, but linear algebra operations are managed by the Intel MKL libraries
-
Complete: all the optional libraries are linked to the code, in order to have the possibility to work with different ODE solvers
-
Instructions to compile the Minimalist version
-
Open the
mybashrc.minimalist
, choose the version of OpenFOAM you are using (2.2, 2.3, 2.4, 3.0) and adjust the paths to the compulsory external libraries -
Type:
source mybashrc.minimalist
-
Go to Section 4
-
Instructions to compile the Minimalist+MKL version
-
Open the
mybashrc.minimalist.mkl
, choose the version of OpenFOAM you are using (2.2, 2.3, 2.4, 3.0) and adjust the paths to the compulsory external libraries and the paths to the Intel MKL library -
Type:
source mybashrc.minimalist.mkl
-
Go to Section 4
-
Instructions to compile the Complete version
-
Open the
mybashrc.complete
, choose the version of OpenFOAM you are using (2.2, 2.3, 2.4, 3.0) and adjust the paths to the compulsory external libraries and the Intel MKL library. You can choose the additional external libraries you want to add to laminarSMOKE, by modifying theEXTERNAL_ODE_SOLVERS
variable: in particular1
means that the support is requested, while0
means that no support is requested. Obviously, for each requested library, you need to provide the correct path. -
Type:
source mybashrc.complete
-
Go to Section 4
-
Compile the solvers
-
Compile the user-defined boundary condition library: from the
libs/boundaryConditionsOpenSMOKE++
folder typewmake
-
Compile the customized radiation library: from the
libs/radiationOpenSMOKE++
folder typewmake
-
Go to Section 5
-
Compile the solvers
- Compile the steady-state (accounting for buoyancy) solver: from the
solver/laminarBuoyantSimpleSMOKE
folder typewmake
- Compile the steady-state solver: from the
solver/laminarSimpleSMOKE
folder typewmake
- Compile the unsteady (accounting for buoyancy) solver: from the
solver/laminarBuoyantPimpleSMOKE
folder typewmake
- Compile the unsteady solver: from the
solver/laminarPimpleSMOKE
folder typewmake
- Compile the post-processor: from the
solver/laminarSMOKEpostProcessor
folder typewmake
The folder run/tutorials/ToroFlames/F3/
contains a simple test case (laminar coflow diffusion flame fed with hydrogen).
-
Unsteady simulation: Open the
laminarBuoyantPimpleSMOKE-Global
folder, build the mesh using theblockMesh
utility, and run the case using thelaminarBuoyantPimpleSMOKE
solver. Even if you are interested in steady state conditions, we strongly suggest to always start with unsteady calculations to create a reasonable first-guess solution for the application of the steady state solver. In this case, you can stop the unsteady simulation after 50 ms of physical time. -
Steady state simuation: you can now move to the
laminarBuoyantSimpleSMOKE-Global
folder. Copy the last time folder calculated by the unsteady solver (point 1 above), build the mesh using theblockMesh
utility, and run the case using thelaminarBuoyantSimpleSMOKE
solver. In order to reach the steady state conditions, 5000-6000 iterations are enough.