Search NASA⌕ Search

SEARCH · Search NASA

Results for “Complete linearization method”

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.

503 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Development of an All-Aqueous Thermally Regenerative Redox Flow Battery to Support Fossil Fuel Assets

Low-temperature thermal energy, a largely untapped resource, holds significant promise for large-scale electrical power generation globally. Various stationary sources, including industrial entities and thermal power plants, emit considerable low-temperature heat that currently remains unutilized. This energy is often overlooked because its low temperature renders it unsuitable for efficient power generation using conventional methods. However, current research is exploring diverse technologies capable of converting low-temperature heat into grid-scale power, aiming to enhance grid efficiency, further decarbonization initiatives, and facilitate a shift toward more decentralized power systems. One such innovative technology is the thermally regenerative battery (TRB), noted for its high power and energy densities compared to similar technologies, positioning it as a potential game-changer in power generation. TRBs integrate two scalable and well-established unit operations: a redox flow battery and a distillation column. This integration suggests that once an effective TRB chemistry is established, the pathway to commercialization could be expedited. The copper-based thermally regenerative ammonia battery (Cu aq -TRAB) stands out as the first TRB that circumvents electrodeposition/dissolution reactions, stabilizing Cu(I) and Cu(II) within the electrolyte and maintaining stability of all electroactive species in an aqueous phase. This stabilization has led to improvements in coulombic efficiency, open circuit potential, and copper solubility, thereby enhancing power density, energy density, and overall energy efficiency. Preliminary tests were conducted to determine the effects of various electrolyte species on the performance metrics of the battery, both theoretically and experimentally. These tests revealed that the solubility of copper in the Cu aq -TRAB electrolyte was constrained by the Cu(I)-NH 3 complex. Adjusting the background electrolyte to 5 M NH4Br and the ligand concentration to 4 M NH 3 enabled the copper concentration to reach a maximum of 0.6 M. This modification led to an estimated theoretical maximum energy density of 9.5 Wh L -1 for the Cu aq -TRAB. Additionally, full cell testing indicated a tradeoff between peak power and energy density with varying copper and ammonia concentrations. Increasing the applied current density during discharge linearly raised the average power output, with a minimal reduction in energy density due to a balance between higher ohmic overpotential and reduced time for undesirable ammonia crossover. Furthermore, a comprehensive numerical sensitivity analysis of the complete Cu aq -TRAB system was performed. This analysis aimed to assess how the battery and the distillation column responded to changes in system input parameters, providing insights into optimal research directions for enhancing system performance. The analysis revealed that at room temperature, battery power was significantly more sensitive to ohmic losses than to mass transfer, with reaction rates having minimal impact. This trend continued even at higher temperatures. Also, the thermal energy required for ammonia separation was studied, showing that increased temperatures generally reduced energy requirements, except in low-pressure scenarios above 65 °C. An investigation into membrane performance in the Cu aq -TRAB was undertaken, given the significant impact of ammonia transport control and ohmic losses on system performance. Various membranes were evaluated to identify key performance metrics. Among the tested membranes, Selemion CMVN exhibited the highest performance, with a peak power density of 84 mW cm -2 and average values of 26 ± 6.8 mW cm -2 for power density and 2.9 Wh L -1 for energy density at an applied current density of 50 mA cm -2 . An economic assessment indicated a levelized cost of storage at $410 per MWh under optimal conditions, highlighting the commercial potential of the Cu aq -TRAB when utilizing cost-effective, readily available materials.

25 ENERGY STORAGE↗
DOE OSTI2022

A 1D–0D–3D coupled model for simulating blood flow and transport processes in breast tissue

In this work, we present mixed dimensional models for simulating blood flow and transport processes in breast tissue and the vascular tree supplying it. These processes are considered, to start from the aortic inlet to the capillaries and tissue of the breast. Large variations in biophysical properties and flow conditions exist in this system necessitating the use of different flow models for different geometries and flow regimes. In total, we consider four different model types. First, a system of 1D nonlinear hyperbolic partial differential equations (PDEs) is considered to simulate blood flow in larger arteries with highly elastic vessel walls. Second, we assign 1D linearized hyperbolic PDEs to model the smaller arteries with stiffer vessel walls. The third model type consists of ODE systems (0D models). It is used to model the arterioles and peripheral circulation. Finally, homogenized 3D porous media models are considered to simulate flow and transport in capillaries and tissue within the breast volume. Sink terms are used to account for the influence of the venous and lymphatic systems. Combining the four model types, we obtain two different 1D–0D–3D coupled models for simulating blood flow and transport processes: The first model results in a fully coupled 1D–0D–3D model covering the complete path from the aorta to the breast combining a generic arterial network with a patient specific breast network and geometry. The second model is a reduced one based on the separation of the generic and patient specific parts. The information from a calibrated fully coupled model is used as inflow condition for the patient specific sub-model allowing a significant computational cost reduction. Several numerical experiments are conducted to calibrate the generic model parameters and to demonstrate realistic flow simulations compared to existing data on blood flow in the human breast and vascular system. Moreover, we use two different breast vasculature and tissue data sets to illustrate the robustness of our reduced sub-model approach.

multiscale blood flow models↗
NASA NTRS2013

Robust Platinum Resistor Thermometer (PRT) Sensors and Reliable Bonding for Space Missions

