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 343 records · Page 19

Investigating the Durability of Iodine Waste Forms in Dilute Conditions

To prevent the release of radioiodine during the reprocessing of used nuclear fuel or in the management of other wastes, many technologies have been developed for iodine capture. The capture is only part of the challenge as a durable waste form is required to ensure safe disposal of the radioiodine. This work presents the first durability studies in dilute conditions of two AgI-containing waste forms: hot-isostatically pressed silver mordenite (AgZ) and spark plasma sintered silver-functionalized silica aerogel (SFA) iodine waste forms (IWF). Using the single-pass flow-through (SPFT) test method, the dissolution rates respective to Si, Al, Ag and I were measured for variants of the IWFs. By combining solution and solid analysis information on the corrosion mechanism neutral-to-alkaline conditions was elucidated. The AgZ samples were observed to have corrosion preferentially occur at secondary phases with higher Al and alkali content. These phases contained a lower proportion of I compared with the matrix. The SFA samples experienced a higher extent of corrosion at Si-rich particles, but an increased addition of Si to the waste led to an improvement in corrosion resistance. The dissolution rates for the IWF types are of similar magnitude to other Si-based waste form materials measured using SPFT.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DFT and microkinetic comparison of Pt, Pd and Rh-catalyzed ammonia oxidation

Ammonia oxidation is the heart of the Ostwald process and is important in emissions control. Catalytic behaviors are a function of conditions and are observed to vary across the platinum group metals (PGMs) Pt, Pd, Rh. Here, we combine density functional theory computations and microkinetic modeling to rationalize these dependencies. We compute reactions over model (2 1 1) and (1 1 1) surfaces of PGMs. Binding energies are similar on Pd and Pt and generally greater on Rh, while activation energies vary across all metals. Rates on (2 1 1) surfaces are greater than (1 1 1) surfaces. Additionally, the stepped Pt is most active and stepped Rh most selective to N 2 at ammonia slip conditions, while at Ostwald process conditions, stepped Pd is most active and stepped Pt most selective to NO. Degree of rate and selectivity control analysis provides insights into the reactions limiting performance of PGMs. Both activation barriers and surface coverages influence rates and selectivities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Mechanical creep testing and characterization of Cr coated Zircaloy

Zirconium based cladding has been developed for use in light water reactor fuel applications due to the material’s desirable low neutron scattering cross section. To help improve corrosion resistance in accident conditions, such as a loss of coolant accident (LOCA), a chromium coating can be applied to the outer diameter of the Zircaloy tube. Cladding strength at reactor relevant temperatures is important for the structural performance. Stress relaxation creep testing was performed on Cr coated Zircaloy-4 to measure strain rate as a function of stress for 350 and 650°C. These were performed at both 5 ± 1 % and 7.5 ± 1 % strain. The lower temperature provides results near normal operating temperature, while the higher temperature provides results relevant to off-normal conditions where recrystallization will take place. Characterization was also performed to analyze the structure of the Cr coating and Zircaloy-4 as well as the interface between the two.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Absence of arterial baroreflex modulation of skin sympathetic activity and sweat rate during whole-body heating in humans

1. Prior findings suggest that baroreflexes are capable of modulating skin blood flow, but the effects of baroreceptor loading/unloading on sweating are less clear. Therefore, this project tested the hypothesis that pharmacologically induced alterations in arterial blood pressure in heated humans would lead to baroreflex-mediated changes in both skin sympathetic nerve activity (SSNA) and sweat rate. 2. In seven subjects mean arterial blood pressure was lowered (approximately 8 mmHg) and then raised (approximately 13 mmHg) by bolus injections of sodium nitroprusside and phenylephrine, respectively. Moreover, in a separate protocol, arterial blood pressure was reduced via steady-state administration of sodium nitroprusside. In both normothermia and heat-stress conditions the following responses were monitored: sublingual and mean skin temperatures, heart rate, beat-by-beat blood pressure, skin blood flow (laser-Doppler flowmetry), local sweat rate and SSNA (microneurography from peroneal nerve). 3. Whole-body heating increased skin and sublingual temperatures, heart rate, cutaneous blood flow, sweat rate and SSNA, but did not change arterial blood pressure. Heart rate was significantly elevated (from 74 +/- 3 to 92 +/- 4 beats x min(-1); P < 0.001) during bolus sodium nitroprusside-induced reductions in blood pressure, and significantly reduced (from 92 +/- 4 to 68 +/- 4 beats x min(-1); P < 0.001) during bolus phenylephrine-induced elevations in blood pressure, thereby demonstrating normal baroreflex function in these subjects. 4. Neither SSNA nor sweat rate was altered by rapid (bolus infusion) or sustained (steady-state infusion) changes in blood pressure regardless of the thermal condition. 5. These data suggest that SSNA and sweat rate are not modulated by arterial baroreflexes in normothermic or moderately heated individuals.

