Deposition Characteristics of Firebrands on and Around Rectangular Cubic Structures

The focus of the present work is on the deposition of firebrands in a flow over a rectangular cubic block representative of a structure in wildland-urban interface (WUI). The study was carried out by physics based modeling where the wind flow turbulence was dealt with by large eddy simulation (LES) and firebrands were treated by Lagrangian tracking. The Lagrangian equations coupled with the flow solver, accounted for both translational and rotational motions as well as thermochemical degradation of firebrands, assumed to be cylindrical. The dimensions of the structure were varied from 3 to 9 m in the simulations for a parametric study. The simulations were carried out by tracking many firebrands randomly released with a uniform distribution from a horizontal plane 35 m above the ground into the computational domain. The coordinates of the deposited firebrands were used to calculate their normalized number density (number of landed firebrands per unit surface area) to quantify their deposition pattern. On the leewardside of the block, an area, referred to as the safe zone, was identified right behind the structure where firebrands never deposit. The size of the safe zone in the direction perpendicular to the wind was nearly identical to the width of the structure. The length of the safe zone in the wind direction was proportional to the height of the structure. The leeward face of the blocks was never hit by a firebrand. The windward face was hit by many more firebrands than the lateral faces but much less than the top face. The distribution of the number density of the deposited firebrands on the top face was found to be correlated with the flow separation and reattachment on this face.


