SEARCH · Engineering Papers

Results for “thermal & statistical models”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 415 records · Page 23

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Disordered hyperuniform solid state materials

Disordered hyperuniform (DHU) states are recently discovered exotic states of condensed matter. DHU systems are similar to liquids or glasses in that they are statistically isotropic and lack conventional long-range translational and orientational order. On the other hand, they completely suppress normalized infinite-wavelength density fluctuations like crystals and, in this sense, possess a hidden long-range correlation. Very recently, there have been several exciting discoveries of disordered hyperuniformity in solid-state materials, including amorphous carbon nanotubes, amorphous 2D silica, amorphous graphene, defected transition metal dichalcogenides, defected pentagonal 2D materials, and medium/high-entropy alloys. It has been found that the DHU states of these materials often possess a significantly lower energy than other disorder models and can lead to unique electronic and thermal transport properties, which results from mechanisms distinct from those identified for their crystalline counterparts. For example, DHU states can enhance electronic transport in 2D amorphous silica; DHU medium/high-entropy alloys realize the Vegard's law and possess enhanced electronic bandgaps and thermal transport at low temperatures. These unique properties open up many promising potential device applications in optoelectronics and thermoelectrics. Here, we provide a focused review on these important new developments of hyperuniformity in solid-state materials, taking an applied and “materials” perspective, which complements the existing reviews on hyperuniformity in physical systems and photonic materials. As a result, future directions and outlook are also provided, with a focus on the design and discovery of DHU quantum materials for quantum information science and engineering.

2D materials↗

NASA NTRS2012

Assessment of Mars Atmospheric Temperature Retrievals from the Thermal Emission Spectrometer Radiances

Motivated by the needs of Mars data assimilation. particularly quantification of measurement errors and generation of averaging kernels. we have evaluated atmospheric temperature retrievals from Mars Global Surveyor (MGS) Thermal Emission Spectrometer (TES) radiances. Multiple sets of retrievals have been considered in this study; (1) retrievals available from the Planetary Data System (PDS), (2) retrievals based on variants of the retrieval algorithm used to generate the PDS retrievals, and (3) retrievals produced using the Mars 1-Dimensional Retrieval (M1R) algorithm based on the Optimal Spectral Sampling (OSS ) forward model. The retrieved temperature profiles are compared to the MGS Radio Science (RS) temperature profiles. For the samples tested, the M1R temperature profiles can be made to agree within 2 K with the RS temperature profiles, but only after tuning the prior and error statistics. Use of a global prior that does not take into account the seasonal dependence leads errors of up 6 K. In polar samples. errors relative to the RS temperature profiles are even larger. In these samples, the PDS temperature profiles also exhibit a poor fit with RS temperatures. This fit is worse than reported in previous studies, indicating that the lack of fit is due to a bias correction to TES radiances implemented after 2004. To explain the differences between the PDS and Ml R temperatures, the algorithms are compared directly, with the OSS forward model inserted into the PDS algorithm. Factors such as the filtering parameter, the use of linear versus nonlinear constrained inversion, and the choice of the forward model, are found to contribute heavily to the differences in the temperature profiles retrieved in the polar regions, resulting in uncertainties of up to 6 K. Even outside the poles, changes in the a priori statistics result in different profile shapes which all fit the radiances within the specified error. The importance of the a priori statistics prevents reliable global retrievals based a single a priori and strongly implies that a robust science analysis must instead rely on retrievals employing localized a priori information, for example from an ensemble based data assimilation system such as the Local Ensemble Transform Kalman Filter (LETKF).

Hoffman, Matthew J.↗

DOE OSTI2023

Chemical controls on the propagation and healing of subcritical fractures