Clinical Trial↗

Proton Dynamics Scenarios in the Integrable Optics Test Accelerator (IOTA) at Fermilab

The Integrable Optics Test Accelerator (IOTA) at Fermilab provides a versatile platform for studying the interplay of space-charge, impedance, and non-linear optics in high-intensity hadron beams within synchrotrons and storage rings. This report examines the parameters and dynamics of 2.5~MeV proton beam operations in two configurations of the bare IOTA lattice: one for demonstrating Non-linear Integrable Optics with the Danilov-Nagaitsev magnet, and the other for use with electron cooling. We offer order-of-magnitude estimates of the transverse emittance growth rate as a function of beam intensity, highlighting contributions from residual gas scattering, intra-beam scattering, and space-charge effects. Under nominal conditions, the beam lifetime is projected to be less than 7~minutes at low intensity with the current vacuum quality, and fewer than 100,000~turns at high intensity due to strong space-charge effects. The calculations presented here will guide strategies to mitigate emittance growth and inform future IOTA experiments.

43 PARTICLE ACCELERATORS↗

Control of thermal balance by a liquid circulating garment based on a mathematical representation of the human thermoregulatory system

Test data and a mathematical model of the human thermoregulatory system were used to investigate control of thermal balance by means of a liquid circulating garment (LCG). The test data were derived from five series of experiments in which environmental and metabolic conditions were varied parametrically as a function of several independent variables, including LCG flowrate, LCG inlet temperature, net environmental heat exchange, surrounding gas ventilation rate, ambient pressure, metabolic rate, and subjective/obligatory cooling control. The resultant data were used to relate skin temperature to LCG water temperature and flowrate, to assess a thermal comfort band, to demonstrate the relationship between metabolic rate and LCG heat dissipation, and so forth. The usefulness of the mathematical model as a tool for data interpretation and for generation of trends and relationships among the various physiological parameters was also investigated and verified.

Kuznetz, L. H.↗

Approach guidance logic for a tilt-rotor aircraft

The distinctive feature of a tilt-rotor aircraft is that the pilot can change the rotor mast angles to go from a helicopter configuration for take-off and landing to an airplane configuration for high cruise speeds and vice-versa. An approach path for such an aircraft is proposed and the logic required to fly along this path in the presence of wind is determined. The main contribution of this work is an efficient and, to my knowledge, new method for generating the nominal state and control histories taking into account an estimate of the mean wind velocity and direction. The method requires the solution of algebraic (mostly linear) equations to generate a 'universal nominal', and feedforward and feedback gains. Then, in flight the additional state and control corrections due to deviation in descent rate, deceleration, and flight in a steady wind are obtained by multiplying simple precalculated functions of time by descent rate, deceleration or sine and cosine components of the mean wind vector. Simulations of approach flights for different wind conditions, assuming perfect state information in the feedback signal, indicated satisfactory performance.

Beser, J.↗

Optimal control of the Shuttle-Tethered-Subsatellite system

Alternate control laws based on an application of the linear regulator problem are developed for possible use in the operation of the Shuttle-Tethered-Subsatellite system. Control is assumed to be provided only by modulating the tension level in the tether as a function of the difference between actual and commanded tether line length, length rate, in (orbital) plane swing angle and swing angle rate. Necessary and sufficient conditions for stability of the (linear) system motion in the vicinity of its nominal local vertical orientation are developed. By proper selection of the state and control penalty matrices it is possible to obtain faster responses with no increase in maximum power levels for use in station keeping when compared with alternate control strategies. The weighting matrices are adjusted in a piecewise adaptive manner to provide control law gains in order to achieve a smooth deployment history. Successful retrieval is dependent mainly on the initial conditions and the rate at which the commanded length is reduced.

