Search

A model examining underground smouldering peat combustion is presented. A one-step chemical reaction is considered where the gas and solid are assumed to be in thermal equilibrium. The full model allows porosity, permeability and gas density to vary and considers a buoyant velocity field determined by Darcy's law. Due to the low bulk thermal conductivity of peat, the diffusion of oxygen through it is characterised by a Lewis number much less than one. This results in thermal-diffusive instabilities. These instabilities can cause flame balls to arise in gaseous combustion and a fingering regime to arise in solid combustion. Analytical solutions to simplified spherically symmetrical equations are derived. These equations assume diffusion to be the dominant transport mechanism as well as taking that the porosity, gas molecular weight and gas density all remain constant. The underlying structure of the combustion region is found to be analogous to that of a flame ball. When studied in cylindrical symmetry a single, stable finger can be modelled propagating against an imposed air flow. The effects of heat losses, velocity magnitude and the Lewis number can be studied and results are compared to existing experimental smouldering combustion data. Although no detailed experiments have studied this phenomenon in peat, predicted results capture key qualitative trends found in both filtration combustion of polyurethane foam and in the fingering combustion of paper. In addition to this, when the imposed air flow is reduced to zero a propagating combustion front is predicted, analogous to a self-travelling flame ball. When the velocity field is determined by Darcy's law the dimensionless permeability of the peat plays a key role in determining the range of values over which fingering combustion can occur. Whilst there is little impact of taking the gas molecular weight to be constant, when porosity is allowed to vary and a relationship between porosity and permeability is included an over-blowing extinction limit is identified. This limit is not found in the constant-porosity model where a low-fuel extinction limit is predicted. Peats of differing ages and locations can possess significantly different characteristics. However, the fingering regime is predicted to occur within the range of parameters in which peat soils lie. Experiments suggest that fingering combustion can take the form of both sparse fingers and a complex fingering regime. The cylindrically symmetrical model can not capture tip-splitting. Hence the model does not explicitly account for the distance between two neighbouring fingers. However, an estimate for this value can be made if peat smouldering were to occur in a regime of multiple fingering. An averaged continuum model describing the spread of an ember storm is also presented. The dominant mechanism determining the spread-rate of the fire is the lofting and landing of embers and individual fires are taken to grow in an elliptical manner under the influence of the wind. When an ember storm is spreading at a steady speed, its spread rate is found to be described by a single similarity solution.

Large eddy simulation of spray combustion in an HSDI engine is carried out in this thesis. The implementation of the code was performed in logical steps that allowed both assessment of the performance of the existing KIVA-LES and later development. The analysis of the liquid annular jet confirmed existence of typical, annular jet exclusive structures like head vortices, stagnation point and recirculation in the inner zone. The influence of the swirl in the ambient domain was found to have profound impact on the development, penetration and radial spreading of the jet. Detailed results were reported in Jagus et al. (2008). The code was further validated by performing an extensive study of large eddy simulation of diesel fuel mixing in an engine environment. The reaction models were switched off in order to isolate the effects of both swirl and the different numerical treatment of LES. Reference RANS simulation was performed and significant differences were found. LES was found to capture much better the influence of the swirl on the liquid and vapour jets, a feature essentially absent in RANS results. Moreover, the predicted penetration of the liquid was much higher for the LES case and more in accordance with experimental measurements. Liquid penetration and subsequent evaporation are very important for prediction of heat release rates and encouraging results formed a good basis to performing a full simulation with models for ignition and combustion employed. The findings were analyzed in the paper by Jagus et al. (2009). Further modifications were introduced into the LES code, among them changes to the combustion model that was originally developed for RANS and calculation of the filter width. A new way of estimating the turbulent timescale (eddy turnover time) assured that the incompatibilities in the numerical treatment were eliminated and benefits of LES maximized. The new combustion model proved to give a very good agreement with experimental data, especially with regard to pressure and accumulated heat release rates. Both qualitative and quantitative results presented a significant improvement with respect to RANS results and old LES formulation. The new LES model was proved to give a very good performance on a spectrum of mesh resolutions. Encouraging results were obtained on a coarse mesh sets therefore proving that the new LES code is able to give good prediction even on mesh sizes more suitable for RANS. Overall, LES was found to be a worthy alternative to the well established RANS methods, surpassing it in many areas, such as liquid penetration prediction, temperature-turbulence coupling and prediction of volume-averaged data. It was also discovered that the improved LES code is capable of producing very good results on under-resolved mesh resolutions, a feature that is especially important in industrial applications and on serial code structure.