Engineering Papers⌕ Search

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.

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

At least 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Self-pressurization of a flightweight liquid hydrogen tank: Effects of fill level at low wall heat flux

Experimental results are presented for the self pressurization and thermal stratification of a 4.89 cu m liquid hydrogen storage tank subjected to low heat flux (2.0 and 3.5 W/sq m) in normal gravity. The test tank was representative of future spacecraft tankage, having a low mass to volume ratio and high performance multilayer thermal insulation. Tests were performed at fill levels of 29 and 49 pcts. (by volume) and complement previous tests at 83 pct. fill. As the heat flux increases, the pressure rise rate at each fill level exceeds the homogeneous rate by an increasing ratio. Herein, this ratio did not exceed a value of 2. The slowest pressure rise rate was observed for the 49 pct. fill level at both heat fluxes. This result is attributed to the oblate spheroidal tank geometry which introduces the variables of wetted wall area, liquid-vapor interfacial area, and ratio of side wall to bottom heating as a function of fill level or liquid depth. Initial tank thermal conditions were found to affect the initial pressure rise rate. Quasi steady pressure rise rates are independent of starting conditions.

Vandresar, N. T.↗
NASA NTRS1992

Self-pressurization of a flightweight liquid hydrogen tank - Effects of fill level at low wall heat flux

Experimental results are presented for the self pressurization and thermal stratification of a 4.89 cu m liquid hydrogen storage tank subjected to low heat flux (2.0 and 3.5 W/sq m) in normal gravity. The test tank was representative of future spacecraft tankage, having a low mass to volume ratio and high performance multilayer thermal insulation. Tests were performed at fill levels of 29 and 49 pcts. (by volume) and complement previous tests at 83 pct. fill. As the heat flux increases, the pressure rise rate at each fill level exceeds the homogeneous rate by an increasing ratio. Herein, this ratio did not exceed a value of 2. The slowest pressure rise rate was observed for the 49 pct. fill level at both heat fluxes. This result is attributed to the oblate spheroidal tank geometry which introduces the variables of wetted wall area, liquid-vapor interfacial area, and ratio of side wall to bottom heating as a function of fill level or liquid depth. Initial tank thermal conditions were found to affect the initial pressure rise rate. Quasi steady pressure rise rates are independent of starting conditions.

Van Dresar, N. T.↗
DOE OSTI2018

A Physical Model Study of Two-Phase Gas–Liquid Flows in a Ladle Shroud

Argon-steel flows inside a ladle shroud during teeming from a ladle to a tundish have been modelled physically. To this end, full-scale Perspex models of bloom as well as slab casting shrouds (BCS and SCS), operating with air and water, have been applied. Both open to air as well as immersed conditions were investigated with and without gas injection. Flows inside a ladle shroud under open to air and immersed conditions were found to be substantially different with a strong function of gas and liquid flow rates, collector nozzle and shroud diameters. Depending on the volumetric gas injection rate relative to liquid flow rate, different flow regimes have been observed in an immersed shroud [i.e., 0<(((d s )/(L s )))≤0.24]. At extremely low gas flow rates, [i.e., (((Q g )/(Q L )))≤0.02], injected gas is completely entrained as bubbles by the down-flowing liquid resulting in a bubbly two-phase flow over the entire length of a shroud. However, with an increasing gas flow rate, two distinctly different regions start to develop within the shroud body: a free liquid jet in the upper part and a gas–liquid mixing zone below. The length of the free jet increases with an increasing gas flow rate and at significantly higher gas to liquid flow rates [viz., (((Q g )/(Q L ))) BCS ≥0.42] and [viz., (((Q g )/(Q L ))) SCS ≥0.30], and the free jet is found to prevail over the entire length of the shroud. Within the range of conditions studied, it is observed that the free jet length or the line of demarcation between the jetting and two-phase mixing zone depends on gas and liquid flow rates and is specific to a particular shroud-collector nozzle system. Physical model results further indicate that a sufficiently large free jet length (~ shroud length) tends to create a high pressure region inside a shroud and prevent ingression of air. Possible implications of the present findings with reference to industrial teeming practices are also discussed in the text.

42 ENGINEERING↗
DOE OSTI2017

Flame thickness and conditional scalar dissipation rate in a premixed temporal turbulent reacting jet

