Search NASA⌕ Search

SEARCH · Search NASA

Results for “Time dependent fluctuations”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
NASA NTRS2000

Section on Observed Impacts on El Nino

Agricultural applications of El Nino forecasts are already underway in some countries and need to be evaluated or re-evaluated. For example, in Peru, El Nino forecasts have been incorporated into national planning for the agricultural sector, and areas planted with rice and cotton (cotton being the more drought-tolerant crop) are adjusted accordingly. How well are this and other such programs working? Such evaluations will contribute to the governmental and intergovernmental institutions, including the Inter-American Institute for Global Change Research and the US National Ocean and Atmospheric Agency that are fostering programs to aid the effective use of forecasts. As El Nino climate forecasting grows out of the research mode into operational mode, the research focus shifts to include the design of appropriate modes of utilization. Awareness of and sensitivity to the costs of prediction errors also grow. For example, one major forecasting model failed to predict the very large El Nino event of 1997, when Pacific sea-surface temperatures were the highest on record. Although simple correlations between El Nino events and crop yields may be suggestive, more sophisticated work is needed to understand the subtleties of the interplay among the global climate system, regional climate patterns, and local agricultural systems. Honesty about the limitations of an forecast is essential, especially when human livelihoods are at stake. An end-to-end analysis links tools and expertise from the full sequence of ENSO cause-and-effect processes. Representatives from many disciplines are needed to achieve insights, e.g, oceanographers and atmospheric scientists who predict El Nino events, climatologists who drive global climate models with sea-surface temperature predictions, agronomists who translate regional climate connections in to crop yield forecasts, and economists who analyze market adjustments to the vagaries of climate and determine the value of climate forecasts. Methods include historical studies to understand past patterns and to test hindcasts of the prediction tools, crop modeling, spatial analysis and remote sensing. This research involves expanding, deepening, and applying the understanding of physical climate to the fields of agronomy and social science; and the reciprocal understanding of crop growth and farm economics to climatology. Delivery of a regional climate forecast with no information about how the climate forecast was derived limits its effectiveness. Explanation of a region's major climate driving forces helps to place a seasonal forecast in context. Then, a useful approach is to show historical responses to previous El Nino events, and projections, with uncertainty intervals, of crop response from dynamic process crop growth models. Regional ID forecasts should be updated with real-time weather conditions. Since every El Nino event is different, it is important to track, report and advise on each new event as it unfolds. The stability of human enterprises depends on understanding both the potentialities and the limits of predictability. Farmers rely on past experience to anticipate and respond to fluctuations in the biophysical systems on which their livelihoods depend. Now scientists are improving their ability to predict some major elements of climate variability. The improvements in the reliability of El Nino forecasts are encouraging, but seasonal forecasts for agriculture are not, and will probably never be completely infallible, due to the chaotic nature of the climate system. Uncertainties proliferate as we extend beyond Pacific sea-surface temperatures to climate teleconnections and agricultural outcomes. The goal of this research is to shed as a clear light as possible on these inherent uncertainties and thus to contribute to the development of appropriate responses to El Nino and other seasonal forecasts for a range of stakeholders, which, ultimately, includes food consumers everywhere.

Rosenzweig, Cynthia↗
DOE OSTI2023

Rapid AI-based Dissection of Ashes using Raman and XRF Spectroscopy (RADAR-X)

