Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “difference scaling”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Development of a Global Reference Surface Reflectance and BRDF Datasets from Geostationary Satellite Observations and AERONET Measurements

Surface reflectances and their dependency on illumination-view geometries (i.e., BRDF) are the foundation of many high-level satellite products for land and water monitoring. Yet it is difficult to evaluate the quality of satellite-based surface reflectances with ground-based measurements due to the spatial scale differences. In order to fill the gap, here we develop a reference dataset of surface reflectance and BRDF at the global AERONET sites with data streams from operational geostationary sensors including Himawari 8/9 AHI, GK-2A AMI, and GOES 16/17 ABI. Taking the top-of-atmosphere (TOA) reflectance and the site measured atmospheric aerosol optical depth (AOD) as the main inputs, we apply the GeoNEX-AC algorithm to performance accurate atmospheric correction and derive 10-minute surface reflectance and daily Ross-Thick-Li-Sparse (RTLS) BRDF parameters at AERONET sites where coincident AOD measurements and TOA observations are available from 2016 (for Himawari) or 2018 (for GOES) onwards. The algorithm ensures that the retrieved surface BRDF parameters, along with the site-measured AOD, allow the atmospheric radiative transfer model, SHARM, accurately simulate the observed TOA reflectance at diurnal and longer time scales. They are our best estimates of the surface optical properties and thus can serve as the “reference” to evaluate the performance of operational atmospheric correction algorithms (where AOD is assumed unknown and needs to be retrieved). The reference BRDF also allow us to evaluate the spectral band ratios between the SWIR (e.g., 2200 nm) and the visible (e.g., 650 nm) regions, which are commonly used in operational atmospheric correction algorithms. Finally, we demonstrate that the reference dataset can be used to develop potential data synergies between different GEO satellites as well as GEO-LEO sensors.

Weile Wang↗
DOE OSTI2021

Computationally Efficient Multiscale Neural Networks Applied to Fluid Flow in Complex 3D Porous Media

Abstract The permeability of complex porous materials is of interest to many engineering disciplines. This quantity can be obtained via direct flow simulation, which provides the most accurate results, but is very computationally expensive. In particular, the simulation convergence time scales poorly as the simulation domains become less porous or more heterogeneous. Semi-analytical models that rely on averaged structural properties (i.e., porosity and tortuosity) have been proposed, but these features only partly summarize the domain, resulting in limited applicability. On the other hand, data-driven machine learning approaches have shown great promise for building more general models by virtue of accounting for the spatial arrangement of the domains’ solid boundaries. However, prior approaches building on the convolutional neural network (ConvNet) literature concerning 2D image recognition problems do not scale well to the large 3D domains required to obtain a representative elementary volume (REV). As such, most prior work focused on homogeneous samples, where a small REV entails that the global nature of fluid flow could be mostly neglected, and accordingly, the memory bottleneck of addressing 3D domains with ConvNets was side-stepped. Therefore, important geometries such as fractures and vuggy domains could not be modeled properly. In this work, we address this limitation with a general multiscale deep learning model that is able to learn from porous media simulation data. By using a coupled set of neural networks that view the domain on different scales, we enable the evaluation of large ( $$>512^3$$ > 512 3 ) images in approximately one second on a single graphics processing unit. This model architecture opens up the possibility of modeling domain sizes that would not be feasible using traditional direct simulation tools on a desktop computer. We validate our method with a laminar fluid flow case using vuggy samples and fractures. As a result of viewing the entire domain at once, our model is able to perform accurate prediction on domains exhibiting a large degree of heterogeneity. We expect the methodology to be applicable to many other transport problems where complex geometries play a central role.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Implementation of a High-Fidelity Interface Resolving Method in Nek5000