INTRODUCTION
A critical mechanism for the spread of large outdoor fires, e.g., wildland-urban interface (WUI) fires is spotting. Spotting is the creation of the secondary (spot) fires by firebrands that are generated by the primary fires. Firebrands can be lofted up into the atmosphere and carried away by the ambient wind to short/long distances (Tarifa et al., 1967;Sardoy et al., 2008). In the presence of strong ambient winds, firebrands can cross distances from a few 100 m to a few kilometers, thus capable of spreading fires over barriers such as rivers, lakes, hills, etc. Spotting is seen frequently in WUI fires and can burn down many WUI structures under extreme conditions such as an ember shower (Manzello, 2014). This motivated the present computational study with a focus on characterizing the deposition pattern of firebrands carried by the wind on top and in the vicinity of a structure shaped as a rectangular cuboid mounted on the ground. The computational configuration here can be considered as a simplified representation of a single isolated WUI structure.
There have been several studies on the role of firebrands in the spread of wildland and WUI fires. Manzello et al. (2007) performed experiments by burning two Douglas-fir trees with 2.6 and 5.2 m heights. They found that the generated firebrands were predominately cylindrical in shape with an average diameter of 3 mm and length of 40 mm for the shorter tree and 4 and 53 mm for the taller tree. Manzello et al. (2008) constructed an apparatus capable of generating glowing firebrands and used it to release firebrand in a wind tunnel. The firebrands released in the wind-tunnel at 9 m/s experienced a mass loss of 20-40% when compared to firebrands released in no wind condition. Tohidi and Kaye (2017b), Tohidi and Kaye (2017a) experimentally and computationally studied the lofting of firebrands in a wind tunnel where in addition to wind, a convective plume was included. They observed that for higher wind speeds, the change in the initial vertical velocity of the convective column did not affect the mean or standard deviation of the heights where the firebrands lofted or the distances they traveled to land. Yin et al. (2003), Oliveira et al. (2014) developed numerical models for the firebrand transport accounting for the drag, lift and gravitational forces and their effect on the rotation of firebrands to model both translational and rotational motions of cylindrical firebrands. To validate their model, Oliveira et al. (2014) performed computations and experiments for a cylindrical firebrand (balsa wood) falling from an elevated point under a no ambient flow condition. The influence of different formulations for the distance between center of pressure and center of mass of a cylindrical object in motion was explored in the modeling by Rayleigh (1876), Marchildon et al. (1964), Rosendahl (2000), Yin et al. (2003). Anand et al. (2018) preformed simulations to investigate the deposition of cylindrical firebrands released in a turbulent wind environment from a fixed elevated point. They assumed for firebrands to retain their mass from release to landing. They reported a bivariate Gaussian function like distribution for the landed firebrand position with a larger variance in the streamwise direction, compared to the spanwise direction. Anand (2018) performed similar simulations while allowing firebrands to experience mass loss due to thermal degradation, taking into account the effect of burning. They observed that, firebrands with a higher mass density (570 kg/m 3 ) experienced a higher mass loss, as compared to lower density (230 kg/m 3 ) firebrands. The lower density firebrands cooled rapidly and reached ambient temperature before landing. On the other hand, the higher density firebrands retained more thermal energy while flying, thus had higher temperatures at landing. Song et al. (2017) performed wind tunnel experiments with disc-shape firebrands and showed the deposited firebrands had uni-modal distribution except for certain wind speed and firebrand conditions where they displayed a bimodal distribution.
Properties of the flow over a cubic obstacle mounted on the ground have been studied in the past (Murakami et al., 1987;Werner and Wengle, 1993;Lee and Bienkiewicz, 1997;Rodi, 1998). One of the earliest works is due to Murakami et al. (1987) who simulated a cube submerged in a boundary layer using largeeddy simulation (LES). Werner and Wengle (1993), Rodi (1998) computationally studied a cube mounted on a surface in a channel flow with a Reynolds number of O(10 4 ) based on the velocity at the height of the cube. Werner and Wengle (1993), Rodi (1998) showed the existence of a horse-shoe vortex on the windward side of the cube and flow separation and reattachment on the top face of the cube. Rodi (1998), using different turbulence models, reported two counter rotating recirculation region on the leeward side of the cube. Vortex shedding was observed originating from the lateral faces with a pair of re-circulation region closer to these faces. More recently, Richards et al. (2001) claimed the pressure coefficient on the surface of the cube is independent of the Reynolds number via a field test. Later in wind-tunnel experiments (Richards et al., 2007), they observed a drop in the pressure coefficient on the windward and leeward faces of the cube as the wind direction changed from 90 + to 45 + with respect to the windward face of the cube. Lim et al. (2009) performed experiments and simulations for a flow around a cube submerged in a turbulent atmospheric surface layer (ASL) and showed that the mean profiles of pressure coefficient and velocity components are independent of the Reynolds number.
The present work is a modeling study focused on deposition of firebrands in a flow over a cubic block representative of a structure in WUI. The flow is dealt with by LES while the deposition of firebrands is treated in the Lagrangian framework. In Section 2, modeling approaches are illustrated for both firebrands and the flow. In Section 3, results are presented with the model validation results included. Concluding remarks are made in Section 4.

Firebrand Equations
The firebrand equations are expressed and solved in the Lagrangian framework. Firebrands are assumed to be cylinders with a large ratio of length to diameter, undergoing both translational and rotational motions (Yin et al., 2003;Oliveira et al., 2014;Anand et al., 2018) and thermal degradation as a result of pyrolysis and charring (Morvan and Dupuy, 2004;Anand, 2018).

Translational Motion
The position and velocity of the center of mass of the firebrand are denoted by x → p and V → p , respectively, which are governed by where d/dt is the time derivative calculated in the Lagrangian framework. Eq. 2 is an expression of Newton's second law where the forces are due to gravity combined with buoyancy F → G , drag F → D and lift F → L , which are calculated (Hoerner, 1965) by: Re 0.778 α 1 + 0.1076Re 0.778 α for 0.1 < Re α ≤ 6 × 10 3 , 1.1 for 6 × 10 3 < Re α ≤ 2 × 10 5 .
Here, ρ p , D p , l, V, m p and V → rel are the firebrand density, diameter, length, volume, mass and velocity relative to the flow at the center of mass of the particle, respectively. It is calculated by where the first term indicates the flow velocity at the position of the center of mass of the firebrand. The drag coefficient C D is calculated, using the particle Reynolds number Re α D p ρ V → rel sinα/μ (Kelbaliyev, 2011), where ρ gas is the density of air and α is the incidence angle between the relative velocity and the major axis of the cylindrical firebrand z r .