The flame structure corresponding to lean hydrogen–air premixed flames in intense sheared turbulence in the thin reaction zone regime is quantified from flame thickness and conditional scalar dissipation rate statistics, obtained from recent direct numerical simulation data of premixed temporally-evolving turbulent slot jet flames. It is found that, on average, these sheared turbulent flames are thinner than their corresponding planar laminar flames. Extensive analysis is performed to identify the reason for this counter-intuitive thinning effect. The factors controlling the flame thickness are analyzed through two different routes i.e., the kinematic route, and the transport and chemical kinetics route. The kinematic route is examined by comparing the statistics of the normal strain rate due to fluid motion with the statistics of the normal strain rate due to varying flame displacement speed or self-propagation. It is found that while the fluid normal straining is positive and tends to separate iso-scalar surfaces, the dominating normal strain rate due to self-propagation is negative and tends to bring the iso-scalar surfaces closer resulting in overall thinning of the flame. The transport and chemical kinetics route is examined by studying the non-unity Lewis number effect on the premixed flames. The effects from the kinematic route are found to couple with the transport and chemical kinetics route. In addition, the intermittency of the conditional scalar dissipation rate is also examined. It is found to exhibit a unique non-monotonicity of the exponent of the stretched exponential function, conventionally used to describe probability density function tails of such variables. As a result, the non-monotonicity is attributed to the detailed chemical structure of hydrogen-air flames in which heat release occurs close to the unburnt reactants at near free-stream temperatures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1998

Human heart rate variability relation is unchanged during motion sickness

In a study of 18 human subjects, we applied a new technique, estimation of the transfer function between instantaneous lung volume (ILV) and instantaneous heart rate (HR), to assess autonomic activity during motion sickness. Two control recordings of ILV and electrocardiogram (ECG) were made prior to the development of motion sickness. During the first, subjects were seated motionless, and during the second they were seated rotating sinusoidally about an earth vertical axis. Subjects then wore prism goggles that reverse the left-right visual field and performed manual tasks until they developed moderate motion sickness. Finally, ILV and ECG were recorded while subjects maintained a relatively constant level of sickness by intermittent eye closure during rotation with the goggles. Based on analyses of ILV to HR transfer functions from the three conditions, we were unable to demonstrate a change in autonomic control of heart rate due to rotation alone or due to motion sickness. These findings do not support the notion that moderate motion sickness is manifested as a generalized autonomic response.

NASA Discipline Neuroscience↗
DOE OSTI2026

FunDiff: diffusion models over function spaces for physics-informed generative modeling

Recent advances in generative modeling-particularly diffusion models and flow matching-have been widely used for synthesizing discrete data such as images and videos. However, adapting these models to physical applications remains challenging, as the quantities of interest are continuous functions governed by complex physical laws. To address this, we introduce FunDiff, an efficient and robust framework for generative modeling in function spaces. FunDiff combines a latent diffusion process with a function autoencoder architecture to handle input functions with varying discretizations, generates continuous functions that can be evaluated at arbitrary locations, and seamlessly incorporate physical priors. These priors are enforced through architectural constraints or physics-informed loss functions, ensuring that generated samples satisfy fundamental physical laws. We theoretically establish minimax optimality guarantees for density estimation in function spaces, demonstrating that diffusion-based estimators achieve optimal convergence rates under suitable regularity conditions. We further demonstrate the practical effectiveness of FunDiff across diverse applications in fluid dynamics and solid mechanics. Empirical results indicate that our method can generate physically consistent samples with high fidelity to the target distribution, and exhibit robustness to noisy and low-resolution data.