The development and utilization of computational fluid dynamics (CFD) models for large, high-temperature electric melters in the Waste Treatment and Immobilization Plant (WTP) in eastern Washington State have proven to be valuable for various purposes. These models allow for a better understanding of the physio-chemical processes occurring within the melter vessels and can contribute to improving operational efficiency, throughput, and addressing operational issues related to vitrification. The CFD models employed for these melter vessels incorporate multiphase fluid flow and heat transfer simulations in different regions, including the plenum, cold cap, and molten glass regions. As the tank waste and glass formers are introduced into the melter, a reacting batch layer known as the cold cap forms on top of the molten glass. To enhance the melt rate, forced convection bubblers located at the bottom of the melters generate convection currents that help homogenize the molten glass and provide heat to the cold cap. As the bubbles rise through the highly viscous glass, they adopt a spherical-cap shape [1]. Meanwhile, the conversion of the batch to glass generates significant amounts of gases (such as water vapor, carbon dioxide, sulfur dioxide, and NOx) due to thermal decomposition [2]. These gases become trapped between the cold cap and molten glass, forming a foam layer [3]. For modeling multiphase flow in CFD and heat transfer simulations of waste glass melters at different scales, efforts are underway to augment the capabilities of the Nek5000 [4] and NekRS [5] open-source codes [6]. Nek5000/NekRS is a scalable and efficient spectral element code that has been successfully applied to a wide range of fluid dynamics problems. By leveraging the Nek5000/NekRS software, it becomes possible to model the melter systems more affordably and with lower computational requirements compared to currently utilized commercial CFD software. The specific objectives of this ongoing effort include: 1. Implementation of a level set method in Nek5000/NekRS: The level set method is a numerical technique commonly employed in CFD simulations to track and represent the interface between different phases or materials accurately. By incorporating this method into Nek5000/NekRS, the ability to simulate multiphase flows in waste glass melters at a high level of fidelity can be achieved. 2. Demonstration of capability for air bubbling through molten glass: As part of the development process, a specific case of air bubbling through molten glass will be simulated using the augmented Nek5000/NekRS code. This demonstration aims to showcase the ability of the software to accurately capture and analyze the complex phenomena involved in the multiphase flow within waste glass melters. By achieving these objectives, the improved Nek5000/NekRS code will offer a powerful computational tool for simulating and analyzing waste glass melter systems, enabling better understanding, optimization, and troubleshooting of these vitrification processes. The ability to accurately model and simulate multiphase flows has broad relevance across many industries and scientific domains, and the improved functionality can contribute to advancements in various fields beyond waste glass melter simulations.

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

Trigger of Standard and Blowout Solar Jets: I. Parametric Studies of the Magnetic Field Geometry

Context. Jets are dynamic, impulsive, well collimated plasma events developing at many different scales and in differentlayers of the solar atmosphere.Aims. Jets are believed to be induced by magnetic reconnection a process central to many astrophysical phenomena.Studying their dynamics can help us to better understand the processes acting in larger eruptive events (e.g. flares andcoronal mass ejection) as well as mass, magnetic helicity and energy transfer at all scale in the solar atmosphere. Therelative simplicity of their magnetic geometry and topology, compared with larger solar active events, makes jets idealcandidate for studying the fundamental role of reconnection in energetic events. Methods. In this study, using our state-of-the-art numerical solver ARMS, we present several parametric studies of anumerical three-dimensional magneto-hydrodynamic model of solar jet-like events. We study the impact of the magneticfield inclination and photospheric field distribution on the generation and properties of two types of observed solar jets:the so-called standard and blowout jets.Results. The present parametric studies validate our model of jets for different geometric properties of the magneticconfiguration. On one hand, a blowout jets is eventually always trigger for the range of parameters that we tested. The3D magnetic null-point configuration is a very robust structure for the energy storage and impulsive release characteristicof blowout jets. On the other hand, the existence of the standard jet depends on the magnetic geometry. We show thatthe reconnection occurring during the standard jet phase influences the properties of the trigger of the blowout jet.Conclusions. Our results allow us to better understand the energization, triggering and driving processes of standardand blowout jets. Our model allows us to predict the impulsiveness and energetics of jets in terms of the surroundingmagnetic field configuration. Finally we discuss the interpretation of standard jets and blowout jets and the physicalfactors that determine which type of jets will occur.

Pariat, E.↗
NASA NTRS1994

Lithospheric-scale buckling and thrust structures on Mars: The Coprates rise and south Tharsis ridge belt

Detailed photogeologic mapping documents for the first time the deformation sequence of the Coprates rise. Materials of Middle Noachian through perhaps Early Hesperian age were involved in the horizontal shortening of crustal and lithospheric rocks that produced uplift and asymmetric, east-vergent folding of these layers during the Late Noachian and perhaps Early Hesperian. Previous suggestions of normal faulting as the deforming mechanism are not supported by the new mapping or by dislocation models of rise topography. The most likely mechanism for producing the Coprates rise is folding due to buckling instability and perhaps thrust faulting. An extensive ridge belt concentric to south Tharsis is defined by a newly recognized system of approximately coeval long-wavelength topogtaphic ridges and swales that include the Coprates rise. Horizontal shortening of crust and lithosphere is manifested by periodically spaced structures at two different scales: wrinkle ridges spaced tens of kilometers apart and the south Tharsis system of ridges spaced hundreds of kilometers apart. Both types of structures may be associated with early volcanotectonic activity and crustal thickening in south Tharsis. These and other large ridges indicate that lithospheric buckling and thrust faulting may have been common on early Mars.