Waste-to-energy (WTE) facilities incinerate ~35 million tons of municipal solid waste annually in the United States. The incineration process reduces the mass and the volume of the waste fraction by over 75 and 95%, respectively. The fraction left after incineration remains as ash residues and is referred to as WTE ash, compromising of bottom and fly ash. In the United States, ~10 million tons of WTE ashes are generated annually and predominantly landfilled because of the lack of secondary end-use pathways. This incurs a significant financial burden (landfilling costs) on U.S. WTE facilities and also results in the loss of materials to the landfill. The primary objective of this research is to understand better the elemental and mineralogical composition of WTE ashes from diverse origins and find composition dependent upcycling pathways for diverting these ashes from landfills. This primary objective was addressed through three research tasks –(Task I) An AI-led Multi-Modal Approach for Compositional Analysis, (Task II) Developing a Dissolution-Based Test for Real Time Analysis, and (Task III) Establishing composition-dependent end uses. The chemical composition of WTE ash is dependent on two factors, i.e., the input waste composition and the operational parameters of a WTE facility (combustion conditions). Amongst these two factors, the input waste composition will likely show spatial and temporal variation. As a result, the chemical composition of WTE ash will also fluctuate. To understand the spatial and temporal variation in WTE ash composition, in Task I, we collected 128 ash samples (62 bottom ash and 66 fly ash samples) from 11 WTE facilities located in 11 U.S. states and characterized them via X-ray Fluorescence, powder X-ray Diffraction, and Raman Spectroscopy. The findings from this extensive characterization work indicated that the key elements in WTE fly ashes are Ca, Cl (greater than 10 wt. %), Si, S, K, Zn ( between 1 and 10 wt. %), Mg, Al, P, Ti, Fe, Cu, Br, and Pb (between 0.1 and 1 wt. %). Similarly, the key elements in WTE bottom ash fraction finer than 45μm are Ca (greater than 10 wt. %), Mg, Al, Si, S, Cl, K, Ti, Fe, Zn (between 1 and 10 wt. %), P, V, Cr, Mn, Cu, Br, and Pb (between 0.1 and 1 wt. %). Here, we note that the dominant fraction of WTE bottom ash is the coarse fraction. The coarse WTE bottom ash fraction (rich in silicon) was not characterized in this study because of excessive grinding requirements and their unsuitability as a supplementary cementitious material due to their coarse nature. The elements in WTE bottom ashes are present as calcite, anhydrite, vaterite, hydroxyapatite, quartz, bassanite, gehlenite, akermanite, hydrocalumite, and portlandite. Similarly, the mineralogical species present in WTE fly ashes are calcium chloride hydroxide, halite, calcite, anhydrite, sylvite, hydrocalumite, vaterite, hannebachite, quartz, and bassanite. Temporal variation in ash composition may also result in significant fluctuations in chemical compositions. Therefore, a WTE facility may need to monitor the ash composition (elemental and mineralogical composition) in real time. In Task I, we evaluated the possibility of using a portable X-ray fluorescence (XRF) spectrometer to monitor the elemental composition in real-time. Specifically, we collected XRF data on identical specimens via a portable XRF spectrometer (low-end) and a lab-based benchtop XRF spectrometer (high-end). The collected data was used to train an A.I. algorithm (portable XRF data as an input and benchtop XRF data as an output) to predict accurate elemental composition using portable XRF data. Finally, we developed a 2-minute photobleaching protocol to monitor the mineralogical characteristics of WTE ashes via Raman spectroscopy. Overall, the activities in Task I improved our understanding of ash composition and developed techniques to monitor elemental and mineralogical composition in near real-time. Based on the findings of Task I, we find that WTE ashes exhibit wide variability in mineralogy. For ICP-based elemental analysis, all the mineralogical species in WTE ashes must be brought into solution. This is traditionally accomplished with acid digestion using a combination of multiple acids. However, acid digestion with multiple acids is time-consuming and often fails to ensure complete digestion of the ash matrix. To address this limitation, in Task II, we developed an alkali-fusion-based digestion protocol using a combination of lithium tetraborate, lithium metaborate, and their combinations as possible alkali fluxes for digesting WTE ashes entirely and rapidly. The validity of the developed method was evaluated on two standard ash specimens, i.e., SRM 1633c coal fly ash and BCR-176R incineration fly ash specimen. The findings suggest that the developed protocol can ensure complete digestion of elements such as Al, Ba, Ca, Cr, Cu, Mg, Mn, P, Sr, V, Zn, Be, K, and rare earth elements. The recent changes in the energy market towards renewables and increased metal recycling have resulted in reduced supplies of supplementary cementitious materials (coal fly ash and slag). Therefore, in Task III, we evaluated the possibility of employing WTE ashes as SCMs. As the chemical composition of WTE ashes varies temporally (on an hourly basis), there was also a need to develop tests that can evaluate the suitability of material to act as supplementary cementitious material rapidly, i.e., in a few minutes. Therefore, in Task II, we also developed a rapid test to assess the suitability of a material to act as an SCM in 5 minutes. This represents a significant advance over the state-of-the-art R 3 test, which takes ~144 hours. This test was initially validated on amorphous aluminosilicates, such as calcined clays, and could be extended to evaluate other industrial by-products, such as WTE ashes. In Task III, we evaluated the possibility of employing WTE ashes for two applications, i.e., as an SCM and a lime substitute for clay stabilization. The findings from Task I indicated that WTE ashes are enriched in chlorine and, therefore, cannot be used directly as an SCM due to corrosion-related risks and altered hydration kinetics. Accordingly, we developed an ash treatment protocol to reduce the solubility of chlorine-containing species in WTE ashes. The developed treatment protocol also immobilized lead in certain mineral forms. As a result of the treatment, WTE ashes can be used as SCMs without any corrosion or heavy metal leaching concerns. The second application examined in this study was clay stabilization. WTE ashes are calcium-rich and can be an adequate lime replacement for clay stabilization. Our findings reveal that the sum of the concentrations of Ca(OH) 2 and CaClOH controls the clay stabilization capability of WTE ashes. In summary, in this work, we evaluated the elemental and mineralogical characteristics of U.S. WTE ashes from diverse origins and developed tests to evaluate the chemical characteristics of these ashes in real time through a portable XRF and a benchtop Raman spectrometer. Based on the chemical characteristics of these ashes, we developed an ash treatment process to enable the use of WTE ashes as an SCM and also evaluated the possibility of employing these ashes for clay stabilization. Overall, the findings from this work enables the diversion of WTE ashes from landfills for multiple end-uses, i.e., as an SCM or a lime substitute for clay stabilization.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Faster network disruption from layered oscillatory dynamics