Bainum, P. M.↗

Axisymmetric confined turbulent jet directed towards the liquid surface from below

A numerical simulation is presented of an axisymmetric turbulent jet discharging axially from below into a cylindrical tank and directed towards the liquid vapor interface. The liquid vapor interface is assumed to be flat and shear free. The k-epsilon turbulence model is used to calculate the eddy viscosity. The turbulence intensity distribution and the length scale associated with the k-epsilon model are calculated as functions of jet flow rates and systems parameters. Numerical results are compared with appropriate experimental data. The problems associated with the free surface boundary conditions for the turbulent quantities are discussed.

Hasan, Mohammad M.↗

Axisymmetric confined turbulent jet directed towards the liquid surface from below

A numerical simulation is presented of an axisymmetric turbulent jet discharging axially from below into a cylindrical tank and directed towards the liquid vapor interface. The liquid vapor interface is assumed to be flat and shear free. The k-epsilon turbulence model is used to calculate the eddy viscosity. The turbulence intensity distribution and the length scale associated with the k-epsilon model are calculated as functions of jet flow rates and systems parameters. Numerical results are compared with appropriate experimental data. The problems associated with the free surface boundary conditions for the turbulent quantities are discussed.

Hasan, Mohammad M.↗

A general kinetic framework for dislocation mobility derived from probabilistic cellular automaton simulations of the kink-pair mechanism

Dislocation mobility laws are essential components of dislocation-density-based crystal plasticity models. For dislocations governed by the kink-pair mechanism, however, existing formulations are typically limited to specific regimes due to the com plex interplay between stochastic kink-pair nucleation and lateral kink migration. In this work, we develop a general kinetic framework that expresses the average dislocation velocity as a function of mechanism-level variables: positive/negative kink pair nucleation rates, kink migration velocity, dislocation segment length, critical kink-pair width, and kink height. Probabilis tic cellular automaton simulations are used to capture the behavior of conceptual dislocation segments between the limiting conditions of migration outpacing nucleation on the one end and nucleation outpacing migration on the other. An elemen tary functional form that captures the system dynamics is suggested and fitted against the simulation results. This framework remains valid for arbitrary combinations of the six variables and is, therefore, compatible with any admissible constitutive re lations that describe their stress and temperature dependence. Comparisons with established approaches and experimental results confirm the robustness and physical consistency of the formulation, making it broadly applicable to material systems in which dislocation motion is governed by the kink-pair mechanism.

36 MATERIALS SCIENCE↗

Numerical identification of boundary conditions on nonlinearly radiating inverse heat conduction problems

An explicit and unconditionally stable finite difference method for the solution of the transient inverse heat conduction problem in a semi-infinite or finite slab mediums subject to nonlinear radiation boundary conditions is presented. After measuring two interior temperature histories, the mollification method is used to determine the surface transient heat source if the energy radiation law is known. Alternatively, if the active surface is heated by a source at a rate proportional to a given function, the nonlinear surface radiation law is then recovered as a function of the interface temperature when the problem is feasible. Two typical examples corresponding to Newton cooling law and Stefan-Boltzmann radiation law respectively are illustrated. In all cases, the method predicts the surface conditions with an accuracy suitable for many practical purposes.

Murio, Diego A.↗

Computationally efficient and error aware surrogate construction for numerical solutions of subsurface flow through porous media

