Numerical study of a crown fire spreading toward a fuel break using a multiphase physical model

2005 ◽  
Vol 14 (2) ◽  
pp. 141 ◽  
Author(s):  
Jean-Luc Dupuy ◽  
Dominique Morvan
Keyword(s):  
Physical Model ◽  
Numerical Solutions ◽  
Numerical Study ◽  
Vertical Plane ◽  
Solid Particles ◽  
Crown Fire ◽  
Fire Behaviour ◽  
Canopy Layer ◽  
Mediterranean Pine ◽  
Fire Spreading

The propagation of a wildfire through a Mediterranean pine stand was simulated using a multiphase physical model of fire behaviour. The heterogeneous character of the vegetation was taken into account using families of solid particles, i.e. the solid phases (foliage, twigs, grass). The thermal decomposition of the solid fuel by drying and pyrolysis, and the combustion of chars were considered, as well as the radiative and convective heat transfer between the gas and the vegetation. In the gaseous phase, turbulence was modelled using a two transport equations model (RNG k–ϵ) and the rate of combustion, which was assumed to be controlled by the turbulent mixing of fuel and oxygen, was calculated using an eddy dissipation concept. The radiation transfer equation, which includes absorption and emission of both the gas–soot mixture and the vegetation, was solved to calculate the contribution of radiation to the energy balance equations. Numerical solutions were calculated in a two-dimensional domain (vertical plane). Results showed the ability of this approach to simulate the propagation of a crown fire and to test the efficiency of a fuel break with success. The effects of the terrain slope were also tested. Some effects on fire behaviour of vortices resulting from the interaction of the wind flow with the canopy layer are shown.

Variation in wind and crown fire behaviour in a northern jack pine – black spruce forest

2004 ◽  
Vol 34 (8) ◽  
pp. 1561-1576 ◽  
Author(s):  
S W Taylor ◽  
B M Wotton ◽  
M E Alexander ◽  
G N Dalrymple
Keyword(s):  
Black Spruce ◽  
Flame Temperature ◽  
Fire Spread ◽  
Crown Fire ◽  
Fire Behaviour ◽  
Canopy Layer ◽  
Crown Fires ◽  
Video Observations ◽  
Before And After ◽  
Fire Front

Fire spread and flame temperature were examined in a series of nine experimental crown fires conducted in the Northwest Territories, Canada. Average rates of spread were 17.8–66.8 m·min–1 (0.3–1.1 m·s–1) over burning periods from about 1.5–10 min across 75 m × 75 m to 150 m × 150 m plots. Detailed maps of fire front progression revealed areas with higher rates of spread in the order of tens of metres in horizontal dimension and tens of seconds in duration in several of the fires, which is consistent with the influence of coherent wind gusts. Comparison of open and in-stand wind speed before and after burning suggests that defoliation in the canopy layer during burning would result in the flaming zone having greater exposure to the ambient wind. Estimates of flame front residence from video observations at the surface averaged 34 s; estimates from temperature measurements decreased significantly with height from 74 s at the surface to 31 s below the canopy.

A numerical study on combined baffles quick-separation device

2021 ◽  
Vol 0 (0) ◽  
Author(s):  
Ehsan Dehdarinejad ◽  
Morteza Bayareh ◽  
Mahmud Ashrafizaadeh
Keyword(s):  
Turbulent Flows ◽  
Separation Efficiency ◽  
Numerical Study ◽  
Turbulence Models ◽  
Complete Separation ◽  
Solid Particles ◽  
Flow Rates ◽  
Separation Device ◽  
Coupling Technique ◽  
Laminar And Turbulent Flows