Nonlinear complex network-coupled systems typically have multiple stable equilibrium states. Following perturbations or due to ambient noise, the system is pushed away from its initial equilibrium, and, depending on the direction and the amplitude of the excursion, it might undergo a transition to another equilibrium. It was recently demonstrated [M. Tyloo, J. Phys. Complex. 3 03LT01 (2022)] that layered complex networks may exhibit amplified fluctuations. Here, I investigate how noise with system-specific correlations impacts the first escape time of nonlinearly coupled oscillators. Interestingly, I show that, not only the strong amplification of the fluctuations is a threat to the good functioning of the network but also the spatial and temporal correlations of the noise along the lowest-lying eigenmodes of the Laplacian matrix. Finally, I analyze first escape times on synthetic networks and compare noise originating from layered dynamics to uncorrelated noise.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Probing iron in Earth’s core with molecular-spin dynamics

Dynamic compression of iron to Earth-core conditions is one of the few ways to gather important elastic and transport properties needed to uncover key mechanisms surrounding the geodynamo effect. Herein, a machine-learned ab initio derived molecular-spin dynamics (MSD) methodology with explicit treatment for longitudinal spin-fluctuations is utilized to probe the dynamic phase-diagram of iron. This framework uniquely enables an accurate resolution of the phase-transition kinetics and Earth-core elastic properties, as highlighted by compressional wave velocity and adiabatic bulk moduli measurements. In addition, a unique coupling of MSD with time-dependent density functional theory enables gauging electronic transport properties, critically important for resolving geodynamo dynamics.

58 GEOSCIENCES↗
DOE OSTI2022

Ab initio calculations for void swelling bias in α - and δ -plutonium

Void swelling can develop in materials under persistent irradiation when nonequilibrium vacancy and self-interstitial populations migrate under sufficiently asymmetric interaction biases. In conventional metals, the propensity is determined to the first approximation by comparing point-defect relaxation strains. We thus present DFT-based calculations of structures and formation energies and volumes of point defects in the α and the δ phases of plutonium. Here we discuss the pros and cons of various levels of electronic structure theory: spin polarization, spin-orbit coupling, and orbital polarization. Our results show that lattice defects in δ-Pu, in contrast to most fcc metals, have surprisingly small formation volumes. Equally unexpected are the large defect formation volumes found in the low-symmetry α-Pu phase. Both these unusual properties can be satisfactorily explained from defect-induced spin/orbital moment formation and destruction in the Pu phases. Surprisingly, the point defects in α-Pu are found to induce far larger transformation of the local electronic structure than in δ-Pu. When we use the calculated defect properties to estimate the classic void swelling bias in each of the phases, we find it to be unusually small in δ-Pu but likely much larger in α-Pu. Hence, swelling rates and mechanisms can diverge dramatically between the different phases of Pu. Especially in the transient regime before the formation of large defect clusters, the swelling rate of α-Pu can reliably be expected to be much larger than δ-Pu. However, accurate forecasts over longer times will require the conventional void-swelling theory to be modified to handle the complexities presented by the different Pu phases. As a case in point, we show the possible anomalous temperature dependence of vacancy properties in δ-Pu, caused by entropic contributions from defect-induced spin-lattice fluctuations. Such complications may affect defect-defect interactions and thus alter the void swelling bias.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