Limiting the injection rate to restrict the pressure below a threshold at a critical location can be an important goal of simulations that model the subsurface pressure between injection and extraction wells. The pressure is approximated by the solution of Darcy’s partial differential equation for a given permeability field. The subsurface permeability is modeled as a random field since it is known only up to statistical properties. This induces uncertainty in the computed pressure. Solving the partial differential equation for an ensemble of random permeability simulations enables estimating a probability distribution for the pressure at the critical location. These simulations are computationally expensive, and practitioners often need rapid online guidance for real-time pressure management. An ensemble of numerical partial differential equation solutions is used to construct a Gaussian process regression model that can quickly predict the pressure at the critical location as a function of the extraction rate and permeability realization. The Gaussian process surrogate analyzes the ensemble of numerical pressure solutions at the critical location as noisy observations of the true pressure solution, enabling robust inference using the conditional Gaussian process distribution. Our first novel contribution is to identify a sampling methodology for the random environment and matching kernel technology for which fitting the Gaussian process regression model scales as O ( n log n ) instead of the typical O ( n 3 ) rate in the number of samples n used to fit the surrogate. The surrogate model allows almost instantaneous predictions for the pressure at the critical location as a function of the extraction rate and permeability realization. Our second contribution is a novel algorithm to calibrate the uncertainty in the surrogate model to the discrepancy between the true pressure solution of Darcy’s equation and the numerical solution. Finally, although our method is derived for building a surrogate for the solution of Darcy’s equation with a random permeability field, the framework broadly applies to solutions of other partial differential equations with random coefficients.

54 ENVIRONMENTAL SCIENCES↗

Growth rate calculations of auroral kilometric radiation using the relativistic resonance condition

The relativistic cyclotron resonance condition for right-handed extraordinary mode waves defines an ellipse in velocity space. The position of the center and size of the semiminor axis of this ellipse are functions of the plasma frequency, gyrofrequency, wave frequency, and wave normal angle. The effect of varying these parameters on the position and size of the resonance contour is analyzed. The results show that as the wave normal angle decreases, the semiminor axis increases in size and as the plasma frequency to gyrofrequency ratio decreases, the minimum energy for resonating electrons decreases and the maximum wave normal angle allowed by the resonance condition increases. Also, as the wave frequency to gyrofrequency ratio increases, the center of the resonance ellipse moves away from the origin. The relativistic resonance condition and the electron distribution in velocity space obtained by the S3-3 satellite are used to calculate numerically growth rates for the terrestrial auroral kilometric radiation. It is shown that the loss cone region of the electron distribution can give rise to growth rates for the extraordinary mode that are sufficiently large to account for the observed radio emission intensities.

Omidi, N.↗

Aerothermal analysis of a wing-elevon cove with variable leakage

A mathematical model of the heating associated with leakage within the space shuttle orbiter wing elevon cove was used to analyze the aerothermal response of the structure. The analysis was applied to wind tunnel results and to shuttle entry conditions. A parametric study described the mechanism by which the energy of the ingested mass is distributed as a function of the Reynolds number, wall conductivity, surface emissivity, and cove height and length. The rate at which energy was transferred into the cove interior was primarily determined by the combined effects of the convection, wall capacitance, and internal radiation terms. Correlation with wind tunnel results indicated that the cold wall convective heating rates, wall temperature, and gas temperature of the cove were predicted by the analytical model for short exposure times. Predicted transient thermal response of the elevon cove subjected to shuttle entry conditions indicated that the allowable leak area at the cove seal was 20 percent less than that previously indicated by experiment.

Hunt, R. L.↗

Retention of Cr in Meta-kaolin-based geo-polymeric Materials under Oxidation and Carbonation - 17453

