Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “RADIATION FLUX”

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 109 records · Page 6

Small-Scale Variability of Solar Radiation (S2VSR) Field Campaign Report

Clouds introduce significant variability in solar radiation, in particular for shallow cumulus and other convective cloud types. Variability at small spatio-temporal scales (specifically below 1 min and 1 km) is, however, not properly resolved by current satellite observations nor atmospheric models. Additionally, 3D photon transport is commonly neglected in analyses and causes significant deviations in 1D radiative transfer simulations at such small scales. A specific challenge is to account for the point-like nature of traditional surface radiative flux measurements in comparison to spatially averaged quantities such as satellite pixel radiances or fluxes representing model grid boxes.

54 ENVIRONMENTAL SCIENCES↗

The Preservation of Super-Earths and the Emergence of Gas Giants after Their Progenitor Cores Have Entered the Pebble-isolation Phase

The omnipresence of super-Earths suggests that they are able to be retained in natal disks around low-mass stars, whereas exoplanets' mass distributions indicate that some cores have transformed into gas giants through runaway gas accretion at ≳1 au from solar-type stars. In this paper, we show that transition to runaway gas accretion by cores may be self-impeded by an increase of the grain opacity in their envelope after they have acquired sufficient mass (typically ~ 10M ⊕ ) to enter a pebble-isolation phase. The accumulation of approximately millimeter- to meter-sized pebbles in their migration barriers enhances their local fragmentation rates. The freshly produced submillimeter grains pass through the barrier, elevate the effective dust opacity, and reduce the radiative flux in the core's envelope. These effects alone are adequate to suppress the transition to runaway accretion and preserve super-Earths in the stellar proximity (~0.1 au), albeit entropy advection between the envelope and the disk can further reduce the accretion rate. At an intermediate distance (~1 au) from their host stars, the escalation in the dust opacity dominates over entropy advection in stalling the transition to runaway accretion for marginally pebble-isolated cores. Beyond a few astronomical unit, the transformation of more-massive cores to gas giants is reachable before severe depletion of disk gas. This requirement can be satisfied either in extended disks with large scale height via orderly accretion of migrating pebbles or through the mergers of oligarchic protoplanetary embryos, and can account for the correlated occurrence of long-period gas giants and close-in super-Earths.

79 ASTRONOMY AND ASTROPHYSICS↗

Aerosols from anthropogenic and biogenic sources and their interactions – modeling aerosol formation, optical properties, and impacts over the central Amazon basin

The Green Ocean Amazon experiment – GoAmazon 2014–2015 – explored the interactions between natural biogenic forest emissions from central Amazonia and urban air pollution from Manaus. Previous GoAmazon 2014–2015 studies showed that nitrogen oxide (NO x = NO + NO 2 ) and sulfur oxide (SO x ) emissions from Manaus strongly interact with biogenic volatile organic compounds (BVOCs), affecting secondary organic aerosol (SOA) formation. In previous studies, ground-based and aircraft measurements provided evidence of SOA formation and strong changes in aerosol composition and properties. Aerosol optical properties also evolve, and their impacts on the Amazonian ecosystem can be significant. As particles age, some processes, such as SOA production, black carbon (BC) deposition, particle growth and the BC lensing effect change the aerosol optical properties, affecting the solar radiation flux at the surface. This study analyzes data and models SOA formation using the Weather Research and Forecasting with Chemistry (WRF-Chem) model to assess the spatial variability in aerosol optical properties as the Manaus plumes interact with the natural atmosphere. The following aerosol optical properties are investigated: single scattering albedo (SSA), asymmetry parameter (g aer ), absorption Ångström exponent (AAE) and scattering Ångström exponent (SAE). These simulations were validated using ground-based measurements at three experimental sites, namely the Amazon Tall Tower Observatory – ATTO (T0a), downtown Manaus (T1), Tiwa Hotel (T2) and Manacapuru (T3), as well as the U.S. Department of Energy (DOE) Gulfstream 1 (G-1) aircraft flights. WRF-Chem simulations were performed over 7 d during March 2014. Results show a mean biogenic SOA (BSOA) mass enrichment of 512 % at the T1 site, 450 % in regions downwind of Manaus, such as the T3 site, and 850 % in areas north of the T3 site in simulations with anthropogenic emissions. The SOA formation is rather fast, with about 80 % of the SOA mass produced in 3–4 h. Comparing the plume from simulations with and without anthropogenic emissions, SSA shows a downwind reduction of approximately 10 %, 11 % and 6 % at the T1, T2 and T3 sites, respectively. Other regions, such as those further downwind of the T3 site, are also affected. The g aer values increased from 0.62 to 0.74 at the T1 site and from 0.67 to 0.72 at the T3 site when anthropogenic emissions are active. During the Manaus plume-aging process, a plume tracking analysis shows an increase in SSA from 0.91 close to Manaus to 0.98 160 km downwind of Manaus as a result of SOA production and BC deposition.