Human activities involving subsurface reservoirs—resource extraction, carbon and nuclear waste storage—alter thermal, mechanical, and chemical steady-state conditions in these systems. Because these systems exist at lithostatic pressures, even minor chemical changes can cause chemically assisted deformation. Therefore, understanding how chemical effects control geomechanical properties is critical to optimizing engineering activities. The grand challenge in predicting the effect of chemical processes on mechanical properties lays in the fact that these phenomena take place at molecular scales, while they manifest all the way to reservoir scales. To address this fundamental challenge, we investigated chemical effects on deformation in model and real systems spanning molecular- to centimeter scales. We used theory, experiment, molecular dynamics simulation, and statistical analysis to (1) identify the effect of simple reactions, such as hydrolysis, on molecular structures in interfacial regions of stressed geomaterials; (2) quantify chemical effects on the bulk mechanical properties, fracture and displacement for granular rocks and single crystals; (3) develop initial understanding of universal scaling for individual displacement events in layered geomaterials; and (4) develop analytic approximations for the single-chain mechanical response utilizing asymptotically correct statistical thermodynamic theory. Taken together, these findings advance the challenging field of chemo-mechanics.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Multiscale simulation of integrated energy system and electricity market interactions

Accelerating the deep decarbonization of the world's electric grids requires the coordination of complex energy systems and infrastructures across timescales from seconds to decades. Here, we present a new multiscale simulation framework that integrates process- and grid-centric modeling paradigms to better design, operate, and control integrated energy systems (IESs), which combine multiple technologies, in wholesale energy markets. Traditionally, IESs are analyzed with a process-centric paradigm such as levelized cost of electricity (LCOE) or annualized net revenue, ignoring important interactions with electricity markets. This framework explicitly models the complex interactions between an IES's bidding, scheduling, and control decisions and the energy market's clearing and settlement processes, while incorporating operational uncertainties. Through two case studies, we show the importance of understanding and quantifying complex resource-grid interactions. In case study 1, we demonstrate that optimized bidding from one resource shifts the profit distribution for all energy systems in the market. This result suggests new and more flexible IES technologies can disrupt the economics of all market participants, possibly leading to accelerated retirements of less flexible resources. Interestingly, the optimized bidding has little impact on grid-level aggregate statistics, such as total generation costs and renewable penetration rate. While aggregate modeling strategies may remain valid under some IES adoption scenarios for analysis focused on regional outcomes, direct comparisons of IES technologies at specific locations without considering these interactions may lead to misleading or incorrect conclusions. In case study 2, we consider the design and flexible operation of IESs that hybridize conventional generators with energy storage. Through a sensitivity analysis, we find that as the size of the storage system increases, the total number of start-ups for coal- and natural gas-based IESs reduced by 25% and 33.6%, and the total thermal generator ramping (i.e., mileage) reduced by 86.5% and 62.5%, respectively. This shows the primary benefit of storage may not be reduced operational costs (which do not change significantly) but fewer start-ups and less ramping, which may greatly simplify the design, operation, and control of carbon capture systems. The new modeling and optimization capabilities from this work enable the coupling of rigorous, dynamic process models with grid-level production cost models to quantitatively identify the nuanced interdependencies across these vast timescales that must be addressed to realize clean, safe, and secure energy production. Moreover, the proposed general multiscale simulation framework is applicable to all IES technologies and can be easily extended to consider other energy carriers (e.g., hydrogen, ammonia) and energy infrastructures (e.g., natural gas pipelines).

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2022

Worldwide Physics-Based Lifetime Prediction of c-Si Modules Due to Solder-Bond Failure

Lifetime prediction of the fielded c-Si solar modules due to location-specific weather conditions has been an important topic of photovoltaic research and the economic viability of solar energy. Data analytic techniques such as the performance ratio method, Statistical clear sky model, and Suns-Vmp methods quantify the degradation from measured data of a solar farm, however, the nonlinear time-dependence and correlated degradations make it difficult to use the empirical degradation rates for ultimate lifetime projection. In this article, we propose a complementary physics-based model to predict the solder bond failure caused by mechanical stress associated with the variations of the temperature. Integrating the worldwide weather information from NASA/NSRDB databases, the model predicts the location-specific output-power degradation and the lifetime of a module due to solder bond failure. The model parameters are calibrated against qualification tests involving thermal cycling of specific batches of modules from a specific technology/manufacturer. The results may be summarized as: 1) Modules installed at higher latitudes show a longer lifetime due to reduced damage accumulation. 2) The reduction of temperature fluctuation close to large bodies of water, such as seashores, increases solder bond lifetime significantly. 3) Relatively speaking, modules installed close to the Tropic of Cancer/Capricorn (23.5 degrees North/South) suffer from a higher solder bond damage and have a shorter lifetime, suggesting a conservative design. This model should serve as a building block of a comprehensive reliability framework that can predict the lifetime of a module that experiences simultaneous and correlated degradation mechanisms involving yellowing, corrosion, and potential-induced degradation.

