Application of a laminar flamelet model to confined explosion hazards

The majority of computational studies of confined explosion hazards apply simple and inaccurate combustion models, requiring ad hoc corrections to obtain realistic flame shapes and often predicting an order of magnitude error in the overpressures. This work describes the application of a laminar flamelet model to a series of two-dimensional test cases. The model is computationally efficient applying an algebraic expression to calculate the flame surface area, an empirical correlation for the laminar flame speed and a novel unstructured, solution adaptive numerical grid system which allows important features of the solution to be resolved close to the flame. Accurate flame shapes are predicted, the correct burning rate is predicted near the walls, and an improvement in the predicted overpressures is obtained. However, in these fully turbulent calculations the overpressures are still too high and the flame arrival times too low, indicating the need for a model for the early laminar burning phase. Due to the computational expense, it is unrealistic to model a laminar flame in the complex geometries involved and therefore a pragmatic approach is employed which constrains the flame to propagate at the laminar flame speed. Transition to turbulent burning occurs at a specified turbulent Reynolds number. With the laminar phase model included, the predicted flame arrival times increase significantly, but are still too low. However, this has no significant effect on the overpressures, which are predicted accurately for a baffled channel test case where rapid transition occurs once the flame reaches the first pair of baffles. In a channel with obstacles on the centreline, transition is more gradual and the accuracy of the predicted overpressures is reduced. However, although the accuracy is still less than desirable in some cases, it is much better than the order of magnitude error previously expected.

Global linear and nonlinear stability of viscous confined plane wakes with co-flow

The global stability of confined uniform density wakes is studied numerically, using two-dimensional linear global modes and nonlinear direct numerical simulations. The wake inflow velocity is varied between different amounts of co-flow (base bleed). In accordance with previous studies, we find that the frequencies of both the most unstable linear and the saturated nonlinear global mode increase with confinement. For wake Reynolds number Re = 100 we find the confinement to be stabilising, decreasing the growth rate of the linear and the saturation amplitude of the nonlinear modes. The dampening effect is connected to the streamwise development of the base flow, and decreases for more parallel flows at higher Re. The linear analysis reveals that the critical wake velocities are almost identical for unconfined and confined wakes at Re ≈ 400. Further, the results are compared with literature data for an inviscid parallel wake. The confined wake is found to be more stable than its inviscid counterpart, whereas the unconfined wake is more unstable than the inviscid wake. The main reason for both is the base flow development. A detailed comparison of the linear and nonlinear results reveals that the most unstable linear global mode gives in all cases an excellent prediction of the initial nonlinear behaviour and therefore the stability boundary. However, the nonlinear saturated state is different, mainly for higher Re. For Re = 100, the saturated frequency differs less than 5% from the linear frequency, and trends regarding confinement observed in the linear analysis are confirmed.

The local and global stability of confined planar wakes at intermediate Reynolds number

At high Reynolds numbers, wake flows become more globally unstable when they are confined within a duct or between two flat plates. At Reynolds numbers around 100, however, global analyses suggest that such flows become more stable when confined, while local analyses suggest that they become more unstable. The aim of this paper is to resolve this apparent contradiction by examining a set of obstacle-free wakes. In this theoretical and numerical study, we combine global and local stability analyses of planar wake flows at $\mathit{Re}= 100$ to determine the effect of confinement. We find that confinement acts in three ways: it modifies the length of the recirculation zone if one exists, it brings the boundary layers closer to the shear layers, and it can make the flow more locally absolutely unstable. Depending on the flow parameters, these effects work with or against each other to destabilize or stabilize the flow. In wake flows at $\mathit{Re}= 100$ with free-slip boundaries, flows are most globally unstable when the outer flows are 50 % wider than the half-width of the inner flow because the first and third effects work together. In wake flows at $\mathit{Re}= 100$ with no-slip boundaries, confinement has little overall effect when the flows are weakly confined because the first two effects work against the third. Confinement has a strong stabilizing effect, however, when the flows are strongly confined because all three effects work together. By combining local and global analyses, we have been able to isolate these three effects and resolve the apparent contradictions in previous work.