Abstract The transfer of particles in laminar and turbulent flows has many applications in combustion systems, biological, environmental, nanotechnology. In the present study, a Combined Baffles Quick-Separation Device (CBQSD) is simulated numerically using the Eulerian-Lagrangian method and different turbulence models of RNG k-ε, k-ω, and RSM for 1–140 μm particles. A two-way coupling technique is employed to solve the particles’ flow. The effect of inlet flow velocity, the diameter of the splitter plane, and solid particles’ flow rate on the separation efficiency of the device is examined. The results demonstrate that the RSM turbulence model provides more appropriate results compared to RNG k-ε and k-ω models. Four thousand two hundred particles with the size distribution of 1–140 µm enter the device and 3820 particles are trapped and 380 particles leave the device. The efficiency for particles with a diameter greater than 28 µm is 100%. The complete separation of 22–28 μm particles occurs for flow rates of 10–23.5 g/s, respectively. The results reveal that the separation efficiency increases by increasing the inlet velocity, the device diameter, and the diameter of the particles.

A numerical study of the interaction between unsteay free-stream disturbances and localized variations in surface geometry

1989 ◽  
Vol 209 ◽  
pp. 285-308 ◽  
Author(s):  
R. J. Bodonyi ◽  
W. J. C. Welch ◽  
P. W. Duck ◽  
M. Tadjfar
Keyword(s):  
Numerical Solutions ◽  
Free Stream ◽  
Stokes Equations ◽  
Numerical Study ◽  
Navier Stokes ◽  
Surface Geometry ◽  
Navier Stokes Equations ◽  
S Waves ◽  
Tollmien Schlichting

A numerical study of the generation of Tollmien-Schlichting (T–S) waves due to the interaction between a small free-stream disturbance and a small localized variation of the surface geometry has been carried out using both finite–difference and spectral methods. The nonlinear steady flow is of the viscous–inviscid interactive type while the unsteady disturbed flow is assumed to be governed by the Navier–Stokes equations linearized about this flow. Numerical solutions illustrate the growth or decay of the T–S waves generated by the interaction between the free-stream disturbance and the surface distortion, depending on the value of the scaled Strouhal number. An important result of this receptivity problem is the numerical determination of the amplitude of the T–S waves.

A Numerical and Experimental Investigation of Turbulent Heat Transport Downstream From an Abrupt Pipe Expansion

1983 ◽  
Vol 105 (4) ◽  
pp. 862-869 ◽  
Author(s):  
R. S. Amano ◽  
M. K. Jensen ◽  
P. Goel
Keyword(s):  
Heat Transfer ◽  
Numerical Solutions ◽  
Numerical Study ◽  
Separated Flow ◽  
Flow Region ◽  
Reynolds Numbers ◽  
Turbulent Heat ◽  
Transfer Coefficients ◽  
Heat Transfer Coefficients ◽  
Pipe Expansion

An experimental and numerical study is reported on heat transfer in the separated flow region created by an abrupt circular pipe expansion. Heat transfer coefficients were measured along the pipe wall downstream from an expansion for three different expansion ratios of d/D = 0.195, 0.391, and 0.586 for Reynolds numbers ranging from 104 to 1.5 × 105. The results are compared with the numerical solutions obtained with the k ∼ ε turbulence model. In this computation a new finite difference scheme is developed which shows several advantages over the ordinary hybrid scheme. The study also covers the derivation of a new wall function model. Generally good agreement between the measured and the computed results is shown.

Computational Fluid Dynamics Evaluation of Heat Transfer Correlations for Sodium Flows in a Heat Exchanger

2010 ◽  
Vol 132 (5) ◽  
Author(s):  
Seok-Ki Choi ◽  
Seong-O Kim ◽  
Hoon-Ki Choi
Keyword(s):  
Heat Transfer ◽  
Heat Exchanger ◽  
Turbulence Model ◽  
Numerical Solutions ◽  
Numerical Study ◽  
Cross Flow ◽  
Turbulence Models ◽  
Parallel Flow ◽  
Sst Model ◽  
Heat Transfer Correlations