Platinum resistance thermometers (PRTs) provide accurate temperature measurements over a wide temperature range and are used extensively on space missions due to their simplicity and linearity. A standard on spacecraft, PRTs are used to provide precision temperature control and vehicle health assessment. This paper reviews the extensive reliability testing of platinum resistor thermometer sensors (PRTs) and bonding methods used on the Mars Science Laboratory (MSL) mission and for the upcoming Soil Moisture Active Passive (SMAP) mission. During the Mars Exploration Rover (MER) mission, several key, JPL-packaged PRTs failed on those rovers prior to and within 1-Sol of landing due to thermally induced stresses. Similar failures can be traced back to other JPL missions dating back thirty years. As a result, MSL sought out a PRT more forgiving to the packaging configurations used at JPL, and extensively tested the Honeywell HRTS-5760-B-U-0-12 sensor to successfully demonstrate suitable robustness to thermal cycling. Specifically, this PRT was cycled 2,000 times, simulating three Martian winters and summers. The PRTs were bonded to six substrate materials (Aluminum 7050, treated Magnesium AZ231-B, Stainless Steel 304, Albemet, Titanium 6AL4V, and G-10), using four different aerospace adhesives--two epoxies and two silicones--that conformed to MSL's low out-gassing requirements. An additional epoxy was tested in a shorter environmental cycling test, when the need for a different temperature range adhesive was necessary for mobility and actuator hardware late in the fabrication process. All of this testing, along with electrostatic discharge (ESD) and destructive part analyses, demonstrate that this PRT is highly robust, and not subject to the failure of PRTs on previous missions. While there were two PRTs that failed during fabrication, to date there have been no in-flight PRT failures on MSL, including those on the Curiosity rover. Since MSL, the sensor has gone through a change in construction such that the manufacturer significantly restricts the minimum temperature. However, significant subsequent testing was performed with this new version of the part to show that it indeed is still robust to at least Mars minimum temperatures of -135 degrees Centigrade. The additional completed testing will be described. This work has resulted in a successful sensor package qualification and a reliable bonding method suitable for use over large temperature extremes

Platinum Resistor Thermometer (PRT)↗
NASA NTRS2013

Robust Platinum Resistor Thermometer (PRT) Sensors and Reliable Bonding for Space Missions

Platinum resistance thermometers (PRTs) provide accurate temperature measurements over a wide temperature range and are used extensively on space missions due to their simplicity and linearity. A standard on spacecraft, PRTs are used to provide precision temperature control and vehicle health assessment. This paper reviews the extensive reliability testing of platinum resistor thermometer sensors (PRTs) and bonding methods used on the Mars Science Laboratory (MSL) mission and for the upcoming Soil Moisture Active Passive (SMAP) mission. During the Mars Exploration Rover (MER) mission, several key, JPL-packaged PRTs failed on those rovers prior to and within 1-Sol of landing due to thermally induced stresses. Similar failures can be traced back to other JPL missions dating back thirty years. As a result, MSL sought out a PRT more forgiving to the packaging configurations used at JPL, and extensively tested the Honeywell HRTS-5760-B-U-0-12 sensor to successfully demonstrate suitable robustness to thermal cycling. Specifically, this PRT was cycled 2,000 times, simulating three Martian winters and summers. The PRTs were bonded to six substrate materials (Aluminum 7050, treated Magnesium AZ231-B, Stainless Steel 304, Albemet, Titanium 6AL4V, and G-10), using four different aerospace adhesives--two epoxies and two silicones--that conformed to MSL's low out-gassing requirements. An additional epoxy was tested in a shorter environmental cycling test, when the need for a different temperature range adhesive was necessary for mobility and actuator hardware late in the fabrication process. All of this testing, along with electrostatic discharge (ESD) and destructive part analyses, demonstrate that this PRT is highly robust, and not subject to the failure of PRTs on previous missions. While there were two PRTs that failed during fabrication, to date there have been no in-flight PRT failures on MSL, including those on the Curiosity rover. Since MSL, the sensor has gone through a change in construction such that the manufacturer significantly restricts the minimum temperature. However, significant subsequent testing was performed with this new version of the part to show that it indeed is still robust to at least Mars minimum temperatures of -135(sup o)C. The additional completed testing will be described. This work has resulted in a successful sensor package qualification and a reliable bonding method suitable for use over large temperature extremes.

dectector↗
NASA NTRS2004

Strain-Gauge Measurement of Weight of Fluid in a Tank

A method of determining the amount of fluid in a tank is based on measurement of strains induced in tank supports by the weight of the fluid. Unlike most prior methods, this method is nonintrusive: there is no need to insert instrumentation in the tank and, hence, no need to run wires, cables, or tubes through the tank wall. Also unlike most prior methods, this method is applicable even if the fluid in the tank is at supercritical pressure and temperature, because it does not depend on the presence of a liquid/gas interface (as in liquid-level-measuring methods). The strain gauges used in this method may be of two types: foil and fiber-optic. Four foil gauges (full bridge) are mounted on each of the tank-supporting legs. As the tank is filled or emptied, the deformation in each leg increases or decreases, respectively. Measured deformations of all legs are added to obtain a composite deformation indicative of the change in weight of the tank plus fluid. An initial calibration is performed by recording data at two points (usually, empty and full) for which the mass or weight of fluid is known. It is assumed that the deformations are elastic, so that the line passing through the two points can be used as a calibration curve of mass (or weight) of fluid versus deformation. One or more fiber-optic gauges may be used instead of the foil gauges. The resolution of the fiber-optic and foil gauges is approximately the same, but the fiber-optic gauges are immune to EMI (electromagnetic interference), are linear with respect to temperature over their entire dynamic range (as defined by the behavior of the sample), and measure thermally induced deformations as predictable signals. Conversely, long term testing has demonstrated that the foil gauges exhibit an erratic behavior whenever subjected to direct sun radiation (even if protected with a rubberized cover). Henceforth, for deployment in outdoor conditions, fiber-optic gauges are the only option if one is to rely on the system for an extended period of time when a recalibration procedure may not be acceptable. A set of foil gauges had been tested on the supports of a 500-gallon (1,900-liter) tank. The gauges were found to be capable of measuring the deformations (up to 22 micro-strain) that occurred during filling and emptying of the tank. The fluid masses calculated from the gauge readings were found to be accurate within 4.5 percent. However, the reliability of the foil gauges over a few hours was not acceptable. Therefore, the foil sensor system is acceptable for use only in controlled environments (complete shade, or indoors).

