SEARCH · Engineering Papers

Results for “Conditional rate functions”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 415 records · Page 23

All sources NASA NTRS DOE OSTI

NASA NTRSDate not supplied

ROCKE-3D 2.0: An Updated General Circulation Model for Simulating the Climates of Rocky Planets

We present the second generation of ROCKE-3D (Resolving Orbital and Climate Keys of Earth and Extraterrestrial Environments with Dynamics), the generalized 3-dimensional General Circulation Model (GCM) for use in Solar System and exoplanetary simulations of rocky planet climates. ROCKE-3D version 2.0 is a descendant of GISS ModelE2.1, the flagship GCM of the NASA Goddard Institute for Space Studies (GISS) used in the most recent Intergovernmental Panel for Climate Change (IPCC) assessment. ROCKE-3D is a continuous effort to expand the capabilities of GISS ModelE to handle a broader range of atmospheric conditions, including different atmospheric planet sizes, gravities, pressures, rotation rates, more diverse chemistries and atmospheric compositions, diverse ocean and land distributions and topographies, and potential basic biosphere functions. We present in this release many more supported configurations which can serve as starting points to simulate terrestrial planets of interest. Two different radiation schemes are supported, the GISS radiation, valid for Earth-like planets, and SOCRATES, which is more generalized but requires more CPU resources. In addition, three different atmospheric compositions are supported, three ocean configurations (prescribed, q-flux, and dynamic), and two resolutions: the medium resolution (4 x 5 degrees in latitude x longitude) used in ROCKE-3D 1.0, and the fine resolution, which has double the resolution in the atmosphere and 4 times the resolution in the ocean. We quantify how the different component choices affect model results, and discuss strengths and limitations of using each component, together with how one can select which component to use. ROCKE-3D is publicly available and tutorial sessions are available for the community, greatly facilitating its use by any interested group.

ROCKE-3D↗

DOE OSTI2020

Calculation of electron-ion temperature equilibration rates and friction coefficients in plasmas and liquid metals using quantum molecular dynamics

In this work, we discuss a method to calculate with quantum molecular dynamics simulations the rate of energy exchanges between electrons and ions in two-temperature plasmas, liquid metals, and hot solids. Promising results from this method were recently reported for various materials and physical conditions. Like other ab initio calculations, the approach offers a very useful comparison with the experimental measurements and permits an extension into conditions not covered by the experiments. The energy relaxation rate is related to the friction coefficients felt by individual ions due to their nonadiabatic interactions with electrons. Each coefficient satisfies a Kubo relation given by the time integral of the autocorrelation function of the interaction force between an ion and the electrons. These Kubo relations are evaluated using the output of quantum molecular dynamics calculations in which electrons are treated in the framework of finite-temperature density functional theory. The calculation presents difficulties that are unlike those encountered with the Kubo formulas for the electrical and thermal conductivities. In particular, the widely used Kubo-Greenwood approximation is inapplicable here. Indeed, the friction coefficients and the energy relaxation rate diverge in this approximation since it does not properly account for the electronic screening of electron-ion interactions. The inclusion of screening effects considerably complicates the calculations. We discuss the physically motivated approximations we applied to deal with these complications in order to investigate a widest range of materials and physical conditions. Unlike the standard method used for the electronic conductivities, the Kubo formulas are evaluated directly in the time domain and not in the energy domain, which spares one from needing to introduce an extraneous undetermined numerical parameter to account for the discrete character of the numerical density of states. We highlight interesting properties of the energy relaxation rate not shared by other electronic properties, in particular its self-averaging character. We then present a detailed parametric and convergence study with the numerical parameters, including the system size, the number of bands and k-points, and the physical approximations for the dielectric function and the exchange-correlation energy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Acidic Oxygen Evolution Reaction Activity–Stability Relationships in Ru-Based Pyrochlores

