Search

International audience

International audience; Eruptive fires are one category of extreme fire behaviour. They are characterized by a sudden and unpredictable change in the fire behaviour which represents an extreme danger for people involved in firefighting. The major point is about the mechanism that turns a usual fire behaviour into an eruptive fire behaviour. Among the different explanations found in the literature, the pioneering interpretation consisting in a feedback effect caused by the convective flow induced by the fire under wind and/or slope conditions, has never been disproved with an example of fire accident. The main goal of this work lies in proposing a physical modelling of this fire induced wind. This modelling attempt is derived from the brand-new version of the Balbi model, which is a simplified physical model for surface fires at the field scale that explicitly depends on the triangle of fire (fuel bed, wind and slope). This work is a first step to the modelling of fire eruption. The model tries to represent accurately the acceleration of the fire rate of spread propagating on different sloped terrain under no-wind or weak wind conditions. It is tested against three sets of experiments carried out at the laboratory scale without external wind and against a high intensity experimental fire spreading on a steep sloped terrain and conducted under weak wind conditions in the north-western of Corsica. Some statistical tools are used to compare predicted and observed rate of spread (NMSE, Normalized Mean Square Error and MAPE, Mean Absolute Percentage Error) and to understand the model’s under-predictions or over-predictions trends (FB, Fractional Bias).

Dimensioning a fuelbreak remains always a challenging problem. For a long time, this problem was tackled using an empirical approach from the experience of operational users such as the fire fighters and the foresters. During the last decades, new approaches coming from fire safety engineering have completed the set of tools adapted to study this problem. These tools are all based on physical considerations, more or- less sophisticated. The simplest ones, consist in assimilating the flame as a radiant panel, calculating the distribution of radiant heat flux as a function of the distance separating the flame to a potential target and defining at what distance this heat flux reached a critical threshold level susceptible to produce damages on this target (pain for people or ignition for materials). The most complex ones, consist in solving the conservation equations (mass, momentum, energy ...) governing the behaviour of complex coupled problem formed by the vegetation, the flame front and the surrounding atmosphere. This new generation of engineering tool, based on CFD approach allows to directly predict the behaviour of a fire front propagating toward a fuelbreak, in order to evaluate its efficiency as a function of the amount of surface fuel (grass, shrubs) removed to reduce locally the fuel load and therefore the intensity of an incoming fire. These two approaches are fully complementary, only the first one has the potentiality to be spread operationally on the field, whereas the second one can contribute to improve the first one and to study with more detail some very sensitive situations such as those encountered in the wildland urban interface (WUI). The main part of this study concerns numerical simulations of the propagation of a fire front through a homogeneous vegetation layer (a grassland) in the vicinity of a fuelbreak represented by a band more or less wide inside which all the fuel was removed. The simulations were performed using a fully physical wildfire model (FIRESTAR3D), three variable parameters were considered in this study: the 1m open wind speed (U1 ranged between 3 and 10 m/s), the fuel height (HFuel ranged between 0.25 and 1m) and the fuelbreak width (LFB). With these conditions, the simulations covered a large range of values of the Byram’s convective number NC (0.3 < NC < 60) in order to explore wind as well driven fires (NC < 2) and plume dominated fires (NC > 10). The 72 simulations carried out in this study have been classified in three categories: 1/ Propagation (if the fire has crossed the fuelbreak with a propagation after); 2/ Overshooting or Marginal (if the fire has crossed the fuelbreak without a propagation after); 3/ No-propagation (if the fuelbreak has stopped the fire). The main objective of this study was to determine the optimal fuelbreak width LFBx separating between the Propagation and the No-propagation regimes, in order to generalize the conclusion, the results have been presented in dimensionless form (similitude theory) in representing as an example the ratio LFBx/HFuel versus the Byram’s convective number NC.

International audience; The effectiveness of a fuelbreak, created in a homogeneous grassland on a flat terrain, was studied numerically. The analysis relies on 3D numerical simulations that were performed using a detailed physical-fire-model (FIRESTAR3D) based on a multiphase formulation. To avoid border effects, calculations were carried out by imposing periodic boundary conditions along the two lateral sides of the computational domain, reproducing that way a quasi-infinitely long fire front. A total of 72 simulations were carried out for various wind speeds, fuel heights, and fuelbreak widths, which allowed to cover a large spectrum of fire behaviour, ranging from plume-dominated fires to wind-driven fires. The results were classified in three main categories: 1- “Propagation” if fire crossed the fuelbreak with a continuous fire front, 2- “Overshooting” and “Marginal” if fire marginally crosses the fuelbreak with the formation of burning pockets, and 3- “No propagation” if fire does not cross at all the fuelbreak. The ratio of fuelbreak width to fuel height, marking the “Propagation”/“No propagation” transition, was found to be scaled with Byram's convection number Nc as 75.07 × Nc−0.46. The numerical results were also compared to an operational wildfire engineering tool (DIMZAL) dedicated to fuelbreaks dimensioning.

