Search NASA⌕ Search

SEARCH · Search NASA

Results for “Code assessment”

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 1,783 records · Page 99

Medium-Scale Methanol Pool Fire Model Validation

In this work, medium scale (30 cm diameter) methanol pool fires were simulated using the latest fire modeling suite implemented in Sierra/Fuego, a low Mach number multiphysics reacting flow code. The sensitivity of model outputs to various model parameters was studied with the objective of providing model validation. This work also assesses model performance relative to other recently published large eddy simulations (LES) of the same validation case. Two pool surface boundary conditions were simulated. The first was a prescribed fuel mass flux and the second used an algorithm to predict mass flux based on a mass and energy balance at the fuel surface. Gray gas radiation model parameters (absorption coefficients and gas radiation sources) were varied to assess radiant heat losses to the surroundings and pool surface. The radiation model was calibrated by comparing the simulated radiant fraction of the plume to experimental data. The effects of mesh resolution were also quantified starting with a grid resolution representative of engineering type fire calculations and then uniformly refining that mesh in the plume region. Simulation data were compared to experimental data collected at the University of Waterloo and the National Institute of Standards and Technology (NIST). Validation data included plume temperature, radial and axial velocities, velocity temperature turbulent correlations, velocity velocity turbulent correlations, radiant and convective heat fluxes to the pool surface, and plume radiant fraction. Additional analyses were performed in the pool boundary layer to assess simulated flame anchoring and the effect on convective heat fluxes. This work assesses the capability of the latest Fuego physics and chemistry model suite and provides additional insight into pool fire modeling for nonluminous, non-sooting flames.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

PCR Amplicon Prediction from Multiplex Degenerate Primer and Probe Sets

Assessing primer specificity and predicting both desired and off-target amplification products is an essential step for robust PCR assay design. Code is described to predict potential polymerase chain reaction (PCR) amplicons in a large sequence database such as NCBI nt from either singleplex or a large multiplexed set of primers, allowing degenerate primer and probe bases, with target mismatch annotates amplicons with gene information automatically downloaded from NCBI, and optionally it can predict whether there are also TaqMan/Luminex probe matches within predicted amplicons.

Gardner, S.N.↗

Preliminary Study of Iodine Gas Removal in Sodium Pools

Potential iodine gas release from failed fuel pins is a critical factor in the source term analysis of oxide fuel-loaded sodium fast reactors (SFRs). The accumulated iodine-containing gas mixtures inside pin plenums are expected to be ejected during pin failures and rise through sodium pool, with potential release of gaseous iodine to the cover gas region. Due to its potential radiological impacts, a proper assessment of iodine behavior is necessary for an accurate source term assessment. Throughout the bubble rise trajectory in the sodium pool, iodine gas is continuously removed or transformed at the bubble interface by diffusion, as the combining reaction between the iodine and sodium to form sodium iodide (NaI) is a chemically preferred process. As the final amount of iodine released from the facility is strongly influenced by the removal phenomenon inside the sodium pool, experiments were previously performed by PNC (Power Reactor and Nuclear Fuel Development Corporation) to provide insight into this phenomenon. To assess the accuracy of present methods for predicting iodine gas removal within sodium pools, several candidate approaches, available in source term analysis codes, have been summarized and evaluated in this study. Spherical cap bubbles and spherical bubbles are considered in accordance with the methods adopted in each approach, and different forms of correlations for major parameters have been implemented in accordance with the original adoptions. Based on the summarized results, important aspects to be considered have been derived.

Decontamination↗

Initial Comparison of Direct and Legacy Modeling Approaches for Radial Core Expansion Analysis