Ru-based oxygen evolution reaction (OER) catalysts show significant promise for efficient water electrolysis, but rapid degradation poses a major challenge for commercial applications. In this work, we explore several Ru-based pyrochlores (A 2 Ru 2 O 7 , A = Y, Nd, Gd, Bi) as OER catalysts and demonstrate improved activity and stability of catalytic Ru sites relative to RuO 2 . Furthermore, we combine complementary experimental and theoretical analysis to understand how the A-site element impacts activity and stability under acidic OER conditions. Among the A 2 Ru 2 O 7 studied herein, we find that a longer Ru–O bond and a weaker interaction of the Ru 4d and O 2p orbitals compared with RuO 2 results in enhanced initial activity. We observe that the OER activity of the catalysts changes over time and is accompanied by both A-site and Ru dissolution at different relative rates depending on the identity of the A-site. Pourbaix diagrams constructed using density functional theory (DFT) calculations reveal a driving force for this experimentally observed dissolution, indicating that all compositions studied herein are thermodynamically unstable in acidic OER conditions. Theoretical activity predictions show consistent trends between A-site cation leaching and OER activity. These trends coupled with Bader charge analysis suggest that dissolution exposes highly oxidized Ru sites that exhibit enhanced activity. Overall, using the stability number (molO2 evolved/molRu dissolved) as a comparative metric, the A 2 Ru 2 O 7 materials studied in this work show substantially greater stability than a standard RuO2 and commensurate stability to some Ir mixed metal oxides. Finally, the insights described herein provide a pathway to enhanced Ru catalyst activity and durability, ultimately improving the efficiency of water electrolyzers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2015

Rate Theory Modeling and Simulations of Silicide Fuel at LWR Conditions

Uranium silicide (U 3 Si 2 ) fuel has higher thermal conductivity and higher uranium density, making it a promising candidate for the accident-tolerant fuel (ATF) used in light water reactors (LWRs). However, previous studies on the fuel performance of U 3 Si 2 , including both experimental and computational approaches, have been focusing on the irradiation conditions in research reactors, which usually involve low operation temperatures and high fuel burnups. Thus, it is important to examine the fuel performance of U 3 Si 2 at typical LWR conditions so as to evaluate the feasibility of replacing conventional uranium dioxide fuel with this silicide fuel material. As in-reactor irradiation experiments involve significant time and financial cost, it is appropriate to utilize modeling tools to estimate the behavior of U 3 Si 2 in LWRs based on all those available research reactor experimental references and state-of-the-art density functional theory (DFT) calculation capabilities at the early development stage. Hence, in this report, a comprehensive investigation of the fission gas swelling behavior of U 3 Si 2 at LWR conditions is introduced. The modeling efforts mentioned in this report was based on the rate theory (RT) model of fission gas bubble evolution that has been successfully applied for a variety of fuel materials at devious reactor conditions. Both existing experimental data and DFT-calculated results were used for the optimization of the parameters adopted by the RT model. Meanwhile, the fuel-cladding interaction was captured by the coupling of the RT model with simplified mechanical correlations. Therefore, the swelling behavior of U 3 Si 2 fuel and its consequent interaction with cladding in LWRs was predicted by the rate theory modeling, providing valuable information for the development of U 3 Si 2 fuel as an accident-tolerant alternative for uranium dioxide.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

NASA NTRS1978

Influence of hydrogen environments on crack growth in Inconel 718

The effect of hydrogen environments on sustained-load and cyclic-load crack growth in Inconel 718 was investigated using fracture-mechanics-type specimens. The sustained-load crack growth was determined to be a function of heat-treatment condition, temperature, and hydrogen pressure. The threshold stress intensity for subcritical crack growth was independent of hydrogen pressure at pressures greater than 21 MN/sq m (3000 psi). The cyclic-load crack growth rate in Inconel 718 at low and moderate stress intensity ranges increased with increasing hydrogen pressure. Decreasing the cyclic frequency from 1.0 to 0.1 Hz considerably increased the cyclic crack growth rate at 0.069 MN/sq m (10 psi) and at 68.9 MN/sq m (10,000 psi) hydrogen pressure. A series of measurements performed between 1.0 and 0.1 Hz showed that the crack growth rate increased as a complex function of the time per cycle.