Figueroa, Jorge↗
NASA NTRS1994

Scaling, elasticity, and CLPT

The first few viewgraphs describe the general solution properties of linear elasticity theory which are given by the following two statements: (1) for stress B.C. on S(sub sigma) and zero displacement B.C. on S(sub u) the altered displacements u(sub i)(*) and the actual stresses tau(sub ij) are elastically dependent on Poisson's ratio nu alone: thus the actual displacements are given by u(sub i) = mu(exp -1)u(sub i)(*); and (2) for zero stress B.C. on S(sub sigma) and displacement B.C. on S(sub u) the actual displacements u(sub i) and the altered stresses tau(sub ij)(*) are elastically dependent on Poisson's ratio nu alone: thus the actual stresses are given by tau(sub ij) = E tau(sub ij)(*). The remaining viewgraphs describe the minimum parameter formulation of the general classical laminate theory plate problem as follows: The general CLT plate problem is expressed as a 3 x 3 system of differential equations in the displacements u, v, and w. The eighteen (six each) A(sub ij), B(sub ij), and D(sub ij) system coefficients are ply-weighted sums of the transformed reduced stiffnesses (bar-Q(sub ij))(sub k); the (bar-Q(sub ij))(sub k) in turn depend on six reduced stiffnesses (Q(sub ij))(sub k) and the material and geometry properties of the k(sup th) layer. This paper develops a method for redefining the system coefficients, the displacement components (u,v,w), and the position components (x,y) such that a minimum parameter formulation is possible. The pivotal steps in this method are (1) the reduction of (bar-Q(sub ij))(sub k) dependencies to just two constants Q(*) = (Q(12) + 2Q(66))/(Q(11)Q(22))(exp 1/2) and F(*) - (Q(22)/Q(11))(exp 1/2) in terms of ply-independent reference values Q(sub ij); (2) the reduction of the remaining portions of the A, B, and D coefficients to nondimensional ply-weighted sums (with 0 to 1 ranges) that are independent of Q(*) and F(*); and (3) the introduction of simple coordinate stretchings for u, v, w and x,y such that the process is neatly completed.

Brunelle, Eugene J.↗
NASA NTRS2016

Retrieval of Daily Maps of PM2.5 Aerosol in the Problematic California Valleys: Bright, Speckled Reflectances, Thin AOT, but High Pollution