FEED Study of Carbon Capture Inc. DAC and CarbonCure Utilization Using United States Steel's Gary Works Waste Plant

Direct Air Capture (DAC) has been proposed as a means of reducing atmospheric concentrations of CO2. While DAC has been evaluated through lab-scale, bench-scale, and small-scale pilot units, large-scale deployment has not been achieved. Feasibility studies are one tool to understand the potential design, operation, performance, and impact of commercial-scale DAC. This paper presents the results of a feasibility study for a passive DAC system deployed at >100,000 tpy scale in three different regions across the U.S. and awarded to Carbon Collect by the U.S. Department of Energy National Energy Technology Laboratory (DOE NETL). The MechanicalTree™technology has been designed and engineered by Carbon Collect based on the concept initially developed at Arizona State University Center for Negative Carbon Emissions. It uses a tower of stacked, sorbentcontaining disks supported by a lifting mechanism, exposed to the air to capture CO2 during the adsorption phase. The disks are then lowered into a regeneration chamber for vacuum and steam regeneration to produce CO2 product during the desorption step. The modular tree design allows large installations with repeatable, mass-manufactured units connected to common utilities such as steam supply, vacuum, and CO2 processing for compression and geologic storage. The passive DAC system eliminates the equipment and energy of forced air fans by using natural air circulation to contact the sorbent with the CO2 in the air. Because of this, the performance is dependent on the wind speed in addition to the temperature and relative humidity. Performance and flow rate fluctuations were incorporated into equipment and facility design with considerations for turndown to 10% of maximum rated flow to allow operation in all seasons. The feasibility study was undertaken to evaluate the technology in different regions and climates. Three locations were selected for this study representing different climates: Alabama (hot and humid), California (hot and arid), and Wyoming (continental). For each region, the adsorption/desorption cycle was optimized including tuning the heat integration and cycle timings. The performance of individual trees was then scaled to the full facility with the same design of more than 20,000 trees common between all regions. The installations had expected average CO2 capture rates of between 330,000 and 485,000 tonnes of CO2 per year depending on the climate. The initial engineering design for the facility at each location was performed with cost and performance estimates for the trees, carbon purification and compression, and the balance of plant. To supply thermal and electrical energy for the facility, carbon-free or low-carbon power must be considered. Options for low-carbon, continuous thermal and electrical energy were considered. The best-performing option from those considered was identified to be an electrically-heated molten salt energy storage system, powered by an on-site photovoltaic field. During molten salt thermal discharge, steam is produced from heat exchange with the molten salt and used to generate power in a steam turbine as well as steam for regenerating the DAC carbon trees. The thermal and electrical supply and analysis is presented in the context of low-carbon power for carbon removal. Modelling and economics of transportation and geologic storage of the CO2 is considered and presented for each location. The results of the feasibility study are incorporated into the presented techno-economic and life-cycle assessments. This work is intended to provide an understanding of the performance, cost, and impact of capturing CO2 at the commercial scale and the impact of climate and regional siting on the considered passive DAC system.

42 ENGINEERING↗
DOE OSTI2024

Statistical analysis of magnetic divertor configuration influence on H-mode transitions