Radial core expansion in sodium-cooled fast reactors provides an important reactivity feedback effect. As the reactor power increases due to normal start up conditions or accident scenarios, the core and surrounding materials heat up, causing both grid plate expansion and bowing of the assembly ducts. When the core restraint system is designed correctly, the resulting structural deformations introduce negative reactivity which decreases the reactor power. Historically, an indirect procedure has been used to estimate the reactivity feedback due to structural deformation which relies upon perturbation theory and coupling legacy physics codes with limited geometry capabilities. With advancements in modeling and simulation, radial core expansion phenomena can now be modeled directly, providing an assessment of the accuracy of the reactivity feedback coefficients generated by indirect legacy methods. Recently a new capability was added to the PROTEUS-SN unstructured geometry neutron transport solver to analyze deformed meshes quickly and directly. By supplying the deformed mesh in addition to the base configuration input files, PROTEUS-SN automatically processes material adjustments including calculation of region densities to conserve mass, calculation of isotopic densities according to material models (for example, sodium density as a function of temperature), and subsequent re-homogenization of materials. To verify the new capability of directly simulating deformed meshes, PROTEUS-SN was used to compute reactivity feedback for a series of contrived yet representative deformed configurations for the Advanced Burner Test Reactor design. The indirect legacy procedure was also performed to generate reactivity feedback coefficients for the same deformed configurations. Interestingly, the legacy procedure consistently overestimated reactivity feedbacks by 35% compared to direct simulations by PROTEUS-SN. This overestimation indicates that the legacy procedures are in fact not conservative and could be overestimating reactivity feedback effects that are closely tied to reactor safety. We conclude that there is indeed value in performing direct simulation of deformed meshes despite the increased computational expense. PROTEUS-SN is already part of the SHARP multi-physics toolkit where both thermal hydraulics and structural mechanical feedback modeling can be applied but this is the first comparison of direct simulation to legacy techniques for radial core expansion.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Case Study: Nuclear-Renewable-Water Integration in Arizona

This document reports the application of the N-R HES software framework to a case study for Arizona Public Service (APS); the manager and part owner of the Palo Verde (PV) nuclear power plant. The case study is a work in progress of which this report presents a detailed description of the current model input data, assumptions and the corresponding results produced by the developed software framework. The goal of the report is, that with this information APS (together with INL) should be able to identify where more detailed data and more accurate model assumptions are necessary to reproduce all driving physical and economic phenomena in order to capture enough complexity of the reality to allow APS to use the results for their strategic decisions. The N-R HES software framework started being developed at INL two years ago [ref]. The framework has reached some level of maturity so that it can be applied to more than simple demonstration cases, e.g. real industry problems. Nevertheless, more capabilities are constantly added to accommodate the special needs of these challenging real life problems. The N-R HES framework is built on top of the RAVEN code [ref] and uses it as driver and workflow manager for all calculations. The framework has specifically been developed for the economic assessment of N-R HES systems. There are four main cornerstones of the N-R HES simulation framework: 1) generation of stochastic time series, 2) a probabilistic analysis and optimization set of algorithms available in RAVEN, 3) a set of models for representation of the physical behavior of N-R HES, and 4) a RAVEN plug-in that maps physical performance into economic performance called CashFlow. Within this framework, a broad spectrum of questions related to N-R HES can be addressed. One of the challenges currently of high interest is that the increasing penetration of variable renewables is altering the profile of the net demand (demand after removing all non curtailable renewable energy sources), with which the other generators on the grid have to cope. The N-R HES software framework is capable to analyze the potential feasibility of mitigating the resultant volatility in the net electricity demand. One possible solution to that problem currently intensively studied by the energy industry is adding stabilizing loads to the grid. These loads can be external industrial processes that are able to ramp production up and down quickly or adding variable internal loads that will lift the base load of the power plant. The latter is a possible solution currently studied by APS. APS if currently anticipating several challenges that will likely arise in the near future. One of them is, as mentioned, coping with the rapid growth of variable renewable energy (VRE) sources put on the grid. Although APS can sell its electricity for a fixed retail price to cover the local demand first, it also sells electricity in excess of the local demand at the Palo Verde energy hub. Due to the growing VRE penetration, the Palo Verde hub electricity spot price is negative more and more frequently. To mitigate this price volatility, APS is seeking to add more base load, so that the quantity of excess energy to be sold for prices potentially less than the internal retail price or even negative prices at the hub is reduced. A second challenge APS is facing is that their cooling water acquisition contract with the Sub Regional Operating Group (SROG) will expire soon and a renewal can only be done for a significantly higher price of the water. Therefore, APS is also seeking for alternative sources for their cooling water. An opportunity currently under investigation is to pump a limited quantity of brackish water from the regional ground water. Although much cheaper than the water from the new SROG contract, the salinity of the brackish water is so high that a blend of brackish and SROG water would need additional treatment to improve its quality for use in the cooling towers. An on-site reverse osmosis (RO) desalination plant is envisaged to reduce the salinity to an acceptable level. This report investigates the economic impact of such an RO plant. As one can see, the RO plant would help with both problems APS is facing: First, it would rise the APS base load helping to mitigate VRE induced hub price volatility. From this perspective, the RO plant can be seen as a stabilizing load on the grid as would be any other industrial process. Second, the RO’s goods, e.g. the clean water produced can be used by APS itself for a (hopefully) lower cost than 100% SROG water would be. The report will first describe the APS water procurement strategy in detail and lay out the possible alternative scenarios studied in this document. This is covered in Chapter 2. Next, in Chapter 3, the APS case model developed within the N-R HES software framework is described. This chapter includes also a description of the needed framework developments to accommodate the APS study.

