|FPE Extra Issue 10, October 2016|
Issue 10, October 2016
Towards Faster Wildland Fire-Spread Models
By: Alan Lattimer, Ph.D. and Brian Lattimer, Ph.D. Jensen Hughes, Inc.
From firefighter safety to the mitigation of damages in the wildland-urban interface, situational awareness is critical to the management of wildland fires. Accurate predictions of the fire spread generated in better than real-time are critical to maintaining this awareness. However, wildland fires are large, multi-scale, complex phenomena that are affected by marginally under- stood or known factors such as the weather and fuel loading. Due to the scale of the problem, fire-spread models must simplify the physics or the scale to meet the real-time requirements for generating a prediction. A third option for reducing the computation time is reduced-order modeling. This article introduces the basic reduced-order modeling technique for fire-spread models and discusses the challenges and opportunities moving forward.
The modeling and simulation of wildland fires is a complex, multi-scale problem. Several fire-spread models (such as FARSITE)1 have been developed to provide operational resource management during fire outbreaks.2-4 Physics-based, large-scale models for predicting wildland fires, such as FIRETEC and WFDS,5-11 are useful for understanding the processes that drive the spread of wildland fires, particularly into the wildland-urban interface (WUI). However, the real-time requirements and limited on-scene computational resources often preclude their use as an operational tool. To reduce the time, level-set methods that track the fire front, such as WRF-Fire,12-13 are often used to produce real-time predictions. The spread of the fire in these models is typically based on the quasi-empirical mathematical fire-spread model given by Rothermel.14
To demonstrate the efficacy of reduced-order modeling, proper orthogonal decomposition (POD) in conjunction with the discrete empirical interpolation method (DEIM)15-17 was used to predict the solution to the advection-reaction-diffusion equation given by Mandel et al.18 This technique also projects the nonlinearity, resulting in ROMs that decrease the online computational cost by 100-1000 times. This ROM technique can be applied in other large-scale fire applications, such as ventilation for mine fires. POD is the most commonly used model-reduction technique for nonlinear dynamics.19-21 This technique has previously been used for wildland fire-spread models with limited success with a speed up of approximately five times over the full-order model (FOM).22-23 By addressing the nonlinearity of the system using DEIM alongside POD, the ROM is significantly faster with little degradation in accuracy.
The details of how to create a reduced-order model are beyond the scope of this article, but the specifics of the mathematical model and how it is reduced can be found in Lattimer.24 The reader is also referred to several good resources in addition to the ones above on the specifics of
POD.25-28 The focus of this article is on POD, but Benner et al.29 and the references therein, are recommended as excellent resources for other reduced-order modeling strategies.
The phenomenological model studied in this article is suggested by Mandel et al.18 as a model to predict flame-front propagation in wildland fires. The one-dimensional (1D) and two- dimensional (2D) versions of this coupled system were studied to assess the benefits of POD with DEIM over the standard POD approach. These approaches look to project the original system into a smaller system that captures the important features of the original system. Once the smaller system has been solved, we can project it back up to the original system size.
Table 1: Results for the ROM. Solution time for the FOM was 99.1 seconds
Figure 1: FOM versus POD and POD/DEIM ROM
where rT = 200, rS = 150, and rDEIM = 250.
Figure 2: Solution times for the FOM and POD/DEIM
models when running ten simulations.
As can be seen in Figure 1, the results for each of these models was almost identical. However, as Table 1 shows, the solution time for the POD model is just over five times faster than the FOM, and the POD/DEIM model was over a hundred times faster than the FOM for those shown in Figure 1. There is an off-line computational cost to develop the ROM. The advantage is that this cost must only be paid once, and after that all models with similar input conditions can be run from the previously built ROM. To see how these time savings accrue, a graph showing the computational costs for running ten ROMs versus ten full-order models is given in Figure 2.
Figure 3: FOM versus POD/DEIM ROM for the 2D fire-spread model.
The results for the 2D model are even better. The same methodology was followed for the 2D models as was done for the 1D models, the difference being that our domain is now a 1000 meter by 1000 meter grass field. The initial conditions have a fire started in the center of the field, and the model was again run for a 3000 second simulation. In this case, the full-order model required 226 seconds to complete while the POD/DEIM ROM took 0.86 seconds. Thus, the ROM is approximately 260 times faster than the original model. Further, the original model required significantly more computational resources. Figure 3 shows a comparison of the two solutions which appear to be well within the operational needs for accuracy.
Nonlinear wildland fire-spread models such as the one given by Mandel et al.18 can be challenging to reduce due to the nonlinearity of the model. By addressing this issue using DEIM, the computational time was reduced by another order of magnitude over POD alone without sacrificing accuracy. While this model is simplistic when compared to the state-of-the-art fire-spread models, it does serve to show that reduced-order modeling can be an effective tool to reduce the computational costs without sacrificing the underlying physics and accuracy of the full-order model. It is important again to note that reduced-order modeling is not a new fire-spread model. Rather, it is mathematical technique that can make existing models faster. In this way, reduced order modeling can be the tool that moves the existing models from running on high-performance computers to laptops, or, put another way, from the laboratory to the field. There is much research to still be done to ensure the efficacy of this technique on existing fire-spread solutions, but the initial results are encouraging.
This research was partially developed under Grant No. 200-2014-59669, awarded by NIOSH. The findings and conclusions in this report are those of the authors and do not reflect the official policies of the Department of Health and Human Services; nor does mention of trade names, commercial practices, or organizations imply endorsement by the U.S. Government.
Alan Lattimer, Ph.D. and Brian Lattimer, Ph.D. are with Jensen Hughes, Inc.
 M. Finney. FARSITE: Fire area simulator: model development and evaluation. 2004.
 P. Benner, S. Gugercin, and K. Willcox. A survey of projection-based model reduction meth- ods for parametric dynamical systems. SIAM Review, 57(4):483–531, 2015.
Subscribe to FPE Extra
To subscribe to receive FPE Extra via email each month, sign up here.