DIII-D plasmas are compared for two upper divertor configurations: with the outer strike point on the small angle slot (SAS) divertor target and with the outer strike point on the horizontal divertor target (HT). Scanning the vertical distance between the magnetic null point and the divertor target over a range 0.10–0.16 m is shown to increase the threshold power, $P_\mathrm{th}$, and edge plasma power, $P_\mathrm{Loss},$ for the low-to-high confinement (L–H) and H–L transitions respectively, by up to a factor of 1.4. The X-point height scans were performed at three L-mode core plasma line average electron densities, $\bar{n}_\mathrm{e} = $ 1.2, 2.2 and 3.6 $\times 10^{19}\,\mathrm{m}^{-3}$, to investigate the density dependence of divertor magnetic configuration influence on $P_\mathrm{th}$. The X-point height, $Z_\textrm{x-pt}$, was further extended across the range 0.16–0.22 m with the more open HT divertor configuration, for which a clear decrease in $P_\mathrm{th}$ with increasing $Z_\textrm{x-pt}$ is observed. The dependence of $P_\textrm{th}$ on divertor magnetic geometry is further investigated using a time-dependent probability density function (PDF) model and information geometry to elucidate the roles played by pedestal plasma turbulence and perpendicular velocity flows. The degree of stochasticity of the plasma turbulence is observed to be sensitive to the plasma heating rate. The calculated square of the information rate shows changes in the relative density fluctuations and perpendicular velocity PDFs begin 2–5 ms prior to the L–H transition for three plasmas; providing a crucial measurement of the dynamic timescale of external transport barrier formation. Additionally, both information length and rate provide potential predictors of the L–H transition for these plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1978

Simulation of turbulent wall pressure

A Monte Carlo procedure was developed to simulate turbulent boundary layer wall pressure fluctuations. The approach utilizes much of the newly available conditional sampling information to construct the required distribution functions. Various disturbance wave forms were examined, as well as the effect of frequency-dependent decay. Good agreement between the simulation and experimental data was achieved for root mean square pressure level, power spectrum, and space time correlation.

Ash, R. L.↗
NASA NTRS2021

Physics and Chemistry of Te and HgTe-based Ternary Melts

Historically, the theoretical treatment of the liquid phases has always been more difficult and complicated than that for solid and gas phases. A liquid has no lattice structure as crystalline solids and the atoms/molecules in the liquid can migrate through it relatively rapidly. On the other hand, it is also interacting with many other atoms/molecules so that the simplifications of the kinetic theory of gases cannot be employed. For more complicated liquids, such as the liquids of high ionicity and those containing hydrogen bonds and electric dipoles, the understanding is far from complete. At the same time, accurate information on the physics and chemistry of semiconductor melts is needed for the quantitative descriptions of the process of crystal growth from melt. The pre-crystallization phenomena in the liquid phase are critical because the properties of the grown crystals depend on the state and structure of the melt as well as the thermal history of the melt during solidification process. However, the data on the liquid phase, such as thermophysical properties of semiconductor melts are scarce, especially for the HgTe-based II-VI ternary compound semiconductors because of their high vapor pressure and extreme toxicity. Analysis of the thermophysical properties of the melt can provide information about structural transitions of the melt during the solidification process. From a broader point of view, the structure of liquids is much more complicated than the crystalline solids, especially the relaxation behavior through different thermal histories. The theory of hetero-phase fluctuations of liquids is applicable to any many-body systems including condensed-matter physics, field theory, physics of nuclear-matter, cosmology, biology and even sociology. This book summarizes the physics and chemistry from the experimental measurements and theoretical analyses of phase diagram, thermodynamic properties, density, thermal conductivity, viscosity, and electrical conductivity on the binary, pseudo-binary and ternary melts of the most advanced IR-detector material systems of HgCdTe and HgZnTe as well as the analyses of these results. The main objectives of this study are: (1) to provide the phase diagrams and thermodynamic properties of Hg-Cd-Te and Hg-Zn-Te systems through quantitatively fitting the experimental data by assuming an associated solution model for the liquid phase, (2) to experimentally measure the thermophysical properties of the Hg-Cd-Te and Hg-Zn-Te melts, including density, viscosity, electrical conductivity and thermal conductivity as functions of temperature and composition and (3) to enhance the fundamental knowledge of hetero-phase fluctuations and relaxation phenomena in the melts and extend our understanding of the solidification process in order to interpret the experimental results of crystal growth so as to improve the melt growth processes of the compound semiconductor. The physics and chemistry of Te and HgTe-based ternary melts were explored through the studies of the structural transformation during melting, the supercooling during solidification, the relaxation phenomena after rapid cooling of the melts and the metal-semiconductor transition in the melts through the analyses of electrical conductivity and Lorenz number. An in-depth study on the thermophysical properties and their time-dependent structural dynamic processes taking place in the vicinity of the solid-liquid phase transition of the narrow homogeneity range HgTe-based ternary semiconductors as well as the analysis of the homogenization process in the melt will also be presented.