97 MATHEMATICS AND COMPUTING↗

Quantitative Performance Assessment of Proxy Apps and Parents (ECP Proxy App Project Milestone ADCD-504-9)

This report presents highlights of these efforts. Section 2 describes work that has been done to compare the performance of proxy applications on AMD MI60 vs. Nvidia V100 GPUs. So far only a small set of ECP proxies are running on AMD GPUs, but we will continue to expand this analysis as additional proxies become available. We find that although the MI60 and V100 have nearly the same measured memory bandwidth, memory bound proxy app kernels perform 20-30% worse on the MI60. Further work is needed to refine these comparisons to determine whether the root cause is due to differences in the hardware, software stack, platform specific optimization, or some combination of the three. Section 3 describes our continuing effort to find methods to accurately assess the similarity of proxies and parents. We have recently seen very encouraging results using a cosine similarity metric. This technique uses the angle between two vectors of hardware performance counters to characterize the similarity (or difference) between two applications or proxies. We show not only that several widely used proxies are highly similar to their parents, but also that they differ from non-related codes. We also show that cosine similarity can be used to identify gaps and redundancies in suites and even to gain insight into the effects of architectural differences between platforms. Our work on assessing the Exascale toolchain is ongoing. Our successes with performance measurement tools are evident from the data provided in this report. However, our assessments across the broader tool chain are still too incomplete to provide a meaningful report at this time. We will continue to assess tools and work with vendors and third party developers as issues are identified.

97 MATHEMATICS AND COMPUTING↗

HED Hydrodynamics in the Common Modeling Framework

Within the Common Modeling Framework (CMF), the HED Hydrodynamics Project has implemented a common methodology, or "Authority", for modeling high energy density (HED) experiments by extending development of the Eulerian Applications Project code Authority, "EAPA," to include relevant physics models, customization tools and templates, and in situ and post-processing capabilities. Additionally, several models for small-scale HED experiments were successfully added to the CMF and simulations reproduce the results of preexisting, benchmarked input decks. We expect use of the CMF to improve our effectiveness in developing predictive capabilities for experiments ranging from small-scale planar single-interface, single-shock to multi-interface, multi-shock configurations and up through multi-shell inertial confinement fusion (ICF) implosions. Additionally, this will strengthen our ability to design future targeted experiments. Using the CMF infrastructure and Authorities for HED modeling has three key impacts. First, it enables a version-controlled, pedigreed and archived base model for each experiment with tools that allow easy setup, execution and analysis of simulations to assess and/or expose sensitivities to various parameters. Second, having a common base set of inputs and physics definitions helps identify systematic differences between models in the codes and experimental data to support code validation, and build confidence in our predictive capabilities. Finally, designed to work in concert with multiple programs and projects including PEM, IC, OES, V&V and DSW, the CMF supports improved collaboration and integration through shared model definitions, data, simulations and results, enhancing the workflow required to underwrite the knowledge base and capabilities ultimately required for stockpile stewardship

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Automation of FARM from Alpha Phase to Beta Phase