Walter, R. J.↗

DOE OSTI2026

Transient Multiphysics Simulations with Pin Power Reconstruction in the Griffin Reactor Physics Code

This work introduces the pin power reconstruction capability available in the Griffin reactor physics code. This capability is implemented in an unstructured mesh framework, and the methods introduced are applied to the 2D SIMBA reactor core, which has assemblies and pins arranged in a hexagonal lattice. Since this reactor has a non-Cartesian geometry and also operates in the thermal spectrum, a general approach to pin power reconstruction is adopted, where SPH-based equivalence is leveraged to preserve assembly-wise reaction rates, while computing full-core form functions to preserve pin-wise fission production rates within the fuel pins of the reactor core. In a 2D microreactor benchmark problem, this pin power reconstruction approach was shown to reproduce pin powers compared to the Serpent2 Monte Carlo code for fixed temperature conditions and control drum rotation angles, yielding a core-wide RMS error level of 0.6\% and a maximum absolute pin error of 2.3\%. In addition, a tabulated library of multigroup cross sections, SPH factors, and form functions was generated to demonstrate the applicability of pin power reconstruction to a thermal feedback problem. Finally, a control drum transient was successfully simulated, showcasing the application of pin power reconstruction in a transient multiphysics feedback problem.

97 - MATHEMATICS AND COMPUTING↗

DOE OSTI2020

Radiolytic degradation of formic acid and formate in aqueous solution: modeling the final stages of organic mineralization under advanced oxidation process conditions

The successful use of advanced oxidation processes to treat aqueous solutions containing undesirable organic species requires the degradation of these species to lower molecular weight, lower hazard compounds. Safe application of this technology requires a thorough understanding of the mechanisms of degradation. These oxidative transformations are mainly initiated by the reactions of reactive oxygen species, particularly hydroxyl radicals. These react with organic molecules to generate carbon-centered radicals. In the presence of dissolved oxygen, the carbon-centered radicals are next converted to peroxyl radicals, which then decay to lower molecular weight species by multiple mechanistic pathways. Formic acid and its conjugate base formate are the last stable chemical species produced immediately before the complete mineralization of any organic molecule undergoing oxidative degradation in aqueous solution. Once understood, the radical-induced chemistry of formic acid/formate under these conditions has wide applicability in all advanced oxidation technologies. To develop this quantitative knowledge, here we have performed a series of 60 Co gamma irradiation studies on aqueous formic acid/formate over different pH and solution conditions. The measured species concentration changes, as a function of applied dose, are compared with the predictions of a kinetic computer model constructed from literature reactions and reported rate coefficients. The excellent agreement found between the results and modeling gives confidence in the mechanism presented here and provide the first complete computer model for the radiolytic degradation of formic acid in water.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Evolutionary reinforcement learning of dynamical large deviations

In this work, we show how to bound and calculate the likelihood of dynamical large deviations using evolutionary reinforcement learning. An agent, a stochastic model, propagates a continuous-time Monte Carlo trajectory and receives a reward conditioned upon the values of certain path-extensive quantities. Evolution produces progressively fitter agents, potentially allowing the calculation of a piece of a large-deviation rate function for a particular model and path-extensive quantity. For models with small state spaces, the evolutionary process acts directly on rates, and for models with large state spaces, the process acts on the weights of a neural network that parameterizes the model's rates. This approach shows how path-extensive physics problems can be considered within a framework widely used in machine learning.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Modeling pre-Exascale AMR Parallel I/O Workloads via Proxy Applications