A numerical study for the evaluation of heat transfer correlations for sodium flows in a heat exchanger of a fast breeder nuclear reactor is performed. Three different types of flows such as parallel flow, cross flow, and two inclined flows are considered. Calculations are performed for these three typical flows in a heat exchanger changing turbulence models. The tested turbulence models are the shear stress transport (SST) model and the SSG-Reynolds stress turbulence model by Speziale, Sarkar, and Gaski (1991, “Modelling the Pressure-Strain Correlation of Turbulence: An Invariant Dynamical System Approach,” J. Fluid Mech., 227, pp. 245–272). The computational model for parallel flow is a flow past tubes inside a circular cylinder and those for the cross flow and inclined flows are flows past the perpendicular and inclined tube banks enclosed by a rectangular duct. The computational results show that the SST model produces the most reliable results that can distinguish the best heat transfer correlation from other correlations for the three different flows. It was also shown that the SSG-RSTM high-Reynolds number turbulence model does not deal with the low-Prandtl number effect properly when the Peclet number is small. According to the present calculations for a parallel flow, all the old correlations do not match with the present numerical solutions and a new correlation is proposed. The correlations by Dwyer (1966, “Recent Developments in Liquid-Metal Heat Transfer,” At. Energy Rev., 4, pp. 3–92) for a cross flow and its modified correlation that takes into account of flow inclination for inclined flows work best and are accurate enough to be used for the design of the heat exchanger.

scholarly journals Assymmetric eddy populations in adjacent basins – a high resolution numerical study of the Tyrrhenian and Ligurian Seas

2012 ◽  
Vol 9 (6) ◽  
pp. 3521-3566 ◽  
Author(s):  
R. M. A. Caldeira ◽  
X. Couvelard ◽  
E. Casella ◽  
A. Vetrano
Keyword(s):  
High Resolution ◽  
Ocean Circulation ◽  
Numerical Solutions ◽  
Numerical Study ◽  
Mesoscale Eddy ◽  
Leeward Side ◽  
Boundary Current ◽  
High Concentration ◽  
Eddy Activity ◽  
Tyrrhenian Basin

Abstract. A high-resolution ocean circulation modelling system forced with a high-resolution numerical wind product was used to study the mesoscale and sub-mesoscale eddy population of the North-Western Mediterranean Sea, contrasting eddy-activity between the Tyrrhenian and Ligurian sub-basins. Numerical solutions reproduced some of the known regional dynamics, namely the occurrence and oceanic implications of Mistral events, the convective cell leeward of the Gulf of Lion, as well as the Balearic frontal system. Calculated transport across the Corsica Channel followed a similar trend, when compared to the transport computed from a moored current meter. The analysis of the results showed that surface eddy activity is mostly confined to the boundary-currents, whereas in the deeper layers most eddies are concentrated on the central-deeper part of the basins. The Liguro-Provençal basin shows a much higher concentration of intermediate and deep-water eddies, when compared to the Tyrrhenian basin. Sub-mesoscale surface eddies tend to merge and migrate vertically onto intermediate waters. Intense eddy activity in the boundary-current surrounding the Liguro-Provençal Gyre, concentrate high-productivity, manifested by higher concentrations of mean sea surface chlorophyll, in the central part of the gyre, defined herein as the Ligurian Productive Pool (LPP). On average, the Tyrrhenian was mostly oligotrophic except for a small productive vortice in the south-eastern (leeward) side of Corsica. The transport in the Tyrrhenian Gyre, and across the basin is one order of magnitude higher than the transport calculated for the Liguro-Provençal basin. A high concentration of eddies in the passage between the Balearic Archipelago and Sardinia suggests retention and longer residence times of nutrient rich water in the "Ligurian pool", compared to a "fast draining" Tyrrhenian basin. Previous studies support the cyclonic gyre circulation generated in the Liguro-Provençal basin but more studies are needed to address the surface and deep mesoscale activity of the Tyrrhenian basin.