The San Joaquin Valley suffers from severe episodes of respirable aerosol (PM2.5) in wintertime.We provide maps of aerosol episodes using daily snapshots of PM2.5 and its changing features despite numerous difficulties inherent to sampling the region. Linear relationships relating aerosol optical thickness (AOT) to PM2.5 give an explained variance of approximately 3.The GEO-CAPE mission has as a goal the provision of relevant measures of respirable aerosol to the community,but has not formulated a science goal beyond the limited goal of retrieval of AOT, bringing the usefulness of GEO-CAPE into doubt.Our special focus was on the DISCOVER-AQ period, Jan-Feb 2013, which had many supporting measurements.Both high pollution and retrieval difficulties tend to occur in many Mediterranean agricultural regions like the San Joauin. One difficulty is the relatively bright surfaces with considerable exposed soil. NASAs MAIAC and MODIS Deep Blue retrieval techniques are shown to have considerable skill even at low aerosol optical thickness (AOT) values, as evaluated by concurrent AERONET sunphotometer measurements.More significantly, these AOT values can correspond to high daytime PM2.5 since aerosol mixed layer depth is thin and variable, 200m 600 m. The thin layers derive from typical subsidence of dry air between more stormy periods. This situation provides an advantage: water vapor column is also almost completely limited to a similar mixed layer depth, and can thus serve as a measure of aerosol dilution.Using the MAIAC Water Vapor Column:In order to make the maps below, we used the MAIAC data but subtracted partial water-vapor columns estimated from MERRA Reanalysis Data availabe from the GSFC GMAO using kriging. We did not use the mixed-layer estimates from MERRA, since such analyses were found problematic during our forecasting exercises for DISCOVER-AQ. Observations from the aircraft soundings suggested that this overlying moisture was mostly due to larger scale flows, not ML venting.However, the specific humidity at the surface and a nearly well-mixed ML was analyzed by kriging from the surface network (MesoWest, University of Utah). These were thought to be truer, uninfluenced by physical process modeling that combines with data observations. (TBD: How different are they?) Procedure: Subtract overlying partial water columns from MAIAC column water and divide this by a surface value of water vapor. (MAIAC column is expressed in cm of water, i.e., water vapor at surface conditions.This method appears to bring out useful details in the distribution of submicron particles in the very problematic Wintertime San Joaquin Valley, and allow analysis of pollution episodes throughout the valley, rather than long-term averages.

san joaquin valley↗
DOE OSTI2022

The Mira–Titan Universe – IV. High-precision power spectrum emulation

Modern cosmological surveys are delivering data sets characterized by unprecedented quality and statistical completeness; this trend is expected to continue in the future as new ground- and space-based surveys come online. In order to maximally extract cosmological information from these observations, matching theoretical predictions are needed. At low redshifts, the surveys probe the non-linear regime of structure formation where cosmological simulations are the primary means of obtaining the required information. The computational cost of sufficiently resolved large-volume simulations makes it prohibitive to run very large ensembles. Nevertheless, precision emulators built on a tractable number of high-quality simulations can be used to build very fast prediction schemes to enable a variety of cosmological inference studies. We have recently introduced the Mira–Titan Universe simulation suite designed to construct emulators for a range of cosmological probes. This gravity-only set of simulations covers the standard six cosmological parameters {ω m , ω b , σ 8 , $h, n_s, w_0$} and, in addition, includes massive neutrinos and a dynamical dark energy equation of state {ω ν , $w_a$}. In this paper, we present the final emulator for the matter power spectrum based on 111 cosmological simulations, each covering a (2.1 Gpc) 3 volume and evolving 3200 3 particles. In this work, an additional set of 1776 lower resolution simulations and TimeRG perturbation theory results for the power spectrum are used to cover scales straddling the linear to mildly non-linear regimes (maximum wavenumber $\textit{k}$ = 5 Mpc –1 ). The emulator provides predictions at the 2–3 percent level of accuracy over a wide range of cosmological parameters and is publicly released as part of this paper.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

The 3D Lyman- α forest power spectrum from eBOSS DR16

We measure the three-dimensional power spectrum (P3D) of the transmitted flux in the Lyman-α (Ly α) forest using the complete extended Baryon Oscillation Spectroscopic Survey data release 16 (eBOSS DR16). This sample consists of ~205 000 quasar spectra in the redshift range 2 ≤ z ≤ 4 at an effective redshift z = 2.334. We propose a pair-count spectral estimator in configuration space, weighting each pair by exp( i k ∙ r), for wave vector k and pixel pair separation r, effectively measuring the anisotropic power spectrum without the need for fast Fourier transforms. This accounts for the window matrix in a tractable way, avoiding artefacts found in Fourier-transform based power spectrum estimators due to the sparse sampling transverse to the line of sight of Ly α skewers. We extensively test our pipeline on two sets of mocks: (i) idealized Gaussian random fields with a sparse sampling of Ly α skewers, and (ii) log-normal LyaCoLoRe mocks including realistic noise levels, the eBOSS survey geometry and contaminants. On eBOSS DR16 data, the Kaiser formula with a non-linear correction term obtained from hydrodynamic simulations yields a good fit to the power spectrum data in the range $(0.02 ≤ k ≤ 0.35)$ h Mpc -1 at the 1–2σ level with a covariance matrix derived from LyaCoLoRe mocks. We demonstrate a promising new approach for full-shape cosmological analyses of Ly α forest data from cosmological surveys such as eBOSS, the currently observing Dark Energy Spectroscopic Instrument and future surveys such as the Prime Focus Spectrograph, WEAVE-QSO, and 4MOST.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Development of Tunable, Precision Permanent Magnet Phase Shifter for X-Ray FEL

In an X-ray free-electron laser (XFEL), a high quality relativistic electron beam from a linear accelerator traverses an array of undulator magnets to generate a high power X-ray beam by allowing the FEL to reach saturation. In many XFELs, phase shifters haven’t been commonly used or included in their designs. The importance of the phase shifter was recognized only recently. A phase shifter can be constructed using either an electromagnet (EM) or a permanent magnet (PM). EM phase shifters are heavy, complex to operate, and take more valuable longitudinal space between the undulators. Compared to EMs, the PMs have several outstanding advantages, including compactness, lighter weights, simple operation, and low maintenance. This project has developed an innovative design of fixed-gap, tunable PM-based compact phase shifter for the XFEL. While keeping the poles and gap fixed, the new design incorporates proven features from the prior work, along with novel field tuning methods, which allow the field strength to be adjusted within a reasonably large range to achieve the desired range of the phase delay control. The PM design without electric currents and power electronics associated with EMs improves reliability by reducing the number of possible points of failure and lowers the operation cost. The compactness makes it possible to use PM phase shifters in the existing XFELs with small gaps between adjacent undulators, as well as in the new XFELs without significantly increasing the FEL length. In Phase I we have accomplished the followings; 1) determined the optimum physics parameters for both weak and strong phase shifters; 2) conducted magnetic design with concrete magnetic and geometrical parameters of actual produced PM materials; 3) provided solutions to achieve precision fine tuning of the magnetic fields with unique methods; 4) elaborated a design based upon ease of manufacturing and cost analysis; and 5) produced rare-earth PM materials, which will be used for prototyping in Phase II. The phase I results provide a foundation for the research and development in the Phase II project, which will aim at building two full-size PM phase shifters and testing at SLAC National Accelerator Laboratory. This project will significantly provide the potential to bring a paradigm change in accelerator technology, replacing conventional electromagnets for many existing accelerator applications and allowing the development of new applications in completely new areas. For example, the PM FEL phase shifters can help significantly improve coherence and power of an XFEL. The proposed use of PM phase shifters will also bring significant cost savings to the construction, operation and maintenance of FEL facilities.