Wang, Sifan [Yale University, New Haven, CT (Unite↗
DOE OSTI2021

Predicting nonresonant pressure-driven MHD modes in equilibria with low magnetic shear

Nonresonant internal modes can be difficult to anticipate as there is no resonant surface in the plasma. However, equilibria that are unstable to multiple nonresonant magnetohydrodynamic (MHD) modes may be more prone to global loss of confinement since these instabilities generate spatially extended linear displacements, potentially enhancing magnetic field line chaos via nonlinear interactions. Here, we successfully predict the unstable nonresonant pressure-driven modes for equilibria with zero shear in the plasma core, irrational q on axis, and a central pressure gradient, which is consistent with pre-crash profiles in sawtoothing tokamak plasmas in the large-aspect-ratio limit. A criterion for identifying nonresonant modes most likely to be unstable is developed from the convergents of the continued fraction representation of q 0 . A higher-order analysis of the standard Energy Principle reveals the conditions under which these modes are expected to dominate. Linear growth rate spectra, as a function of toroidal mode number (up to n = 30), calculated using the initial-value extended-MHD code, M3D-C1, recover the characteristic dependence observed for ideal infernal modes. Nonresonant modes have also been invoked in some ideal sawtooth crash models. Finally, this work provides a mechanism to predict the mode numbers of infernal modes and, potentially, the width of some post-sawtooth-crash profiles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Reconnection and interchange instability in the near magnetotail

This paper provides insights into the possible coupling between reconnection and interchange/ballooning in the magnetotail related to substorms and flow bursts. The results presented are largely based on recent simulations of magnetotail dynamics, exploring onset and progression of reconnection. 2.5-dimensional particle-in-cell (PIC) simulations with different tail deformation demonstrate a clear boundary between stable and unstable cases depending on the amount of deformation, explored up to the real proton/electron mass ratio. The evolution prior to onset, as well as the evolution of stable cases, are governed by the conservation of integral flux tube entropy S as imposed in ideal MHD, maintaining a monotonic increase with distance downtail. This suggests that ballooning instability in the tail should not be expected prior to the onset of tearing and reconnection. 3-D MHD simulations confirm this conclusion, showing no indication of ballooning prior to reconnection, if the initial state is ballooning stable. The simulation also shows that, after imposing resistivity necessary to initiate reconnection, the reconnection rate and energy release initially remain slow. However, when S becomes reduced from plasmoid ejection and lobe reconnection, forming a negative slope in S as a function of distance from Earth, the reconnection rate and energy release increase drastically. The latter condition has been shown to be necessary for ballooning/interchange instability, and the cross-tail structures that develop subsequently in the MHD simulation are consistent with such modes. The simulations support a concept in which tail activity is initiated by tearing instability but significantly enhanced by the interaction with ballooning/interchange enabled by plasmoid loss and lobe reconnection.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Analysis of Monolith Cores from an Engineering Scale Demonstration of a Prospective Cast Stone Process

The primary disposition path of Low Activity Waste (LAW) at the DOE Hanford Site is vitrification. A cementitious waste form is one of the alternatives being considered for the supplemental immobilization of the LAW that will not be treated by the primary vitrification facility. Washington River Protection Solutions (WRPS) has been directed to generate and collect data on cementitious or pozzolanic waste forms such as Cast Stone. This report documents the coring and leach testing of monolithic samples cored from an engineering-scale demonstration (ES Demo) with non-radioactive simulants. The ES Demo was performed at SRNL in October of 2013 using the Scaled Continuous Processing Facility (SCPF) to fill an 8.5 ft. diameter x 3.25 ft. high container with simulated Cast Stone grout. The Cast Stone formulation was chosen from the previous screening tests. Legacy salt solution from previous Hanford salt waste testing was adjusted to correspond to the average LAW composition generated from the Hanford Tank Waste Operation Simulator (HTWOS). The dry blend materials, ordinary portland cement (OPC), Class F fly ash, and ground granulated blast furnace slag (GGBFS or BFS), were obtained from Lafarge North America in Pasco, WA. In 2014 core samples originally obtained approximately six months after filling the ES Demo were tested along with bench scale molded samples that were collected during the original pour. A latter set of core samples were obtained in late March of 2015, eighteen months after completion of the original ES Demo. Core samples were obtained using a 2” diameter x 11” long coring bit. The ES Demo was sampled in three different regions consisting of an outer ring, a middle ring and an inner core zone. Cores from these three lateral zones were further segregated into upper, middle and lower vertical segments. Monolithic core samples were tested using the Environmental Protection Agency (EPA) Method 1315, which is designed to provide mass transfer rates (release rates) of inorganic analytes contained in monolithic material under diffusion controlled release conditions as a function of leaching time. Compressive strength measurements and drying tests were also performed on the 2015 samples. Leachability indices reported are based on analyte concentrations determined from dissolution of the dried samples.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS1994

The conductive propagation of nuclear flames. 2: Convectively bounded flames in C + O and O + Ne + Mg cores

We determine the speeds, and many other physical properties, of flame fronts that propagate inward into degenerate and semidegenerate cores of carbon and oxygen (CO) and neon and oxygen (NeOMg) white dwarfs when such flames are bounded on their exterior by a convective region. Combustion in such fronts, per se, is incomplete, with only a small part of the initial mass function burned. A condition of balanced power is set up in the star where the rate of energy emitted as neutrinos from the convective region equals the power available from the unburned fuel that crosses the burning front. The propagation of the burning front itself is in turn limited by the temperature at the base of the convective shell, while cannot greatly exceed the adiabatic value. Solving for consistency between these two conditions gives a unique speed for the flame. Typical values for CO white dwarfs are a few hundredths of a centimeter per second. Flames in NeOMg mixtures are slower. Tables are presented in a form that can easily be implemented in stellar evolution codes and yield the rate at which the convective shell advances into the interior. Combining these velocities with the local equations for stellar structure, we find a minimum density for each gravitational potential below with the local equations for stellar structure, we find a minimum density for each gravitational potential below which the flame cannot propagate, and must die. Although detailed stellar models will have to be constructed to reslove some issues conclusively, our results that a CO white dwarf inginted at its edge will not burn carbon all the way to its center unless the mass of the white dwarf exceeds 0.8 solar mass. On the other hand, it is difficult to ignite carbon burning by compression alone anywhere in a white dwarf whose mass does not exceed 1.0 solar mass. Thus, compressionally ignited shell carbon burning in an accerting CO dwarf almost certainly propagates all the way to the center of the star. Implications for neutron star formation, and Type Ia supernova models, are briefly discussed. These are also applicable to massive stars in the about 10-12 solar mass range which ignite neon burning off center.

Timmes, F. X.↗
DOE OSTI2021

Identification of Life Models for Li-Ion Batteries Using Penalized Regression and Bilevel Optimization

Reduced-order physics-based life models are extremely useful for rapidly predicting battery state-of-health and for simulating battery lifetime in arbitrary aging conditions. However, identification of well-parameterized models is difficult. This is because, for maximum usefulness in predicting lifetime under a variety of conditions, aging test data exhibits many degradation mechanisms, which all need to be accurately modeled. However, because aging tests are time-consuming and expensive, especially for large-format batteries, a minimum of tests are conducted while probing many stress factors. Building a well-parameterized model is then very challenging: an under-parameterized model will neglect critical degradation modes, and an over-parameterized model will extrapolate poorly to new testing conditions. To complicate this matter, the functional form of the model for any individual degradation rate can be very difficult to identify. In this work, the statistical tools of penalized regression and bilevel optimization are used to help identify both the functional forms of and optimize the parameters of reduced-order life models, accelerating identification of robust models. Model robustness is demonstrated through traditional statistics methods of cross-validation and sensitivity analysis, uncertainty quantification through bootstrap resampling and Monte-Carlo simulation, and simulation of real-world use cases.

47 OTHER INSTRUMENTATION↗
DOE OSTI2018

Rate theory scenarios study on fission gas behavior of U 3 Si 2 under LOCA conditions in LWRs

Fission gas behavior of U3Si2 under various loss-of-coolant accident (LOCA) conditions in light water reactors (LWRs) was simulated using rate theory. A rate theory model for U3Si2 that covers both steady-state operation and power transients was developed for the GRASS-SST code based on existing research reactor/ion irradiation experimental data and theoretical predictions of density functional theory (DFT) calculations. The steady-state and LOCA condition parameters were either directly provided or inspired by BISON simulations. Due to the absence of in-pile experiment data for U3Si2's fuel performance under LWR conditions at this stage of accident tolerant fuel (ATF) development, a variety of LOCA scenarios were taken into consideration to comprehensively and conservatively evaluate the fission gas behavior of U3Si2 during a LOCA.

Miao, Yinbin (ORCID:0000000231284275)↗
DOE OSTI2021

Consequences of adsorbate-adsorbate interactions for apparent kinetics of surface catalytic reactions

Lateral adsorbate interactions at catalyst surfaces are known to influence adsorption energies and reaction rates. Lattice-based kinetic Monte Carlo (kMC) simulations are able to capture these influences, but such models are typically parameterized for a specific reaction network and catalyst surface. Here we report kMC simulations to probe the influence of lateral interactions on simulated rates, rate orders, apparent activation energies, and Sabatier plots. We construct a simple, two-step reaction network involving a single adsorbate and rate- limiting diatomic dissociation, employ a generic repulsive lateral interaction model consistent with known adsorbate-adsorbate interactions on metal surfaces, and a rate model consistent with known Brønsted-Evans-Polyani relationships for diatomic dissocations. Furthermore, we juxtapose reaction kinetics over a wide range of reaction conditions and catalyst binding energies, as a function of interaction strength. We find that at a given zero-coverage binding energy and external conditions, adsorbate coverage decreases monotonically with increasing interaction strength, but absolute rates can vary linearly or nonlinearly. Interactions flatten the Sabatier volcano and shift the maximum towards stronger binding. Influences on apparent rate orders and activation energies are modest and are sensitive to interaction-induced adsorbate ordering. Model predictions are sensitive to lattice size effects at high coverages. The results, modelled for a simple reaction system, highlight the generic consequences of lateral interactions and guidance for identifying their signatures in observed kinetics.

kinetic Monte Carlo↗
DOE OSTI2015

Enhanced Dry Reforming of Methane on Ni and Ni-Pt Catalysts Synthesized by Atomic Layer Deposition

Atomic layer deposition (ALD) was used to deposit Ni and Pt on alumina supports to form monometallic and bimetallic catalysts with initial particle sizes of 1–2.4 nm. The ALD catalysts were more active (per mass of metal) than catalysts prepared by incipient wetness (IW) for dry reforming of methane (DRM), and they did not form carbon whiskers during reaction due to their sufficiently small size. Catalysts modified by Pt ALD had higher rates of reaction per mass of metal and inhibited coking, whereas NiPt catalysts synthesized by IW still formed carbon whiskers. Temperature-programmed reduction of Ni catalysts modified by Pt ALD indicated the presence of bimetallic interaction. Density functional theory calculations suggested that under reaction conditions, the NiPt surfaces form Ni-terminated surfaces that are associated with higher DRM rates (due to their C and O adsorption energies, as well as the CO formation and CH4 dissociation energies).

Environmental Molecular Sciences Laboratory↗
DOE OSTI2016

Optimizing conditions for utilization of an H 2 oxidation catalyst with outer coordination sphere functionalities

[Ni(PCy2NArginine2)2]2+ ( CyArg ) or [Ni(PCy2NBenzyl2)2]2+ ( CyBn ) were evaluated for H 2 oxidation as a function of temperature, pressure, and solvent. 70 °C and 100 atm H 2 result in a TOF of 1.1 × 106 s -1 and an overpotential of 240 mV for CyArg in water. In methanol the rates were 280 s -1 for CyArg and 80 s -1 for CyBn , demonstrating the importance of water and the outer coordination sphere ( OCS ).

bio-inspired↗
NASA NTRS1993

Kinematics and aerodynamics of the velocity vector roll

The velocity vector roll is an angular rotation of an aircraft about its instantaneous velocity vector, constrained to be performed at constant angle-of-attack (AOA), no sideslip, and constant velocity. Consideration of the aerodynamic force equations leads to requirements for body-axis yawing and pitching rotations that satisfy these constraints. Here, the body axis rotations, and the constraints, are used in the moment equations to determine the aerodynamic moments required to perform the velocity vector roll. For representative tactical aircraft, the conditions for maximum pitching moment are a function of orientation, occurring at about 90 deg of bank in a level trajectory. Maximum required pitching moment occurs at peak roll rate, and is achieved at AOA above 45 deg. The conditions for maximum rolling moment depend on the value of the roll mode time constant. For a small time constant (fast response) the maximum rolling moment occurs at maximum roll acceleration and zero AOA, largely independent of aircraft orientation; for a large time constant, maximum required rolling moment occurs at maximum roll rate, at maximum AOA, and at 180 deg of bank in level flight. Maximum yawing moment occurs at maximum roll acceleration, maximum AOA, and is largely independent of airplane orientation.

Durham, Wayne C.↗
DOE OSTI2019

Methodology to estimate vapor quality of refrigerant leaks

The US Department of Energy (DOE) Building Technology Office (BTO) and the Air-conditioning, Heating & Refrigeration Technology Institute (AHRTI) are sponsoring an experimental study of refrigerant leak characteristics. The objectives are to conduct refrigerant leak tests on several heating, ventilating, air conditioning and refrigeration (HVAC&R) systems under operating conditions representative of actual applications to document the pressure decay rate and the mass flow rate of leaked refrigerant as a function of time. Another aspect of the study was to formulate and document a procedure for estimating the quality (vapor versus liquid) of the leaking refrigerant and oil escaping from the high- and low-pressure sides with the unit energized and with the unit turned off. Refrigerant quality was estimated by tracking total system weight of a packaged air-conditioning (A/C) unit and the total weight of a large container used to capture the refrigerant leakage escaping the packaged unit.

Miller, William A.↗