14 SOLAR ENERGY↗

NASA NTRS1993

Plasmasphere dynamics in the duskside bulge region: A new look at old topic

Data acquired during several multiday periods in 1982 at ground stations Siple, Halley, and Kerguelen and on satellites Dynamics Explorer 1, International Sun Earth Explorer 1, and GEOS 2 have been used to investigate thermal plasma structure and dynamics in the duskside plasmasphere bulge region of the Earth. The distribution of thermal plasma in the dusk bulge sector is difficult to describe realistically, in part because of the time integral manner in which the thermal plasma distribution depends upon on the effects of bulk cross-B flow and interchange plasma flows along B. While relatively simple MHD models can be useful for qualitatively predicting certain effects of enhanced convection on a quiet plasmasphere, such as an initial sunward entrainment of the outer regions, they are of limited value in predicting the duskside thermal plasma structures that are observed. Furthermore, use of such models can be misleading if one fails to realize that they do not address the question of the formation of the steep plasmapause profile or provide for a possible role of instabilities or other irreversible processes in plasmapause formation. Our specific findings, which are based both upon the present case studies and upon earlier work, include the following: (1) during active periods the plasmasphere appears to become divided into two entities, a main plasmasphere and a duskside bulge region. (2) in the aftermath of an increase in convection activity, the main plasmasphere tends (from a statistical point of view) to become roughly circular in equatorial cross section, with only a slight bulge at dusk; (3) the abrupt westward edge of the duskside bulge observed from whistlers represents a state in the evolution of sunward extending streamers; (4) in the aftermath of a weak magnetic storm, 10 to 30% of the plasma 'removed' from the outer plasmasphere appears to remain in the afternoon-dusk sector beyond the main plasmasphere. (5) outlying dense plasma structures may circulate in the outer duskside magetosphere for many days following an increase in convection, unless there is extremely deep quieting; (6) a day-night plasmatrough boundary may be identified in equatorial satellite data; (7) factor-of-2-to-10 density irregularities appear near the plasmatrough from the ionosphere at L = 4.6, predominantly bidirectional field aligned and equatorially trapped light ion pitch angle distributions give away to a predominantly isotropic distribution (as seen by DE 1) when the plasma density reaches a level a factor of about 3 below the satured plasmasphere level; (9) some outlying dense plasma structures are effectively detached from the main plasmasphere, while others appear to be connected to that body.

Carpenter, D. L.↗

NASA NTRS2000

Integration of the Empirical Exospheric GCPM Plasma Model into IRI

The Global Core Plasma Model (GCPM) is an empirical description of typical thermal magnetospheric plasma densities in the plasmasphere, plasma trough, and polar cap. The GCPM makes use of the International Reference Ionosphere (IRI) for low altitudes. Densities are continuous and smooth. Plasmaspheric ion composition is also included in the GCPM. For the purpose of supporting the ionospheric community, the densities derived in the GCPM for high altitudes will be expressed as an extension of the IRI. The GCPM exospheric extension of IRI improves the topside densities and provides typical thermal plasma densities for the plasmasphere, trough, and polar cap. The GCPM is modular, having been designed for continued improvement as statistical density and composition measurements become available.

Gallagher, Dennis L.↗

NASA NTRS2015

Phase-Angle Dependence of Determinations of Diameter, Albedo, and Taxonomy: A Case Study of NEO 3691 Bede