43 PARTICLE ACCELERATORS↗
NASA NTRS1958

Three-Dimensional Orbits of Earth Satellites, Including Effects of Earth Oblateness and Atmospheric Rotation

The principal purpose of the present paper is to present sets of equations which may be used for calculating complete trajectories of earth satellites from outer space to the ground under the influence of air drag and gravity, including oblateness effects, and to apply these to several examples of entry trajectories starting from a circular orbit. Equations of motion, based on an "instantaneous ellipse" technique, with polar angle as independent variable, were found suitable for automatic computation of orbits in which the trajectory consists of a number of revolutions. This method is suitable as long as the trajectory does not become nearly vertical. In the terminal phase of the trajectories, which are nearly vertical, equations of motion in spherical polar coordinates with time as the independent variable were found to be more suitable. In the first illustrative example the effects of the oblateness component of the earth's gravitational field and of atmospheric rotation were studied for equatorial orbits. The satellites were launched into circular orbits at a height of 120 miles, an altitude sufficiently high that a number of revolutions could be studied. The importance of the oblateness component of the earth's gravitational field is shown by the fact that a satellite launched at circular orbital speed, neglecting oblateness, has a perigee some 67,000 feet lower when oblateness forces are included in the equations of motion than when they are not included. Also, the loss in altitude per revolution is double that of a satellite following an orbit not subject to oblateness. The effect of atmospheric rotation on the loss of altitude per revolution was small. As might be surmised, the regression of the line of nodes as predicted by celestial mechanics is unchanged when drag is included. It is clear that the inclination of the orbital plane to the equator will be relatively unaffected by drag for no atmospheric rotation since the drag lies in the orbital plane in this case. With the inclusion of atmospheric rotation it was found that the inclination of the plane changed about one-millionth of a radian per revolution. Thus the prediction of the position of the orbital plane of an earth satellite is not complicated by the introduction of drag. The line of apsides, which without drag but with oblateness moves slowly in space, tends to move with the satellite when drag is included in the calculations. As a results, the usual linearized solutions based on oblateness alone must be basically altered when drag is included to take into account the rapid movement of the line of apsides. In the second illustrative example the final revolution was calculated to impact for a number of trajectories in an orbital plane inclined at 650 to the equator. Of particular interest is the large effect the oblateness gravitational field and atmospheric rotation can have on the impact point. For a value of CDA/m of unity, and for an initial downward angle at 80 miles altitude of 0.01 radian, such as might be utilized for manned re-entry, oblateness had an influence of about 300 miles in the impact point, and atmospheric rotation had about a 150-mile influence.

Nielsen, Jack N.↗
DOE OSTI2022

MAGNETICALLY SUSPENDED CANNED ROTOR PUMPS FOR THE INTEGRAL MOLTEN SALT REACTOR

This project developed a molten salt pump design for small modular reactors, concentrated solar, and Gen IV nuclear reactors with a magnetically levitated rotor that can operate at temperatures of up to 700 °C. This eliminates the need for rotating seals and roller element bearings which require maintenance and are prone to failure. The pump design is also more compact than existing molten salt pump designs which require long shafts to thermally isolate the electric motor. The aims of the project are to: 1) design and test the molten salt pump, 2) develop and test a high temperature position sensor needed for magnetic bearing control, and 3) develop and test high temperature coil fabrication methods for electric motor and magnetic bearing fabrication. To develop the high temperature coils, conductor materials were analyzed for oxidation, diffusion, cladding, and resistivity. Of the candidate conductor materials that met the requirements, silver was chosen. The long-term impact of radiation on silver resistivity was analyzed by developing a model of a generic molten salt fast reactor to model the conductor irradiation in worst case scenarios. The analysis showed that increased resistivity due to transmutation of silver decreases exponentially with the reactor reflector thickness. Different electrical insulation materials were tested and the final high temperature coils developed were able to withstand conductor to ground voltages up to 1800 V before insulation breakdown. Magnetic bearings require high-precision high-speed measurements of the rotor position to operate. Current commercial position sensors can only operate up to 550 °C so a novel position sensor was developed that can operate up to 800 °C. These sensors were also designed to measure the position of the metallic rotor through a thin metallic corrosion barrier that protects that coils, insulation, and magnetic cores so that the salt containment barrier does not have any penetrations. The high-temperature sensor was tested and had comparable performance to commercially available inductive position sensors. The sensor design and signal processing has led to several invention disclosures. The magnetically levitated molten salt pump design was completed and analyzed. The chloride salt test loop at Oak Ridge National Laboratory was chosen to test the pump and the pump hydraulics and power were designed for the test loop flow and pressure requirements. Classes of motor designs that do not utilize permanent magnetic were considered for the high-temperature motor. After extensive analysis of their performance, a synchronous reluctance motor was chosen for the pump. The design was refined and analyzed to optimize performance. The magnetic bearings were designed to meet the force and frequency requirements necessary to levitate the rotor and reject disturbances from rotor imbalance, motor forces, and hydraulic forces. Compliant mounts for the position sensors were designed that would not damage the ceramic sensor components and maintain sensor alignment over a wide temperature range. A thermal hydraulic analysis of the molten salt coolant flow used to cool the rotor and stator was performed to ensure that the maximum temperature in the rotor and stator will not exceed 725 °C. The corrosion barrier and casing were designed to maintain the salt fluid boundary by using compliant features to absorb large expansion mismatches due to different coefficients of thermal expansion (CTE). Finally, the salt drain and fill systems were designed along with the pump support structure, integration with the test loop, and thermal insulation and heating. Corrosion barrier fabrication experiments were performed to assess using laser cladding to apply a sub-millimeter corrosion barrier to the surface of a Fe-Co-V magnetic core material. The cladding method showed minimal mixing between materials at their interface and withstood temperature cycling without damage to the cladding. The technology commercialization is focused on the high-temperature position sensors and there has been significant interest from the nuclear industry in utilizing these sensors along with some interest from the aerospace industry. The project also led to a research collaboration with a fusion energy industrial partner to study barriers to scaling the pump design from the kW size to the MW size. This research collaboration will also study methods to create high-temperature coils that have conductor to ground breakdown voltages above 5 kV. The successful development of high-temperature electromagnetic coils and high-temperature position sensors greatly increases the temperature limits for electromagnetic devices that don’t rely on permanent magnets. This includes electric motors, magnetic bearings, linear actuators, rotary position and velocity sensing, actuated valves, and generators to name a few. This will fill a current need in molten salt reactors, concentrated solar, and fusion energy for these critical peripheral devices needed for practical reactor designs that do not currently have commercial solutions.