Numerical Solution for Boundary Layer Flow of a Dusty Micropolar Fluid Due to a Stretching Sheet with Constant Wall Temperature

2021 ◽  
Vol 87 (1) ◽  
pp. 30-40
Author(s):  
Anisah Dasman ◽  
Abdul Rahman Mohd Kasim ◽  
Iskandar Waini ◽  
Najiyah Safwa Khashi’ie
Keyword(s):  
Differential Equations ◽  
Ordinary Differential Equations ◽  
Wall Temperature ◽  
Micropolar Fluid ◽  
Stretching Sheet ◽  
Numerical Solutions ◽  
Numerical Study ◽  
Particle Interaction ◽  
Dust Particles ◽  
Constant Wall Temperature

This paper aims to present the numerical study of a dusty micropolar fluid due to a stretching sheet with constant wall temperature. Using the suitable similarity transformation, the governing partial differential equations for two-phase flows of the fluid and the dust particles are reduced to the form of ordinary differential equations. The ordinary differential equations are then numerically analysed using the bvp4c function in the Matlab software. The validity of present numerical results was checked by comparing them with the previous study. The results graphically show the numerical solutions of velocity, temperature and microrotation distributions for several values of the material parameter K, fluid-particle interaction parameter and Prandtl number for both fluid and dust phase. The effect of microrotation is investigated and analysed as well. It is found that the distributions are significantly influenced by the investigated parameters for both phases.

Implementation of DTM as a numerical study for the Casson fluid flow past an exponentially variable stretching sheet with thermal radiation

2021 ◽  
pp. 2150111
Author(s):  
Khadijah M. Abualnaja
Keyword(s):  
Fluid Flow ◽  
Thermal Radiation ◽  
Numerical Method ◽  
Stretching Sheet ◽  
Numerical Solutions ◽  
Numerical Study ◽  
Transformation Method ◽  
Casson Fluid ◽  
Physical Parameters ◽  
Special Cases

This paper introduces a theoretical and numerical study for the problem of Casson fluid flow and heat transfer over an exponentially variable stretching sheet. Our contribution in this work can be observed in the presence of thermal radiation and the assumption of dependence of the fluid thermal conductivity on the heat. This physical problem is governed by a system of ordinary differential equations (ODEs), which is solved numerically by using the differential transformation method (DTM). This numerical method enables us to plot figures of the velocity and temperature distribution through the boundary layer region for different physical parameters. Apart from numerical solutions with the DTM, solutions to our proposed problem are also connected with studying the skin-friction coefficient. Estimates for the local Nusselt number are studied as well. The comparison of our numerical method with previously published results on similar special cases shows excellent agreement.

A numerical study of the nucleation, growth and settling of crystals from a turbulent convecting fluid

2021 ◽  
Author(s):  
Vojtech Patocka ◽  
Nicola Tosi ◽  
Enrico Calzavarini
Keyword(s):  
Magma Chamber ◽  
Solid Fraction ◽  
Large Scale ◽  
Solid Phase ◽  
Numerical Study ◽  
Equilibrium Concentration ◽  
Solid Particles ◽  
Settling Rate ◽  
Carrier Fluid ◽  
Nucleation Growth