Rotational Motion
The rotational motion is described by the Euler rotation equation: where I x′ , I y′ and I z′ are the moments of inertia with respect to the Cartesian frame of reference x ′ -y ′ -z ′ attached to the cylindrical firebrand with the origin at the cylinder center and the z′ axis constituting the cylinder axis. The total torque is the addition of the torque T → ′ hydro due to the hydrodynamic forces and the torque T → ′ resist (Oliveira et al., 2014) due to the frictional air resistance experienced by the firebrand T x′ resist −ρ gas ω x′ ab 4 ⎡ ⎣ 0.538 + 3.62 ρ gas a ω x′ b T y′ resist −ρ gas ω y′ ab 4 ⎡ ⎣ 0.538 + 3.62 ρ gas a ω y′ b x cp l(90 − α) 480.
Here, a D p /2 is the radius of the firebrand, b l/2 is the half length and x cp is the distance between the center of pressure and the center of mass (Marchildon et al., 1964), and A is the transformation matrix expressed in terms of quaternions ϵ 1 , ϵ 2 , ϵ 3 and η (Yin et al., 2003): Quaternions are governed by The quaternions are correlated with Euler angles (ϕ, ψ, θ) through the following equations, which are used here to find initial values of the quaternions:

Mass and Temperature
Heat is transfered from the firebrand to the surrounding gas through thermal radiation and convection. The firebrand undergoes thermal degradation and loses mass as a result of pyrolysis and char oxidation. To take this effect into account, the firebrand model assumes for the firebrand to be thermally thin (i.e. temperature throughout the firebrand is spatially uniform) with a mass governed by: where _ m pyr and _ m char are the mass loss rates due to pyrolysis and char oxidation, respectively, which are modeled by the Arrhenius equation: where m i represents the mass of the solid constituent, namely i pyr for the charring of the fuel and char for char oxidation, A i is the pre-exponential factor, T p is the temperature of the firebrand and T i E i /R is the activation temperature where E i is the activation energy. The pre-exponential factor and activation temperature for pyrolysis are A pyr 725 s −1 , T pyr 6899 K (for Pinus) and for char oxidation are A char 430 m/s, T char 9000 K ( Morvan and Dupuy, 2004;Sardoy et al., 2007;Anand, 2018). The firebrand temperature is governed by where Δh pyr 418 kJ/kg and Δh char 12 × 10 3 kJ/kg are the enthalpy of pyrolysis and char oxidation, respectively (Sardoy et al., 2007;Mell et al., 2009;Anand, 2018). Here, _ q c and _ q r are the rates of the convective and radiative heat transfer, respectively: where A is the surface area of the firebrand, T ∞ is the ambient temperature, h is the heat transfer coefficient, σ is the Stefan-Boltzmann constant and ϵ is the emmisitivity of the firebrand set to 0.9. It is noted that for improved modeling of the mass loss and thermal energy, combustion models are needed in addition to the char oxidation representation here to more accurately represent the burning effect.

Computational Approach
Our group developed a model that handles the transport and burning of firebrands, according to Eqs. 1-26, in the framework of Fire Dynamic Simulator (FDS, version 6.7.0) (McGrattan et al., 2018). FDS is computational fluid dynamics (CFD) based software capable of modeling the fire dynamics while representing significant thermal, chemical and physical processes such as combustion, turbulence, radiation, etc. In the present study, only the fluid dynamical features of FDS are relevant. Turbulence is dealt with by LES in FDS with the default option of Deardoff model (Deardorff, 1980) set to represent the subgrid-scale (SGS) terms here. FDS uses Wall-Adapting Local Eddy-viscosity model (WALE) (Nicoud and Ducros, 1999) as the near-wall model by default. The firebrand equations are solved by a second-order Admas-Bashforth time integration method, as described by Anand et al. (2018) and Anand (2018).
, t] defined in Section 2.1.1 is calculated via a trilinear interpolation of the flow velocities at cell faces to the location of center of mass of the firebrand. The coupling of firebrands to the flow solver is one-way, as the influence of firebrands on the flow is assumed negligible. The firebrands deposited on the solid surfaces, i.e., ground and faces of the block, are removed from the simulation after their deposition coordinates are recorded.