Integrated energy systems (IES) combine different energy technologies in synergistic ways to achieve a more secure and economical energy supply. The RAVEN-based HYBRID framework and the RAVEN plugin for grid and capacity optimization (HERON) are used to find the optimal installed capacity and the optimal economical dispatch of each component of the IES, by respecting the limits on the production variables and the corresponding rates of variation (explicit constraints). Besides, there are other process variables whose evolution needs to be bounded to avoid damaging the components (e.g., condensers, heat exchangers, steam generators, etc.) or degrading the process efficiency (e.g., electrolysis in the hydrogen production process). To avoid violating these latter limits (implicit constraints), a proof-of-concept HERON validator based on Feasible Actuator Range Modifier (FARM-Alpha) was developed by Argonne National Laboratory in January 2021. This FARM-Alpha validator calculates the evolution of process variables on whom the implicit constraints are placed, and then provides feedback to HERON dispatcher to adjust the power setpoints of three IES components, i.e., Balance of Plant, Secondary Energy Source, and Thermal Energy Storage, so as to meet both the explicit and implicit constraints. FARM-Alpha was designed to assess the performance of FARM as a HERON validator only, i.e., the list of components and implicit operational constraints were hard-coded within the source code. The lack of flexibility of the corresponding software structure does not allow the deployment in production environment. This report describes the development and the implementation of an enhanced version of the FARM-based validator (FARM-Beta), which ensures more flexibility for the end user in modeling multiple IES configurations and scenarios. Several test cases of the power dispatch problem were then selected to demonstrate the capabilities offered by FARM-beta. The test cases illustrate the efficiency of the closed-loop optimization scheme and the capability to calculate set-point trajectories satisfying both explicit and implicit constraints.

24 POWER TRANSMISSION AND DISTRIBUTION↗

MELCOR Code Change History: Revision 18019 to 21402

This document summarily provides brief descriptions of the MELCOR code enhancement made between code revision number 18019and 21440. Revision 18019 represents the previous official code release; therefore, the modeling features described within this document are provided to assist users that update to the newest official MELCOR code release, 21440. Along with the newly updated MELCOR Users’ Guide [2] and Reference Manual [3], users are aware and able to assess the new capabilities for their modeling and analysis applications.

97 MATHEMATICS AND COMPUTING↗

Assessing the Impact of Large Removable Beryllium Reflector Experiments on HFIR Performance and Safety Metrics

This study investigated the effect of various removable beryllium (RB) experiment configurations on High Flux Isotope Reactor (HFIR) metrics to address increasing interest in these facilities for materials and fuels irradiation research. The RB reflector contains eight large and four small irradiation experiment facilities that offer excellent neutron flux conditions to perform fission and fusion reactor materials and fuels irradiation research. This work aims to outline acceptable RB configurations based on safety, performance, and programmatic metrics. This study assessed experiment effects on reactivity, cycle length, fission rate density distributions, and neutron flux distributions using the Shift, HFIRCON, and SCALE ORIGEN codes. Initial evaluations focused on generic experiment materials including aluminum plugs, stainless steel plugs, molybdenum plugs, and aluminum plugs with gadolinium shields. Subsequent analyses of MiniFuel experiments were performed to evaluate a heterogenous experiment, consisting of a more complex geometry and bearing several materials, and to test the correlations developed with the generic materials on a real experiment. Furthermore, the effects of neutron poison concentrations in the standard beryllium plugs were evaluated. When assuming a reference RB configuration with eight large fresh beryllium plugs, perturbed configurations with three aluminum plugs, one stainless steel plug, one molybdenum plug, one gadolinium-shielded aluminum plug, and one MiniFuel experiment resulted in a cycle length reduction of less than 1.3 days, the current threshold before requiring additional approvals. Configurations with five aluminum plugs, one stainless steel plug, one molybdenum plug, one gadolinium-shielded aluminum plug, and two MiniFuel experiments meet the 1.3 day limit if irradiated beryllium plugs are considered. Fuel element fission rate density distributions remained within safety limits, with maximum local increases under 9%. Neutron flux calculations revealed large thermal flux depressions inside and around the perturbed RB facilities, while epithermal and fast neutron fluxes increased because of reduced neutron moderation by the perturbed materials. These findings provide valuable guidance for optimizing RB configurations to balance safety and performance at HFIR.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Anti-Ultralocality and Plateau Models of Inflation