Hines, J Wesley↗
DOE OSTI2024

Comparing pion production in transport simulations of heavy-ion collisions at 270⁢𝐴 MeV under controlled conditions

Within the TMEP, we present a detailed study of the performance of different transport models in Sn+Sn collisions at 270A MeV, and put particular emphasis on the production of pions and Delta resonances, which have been used as probes of the nuclear symmetry energy. We prescribe a common and rather simple physics model, and follow in detail the results of 4 BUU models and 6 QMD models. The nucleonic evolution of the collision and the nucleonic observables in these codes do not completely converge, but the differences among the codes can be understood as being due to several reasons: the basic differences between BUU and QMD models in the representation of the phase-space distributions, computational differences in the mean-field evaluation, and differences in the adopted strategies for the Pauli blocking in the collision integrals. For pionic observables, we find that a higher maximum density leads to an enhanced pion yield and a reduced pi- pi+ yield ratio, while a more effective Pauli blocking generally leads to a slightly suppressed pion yield and an enhanced pi- pi+ yield ratio. We specifically investigate the effect of the Coulomb force, and find that it increases the total pi- pi+ yield ratio but reduces the ratio at high pion energies, although differences in its implementations do not have a dominating role in the differences among the codes. Taking into account only the results of codes that strictly follow the homework specifications, we find a convergence of the codes in the final charged pion yield ratio to a 1 sigma deviation of about 5%. Furthermore, the uncertainty is expected to be reduced to about 1.6% if the same or similar strategies and ingredients, i.e., an improved Pauli blocking and calculation of the non-linear term in the mean-field potential, are similarly used in all codes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Beyond fragmented dopant dynamics in quantum spin lattices: Robust localization and non-Gaussian diffusion

The motion of dopants in magnetic spin lattices has received tremendous attention for at least four decades due to its connection to high-temperature superconductivity. Despite these efforts, we lack a complete understanding of their behavior, especially out of the equilibrium and at nonzero temperatures. In this paper, we take a significant step towards a much deeper understanding based on state-of-the-art matrix-product-state calculations. In particular, we investigate the nonequilibrium dynamics of a dopant in two-leg 𝑡−𝐽 ladders with antiferromagnetic XXZ spin interactions. In the Ising limit, we find that the dopant is localized for all investigated nonzero temperatures due to an emergent disordered potential, with a localization length controlled by the underlying correlation length of the spin lattice, which increases exponentially with decreasing temperature. The dopant, hereby, only delocalizes asymptotically in the zero temperature limit. This greatly generalizes the localization effect discovered recently in Hilbert space fragmented models [Phys. Rev. Res. 6, 023325 (2024), SciPost Phys. Core 7, 054 (2024)]. In the presence of spin-exchange processes at rate 𝛼, the dopant diffuses with a diffusion coefficient, 𝐷 ℎ , depending nonmonotonically on 𝛼. It initially increases linearly as 𝐷 ℎ ∝ 𝛼 for 𝛼 ≪ 1 before dropping off as 𝛼 −1 for 𝛼 > 1. Moreover, we show that the underlying spin dynamics at infinite temperature behaves qualitatively the same, albeit with important quantitative differences. We substantiate these findings by showing that the dynamics features self-similar scaling behavior, which strongly deviates from the Gaussian behavior of regular diffusion, especially for weak spin exchange. Finally, we show that the diffusion coefficient 𝐷 ℎ follows an Arrhenius relation at high temperatures, whereby it is exponentially suppressed upon cooling.

Anomalous diffusion↗
NASA NTRS2010

Performance Bounds on Two Concatenated, Interleaved Codes