Firebrand Model Validation
To validate the firebrand model, first, a firebrand drop test previously investigated both experimentally and computationally (Oliveira et al., 2014) was considered. The exercise involved a non-burning cylindrical firebrand made from balsa wood with diameter 10 mm and length 80 mm, which was released from the height 8.7 m in a no-wind condition. At the release point, the firebrand had zero velocities and made an angle of 60 + with the vertical axis. The firebrand mass density was reported ρ p 215.5 kg/m 3 . Using the firebrand model illustrated in Section 2.1.1 and Section 2.1.2, the drop test was simulated here in a computational domain 1.5 × 1.5 × 9 m (length × width × height). In lieu of Eq. 15 (Marchildon et al., 1964), other formulas ( Table 1) have been also reported in the literature (Rayleigh, 1876;Rosendahl, 2000;Yin et al., 2003) for calculation of x cp . This motivated a sensitivity study of the model to these formulas to be a part of this validation exercise. Table 1 tabulates the landing times calculated in the current study using various x cp formulas and compares them against those obtained in the modeling and measurement of Oliveira et al. (2014). Corresponding trajectories of the firebrand from release to landing are shown in Figure 1. Both table and figure suggest the significance of the x cp formula in the firebrand landing time and trajectory. Discussed by Oliveira et al. (2014) was the notable difference between the amplitudes of the trajectory oscillation in their model (panel E in Figure 1) and their measurement (panel F). They additionally argued that this difference was correlated with the difference between their corresponding calculated and measured landing times, as tabulated in Table 1. On the other hand, Figure 1 suggests that the amplitude obtained in the current simulations, regardless of the formula used x cp , was significantly more consistent with the experimental data of Oliveira et al. (2014). When the x cp formulas of Rayleigh (1876) (panel A) and Marchildon et al. (1964) (panel B) were used, the amplitudes of the trajectories were slightly larger than those observed in the experiment and accordingly, the calculated landing times were slightly smaller than the measured landing time. When the formula of Rosendahl (2000) (panel C) was used, the calculated amplitude seemed to be more consistent with the amplitude in the experiment. However, the calculated landing time was greater than the measured landing time by a larger amount. When the formula of Yin et al. (2003) (panel D) was used in the calculations, the resulting amplitude was larger than both the measured amplitude and the amplitude measured by other formals. However, the landing time was closer to the measured landing time.