54 ENVIRONMENTAL SCIENCES↗

Relative importance of high-latitude local and long-range-transported dust for Arctic ice-nucleating particles and impacts on Arctic mixed-phase clouds

Dust particles, serving as ice-nucleating particles (INPs), may impact the Arctic surface energy budget and regional climate by modulating the mixed-phase cloud properties and lifetime. In addition to long-range transport from low-latitude deserts, dust particles in the Arctic can originate from local sources. However, the importance of high-latitude dust (HLD) as a source of Arctic INPs (compared to low-latitude dust, LLD) and its effects on Arctic mixed-phase clouds are overlooked. In this study, we evaluate the contribution to Arctic dust loading and INP population from HLD and six LLD source regions by implementing a source-tagging technique for dust aerosols in version 1 of the US Department of Energy's Energy Exascale Earth System Model (E3SMv1). Our results show that HLD is responsible for 30.7 % of the total dust burden in the Arctic, whereas LLD from Asia and North Africa contributes 44.2 % and 24.2 %, respectively. Due to its limited vertical transport as a result of stable boundary layers, HLD contributes more in the lower troposphere, especially in boreal summer and autumn when the HLD emissions are stronger. LLD from North Africa and East Asia dominates the dust loading in the upper troposphere with peak contributions in boreal spring and winter. The modeled INP concentrations show better agreement with both ground and aircraft INP measurements in the Arctic when including HLD INPs. The HLD INPs are found to induce a net cooling effect (-0.24 W m -2 above 60 °N) on the Arctic surface downwelling radiative flux by changing the cloud phase of the Arctic mixed-phase clouds. The magnitude of this cooling is larger than that induced by North African and East Asian dust (0.08 and -0.06 W m -2 , respectively), mainly due to different seasonalities of HLD and LLD. Uncertainties of this study are discussed, which highlights the importance of further constraining the HLD emissions.

54 ENVIRONMENTAL SCIENCES↗

Size-resolved process understanding of stratospheric sulfate aerosol following the Pinatubo eruption

Stratospheric sulfate aerosol produced by volcanic eruptions plays important roles in atmospheric chemistry and the global radiative balance of the atmosphere. The simulation of stratospheric sulfate concentrations and optical properties is highly dependent on the chemistry scheme and microphysical treatment. In this work, we implemented a sophisticated gas-phase chemistry scheme (full chemistry, FC) and a 5-mode version of the Modal Aerosol Module with Prognostic Stratospheric Aerosol (MAM5-PSA) for the interactive treatment of stratospheric sulfate aerosol in the Department of Energy's Energy Exascale Earth System Model version 2 (E3SMv2) model to better simulate the chemistry-aerosol feedback following the Pinatubo eruption, and to compare it against a simulation using simplified chemistry (SC) and the default 4-mode version of the Modal Aerosol Module (MAM4). MAM5-PSA experiments were found to better capture the stratospheric sulfate burden from the eruption of the volcano to the end of 1992 as compared to the High-resolution Infrared Sounder (HIRS) observations, and the formation of sulfate in MAM5-PSA with FC (with an additional OH replenishment reaction) was significantly faster than in MAM4 with FC. Analyses of microphysical processes indicate that more sulfate aerosol mass was generated in total in FC experiments than in SC experiments. MAM5-PSA performs better than MAM4 in simulation of aerosol optical depth (AOD); AOD anomalies from the MAM5-PSA experiment have better agreement with observations. The simulated largest changes in global mean net radiative flux at the top of the atmosphere following the eruption were about −3 W m −2 in MAM5-PSA experiments and roughly −1.5 W m −2 in MAM4 experiments.

