QUIC-Fire_GUI
O5080 - QUIC-Fire GUI ment to run QUIC-Fire and other applications
Engineering topics
Publications and source records attributed to Linn, Rodman.
O5080 - QUIC-Fire GUI ment to run QUIC-Fire and other applications
This is the public-facing repository for the LANL-developed and maintained HIGRAD/FIRETEC software, a CFD software designed to simulate the coupled fire-atmospheric processes at a landscape-scale of resolution. Distribution is possible under OSS license O5086.
Active management of fire-dependent ecosystems for specific species leads to complex tradeoffs, which affect conservation outcomes to other species. Therefore a multi-species evaluation of management actions is required. Habitat Suitability Models (HSMs) can help in predicting the likelihood of species occurrence using corresponding environmental variables and empirical relationships that link occurrence with specific environmental conditions. Incorporating multiple species into HSMs and relating them to habitat dynamics is crucial for ecosystems that require active management with prescribed fire. To address this issue, we developed multi-species HSM driven within an existing population model of the longleaf pine-hardwood ecosystem to assess the suitability of an ecosystem given different fire management strategies and environmental conditions. The population model used in this study provides spatial and temporal changes of longleaf pine-hardwood habitat structure in response to fire. These habitat values are used by the HSM to calculate habitat suitability for three threatened and endangered faunal species of this ecosystem, which all thrive with frequent fire, but have unique habitat requirements. Transient habitat conditions are traced to predict longleaf pine ecosystem trajectories under various management strategies, thereby evaluating current land management actions, such as thinning or prescribed fire frequencies. We tested a suite of environmental conditions to emphasize the sensitivity of the species to different fire management actions. The results of our modeling suggest that maximum suitable habitat for all three species can be achieved with fire frequency occurring at approximately once every three years. The modeling results support current management actions and provide a new habitat assessment tool that incorporates ecological factors for multiple species, thus providing for habitat optimization.
Data accompanying Jonko, Yedinak, Conley, and Linn (2021): "Sensitivity of grass fires burning in marginal conditions to atmospheric turbulence", JGR-Atmospheres. LA-UR-21-24832. The following files were generated from simulations performed with Los Alamos National Laboratory coupled atmosphere-fire behavior model FIRETEC v1.3: -area_burned.csv contains the area burned per second in m^2 for 45 CORE simulations, 20 INIT simulations and 20 BOUND simulations (85 columns), and 600 seconds of simulation (600 rows). -fire_perimeter.csv containing the fire perimeters in m for 45 CORE simulations, 20 INIT simulations and 20 BOUND simulations (85 columns), and 600 seconds of simulation (600 rows). -consumption_per_second.csv containing consumption in kg/s for 45 CORE simulations, 20 INIT simulations and 20 BOUND simulations (85 columns), and 600 seconds of simulation (600 rows). -rate_of_spread.csv containing the distance of the fire front from the ignition line in m for 45 CORE simulations, 20 INIT simulations and 20 BOUND simulations (85 columns), and 600 seconds of simulation (600 rows). -CORE1_u_10m_300_sec.hdf, INIT3_u_10m_300_sec.hdf, and BOUND3_u_10m_300_sec.hdf contain u wind speeds at 10 m height for 300 seconds (3D: time, x, y) of simulation for simulations CORE1, INIT3, and BOUND3, respectively. -CORE.final.fuel.density.hdf contains final fuel densities in the surface layer of the model (3D: ensemble member, x, y) for all 45 CORE simulations. -CORE{1-45}_3D_u_120_sec.hdf: 45 files containing vertically and horizontally resolved u winds for 120 seconds of simulation (4D: time, x, y, z) for the 45 CORE ensemble members. For additional information about this data, please contact the corresponding author at ajonko@lanl.gov.
Wildfire ignition, intensity, and spread rates are tightly linked with water cycle extremes. The science of wildfire prediction has traditionally encompassed the use of physical and empirical models to quantify the direction and speed of fire spread, plume injection and fire-aerosol impacts on atmospheric composition, predictions of fire season severity on subseasonal-to-seasonal (S2S) time scales, and assessment of the spatial and temporal patterns of fire risk across landscapes. Together with expanding observation networks, machine learning and artificial intelligence (AI) have the potential to revolutionize the application of such models for fire science, saving lives, protecting critical infrastructure, and providing more accurate estimates of wildfire-climate feedbacks.
Wildfire behavior is dictated by the complex interaction of numerous physical phenomena including dynamic ambient and fire-induced winds, heat transfer, aerodynamic drag on the wind by the fuel and combustion. These phenomena create complex feedback effects between the fire and its surroundings. In this study, we aim to study the mechanisms by which buoyant flame dynamics along with vortical motions and instabilities control wildfire propagation. Specifically, this study employs a suite of simulations conducted with the physics-based coupled fire-atmosphere behavior model (FIRETEC). The simulations are initialized with a fire line and the fires are allowed to propagate on a grass bed, where the fuel heights and wind conditions are varied systematically. Flow variables are extracted to identify the characteristics of the alternating counter-rotational vortices, called towers and troughs, that drive convective heat transfer and fire spread. These vortices have previously been observed in wildfires and laboratory fires, and have also been observed to arise spontaneously in FIRETEC due to the fundamental physics incorporated in the model. However, these past observations have been qualitative in nature and no quantitative studies can be found in the literature which connected these coherent structures fundamental to fire behavior with the constitutive flow variables. To that end, a variety of state variables are examined in the context of these coherent structures under various wind profile and grass height conditions. Identification of various correlated signatures and fire-atmosphere feedbacks in simulations provides a hypothesis that can be tested in future observational or experimental efforts, potentially assisting experimental design, and can aid in the interpretation of data from in situ detectors.
Here, we performed a suite of flow simulations for a 12-wind-turbine array with varying inflow conditions and lateral spacings, and compared the impacts of the flow on velocity deficit and wake recovery. We imposed both laminar inflow and turbulent inflows, which contain turbulence for the Ekman layer and a low-level jet (LLJ) in the stable boundary layer. To solve the flow through the wind turbines and their wakes, we used a large-eddy simulation technique with an actuator-line method. We compared the time series for the velocity deficit at the first and rear columns to observe the temporal change in velocity deficit for the entire wind farm. The velocity deficit at the first column for LLJ inflow was similar to that for laminar inflow. However, the magnitude of velocity deficit at the rear columns for the case with LLJ inflow was 11.9% greater because of strong wake recovery, which was enhanced by the vertical flux of kinetic energy associated with the LLJ. In order to observe the spatial transition and characteristics of wake recovery, we performed statistical analyses of the velocity at different locations for both the laminar and LLJ inflows. These studies indicated that strong wake recovery was present, and a kurtosis analysis showed that the probability density function for the streamwise velocity followed a Gaussian distribution. In a quadrant analysis of the Reynolds stress, we found that the ejection and sweep motions for the LLJ inflow case were greater than those for the laminar inflow case.
Studying turbulence in vegetation canopies is vital in the context of a number of micrometeorological and hydrological applications. While recent focus has shifted more towards exploring different kinds of canopy heterogeneities, there are still gaps in the existing knowledge on the multiple types of dynamics involved in the case of horizontally homogeneous canopies. For example, experimental studies have indicated that turbulence in the canopy sublayer (CSL) can be divided into three regimes. In the deep-zone, the flow-field is dominated by von Kármán vortex streets and interrupted by strong sweep events. The second zone near the canopy top is dominated by attached eddies and Kelvin-Helmholtz waves associated with the velocity inflection point in the mean longitudinal velocity profile. Above the canopy, the flow resembles classical boundary layer flow. In this study, these different kinds of dynamics are studied together by means of a large eddy simulation (LES). The main theme of this work is to address the question whether the parametrization of the canopy by a distributed drag force in numerical simulations instead of placing real solid obstacles is consistent with the three layer conceptual model. Unique techniques such as measures from information theory and coupled oscillator analysis are used to extract the coherent structures associated with the two motions. It can be said that a better understanding of the rich dynamics associated with the simplest case of canopy turbulence can lead to more efficient simulations and more importantly improve the interpretation of more complex scenarios.
This is a computing progress report for w14_firetec. FIRETEC simulations will explore the prescribed fire ignition methods to achieve burning objectives (understory reduction and ecosystem health) but at the same time minimize the risk of escaped fire.
High-resolution large-eddy simulation of the flow over a large wind farm (64 wind turbines) is performed using the HIGRAD/FIRETEC-WindBlade model, which is a high-performance computing wind turbine–atmosphere interaction model that uses the Lagrangian actuator line method to represent rotating turbine blades. These high-resolution large-eddy simulation results are used to parameterize the thrust and power coefficients that contain information about turbine interference effects within the wind farm. Those coefficients are then incorporated into the WRF (Weather Research and Forecasting) model in order to evaluate interference effects in larger-scale models. In the high-resolution WindBlade wind farm simulation, insufficient distance between turbines creates the interference between turbines, including significant vertical variations in momentum and turbulent intensity. The characteristics of the wake are further investigated by analyzing the distribution of the vorticity and turbulent intensity. Quadrant analysis in the turbine and post-turbine areas reveals that the ejection motion induced by the presence of the wind turbines is dominant compared to that in the other quadrants, indicating that the sweep motion is increased at the location where strong wake recovery occurs. Regional-scale WRF simulations reveal that although the turbulent mixing induced by the wind farm is partly diffused to the upper region, there is no significant change in the boundary layer depth. The velocity deficit does not appear to be very sensitive to the local distribution of turbine coefficients. However, differences of about 5% on parameterized turbulent kinetic energy were found depending on the turbine coefficient distribution. Furthermore, turbine coefficients that consider interference in the wind farm should be used in wind farm parameterization for larger-scale models to better describe sub-grid scale turbulent processes.
We validate our proposed mechanisms with the aid of numerical simulations using a mesoscale atmospheric model and LANL's local-scale dispersion model.