The present work investigates the modeling of preexascale input/output (I/O) workloads of Adaptive Mesh Refinement (AMR) simulations through a simple proxy application. We collect data from the AMReX Castro framework running on the Summit supercomputer for a wide range of scales and mesh partitions for the hydrodynamic Sedov case as a baseline to provide sufficient coverage to the formulated proxy model. The non-linear analysis data production rates are quantified as a function of a set of input parameters such as output frequency, grid size, number of levels, and the Courant-Friedrichs-Lewy (CFL) condition number for each rank, mesh level and simulation time step. Linear regression is then applied to formulate a simple analytical model which allows to translate AMReX inputs into MACSio proxy I/O application parameters, resulting in a simple “kernel” approximation for data production at each time step. Results show that MACSio can simulate actual AMReX nonlinear “static” I/O workloads to a certain degree of confidence on the Summit supercomputer using the present methodology. The goal is to provide an initial level of understanding of AMR I/O workloads via lightweight proxy applications models to facilitate autotune data management strategies in anticipation of exascale systems.

Godoy, William↗

DOE OSTI2019

On Gas Ingression of Hermetic Packages

In this paper, a closed-form mathematic equation that governs gas ingression of hermetic packages is derived from first principles and applied to moist air and water vapor ingression conditions. The equation models internal gas partial pressure change as a function of time, external conditions, and package characteristics. The equation provides the theoretical basis for direct comparisons of ingression behaviors of different gases into hermetic packages. Comparing the rates of internal air pressure increase due to air ingression and water vapor partial pressure buildup due to water vapor ingression, the authors theorize that vacuum decay may present a greater challenge to the performance of microelectromechanical systems (MEMS) devices within hermetic packages than that of water vapor content induced corrosion failures. Furthermore, this paper also examines gas ingression of hermetic enclosures with multiple layers of seals.

42 ENGINEERING↗

DOE OSTI2016

Equivalent radiation damage in zirconium irradiated in various reactors - 2016-0101

The physical properties and in-reactor performance of core components made from zirconium alloys are a function of the radiation-induced atomic displacements that occur during service. The vacancy and interstitial point defects and their clusters created by the atomic displacements modify the microstructure, and therefore physical properties such as strength and ductility, by forming dislocation loops. The rate of production of point defects also affects time-dependent properties such as creep and diffusion (mass transport). To simulate the effect of irradiation on core components, researchers have used materials test reactors for many years to irradiate and test materials. Such reactor tests can target an accumulated displacement damage dose to simulate the state of a material at a given stage in reactor life but at an accelerated rate. Each nuclear reactor has a specific neutron energy spectrum that also depends on the in-reactor location. Therefore, to use test reactor data to determine the effect of radiation damage on power reactor materials, one must convert the accumulated dose into a unit that is common in its effect on the material properties. For many property changes in nuclear reactor cores, this unit is displacements per atom. This paper compares the damage production for zirconium irradiated under differing neutron energy spectra in various power and test reactors. The damage rates as a function of different measures of neutron flux were used to relate atomic displacements in various materials test reactors to CANDU and light-water reactor conditions. (authors)

36 MATERIALS SCIENCE↗

DOE OSTI2021

Assessing climate change impacts on live fuel moisture and wildfire risk using a hydrodynamic vegetation model

Abstract. Live fuel moisture content (LFMC) plays a critical role in wildfire dynamics, but little is known about responses of LFMC to multivariate climate change, e.g., warming temperature, CO2 fertilization, and altered precipitation patterns, leading to a limited prediction ability of future wildfire risks. Here, we use a hydrodynamic demographic vegetation model to estimate LFMC dynamics of chaparral shrubs, a dominant vegetation type in fire-prone southern California. We parameterize the model based on observed shrub allometry and hydraulic traits and evaluate the model's accuracy through comparisons between observed and simulated LFMC of three plant functional types (PFTs) under current climate conditions. Moreover, we estimate the number of days per year of LFMC below 79 % (which is a critical threshold for wildfire danger rating of southern California chaparral shrubs) from 1960 to 2099 for each PFT and compare the number of days below the threshold for medium and high greenhouse gas emission scenarios (RCP4.5 and 8.5). We find that climate change could lead to more days per year (5.2 %–14.8 % increase) with LFMC below 79 % between the historical (1960–1999) and future (2080–2099) periods, implying an increase in wildfire danger for chaparral shrubs in southern California. Under the high greenhouse gas emission scenario during the dry season, we find that the future LFMC reductions mainly result from a warming temperature, which leads to 9.1 %–18.6 % reduction in LFMC. Lower precipitation in the spring leads to a 6.3 %–8.1 % reduction in LFMC. The combined impacts of warming and precipitation change on fire season length are equal to the additive impacts of warming and precipitation change individually. Our results show that the CO2 fertilization will mitigate fire risk by causing a 3.5 %–4.8 % increase in LFMC. Our results suggest that multivariate climate change could cause a significant net reduction in LFMC and thus exacerbate future wildfire danger in chaparral shrub systems.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1976