Flow Model Validation
The flow model used here was first validated against the previous experimental and modeling data obtained in a wind tunnel for a flow over a cubic block (Lim et al., 2009). The test section of the wind tunnel had dimensions of 4.5 × 0.9 × 0.6 m (length, width and height, respectively) with a cube of height of 0.08 m situated 2.36 m from the inlet of the tunnel.  Lim et al. (2009) reported that they conducted their experiments for Reynolds numbers in the range between 18, 600 and 73, 100 but did not find the mean and variance of measured velocities to significantly change at this range of Reynolds numbers. The lateral and top boundaries were set to be free slip and the outflow boundary condition was set to be open. At the inlet, turbulence with the intensity of 5% was introduced.
Flow turbulence was dealt with by LES with the Deardoff SGS (Deardorff, 1980) and near-wall models, as discussed in Section 2.2. However, the simulations were repeated with other SGS models including constant Smagorinsky (Smagorinsky, 1963), dynamic Smagorinsky (Germano et al., 1991;Moin et al., 1991), Vreman (Vreman, 2004) and RNG (Yakhot et al., 1989) available in FDS. It was determined that the results were negligibly sensitive to the SGS models. Hence, only the results of the Deardoff model are presented here. Figure 3 shows the mean velocity streamlines at a slice y 0 and z 0.5h obtained from present simulations. This figure shows the key flow structures around the cube, viz. the center of the horseshoe vortex, the flow separation and reattachment on the top and lateral faces, flow reattachment on the leeward side of the cube, the two counter rotating re-circulation region and the stagnation point of the windward face of the cube. Table 2 compares the locations of these points of interest obtained in the current study with those obtained in the simulation of Lim et al. (2009). The center of the horseshoe vortex obtained here is a little further away from the windward face of the cube when compared to the previous simulation (Lim et al., 2009). On the other hand, the locations of the stagnation point on the windward face of the cube, the reattachment length on the top face of the cube and the reattachment length on the leeward side of the cube obtained here closely match those in the simulation of Lim et al. (2009), as seen in Table 2. Figure 4 shows the pressure coefficient C p on the axial (i.e. y/h 0) and transverse (i.e. x/h 0.5) center-lines on the faces of the block as indicated in Figure 2. As could be seen Figure 4, the pressure coefficient calculated here for the top 1 | Landing time of a cylindrical firebrand released in a still air in the previous experiment and simulation (Oliveira et al., 2014), and present simulations using different formula for the center of pressure x cp (Rayleigh, 1876;Marchildon et al., 1964;Rosendahl, 2000;Yin et al., 2003).

Landing time (s)
Formula of x cp /l References  face of the cube compares very well against the experimental and simulation data of Lim et al. (2009). The agreement between the current simulation and the previous works for this coefficient is reasonable for the rest of the faces. The experimental data of Richards et al. (2001) is also shown here for a comparison albeit they were obtained for a different Reynolds number of 4.1 × 10 6 . The profile in Figure 4A shows the largest positive pressure on the windward face of the cube closer to the leading edge which is a result of the cube blocking the flow. On the top face, the largest negative pressure right after the leading edge is associated with the flow separation at the leading edge which is followed pressure recovery corresponding to the flow reattachment.
In Figure 5, the mean streamwise and vertical velocities are plotted vs. z on the axial centerline of the top face of the block at various x's. The agreement between the current simulation and the previous experiment and simulation (Lim et al., 2009) is very good. The change of the velocity profile in the x direction is attributed to the flow separation on the top face. Figure 6 shows the profiles of the root mean square (rms) of the streamwise and vertical velocities as well as the Reynolds shear stress at various x's on the axial center-line of the top face of the block. As could be seen in Figures 6A,B, the current simulation substantially over-predicts the rms values obtained in the previous experiments and the simulation (Lim et al., 2009). On the other hand, the Reynolds shear stress in the simulation is in reasonably good agreement with the previous experimental and simulation data. Figure 7 shows the computational configuration used in the simulation of firebrand deposition in a flow over a single cubic structure. The length, width and height of the structure are indicated by L, W and H, which are its dimensions in the x, y and z directions, respectively. Simulations were carried out for structures with various lengths, widths and heights. The domain size is 75 × 36 × 36 m in the x, y and z directions, respectively. The domain is divided into two sub-domains with a finer gird (0.15 × 0.15 × 0.15 m) between heights 0-12 m and a coarse grid (0.3 × 0.3 × 0.3 m) between heights 12-36 m. The inlet flow velocity was specified by a power law with an exponent of 0.18 with a velocity of 6 m/s at a reference height h 3 m which resulted in Re h 1 × 10 6 . The turbulent intensity at the inlet was set to 20%. This inlet boundary condition is an approximate   2 | The locations of the center of the horseshoe vortex (HSV) (x HVC , y HVC ); the stagnation point on the windward face y stag ; the flow reattachment point on the top face x top and the flow reattachment point on the leeward side of the structure x lee in the previous (Lim et al., 2009)  representation of a neutrally stable ASL. The modeling approaches such as SGS turbulent closure model and the nearwall models are the same described in Section 3.2. The dimension and velocity scales the structures are selected here to be relevant to WUI. The firebrands were released every second from positions with coordinates randomly selected with a uniform distribution from a horizontal plane passing z 35 m, as shown in Figure 7, after the flow reached a statistically stationary state. At the release points, firebrands had a zero velocity with the orientation of 60 + with respect to the vertical axis and the initial firebrand temperature T p 773 K. The initial firebrand mass density was 570 kg/m 3 , and the firebrand diameter and length of 3 mm and 40 mm, respectively (Manzello et al., 2007). Considering the flow and firebrand release conditions, the simulations here will be relevant to long range spotting. The random initial distribution of firebrands is to account for the uncertainty of the firebrand release point.