AEROSOL↗

Biomass burning aerosol radiative effects in the Southeast Atlantic depend strongly on meteorological forcing method

Biomass burning aerosols (BBAs) from African fires may strongly impact Earth’s radiation budget in the southeast Atlantic (SEA), but the sign and magnitude of the overall radiative effect (RE) remain uncertain. Aerosol–climate models are needed to separately quantify direct, indirect, and semi-direct REs. Here, we evaluate improved simulations with the UK Met Office's Unified Model and explore REs resulting from various methods used to match observed meteorology (nudging or running forecasts reinitialized at different frequencies). REs are calculated as differences in radiative fluxes between simulations with and without smoke emissions and with and without aerosol absorption. All model setups agree on net warming for the SEA dominated by the direct effect. Simulated smoke, clouds, and the direct effect agree better with observations than previous studies using the same model, though biases in aerosol extinction and liquid water path remain. Changes in cloud droplet number concentration due to BBA self-lofting influence how cleanly we can separate cloud effects into semi-direct and indirect effects. Total RE, which remains unaffected, ranges from +3.0 to +7.9 W m −2 . The 4.9 W m −2 spread arises mainly from simulated semi-direct effects. Forecasts three days long or less probably do not allow time for plausible differences in boundary layer properties due to semi-direct effects to accumulate. Free running simulations with and without smoke accumulate differences in meteorology that are likely spurious “butterfly effects”. We recommend future research quantifying BBA REs over weeks to months to use meteorological forcing techniques that allow aerosol absorption to affect the boundary layer.