Anti-ultralocality refers to the growth of spatial gradient terms relative to velocity terms in the coupled Einstein--scalar field equations. It is a characteristic feature of decelerated expansion before the onset of inflation. Previous numerical relativity studies have shown that anti-ultralocality prevents the onset of inflation in models with power-law inflaton potentials. In this paper, we show that models with plateau-shaped inflaton potentials, which are considered to be the simplest way to generate a tensor-to-scalar ratio below current observational upper limits, are especially vulnerable to anti-ultralocality effects. The reasons are the flatness of the plateau and the energy density gap of $\sim 10$ orders of magnitude between the Planck density and the plateau potential energy. To study the problem, we develop a protocol for assessing the viability of inflationary models in general, and we apply it to a plateau potential using a previously validated numerical relativity code. We find that, starting from generic initial conditions, the growth of gradient terms in the Einstein equations relative to non-gradient terms either prevents inflation from lasting for enough $e$-folds or triggers a phase of quantum runaway. We show that the fine-tuning of initial conditions necessary to avoid these issues becomes more severe as the energy scale of inflation is made smaller, disfavoring common approaches for reducing the tensor-to-scalar ratio.

FOS: Physical sciences↗

Measurement of recovery temperature on an airfoil in the Langley 0.3-m transonic cryogenic tunnel

Experimental measurements of recovery temperature were made on an airfoil in the Langley 0.3-m Transonic Cryogenic Tunnel at Mach numbers of 0.60 and 0.84 over a Reynolds number per meter range from about 15,000,000 to about 335,000,000. The measured recovery temperatures were considerably below those associated with ideal-gas ambient temperature wind tunnels. This difference was accentuated as the stagnation pressure increased and the total temperature decreased. A boundary-layer code modified for use with cryogenic nitrogen adequately predicted the measured adiabatic wall temperature at all conditions. A quantitative, on-line assessment of the nonadiabatic condition of a model can be made during the operation of a cryogenic wind tunnel by using a correlation for the adiabatic wall temperature which is only a function of total temperature, total pressure, and local Mach number on the model.

Johnson, C. B.↗

A computer simulation of aircraft evacuation with fire

A computer simulation was developed to assess passenger survival during the post-crash evacuation of a transport category aircraft when fire is a major threat. The computer code, FIREVAC, computes individual passenger exit paths and times to exit, taking into account delays and congestion caused by the interaction among the passengers and changing cabin conditions. Simple models for the physiological effects of the toxic cabin atmosphere are included with provision for including more sophisticated models as they become available. Both wide-body and standard-body aircraft may be simulated. Passenger characteristics are assigned stochastically from experimentally derived distributions. Results of simulations of evacuation trials and hypothetical evacuations under fire conditions are presented.

Middleton, V. E.↗

Analysis of turbulent underexpanded jets. II - Shock noise features using SCIPVIS

SCIPVIS, the computational model discussed by Dash et al. (1985), is assessed in predicting the complicated flow structure associated with shock-containing plumes. In addition, the analysis in this study examines this code's applicability as a basic part of a program for estimating broadband shock noise radiation. The results of this study show that excellent agreement exists between predicted and measured static pressure distributions for both underexpanded and overexpanded flow cases considered. Of the three turbulence closure models incorporated in the SCIPVIS code, the kW model of Spalding produces the most uniform agreement with measurement. The k-epsilon-2 model of Launder consistently overestimates plume spreading for supersonic jets with exit Mach numbers in the 1-2 range. Dash's (1983) k-epsilon-2-cc, compressibility-corrected version of Launder's model underestimates plume spreading. Good qualitative agreement was also obtained between the measured longitudinal turbulence intensity and that predicted by the code for the same trial case. Comparison of measured and predicted broadband shock noise spectrum peak values were found to be in excellent agreement. This utilized a variant of the Harper-Bourne and Fisher (1973) phase-array model: the effective shock spacing was reinterpreted as the value of the end of the plume potential core, determined herein by the SCIPVIS code.