Parameters important for NEO risk assessment and mitigation include Near-Earth Object diameter and taxonomic classification, which translates to surface composition. Diameters of NEOs are derived from the thermal fluxes measured by WISE, NEOWISE, Spitzer Warm Mission and ground-based telescopes including the IRTF and UKIRT. Diameter and its coupled parameters Albedo and IR beaming parameter (a proxy for thermal inertia and/or surface roughness) are dependent upon the phase angle, which is the Sun-target-observer angle. Orbit geometries of NEOs, however, typically provide for observations at phase angles greater than 20 degrees. At higher phase angles, the observed thermal emission is sampling both the day and night sides of the NEO. We compare thermal models for NEOs that exclude (NEATM) and include (NESTM) night-side emission. We present a case study of NEO 3691 Bede, which is a higher albedo object, X (Ec) or Cgh taxonomy, to highlight the range of H magnitudes for this object (depending on the albedo and phase function slope parameter G), and to examine at different phase angles the taxonomy and thermal model fits for this NEO. Observations of 3691 Bede include our observations with IRTF+SpeX and with the 10 micrometer UKIRT+Michelle instrument, as well as WISE and Spitzer Warm mission data. By examining 3691 Bede as a case study, we highlight the interplay between the derivation of basic physical parameters and observing geometry, and we discuss the uncertainties in H magnitude, taxonomy assignment amongst the X-class (P, M, E), and diameter determinations. Systematic dependencies in the derivation of basic characterization parameters of H-magnitude, diameter, albedo and taxonomy with observing geometry are important to understand. These basic characterization parameters affect the statistical assessments of the NEO population, which in turn, affects the assignment of statistically-assessed basic parameters to discovered but yet-to-be-fully-characterized NEOs.

Albedo and taxonomy↗

DOE OSTI2021

Specifics of Thermal Transport in Graphene Composites: Effect of Lateral Dimensions of Graphene Fillers

In this work, we report on the investigation of thermal transport in noncured silicone composites with graphene fillers of different lateral dimensions. Graphene fillers are comprised of few-layer graphene flakes with lateral sizes in the range from 400 to 1200 nm and the number of atomic planes from 1 to ~100. The distribution of the lateral dimensions and thicknesses of graphene fillers has been determined via atomic force microscopy statistics. It was found that in the examined range of the lateral dimensions, the thermal conductivity of the composites increases with increasing size of the graphene fillers. The observed difference in thermal properties can be related to the average gray phonon mean free path in graphene, which has been estimated to be around ~800 nm at room temperature. The thermal contact resistance of composites with graphene fillers of 1200 nm lateral dimensions was also smaller than that of composites with graphene fillers of 400 nm lateral dimensions. The effects of the filler loading fraction and the filler size on the thermal conductivity of the composites were rationalized within the Kanari model. The obtained results are important for the optimization of graphene fillers for applications in thermal interface materials for heat removal from high-power-density electronics.

36 MATERIALS SCIENCE↗

DOE OSTI2026

SALSAA: a statistical approach to line shapes from an average atom

Ion-Stark line broadening is a key density diagnostic for hot dense plasmas relevant to inertial fusion and astrophysics. It is caused by interactions of a radiating ion with nearby perturbing ions, whose electric microfields lead to changes in bound-bound transition energies. Ion-Stark broadening becomes increasingly difficult to compute for complex, many-electron ions with myriad transitions. In this paper, we propose a simplified approach to ion-Stark broadening based on self-consistent ion distributions and electronic structure from an average-atom model. We find that this approach reproduces the line shape predictions of one traditional method for high-n K-shell emission lines in aluminum ions with accuracy sufficient for density diagnostics in thermal plasmas with equal ion and electron temperatures. We expect that this approach can be extended to provide a reasonable picture of ion-Stark broadening in many-electron ions, enabling rapid calculations of line profiles in complex spectra.

average-atom↗

NASA NTRS2000

Comments on "The Sensitivity Study of Radiative-Convective Equilibrium in the Tropics with a Convective Resolving Model"