Giuffrida, Eric [Carnegie Mellon University, Pitts↗

ClimGen: Learning the Forcing-Response Relationship in Climate System

Solar Radiation Management (SRM) is emerging as a potential geoengineering strategy to address the anthropogenic impact on climate, but its effective implementation requires an iterative and large ensemble of highly accurate and efficient climate projections. Traditional climate projections rely on executing computationally demanding and time-consuming numerical climate models. Recent advances in machine learning (ML) aim to enhance these approaches by emulating traditional methods. In this work, we propose a novel framework for directly learning the relationship between solar radiation flux at the top of the atmosphere and the corresponding surface temperature response. To evaluate the feasibility of this direct ML-based projection, we developed a dataset using an intermediate complexity model, incorporating a comprehensive suite of different forcing patterns and evaluation metrics to rigorously assess the ML model’s performance. We introduce a Conditional Denoising Diffusion Probabilistic Model (cDDPM) for this task, which demonstrates encouraging skill in representing climate statistics under previously unseen forcing patterns. This approach provides a promising pathway for direct climate projections by accurately learning the forcing-response relationship, with a wide range of applications in impact mitigation, emissions policy design, and SRM strategies.

Chen, Tse-Chun [BATTELLE (PACIFIC NW LAB)] (ORCID:↗

A network of soil moisture, soil temperature, air temperature, net radiation, ground heat flux and ground water for Chicago, Illinois

This dataset contains environmental monitoring data collected using solar-powered Multi-Function Research (MFR) Long Range Wide Area (LoRaWAN)-enabled nodes at 11 sites in Chicago, Illinois, as part of the DOE Urban Integrated Field Lab CROCUS project. The MFR node system consists of an Input/Output Digital Input Module (IB8) interface box (ICT International) providing wired connections for environmental sensors and an MFR-Node-L data logger that manages power, data processing, and LoRaWAN communication. The wireless data are ingested via Sage network (https://sagecontinuum.org/) nodes that contain LoRaWAN antennae. Measurements were collected from 11 MFR nodes deployed across Chicago State University (CSU), Northeastern Illinois University (NEIU), Northwestern University (NU), University of Illinois Chicago (UIC), West Woodlawn "Blacks in Green" (BIG), and Indian Boundary Prairies (IBP). Each MFR node supports a consistent suite of sensors measuring atmospheric, soil, and hydrological variables. Atmospheric measurements include 2m air temperature (°C), 2m vapor pressure deficit (kPa), and 2m shortwave/longwave radiation (incoming and outgoing, W/m²) measured using ATH-VPD and Apogee SN500 sensors. Soil measurements include volumetric water content (VWC, %) and temperature (°C) at four depths (15, 30, 45, and 60 cm below surface) using Meter Teros54 sensors, and heat flux (W/m²) at 10 cm depth using Huske HFP01-05 sensors. At selected locations, Meter Hydros21 sensors measure groundwater depth (mm), specific conductivity (dS/m), and temperature (°C). The dataset includes timestamps, site identifiers with location names, device IDs, Global Positioning System (GPS) coordinates, variable names with units, measurement depths, values, sensor names, and Sage node identifiers. All timestamps are in local Chicago time (CDT/CST). Quality control flags are provided using a 6-bit binary system indicating physical range violations, step spikes, 24-hour flat-line conditions, 6-hour jitter, 7-day ultra-low variance, and persistent high offset. Data is provided in CSV and CF-compliant NetCDF formats. This dataset is part of a larger collection of CROCUS environmental monitoring data, including linked datasets from Air Quality Transmitter (AQT) sensors, Weather Transmitter (WXT) sensors, and Sap Flow Meter (SFM1x) sensors.

Chicago↗

A More La Niña–Like Response to Radiative Forcing after Flux Adjustment in CESM2

Abstract In response to greenhouse gas forcing, most coupled global climate models project the tropical Pacific SST trend toward an “El Niño–like” state, with a reduced zonal SST gradient and a weakened Walker circulation. However, observations over the last five decades reveal a trend toward a more “La Niña–like” state with a strengthening zonal SST gradient. Recent research indicates that the identified trend differences are unlikely to be entirely due to internal variability and probably result, at least in part, from systematic model biases. In this study, Community Earth System Model, version 2 (CESM2), is used to explore how mean-state biases within the model may influence its forced response to radiative forcing in the tropical Pacific. The results show that using flux adjustment to reduce the mean-state bias in CESM2 over the tropical regions results in a more La Niña–like trend pattern in the tropical Pacific, with a strengthening of the tropical Pacific zonal SST gradient and a relatively enhanced Walker circulation, as hypothesized to occur if the ocean thermostat mechanism is stronger than the atmospheric mechanisms which by themselves would weaken the Walker circulation. We also find that the historical strengthening of the tropical Pacific zonal gradient is transient but persists into the near term in a high-emissions future warming scenario. These results suggest the potential of flux adjustment as a method for developing alternative projections that represent a wider range of possible future tropical Pacific warming scenarios, especially for a better understanding of regional patterns of climate risk in the near term.

Zhuo, Jing-Yi [Lamont-Doherty Earth Observatory, C↗

Over three decades, and counting, of near-surface turbulent flux measurements from the Atmospheric Radiation Measurement (ARM) user facility

Processes mediating the coupling of terrestrial, aquatic, biospheric, and atmospheric systems influence weather, climate, and ecosystem dynamics via transfer of energy, momentum, water, and carbon (or other species). These exchange processes are quantified by measurements of near-surface turbulent fluxes. Understanding processes at these interfaces provides insight toward understanding and predicting current and future states within the Earth system. The Atmospheric Radiation Measurement (ARM) user facility has been conducting measurements of near-surface turbulent fluxes since the early 1990s at long-term fixed locations and shorter-term mobile deployments across the Earth. ARM has utilized two established methods for conducting these measurements: energy balance Bowen ratio (EBBR) and eddy covariance (EC). Primary measurements from the former include sensible and latent heat flux, while the latter also measures fluxes of momentum and carbon (primarily carbon dioxide, with methane fluxes measured at two locations to date). The EBBR systems have been deployed at 22 locations, and, to date, the EC systems have been deployed at over 50 sites, with plans for additional novel site locations in the future. Herein, the history, evolution, and key aspects of these instrument systems are documented, along with information on data quality assurance and post-processing, as well as best use practices. Additionally, three data validation experiments were recently conducted, and their key findings are summarized. Finally, ancillary datasets acquired by ARM, which can contextualize and aid interpretation of the near-surface turbulent flux measurements, are discussed. The datasets described herein include the eddy correlation flux measurement system: 30ECOR (https://doi.org/10.5439/1879993, Sullivan et al., 1997), 30QCECOR (https://doi.org/10.5439/1097546, Gaustad, 2003), ECORSF (https://doi.org/10.5439/1494128, Sullivan et al., 2019a), and associated AmeriFlux and Methane Value-Added Product, AMCMETHANE (https://doi.org/10.5439/1508268, Billesbach, 2011); the energy balance Bowen ratio system: 30EBBR (https://doi.org/10.5439/1023895, Sullivan et al., 1993) and 30BAEBBR (https://doi.org/10.5439/1027268, Gaustad and Xie, 1993); and the carbon dioxide flux measurement system: CO2FLX (https://doi.org/10.5439/1287574, https://doi.org/10.5439/1287575, https://doi.org/10.5439/1287576, Koontz et al., 2015a, b, c; https://doi.org/10.5439/1989774, https://doi.org/10.5439/1989776, https://doi.org/10.5439/1992202, Biraud and Chan, 2002a, b, c). These data can be found by searching the above data stream names at https://adc.arm.gov/discovery/#/results/ (last access: 8 September 2025).

Sullivan, Ryan C. [Argonne National Laboratory (AN↗

Effect of soot concentration on radiative transport in a 2 m JP-8 pool fire

Hubbard et al. (2024) investigated conjugate heat transfer in a 7.9 m JP-8 pool fire, utilizing both heavily sooting and non-sooting fuel models. Their findings indicated that all models produced pool radiative heat fluxes within the measurement uncertainty for large-scale outdoor experiments. This study aims to further investigate the influence of soot concentration on radiative transport. To this end, Large Eddy Simulations (LES) of a 2 m diameter JP-8 pool fire were conducted, employing four turbulent combustion chemistry models: the Eddy Dissipation Concept (EDC) and three variants of flamelet model. Additionally, three variations of the one-equation subgrid scale turbulence kinetic energy model, and multiple mesh resolutions were utilized. The primary validation metrics were pool radiative heat flux and plume radiative fraction, with simulation data compared to experimental results published by Blanchat et al. (2010). The EDC model produced high soot concentrations near the pool surface, resulting in a significant shielding effect. In contrast, the flamelet combustion model has different kinetic limiters for soot production, leading to soot formation higher in the flame, where it can be more readily oxidized. Flamelet models incorporating both acetylene and benzene contributions to soot growth demonstrated better agreement with experimental data. To assess the sensitivity of pool radiative heat flux to soot mass concentration, we used a one-dimensional approximation over a wide range of soot concentrations and compared it to three-dimensional solutions obtained using Discrete Ordinates. Furthermore, this analysis illustrates the impacts of soot shielding and Turbulence Radiation Interactions (TRI), with TRI emission enhancement factors predicted to be in the range of 3–5 near the edges of the plume.

Computational fluid dynamics↗

Proper orthogonal decomposition based reduced-order modeling of flux-Limited gray thermal radiation

Here, in this work, a proper orthogonal decomposition (POD) based reduced-order model (ROM) is developed to solve gray, flux-limited thermal radiation diffusion. We focus on the variable opacity radiation penetration benchmark posed by Olson, Auer, and Hall. The T -3 relationship for opacity in conjunction with high-temperature radiation penetrating an initially cold material produces a strong thermal radiation shock. This class of problems is particularly challenging for standard POD-based reduced-order modeling due to the nonlinearities presented by 1) the T 4 source term, and 2) flux-limited diffusion operator. To address these challenges and develop a cost competitive ROM, we employ a “hyper-reduction” technique through discrete empirical interpolation (DEIM) and allow for adaptive reduced-order projections through principal interval decomposition (PID). Performance of the proposed methodology is quantified by comparing the cost savings and accuracy relative to a full-order computation. Reference solutions and snapshot data are obtained through a full-order calculation performed by the University of Chicago maintained astrophysics code, FLASH. For consistency and potential extensibility, the developed ROM is also implemented in FLASH. We find that in the initialization regime, where the thermal radiation wave is initially created by the warming of the material, this class of problems is highly reducible and suitable for POD-based ROMs. However, the strong convective nature of the wave propagation regime is less reducible and more challenging to create an efficient ROM.

42 ENGINEERING↗

Spontaneous Cyclogenesis without Radiative and Surface-Flux Feedbacks

Tropical cyclones (TCs) are among the most intense and feared storms in the world. What physical processes lead to cyclogenesis remains the most mysterious aspect of TC physics. Here, we study spontaneous TC genesis in rotating radiative–convective equilibrium using cloud-resolving simulations over an f plane with constant sea surface temperature. Previous studies proposed that spontaneous TC genesis requires either radiative or surface-flux feedbacks. To test this hypothesis, we perform mechanism-denial experiments, in which we switch off both feedback processes in numerical simulations. We find that TCs can self-emerge even without radiative and surface-flux feedbacks. Although these feedbacks accelerate the genesis and impact the size of the TCs, TCs in the experiments without them can reach similar intensities as those in the control experiment. We show that TC genesis is associated with an increase in the available potential energy (APE) and that convective heating dominates the APE production. Overall, our result suggests that spontaneous TC genesis may result from a cooperative interaction between convection and circulation and that radiative and surface-flux feedbacks accelerate the process. Furthermore, we find that increasing the planetary rotation favors spontaneous TC genesis.

54 ENVIRONMENTAL SCIENCES↗

ARM Data-Oriented Metrics and Diagnostics Package for Climate Model Evaluation

A Python-based metrics and diagnostics package is currently being developed by the U.S. Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) Infrastructure Team at Lawrence Livermore National Laboratory (LLNL) to facilitate the use of long-term, high-frequency measurements from the ARM Facility in evaluating the regional climate simulation of clouds, radiation, and precipitation. This metrics and diagnostics package computes climatological means of targeted climate model simulation and generates tables and plots for comparing the model simulation with ARM observational data. The Coupled Model Intercomparison Project (CMIP) model data sets are also included in the package to enable model intercomparison as demonstrated in Zhang et al. (2017). The mean of the CMIP model can serve as a reference for individual models. Basic performance metrics are computed to measure the accuracy of mean state and variability of climate models. The evaluated physical quantities include cloud fraction, temperature, relative humidity, cloud liquid water path, total column water vapor, precipitation, sensible and latent heat fluxes, and radiative fluxes, with plan to extend to more fields, such as aerosol and microphysics properties. Process-oriented diagnostics focusing on individual cloud- and precipitation-related phenomena are also being developed for the evaluation and development of specific model physical parameterizations. The version 1.0 package is designed based on data collected at ARM’s Southern Great Plains (SGP) Research Facility, with the plan to extend to other ARM sites. The metrics and diagnostics package is currently built upon standard Python libraries and additional Python packages developed by DOE (such as CDMS and CDAT). The ARM metrics and diagnostic package is available publicly with the hope that it can serve as an easy entry point for climate modelers to compare their models with ARM data. In this report, we first present the input data, which constitutes the core content of the metrics and diagnostics package in section 2, and a user's guide documenting the workflow/structure of the version 1.0 codes, and including step-by-step instruction for running the package in section 3.

54 ENVIRONMENTAL SCIENCES↗

Near-Field Radiative Heat Transfer Between Two α-MoO3 Biaxial Crystals

The near-field radiative heat transfer (NFRHT) between two semi-infinite α-MoO3 biaxial crystals is investigated numerically based on the fluctuation–dissipation theorem combined with the modified 4 × 4 transfer matrix method in this paper. In the calculations, the near-field radiative heat flux (NFRHF) along each of the crystalline directions of α-MoO3 is obtained by controlling the orientation of the biaxial crystals. The results show that much larger heat flux than that between two semi-infinite hexagonal boron nitride can be achieved in the near-field regime, and the maximum heat flux is along the [001] crystalline direction. The mechanisms for the large radiative heat flux are explained as due to existence of hyperbolic phonon polaritons (HPPs) inside α-MoO3 and excitation of hyperbolic surface phonon polaritons (HSPhPs) at the vacuum/α-MoO3 interfaces. The effect of relative rotation between the emitter and the receiver on the heat flux is also investigated. It is found that the heat flux varies significantly with the relative rotation angle. The modulation contrast can be as large as two when the heat flux is along the [010] direction. We attribute the large modulation contrast mainly to the misalignment of HSPhPs and HPPs between the emitter and the receiver. Hence, the results obtained in this work may provide a promising way for manipulating near-field radiative heat transfer between anisotropic materials.

Engineering↗

Validation of the HyRAM+ physics models for use with propane

The Hydrogen Risk Assessment Model Plus (HyRAM+) toolkit combines quantitative risk assessment with simulations of unignited dispersion, ignited turbulent diffusion flames, and indoor accumulation with delayed ignition of fuels. HyRAM+ is differentiated from HyRAM in that it includes models and leak data for other alternate fuels. The models of the physical phenomena need to be validated for each of the fuels in the toolkit. This report shows the validation for propane which is being used as a surrogate for autogas, which is a mixture of propane and butane and used in internal combustion engines in vehicles. For flame length comparisons, five previously published experiments from peer reviewed journals were used to validate our models. The validation looked at flame lengths and flame widths with respect to different leak diameters, mass flow rates, and source pressures. Most of the sources included more than one set of experimental data, which were collected using different methods (CCD cameras, IR visualization etc.). In general, HyRAM+ overpredicts the flame lengths by around 65%. For heat and radiation models, we compared the heat flux and radiation data reported from two different sources to the values calculated by HyRAM+. For higher mass flow rates, the HyRAM+ calculated flame length results gave a better estimate of what is found in the experiments (65% error), but a higher error (85%) is observed between the HyRAM+ calculated lengths and the experimental flame lengthsfor lower mass flows. Some differences can be attributed to outdoor environmental effects (i.e. wind speed) and uncertainties in jet flame shapes. The propane flame trajectory is predicted for a high Reynolds number case with Re = 12,500 and a low Reynolds number case where Re = 2,000. The Re=12,500 case which is momentum dominated matches well with the experimental flame trajectory, but the agreement for the bouancy driven low Reynolds number case is not as good. Dispersion modeling for unignited propane was also analyzed. We compared the mole fraction, mixture fraction, mean velocity, concentration half width, and inverse mass concentration over an axial distance from different credible journals to the values calculated by HyRAM+. The results display good agreement but generally, HyRAM+ predicts a wider profile for mole fraction and mixture fraction experiments. Overall, HyRAM+’s results are reasonable for predicting the flame length, heat flux, flame trajectory, and dispersion for propane and can be used in risk analyses

08 HYDROGEN↗