Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Leading Edge”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

Phenomena Identification and Ranking Tables (PIRT) analysis of wind turbine blade leading edge erosion

Leading edge erosion (LEE) of wind turbine blades has been identified as a major factor in decreased wind turbine blade lifetimes and energy output over time. Accordingly, the International Energy Agency Wind Technology Collaboration Programme (IEA Wind TCP) has created the Task 46 Phase 2 to undertake cooperative research in the key topic of blade erosion. The purpose of Task 46 Erosion Phase 2 is to further improve understanding of erosion driving factors, develop datasets and model tools to enhance prediction of leading-edge erosion likelihood, identify damage at the earliest possible stage and advance potential solutions. The scope of work covers several technical areas, reflecting the multidisciplinary nature of the challenge. Participants in the task are given in Table 1.

17 WIND ENERGY↗

Spatiotemporal Variability in Wind Turbine Blade Leading Edge Erosion

Wind turbine blade leading edge erosion (LEE) reduces energy production and increases wind energy operation and maintenance costs. Degradation of the blade coating and ultimately damage to the underlying blade structure are caused by collisions of falling hydrometeors with rotating blades. The selection of optimal methods to mitigate/reduce LEE are critically dependent on the rates of coating fatigue accumulation at a given location and the time variance in the accumulation of material stresses. However, no such assessment currently exists for the United States of America (USA). To address this research gap, blade coating lifetimes at 883 sites across the USA are generated based on high-frequency (5-min) estimates of material fatigue derived using a mechanistic model and robust meteorological measurements. Results indicate blade coating failure at some sites in as few as 4 years, and that the frequency and intensity of material stresses are both highly episodic and spatially varying. Time series analyses indicate that up to one-third of blade coating lifetime is exhausted in just 360 5-min periods in the Southern Great Plains (SGP). Conversely, sites in the Pacific Northwest (PNW) exhibit the same level of coating lifetime depletion in over three times as many time periods. Thus, it may be more cost-effective to use wind turbine deregulation (erosion-safe mode) for damage reduction and blade lifetime extension in the SGP, while the application of blade leading edge protective measures may be more appropriate in the PNW. Annual total precipitation and mean wind speed are shown to be poor predictors of blade coating lifetime, re-emphasizing the need for detailed modeling studies such as that presented herein.

Pryor, Sara C. (ORCID:0000000348473440)↗

Healable Coatings as a Mechanism to Repair Leading Edge Erosion in Wind Energy