<p>We evaluate the equilibrium concentration of a thermally convecting suspension that is cooled from above and in which<br>solid crystals are self-consistently generated in the thermal boundary layer near the top. In a previous study (Patočka et<br>al., 2020), we investigated the settling rate of solid particles suspended in a highly vigorous (Ra = 10<sup>8</sup> , 10<sup>10</sup>, and 10<sup>12</sup> ),<br>finite Prandtl number (Pr = 10, 50) convection. In this follow-up study we additionally employ the model of crystal<br>generation and growth of Jarvis and Woods (1994), instead of using particles with a predefined size and density that are<br>uniformly injected into the carrier fluid.</p><p>We perform a series of numerical experiments of particle-laden thermal convection in 2D and 3D Cartesian geometry<br>using the freely available code CH4 (Calzavarini, 2019). Starting from a purely liquid phase, the solid fraction gradually<br>grows until an equilibrium is reached in which the generation of the solid phase balances the loss of crystals due to<br>sedimentation at the bottom of the fluid. For a range of predefined density contrasts of the solid phase with respect to<br>the density of the fluid (ρ<sub>p</sub> /ρ<sub>f</sub> = [0, 2]), we measure the time it takes to reach such equilibrium. Both this time and<br>the equilibrium concentration depend on the average settling rate of the particles and are thus non-trival to compute for<br>particle types that interact with the large-scale circulation of the fluid (see Patočka et al., 2020).</p><p>We apply our results to the cooling of a large volume of magma, spanning from a large magma chamber up to a<br>global magma ocean. Preliminary results indicate that, as long as particle re-entrainment is not a dominant process, the<br>separation of crystals from the fluid is an efficient process. Fractional crystallization is thus expected and the suspended<br>solid fraction is typically small, prohibiting phenomena in which the feedback of crystals on the fluid begins to govern the<br>physics of the system (e.g. Sparks et al, 1993).</p><p>References<br>Patočka V., Calzavarini E., and Tosi N.(2020). Settling of inertial particles in turbulent Rayleigh-Bénard convection.<br>Physical Review Fluids, 26(4) 883-889.</p><p>Jarvis, R. A. and Woods, A. W.(1994). The nucleation, growth and settling of crystals from a turbulently convecting<br>fluid. J. Fluid. Mech, 273 83-107.</p><p>Sparks, R., Huppert, H., Koyaguchi, T. et al (1993). Origin of modal and rhythmic igneous layering by sedimentation in<br>a convecting magma chamber. Nature, 361, 246-249.</p><p>Calzavarini, E (2019). Eulerian–Lagrangian fluid dynamics platform: The ch4-project. Software Impacts, 1, 100002.</p>

scholarly journals Numerical Study of Branched Turboprop Inlet Ducts Using a Multiple Block Grid Procedure

1991 ◽  
Author(s):  
Anil K. Tolpadi ◽  
Mark E. Braaten
Keyword(s):  
Numerical Solutions ◽  
Numerical Study ◽  
Three Dimensional ◽  
Separated Flow ◽  
Experimental Tests ◽  
Spherical Particles ◽  
Branch Duct ◽  
Pressure Losses ◽  
Overlapping Grids ◽  
Inlet Duct

An important requirement in the design of an inlet duct of a turboprop engine is the ability to provide foreign object damage protection. A possible method for providing this protection is to include a bypass branch duct as an integral part of the main inlet duct. This arrangement would divert ingested debris away from the engine through the bypass. However, such an arrangement could raise the possibility of separated flow in the inlet, which in turn can increase pressure losses if not properly accounted for during the design. A fully elliptic three-dimensional body-fitted computational fluid dynamics (CFD) code based on pressure correction techniques has been developed that has the capability of performing multiple block grid calculations compatible with present day turboshaft and turboprop branched inlet ducts. Calculations are iteratively performed between sets of overlapping grids with one grid representing the main duct and a second grid representing the branch duct. Both the grid generator and the flow solver have been suitably developed to achieve this capability. The code can handle multiple branches in the flow. Using the converged flow field from this code, another program was written to perform a particle trajectory analysis. Numerical solutions were obtained on a supercomputer for a typical branched duct for which experimental flow and pressure measurements were also made. The flow separation zones predicted by the calculations were found to be in good agreement with those observed in the experimental tests. The total pressure recovery factors measured in the experiments were also compared with those obtained numerically. Within the limits of the grid resolution and the turbulence model, the agreement was found to be fairly good. In order to simulate the path of debris entering the duct, the trajectories of spherical particles of different sizes introduced at the inlet were determined.

Sign in / Sign up

Export Citation Format

Share Document