Seiner, J. M.↗

A methodology for evaluating the reliability and risk of structures under complex service environments

The theoretical basis and numerical implementation of NESSUS (Numerical Evaluation of Stochastic Structures Under Stress; Cruse et al., 1988), a computer code for probabilistic structural analyses of aerospace components, are described, with an emphasis on the use of NESSUS for reliability and risk assessment. Topics addressed include the structure of NESSUS, the material-properties model, computational procedures, probabilistic models of fatigue-crack initiation, risk/cost evaluation, fatigue-fracture analyses, and fatigue-crack propagation. Numerical results from typical applications are presented in extensive graphs and briefly characterized; the usefulness of NESSUS predictions for establishing inspection and retirement schedules and for component certification is indicated.

Shiao, Michael C.↗

A methodology for evaluating the reliability and risk of structures under complex service environments

The theoretical basis and numerical implementation of NESSUS (Numerical Evaluation of Stochastic Structures Under Stress), a computer code for probabilistic structural analysis of aerospace components, are described, with an emphasis on the use of NESSUS for reliability and risk assessment. Topics addressed include the structure of probabilistic models of fatigue-crack initiation, risk/cost evaluation, fatigue-fracture analysis, and fatigue-crack initiation. Numerical results from typical applications are presented in graphs and briefly characterized. The usefulness of NESSUS predictions for establishing inspection and retirement schedules and for component certification is indicated.

Shiao, Michael C.↗

Assessment of antipodal-impact terrains on Mars

The regions anitpodal to Mars' three largest impact basins, Hellas, Isidis, and Argyre, were assessed for evidence of impact-induced disrupted terrains. Photogeology and computer modeling using the Simplified Arbitrary Lagrangian Eulerian (SALE) finite element code suggest that such terrains could have been found by the Hellas impact. Maximum antipodal pressures are 1100 MPa for Hellas, 520 MPa for Isidis, and 150 MPa for Argyre. The results suggest that if antipodal fracturing were associated with later volcanism, then Alba Patera may be related to the Hellas event, as proposed by Peterson (1978). Alba Patera is a unique volcano in the solar system, being a shield volcano which emitted large volume lava flows. This volcanism could be the result of the focusing of seismic energy which created a fractured region that served as a volcanic conduit for the future release of large volumes of magma. No disrupted terrain features are observed antipodal to the Isidis or Argyre basins, although some of the old fractures in Noctis Labyrinthus could have originated in response to the Isidis impact, and later have been reactivated by the Tharsis tectonics assumed to have produced Noctis. If the lower calculated antipodal pressures for Argyre were capable of producing disrupted terrains, then the terrains have been covered subsequently by volcanic or aeolian material, or modified beyond recognition.

Williams, David A.↗

A Performance Map for Ideal Air Breathing Pulse Detonation Engines

The performance of an ideal, air breathing Pulse Detonation Engine is described in a manner that is useful for application studies (e.g., as a stand-alone, propulsion system, in combined cycles, or in hybrid turbomachinery cycles). It is shown that the Pulse Detonation Engine may be characterized by an averaged total pressure ratio, which is a unique function of the inlet temperature, the fraction of the inlet flow containing a reacting mixture, and the stoichiometry of the mixture. The inlet temperature and stoichiometry (equivalence ratio) may in turn be combined to form a nondimensional heat addition parameter. For each value of this parameter, the average total enthalpy ratio and total pressure ratio across the device are functions of only the reactant fill fraction. Performance over the entire operating envelope can thus be presented on a single plot of total pressure ratio versus total enthalpy ratio for families of the heat addition parameter. Total pressure ratios are derived from thrust calculations obtained from an experimentally validated, reactive Euler code capable of computing complete Pulse Detonation Engine limit cycles. Results are presented which demonstrate the utility of the described method for assessing performance of the Pulse Detonation Engine in several potential applications. Limitations and assumptions of the analysis are discussed. Details of the particular detonative cycle used for the computations are described.

Paxson, Daniel E.↗