Thermophysical properties↗
NASA NTRS2019

Comparing Flood-Producing Catastrophic Hurricanes in GEOS: A Study of Juan 1985 and Harvey 2017

Hurricanes Juan (1985) and Harvey (2017) bear remarkable similarities because of the exceptional amount of accumulated precipitation over coastal regions in the Gulf, consequent to the peculiar and rare situation in which a hurricane stalls for several days on a coastal area, displaying minimal progress. This is particularly dangerous because: a) terrestrial drag affect only part of the circulation so that weakening is not as fast as when there is rapid motion inland; b) a part of the circulation remains over water and as such low-level moisture flux convergence is maintained; c) precipitation lasts for a very long time over the same area; d) moisture recycling from flooded areas further contributes to maintaining flood conditions. From the forecast perspective, hurricanes with track singularities such as stalling points, loops, or multiple landfalls over the same region are exceptionally difficult to predict because they result from small variations of opposing mechanisms and are therefore intrinsically unstable. Intensity fluctuations and track singularities are strongly coupled. The representation of fluxes undergoing differential drag (ocean vs. land) inside the same circulation is critically dependent upon resolution. Recent modeling development efforts within the NASA Global Modeling and Assimilation Office (GMAO) have brought the horizontal resolution of the global cubed-sphere Goddard Earth Observing System (GEOS) atmospheric model down to 12 km. This has led to innovative findings in the field of tropical cyclone (TC) modeling within a global framework. The GEOS atmospheric model and data assimilation system (DAS), used to produce the MERRA-2 (Modern Era Retrospective analysis for Research and Applications) reanalysis, has a "Replay" capability that adds a forcing term to the model equations to constrain it to follow a specified trajectory. This capability is used to produce global, downscaled versions of the 50 km MERRA-2 reanalysis at a 12.5 km resolution. The full DAS and replay modes both use the Incremental Analysis Update or IAU for applying an analysis increment to the model's background state four times each day at 00Z, 06Z, 12Z and 18Z. In a full DAS mode, a six-hour free run of the model is made (Predictor Step), writing states that will be used as backgrounds or first-guesses for the analysis centered on 00Z. The Analysis system then produces its best estimate of the 00Z state by blending backgrounds with observations. The difference between this best estimate and the model's original first guess at 00Z is treated as an analysis increment and added uniformly as a forcing term to the model during a second, now forced, run of the 12Z to 03Z window, which produces the assimilated results. Replay mode is very similar, except only the 00Z background is needed from the predictor segment, which is thus half as long, and instead of doing an analysis, a previously analyzed state at 00Z is read in. The increment is then computed and applied to the corrector step in the same way as in the original DAS. When downscaling horizontal spatial filters can be applied to the "analysis" increment used to force the Corrector step, so as to reproduce the full meteorological state or selectively choose or discard certain forcings. In this presentation, full global hindcast simulations of Juan and operational forecasts of Harvey are compared.

Flooding↗
DOE OSTI2026

Data-driven particle dynamics: Structure-preserving coarse-graining for emergent behavior in non-equilibrium systems

Multiscale systems are ubiquitous in science and technology, but are notoriously challenging to simulate as short spatiotemporal scales must be appropriately linked to emergent bulk physics. When expensive high-dimensional dynamical systems are coarse-grained into low-dimensional models, the entropic loss of information leads to emergent physics which are dissipative, history-dependent, and stochastic. To machine learn coarse-grained dynamics from time-series observations of particle trajectories, we propose a framework using the metriplectic bracket formalism that preserves these properties by construction; most notably, the framework guarantees discrete notions of the first and second laws of thermodynamics, conservation of momentum, and a discrete fluctuation-dissipation balance crucial for capturing non-equilibrium statistics. We introduce the mathematical framework abstractly before specializing to a particle discretization. As labels are generally unavailable for entropic state variables, we introduce a novel self-supervised learning strategy to identify emergent structural variables. We validate the method on benchmark systems and demonstrate its utility on two challenging examples: (1) coarse-graining star polymers at challenging levels of coarse-graining while preserving non-equilibrium statistics, and (2) learning models from high-speed video of colloidal suspensions that capture coupling between local rearrangement events and emergent stochastic dynamics. We provide open-source implementations in both PyTorch and LAMMPS, enabling large-scale inference and extensibility to diverse particle-based systems.