Schultz, Richard A.↗
DOE OSTI2020

Superdiffusive transport of energy in one-dimensional metals

Metals in one spatial dimension are described at the lowest energy scales by the Luttinger liquid theory. It is well understood that this free theory, and even interacting integrable models, can support ballistic transport of conserved quantities including energy. In contrast, realistic one-dimensional metals, even without disorder, contain integrability-breaking interactions that are expected to lead to thermalization and conventional diffusive linear response. We argue that the expansion of energy when such a nonintegrable Luttinger liquid is locally heated above its ground state shows superdiffusive behavior (i.e., spreading of energy that is intermediate between diffusion and ballistic propagation), by combining an analytical anomalous diffusion model with numerical matrix-product–state calculations on a specific perturbed spinless fermion chain. Different metals will have different scaling exponents and shapes in their energy spreading, but the superdiffusive behavior is stable and should be visible in time-resolved experiments.

Bulchandani, Vir B.↗
DOE OSTI2025

Evolution of Magnetohydrodynamic Turbulence in the Expanding Solar Wind: Residual Energy and Intermittency

We conduct 3D magnetohydrodynamic simulations of decaying turbulence in the context of the solar wind. To account for the spherical expansion of the solar wind, we implement the expanding box model. The initial turbulence comprises uncorrelated counterpropagating Alfvén waves and exhibits an isotropic power spectrum. Our findings reveal the consistent generation of negative residual energy whenever nonlinear interactions are present, independent of the normalized cross helicity σc and compressibility. The spherical expansion facilitates this process. The resulting residual energy is primarily distributed in the perpendicular direction, with S 2 (b) − S 2 (u) ∝ l ⊥ or equivalently –E r ∝ k$^{–2}_{⟂}$. Here S 2 (b) and S 2 (u) are second-order structure functions of magnetic field and velocity respectively. In most runs, S 2 (b) develops a scaling relation S 2 (b) ∝ $l$$^{1/2}_{⊥}$ (E b ∝ k$^{–3/2}_{⟂}$). In contrast, S 2 (u) is consistently shallower than S 2 (b), which aligns with in situ observations of the solar wind. We observe that the higher-order statistics of the turbulence, which act as a proxy for intermittency, depend on the initial σc and are strongly affected by the expansion effect. Generally, the intermittency is more pronounced when the expansion effect is present. Finally, we find that in our simulations, although the negative residual energy and intermittency grow simultaneously as the turbulence evolves, the causal relation between them seems to be weak, possibly because they are generated on different scales.

Interplanetary turbulence↗
DOE OSTI2020

Thermal Energy Storage Using Solid Particles for Long-Duration Energy Storage

The rapid growth of renewable energy increases the importance of economically firming the electricity supply from variable solar photovoltaic- and wind-power generators. Energy storage will be the key to manage variability and to bridge the generation gap over time scales of hours or days for high renewable grid integration. The integration of renewable power and storage of excess electricity has several significant and positive impacts including: 1) expanding the renewable energy portion of total electricity generation, 2) improving the peak-load response, and 3) coordinating the electricity supply and demand over the grid. Long-duration energy storage can potentially complement the reduction of fossil-fuel baseload generation that otherwise would risk grid security when a large portion of grid power comes from variable renewable sources. Several energy storage methods are deployed or under development, including mechanical, chemical or electrochemical, and thermal energy storage (TES). Comparing their economic potential for different scales and applications helps identify suitable technology to support high renewable grid integration. Despite the progress of TES technologies developed and deployed with concentrating solar power (CSP) systems, TES has been undervalued for its potential role in electric energy storage. This paper introduces TES methods applicable to grid energy storage and particularly focuses on solid-particle-based TES to serve the purpose of long-duration energy storage (LDES). The objective of this paper is to present a standalone particle-based TES system for electric storage and to show the potential of TES systems for LDES applications over other energy storage methods such as batteries, compressed-air energy storage, or pumped-storage hydropower.

27 ARPA - Advanced Research Projects Agency-Energy↗
DOE OSTI2023

Development of Analysis Methods that Integrate Numeric and Textual Equipment Reliability Data