A method has been developed of computing bounds on the performance of a code comprised of two linear binary codes generated by two encoders serially concatenated through an interleaver. Originally intended for use in evaluating the performances of some codes proposed for deep-space communication links, the method can also be used in evaluating the performances of short-block-length codes in other applications. The method applies, more specifically, to a communication system in which following processes take place: At the transmitter, the original binary information that one seeks to transmit is first processed by an encoder into an outer code (Co) characterized by, among other things, a pair of numbers (n,k), where n (n > k)is the total number of code bits associated with k information bits and n k bits are used for correcting or at least detecting errors. Next, the outer code is processed through either a block or a convolutional interleaver. In the block interleaver, the words of the outer code are processed in blocks of I words. In the convolutional interleaver, the interleaving operation is performed bit-wise in N rows with delays that are multiples of B bits. The output of the interleaver is processed through a second encoder to obtain an inner code (Ci) characterized by (ni,ki). The output of the inner code is transmitted over an additive-white-Gaussian- noise channel characterized by a symbol signal-to-noise ratio (SNR) Es/No and a bit SNR Eb/No. At the receiver, an inner decoder generates estimates of bits. Depending on whether a block or a convolutional interleaver is used at the transmitter, the sequence of estimated bits is processed through a block or a convolutional de-interleaver, respectively, to obtain estimates of code words. Then the estimates of the code words are processed through an outer decoder, which generates estimates of the original information along with flags indicating which estimates are presumed to be correct and which are found to be erroneous. From the perspective of the present method, the topic of major interest is the performance of the communication system as quantified in the word-error rate and the undetected-error rate as functions of the SNRs and the total latency of the interleaver and inner code. The method is embodied in equations that describe bounds on these functions. Throughout the derivation of the equations that embody the method, it is assumed that the decoder for the outer code corrects any error pattern of t or fewer errors, detects any error pattern of s or fewer errors, may detect some error patterns of more than s errors, and does not correct any patterns of more than t errors. Because a mathematically complete description of the equations that embody the method and of the derivation of the equations would greatly exceed the space available for this article, it must suffice to summarize by reporting that the derivation includes consideration of several complex issues, including relationships between latency and memory requirements for block and convolutional codes, burst error statistics, enumeration of error-event intersections, and effects of different interleaving depths. In a demonstration, the method was used to calculate bounds on the performances of several communication systems, each based on serial concatenation of a (63,56) expurgated Hamming code with a convolutional inner code through a convolutional interleaver. The bounds calculated by use of the method were compared with results of numerical simulations of performances of the systems to show the regions where the bounds are tight (see figure).

Moision, Bruce↗
DOE OSTI2023

X-ray absorption spectroscopy study of Mn reference compounds for Mn speciation in terrestrial surface environments

Abstract X-ray absorption spectroscopy (XAS) offers great potential to identify and quantify Mn species in surface environments by means of linear combination fit (LCF), fingerprint, and shell-fit analyses of bulk Mn XAS spectra. However, these approaches are complicated by the lack of a comprehensive and accessible spectrum library. Additionally, molecular-level information on Mn coordination in some potentially important Mn species occurring in soils and sediments is missing. Therefore, we investigated a suite of 32 natural and synthetic Mn reference compounds, including Mn oxide, oxyhydroxide, carbonate, phosphate, and silicate minerals, as well as organic and adsorbed Mn species, by Mn K-edge X-ray absorption near edge structure (XANES) and extended X-ray absorption fine structure (EXAFS) spectroscopy. The ability of XAS to infer the average oxidation state (AOS) of Mn was assessed by comparing XANES-derived AOS with the AOS obtained from redox titrations. All reference compounds were studied for their local (<5 Å) Mn coordination environment using EXAFS shell-fit analysis. Statistical analyses were employed to clarify how well and to what extent individual Mn species (groups) can be distinguished by XAS based on spectral uniqueness. Our results show that LCF analysis of normalized XANES spectra can reliably quantify the Mn AOS within ~0.1 v.u. in the range +2 to +4. These spectra are diagnostic for most Mn species investigated, but unsuitable to identify and quantify members of the manganate and Mn(III)-oxyhydroxide groups. First-derivative XANES fingerprinting allows the unique identification of pyrolusite, ramsdellite, and potentially lithiophorite within the manganate group. However, XANES spectra of individual Mn compounds can vary significantly depending on chemical composition and/or crystallinity, which limits the accuracy of XANES-based speciation analyses. In contrast, EXAFS spectra provide a much better discriminatory power to identify and quantify Mn species. Principal component and cluster analyses of k2-weighted EXAFS spectra of Mn reference compounds implied that EXAFS LCF analysis of environmental samples can identify and quantify at least the following primary Mn species groups: (1) Phyllo- and tectomanganates with large tunnel sizes (2 × 2 and larger; hollandite sensu stricto, romanèchite, todorokite); (2) tectomanganates with small tunnel sizes (2 × 2 and smaller; cryptomelane, pyrolusite, ramsdellite); (3) Mn(III)-dominated species (nesosilicates, oxyhydroxides, organic compounds, spinels); (4) Mn(II) species (carbonate, phosphate, and phyllosilicate minerals, adsorbed and organic species); and (5) manganosite. All Mn compounds, except for members of the manganate group (excluding pyrolusite) and adsorbed Mn(II) species, exhibit unique EXAFS spectra that would allow their identification and quantification in mixtures. Therefore, our results highlight the potential of Mn K-edge EXAFS spectroscopy to assess bulk Mn speciation in soils and sediments. A complete XAS-based speciation analysis of bulk Mn in environmental samples should preferably include the determination of Mn valences following the “Combo” method of Manceau et al. (2012), EXAFS LCF analyses based on principal component and target transformation results, as well as EXAFS shell-fit analyses for the validation of LCF results. For this purpose, all 32 XAS reference spectra are provided in the Online Materials1 for further use by the scientific community.