Alkali activated alumino-silicates, often referred to as geo-polymers, are considered promising immobilization matrices for contaminants due to their high sorption capacity, durability and low leachability. For some contaminants they may prove superior to cement pastes that are conventional immobilization materials in waste management applications. However, estimating the retention capacity of some contaminants under environmental conditions that include oxidation and carbonation and may induce a valance change, is complicated. These contaminants tend to be soluble under oxidizing condition and form precipitates under reducing conditions. Thus, contaminant retention is highly dependent on the initial oxidation state at which the metal was introduced, as well as the matrix sorption ability, ability to resist oxidation and carbonation, and the contaminant leachability. Moreover, carbonation in considered one of the most important degradation mechanisms in cementitious material, causing porosity reduction and significant changes in contaminant leachability, however, geo-polymers have considerably lower Ca content relative to cements, or none, and the effect of carbonation on these pastes is less clear. Thus, contaminant release mechanisms and rates must be estimated based upon exposure to various atmospheres appropriate to waste management applications. In this study we aim to estimate the sorption capacity, and leachability of Cr from a meta-kaolin based geo-polymer, doped with Cr, as a function of its oxidation state and carbonation. Cr has a high valance range from +6 to 0 similar to Tc, and is used as its simulant in the current study. Sodium dichromate was used as an oxidized source metal material. To induce reducing conditions, ascorbic acid was added except for experimental control cases. The samples were aged under N{sub 2}, 98% N{sub 2} + 2% CO{sub 2}, CO{sub 2} deficient air and air environments under constant relative humidity. Following aging, samples were subject to equilibrium leaching as a function of pH (EPA Method 1313) and monolith mass transport rate testing (EPA Method 1315). Mass balance calculations, based on the results of a pH dependent leaching test, have shown that oxidized Cr samples lost 85-95% of their Cr content, whereas, the reduced Cr samples lost only 20% of their Cr content. On the other hand, the leachability of Si and Al is much higher in the reduced-Cr samples then in the samples with oxidized Cr, and their natural pH is up to 2.5 pH units lower than that of the oxidizing cases. Thus, we conclude that the combination of the reducing agent and the spike chemistry influence the matrix formation. Sample aging have not resulted in any significant change in the Cr-leaching behavior of the oxidizing matrices in the pH dependent leaching test, although a significant decrease in the natural pH, of 1.5 pH unit, was observed upon exposure to elevated CO{sub 2} levels (2%), probably due to the formation of Na-carbonate species. The pH and leaching behavior of the reduced Cr samples under the various aging conditions have not yielded significant variation either in their natural pH or in their leaching behavior. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Convergence behavior of the random phase approximation renormalized correlation energy

Based on the random phase approximation (RPA), RPA renormalization is a robust many-body perturbation theory that works for molecules and materials because it does not diverge as the Kohn-Sham gap approaches zero. Additionally, RPA renormalization enables the simultaneous calculation of RPA and beyond-RPA correlation energies since the total correlation energy is the sum of a series of independent contributions. The first-order approximation (RPAr1) yields the dominant beyond-RPA contribution to the correlation energy for a given exchange-correlation kernel, but systematically underestimates the total beyond-RPA correction. For both the homogeneous electron gas model and real systems, we demonstrate numerically that RPA renormalization beyond first order converges monotonically to the infinite-order beyond-RPA correlation energy for several model exchange-correlation kernels and that the rate of convergence is principally determined by the choice of the kernel and spin polarization of the ground state. The monotonic convergence is rationalized from an analysis of the RPA renormalized correlation energy corrections, assuming the exchange-correlation kernel and response functions satisfy some reasonable conditions. In this work, for spin-unpolarized atoms, molecules, and bulk solids, we find that RPA renormalization is typically converged to 1 meV error or less by fourth order regardless of the band gap or dimensionality. Most spin-polarized systems converge at a slightly slower rate, with errors on the order of 10 meV at fourth order and typically requiring up to sixth order to reach 1 meV error or less. Slowest to converge, however, open-shell atoms present the most challenging case and require many higher orders to converge.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Evaluation of a locally homogeneous model of spray evaporation

Measurements were conducted on an evaporating spray in a stagnant environment. The spray was formed using an air-atomizing injector to yield a Sauter mean diameter of the order of 30 microns. The region where evaporation occurred extended approximately 1 m from the injector for the test conditions. Profiles of mean velocity, temperature, composition, and drop size distribution, as well as velocity fluctuations and Reynolds stress, were measured. The results are compared with a locally homogeneous two-phase flow model which implies no velocity difference and thermodynamic equilibrium between the phases. The flow was represented by a k-epsilon-g turbulence model employing a clipped Gaussian probability density function for mixture fraction fluctuations. The model provides a good representation of earlier single-phase jet measurements, but generally overestimates the rate of development of the spray. Using the model predictions to represent conditions along the centerline of the spray, drop life-history calculations were conducted which indicate that these discrepancies are due to slip and loss of thermodynamic equilibrium between the phases.

Shearer, A. J.↗