Within the Light Water Reactor Sustainability (LWRS) program, the Risk-Informed Systems Analysis (RISA) Pathway is performing collaborative research on the development and deployment of technologies designed to assist operating nuclear power plants (NPPs) to reduce operating costs improve plant reliability and availability. One of the RISA research areas is focusing on the development of methods and tools designed to optimize plant operations (e.g., maintenance/replacement schedules, optimal maintenance postures for plant structures, systems, and components [SSCs]) in a manner that is more cost effective than current approaches and makes better use of available SSC health data. The Risk-Informed Asset Management (RIAM) project targets this research area by creating a direct bridge between component equipment reliability (ER) data and system engineer decision making regarding maintenance activity scheduling and component aging management. In this respect, one challenge that NPP system engineers are facing is that the amount of ER data being continuously generated is not only extremely large in size, but it comes in different forms: textual (e.g., condition or maintenance reports) and numeric (e.g., generated by monitoring systems). All these data elements provide them with valuable insights and information regarding: 1) the discovery of anomalous behaviors or degradation trends, 2) the identification of the possible causes behind such behaviors/trends, and 3) the prediction of their direct consequences. However, several challenges have proved to be roadblocks to this process. While some of these challenges are technical in nature (i.e., data are often distributed over several physical servers/databases), others are conceptual in nature: data elements come in different formats (e.g., numeric or textual), and measured values have different scales (e.g., vibration spectra and oil temperature). The activities performed by the RIAM project during FY23 directly tackles the need to simultaneously integrate the analysis of ER data in all its forms, numeric and textual. Note that such task has never been performed before due to the complexity of the systems under consideration but, most importantly, because of the technical challenges behind the harmonization of ER data formats and the lack of adequate computational methods to analyze them. Our approach borrows ideas and concepts from the medical field where integration of several data sources is vital to assist medical practitioners to perform correct diagnosis and indicate optimal treatments. In our view a NPP asset is equivalent to a patient in a medical context. The main difference is the complexity of a human body is a magnitude more complex when compared to typical assets commonly present in NPPs (e.g., centrifugal pumps, or motor operated valves). This simplifies our first requirement when analyzing heterogenous ER data formats: to put data into “context”. Context is here intended as the additional piece of information that is needed by ER data analysis tools to understand what these data elements are referring to, i.e., which king of knowledge they are generating. In our context, this knowledge can be translated into models that capture the form and functional architecture of assets/systems, their dependencies, and how they interact. These models actually emulate the knowledge that that NPP system engineers possess about assets and systems; this is their key of success when analyzing ER data, their challenge is ability to handle large amount of data. Here, we employ model-based system engineering (MBSE) models of systems and assets to represent and capture their architecture and functional, i.e. cause-effect, relations. Then, ER data elements are processed by identifying first of all which elements of the developed MBSE elements they are referring to. For numeric ER data this task is fairly easy since it is possible to precisely pinpoint what MBSE elements the corresponding sensor are observing (e.g., bearing temperature of a centrifugal pump). Task is much harder for textual data since the information contained in issue or maintenance reports needs to “be understood” by a computational tool. Here we called this process as “knowledge extraction”. Once again, we borrow the experience in the medical field where methods to extract knowledge from textual data have been developed in the past decade. The missing element for us is the availability of a complete dictionary of NPP related concepts (in addition to the MBSE models presented earlier) that can put “text into context”. In FY23, such dictionary has been developed along with all the computational elements required for knowledge extraction. Lastly, once numeric and textual ER data elements have been processed and “understood”, then the last step is the discovery of possible cause-effect relations among them. This is performed by observing if a logical connection through the MBSE models exists, and if the

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Modularity-at-scale for advanced reactors presentation

Modularity options have been limited for traditional nuclear energy deployment due to the conventional light water reactor (LWR) safety requirements, such as high pressure retaining heavy and robust containment structures. However, a relatively new regulatory approach called ‘Functional Containment’ has potential to allow less expensive and more flexible designs for non-LWRs. Functional containment provides flexibility in design and deployment based on risk-informed and performance-based criteria, so that reactors are not over-designed. Non-nuclear industry has successfully used modular design approaches in automotive, aerospace, chemical processing, building construction, and ship building. These industries have shown that modular construction reduces construction time by around 30% - 50% compared to the conventional stick-built approach. The nuclear industry can use similar approaches to reduce construction time and costs, balanced with safety requirements, using the functional containment approach. This paper discussing the background of modularity in nuclear energy, examples of less learned, modularity approaches in non-nuclear industries and the potential of cost and schedule savings through the emerging regulatory design flexibilities potentially enabling combination of modular deployment at different scales.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

MODULARITY-AT-SCALE FOR COST COMPETITIVE DEPLOYMENT OF NUCLEAR ENERGY