In general, there are two broad scientific objectives when using cloud resolving models (CRMs or cloud ensemble models-CEMs) to study tropical convection. The first one is to use them as a physics resolving models to understand the dynamic and microphysical processes associated with the tropical water and energy cycles and their role in the climate system. The second approach is to use the CRMs to improve the representation of moist processes and their interaction with radiation in large-scale models. In order to improve the credibility of the CRMs and achieve the above goals, CRMs using identical initial conditions and large-scale influences need to produce very similar results. Two CRMs produced different statistical equilibrium (SE) states even though both used the same initial thermodynamic and wind conditions. Sensitivity tests to identify the major physical processes that determine the SE states for the different CRM simulations were performed. Their results indicated that atmospheric horizontal wind is treated quite differently in these two CRMs. The model that had stronger surface winds and consequently larger latent and sensible heat fluxes from the ocean produced a warmer and more humid modeled thermodynamic SE state. In addition, the domain mean thermodynamic state is more unstable for those experiments that produced a warmer and more humid SE state. Their simulated wet (warm and humid) SE states are thermally more stable in the lower troposphere (from the surface to 4-5 km in altitude). The large-scale horizontal advective effects on temperature and water vapor mixing ratio are needed when using CRMs to perform long-term integrations to study convective feedback under specified large-scale environments. In addition, it is suggested that the dry and cold SE state simulated was caused by enhanced precipitation but not enough surface evaporation. We find some problems with the interpretation of these three phenomena.

Tao, W.-K.↗

DOE OSTI2025

One Million Open-source Cislunar Orbits

Cislunar space, encompassing the region from geosynchronous orbit to beyond the Moon, is poised to become a cornerstone for future exploration, scientific discovery, and national security. Missions in this region, spanning durations from weeks to decades, require robust infrastructure and reliable transit capabilities. The complex gravitational influences of the Moon, Sun, and planets, along with thermal radiation from Earth and the Sun, lead to significant trajectory deviations, resulting in kilometer-scale errors within days. Leveraging the high-performance computing resources at Lawrence Livermore National Laboratory (LLNL), we have simulated one million high-fidelity cislunar trajectories, now publicly available via LLNL’s Green Data Oasis and the Unified Data Library. Generated using the open-source Space Situational Awareness Python package, these trajectories match the precision of commercial tools such as AGI’s Systems Tool Kit and NASA’s General Mission Analysis Tool. This data set is a valuable resource for reference, statistical analysis of cislunar orbit populations, and training machine learning models for rapid orbit classification with minimal observational input. Preliminary analysis reveals stable bands in Keplerian element space, particularly around five geosynchronous radii across a range of inclinations and eccentricities. Beyond this threshold, the Moon’s influence disrupts most unassisted orbits, though co-orbiting L4/L5 Lunar Trojans persist throughout the six-year simulation.

Astronomy and AstroPhysics↗

DOE OSTI2020

Temperature statistics in a radiatively heated particle-laden turbulent square duct flow

Radiation absorption by preferentially concentrated particles in a turbulent square duct flow is studied experimentally. The particle-laden flow is exposed to near-infrared radiation, and the gas phase temperature statistics are measured along the wall bisector of the duct. It is found that the instantaneous temperature fluctuations are comparable to the overall mean temperature rise. The temperature statistics at the duct centerline and near the wall are qualitatively different. The former reflects preferential concentration in isotropic flows while the latter displays evidence of particle clustering into streamwise elongated streaks. Comparison of the experimental data to a simplified heat transfer model suggests that the Lagrangian evolution of particle clusters and voids, and turbulent mixing in the vicinity of particle clusters, are important. This work was motivated by particle solar receiver technology, but the findings are also relevant to systems where there is localized heat release or mass transfer from disperse particles or droplets. It shows that obtaining Lagrangian histories of particle trajectories is an important next step towards understanding thermal transport phenomena in particle-laden turbulent flows.

42 ENGINEERING↗

NASA NTRS1981

The distribution of free electrons in the inner galaxy from pulsar dispersion measures

The dispersion measures of a sample of 149 pulsars in the inner Galaxy (absolute value of l 50 deg) were statistically analyzed to deduce the large-scale distribution of free thermal electrons in this region. The dispersion measure distribution of these pulsars shows significant evidence for a decrease in the electron scale height from a local value greater than the pulsar scale height to a value less than the pulsar scale height at galactocentric radii inside of approximately 7 kpc. An increase in the electron density (to a value around .15/cu cm at 4 to 5 kpc) must accompany such a decrease in scale height. There is also evidence for a large-scale warp in the electron distribution below the b + 0 deg plane inside the Solar circle. A model is proposed for the electron distribution which incorporates these features and Monte Carlo generated dispersion measure distributions are presented for parameters which best reproduce the observed pulsar distributions.

