Developing and evaluating approaches for direct geologic disposal of commercial spent nuclear fuel (SNF) in dual-purpose canisters (DPCs) is a cross-cutting multi-disciplinary activity that is directly tied to the implementation of DPCs by the nuclear industry. The ultimate goal of the DPC direct disposal R&D program is to facilitate and maximize safe, cost-effective, licensed direct disposal. Independent Technical Review (ITR) is needed to maximize the impact of the R&D program on future implementation. The review will involve a team of experts representing the nuclear industry, repository sciences, and licensing. The team will be charged to review a set of representative technical reports and other information, and answer a set of questions that focus on R&D steering.
Macroscopic simulations of dense plasmas rely on detailed microscopic information that can be computationally expensive and is difficult to verify experimentally. In this work, we delineate the accuracy boundary between microscale simulation methods by comparing Kohn-Sham density functional theory molecular dynamics (KS-MD) and radial pair potential molecular dynamics (RPP-MD) for a range of elements, temperature, and density. By extracting the optimal RPP from KS-MD data using force matching, we constrain its functional form and dismiss classes of potentials that assume a constant power law for small interparticle distances. Our results show excellent agreement between RPP-MD and KS-MD for multiple metrics of accuracy at temperatures of only a few electron volts. The use of RPPs offers orders of magnitude decrease in computational cost and indicates that three-body potentials are not required beyond temperatures of a few eV. Due to its efficiency, the validated RPP-MD provides an avenue for reducing errors due to finite-size effects that can be on the order of ∼ 20 %.
Broderick, Robert J.; Reno, Matthew J.; Lave, Matthew S.; Azzolini, Joseph A.; Blakely, Logan; Galtieri, Jason; Mather, Barry; Weekley, Andrew; Hunsberger, Randolph; Chamana, Manohar; Li, Qinmiao; Zhang, Wenqi; Latif, Aadil; Zhu, Xiangqi; Grijalva, Santiago; Zhang, Xiaochen; Deboever, Jeremiah; Qureshi, Muhammad U.; Therrien, Francis; Lacroix, Jean-Sebastien; Li, Feng; Belletete, Marc; Hebert, Guillaume; Montenegro, Davis; Dugan, Roger
The rapid increase in penetration of distributed energy resources on the electric power distribution system has created a need for more comprehensive interconnection modeling and impact analysis. Unlike conventional scenario-based studies, quasi-static time-series (QSTS) simulations can realistically model time-dependent voltage controllers and the diversity of potential impacts that can occur at different times of year. However, to accurately model a distribution system with all its controllable devices, a yearlong simulation at 1-second resolution is often required, which could take conventional computers a computational time of 10 to 120 hours when an actual unbalanced distribution feeder is modeled. This computational burden is a clear limitation to the adoption of QSTS simulations in interconnection studies and for determining optimal control solutions for utility operations. The solutions we developed include accurate and computationally efficient QSTS methods that could be implemented in existing open-source and commercial software used by utilities and the development of methods to create high-resolution proxy data sets. This project demonstrated multiple pathways for speeding up the QSTS computation using new and innovative methods for advanced time-series analysis, faster power flow solvers, parallel processing of power flow solutions and circuit reduction. The target performance level for this project was achieved with year-long high-resolution time series solutions run in less than 5 minutes within an acceptable error.
Accurate and efficient constitutive modeling remains a cornerstone issue for solid mechanics analysis. Over the years, the LAMÉ advanced material model library has grown to address this challenge by implementing models capable of describing material systems spanning soft polymers to stiff ceramics including both isotropic and anisotropic responses. Inelastic behaviors including (visco)plasticity, damage, and fracture have all incorporated for use in various analyses. This multitude of options and flexibility, however, comes at the cost of many capabilities, features, and responses and the ensuing complexity in the resulting implementation. Therefore, to enhance confidence and enable the utilization of the LAMÉ library in application, this effort seeks to document and verify the various models in the LAMÉ library. Specifically, the broader strategy, organization, and interface of the library itself is first presented. The physical theory, numerical implementation, and user guide for a large set of models is then discussed. Importantly, a number of verification tests are performed with each model to not only have confidence in the model itself but also highlight some important response characteristics and features that may be of interest to end-users. Finally, in looking ahead to the future, approaches to add material models to this library and further expand the capabilities are presented
While studies of urban acoustics are typically restricted to the audio range, anthropogenic activity also generates infrasound (<20 Hz, roughly at the lower end of the range of human hearing). Shutdowns related to the COVID-19 pandemic unintentionally created ideal conditions for the study of urban infrasound and low frequency audio (20-500 Hz), as closures reduced human-generated ambient noise, while natural signals remained relatively unaffected. An array of infrasound sensors deployed in Las Vegas, NV, provides data for a case study in monitoring human activity during the pandemic through urban acoustics. The array records a sharp decline in acoustic power following the temporary shutdown of businesses deemed nonessential by the state of Nevada. This decline varies spatially across the array, with stations close to McCarran International Airport generally recording the greatest declines in acoustic power. Further, declines in acoustic power fluctuate with the time of day. As only signals associated with anthropogenic activity are expected to decline, this gives a rough indication of periodicities in urban acoustics throughout Las Vegas. The results of this study reflect the city's response to the pandemic and suggest spatiotemporal trends in acoustics outside of shutdowns.
Narrowband and broadband low frequency magnetic field sensors have applications in EMI/EMC measurements and atmospheric and space research. This report gives analytical development of such sensors.
The explosion of both sensors and GPS-enabled devices has resulted in position/time data being the next big frontier for data analytics. However, many of the problems associated with large numbers of trajectories do not necessarily have an analog with many of the historic big-data applications such as text and image analysis. Modern trajectory analytics exploits much of the cutting-edge research in machine-learning, statistics, computational geometry and other disciplines. We will show that for doing trajectory analytics at scale, it is necessary to fundamentally change the way the information is represented through a feature-vector approach. We then demonstrate the ability to solve large trajectory analytics problems using this representation.
The objective of the Photovoltaic Collaborative to Advance Multi-climate and Performance Research (PVCAMPER) is to: 1) Build and maintain a multi-climate research platform to enable pioneering photovoltaic research; 2) Validate the performance of emerging technologies in specific climates; 3) Help accelerate the world’s transition to a solar-intensive economy. Our focus in achieving those goals is to foster collaborative research and to build an international organization dedicated to improving data quality, minimizing measurement uncertainty and exchanging best practices related to PV performance.
Three M-MOF-74 (M = Co, Mg, Ni) metal-organic framework (MOF) thin film membranes have been synthesized through a sensor functionalization method for the direct electrical detection of NO2. The two-step surface functionalization procedure on the glass/Pt interdigitated electrodes resulted in a terminal carboxylate group, with both steps confirmed through infrared spectroscopic analysis. This surface functionalization allowed the MOF materials to grow largely in a uniform manner over the surface of the electrode forming a thin film membrane over the Pt sensing elec-trodes. The growth of each membrane was confirmed through scanning electron microscopy (SEM) and X-ray diffraction analysis. The Ni and Mg MOFs grew as a continuous but non-defect free membrane with overlapping polycrystallites across the glass surface, whereas the Co-MOF-74 grew dis-continuously. To demonstrate the use of these MOF membranes as an NO2 gas sensor, Ni-MOF-74 was chosen as it was consistently fabricated as the best thin and homogenous membrane, as confirmed by SEM. The membrane was exposed to 5 ppm NO2 and the impedance magnitude was observed to decrease 123× in 4 h, with a larger change in impedance and a faster response than the bulk material. Importantly, the use of these membranes as a sensor for NO2 does not require them to be defect-free, but solely continuous and overlapping growth.
The “how to” document guides the user through complicated aspects of software usage. It should supplement both the User’s manual and the Theory document, by providing examples and detailed discussion that reduce learning time for complex set ups. These documents are intended to be used together. We will not formally list all parameters for an input here – see the User’s manual for this. All the examples in the “How To” document are part of the Sierra/SD test suite, and each will run with no modification. The nature of this document casts together a number of rather unrelated procedures. Grouping them is difficult. Please try to use the table of contents and the index as a guide in finding the analyses of interest.
Garg, Raveesh; Qin, Eric; Martinez, Francisco M.; Guirado, Robert; Jain, Akshay; Abadal, Sergi; Abellan, Jose L.; Acacio, Manuel E.; Alarcon, Eduard; Rajamanickam, Sivasankaran R.; Krishna, Tushar
Recently, Graph Neural Networks (GNNs) have received a lot of interest because of their success in learning representations from graph structured data. However, GNNs exhibit different compute and memory characteristics compared to traditional Deep Neural Networks (DNNs). Graph convolutions require feature aggregations from neighboring nodes (known as the aggregation phase), which leads to highly irregular data accesses. GNNs also have a very regular compute phase that can be broken down to matrix multiplications (known as the combination phase). All recently proposed GNN accelerators utilize different dataflows and microarchitecture optimizations for these two phases. Different communication strategies between the two phases have been also used. However, as more custom GNN accelerators are proposed, the harder it is to qualitatively classify them and quantitatively contrast them. In this work, we present a taxonomy to describe several diverse dataflows for running GNN inference on accelerators. This provides a structured way to describe and compare the design-space of GNN accelerators.
The U.S. Department of Energy supports an R&D program for evaluating approaches to direct disposal of commercial spent fuel in dual-purpose canisters (DPCs). The major thrusts include alternative measures for treating the possibility of internal criticality events in DPC-based waste packages after thousands of years in a repository. These measures include: 1) injectable fillers, 2) analysis of the consequences of criticality events in a repository should they occur, and 3) options for modifying fuel assemblies or baskets in DPCs at the time they are loaded. This report presents a snapshot of progress in each of these areas drawing on deliverable reports generated during FY18 through FY20. Another aspect of the R&D program is to develop concepts of operations for repositories that would permanently dispose of DPC-based waste packages, considering different generic host media (not site-specific). The idea is to examine whether the disposal of large, heavy, heat-generating waste packages is technically feasible, and to identify the engineering challenges that would arise during implementation of the different disposal concepts. Descriptions of repository features are presented for repositories in salt media, argillite (clay/shale) media, crystalline (e.g., granitic) media, and unsaturated media (considering either alluvium or hard rock). Thermal management criteria for each concept are presented in terms of the maximum waste package thermal power at emplacement, when the repository could be opened, and the duration of repository emplacement operations. The overall message of this report is that direct disposal of commercial spent fuel is technically feasible in different types of geologic host media, but that thermal management and postclosure criticality impose different constraints on each concept. Engineering challenges are recognized and discussed. Treatment of postclosure criticality is identified as an important technical question that receives the majority of attention in the R&D program.
The application of hydrogen as an energy carrier has been expanding into industrial and transportation sectors enabling sustainable energy resources and providing a zero-emission energy infrastructure. The hydrogen supply infrastructure includes processes from production and storage, to transportation and distribution, to end use. Each portion of the hydrogen supply infrastructure is regulated by international, federal, state, and local entities. Regulations are enforced by entities which provide guidance and updates as necessary. While energy sources such as natural gas are currently regulated via the Code of Federal Regulations and United States Code, there might be some ambiguity as to which regulations are applicable to hydrogen and where regulatory gaps may exist. This report contains an overview of the regulations that apply to hydrogen, and those that may indirectly cover hydrogen as an energy carrier participating in a sustainable zero emission global energy system. As part of this effort, the infrastructure of hydrogen systems and regulation enforcement entities are defined, and a visual map and reference table are developed. This regulatory map and table can be used to identify the boundaries of federal oversight for each component of the hydrogen supply value chain which includes production, storage, distribution, and use.
Numerical simulations of metallic structures undergoing rapid loading into the plastic range require material models that accurately represent the response. In general, the material response can be seen as having four interrelated parts: the baseline response under slow loading, the effect of strain rate, the conversion of plastic work into heat and the effect of temperature. In essence, the material behaves in a thermal-mechanical manner if the loading is fast enough so when heat is generated by plastic deformation it raises the temperature and therefore influences the mechanical response. In these cases, appropriate models that can capture the aspects listed above are necessary. The matters of interest here are the elastic-plastic response and ductile failure behavior of 6061-T651 aluminum alloy under the conditions described above. The work was accomplished by first designing and conducting a material test program to provide data for the calibration of a modular $J_2$ plasticity model with isotropic hardening as well as a ductile failure model. Both included modules that accounted for temperature and strain rate dependence. The models were coupled with an adiabatic heating module to calculate the temperature rise due to the conversion of plastic work to heat. The test program included uniaxial tension tests conducted at room temperature, 150 and 300 C and at strain rates between 10–4 and 103 1/s as well as four geometries of notched tension specimens and two tests on specimens with shear-dominated deformations. The test data collected allowed the calibration of both the plasticity and the ductile failure models. Most test specimens were extracted from a single piece of plate to maintain consistency. Notched tension tests came from a possibly different plate, but from the same lot. When using the model in structural finite element calculations, element formulations and sizes different from those used to model the test specimens in the calibration are likely to be used. A brief investigation demonstrated that the failure model can be particularly sensitive to the element selection and provided an initial guide to compensate in a specific example.
Irradiance from a nuclear weapon can be the source of heat on gas infrastructure. This exposure when sufficiently intense can result in failure. An estimation tool for this behavior is the object of this study. A lumped capacity technique is employed to estimate the system temperature rise. The temperature rise is related to three possible outcomes. Two of the outcomes are relatively certain failure and relatively certain lack of effect. A large range of exposures are assessed with the model, and a relatively small number of cases are in the uncertain range. This model is presented as a tool that can be used in conjunction with a structural assessment model to sensitivities to the overpressure and shock to screen potential outcomes from subject events .
The US generates approximately 2,000 Metric Tons of Heavy Metal (MTHM) of commercial spent fuel (SNF) each year and currently stores ~ 85,000 MTHM of commercial SNF at 70+ reactor sites, the disposal of which is the responsibility of the US Department of Energy (DOE). SNF is initially stored in spent nuclear fueld pools (SFPs). SFPs were initially constructed by US utilities for temporary fuel storage, but with no final disposal pathway available, SFPs are reaching capacity. To allow continued operation of the nation's commercial nuclear reactor fleet, utilities started transferring SNF from SFPs (wet storage) to dry cask storage systems, typically using dual-purpose (storage and transportation) canisters (DPCs). And while DPCs were designed, licensed and loaded to meet Nuclear Regulatory Commission (NRC) requirements that preclude the possibility of a critical event during SNF storage and transport, they were not designed or loaded to preclude the possibility of a criticality event during the regulated post-closure period following disposal, which could be up to 1,000,000 years (Price, 2019). DPC filler option criteria are detailed and materials that exhibit these attributes are explored. This document is an update of the SNL Joint Workplan on Filler Investigations for DPCs.
This report provides a summary of notes for building and running the Sandia Computational Engine for Particle Transport for Radiation Effects (SCEPTRE) code. SCEPTRE is a general- purpose C++ code for solving the li near Boltzmann transport equation in serial or parallel using unstructured spatial finite elements, multigroup energy treatment, and a variety of angular treatments including discrete ordinates and spherical harmonics. Either the first-order form of the Boltzmann equation or one of the second-order forms may be solved. SCEPTRE requires a small number of open-source Third Part y Libraries (TPL) to be available, and example scripts for building these TPLs are provided. The TPLs needed by SCEPTRE are Trilinos, boost, and netcdf. SCEPTRE uses an autotools build system , and a sample configure script is provided. Running the SCEPTRE code requires that the user provide a spatial finite-elements mesh in Exodus format and a cross section library in a format that will be described. SCEPTRE uses an xml-based input, and several examples will be provided.
NaI-AlBr3 is a very appealing low melting temperature (<100 C), salt system for use as an electrochemically-active electrolyte. This system was investigated for its electrochemical and physical properties with focus to energy storage considerations. A simple phase diagram was generated; at >100 C, lower NaI concentrations had two partially miscible liquid phases, while higher NaI concentrations had solid particles. Considering the fully molten regime, electrical conductivities were evaluated over 5-25 mol% NaI and 110 C-140 C. Conductivities of 6.8-38.9 mS cm-1 were observed, increasing with temperature and NaI concentration. Effective diffusion coefficients of the I-/I3- redox species were found to decrease with both increasing NaI concentration and increasing applied potential. Regardless, oxidation current density at 3.6 V vs Na/Na+ was observed to increase with increasing NaI concentration over 5-25 mol%. Finally, the critical interface between the molten salt electrolyte and electrode materials was found to significantly affect reaction kinetics. When carbon was used instead of tungsten, an adsorbed species, most likely I2, blocked surface sites and significantly decreased current densities at high potentials. This study shows the NaI-AlBr3 system offers an attractive, low-temperature molten salt electrolyte that could be useful to many applied systems, though composition and electrode material must be considered.
This document includes details of the angular quadrature sets available in SCEPTRE for performing numerical integrations in the angular phase space. The angular dependence of the boundary and fixed-source terms an d initial angular flux are specified by angular index rather than by direction. It is, therefore, necessary to know the mapping from a specific direction to a direction index. This document includes angular quadrature weights and direction cosines for most of the quadrature sets available in SCEPTRE.
This document includes the information pertaining to the practices of the Emergency Response Software Development Team (ERSDT) for the RASCAL software project. The content includes information regarding: 1) Software Quality Assurance. 2) Project Management. 3) Configuration Management. 4) Development Standards. 5) Third-Party Software. 6) Verification and Validation. The information contained in this report is considered living documentation. The information contained in this report was assembled from multiple documents and content management systems.
Learning 3D representations that generalize well to arbitrarily oriented inputs is a challenge of practical importance in applications varying from computer vision to physics and chemistry. We propose a novel multi-resolution convolutional architecture for learning over concentric spherical feature maps, of which the single sphere representation is a special case. Our hierarchical architecture is based on alternatively learning to incorporate both intra-sphere and inter-sphere information. We show the applicability of our method for two different types of 3D inputs, mesh objects, which can be regularly sampled, and point clouds, which are irregularly distributed. We also propose an efficient mapping of point clouds to concentric spherical images, thereby bridging spherical convolutions on grids with general point clouds. We demonstrate the effectiveness of our approach in improving state-of-the-art performance on 3D classification tasks with rotated data.
Alternatives to conventional diesel electric propulsion are currently of interest to rail operators. In the U.S., smaller railroads have implemented natural gas and other railroads are exploring hydrogen technology as a cleaner alternative to diesel. Diesel, battery, hydrogen fuel cell, or track electrification all have trade-offs for operations, economics, safety, and public acceptability. A framework to compare different technologies for specific applications is useful to optimize the desired results. Standards from the Association of American Railroads (AAR) and other industry best practices were reviewed for applicability with hydrogen fuel cell technology. Some technical gaps relate to the physical properties of hydrogen, such as embrittlement of metals, invisible flames, and low liquid temperatures. A reassessment of material selection, leak/flame detection, and thermal insulation methods is required. Hydrogen is less dense and diffuses more easily than natural gas, and liquid hydrogen is colder than liquefied natural gas. Different densities between natural gas and hydrogen require modifications to tank designs and flow rates. Leaked hydrogen will rise rather than pool on the ground like diesel, requiring a modification to the location of hydrogen tanks on rolling stock. Finally, the vibration and shock experienced in the rail environment is higher than light-duty vehicles and stationary applications for which current fuel cell technology has been developed, requiring a modification in tank design requirements and testing.
Even with the advent of additive manufacturing, the vast majority of complex structures are comprised of individual components held together with bolted joints. However, bolted joints present a challenge for mechanical design as they are a source of nonlinearity and increase the uncertainty in the overall behavior of the system in a dynamic environment. While many advances have been made in the ability to accurately model and test bolted joints, it is still an open area of research. Modes of vibration that exercise bolted joints typically exhibit nonlinear behavior where, with increased excitation level, the natural frequency decreases (i.e. softens) and the damping increases. However, the system under study for this work has an axial mode which does not follow this trend; it does soften as expected, but, after an initial increase, the apparent damping decreases with excitation amplitude. At the highest excitation level, the frequency of the mode decreases to that of a nearby bending mode and the response is amplified nearly 500% above that at lower levels. It is unclear whether the decrease in damping is due to the coupling of the two modes or if it is a characteristic of the axial mode. Therefore, the objective of this project is to investigate the coupling between the axial and bending modes and the dynamics leading to the decrease in damping.
Ekanayaka, Thilini K.; Hao, Guanhua; Mosey, Aaron; Dale, Ashley S.; Jiang, Xuanyuan; Yost, Andrew J.; Sapkota, Keshab R.; Wang, George T.; Zhang, Jian; N'Diaye, Alpha T.; Marshall, Andrew; Cheng, Ruihua; Naeemi, Azad; Xu, Xiaoshan; Dowben, Peter A.
Nonvolatile, molecular multiferroic devices have now been demonstrated, but it is worth giving some consideration to the issue of whether such devices could be a competitive alternative for solid‐state nonvolatile memory. For the Fe (II) spin crossover complex [Fe{H2B(pz)2}2(bipy)], where pz = tris(pyrazol‐1‐yl)‐borohydride and bipy = 2,2′‐bipyridine, voltage‐controlled isothermal changes in the electronic structure and spin state have been demonstrated and are accompanied by changes in conductance. Higher conductance is seen with [Fe{H2B(pz)2}2(bipy)] in the high spin state, while lower conductance occurs for the low spin state. Plausibly, there is the potential here for low‐cost molecular solid‐state memory because the essential molecular thin films are easily fabricated. However, successful device fabrication does not mean a device that has a practical value. Here, we discuss the progress and challenges yet facing the fabrication of molecular multiferroic devices, which could be considered competitive to silicon.
Rembold, Randy K.; Merchant, Bion J.; Davis, J.P.; Ebeling, Carl W.; Wilson, David C.; Ringler, Adam T.; Anthony, Robert E.
The Comprehensive Nuclear-Test-Ban Treaty mandates the four verification technologies to be used by the International Monitoring System (IMS) to monitor compliance to the Treaty. These technologies are seismic, hydroacoustic, infrasound, and radionuclide. Operational manuals for each of these technologies specify the requirements that equipment installed at IMS stations must meet for data from said stations to be accepted by the International Data Center. Following a model in which a revised set of infrasound sensor specifications and accompanying test procedures were developed in an international collaboration overseen by the Comprehensive Nuclear-Test-Ban Treaty Organization, the same process was recently followed for seismometers potentially of use by the IMS. This document, a product of that collaboration, defines key concepts and recommends test procedures to be followed to characterize broadband seismometers in a consistent and standardized way. Results from these test procedures can be used to verify that a seismometer meets the manufacturers stated specifications and IMS requirements. ACKNOWLEDGEMENTS The writers of this document would like to acknowledge Charles R. Hutt, John R. Evans, Fred Followill, Robert L. Nigbor, and Erhard Wielandt, the authors of the 2009 Guidelines for Standardized Testing of Broadband Seismometers and Accelerometers funded by the U.S. Department of Interior, on which this document is based. Although this document is limited in scope to seismometers and much of the document has been updated, the structure remains essentially the same.
The atomic cluster expansion is a general polynomial expansion of the atomic energy in multi-atom basis functions. Here we implement the atomic cluster expansion in the performant C++ code PACE that is suitable for use in large scale atomistic simulations. We briefly review the atomic cluster expansion and give detailed expressions for energies and forces as well as efficient algorithms for their evaluation. We demonstrate that the atomic cluster expansion as implemented in PACE shifts a previously established Pareto front for machine learning interatomic potentials towards faster and more accurate calculations. Moreover, general purpose parameterizations are presented for copper and silicon and evaluated in detail. We show that the new Cu and Si potentials significantly improve on the best available potentials for highly accurate large-scale atomistic simulations.
This report provides results from a series of 2-m pool fire experiments performed in the Thermal Test Complex at Sandia National Laboratories testing heptane, Bakken crude oil, and dilbit crude oil. The effect of the presence and placement of a calorimeter, fuel supply temperature, and maintaining a constant fuel level were assessed. Measurements include burn rate, surface emissive power, flame height, heat flux to an engulfed calorimeter, heat flux to external instruments, thermocouple temperatures within the fuel and fire plume, and heat release rate. The results indicate that the presence and placement of the calorimeter has the most effect on the measured quantities for the Bakken crude oil and indicated no effect for the Dilbit crude oil. The fuel feed temperature had a slight effect for the heptane fuel, but not for the crude oils. Allowing the fuel to burn down did not have a significant effect on any of the fuels. The Bakken crude oil resulted in the highest average total heat flux to the calorimeter by a factor of about 1.5 and 1.3 higher compared to heptane and the dilbit crude oil, respectively.
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