Modularity options have been limited for traditional nuclear energy deployment due to the conventional light water reactor (LWR) safety requirements, such as high pressure retaining heavy and robust containment structures. However, a relatively new regulatory approach called ‘Functional Containment’ has potential to allow less expensive and more flexible designs for non LWRs. Functional containment provides flexibility in design and deployment based on risk informed and performance-based criteria, so that reactors are not over-designed. Non-nuclear industry has successfully used modular design approaches in automotive, aerospace, chemical processing, building construction, and ship building. These industries have shown that modular construction reduces construction time by around 30% - 50% compared to the conventional stick built approach. The nuclear industry can use similar approaches to reduce construction time and costs, balanced with safety requirements, using the functional containment approach. This paper discussing the background of modularity in nuclear energy, examples of less learned, modularity approaches in non-nuclear industries and the potential of cost and schedule savings through the emerging regulatory design flexibilities potentially enabling combination of modular deployment at different scales.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
NASA NTRS2008

Plasmaspheric Plumes Observed by the CLUSTER and IMAGE Spacecraft

Global IMAGE/EUV observations have revealed complex changes in plasmaspheric structures as the plasmasphere responds to geomagnetic activity while remaining under varying degrees of influence by co-rotation, depending on the radial distance. The complex plasmaspheric dynamics, with different scales of variability, is clearly far from being well understood. There is now renewed interest in the plasmasphere due to its apparent connections with the development of the ring current and radiation belt, and loss of ionospheric plasmas. Early in the mission, the Cluster spacecraft only crossed the plasmapause (L - 4) occasionally and made measurements of the outer plasmasphere and plasmaspheric drainage plumes. The study by Darrouzet et al. [2006] provided detailed analyses of in situ Cluster observations and IMAGE EUV observations of three plasmaspheric plumes detected in April-June, 2002. Within the next couple of years, Cluster orbit will change, causing perigee to migrate to lower altitudes, and thus providing excellent opportunities to obtain more detailed measurements of the plasmasphere. In this paper, we report our analyses of the earlier Cluster-IMAGE events by incorporating the different perspectives provided by the IMAGE Radio Plasma Imager (RPI) observations. We will discuss our new understanding of the structure and dynamics of the Cluster-IMAGE events.

Fung, S. F.↗
NASA NTRS2013

Dynamics and Control of a Disordered System in Space

