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.

Explosion pressure prediction via polynomial mathematical correlation based on advanced CFD Modelling

Computational Fluid Dynamics CFD can be used as a powerful tool supporting engineers throughout the steps of the design. The combination of CFD with response surface methodology can play an important role in such cases. During the conceptual engineering design phase, a quick response is always a matter of urgency. During this phase even a sketch of the geometrical model is rare. Therefore, the utilisation of typical response surface developed for congested and confined environment rather than CFD can be an important tool to help the decision making process, when the geometrical model is not available, provided that similarities can be considered when taking into account the characteristic of the geometry in which the response surface was developed. The present work investigates how three different types of response surfaces behave when predicting overpressure in accidental scenarios based on CFD input. First order, partial second order and complete second order polynomial expressions are investigated. The predicted results are compared with CFD findings for a classical offshore experiment conducted by British Gas on behalf of Mobil and good agreement is observed for higher order response surfaces. The higher order response surface calculations are also compared with CFD calculations for a typical offshore module and good agreement is also observed. © 2011 Elsevier Ltd.

Initial phase modelling in numerical explosion applied to process safety

The utilisation of computational fluid dynamics (CFD) in process safety has increased significantly in recent years. The modelling of accidental explosion via CFD has in many cases replaced the classical Multi Energy and Brake Strehlow methods. The benefits obtained with CFD modelling can be diminished if proper modelling of the initial phase of explosion is neglected. In the early stages of an explosion, the flame propagates in a quasi-laminar regime. Proper modelling of the initial laminar phase is a key aspect in order to predict the peak pressure and the time to peak pressure. The present work suggests a modelling approach for the initial laminar phase in explosion scenarios. Findings are compared with experimental data for two classical explosion test cases which resemble the common features in chemical process areas (confinement and congestion). A detailed analysis of the threshold for the transition from laminar to turbulent regime is also carried out. The modelling is implemented in a fully 3D Navier-Stokes compressible formulation. Combustion is treated using a laminar flamelet approach based on the Bray, Moss and Libby (BML) formulation. A novel modified porosity approach developed for the unstructured solver is also considered. Results agree satisfactorily with experiments and the modelling is found to be robust. © 2013 The Institution of Chemical Engineers.