International audience; The safety during prescribed burnings could be achieved by conducting these operations under marginal conditions of fire propagation. This type of fire can or cannot propagate on account of small deviations of the burning conditions, mainly the wind speed, the fuel load, and the fuel moisture-content. In this context, numerical simulations of grassland fires were conducted under marginal conditions in order to relate the moisturecontent threshold of propagation success to the wind speed and the fuel load. The simulations were conducted using FireStar2D, a complete physical 2D fire simulator based on a multiphase modelling approach. The 10 mopen wind speed ranged from 0 to 10 m/s and the fuel load varied from 0.1 kg/m 2 to 0.7 kg/m 2. The effects of wind speed and fuel moisture-content on the fire behaviour and on the flame parameters are discussed. The results show that the moisture threshold increases with the fuel load until it reaches a value beyond which there is no dependence. A similar dependence of the moisture threshold on the wind speed is also observed. Finally, empirical formulae were constructed to relate the fuel moisture content threshold to the wind speed and the fuel loading implicitly through Byram's convective number.

Numerical simulation of multi-physical processes requires a lot of processor time, especially when solving ill-conditional linear systems arising in fluid dynamics problems. This paper is devoted to the development of efficient parallel methods for such systems for FireStar3D wildfire modeling code. Two alternative approaches are discussed and analyzed, based on the MILU-preconditioned conjugate gradient method and on the algebraic multigrid, respectively. The main difficulties of parallelizing these methods are considered and solutions are presented: in the first case, nested twisted factorization with a staircase pipelining, and in the second, a multicolor technique for a new smoother for strongly anisotropic grids. A novel quasi-geometric interpolation technique is presented for solving the problem of positive off-diagonal matrix entries in the multigrid. The limits of applicability of the methods are determined depending on their flexibility, robustness and parallelization capabilities. The performance comparison demonstrates the superiority of the new methods over the widely used variants of the traditional conjugate gradient method.

International audience; 12 This study reports 3D numerical simulations of the ignition and the propagation of 13 grassland fires. The mathematical model is based on a multiphase formulation and on a 14 homogenization approach that consists in averaging the conservation equations (mass, 15 momentum, energy …) governing the evolution of variables representing the state of the 16 vegetation/atmosphere system, inside a control volume containing both the solid-17 vegetation phase and the surrounding gaseous phase. This preliminary operation results 18 in the introduction of source/sink additional terms representing the interaction between 19 the gaseous phase and the solid-fuel particles. This study was conducted at large scale in 20 grassland because it represents the scale at which the behavior of the fire front presents 21 most similarities with full scale wildfires and also because of the existence of a large 22 number of relatively well controlled experiments performed in Australia and in the 23 United States. The simulations were performed for a tall grass, on a flat terrain, and for 24 six values of the 10-m open wind speed ranged between 1 and 12 m/s. The results are in 25 fairly good agreement with experimental data, with the predictions of operational 26 empirical and semi-empirical models, such as the McArthur model (MK5) in Australia and 27 the Rothermel model (BEHAVE) in USA, as well as with the predictions of other fully 3D 28 physical fire models (FIRETEC and WFDS). The comparison with the literature was 29 mainly based on the estimation of the rate of fire spread (ROS) and of the fire intensity, 30 as well as on the analysis of the fire-front shape. 31 32