In this paper, we present some ideas regarding the modeling, dynamics and control aspects of granular spacecraft. Granular spacecraft are complex multibody systems composed of a spatially disordered distribution of a large number of elements, for instance a cloud of N grains in orbit, with N greater than 10(exp 3). These grains can be large (Cubesat-size) or small (mm-size), and can be active, i.e., a fully equipped vehicle capable sensing their own position and attitude, and enabled with propulsion means, or entirely passive. The ultimate objective would be to study the behavior of the single grains and of large ensembles of grains in orbit and to identify ways to guide and control the shape of a cloud composed of these grains so that it can perform a useful function in space, for instance, as an element of an optical imaging system for astrophysical applications. This concept, in which the aperture does not need to be continuous and monolithic, would increase the aperture size several times compared to large NASA observatories such as ATLAST, allowing for a true Terrestrial Planet Imager that would be able to resolve exo-planet details and do meaningful spectroscopy on distant world. In the paper, we address the modeling and autonomous operation of a distributed assembly (the cloud) of large numbers of highly miniaturized space-borne elements (the grains). A multi-scale, multi-physics model is proposed of the dynamics of the cloud in orbit, as well as a control law for cloud shape maintenance, and preliminary simulation studies yield an estimate of the computational effort, indicating a scale factor of approximately N(exp 1.4) as a function of the number of grains. A granular spacecraft can be defined as a collection of a large number of space-borne elements (in the 1000s) designed and controlled such that a desirable collective behavior emerges, either from the interactions among neighboring grains, and/or between the grains and the environment. In this paper, each grain is considered to be a highly miniaturized spacecraft which has limited size and mass, hence it has limited actuation, limited propulsive capability, limited power, limited sensing, limited communication, limited computational resources, limited range of motion, limited lifetime, and may be expendable. The modeling and dynamics of clouds of vehicles is more challenging than with conventional vehicles because we are faced with a probabilistic vehicle composed of a large number of physically disconnected vehicles. First, different scales of motion occur simultaneously in a cloud: translations and rotations of the cloud as a whole (macro-dynamics), relative rotation and translation of one cloud member with respect to another (meso-dynamics), and individual cloud member dynamics (micro-dynamics). Second, the control design needs to be tolerant of the system complexity, of the system architecture (centralized vs. decentralized large scale system control) as well as robust to un-modeled dynamics and noise sources. Figure 1, top left, shows the kinematic parameters of a 1000 element cloud in orbit. The motion of the system is described with respect to a local vertical-local horizontal (LV-LH) orbiting reference frame (x,y,z)=F(sub ORF) of origin O(sub ORF) which rotates with mean motion omega and orbital semi-major axis R(sub 0). The orbital geometry at the initial time is defined in terms of its six orbital elements, and the orbital dynamics equation for point O(sub ORF) is propagated forward in time under the influence of the gravitational field of the primary and other external perturbations, described below. The origin of this frame coincides with the initial position of the center of mass of the system, and the coordinate axes are z along the local vertical, x toward the flight direction, and y in the orbit normal direction. The assumptions we used to model the dynamics are as follows: 1) The inertial frame is fixed at Earth's center. 2) The orbiting Frame ORF follows Keplerian orbit. 3) the cloud system dynamics is referred to ORF. 4) the attitude of each grain uses the principal body frame as body fixed frame. 5) the atmosphere is assumed to be rigidly rotating with the Earth. Regarding the grains forming the cloud: 1) each grain is modeled as a rigid body; 2) a simple attitude estimator provides attitude estimates, 3) a simple guidance logic commands the position and attitude of each grain, 4) a simple local feedback controller based on PD control of local states is used to stabilize the attitude of the vehicle. Regarding the cloud: 1) the cloud as a whole is modeled as an equivalent rigid body in orbit, and 2) an associated graph establishes agent connectivity and enables coupling between modes of motion at the micro and macro scales; 3) a simple guidance and estimation logic is modeled to estimate and command the attitude of this equivalent rigid body; 4) a cloud shape maintenance controller is based on the dynamics of a stable virtual truss in the orbiting frame. Regarding the environmental perturbations acting on the cloud: 1) a non-spherical gravity field including JO (Earth's spherical field) zonal component, J2 (Earth's oblateness) and J3 zonal components is implemented; 2) atmospheric drag is modeled with an exponential model; 3) solar pressure is modeled assuming the Sun is inertially fixed; and 4) the Earth's magnetic field is model using an equivalent dipole model. The equations of motion are written in a referential system with respect to the origin of the orbiting frame and the state is propagated forward in time using an incremental predictor-corrector scheme. A representative cloud with varying number of grains is simulated to identify the limitations in computation time as the number of grains grows. We derive a control law to track a desired surface in the ORF (equivalently to maintain a reference cloud shape) by defining an error from a desired surface shape, and designing a control law that is exponentially stable and reduces the tracking error to zero. Figure 1 (top right) shows a comparison of various requirements for simulation of single spacecraft vs. granular spacecraft, indicating the high degree of complexity that needs to be taken into consideration. The ORF components of control force required by one of the grains is, for this particular case, in the micro-Newton range. However, no attempt has been made yet to reconfigure (or re-orient) the cloud configuration internally, for which forces in the milli-Newton level are expected, depending on the time required to do the reconfiguration. Figure 1, bottom, shows the computation time as a function of the number of grains, indicating an order N(exp 1.43) scaling on a 8 Gb, 1067 MHz RAM MacOSX computer with a 3.06 GHz Intel Core 2 Duo processor. With this metric, the same simulation for a system of N=1000 grains would take 5.4 hours, and 146 hours (i.e., 6 days) for a system with N=10,000 grains. Therefore, efficient ways to simulate this complex system, where not only the time scales of natural system dynamics, but also the sampling times of the Guidance, Navigation, and Control are included, remain to be explored. Additional details on the cloud modeling, dynamics, and control will be described in the paper.

simulation↗
DOE OSTI2025

Solute Redistribution During Aging in a Laser-Based Powder Bed Fusion-Produced Fe-18Ni-9Co-5Mo-Ti Maraging Steel

Solute distribution during the laser-based powder bed fusion of metals process is highly complex because of solidification. This study builds on a previous investigation of how solute heterogeneities across different scales influence phase distribution and aging physical phenomena in a Fe-18Ni-9Co-5Mo-Ti maraging steel (Santana et al. in Addit Manuf 104494, 2024). Previous work using energy-dispersive X-ray spectroscopy revealed heterogeneous Ni and Mo enrichment at intercellular regions. Using X-ray diffraction, 8 ± 3% retained austenite was quantified in as-built samples. Additionally, electron backscatter diffraction analysis exposed the location of retained austenite along cell boundaries. In this work, scanning transmission electron microscopy results revealed the presence of non-enriched retained austenite in the as-built condition. With increased aging temperature, Ni and Mo enrichment in austenite increased, especially at the austenite/martensite boundary. Solute diffusion was simulated in DICTRA Thermo-Calc® software to assess solute redistribution between the matrix and second phases and to comprehend whether elemental diffusion processes are interconnected. Three different scenarios were considered: (1) solute diffusion in the matrix; (2) diffusion with enriched austenite in the intercellular region; (3) diffusion with non-enriched austenite in the intercellular region. Results showed that Ni partitioning controls austenite growth and that boundary regions exhibit higher element enrichment than the core.

42 ENGINEERING↗
DOE OSTI2026

Impact of Extreme Heat on Emergency Department Admissions for Childhood and Adult Asthma: An Evaluation of Earth Observations and Heat Wave Definitions

Extreme heat has been associated with adverse health outcomes, yet its impact on asthma exacerbations remains understudied. This is, in part, due to data limitations: research that relies on weather station records and aggregated health statistics cannot resolve fine-scale differences in heat impacts. This study investigates the association between heat wave definitions and summertime asthma-related emergency department visits in Baltimore, Maryland from 2016 to 2022, including 819 adult and 695 pediatric exacerbations. Using geocoded electronic health records and air temperature measurements at several spatial resolutions, we applied a case-crossover design with conditional logistic regressions at the census block group and tract levels. We found strong associations between asthma exacerbations and nighttime heat wave definitions based on relative thresholds of minimum temperatures when census block group or tract level temperature estimates were used. These relationships were significant for both age groups and showed elevated risks in socially vulnerable areas. In contrast, heat wave definitions derived from the city's primary National Weather Service synoptic weather station show associations between asthma and daytime heat extremes, suggesting that the character of the heat hazard depends on the scale at which it is defined. The extreme heat event definition used by Baltimore City's Code Red system showed no significant association with exacerbations. These findings highlight the importance of data resolution in shaping health inferences related to extreme heat in urban environments. Further, this study demonstrates that, regardless of spatial scale, extreme heat is associated with asthma exacerbations in both age groups.

Corpuz, B. [Johns Hopkins University, Baltimore, M↗
NASA NTRS1995

A similarity hypothesis for the two-point correlation tensor in a temporally evolving plane wake

The analysis demonstrated that the governing equations for the two-point velocity correlation tensor in the temporally evolving wake admit similarity solutions, which include the similarity solutions for the single-point moment as a special case. The resulting equations for the similarity solutions include two constants, beta and Re(sub sigma), that are ratios of three characteristic time scales of processes in the flow: a viscous time scale, a time scale characteristic of the spread rate of the flow, and a characteristic time scale of the mean strain rate. The values of these ratios depend on the initial conditions of the flow and are most likely measures of the coherent structures in the initial conditions. The occurrences of these constants in the governing equations for the similarity solutions indicates that these solutions, in general, will only be the same for two flows if these two constants are equal (and hence the coherent structures in the flows are related). The comparisons between the predictions of the similarity hypothesis and the data presented here and elsewhere indicate that the similarity solutions for the two-point correlation tensors provide a good approximation of the measures of those motions that are not significantly affected by the boundary conditions caused by the finite extent of real flows. Thus, the two-point similarity hypothesis provides a useful tool for both numerical and physical experimentalist that can be used to examine how the finite extent of real flows affect the evolution of the different scales of motion in the flow.

Ewing, D. W.↗
DOE OSTI2024

DECOVALEX-2023: Task F1 Final Report

DECOVALEX-2023 Task F is a comparison of models and methods for post-closure performance assessment (PA) of a deep geologic repository for radioactive waste. The general aims of Task F are to build confidence in the models, methods, and software used for PA and to stimulate additional research and development in PA methodologies. The task objectives are to motivate development of PA modelling skills and capabilities, to examine the influence of model choices on calculated repository performance, and to compare the uncertainties introduced by model choices to other sources of uncertainty. Task F involves no actual experiment or site. It is a PA modelling exercise that requires the conceptual development of hypothetical repository designs and geologic settings. Because three of the teams were interested in salt and the rest of the teams were interested in crystalline rock, Task F was split into two branches: Task F1 for crystalline rock and Task F2 for salt. This report is for Task F1, crystalline rock. Teams from seven countries (Canada, Czech Republic, Germany, Korea, Sweden, Taiwan, and United States) participated in Task F1. The teams worked together to define the features, events, and processes of the reference case repository and established a set of performance measures. In addition, they defined a set of benchmark problems designed to test and compare modelling capabilities for fracture flow and transport at different scales. The repository design and benchmark problems are documented in a Task Specification that evolved over time as the group honed the specifications. The benchmark problems verified that each team can aptly model flow and transport in fractured media in 1-, 2-, and 3-dimensions. Two general approaches were used for the 3-dimensional benchmarks: discrete fracture network (DFN) and equivalent continuous porous medium (ECPM). DFN modelling involves explicit meshing of each fracture while ECPM modelling aims to capture the effective porosity and directional permeability of each cell in a space-filling mesh as affected by intersecting fractures. In some models, a combination of the two is used, i.e., DFN for large known fractures and ECPM for the rest of the domain. Transport is solved by using either the advection-dispersion equation or particle tracking. Although some variation is observed among model breakthrough curves in the benchmark problems, there is strong agreement in breakthrough behaviour up to at least the 75 th percentile for all benchmarks. At the 90 th percentile, breakthrough results show larger differences, suggesting several models retain substantially higher fractions of tracer in regions of slower moving water. In addition to the flow and transport benchmarks, several teams completed the source term benchmark, verifying capabilities for modelling radionuclide decay and ingrowth, waste package breach, instant release fractions, fuel matrix degradation rates, and radionuclide solubility limitations. The reference case is conceptualized as a generic spent fuel repository at a depth of 450 m in fractured crystalline rock. The repository has 50 parallel backfilled drifts, each with 50 deposition holes 6 m apart. Each deposition hole contains a 4-PWR waste package and bentonite buffer. The rock domain is 5 km in length, 2 km in width, and 1 km in depth. It has 6 deterministic fractured deformation zones and a multitude of stochastic fractures. Teams generally used the ECPM approach for the entire rock or a hybrid approach in which the deterministic fracture zones are modelled with a DFN and the rest of the rock is modelled by ECPM. Of the reference case problems specified, only the results of the initial reference case problem are compared in this report. The initial problem focuses on transport from the deposition holes to the surface, i.e., it neglects waste package performance. Tracers are released at all waste package locations at time zero and tracked for their releases to the near field and ground surface. The water fluxes calculated at the ground surface entry and exit regions of the domain are similar for all models except for two that have considerably lower fluxes. For tracer transport, large differences are observed among models in the magnitude of tracer transported. Much of the difference appears to be due to how the repository is implemented and hence the different degrees of repository simplification. Models that exclude the drifts, buffer, and backfill from the domain tend to show greater release of tracers and radionuclides from the repository. The initial study presented here indicates that major differences in modelling important processes within the repository (e.g., diffusion through buffer and backfill) can produce broadly different release and transport results, especially when those processes are excluded. Even for the models that included all specified features, events, and processes, the results show significant differences and demonstrate the importance of examining multiple modelling approaches in performance assessment. The differences in results observed in this study are expected to motivate teams to either increase complexity in future versions of the reference case models or to improve methods to account for the effects of simplified features and processes. Either way, future improvements in these models are expected to produce results that more closely agree.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

Out-of-equilibrium interactions and collective locomotion of colloidal spheres with squirming of nematoelastic multipoles

Many living and artificial systems show similar emergent behavior and collective motions on different scales, starting from swarms of bacteria to synthetic active particles, herds of mammals, and crowds of people. What all these systems often have in common is that new collective properties like flocking emerge from interactions between individual self-propelled or driven units. Such systems are naturally out-of-equilibrium and propel at the expense of consumed energy. Mimicking nature by making self-propelled or externally driven particles and studying their individual and collective motility may allow for deeper understanding of physical underpinnings behind collective motion of large groups of interacting objects or beings. Here, in this work, using a soft matter system of colloids immersed into a liquid crystal, we show that resulting so-called nematoelastic multipoles can be set into a bidirectional locomotion by external oscillating electric fields. Out-of-equilibrium elastic interactions between such colloidal objects lead to collective flock-like behaviors emerging from time-varying elasticity-mediated interactions between externally driven propelling particles. Repulsive elastic interactions in the equilibrium state can be turned into attractive interactions in the out-of-equilibrium state under applied external electric fields. We probe this behavior at different number densities of colloidal particles and show that particles in dense dispersions collectively select the same direction of a coherent motion due to elastic interactions between near neighbors. In our experimentally implemented design, their motion is highly ordered and without clustering or jamming often present in other colloidal transport systems, which is promising for technological and fundamental-science applications, like nano-cargo transport, out-of-equilibrium assembly, and microrobotics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