Geochemistry & Geophysics↗
NASA NTRSDate not supplied

An Optimization Approach to Support Science Decision Making for Lunar Surface Exploration

Introduction: Scientific exploration is one of the three pillars of NASA’s Moon2Mars architecture, with crew surface extra vehicular activities (EVA) serving a critical enabling function. Development of surface EVA operational planning and execution, specifically integrating science and flight control teams (FCT), is currently being explored through analog scenarios. This integration, exercised, for example, through the Joint EVA and Hu-man Surface Mobility Test Team (JETT), allows for science input on EVA activities in near real-time through a Science Evaluation Room (SER), or Arte-mis science backroom, which integrates with the broader FCT through the Science Officer. The SER works within the FCT to support dynamic EVA planning in response to changes in operational constraints as well as science opportunities and re-prioritization, increasing the mission science return and accelerating the accomplishment of the Moon2Mars science objectives. The SER works within the FCT to provide recommendations to traverse execution in near real-time. One challenge is the requirement to deliver SER inputs to the FCT on operationally relevant timelines. Failure to do so may result in suboptimal execution of science exploration EVAs or even loss of key science objectives. To close this gap, we present a network optimization tool to allow the SER to provide rapid input to the FCT in response to changes in operational constraints or science opportunities. Inputs are predicated on approved science objectives, and clear rationale must be provided to the FCT for any requested change. Accordingly, this tool incorporates the Science Traceability Matrix (STM), SER prioritization scheme, and station characterization and action planning with operational constraints such as duration, traverse speed, and distance to maximize science objectives based on SER priorities, consistent with FCT operational requirements. Method: As a proof of concept, we used an existing linear programing software package used to simulate optimal routes through cellular metabolism. We built a Demonstrative Model with three STM objectives and four stations on a region of the Moon. The objectives were given an arbitrary prioritization and mapped to the stations through four possible crew actions. (Figs. 1 and 2). This station to STM mapping is consistent with the method used by the JETT5 Science Team to develop analog surface EVA science planning. We used a grid system with the landing site at the origin and the four stations placed across the positive x,y quadrant. Actions were assigned to each station and the accomplishment of those actions resulted in a numerical “reward” based on the ability of that action to achieve science objectives. The aggregate reward from each individual STM objective contributes to a global score (Science Yield), weighted by its priority. Operational constraints included a requirement to start and end at the landing site, 5 minutes each for initial station characterization and “clean up,” and variable total EVA time, traverse rate (fixed to 0.5 meters per second in our example), and time to perform each action (10, 5, 7, and 15 min for actions 1, 2, 3, and 4, respectively). Additional constraints and variables will be added in the future (e.g., sample mass, number of stations, traverse route constraints, illumination). Optimization. We converted the connections (arcs) between these stations (nodes) into a mixed integer linear programming optimization problem (arcs = constraints, nodes = variables) with the objective to maximize Science Yield. For any action, the Science Yield is equal to the relevance of that action to an STM objective [3, 2, and 1 point(s) for High, Med., and Low relevance, respectively], multiplied by the STM Objective Priority [3, 2, and 1 point(s) for High, Med., and Low priority, respectively]. This resulted in a model that computes the optimal station and action combination to maximize the Science Yield. These weightings can be adjusted by the SER as desired. Results: We explored three test cases for the Demonstrative Model. First, we set the maximum EVA duration to 120 minutes and computed the optimal route (Fig. 3A). The model suggested per-forming Actions 1 and 2 at Station P01, followed by Actions 1 and 2 at Station P02, and finally Actions 1 and 3 at Station P04 before returning to the Landing Site. Second, we adjusted the STM Objective Priori-ty order and computed the new optimal route (Fig. 3B). Under this situation, the model suggested per-forming all Actions at Station P02 followed by all Actions at Station P03. The previous test cases were relevant to SER planning activities. Next, we explored providing mid-EVA replanning input to the FCT. Scenario: While executing the Route in Fig. 3A the crew finishes at Station P01 and FCT decides that the EVA needs to finish in 45 minutes back at the Landing Site. FCT asks SER to recommend changes to the plan to accommodate this operation-al change. Using the model and incorporating these new constraints (start at Station P01, max. time of 45 min), the model suggested performing Actions 2 and 4 at Station P03 (Fig. 4), requiring 41 minutes to complete and return to the Landing Site. Interestingly, Station 3 was not part of the original route. Using the model, we determined the EVA would need 66 minutes, instead of 45, in order for the original Station P04 to yield a larger Science Yield than Station P03. The parametrization and simulation was per-formed in less than a minute, demonstrating the operational relevance of the approach. Future Efforts: The results from the Demonstrative Model suggest this tool can accelerate SER decision making on operationally relevant timelines. Use in analog activities, such as JETT5 or follow-ons, which have over a dozen stations for a crew to explore and over a dozen actions per station, will provide needed validation of the utility of this tool for planning EVAs, replanning mid-EVA, or planning follow-on EVAs based on previous results. Further integration with FCT execution monitoring tools may provide additional efficiency gains, al-lowing rapid and iterative exploration of operation-al and science decision space by the FCT and SER.

Science Operations↗