Firebrand Deposition in the Flow Over a Single Structure
To quantify the spatial distribution of the firebrands deposited on the ground and the top face of the block, a criterion proposed by Anand et al. (2018) with the following function, was used: where n is the total number of the deposited firebrands, and κ(·, ·) is the kernel function with f satisfying the normalization condition f dxdy 1. Here, B is the bandwidth, which set to 0.25 m in this study, and x i and y i are the landing coordinates of the firebrand number i. In the simulation, n ∼ 3.8 × 10 6 firebrands were deposited. The reason for release of many firebrands is to generate enough samples for the statistical description of the deposition location of firebrands. A Gaussian function was selected as the kernel function MathWorks (2019a) here. It is noted that f (x, y)   Figure 8 shows the mean velocity streamlines superimposed on the contour plots of mean velocity magnitude on the slice y 0 for varying structure sizes (panels B-H) and no structure (panel A). The streamline features here, when there is a structure, overall resemble the ones seen in Figure 3A, which is for a low Reynolds number. However, the details of these features are different for various displayed cases. In the group of structures (panels C, E, F) with fixed lengths and heights but varying widths, the horseshoe vortex and the length of the wake on the leeward side of the structure increases in size with the increase of the width. It can also be seen that the flow accelerates above the leading edge of the structure. This acceleration in the flow is more prominent for a group of structures (panels B,C,D) with fixed lengths and widths but varying heights, as the height of the structure increases. The length of the wake on the leeward side of the structure decreases slightly as the length of the structure increases as seen in the group of structures (panels C,G,H). Figure 9 shows a top view of the contour plots of NND of the deposited firebrands for cases displayed in Figure 8. As seen in this figure, there is a region of very low NND on the leeward side in panels with structures. Examining the scattered deposited particle data revealed that no firebrands were deposited on this region. This region is hereby referred to as the safe zone. The safe zone is approximately shaped like a rectangle with a  length l r and a width w r (in the spanwise direction), which is almost identical to the width of the structure W. The length l r is calculated as the distance from the leeward face of the structure to where the NND is 3.85 × 10 −4 . As seen in Figure 9, the safe zone length is larger for the structures with larger heights. Figure 10 displays l r vs. H and indicates that for every 3 m increase of the structure height, the safe zone length increases roughly by one meter. The change in width W or length L of the structure barely affected the length of the safe zone. The simulation of the structure size L W H 3 m was repeated with a grid size twice larger in each direction and it was found that l r decreased less than 6%. Figures 11A-C shows the NND of deposited firebrands vs. x at y 0 and Figures 11E,F plots it against y at x 0 for various structure sizes and the case with no structure. Seen in Figures  11A-C, are distinct troughs in cases with a structure, which correspond to the safe zones. It is also seen in these panels that NND overall decreases from the leading to the trailing edge on top of the structures. This feature is associated with the flow separation that occurs on top faces of the block, which is visible in Figure 8. It is seen in Figures 8G,H, which are for the blocks with longer lengths, this separated flow reattaches. It is believed that this reattachment gives rise to the local peaks of NND on the top face of the structure which are more pronounced for L 6 and 9 m in Figure 11C. This could be a result of some firebrands gaining momentum from the accelerated flow above the leading edge of the structure ( Figure 8) and depositing closer to its trailing edge. The curves of the cases with structures in Figures 11D-F show that the NND on top faces overall has smaller values compared to the neighboring areas on the ground. Figures 11A,D shows that an increase in the height of the structure results in a slightly higher NND on the top face of the structure. Table 3 shows the number of firebrands deposited and their temperatures on the top, front and lateral faces of the structure. In none of the cases, a firebrand was deposited on the back face of the structure. This table shows that in the cases with varying height but the same width and length, the number of firebrands deposited on the top face and their average temperature increase with an increasing height. The reason for this average temperature increase is that overall as firebrands descend, their temperatures drop. Figure 12 shows the exact location and temperature of each deposited firebrand on all faces of one of the considered structures but the leeward face. As noted earlier, the leeward face did not receive any firebrands in any of the cases. As evident in this figure, the temperature of the firebrands deposited on the windward face decrease with the decrease of z. A triangular like region with no firebrands on either lateral face of the block is noticeable.

