SEARCH · Search NASA
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.
Well-balanced reconstruction techniques have been developed for stellar hydrodynamics to address the challenges of maintaining hydrostatic equilibrium during evolution. I show how to adapt a simple well-balanced method to the piecewise parabolic method for hydrodynamics. A python implementation of the method is provided.
The random ray method is a recently developed stochastic method for solving neutral particle transport problems based on the method of characteristics. Perhaps surprisingly for a characteristics-based method using flat sources, we note that the random ray method can produce negative fluxes which may be numerically troublesome in several situations. These occur most severely in fixed source problems where the source is in a region with a small cross section. Additionally, we briefly discuss another source of bias which can occur in similar situations, namely a ray missing a mesh with a strong source and small cross section, resulting in the entirety of the source being unphysically deposited locally. This paper describes the mechanism by which negative fluxes may occur and several different methods to mitigate their effects. These fixes are tested on an eigenvalue problem, a ‘fusion-like’ shielding problem, and a shielding problem featuring an adjoint calculation. Even when extremely coarse random ray quadratures are used such that 20%–30% of cells are missed during a given iteration, use of the preferred fix technique ensures local flux tally errors remain trivial (below 1%). The preferred fix is now the default option in SCONE and OpenMC.
Herein, we present a new method for solving the linear Boltzmann transport equation. Two commonly used and well-understood methods for solving partial differential equations are the method of characteristics (MOC) and the finite element method (FEM). We propose a new method that combines the fundamental concept of the FEM with the analytic solution from the MOC to obtain coefficients for the FEM basis function expansion. Traditionally, coefficients for the FEM basis function expansion are obtained via matrix inversion. Instead, we solve for the coefficients with the MOC and represent the underlying fields with the basis function expansion using these coefficients. We provide a convergence study for our method with results from two sets of FEM basis functions: Gauss-Legendre and Gauss-Lobatto sets. We also compare two different variations of our method categorized as short characteristics and intermediate characteristics.
Earthquake nowcasting has been proposed as a means of tracking the change in large earthquake potential in a seismically active area. The method was developed using observable seismic data, in which probabilities of future large earthquakes can be computed using Receiver Operating Characteristic methods. Furthermore, analysis of the Shannon information content of the earthquake catalogs has been used to show that there is information contained in the catalogs, and that it can vary in time. So an important question remains, where does the information originate? In this paper, we examine this question using stochastic simulations of earthquake catalogs. Our catalog simulations are computed using an Earthquake Rescaled Aftershock Seismicity (“ERAS”) stochastic model. This model is similar in many ways to other stochastic seismicity simulations, but has the advantage that the model has only 2 free parameters to be set, one for the aftershock (Omori-Utsu) time decay, and one for the aftershock spatial migration away from the epicenter. Generating a simulation catalog and fitting the two parameters to the observed catalog such as California takes only a few minutes of wall clock time. While clustering can arise from random, Poisson statistics, we show that significant information in the simulation catalogs arises from the “non-Poisson” power-law aftershock clustering, implying that the practice of de-clustering observed catalogs may remove information that would otherwise be useful in forecasting and nowcasting. We also show that the nowcasting method provides similar results with the ERAS model as it does with observed seismicity.
The U.S. nuclear energy expansion goals are driving the demand for manufacturing routes that can rapidly produce large, complex, near net shape components. Powder metallurgy hot isostatic pressing (PM HIP) is an advanced manufacturing technique that can be economically scaled-up, while alleviating the supply chain challenges that forging and casting face in terms of cost and lead time constraints. This makes PM-HIP a viable technology to aid and accelerate large-scale part manufacturing for nuclear applications. However, large-scale qualification and deployment of this technology require a thorough understanding of the influence of powder feedstock quality, powder handling history, hot isostatic pressing (HIP) parameters, and subsequent heat treatment on microstructural evolution and elevated temperature mechanical performance. The present work focusses on 316 austenitic stainless steel (SS316) which is most commonly used in high temperature environments for nuclear applications Results from this study show that PM HIPed SS316 meets ASME tensile requirements at room temperature and at elevated temperature. However, creep performance of PM-HIPed SS316 remains inferior to its wrought counterpart, demonstrating that tensile performance alone is not a reliable metric for long duration high temperature integrity. Further, this report delineates powder derived microstructural features that govern creep damage, with key evidences pointing to “microstructural inheritance” from gas atomized powder feedstocks. Commercial SS316 powders of varying chemical compositions and recycling histories were studies, and the results showed large differences in elemental segregation, oxide surface layers and secondary phase distributions. Multi-scale characterization revealed segregation of chromium, molybdenum, manganese and silicon at the boundaries and the precipitation of manganese-, silicon-, and molybdenum-oxides. During HIP consolidation, these surface oxides transform into decorated prior particle boundaries (PPBs) and grain boundary inclusions that persist through conventional post-HIP solution annealing treatment. The retained oxides in post-HIP microstructures were found to influence grain growth, precipitation behavior, and ultimately creep cavitation and fracture. Such post-HIP heat treatments are therefore limited by a complex trade-off between grain growth, and oxide coarsening which aggravate creep damage by acting as nucleation sites for cavities. The objective of this work is to establish an integrated processing–structure–property framework for PM-HIP 316 stainless steel by investigating the influence of powder feedstock characteristics in pre- and post-HIP processing as well as to understand the significance of post-HIP heat treatment on microstructural evolution and creep properties. The results from this report emphasize the significance of powder feedstock integrity in improving creep performance of PM-HIPed SS316, by highlight the effect of rapid solidification, elemental segregation, oxide formation and powder recycling on microstructural defect inheritance following HIP consolidation. Rather than considering HIP processing, solution annealing, and mechanical performance independently, this report treats powder production, HIP consolidation, post-HIP thermal processing, and creep deformation as interconnected stages within a continuous metallurgical process. The resulting framework will provide a scientific basis for developing feedstock engineering strategies capable of improving long-term reliability of PM-HIP stainless steels and accelerating their qualification for advanced nuclear applications.
Whey protein isolates (WPI) are widely used in processed foods for their versatile functional properties. Modifying the structural properties of proteins by assembling them into mesoscale or microscale particles may improve their functionality and broaden their applications. This study aims to manufacture and characterize mesoscale whey protein particles (WPP) derived from WPI. Two types of WPP, WPP1 (0.05 mL/min) and WPP2 (0.25 mL/min), were prepared through a multistep approach involving liquid antisolvent (LAS) precipitation, heat treatment, and microfluidization. Liquid antisolvent precipitation was performed by injecting a 20% (wt/vol) WPI dispersion (pH 7) into an ethanol-glycerol mixture (75:25, vol/vol) under laminar flow, followed by heat treatment at 80°C for 20 min as a particle hardening step. This process produced stable fiber- and ribbon-shaped whey protein assemblies (WPA), which served as precursors to WPP. Subsequent microfluidization (150 MPa, 6 passages) reduced the size of WPA, yielding mesoscale WPP with irregular morphologies and a more uniform size distribution, as revealed by microscopy and dynamic light scattering. ζ-Potential and fluorescence labeling indicated higher surface charge and surface hydrophobicity of WPP compared with untreated WPI. The WPP showed internal mass fractal and surface fractal structures at larger length scales, analyzed using small-angle X-ray scattering. Fourier transform infrared spectroscopy demonstrated an increased fraction of intermolecular β-sheets in WPP, suggesting that hydrogen bonding contributed to their formation. Gel electrophoresis confirmed that disulfide bonds served as the primary cross-links stabilizing the WPP structure. Furthermore, turbidity measurements showed that WPP exhibited superior colloidal phase stability compared with untreated WPI and maintained high colloidal stability under both acidic and neutral pH conditions.
This report, Advancing Building Air Leakage Measurement and Modeling: Measurement of Interior Air Leaks, Experimental Comparison to Conventional Methods, and Building Modeling Characteristics, developed two new building air leakage measurement methods, experimentally evaluated the performance in comparison to conventional methods, experimentally evaluated the behavior of leaks for specific crack structures, and identified testing constraints and model implications.
Explore the source record for details and available documents.
Control blade cusping can introduce significant error in boiling water reactor (BWR) calculations with MPACT when blade tips fall partway within an axial method of characteristics (MOC) plane, requiring homogenization of controlled and uncontrolled regions. This work implements subplane decusping for BWRs in MPACT by enabling BWR-compatible subplane coarse mesh finite difference (CMFD) and extending the existing decusping framework to represent between-assembly control blades that insert from the bottom of the core. The method resolves axial heterogeneity on a refined subplane mesh in the low-order solve and uses the resulting subplane flux shape to form flux-volume homogenized transport cross sections for the partially rodded MOC plane. The capability is evaluated using a single physics General Electric (GE)-14 assembly and a multiphysics Peach Bottom Unit 2 Type 1 assembly (PB2T1A) with thermal hydraulic feedback. In both cases, coarse axial meshes with and without subplane decusping are compared against fine mesh reference solutions over the full range of blade withdrawal positions. Subplane decusping reduces maximum/average eigenvalue errors from 3,369/220 pcm to 303/36 pcm for GE-14 and from 6,689/585 pcm to 1,582/109 pcm for PB2T1A. Additionally, it reduces pin power root mean square errors from 5.2%/0.8% to 1.8%/0.2% for GE-14 and from 6.9%/0.6% to 2.3%/0.1% for PB2T1A. These results demonstrate an effective, practical correction for BWR blade cusping in VERA-MPACT.
The dynamic behaviour of floating offshore wind turbines (FOWTs) involves complex interactions of multivariate loads from wind, waves, and currents, which result in complex motion characteristics. Although methods for analysing global motion responses are well-established, the time- and location-dependent kinematics remain underexplored. This paper investigates the instantaneous centre of rotation (ICR), a point of zero velocity at a time instance of general plane motion. Understanding and strategically positioning the ICR can reduce the dynamic motion in critical structural locations, enhancing the performance and structural robustness of FOWTs. The paper presents a method for computing the ICR using time-domain simulation results and proposes a statistical analysis approach suitable for design studies. Building on prior research, it examines the sensitivity of the ICR to external loading and design features, providing insights into how these factors influence motion response and how the motion response influences the statistics of the ICR, structural loads, and other performance metrics of interest. The study explores two FOWT configurations, a spar and a semisubmersible, identifying design variables that most effectively control the ICR statistics and identifying the ICR statistics most correlated with the responses of interest. Finally, through two case studies, we demonstrate how to apply these new insights in a practical design scenario. By adjusting the design variables most correlated with the ICR (fairlead vertical position and centre of mass for the spar and mooring line length and offset column diameter for the semisubmersible), we successfully modified the designs of the floating support structures to reduce the loads in the mooring lines, tower base, and blade roots, improving the ultimate strength and fatigue characteristics compared to the original designs.
In this paper, we extend existing turret defense differential game formulations involving a turn-constrained turret and mobile agent to include specified final time and a constraint. For the purposes of this analysis, the specified final time may represent some exogenous input, perhaps representing the time at which some other event will take place. As for the constraint, it represents a no-fly zone for the mobile agent. The scenario is formulated as a two-player, zero-sum differential game and solved via the method of characteristics (i.e., back-propagation of equilibrium trajectories). Three different trajectory types make up the solution: trajectories that end with the turret aligned with the mobile agent, trajectories that end with the mobile agent on the constraint boundary, and regular trajectories.
The emergence of renewable deep eutectic solvents (DES) as clean and efficient catalysts and solvents has created new opportunities for lignocellulosic biorefineries and critical material sectors, including chemical, energy, pharmaceutical, textile, and hydrometallurgical industries. This review provides an in-depth overview of DES, covering their chemistry, classifications, preparation methods, processing characteristics, and recyclability, while highlighting their unique attributes and industry relevant applications. Emphasis is placed on the integration of DES into advanced biorefinery systems, focusing on their tunable physicochemical and thermodynamic properties for biomass pretreatment and the production of value-added products. The review explores how DES can be tuned for selective dissolution of biomass components and evaluates production and valorization of DES-derived biochemicals, with attention to lignin extraction mechanisms and conversion of biomass into bioproducts and biofuels. Beyond biorefineries, the scope extends to DES applications in electrochemical energy devices, where they serve as electrolytes, synthesis media for electrode materials, and leaching agents in battery recycling. The multifunctional roles of DES in pharmaceutical, hydrometallurgical, and textile sectors are also explored for contributions to sustainable processing. Finally, the review identifies future research directions, outlining benefits, challenges, and knowledge gaps, for continued industrial development.
Bottom-up urban energy models are crucial for understanding current energy use patterns and informing design strategies. However, accurately characterizing these models to represent different communities remains a challenge due to the extensive data needed for simulating existing energy use behavior. This data includes information related to human activities and building characteristics, all of which correlate with socioeconomic factors. To overcome this challenge, we developed an automated framework that utilizes both top-down and bottom-up data, to predict unknown building and occupant characteristics that are needed for more accurate and equitable modeling and analytics. Our framework, integrated into the URBANopt district energy modeling platform, uses statistical data models from ResStock. URBANopt models co-located buildings and neighborhoods. At this scale there are data gaps in building characteristic data, such as materials, insulation, occupancy, income, and energy usage of the buildings. To address this data gap, we use ResStock data, representative at the census tract scale, and develop machine-learning and deeplearning techniques to disaggregate it to individual buildings. By mapping unique occupant, building and economic properties to URBANopt energy models, we gain detailed insights into the variability of building energy use across different neighborhoods. This insight helps deploy technologies for co-located buildings and supports targeted upgrades for communities with unique economic and demographic characteristics, ensuring energy equity. Accurate characterization of energy models allows us to develop equitable strategies tailored to diverse neighborhoods, whether underserved or affluent. Our automated framework streamlines energy modeling and provides a reliable tool for building energy characterization.
Inverter-based resources (IBRs) such as photovoltaics (PVs), wind turbines, and battery energy storage systems (BESSs) are widely deployed in low-carbon power systems. However, these resources typically do not provide the inertia needed for grid stability, resulting in a low-inertia power system. IBRs and lack of inertia have been known to cause anomalies such as waveform distortions and wideband oscillations in power systems due to the limited inertia level, leading to increased generation trips and load shedding. Here, to achieve effective anomaly identification, this paper proposes a synchro-waveform-based algorithm utilizing real-time synchronized voltage waveform measurements from waveform measurement units (WMUs). In the proposed method, different physical characteristics, as well as statistical features, are extracted from synchronized voltage waveform measurements to filter anomalies. Then, the anomaly identification approach based on the random forest is developed and deployed into the FNET/GridEye system considering trade-offs among accuracy, computational burden, and deployment cost. Moreover, four WMUs are specially designed and deployed on Kauai Island to receive instantaneous synchronized voltage waveform measurements. To verify the performance of the proposed algorithm, different experiments are carried out with collected field test data. The result demonstrates that the performance of the proposed synchro-waveform-based anomaly categorization algorithm can accurately identify anomalies 95.35% of the time, which has comparable performance among benchmarking algorithms.
Item purpose: SnT virtual poster showing LYNM work Internal peer review status: This poster has been peer reviewed HQ review status: this poster has been approved by HQ
Background Age-related changes in adipose tissue affect chronic medical diseases and mobility disability but mechanism remains poorly understood. The goal of this study is to define methods for phenotyping unique characteristics of adipose tissue from older adults. Methods Older adults enrolled in study of muscle, mobility, and aging selected for the adipose tissue ancillary (SOMMA-AT; N = 210, 52.38% women, 76.12 ± 4.37 years) were assessed for regional adiposity by whole-body magnetic resonance (AMRA) and underwent a needle-aspiration biopsy of abdominal subcutaneous adipose tissue (ASAT). ASAT biopsies were flash frozen, fixed, or processed for downstream applications and deposited at the biorepository. Biopsy yields, qualitative features, adipocyte sizes, and concentration of adipokines secreted in ASAT explant conditioned media were measured. Inter-measure Spearman correlations were determined. Results Regional, but not total, adiposity differed by sex: women had greater ASAT mass (8.20 ± 2.73 kg, p < .001) and biopsy yield (3.44 ± 1.81 g, p < .001) than men (ASAT = 5.95 ± 2.30 kg, biopsy = 2.30 ± 1.40 g). ASAT mass correlated with leptin (r = 0.54, p < .001) and not resistin (p = .248) and adiponectin (p = .353). Adipocyte area correlated with ASAT mass (r = 0.34, p < .001), BMI (r = 0.33, p < .001), adiponectin (r = −0.22, p = .005) and leptin (r = 0.18, p = .024) but not with resistin (p = .490). Conclusion In addition to the detailed ASAT biopsy processing in this report, we found that adipocyte area correlated with ASAT mass, and both measures related to some key adipokines in the explant conditioned media. In conclusion, these results, methods, and biological repositories underscore the potential of this unique cohort to impact the understanding of aging adipose biology on disease, disability, and other aging tissues.
This report explores using "generator knowledge" to determine whether a solar photovoltaic (PV) module must be managed as hazardous waste prior to recycling or landfilling. Generator knowledge is a legal term and existing regulatory pathway for making a hazardous waste determination that has been used by other industries but is a relatively unknown option to the PV industry. In the United States, a hazardous waste determination often acts as a pre-requisite to recycle or landfill a PV module. The results of the hazardous waste determination dictate whether the PV module must be managed as hazardous waste or nonhazardous solid waste. Managing a PV module as hazardous waste requires compliance with stringent U.S. federal and state hazardous waste law. In addition to increased management costs, which can be ten times higher, legal liability for PV modules regulated as hazardous waste is also heightened with both civil and criminal penalties for noncompliance which includes making an inaccurate or faulty hazardous waste determination. The most common reason a PV module would be regulated as hazardous is if it contains a regulated metal in an amount that equals or exceeds the toxicity characteristic limits. To determine whether a PV module exhibits a hazardous characteristic, the regulated person/entity must "apply knowledge...in light of the materials and processes used." In the absence of adequate knowledge to determine whether the PV module is hazardous, it must be tested using Test Method 1311 Toxicity Characteristic Leaching Procedure (TCLP) or an equivalent EPA-approved method. Although TCLP is the predominant method used to today to make a hazardous waste determination for PV modules, evidence from this study concludes it is not practical to TCLP test every PV module even in a single utility-scale installation, and a scalable solution is needed. This study finds that knowledge-based hazardous waste determinations may allow a regulated person/entity to make a hazardous waste determination for more than one PV module at one time - making this regulatory pathway a potential scalable solution. Through legal analysis and interviews with 44 experts, the authors explore what it means to make an accurate knowledge-based hazardous determination for PV modules considering sources and forms of information as well as potential limitations. The work aims to provide a foundation for building consensus on whether knowledge-based hazardous determinations are a viable, scalable industry approach for solar.
Introduction: The distribution of sediment grain size in streams and rivers is often quantified by the median grain size (D50), a key metric for understanding and predicting hydrologic and biogeochemical function of streams and rivers. Manual D50 measurements are time-consuming and ignore larger grains, while approaches to model D50 based on catchment characteristics may over-generalize and miss site-scale heterogeneity. Machine learning-enabled object detection methods like You Only Look Once (YOLO) provides an alternative that enables estimation of D50 that is faster than manual measurements and more site-specific than predictions based on catchment characteristics. Methods: To understand the potential role of object detection methods for improving understanding of D50, we compared D50 estimates made manually, predicted from catchment characteristics, and using a YOLO-enabled approach across the Yakima River Basin. Results: We found distinct differences between methods for D50 averages and variability, and relationships between D50 estimates and basin characteristics. Discussion: We discuss the advantages and limitations of object detection methods versus current methods, and explore potential future directions to combine D50 methods to better estimate spatiotemporal variation of D50, and improve incorporation into basin-scale models.