Wind turbine blades are highly engineered structures designed to face temperature extremes and high winds. However, erosion of the blade's leading edge and subsequent repair remains a significant and costly challenge for the wind energy industry. Repair of these leading edges can lead to large amounts of downtime for the turbine and significant operational inefficiencies. In this work, the strength of adhesion and healing ability of a commercially available vitrimer (Mallinda's VITRIMAX) was compared to that of a thermoplastic resin, which has previously been demonstrated in wind energy applications (Arkema's Elium) to evaluate their efficacy as surface coatings for wind turbine blades, particularly their leading edges. Vitrimers are a class of inherently reprocessable thermosets, and it was theorized that vitrimer-based leading edge coatings could enable more robust and efficient wind turbine blades with decreased operational downtime and safer maintenance practices. It was found that the VITRIMAX adhered better to the wind blades' surfaces than both the manufacturer's paint and Elium, with increases in pull-off strength of adhesion ranging from 24% to 83% above that of the original paint. Furthermore, the VITRIMAX adhered strongly to the underlying composite of each blade with strength of adhesion values increasing in ranges from 42% to 97% above that of the original paint. Finally, the vitrimer coating showed an 88% decrease in surface roughness compared to end-of-life blade materials, and initial healing demonstrations in which coatings were manually scratched and subsequently healed exhibited an ~84.5% decrease in scratch depths.

Hubbard, Amber [ORNL]↗

A variational method for the sheath potential of hypersonic leading edges with space-charge limitations

Electron transpiration cooling for the leading edges (LE) of hypersonic aircraft utilizes thermionic emission; however, space-charge effects limit the electron emission rate, potentially diminishing the efficiency of this cooling mechanism. We develop a variational weak form of the Poisson equation that describes the sheath potential and then numerically solve it using the finite element method. This formulation has two main benefits: (1) the space-charge limit condition can be incorporated as a constraint and (2) it allows for the analysis of three-dimensional geometries with complex boundary conditions. We demonstrate that the current emitted from the surface of an LE is generally a small fraction of the Child–Langmuir limit due to space charge. We then propose several methods to enhance the emitted current from the surface and to boost the cooling effect of thermionic emission. These include increasing the plasma density, applying a negative surface potential, and using fringe fields under suitable geometric conditions. For a LaB6 emitting LE, the total emitted current is shown to be minimal and independent of the temperature of a surface with floating potential. However, when a negative potential is applied and the surface is heated, the emitted current follows the Richardson–Dushman relationship up to a critical temperature, beyond which it remains constant. At an applied surface potential of −5 V, the critical temperature is around 1700 K.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

VitriEdge: Repairable & Durable Vitrimer Coatings for Wind Turbine Blade Leading Edges

The primary goal of this Level 1b incubator project was to explore the use of vitrimer coatings for repair of leading-edge erosion on end-of-life wind turbine blade surfaces, beyond coating strength of adhesion which has previously been demonstrated in the Level 1a project. Uniform vitrimer coatings (thickness: 400 µm) were applied to two end-of-life wind turbine blades for flexural, fatigue, and laminate tensile testing where the addition of the coating did not produce any statistical variation in tensile properties with minor drops in flexural strength for some laminate formulations. However, a <2% variation in storage modulus was measured for laminate structures (i.e., blade samples with vitrimer coatings) across 100,000 flexural cycles and upon laminate tensile failure, the vitrimer coatings displayed no visible signs of delamination. In addition, three methods to heal vitrimer coating damage was displayed: oven heating, addition of hot water, and a forced convection heat gun. All three heating and healing mechanisms demonstrated significant healing with scratch depths decreasing between 79-91% at healing times ranging between 1-min and 10-minutes. Finally, a water jet machine was used to simulate rain erosion for both the blade surfaces and vitrimer-coated blade surfaces where the diameter and depth of the damage was recorded as a function of exposure time, water pressure, height of exposure, and angle of exposure. Of interest, while the vitrimer coating did not significantly lessen the damage experienced during rain erosion, the addition of vitrimer composite coatings(5 wt.% mica addition) did result in a crack-resistant, durable coating capable of self-healing behavior and in all cases the angle of rain exposure was the most critical parameter explored. It is crucial to continue exploring this space where vitrimer coatings are of interest for both their self-healing properties and potential use as reversible adhesives.

17 WIND ENERGY↗

Investigating the effect of B T direction on W source-to-core pathways during the SAS-VW campaign on DIII-D

Experiments using the V-shaped closed slot tungsten (W) coated SAS-VW divertor in DIII-D studied the effects of the BT direction on core contamination of eroded tungsten from a closed slot divertor configuration. Core W content is inferred using soft-X ray tomography (SXR) and vacuum ultraviolet spectroscopy (SPRED), while W divertor erosion is inferred from visible spectroscopy of W emission (400.9 nm) measured by in-slot filterscopes (filtered photo-multipliers). Post-mortem analysis from the campaign discovered tile misalignment leading to suspected pronounced leading-edge erosion in the unfavorable BT direction (ion B ⃗×∇B ⃗ drift away from divertor) likely not captured by diagnostics. However, empirical findings show up to ~2-3x larger core contamination in the favorable BT direction even considering no additional W erosion from leading edges. A “source-to-core efficiency factor” is derived to estimate the effects of leading-edge erosion and compare W contamination for two pairs of H-mode discharges in opposite BT directions. While having differing absolute parameters, similar core impurity density gradients suggest comparable core impurity transport. These results show that favorable BT may have stronger source-to-core pathways for W impurities sourced from the outer divertor region. Possible explanations could include the effects of E ⃗×B ⃗ drifts on W transport in the scrape-off-layer (SOL) as well as previously determined fast SOL inner target directed flows in favorable BT.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Reduced order modeling of a fluidized bed particle receiver for concentrating solar power with thermal energy storage

Oxide particles can serve as both the heat transfer and thermal energy storage (TES) media for next-generation concentrating solar power (CSP) plants where high-temperature TES enables dispatchable electricity from efficient power cycles with firing temperatures above 600 °C. Transferring heat to flowing particles at such high temperatures in a MW-scale central tower receiver remains a challenge for the CSP community. For indirect receivers with external walls to contain the particles, maintaining wall temperatures below the limits of structural metal alloys requires high heat transfer coefficients between the wall and the moving particle stream. Bubbling fluidization of downward-flowing particles can sustain high bed-wall heat transfer coefficients (> 1000 W m -2 K -1 ). Using experimentally calibrated correlations for bed-wall heat transfer and vertical particle dispersion, this study implements an axially discretized zonal model of a counterflow fluidized bed receiver to explore how bubbling fluidization may enable indirect cavity particle receivers. High bed-wall heat transfer coefficients support solar fluxes on angled cavity walls > 200 kW m -2 at peak aperture fluxes of 980 kW m -2 while maintaining external wall temperatures < 950 °C. Lateral particle dispersion enables hotter particles near the receiver leading edge to mix with cooler particles further from the leading edge to lower maximum external wall temperatures. Parametric studies identify how mass fluxes, particle dispersion, and solar concentrations impact indirect receiver thermal efficiency and uniformity for a CSP plant. These studies provide a basis for the design of indirect fluidized-bed cavity receivers that can maintain particle outlet temperatures for TES above 750 °C.

14 SOLAR ENERGY↗

Comparison of juvenile white sturgeon blade strike survival through a conventional and a novel hydropower turbine

Hydropower is a globally important renewable energy source with profound impacts on aquatic organisms, including direct injury and mortality to fish that move through hydropower facilities. While it is expected that the characteristics of turbine design and operation influence the type, frequency, and magnitude of injuries suffered by entrained fish, few direct comparisons exist to illustrate differences in fish survival outcomes for conventional and novel turbine designs. This study evaluated blade strike-associated injury and mortality rates for juvenile white sturgeon (Acipenser transmontanus) passed through a model-scale turbine equipped with a runner having conventional blade profiles (thin, straight leading edges) and the same turbine equipped with a runner designed for improved fish survival (thick, slanted leading edges). In both trials, head and turbine runner rotational speeds were matched to produce five blade peripheral speeds between 15.0 and 27.6 m/s. High-speed video was captured for all turbine passage events. Fish were assessed for injuries following passage, and mortalities were assessed immediately after passage and after 48 h. Conventional runner passage resulted in 42% to 78% survival after 48 h, with approximately one-third of all tested fish killed by severing. Under the same test conditions, immediate and 48 h survival rates through the novel runner were 100% except at the highest speed condition (95.6% survival at 48 h). These results clearly indicate that turbine design has a profound effect on fish survival, and suggest that thoughtful redesign of hydropower equipment could significantly ameliorate the risks of blade strike mortality.

13 HYDRO ENERGY↗

Aerodynamic Characterization of 3D Scanned Wind Turbine Blades Using Experimental and Computational Methods

This study presents an aerodynamic characterization of 3D scanned wind turbine blades using both experimental and computational methods. The research was conducted by Gulf Wind Technology and Sandia National Laboratories. The primary objective was to investigate the aerodynamic impacts of leading-edge manufacturing defects on wind turbine blades. The study utilized the Stratasys NEO 800 3D Printer for high-precision manufacturing and the GWT Accelerator Wind Tunnel for experimental testing. Computational simulations were performed using COMSOL Multiphysics to model the wind tunnel and analyze flow characteristics and OpenFOAM to study the aerodynamic impacts of leading-edge defects. OpenFAST was used to estimate how these defects can lead to revenue losses for wind farm operators as high as 6%. The results demonstrated significant aerodynamic performance variations due to defects, with detailed analysis provided through wind tunnel and CFD data. The findings contribute to the understanding of defect impacts on wind turbine blade performance and offer insights for future design improvements.

17 WIND ENERGY↗

A Flang Plugin for Fortran Feature Characterization

As new compute systems are developed, there is still a need to compile and execute codes authored in Fortran on these leading edge systems. In order to achieve this, development of compilers that support the latest hardware is continuously under development. Though the specification of Fortran is extensive, it is helpful to compiler authors to be able to prioritize the development of key features in order to get certain codes deemed important, e.g., applications of interest to leadership computing facilities, executable on leading edge compute systems. Identifying key features though is largely done through querying software experts or users of the Fortran applications of interest, who then manually report what features are and are not present. This exercise can both time consuming and error prone. To automate this process, we present a compiler plugin to Flang, the Fortran frontend for LLVM. This plugin is a tool that operates on the parse tree representation generated by Flang and detects key features based on walking parse tree nodes that correspond to features of interest. We show the result of our tool on four applications, three of which were manually profiled by software experts. We show the discrepancies between our tool and the manual characterization of the three applications, as well as generate a characterization for an application not yet profiled. We intend to open-source our tool in order to invite the community to benefit from the tool and make contributions for other features.

Cabrera, Anthony [ORNL]↗

DUNE Phase II: scientific opportunities, detector concepts, technological solutions

The international collaboration designing and constructing the Deep Underground Neutrino Experiment (DUNE) at the Long-Baseline Neutrino Facility (LBNF) has developed a two-phase strategy toward the implementation of this leading-edge, large-scale science project. The 2023 report of the US Particle Physics Project Prioritization Panel (P5) reaffirmed this vision and strongly endorsed DUNE Phase I and Phase II, as did the European Strategy for Particle Physics. While the construction of the DUNE Phase I is well underway, this White Paper focuses on DUNE Phase II planning. DUNE Phase-II consists of a third and fourth far detector (FD) module, an upgraded near detector complex, and an enhanced 2.1 MW beam. The fourth FD module is conceived as a “Module of Opportunity”, aimed at expanding the physics opportunities, in addition to supporting the core DUNE science program, with more advanced technologies. This document highlights the increased science opportunities offered by the DUNE Phase II near and far detectors, including long-baseline neutrino oscillation physics, neutrino astrophysics, and physics beyond the standard model. It describes the DUNE Phase II near and far detector technologies and detector design concepts that are currently under consideration. A summary of key R&D goals and prototyping phases needed to realize the Phase II detector technical designs is also provided. DUNE's Phase II detectors, along with the increased beam power, will complete the full scope of DUNE, enabling a multi-decadal program of groundbreaking science with neutrinos.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Bio-inspired alula-based winglet design for enhanced heat transfer in high temperature fin-and-tube heat exchangers

Fin-and-tube heat exchangers (FTHEs) are widely used for high-temperature flue-gas heat recovery, but their performance is often limited by wake regions and non-uniform fin-surface temperatures. This study proposes and numerically evaluates four bio-inspired longitudinal vortex generator (VG) configurations in a high-temperature FTHE with flue-gas inlet temperature ∼1230 K: double-delta, curved double-delta, alula, and a new curved-alula geometry. The reference fin is not hydraulically plain; it already incorporates leading-edge separation columns and convex protrusions, so the alula-type winglets are assessed as downstream add-ons acting on a strongly disturbed flow. In a second step, perforations (one, two and three circular holes) are introduced into the curved-alula VGs to further tailor the flow field. Three-dimensional simulations with the Shear Stress Transpor (SST) $k - ω$ model, temperature-dependent flue-gas properties and conjugate conduction are carried out for gas-side Reynolds numbers $Re_g ≈ 8.0$ x $10^2 - 3.6$ x $10^3$ (mass flow rates 0.5 – 2.5 g/s), and the designs are compared in terms of surface heat flux, Nusselt number, friction factor and hydrothermal performance factor (HTPF). For this already-promoted fin, the additional downstream winglets provide moderate, incremental hydrothermal gains. At the highest Reynolds number, the best non-perforated design (curved-alula) increases surface heat flux from 1630.9 to 1794.7 kW/m² (∼ 10 % gain) and the Nusselt number from 227.6 to 242.6 (∼ 7 % gain), while the friction factor rises from 0.26 to about 0.30, yielding HTPF values close to unity (∼ 0.9 – 1.0). Introducing circular perforations into the curved-alula winglets acts mainly as a wake-bleeding refinement: the three-hole configuration provides a heat flux of 1824.7 kW/m² and a pressure drop of 127.9 Pa, with HTPF in the range ∼ 1.03 – 1.14 and a small (∼ 1 – 3 %) improvement over the solid curved-alula design. Flow-field analysis shows that the perforated curved-alula VGs shrink tube-wake regions, thin the thermal boundary layer and homogenize the fin-surface temperature (outlet-gas temperature ∼ 510 – 520 K and fin-surface temperature ∼ 420 – 421 K for the three-hole case). An optimal flue-gas mass flow rate of ∼ 1 g/s ($Re_g ≈ 1.5$ x $10^3$) is identified, beyond which additional heat-transfer gains are offset by rapidly increasing pressure losses. Overall, the results highlight that initial fin geometry and VG placement are as important as VG shape: alula-based winglets are expected to yield larger relative gains on simpler flat-fin layouts or when positioned closer to the fin leading edge and tube

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

The trailing edge is everywhere: tree rings reveal the transient risk of extinction hidden inside climate envelope forecasts

Given the importance of climate in shaping species’ geographic distributions, climate change poses an existential threat to biodiversity. Climate envelope modeling, the predominant approach used to quantify this threat, presumes that individuals in populations respond to climate variability and change according to species-level responses inferred from spatial occurrence data—such that individuals at the cool edge of a species’ distribution should benefit from warming (the “leading edge”), whereas individuals at the warm edge should suffer (the “trailing edge”). Using 1,558 tree-ring time series of an aridland pine (Pinus edulis) collected at 977 locations across the species’ distribution, we found that trees everywhere grow less in warmer-than-average and drier-than-average years. Ubiquitous negative temperature sensitivity indicates that individuals across the entire distribution should suffer with warming—the entire distribution is a trailing edge. Species-level responses to spatial climate variation are opposite in sign to individual-scale responses to time-varying climate for approximately half the species’ distribution with respect to temperature and the majority of the species’ distribution with respect to precipitation. These findings, added to evidence from the literature for scale-dependent climate responses in hundreds of species, suggest that correlative, equilibrium-based range forecasts may fail to accurately represent how individuals in populations will be impacted by changing climate. A scale-dependent view of the impact of climate change on biodiversity highlights the transient risk of extinction hidden inside climate envelope forecasts and the importance of evolution in rescuing species from extinction whenever local climate variability and change exceeds individual-scale climate tolerances.

54 ENVIRONMENTAL SCIENCES↗

Tree rings reveal the transient risk of extinction hidden inside climate envelope forecasts

Given the importance of climate in shaping species’ geographic distributions, climate change poses an existential threat to biodiversity. Climate envelope modeling, the predominant approach used to quantify this threat, presumes that individuals in populations respond to climate variability and change according to species-level responses inferred from spatial occurrence data—such that individuals at the cool edge of a species’ distribution should benefit from warming (the “leading edge”), whereas individuals at the warm edge should suffer (the “trailing edge”). Using 1,558 tree-ring time series of an aridland pine (Pinus edulis) collected at 977 locations across the species’ distribution, we found that trees everywhere grow less in warmer-than-average and drier-than-average years. Ubiquitous negative temperature sensitivity indicates that individuals across the entire distribution should suffer with warming—the entire distribution is a trailing edge. Species-level responses to spatial climate variation are opposite in sign to individual-scale responses to time-varying climate for approximately half the species’ distribution with respect to temperature and the majority of the species’ distribution with respect to precipitation. These findings, added to evidence from the literature for scale-dependent climate responses in hundreds of species, suggest that correlative, equilibrium-based range forecasts may fail to accurately represent how individuals in populations will be impacted by changing climate. A scale-dependent view of the impact of climate change on biodiversity highlights the transient risk of extinction hidden inside climate envelope forecasts and the importance of evolution in rescuing species from extinction whenever local climate variability and change exceeds individual-scale climate tolerances.

54 ENVIRONMENTAL SCIENCES↗

Atomically resolved edges and defects in lead halide perovskites

Although edges and defects constitute only a small fraction of crystalline materials, they exert an outsized impact on a material′s properties. Organic–inorganic halide perovskites are promising next-generation semiconductor materials with superior cost effectiveness and interesting optoelectronic properties. However, clear images of their edges have remained challenging to obtain owing to their extreme sensitivity. Using truly high-speed ultralow-dose four-dimensional scanning transmission electron microscopy with dose fractionation, we perform ptychography at, to our knowledge, the lowest-dose atomic resolution to date, revealing not only the detailed atomic structure of the edges of a halide perovskite but also their structural dynamics. Here, a majority methylammonium (MA) and iodine (I) edge termination is observed in methylammonium lead iodide (MAPbI 3 ), and the damage rate of its edges and internal defects is found to depend on the concentration and type of vacancies present, with a preponderance of I vacancies in particular correlating with higher rates of damage.

4D-STEM↗

Semiconductor-compatible topological digital alloys

Recently, GeSn alloys have attracted much interest for direct-gap infrared photonics and as potential topological materials which are compatible with the semiconductor industry. However, for photonics, the high-Sn content required leads to low detectivity, associated with poor material quality, and the (>35 %) Sn required for topological properties have been out of reach experimentally. Here, we demonstrate that by patterning the Sn distribution within Ge, the electronic properties have a far greater tunability than is possible with the random alloy. For the GeSn δ-digital alloy (DA) formed by confining Sn atoms in atomic layer(s) along the [111] direction of Ge, we show that ∼ 10 % Sn can lead to a triple-point semimetal. These findings are understood in terms of Sn ordering causing spatial separation of Sn and Ge band edges, leading to band inversion. This mechanism can also lead to a weak topological insulator, Weyl semimetal, and enables tunable direct bandgaps down to 2 meV, covering the entire infrared range. This DA induced topological properties are also identified in compound semiconductors, such as InAs 1-x Sb x , showing the general applicability of the DA design for realizing topological properties on conventional semiconductor platforms. Furthermore, our findings not only point to a new class of currently unexplored topological systems accessible by epitaxy, but also establish the promise of low-Sn GeSn DAs for application as infrared laser diodes and photodetectors in Si photonic integrated circuits and infrared image sensors.

3D topological semimetal↗

Electron-Phonon Renormalization in the Proton-Conducting Electrolyte Ba Zr O 3 and Its Implications for High-Temperature Electrolysis

Electrical leakage is an inherent problem in solid-oxide fuel and electrolyzer cells, limiting their energy-conversion efficiency. High concentrations of electrons or holes can exacerbate this issue. However, it is largely unclear how the typical high operating temperatures of the systems influence their charge-carrier concentrations. In this work, we use first-principles calculations to examine how lattice vibrations impact electrical conductivity in conventional electrolytes, using Ba Zr O 3 as a representative material. Our analysis shows that phonon-induced shifts in the band gap and band edges lead to a dramatic increase in p -type carrier concentrations at temperatures above 600 K, compared to models that neglect temperature effects. Additionally, we reveal the importance of oxygen-ion motion on band-edge positions, which makes valence-band-edge shift dominate the band-gap change. Our study provides a protocol for calculating phonon-induced changes in similar oxides, paving the way for interrogating electrical leakage in electrolytes for high-temperature operation. Published by the American Physical Society 2025

25 ENERGY STORAGE↗

Characterizing Turbulence at a Forest Edge: Comparing Sub-Filter Scale Turbulence Models in Simulations of Flow over a Canopy

In wildfires, atmospheric turbulence plays a major role in the transfer of turbulent kinetic energy. Understanding how turbulence feeds back into a dynamical system is important, down to the varying small scales of fuel structures (i.e. pine needles, grass). Large eddy simulations (LES) are a common way of numerically representing turbulence. The Smagorinsky model (1963) serves as one of the most studied sub-grid scale representations in LES. In this investigation, the Smagorinsky model was implemented in HIGRAD/FIRETEC, LANL’s coupled fire-atmosphere model. This study was motivated by the need to quantitatively investigate the vorticity budget equation in HIGRAD/FIRETEC. The Smagorinsky turbulent kinetic energy (TKE) was compared to FIRETEC’s 1.5-order TKE eddy-viscosity subgrid-scale model, known as the Linn turbulence model. This was done in simulations of flow over flat terrain with a homogeneous, cuboidal canopy in the center of the domain. Examinations of the modeled vertical TKE profile and turbulent statistics at the leading edge, and throughout the canopy, show that the Smagorinsky model provides comparable results to that of the original closure model posed in FIRETEC.

58 GEOSCIENCES↗