Harding, D. S.↗

NASA NTRS2018

Learning Model Structural Uncertainty with Gaussian Processes

The advent of commercially available quantum computers has marked the beginning of quantum computing as a reality. Both quantum gate and annealing computers have been released by major computer hardware companies. In this work, the D-Wave 2XTM quantum annealing computer housed at the NASA Advanced Systems computational facility is investigated to accelerate Machine Learning (ML) for image registration. NASA collects large amounts of images over the globe remotely using space-based monitoring. Images of a fixed areas of the land surface are taken over time. Due to the orbit of the sensors, the viewing angles deviate slightly, and it is necessary to align or register the images precisely to create image time series over the land surface. Unaligned images can lead to substantial analysis errors. These time-series are then used in modeling Earth Systems models such as hydrological, weather, and carbon monitoring models. In this work, we consider the Moderate Resolution Image Spectrometer (MODIS) data collected by the NASA's terra satellite. Artificial Neural Networks (ANNs) is a natural fit for ML modelling of images. Several successes have been reported using machine learning related to image processing. We investigate the use of ML to register MODIS images. ANNs are investigated in combination with a Restricted Boltzmann Machines (RBM) as an auto-encoder. We will present results showing the accuracy and efficiency of this approach.The D-Wave 2XTM quantum annealer samples the ground-state wave-function of a spin-Ising systems with quadratic interactions between qubits and a Chimera connectivity. The system sits in a ~15 mK thermal bath. One can think of the system as being placed in the ground state initially and subject to thermal excitations governed by Boltzmann statistics. If this is assumed true, one can use the statistics from the D-Wave 2XTM to train RBMs. Generating statistics for training Boltzmann machines is an NP-hard problem and constitutes the largest compute cost. We investigate the use of the D-Wave 2XTM to accelerate the training of the RBMs in our ANNs and report on the results.

Kouatchou, Jules↗

DOE OSTI2023

Coupled Decay Heat and Thermal Hydraulic Capability for Loss-of-Coolant Accident Simulations

As the nuclear energy industry considers ways to achieve improved economics in the current fleet of light-water reactors (LWRs), one possible approach is to operate each cycle for longer durations. This causes a greater portion of the fuel to be burned and reduces the frequency of outages, which ultimately reduces the cost to operate the reactor. However, this also leads to higher burnup fuels than has traditionally been allowed in these reactors. Thus, there are concerns about integrity of high-burnup (HBu) fuel, especially during accident conditions such as loss-of-coolant accidents (LOCAs), as shown by Capps et al.. To investigate these concerns, advanced modeling and simulation capabilities are being leveraged to determine the susceptibility of HBu fuel to fuel fragmentation, relocation, and dispersion (FFRD). Improvements have previously been made to fuel performance capabilities to more accurately model these phenomena; multiphysics simulations have also been conducted to determine the power and burnup histories of the HBu fuel, which are needed as inputs for the fuel performance calculations. Most recently, new statistical approaches have been developed to identify a subset of fuel rods that have greater FFRD susceptibility, reducing the total number of fuel performance simulations required. Prior LOCA simulations have relied on the TRACE systems code, which can model the core and primary loop during accident conditions. TRACE includes many models for various aspects of the primary loop, but two sets of models are important for this report. First, TRACE uses a lumped-fuel approach for modeling the core. This approximates the ~50,000 fuel rods in the core with a much smaller number of rods. The rods can be lumped in various ways as determined by the user. For example, one lumped rod may be used to represent all rods in an assembly, sometimes with an additional rod representing the hottest fuel rod. However, due to runtime constraints and complexity of modeling, a more common approach is to group several assemblies or larger regions of the core into single lumped rods. These lumping schemes apply not only to fuel rods but to flow channels as well. Second, TRACE has several different models for treating decay heat, ranging from pregenerated decay heat curves based on an ANSI/ANS-5.1 standard (hereinafter abbreviated simply as ANSI) to explicit time-dependent heat inputs from the user. None of these models account for differences in isotopics between different rods, which is an approximation the work in this report seeks to eliminate. This report focuses on the implementation of coupled decay heat capabilities in the Virtual Environment for Reactor Applications (VERA) code suite to address a gap identified in previous LOCA simulations. This constitutes an improvement for both the lumped-fuel and decay heat models in TRACE. VERA has been developed to perform high-fidelity, whole-core multiphysics simulations for LWRs. Previously, during the Consortium for Advanced Simulation of LWRs (CASL) program, the emphasis was on providing accurate steady-state analysis—with a secondary focus on reactivity insertion accident (RIA) analysis—to address operational challenges in the nuclear energy industry. Under the Department of Energy (DOE) Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, these capabilities are being extended to a broader range of transient analyses with the goal of quantifying the risk of fuel failures such as FFRD. To properly model such conditions with VERA, decay heat calculations have been integrated with the multiphysics to enable rod-by-rod thermal hydraulic (TH) conditions to be driven by the decay heat in long-running accidents such as LOCAs.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2025