Computational Engineering, Finance, and Science (c↗
NASA NTRS1980

Electromagnetic scattering lifetimes for dust in Jupiter's ring

The lifetime of the charged, micron-sized dust particles found to lie in a ring at 1.8 Jupiter radii from Jupiter against electromagnetic scattering is estimated. Interactions of the particles with the Jovian magnetic field dipole are found to spread particle orbits out to a much greater extent than observed for particle charges greater than -10 V, and to result in slightly eccentric, open orbits for the dust particles, with eccentricities dependent on particle size. In addition, fluctuations within the magnetic field are found to lead to significant changes in the eccentricity and inclination of the ring particle orbits. It is concluded that Jovian ring particles, if charged to -10 V as is consistent with Voyager data, will be strongly perturbed and scattered on a short time scale (100 yr), so the ring must be constantly replenished or it is a temporary phenomenon. Finally, the dust cloud observed to extend above, below, and beyond the Jupiter ring is attributed to the electromagnetic scattering of particles smaller than the ring particles.

Consolmagno, G. J.↗
NASA NTRS1988

Inflow Measurements Made with a Laser Velocimeter on a Helicopter Model in Forward Flight. Volume 1: Rectangular Planform Blades at an Advance Ration of 0.15

An experimental investigation was conducted in the 14- by 22-Foot Subsonic Tunnel at NASA Langley to measure the inflow into a scale model helicopter rotor in forward flight (microinf = 0.15). The measurements were made with a two component Laser Velocimeter (LV) one chord above the plane formed by the path of the rotor tips (tip path plane). A conditional sampling technique was employed to determine the azimuthal position of the rotor at the time each velocity measurement was made so that the azimuthal fluctuations in velocity could be determined. Measurements were made at a total of 147 separate locations in order to clearly define the inflow character. This data is presented without analysis. In order to increase the availability of the resulting data, both the mean and azimuthally dependent values are included as part of this report on two 5.25 inch floppy disks in Microsoft Corporation MS-DOS format.

Elliott, Joe W.↗
NASA NTRS1989

Age structure of refractory interstellar dust and isotopic consequences

A sputtering and recycling Monte Carlo model, developed by Liffman and Clayton (1988) is used to calculate the distribution of existence times of the matter in interstellar dust composed of refractory metals. The mean age of each dust particle is defined not as the time it has existed but rather as the mass-weighted existence times of its parts at t = 6 Gyr of the modeled solar system formation. It is shown that Galactic evolution generates a mean correlation, applying to large numbers of particles binned according to size rather than according to individual particles, whose mean ages fluctuate statistically. The cosmochemical consequence is that if interstellar particles can be dynamically sorted into separate size populations during the aggregation history of solar system bodies, the collections of larger grains will constitute matter that is chemically older than collections of smaller grains. The macroscopic age difference generates isotopic anomalies by virtue of the time dependence of the secondary/primary nucleosynthesis yields. Results are compared with three different prescriptions for the sputtering of interstellar dust.

Clayton, Donald D.↗
NASA NTRS1992

Properties of quasi-periodic oscillations in accreting magnetic white dwarfs

Previous studies of time-dependent accretion onto magnetic white dwarfs, in which the cooling was assumed to be due to bremsstrahlung emission, have shown that the accretion shock undergoes oscillations. However, when cyclotron cooling is also included, the oscillations are damped for sufficiently strong magnetic fields. Here we demonstrate that the oscillations can be sustained by accretion-fluctuation-induced excitations. The frequency of the QPOs are shown to increase quadratically with the magnetic field strength. We interpret the oscillations as a two-phase process in which bremsstrahlung cooling dominates in one half-cycle and cyclotron cooling in the other. Such a process may have very different consequences compared to a single-phase process where the functional form of the cooling is essentially the same throughout the cycle. If in the two-phase process damping occurs mainly in the cyclotron cooling half-cycle, there will be a universal effective damping factor which tends to suppress all oscillation modes indiscriminately. The oscillations of the accretion shock also could be a limit cycle process in which the system vacillates between two branches.

Wu, Kinwah↗
DOE OSTI2022

Frequency-domain vs time-domain SP{sub N} equations to simulate neutron noise

Inside the reactor core mechanical vibrations of fuel assemblies can produce high fluctuations around a steady-state configuration, known as neutron noise. This effect can cause the triggering of power reduction measures. Classically, diffusion theory has been used to simulate this behavior. However, this equation has some limitations if the materials of the reactor have strong variations. In this work, we use the diffusive time-dependent simplified spherical harmonics equations that improve the previous results without the necessity of using high computational requirements. In particular, two types of analyses with these equations (SP{sub 3}) are made: a frequency-domain and a time-domain. A numerical neutron noise benchmark tests the methodology and compare both formulations. First, numerical results show a good agreement between the amplitudes and phases of the SP3 equations computed with the frequency-domain and time-domain. Therefore, as the frequency-domain computation only requires to solve a linear system, it is a recommendable option for neutron noise computations. Second, one can conclude that for this type of nuclear systems, where the assemblies are not homogenized, the SP{sub 3} approximation results improve considerably the accuracy of the diffusion theory.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1994

Excitation of solar p-modes

We investigate the rates at which energy is supplied to individual p-modes as a function of their frequencies nu and angular degrees l. The observationally determined rates are compared with those calculated on the hypothesis that the modes are stochastically excited by turbulent convection. The observationally determined excitation rate is assumed to be equal to the product of the mode's energy E and its (radian) line width Gamma. We obtain E from the mode's mean square surface velocity with the aid of its velocity eigenfuction. We assume that Gamma measures the mode's energy decay rate, even though quasi-elastic scattering may dominate true absorption. At fixed l, E(Gamma) arises as nu(exp 7) at low nu, reaches a peak at nu approximately equal 3.5 mHz, and then declines as nu(exp 4.4) at higher nu . At fixed nu, E(Gamma) exhibits a slow decline with increasing l. To calculate energy input rates, P(sub alpha), we rely on the mixing-length model of turbulent convection. We find entropy fluctuations to be about an order of magnitude more effective than the Reynolds stress in exciting p-modes . The calculated P(sub alpha) mimic the nu(exp 7) dependence of E(Gamma) at low nu and the nu(exp -4.4) dependence at high nu. The break of 11.4 powers in the nu-dependence of E(Gamma) across its peak is attributed to a combination of (1) the reflection of high-frequency acoustic waves just below the photosphere where the scale height drops precipitously and (2) the absence of energy-bearing eddies with short enough correlation times to excite high-frequency modes. Two parameters associated with the eddy correlation time are required to match the location and shape of the break. The appropriate values of these parameters, while not unnatural, are poorly constrained by theory. The calculated P(sub alpha) can also be made to fit the magnitude of E(Gamma) with a reasonable value for the eddy aspect ratio. Our resutls suggest a possible explanation for the decline of mode energy with increasing l at fixed nu. Entropy fluctuations couple to changes in volume associated with the oscillation mode. These decrease with decreasing n at fixed nu, becoming almost zero for the f-mode.

Goldreich, Peter↗
DOE OSTI2020

Quantum Versus Classical Switching Dynamics of Driven Dissipative Kerr Resonators

Here, we report a first-principles study of the driven dissipative dynamics for Kerr oscillators in the mesoscopic regime. This regime is characterized by large Kerr nonlinearity, realized here using the nonlinear kinetic inductance of a large array of Josephson junctions. The experimentally measured nonlinear resonance lineshapes of the junction array modes show significant deviations from steady-state numerical predictions, and necessitate time-dependent numerical simulations indicative of strong measurement-induced dephasing in the system arising from the large cross-Kerr effect between array modes. Analytical and numerical calculations of switching rate corroborate this by showing the emergence of a slow time scale, which is much longer than the linear decay rate and is set by fluctuation-induced switching times in the bistable regime. Furthermore, our analysis shows that the usual quantum-activated escape treatment is inadequate for prediction of the switching rates at large frequency shifts caused by strong nonlinearities, necessitating a quantum treatment that utilizes the full system Liouvillian. Based on our analysis, we identify a universal crossover parameter that delineates the regimes of validity of semiclassical and quantum descriptions, respectively. Our work shows how dynamical switching effects in strongly nonlinear systems provide a platform to study quantum-to-classical transitions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