Engineering flowfields and heating rates for highly swept wing leading edges

An engineering flow model is developed (and verified experimentally) which describes the flowfield that arises when a supersonic stream encounters a wedge/cylinder configuration whose angles are such that the flow includes only weak shock waves. A numerical code using the perfect gas relations is used to describe the flow in the plane of symmetry inboard of the shock interaction region. Theoretical surface-pressure and heat-transfer distributions are computed for freestream velocities ranging from 1167 to 7610 m/sec. Nondimensionalization of the heat-transfer rates in terms of local flow parameters produced a correlation of Stanton number as a function of the local Reynolds number, which is independent of the freestream flow conditions and of the surface temperature.

Bertin, J. J.↗

NASA NTRS1976

On the energy deposition of photoelectrons in the atmosphere of Venus

Vertical components of photoelectron fluxes in the atmosphere of Venus are computed by solving an appropriate form of the Boltzmann equation in the cases where there is no flux of either photoelectrons or solar-wind particles across the ionopause and where photoelectrons are free to escape from the atmosphere. It is assumed that Venus has no magnetic field and that the atmosphere is composed of carbon dioxide, atomic oxygen, and helium. The results are plotted as a function of altitude for several energies in the range from 100 eV to a cutoff of the order of 1 eV. Heating rates for the two upper boundary conditions and the case of no spatial transport are determined which show that transport effects dominate at altitudes greater than about 200 km. Electron temperatures are calculated for the adopted model atmosphere and ionosphere by solving the pertinent conservation equation, and excitation rates are computed for the CO Cameron band as well as the CO2(+) A and B bands.

Mccormick, P. T.↗

NASA NTRS1981

Effect of gamma irradiation on the friction and wear of ultrahigh molecular weight polyethylene

The effect of sterilization gamma irradiation on the friction and wear properties of ultrahigh molecular weight polyethylene (UHMWPE) sliding against stainless steel 316L in dry air at 23 C is investigated, the results to be used in the development of artificial joints which are to surgically replace diseased human joints. A pin-on-disk sliding friction apparatus is used, a constant sliding speed in the range 0.061-0.27 m/s is maintained, a normal load of 1 kgf is applied with dead weight, and the irradiation dose levels are: 0, 2.5, and 5.0 Mrad. Wear and friction data and conditions for each of the ten tests are summarized, and include: (1) wear volume as a function of the sliding distance for the irradiation levels, (2) incremental wear rate, and (3) coefficient of friction as a function of the sliding distance. It is shown that (1) the friction and wear properties of UHMWPE are not significantly changed by the irradiation doses of 2.5 and 5.0 Mrad, (2) the irradiation increases the amount of insoluble gel as well as the amount of low molecular weight material, and (3) after run-in the wear rate is either steady or gradually decreases as a function of the sliding distance.

Jones, W. R.↗

NASA NTRS1988

Semi-discrete approximations to nonlinear systems of conservation laws; consistency and L(infinity)-stability imply convergence