Molecular fluctuations inhibit intermittency in compressible turbulence

In the standard picture of fully developed turbulence, highly intermittent hydrodynamic fields are nonlinearly coupled across scales, where local energy cascades from large scales into dissipative vortices and large density gradients. Microscopically, however, constituent fluid molecules are in constant thermal (Brownian) motion, but the role of molecular fluctuations in large-scale turbulence is largely unknown, and with rare exceptions, it has historically been considered irrelevant at scales larger than the molecular mean free path. Recent theoretical and computational investigations have shown that molecular fluctuations can impact energy cascade at Kolmogorov length scales. Here, we show that molecular fluctuations not only modify energy spectrum at wavelengths larger than the Kolmogorov length in compressible turbulence, but also significantly inhibit spatio-temporal intermittency across the entire dissipation range. Using large-scale direct numerical simulations of computational fluctuating hydrodynamics, we demonstrate that the extreme intermittency characteristic of turbulence models is replaced by nearly Gaussian statistics in the dissipation range. These results demonstrate that the compressible Navier–Stokes equations should be augmented with molecular fluctuations to accurately predict turbulence statistics across the dissipation range. Our findings have significant consequences for turbulence modelling in applications such as astrophysics, reactive flows and hypersonic aerodynamics, where dissipation-range turbulence is approximated by closure models.

compressible turbulence↗

DOE OSTI2026

The Thermal and Kinematic Sunyaev–Zeldovich Effect in Galaxy Clusters and Filaments Using Multifrequency Temperature Maps of the Cosmic Microwave Background: A399–A401 Cluster Pair Case Study

We present a multifrequency and multi-instrument methodology to study the physical properties of galaxy clusters and cosmic filaments using cosmic microwave background observations. Our approach enables simultaneous measurement of both the thermal (tSZ) and kinematic Sunyaev–Zeldovich (kSZ) effects, incorporates relativistic corrections, and models astrophysical foregrounds such as thermal dust emission. We do this by jointly fitting a single physical model across multiple maps from multiple instruments at different frequencies, rather than fitting a model to a single Compton-y map. We demonstrate the success of this method by fitting the A399–A401 galaxy cluster pair and filament system using archival data from the Planck satellite and new, targeted deep data from the Atacama Cosmology Telescope, covering 11 different frequencies over 14 maps from 30 GHz to 545 GHz. Our tSZ results are consistent with previous work using Compton-y maps. We measure the line-of-sight peculiar velocities of the cluster–filament system using the kSZ effect and find statistical uncertainties on individual cluster peculiar velocities of ≲600 km s −1 , which are competitive with current state-of-the-art measurements. Additionally, we measure the optical depth of the filament component with a signal-to-noise of 8.5σ and reveal hints of its morphology. This modular approach is well-suited for application to future instruments across a wide range of millimeter and submillimeter wavebands.

Gill, Ajay S. [National Research Council of Canada↗