OpenSMOKE++ is a general framework for numerical simulations of reacting systems with detailed kinetic mechanisms, including thousands of chemical species and reactions. The framework is entirely ...written in object-oriented C++ and can be easily extended and customized by the user for specific systems, without having to modify the core functionality of the program. The OpenSMOKE++ framework can handle simulations of ideal chemical reactors (plug-flow, batch, and jet stirred reactors), shock-tubes, rapid compression machines, and can be easily incorporated into multi-dimensional CFD codes for the modeling of reacting flows. OpenSMOKE++ provides useful numerical tools such as the sensitivity and rate of production analyses, needed to recognize the main chemical paths and to interpret the numerical results from a kinetic point of view. Since simulations involving large kinetic mechanisms are very time consuming, OpenSMOKE++ adopts advanced numerical techniques able to reduce the computational cost, without sacrificing the accuracy and the robustness of the calculations.
In the present paper we give a detailed description of the framework features, the numerical models available, and the implementation of the code. The possibility of coupling the OpenSMOKE++ functionality with existing numerical codes is discussed. The computational performances of the framework are presented, and the capabilities of OpenSMOKE++ in terms of integration of stiff ODE systems are discussed and analyzed with special emphasis. Some examples demonstrating the ability of the OpenSMOKE++ framework to successfully manage large kinetic mechanisms are eventually presented.
Program title: OpenSMOKE++
Catalogue identifier: AEVY_v1_0
Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEVY_v1_0.html
Program obtainable from: CPC Program Library, Queen’s University, Belfast, N. Ireland
Licensing provisions: GNU General Public License, version 3
No. of lines in distributed program, including test data, etc.: 146353
No. of bytes in distributed program, including test data, etc.: 4890534
Distribution format: tar.gz
Programming language: C++.
Computer: Any computer that can run a C++ Compiler.
Operating system: Tested on Microsoft Windows 7, Ubuntu 14.4.
RAM: From a few Mb to several Gb depending on the size of the system being simulated.
Classification: 22.
External routines: Eigen, Boost C++ Libraries, RapidXML
Nature of problem: Evolution of reacting gas mixtures with detailed description of thermodynamic, kinetic and transport data.
Solution method: Stiff systems of Ordinary differential Equations, whose solution is obtained using methods based on the Backward Differentiation Formulas (BDF) (LU factorization of dense matrices is required).
Additional comments: The code was specifically conceived for managing homogeneous, reacting mixtures including thousands of species and reactions.
Running time: Problem-dependent, from seconds (small kinetics) to hours
The primary objective of the present endeavor is to collect, consolidate, and review the vast amount of experimental data on the laminar flame speeds of hydrocarbon and oxygenated fuels that have ...been reported in recent years, analyze them by using a detailed kinetic mechanism for the pyrolysis and combustion of a large variety of fuels at high temperature conditions, and thereby identify aspects of the mechanism that require further revision. The review and assessment was hierarchically conducted, in the sequence of the foundational C0–C4 species; the reference fuels of alkanes (n-heptane, iso-octane, n-decane, n-dodecane), cyclo-alkanes (cyclohexane and methyl-cyclo-hexane) and the aromatics (benzene, toluene, xylene and ethylbenzene); and the oxygenated fuels of alcohols, C3H6O isomers, ethers (dimethyl ether and ethyl tertiary butyl ether), and methyl esters up to methyl decanoate. Mixtures of some of these fuels, including those with hydrogen, were also considered. The comprehensive nature of the present mechanism and effort is emphasized.
The use of simplifying techniques to obtain skeletal kinetic mechanisms with the required accuracy is often a necessary step when computationally demanding simulations are concerned. In this work, a ...novel approach for an automatic mechanism reduction, aimed at retaining accuracy on specific target species, is proposed. Starting from the consolidated coupling between flux analysis and sensitivity analysis, a methodology based on curve matching and functional data analysis was developed, through which the importance of a species in the target accuracy is assessed via a proper metric. The error associated with the removal of uncertain species from the detailed mechanism is quantified in terms of distance and similarity indices before and after such a removal, within a Species-Targeted Sensitivity Analysis (STSA) framework. A species ranking is then generated, and the original mechanism is progressively reduced. The whole algorithm also implements several improvements to enhance a faster convergence, and adds a novel criterion to remove unimportant reactions, based on sensitivity analysis to kinetic parameters.
The capability of this algorithm was tested through two case studies in this work. A kinetic mechanism for a Toluene Reference Fuel (TRF) was first obtained, with the overall reactivity as reduction target. The numerical procedure allowed to obtain a compact skeletal mechanism (115 species and 856 reactions), able to retain good accuracy in ignition delay time and laminar flame speed predictions of both fuel mixture and pure compounds. More important, two skeletal mechanisms for methane combustion, including chemistry of nitrogen oxides (NOx), were developed, with different degrees of reduction. The agreement between the original and the skeletal mechanisms in terms of NO formation was successfully assessed with satisfactory results. Attention was also dedicated to the choice of the type of reactor where undertaking reduction, which turned out to play a major role in the overall process.
•Biomass origin and structure affect composition and yield of pyrolysis biochar.•Operating conditions of thermal treatment affect composition and yield of pyrolysis biochar.•Further pyrolysis ...reactions modify young biochar structure and release gas products.•A predictive kinetic mechanism of biomass pyrolysis accounts for all these effects.•Large collection of literature experimental data supports the validation of the kinetic mechanism.
Biomass is increasingly being recognized as a promising carrier for both heat, energy and chemicals production. However, several aspects still require intense research activity towards a better design and optimization of industrial combustors, gasifiers and pyrolyzer. The objective of this work is to update the CRECK kinetic mechanism of biomass pyrolysis, allowing a better prediction of both yield and composition of the solid residue (biochar). Moreover, further model modifications allow to better describe the variability of hemicellulose in different biomass. To this end, a large set of literature experimental data is collected and organized into a database, which is used to further tune and validate the proposed kinetic mechanism. Although the kinetic model maintains the previous agreement in respect of the rate of biomass pyrolysis, formation and distribution of gas and tar products, the novelty of this work is the greater attention to the predictions of biochar yield and composition, in a wide range of operative conditions. The model describes the solid residue as a mixture of pure carbon together with lumped metaplastic compounds, which represent the whole range of oxygenated and hydrogenated groups bonded to the carbonaceous matrix. These metaplastic species are released to the gas phase with their own kinetics and describe the change of both mass loss and elemental composition of the biochar. These comprehensive predictions of biochar composition are crucial for an accurate description of the successive oxidation and gasification processes.
An interface-resolved simulation of the combustion of a fuel droplet suspended in normal gravity is presented in this work, followed by an extensive analysis on the physical aspects involved. The ...modeling is based on DropletSMOKE++, a multiphase solver developed for the modeling of droplet vaporization and combustion in convective conditions. A wide range of phenomena can be described by the model, including the interface advection, the phase-change, the combustion chemistry, non-ideal thermodynamics and multicomponent mixtures. To our knowledge, this is the most detailed simulation performed on this configuration, providing a useful theoretical and numerical support for the experimental activity on this field. A recent experimental work is used as a reference, in which a methanol droplet is suspended on a quartz fiber and ignited at different oxygen concentrations. The numerical analysis offers a detailed insight into the physics of the problem and a satisfactory agreement with the experiments in terms of diameter decay, radial temperature profiles and sensitivity to the oxygen concentration. The vaporization rate is affected by the thermal conduction from the fiber, due to the high temperatures involved. Moreover, the fiber perturbs the flame itself, providing quenching at its surface. The combustion physics is compared to the one predicted at zero-gravity, evidencing a lower standoff-ratio, a higher flame temperature and an intense internal circulation. The distribution of the species around the droplet shows (i) a local accumulation of intermediate oxidation products at the fiber surface and (ii) water absorption in the liquid phase, affecting the vaporization rate.
Display omitted
•An analysis on the evaporation of acetic acid and ethylene glycol is presented.•A multiregion approach for the fluid-solid heat transfer is adopted.•The fiber transports heat towards ...the liquid, which distribute it by convection.•The agreement with the experiments is excellent for all the cases.•Internal convection influences the preferential vaporization in the mixture.
A detailed analysis on the evaporation of acetic acid and ethylene glycol droplets is performed experimentally and numerically. The isolated droplet is positioned in a combustion chamber, suspended on a thermocouple and evaporated in buoyancy driven convection, following the thermal history throughout the droplet lifetime. The experiments provide quantitative and qualitative data on the evaporation physics of acetic acid, ethylene glycol and their mixture. The data are then modeled adopting the multiphase CFD code DropletSMOKE++, describing the flow field around the droplet, the heating rate and the evaporation process. The main novelty introduced in this work is a multiregion approach to describe the solid fiber, which allows to model the conjugate heat transfer with the liquid and the gas phase, as well as its impact on the droplet evaporation. DropletSMOKE++ results show a good agreement with the experimental data, regarding both the diameter decay and the liquid temperature, whose internal distribution in the liquid is shown to be highly affected by the heat flux from the fiber (which can contribute up to 30–40% in the total heat flux on the droplet). The effect of the thermocouple on the evaporation rate has been highlighted simulating the same experiments considering the solid as adiabatic, showing in this case a large underprediction of the vaporization rate and confirming the need of a detailed model for the fiber to correctly predict the vaporization phenomenon. The mixture evaporation has been investigated, emphasizing the importance of adopting a detailed thermodynamic model (which includes activity coefficients) and the impact of the mixture non-ideality on the evaporation process. The mixture also exhibits preferential vaporization, facilitated by the internal convection in the liquid phase.
The formation of pollutant species in turbulent diffusion flames is strongly affected by turbulence/chemistry interactions. Unsteady counterflow diffusion flames can be conveniently used to address ...the unsteady effects of hydrodynamics on the pollutant chemistry, because they exhibit a larger range of combustion conditions than those observed in steady flames.
In this paper, unsteady effects on the formation of soot (and its main precursors) and nitrogen oxides (NO
x
) are investigated by imposing harmonic oscillations on the strain rate of several counterflow diffusion flames fed with propane. Numerical results confirm that the dynamic response of each species is strongly affected by the strain rate oscillations and the characteristic time governing its chemistry. At low frequencies of imposed oscillations the soot and NO
x
profiles show strong deviations from the steady-state profile. At large frequencies a decoupling between the concentration and the velocity field is evident. In particular, the formation of soot and NO
x
is found less sensitive to velocity fluctuations for flames with large initial strain rate. The significant increase of soot and NO
x
concentrations in unsteady conditions appears to be a function of both forcing frequency and flame global strain rate. Moreover, the cut-off frequency, defined as the minimum frequency above which the strain rate oscillations have negligible effects on the formation of each species, was found to be strongly dependent on the chemical characteristic time and the flame global strain rate, but only marginally affected by the amplitude of imposed oscillations.
In this work, the performances of two recently developed finite-rate dynamic scale similarity (SS) sub-grid scale (SGS) combustion models (named DB and DC) for non-premixed turbulent combustion are a ...priori assessed based on three Direct Numerical Simulation (DNS) databases. These numerical experiments feature temporally evolving syngas jet flames with different Reynolds (Re) numbers (2510, 4487 and 9079), experiencing a high level of local extinction. For comparison purposes, the predicting capability of these models is compared with three classical non-dynamic SS models, namely the scale similarity resolved reaction rate model (SSRRRM or A), the scale similarity filtered reaction rate model (SSFRRM or B), and a SS model derived by the “test filtering” approach (C), as well as an existing dynamic version of SSRRRM (DA). Improvements in the prediction of heat release rates using a new dynamic model DC are observed in high Re flame case. By decreasing Re, dynamic procedures produce results roughly similar to their non-dynamic counterparts. In the lowest Re, the dynamic methods lead to higher errors.
In this work, recently developed finite-rate dynamic scale similarity (SS) sub-grid scale (SGS) combustion models have been a priori assessed and compared with the Eddy Dissipation Concept (EDC) and ...“no model” approaches based on a Direct Numerical Simulation (DNS) database of a temporally evolving non-premixed jet flame. Two different filter widths, one placed in the inertial range and the other in the near dissipation range, have been used. The analyses were carried out in two time instants corresponding to instants of maximum local extinction and re-ignition. Conditional averaged filtered chemical source terms, conditioned on different parameters in the composition space, have been presented. Improvements are observed using the dynamic SS models compared to the two other approaches in the prediction of filtered chemical source terms of individual species while using larger filter widths. However, discrepancies still exists using the dynamic SS model on the turbulent/non-turbulent interfaces of the jet, mainly in the prediction of the oxidizer consumption rate.
This work’s aim is the development of a simplified kinetic model for coal devolatilization, suitable for CFD applications. The detailed model of coal devolatilization, already developed and tested ...against a very large set of experiments and conditions, is too large to be implemented and used inside a CFD code.
An automatic procedure is developed to estimate the kinetic parameters and stoichiometric coefficients of a one step model. An optimization technique manages to minimize an objective function which measures the distance between the simple one step kinetic mechanism and the results of the detailed model of coal volatilization. The results show that rate parameters can be expressed as function of the coal rank, or better of its carbon content. Despite its simplicity, the one step model is able to characterize not only the weight loss, but also the relative yields of gaseous species (CO, H2O, CO2, H2, C2H4, HCN, H2S), tar species and char. The coal database discussed in this work includes 13 coals of different elemental composition, from lignite to anthracite. The agreement of the one step model with the detailed model of coal volatilization is satisfactory for both evolution rates and cumulative values. Finally, it is important to newly underline that the proposed model, based on the detailed model, is predictive and only needs the coal elemental composition (coal rank) as an input. Two sets of correlation for the distribution of the nitrogen and sulfur compounds into solid, tar and gas phases are also proposed.