A convergence theory for semi-discrete approximations to nonlinear systems of conservation laws is developed. It is shown, by a series of scalar counter-examples, that consistency with the conservation law alone does not guarantee convergence. Instead, a notion of consistency which takes into account both the conservation law and its augmenting entropy condition is introduced. In this context it is concluded that consistency and L(infinity)-stability guarantee for a relevant class of admissible entropy functions, that their entropy production rate belongs to a compact subset of H(loc)sup -1 (x,t). One can now use compensated compactness arguments in order to turn this conclusion into a convergence proof. The current state of the art for these arguments includes the scalar and a wide class of 2 x 2 systems of conservation laws. The general framework of the vanishing viscosity method is studied as an effective way to meet the consistency and L(infinity)-stability requirements. How this method is utilized to enforce consistency and stability for scalar conservation laws is shown. In this context we prove, under the appropriate assumptions, the convergence of finite difference approximations (e.g., the high resolution TVD and UNO methods), finite element approximations (e.g., the Streamline-Diffusion methods) and spectral and pseudospectral approximations (e.g., the Spectral Viscosity methods).

Tadmor, Eitan↗

NASA NTRS1992

Stochastic modeling of turbulent reacting flows

Direct numerical simulations of a single-step irreversible chemical reaction with non-premixed reactants in forced isotropic turbulence at R(sub lambda) = 63, Da = 4.0, and Sc = 0.7 were made using 128 Fourier modes to obtain joint probability density functions (pdfs) and other statistical information to parameterize and test a Fokker-Planck turbulent mixing model. Preliminary results indicate that the modeled gradient stretching term for an inert scalar is independent of the initial conditions of the scalar field. The conditional pdf of scalar gradient magnitudes is found to be a function of the scalar until the reaction is largely completed. Alignment of concentration gradients with local strain rate and other features of the flow were also investigated.

Fox, R. O.↗

NASA NTRS2001

Responses of Precipitation and Hydrologic Processes to Tropical SST

The goal of the research is to identify the mechanisms in the response of tropical precipitation and atmospheric hydrologic cycle to sea surface temperature (SST) variability at seasonal-to-interannual time scales, and to utilize the knowledge for better understanding of climate feedback processes relevant to global change. As a first step to achieve the goal, we characterize the inter-relationship among convective/stratiform rain, ice/water clouds water vapor, and SST using TRMM satellite data and a cloud-resolving model. We examined the daily hydrologic variables [column water vapor (PW), cloud liquid water (CW), rainfall rates (RR)] as a function of SST using high-resolution data (0.25 x 0.25, daily) derived from TRMM satellite measurements. Comparing the winter of 97/98 (El Nino condition) against the winter of 99/00 (La Nina condition), area-mean values of all four hydrologic variables in cloudy areas within the tropical Pacific are higher in the El Nino winter than in the La Nina winter. This is consistent with previous observational analyses and SST warming experiments (idealized or ENSO-like) that showed the interaction between hydrologic cyclic and radiation at the seasonal to interannual time scales leads to intensified tropical circulation and hydrologic cycle. However, there is evidence that the enhanced hydrologic cycle over the warm pool is accompanied by an expansion of radiatively -driven subsidence in response to a stronger SST gradient between warm pool and surrounding cold pool. The expanding subsidence effectively reduces cloud amounts over the warm pool. Our analysis of daily variability further indicates a more vigorous water cycle characterized by higher PW, CW, and RR in response to overall warming. This is expected from the Clausius Clapeyron relation as a thermodynamic response to warming. However cloudy areas decrease in response to overall warming. This may be due to factors that are fundamentally different. One possibility is that in a warmer climate, there may be alteration in the microphysical processes in clouds, e.g., increased coalescence of cloud size droplets and increased entrainment of dry air from above due to the invigorated convective motions in response to overall warming. As a result, cloud amounts are less. Another interpretation is that, in colder air masses, more ice particles will be produced and, because of their lower densities, are less likely to precipitate, hence the higher cloud amount. This suggests that, when averaged over a sufficiently large area covering both the rising and sinking branches of the tropical circulation, there maybe an increase in the precipitation efficiency in stratiform clouds, in conjunction with increased subsidence in an enhanced water cycle. This is supported by the residence time at different rainrates.

Sui, C.-H.↗