SUMMARY AND CONCLUSION
A model was developed for simulation of cylindrical firebrand motion and burning in the FDS computational framework. The model was validated against the previous experimental and computational data (Oliveira et al., 2014) for a firebrand falling in a no-wind condition. The current model showed better agreement with the experimental data than the previous computational model. In addition, the previous experimental and CFD data (Lim et al., 2009) for a flow over a mounted 0.08 m height cube in a wind tunnel was used to validate FDS for simulation of flows over obstacles. The pressure coefficients in the simulation was in relatively good agreement with the experimental data. The mean velocity profiles in the streamwise  and vertical directions as well as the Reynolds shear stress in the simulation closely matched the experimental data. On the other hand, the simulation substantially over-predicted the measured rms of the velocities in the streamwise and vertical directions. The developed firebrand model then used with FDS to simulate the deposition of firebrands carried by a flow over a rectangular cubic structure, as a representative of a single structure in an open domain. The Reynolds number in the deposition study was an order of magnitude larger than that in the validation study. A parametric study was conducted where heights, widths and lengths were varied from 3 to 9 m. It revealed an area on the leeward side of the structure on the ground where no firebrands  were deposited. This area was refereed to as the safe zone. The width of the zone was the same as the width of the structure (the dimension of the structure in the spanwise direction). The length of this zone in the streamwise direction was proportional to the height of the structure. No firebrand was deposited on the leeward face of the structure regardless of the size of the structure. The NND on the top face of the structure increases slightly with its height. For structures with longer lengths, the NND dropped near the leading edge and rose back again toward the trailing edge of the structure. This effect was attributed to the flow accelerating above the leading edge of the structure thus imparting extra momentum onto the firebrands and carrying them farther away. Shapes of the structures considered here were simple but fundamental. Understanding the problem in fundamental setups seems an essential first step but considerations should be given to shapes representing more realistic structures. Realistic structures can significantly change from one to another in shape while involving additional geometric parameters, which can hinder the interpretation of the results. It is noted that the dimensions chosen for the structures here ranged from 3 to 9 m which are relevant to the overall dimensions of realistic small structures, e.g., houses. Future work should include sensitivity studies of the wind speed and direction. It should also include heat flux transferred from the deposited firebrands because of its consequence on ignition of the recipient fuel. Calculations of this flux require additional models to represent this phenomenon.

DATA AVAILABILITY STATEMENT
The data that support the findings of this study are available from the corresponding author upon reasonable request.

AUTHOR CONTRIBUTIONS
AM made changes to the code, ran the simulation and post processed the acquired data for this study under the supervision of BS. Both authors contributed to the development of the manuscript.

FUNDING
This work was performed under the following financial assistance award 70NANB17H281 from United States Department of Commerce, National Institute of Standards and Technology.