International audience; A 3D multi-physical model referred to as “FireStar3D” has been developed in order to predict the behavior of wildfires at a local scale (

International audience; Flow simulation models often contain open boundaries in order to avoid computations on too large or infinite domains. The open boundary conditions partially drives the flow and the influence of the outer domain has to be taken into account on such boundaries. " Free " open boundary conditions are then essential to the model. The goal of this work is to introduce a new open boundary condition, adapted both for inlet and outlet, parametrized and driven by simple informations such as flow rate.

International audience; This study reports 3D numerical simulations of the ignition and the propagation of grassland fires. The mathematical model consists in solving the conservation equations governing the behavior of the coupled system formed by the vegetation and the ambient atmosphere. Based on a finite volume discretization, the computation code is parallelized using OpenMP directives and had been the subject of numerous validations. the model was previously tested on small scale in case of homogeneous litter fires; in this study it is proposed to test it on larger scale in case of grassland fires. the results are in fairly agreement with experimental data, with predictions of operational empirical-models developed in Australia (MK5) and in the USA (BEHAVE), as well as with predictions of physical models (FIRETEC, WFDS, FireStar2D). The comparison with the literature is based on the estimation of the rate of fire spread (ROS) and the analysis of the fire-front shape.

Special issue: LES of turbulence aeroacoustics and combustion; International audience; This study deals with the problem of turbulent atmospheric boundary-layer flow over a forest canopy. Numerous previous works showed that this flow presents more similarities with a mixing-layer flow than with the standard boundary-layer flow. In this paper, this problem was studied for homogeneous canopies, using large eddy simulation (LES). The numerical results reproduced correctly the various steps of development of this flow: the appearance of a first generation of coherent structures resulting from the development of a primary Kelvin-Helmholtz instability, the reorganization of these structures, by vortex pairing and kinking, the development of a secondary instability and the formation of horseshoe vortices. Then, the process of transport of a passive scalar from a forest canopy into a clear atmosphere was studied in two cases, i.e., when the passive scalar concentration at the surface foliage is either constant or time-varying. Even though this small difference has little influence on the concentration patterns, the results showed that it can significantly affect the concentration magnitude as well as the dynamics of the total concentration in the atmosphere.

International audience; This paper deals with 3D numerical simulations of two fires fronts (head and backfire) propagating simultaneously through a grassland fuel. The simulations were carried out using a "fully" physical and three-dimensional fire model (namely WFDS). One of the objectives of this work, was to evaluate the potential for fully physical fire model to simulate the interactions between two fire fronts (a head fire and a backfire), in conditions similar to those encountered during suppression fire operations. A set of numerical simulations was first carried out for standalone head fires propagating through grasslands on a flat terrain and for various wind conditions ranging between 1 and 10 m/s. These results were compared with experimental data and numerical results from the literature. The same calculations were then repeated, with a backfire ignited at the downwind side of the plot. The numerical results highlighted that, for these particular conditions, head fire and backfire can interact, mainly, via two mechanisms: - at relatively large distances (greater than 10 m) the head fire acts on backfire as a screen and reduces the direct action of the wind flow on the backfire, - at relative small distances (nearly equal to 10 m) the gas flow (entrainment) generated in the vicinity of the head fire promotes the aspiration of the backfire towards the main fire front.

International audience; One of the objectives of this paper was to simulate numerically the interaction between two line fires ignited in a grassland, on a flat terrain, perpendicularly to the wind direction, in such a way that the two fire fronts (a head fire and a backfire) propagated in opposite directions parallel to the wind. The numerical simulations were conducted in 3D using the new fuel element module recently implemented in WFDS. We chose a grassland for the fuel layer, because it represents certainly one of the most studied ecosystem, in the frame of experimental campaigns to study the behaviour of surface fires at large scale. The aim of this numerical study is to understand what are the physical phenomena and the favorable conditions of ignition of a counter fire, during a fire suppression operation. The preliminary results highlighted that the two fire fronts interacted only at a relative short distance (10-20 m), following this scenario: • The thermal plume above the head fire (representing the main fire front) formed a sort of shelter, protecting the backfire to the direct effect of the wind flow • Before the merging between the two fire fronts, an in-draft flow was observed ahead of the head fire, promoting and accelerating the propagation of the backfire During the last step of the merging of the two fire fronts, a sudden increase of the heat release rate was observed, indicating a significant interaction between the two fires, which can potentially represent a safety problem for people in charge of this kind of operation.

International audience; This paper deals with the numerical simulation (using a LES approach) of the interaction between an atmospheric boundary layer (ABL) and a canopy, representing a forest cover. This problem was studied for a homogeneous configuration representing the situation encountered above a continuous forest cover, and a heterogeneous configuration representing the situation similar to an edge or a clearing in a forest. The numerical results, reproduced correctly all the main characteristics of this flow, as reported in the literature: the formation of a first generation of coherent structures aligned transversally from the wind flow direction, the reorganisation and the deformation of these vortex tubes to horse shoe structures. The results obtained, introducing a discontinuity in the canopy (reproducing a clearing or a fuel break in a forest), were compared with experimental data collected in a wind tunnel. The results confirmed the existence of a strong turbulence activity inside the canopy at a distance equal to 8 times the height of the canopy, referenced in the literature as an Enhance Gust Zone (EGZ) characterized by a local peak of the skewness factor.

The propagation of grassland fires is simulated using a fully physical based model, partially developed during the FIRESTAR European Union programme. This approach, based on a multiphase formulation, includes the calculation of the degradation of the vegetation (by dehydration and pyrolysis) and the turbulent/reactive flow resulting from the mixing between the ambient gas (wind flow) and the pyrolizate. The solid fuel is simulated as homogeneous distribution of solid particles forming a porous media, interacting with the gas flow using a continuous distribution of drag forces. Other source terms representing the interactions between the vegetation and the gas flow are also taken into account, such as the production of water vapour and gaseous fuel, the radiation of soot particles and ashes, and the convective exchange in the energy balance. The model was validated from preliminary calculations carried out at small scale, for a homogeneous fuel bed (pine needles, excelsior, sticks) and compared with experimental results obtained in a wind tunnel. Calculations are then extended to study the propagation of fires through a flat grassland, for various wind speed conditions. The numerical results are compared to empirical and semi-empirical predictions obtained in similar conditions.

International audience; Most of the models proposed in literature for binary diffusion coefficients of solids in supercritical fluids are restricted to infinite dilution; this can be explained by the fact that most of experimental data are performed in the dilute range. However some industrial processes, such as supercritical fluid separation, operate at finite concentration for complex mixtures. In this case, the concentration dependence of diffusion coefficients must be considered, especially near the upper critical endpoint (UCEP) where a strong decrease of diffusion coefficients was experimentally observed. In order to represent this slowing down, a modified version of the Darken equation was proposed in literature for naphthalene in supercritical carbon dioxide. In this paper, the conditions of application of such a modelling are investigated. In particular, we focus on the order of magnitude of the solubility of the solid and on the vicinity of the critical endpoint. Various equations proposed in literature for the modelling of the infinite dilution diffusion coefficients of the solutes are also compared. Ten binary mixtures of solids with supercritical carbon dioxide were considered for this purpose.

The aim of this article is to study the propagation of a wildfire, propagating through two eco-systems characteristics of the Mediterranean region located in south of France: the garrigue and the aleppo pines forest. The behavior of a line fire is analysed using a physical model based on the resolution of conservation equations (mass, momentum, energy) governing the coupled system formed by the vegetation and the surrounding atmosphere. The effects of wind (ranged from 3 to 30 m/s) and slope of the terrain (ranged from 0 to 40°) have been studied for surface fires propagating through shrubs (quercus coccifera) composing the surface vegetation of a garrigue and the understorey vegetation of a young pines stand (pinus halepensis). Then the configuration has been generalized, in adding a solid fuel layer above the first one, representing the canopy of a young pines stand. The numerical results have highlighted the strong connection between the surface fire and the canopy. In all cases (on flat terrain) the a...

International audience; The Sanchez-Lacombe equation of state is known to describe the thermodynamic properties of molecular fluids of arbitrary size, especially polymer-solvent phase behavior. However, it is rarely used for modeling solid-supercritical fluid equilibria. In this work, it is shown that a proper estimation of the EoS characteristic parameters together with a thermodynamically consistent expression of fugacity coefficients allows a satisfactory correlation of the solubility of solids in the supercritical phase. Binary mixtures containing carbon dioxide, ethane, ethylene and xenon were considered for this purpose. In a first step, the consistency of experimental data was checked using variance analysis. Then, different mixing rules were considered and results compared with those obtained with the Peng-Robinson equation. Finally, the lower and upper boundaries of the solid-liquid-vapor regions (LCEP and UCEP) were also determined and compared with experimental values.

Free oscillations of a viscous liquid drop surrounded by a dynamically inactive ambient gas, in zero gravity, are investigated numerically using FIDAP TM package in the axisymmetrical case. The full Navier-Stokes equations with appropriate interfacial conditions are solved by using Galerkin/Finite element technique along with the spine method for the advection of the free boundary are used. The aim of this preliminary study is to demonstrate the ability of the package to accurately solve nonlinear free surface problems. Oscillations of viscous drops released from an initial static deformation of small-to large-amplitude proportional to the second spherical harmonic, without initial internal circulation, are considered