SEARCH · Search NASA

Results for “Numerical analysis”

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 325 records · Page 18

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Analysis of the SiMPL Method for Density-Based Topology Optimization

We present a rigorous convergence analysis of a new method for density-based topology optimization that provides pointwise bound-preserving design updates and faster convergence than other popular first-order topology optimization methods. Due to its strong bound preservation, the method is exceptionally robust, as demonstrated in numerous examples here and in the companion article [D. Kim et al., Struct. Multidiscip. Optim., 68 (2025), 74]. Furthermore, it is easy to implement with clear structure and analytical expressions for the updates. Our analysis covers two versions of the method, characterized by the employed line search strategies. We consider a modified Armijo backtracking line search and a Bregman backtracking line search. For both line search algorithms, our algorithm delivers a strict monotone decrease in the objective function and further intuitive convergence properties, e.g., strong and pointwise convergence of the density variables on the active sets, norm convergence to zero of the increments, convergence of the Lagrange multipliers, and more. In addition, the numerical experiments demonstrate apparent mesh-independent convergence of the algorithm. Here, we refer to the new algorithm as the SiMPL method (pronounced “simple”), which stands for Sigmoidal Mirror descent with a Projected Latent variable.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

Survey of baffle plate samples, specimen machining, and hydrogen measurements

This report outlines preliminary results on the survey and characterization of irradiated baffle plate, baffle former, and flux thimble tube specimens as part of the ongoing investigation into irradiation-induced embrittlement in austenitic stainless steels at ambient (room temperature) conditions. The specimens originated from commercial pressurized water reactor (PWR) components, covering displacement damage doses ranging from approximately 0.065 to 75 dpa. Initial scanning electron microscopy (SEM) surveys revealed that specimen surfaces exhibited fine machining marks and in-service-formed oxide layers on the side surfaces of the analyzed specimens. The oxide layers revealed specific features resembling localized "pitting-like" corrosion. An additional set of specimens representing in-service oxidation is being preserved for future microstructure analysis work. As believed, this will provide additional insights into long-term material degradation in PRWs. Specimen machining challenges emerged due to a complex failure of the electric discharge machining (EDM) system located in hot area. Despite partial restoration, issues persist with the EDM’s secondary power supply, necessitating the exploration of alternative EDM or computer numerical control (CNC) machining approaches to facilitate tensile specimen preparation. Currently, low speed saw cutting is in progress to prepare specimens for hydrogen measurements and general microstructure analysis. The near-term goals include completing tensile specimen machining for advanced mechanical testing and characterizing fracture mechanisms, and stress-corrosion cracking testing, ultimately aiming to identify and mitigate the ambient-condition intergranular cracking through targeted post-irradiation annealing strategies.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2024

A Survey of Open-Source Tools for Transmission and Distribution Systems Research

This work presents a review of open-source electric power transmission and distribution systems analysis tools suitable for use by industry professionals and academic researchers. Due to the high complexity of the electric grid, there exist numerous tools and extensive research pertaining to nearly every aspect of the design, operation, and control of transmission and distribution networks. In addition to the commercial tools, a wide range of free, open-source tools, models, and data usable by the scientific community for related research have been developed by different organizations, including both international and US universities and national laboratories. However, due to the absence of a catalog of available tools, models and data, researchers often lack a knowledge of existing capabilities and may develop duplicative software and tools. Increasing awareness of these available resources seeks to accelerate their broader use, leading to more efficient and standardized grid analysis. This review paper (which is part of a larger survey effort that studied over 400 tools in the transmission, distribution, buildings, and electric vehicles space) outlines selected open-source resources that have been developed in power transmission and distribution systems research. It is anticipated that this work can serve as a guide for industry and academic researchers alike, ensuring that research efforts are well-channeled.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2024

A Global Increase in Nearshore Tropical Cyclone Intensification

Abstract Tropical Cyclones (TCs) inflict substantial coastal damages, making it pertinent to understand changing storm characteristics in the important nearshore region. Past work examined several aspects of TCs relevant for impacts in coastal regions. However, few studies explored nearshore storm intensification and its response to climate change at the global scale. Here, we address this using a suite of observations and numerical model simulations. Over the historical period 1979–2020, observations reveal a global mean TC intensification rate increase of about 3 kt per 24‐hr in regions close to the coast. Analysis of the observed large‐scale environment shows that stronger decreases in vertical wind shear and larger increases in relative humidity relative to the open oceans are responsible. Further, high‐resolution climate model simulations suggest that nearshore TC intensification will continue to rise under global warming. Idealized numerical experiments with an intermediate complexity model reveal that decreasing shear near coastlines, driven by amplified warming in the upper troposphere and changes in heating patterns, is the major pathway for these projected increases in nearshore TC intensification.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Absorbing boundary conditions in material point method adopting perfectly matched layer theory

This study focuses on solving the numerical challenges of imposing absorbing boundary conditions for dynamic simulations in the material point method (MPM). To attenuate elastic waves leaving the computational domain, the current work integrates the Perfectly Matched Layer (PML) theory into the implicit MPM framework. The proposed approach introduces absorbing particles surrounding the computational domain that efficiently absorb outgoing waves and reduce reflections, allowing for accurate modeling of wave propagation and its further impact on geotechnical slope stability analysis. The study also includes several benchmark tests to validate the effectiveness of the proposed method, such as several types of impulse loading and symmetric and asymmetric base shaking. The conducted numerical tests also demonstrate the ability to handle large deformation problems, including the failure of elasto-plastic soils under gravity and dynamic excitations. The findings extend the capability of MPM in simulating continuous analysis of earthquake-induced landslides, from shaking to failure.

58 GEOSCIENCES↗

DOE OSTI2025

A three-dimensional laser ray-tracing methodology for radiation-hydrodynamics simulations

We report on a methodology for performing laser ray-tracing in three spatial dimensions for radiation-hydrodynamics simulation codes. Our method, which is an extension of that developed in Haines et al., Comput. Fluids 201, 104478 (2020), utilizes an automatically generated separate mesh for the laser ray-tracing from the radiation-hydrodynamics mesh. This enables the laser mesh to be tailored to minimize ray noise with significantly fewer rays than would be required when the ray-tracing is performed on the radiation-hydrodynamics mesh, primarily by allowing the use of high-aspect-ratio cells that are not suitable for hydrodynamics solvers. For a planar target, we show that our method provides a ≈ 100× reduction in computational expense to achieve a fixed level of ray noise relative to ray-tracing directly on the radiation-hydrodynamics mesh. The relatively low ray requirement also enables efficient computation of cross-beam energy transfer. Each cell in the logically cubic laser mesh is a non-convex dodecahedron with triangular sides, and numerical integration of the ray trajectories and inverse bremsstrahlung is performed by mapping each cell to the unit cube. We will describe our methodology in detail as well as its implementation in the xRAGE radiation-hydrodynamics code, discuss performance, and present the results from applying the methodology to test problems with analytic solutions for laser ray-tracing through a quadratic density gradient with an analytic solution as well as for a laser-driven heat front. In 3D radiation-hydrodynamics simulations of laser-driven experiments performed on the National Ignition Facility, laser ray-tracing with our methodology uses less than 1% of total computational time while introducing acceptably low levels of ray noise.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

TRUST Sensors in Environments: Fiber Optic Displacement (SE-FOD) Report, Release FY24

The objective of the Delivery Environments (DE) Testbeds to Reduce Uncertainty in Simulations and Tests (TRUST) project is to quantify and help increase confidence in specific areas of computational and experimental capabilities that are applicable to the development, on-target assessment, and qualification of current and future delivery environments. More complete quantification of confidence in experimental and computational capabilities and the sufficient increase of confidence in those capabilities is critical to improving weapons engineering design, qualification, and assessment efforts that are critical to the current and future stockpile. This work uses and provides feedback on analysis tools and experimental results databases for efficient and responsive engineering which are currently under development. Under the TRUST work package, there are five testbeds and their associated engineering analysis baseline models (EABMs): 1. Contact Thermal Conductivity (CTC), 2. Nonlinear Dynamics (ND), 3. Sensors in Environments: Accelerometers (SE-A), 4. Sensors in Environments: Fiber Optic Displacement Gauges (SE-FOD), and 5. Sensors in Environments: Thermocouples (SE-TC). The TRUST project uses single-feature testbeds to quantify uncertainties in specific models and experiments and to identify capability development needs that can help to reduce these uncertainties. Each testbed is designed, configured, and tested in collaboration with groups with design and experimental capability: E-14 and MPA-CINT. The complementary simulations are conducted using W-13 analysis tools and stored in model repositories with plans for incremental progress toward EABM requirements. W-13 extends and exercises the testbed simulations in collaboration with experimentalists for uncertainty quantification of current and future materials, geometries, and environments. Additionally, TRUST is intended to provide engineers in W-13 and E-14 with experience in both numerical simulations and experimental methods through cross-discipline collaborations. Following the introduction to the TRUST project, the remainder of this report focuses on experimental, modeling and simulation, and analysis efforts conducted in FY24 relevant to the TRUST Sensors in Environments: Fiber Optic Displacement (SE-FOD) testbed.

42 ENGINEERING↗

DOE OSTI2024

Wavelength modulation laser-induced fluorescence for plasma characterization

Laser-Induced Fluorescence (LIF) spectroscopy is an essential tool for probing ion and atom velocity distribution functions (VDFs) in complex plasmas. VDFs carry information about the kinetic properties of species that is critical for plasma characterization. Accurate interpretation of these functions is challenging due to factors such as multicomponent distributions, broadening effects, and background emissions. Our research investigates the use of Wavelength Modulation (WM) LIF to enhance the sensitivity of VDF measurements. Unlike standard Amplitude Modulation (AM) methods, WM–LIF measures the derivative of the LIF signal. This approach makes variations in VDF shape more pronounced. VDF measurements with WM–LIF were investigated with both numerical modeling and experimental measurements. The developed model enables the generation of both WM and AM signals, facilitating comparative analysis of fitting outcomes. Experiments were conducted in a weakly collisional argon plasma with magnetized electrons and non-magnetized ions. Measurements of the argon ion VDFs employed a narrow-band tunable diode laser, which scanned the 4p 4 D 7/2 –3d 4 F 9/2 transition centered at 664.553 nm in vacuum. A lock-in amplifier detected the second harmonic WM signal, which was generated by modulating the laser wavelength with an externally controlled piezo-driven mirror of the diode laser. Finally, our findings indicate that the WM–LIF signal is more sensitive to fitting parameters, allowing for better identification of VDF parameters such as the number of distribution components, their temperatures, and velocities. In addition, WM–LIF can serve as an independent method to verify AM measurements and is particularly beneficial in environments with substantial light noise or background emissions, such as those involving thermionic cathodes and reflective surfaces.

47 OTHER INSTRUMENTATION↗

DOE OSTI2025

Kolmogorov-Arnold wavefunctions

Here, this work investigates Kolmogorov-Arnold network-based (KAN) wave-function Ansätz as viable representations for quantum Monte Carlo simulations. Through systematic analysis of one-dimensional model systems, we evaluate their computational efficiency and representational power against established methods. Our numerical experiments suggest some efficient training methods and we explore how the computational cost scales with desired precision, particle number, and system parameters. Roughly speaking, KANs seem to be 10 times cheaper computationally than other neural-network-based Ansätz . We also introduce a novel approach for handling strong short-range potentials—a persistent challenge for many numerical techniques—which generalizes efficiently to higher-dimensional, physically relevant systems with short-ranged strong potentials common in atomic and nuclear physics.

1-dimensional systems↗

DOE OSTI2025

Adopting Code Verification Methodology Based on Model Form

Code verification is an essential part of credibility analysis for computational models. It assesses whether the mathematical model is implemented correctly into the code and whether the numerical methods behave consistently, and is done before solution verification and validation. Robust guidance for code verification exists in the literature. However, there is no known, concise guide for selecting the approach based on the model form that also presents an overview of the common elements. This document was written to address this gap as an accessible reference for beginning a code-verification effort.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

Adapting Code Verification Methodology to Model Form

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

Multi-Entity Simulation with CoSim Toolbox

Co-simulation is an analysis technique for linking multiple software models during runtime by facilitating data exchange and simulation time synchronization. There are numerous challenges when constructing an effective co-simulation including simulation tool installation, data management, and writing new models in a manner compatible with the co-simulation framework of choice. CoSim Toolbox is an integration of multiple pieces of software designed to make assembling such a co-simulation in HELICS easier. This report summarizes the existing capabilities of CoSim Toolbox and outlines future development plans.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Proteogenomic characterization of difficult-to-treat breast cancer with tumor cells enriched through laser microdissection

Abstract Background Breast cancer (BC) is the most commonly diagnosed cancer and the leading cause of cancer death among women globally. Despite advances, there is considerable variation in clinical outcomes for patients with non-luminal A tumors, classified as difficult-to-treat breast cancers (DTBC). This study aims to delineate the proteogenomic landscape of DTBC tumors compared to luminal A (LumA) tumors. Methods We retrospectively collected a total of 117 untreated primary breast tumor specimens, focusing on DTBC subtypes. Breast tumors were processed by laser microdissection (LMD) to enrich tumor cells. DNA, RNA, and protein were simultaneously extracted from each tumor preparation, followed by whole genome sequencing, paired-end RNA sequencing, global proteomics and phosphoproteomics. Differential feature analysis, pathway analysis and survival analysis were performed to better understand DTBC and investigate biomarkers. Results We observed distinct variations in gene mutations, structural variations, and chromosomal alterations between DTBC and LumA breast tumors. DTBC tumors predominantly had more mutations inTP53,PLXNB3, Zinc finger genes, and fewer mutations inSDC2,CDH1,PIK3CA,SVIL, andPTEN. Notably, Cytoband 1q21, which contains numerous cell proliferation-related genes, was significantly amplified in the DTBC tumors. LMD successfully minimized stromal components and increased RNA–protein concordance, as evidenced by stromal score comparisons and proteomic analysis. Distinct DTBC and LumA-enriched clusters were observed by proteomic and phosphoproteomic clustering analysis, some with survival differences. Phosphoproteomics identified two distinct phosphoproteomic profiles for high relapse-risk and low relapse-risk basal-like tumors, involving several genes known to be associated with breast cancer oncogenesis and progression, includingKIAA1522,DCK,FOXO3,MYO9B,ARID1A,EPRS,ZC3HAV1, andRBM14. Lastly, an integrated pathway analysis of multi-omics data highlighted a robust enrichment of proliferation pathways in DTBC tumors. Conclusions This study provides an integrated proteogenomic characterization of DTBC vs LumA with tumor cells enriched through laser microdissection. We identified many common features of DTBC tumors and the phosphopeptides that could serve as potential biomarkers for high/low relapse-risk basal-like BC and possibly guide treatment selections.

Oncology↗

DOE OSTI2024

Cancer Incidence Trends in Successive Social Generations in the US

Importance: The incidence of some cancers in the US is increasing in younger age groups, but underlying trends in cancer patterns by birth year remain unclear. Objective: To estimate cancer incidence trends in successive social generations. Design, Setting, and Participants: In this cohort study, incident invasive cancers were ascertained from the Surveillance, Epidemiology, and End Results (SEER) program’s 13-registry database (November 2020 submission, accessed August 14, 2023). Invasive cancers diagnosed at ages 35 to 84 years during 1992 to 2018 within 152 strata were defined by cancer site, sex, and race and ethnicity. Exposure: Invasive cancer. Main Outcome and Measures: Stratum-specific semiparametric age-period-cohort (SAGE) models were fitted and incidence per 100 000 person-years at the reference age of 60 years was calculated for single-year birth cohorts from 1908 through 1983 (fitted cohort patterns [FCPs]). The FCPs and FCP incidence rate ratios (IRRs) were compared by site for Generation X (born between 1965 and 1980) and Baby Boomers (born between 1946 and 1964). Results: A total of 3.8 million individuals with invasive cancer (51.0% male; 8.6% Asian or Pacific Islander, 9.5% Hispanic, 10.4% non-Hispanic Black, and 71.5% non-Hispanic White) were included in the analysis. In Generation X vs Baby Boomers, FCP IRRs among women increased significantly for thyroid (2.76; 95% CI, 2.41-3.15), kidney (1.99; 95% CI, 1.70-2.32), rectal (1.84; 95% CI, 1.52-2.22), corpus uterine (1.75; 95% CI, 1.40-2.18), colon (1.56; 95% CI, 1.27-1.92), and pancreatic (1.39; 95% CI, 1.07-1.80) cancers; non-Hodgkins lymphoma (1.40; 95% CI, 1.08-1.82); and leukemia (1.27; 95% CI, 1.03-1.58). Among men, IRRs increased for thyroid (2.16; 95% CI, 1.87-2.50), kidney (2.14; 95% CI, 1.86-2.46), rectal (1.80; 95% CI, 1.52-2.12), colon (1.60; 95% CI, 1.32-1.94), and prostate (1.25; 95% CI, 1.03-1.52) cancers and leukemia (1.34; 95% CI, 1.08-1.66). Lung (IRR, 0.60; 95% CI, 0.50-0.72) and cervical (IRR, 0.71; 95% CI, 0.57-0.89) cancer incidence decreased among women, and lung (IRR, 0.51; 95% CI, 0.43-0.60), liver (IRR, 0.76; 95% CI, 0.63-0.91), and gallbladder (IRR, 0.85; 95% CI, 0.72-1.00) cancer and non-Hodgkins lymphoma (IRR, 0.75; 95% CI, 0.61-0.93) incidence decreased among men. For all cancers combined, FCPs were higher in Generation X than for Baby Boomers because gaining cancers numerically overtook falling cancers in all groups except Asian or Pacific Islander men. Conclusions and Relevance: In this model-based cohort analysis of incident invasive cancer in the general population, decreases in lung and cervical cancers in Generation X may be offset by gains at other sites. Generation X may be experiencing larger per-capita increases in the incidence of leading cancers than any prior generation born in 1908 through 1964. On current trajectories, cancer incidence could remain high for decades.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2025

Resolution requirements for numerical modeling of neutrino quantum kinetics

Neutrino quantum kinetics is a rapidly evolving field in computational astrophysics, with a primary focus on collective neutrino oscillations in core-collapse supernovae and postmerger phases of binary neutron star mergers. In recent years, there has been considerable debate concerning resolution dependence in numerical simulations. In this paper, we conduct a comprehensive resolution study in both angular- and spatial directions by using two independent schemes of quantum kinetic neutrino transport: finite volume and pseudospectral methods. We complement our discussion by linear stability analysis including inhomogeneous modes. Our result suggests that decreasing spatial resolutions underestimates the growth of flavor instability, and then leads to wrong asymptotic states of flavor conversions, which potentially has a critical impact on astrophysical consequences. We further delve into numerical results of low resolution simulations, that reveals the underlying mechanism responsible for numerical artifacts caused by insufficient resolutions. Finally, this study settles the debate on requirements of resolutions and serves as a guideline for numerical modeling of quantum kinetic neutrino transport.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2026

Experimental characterization and analysis of phase change material-based thermal energy storage system for refrigerated display case

Refrigerated display cases that are used to store and exhibit food products in supermarkets and retail spaces consume a significant portion of the buildings’ total electricity. Importantly, the refrigeration-related energy cost and demand charges are greatly affected by the time-of-use electricity pricing and demand rates, which are at their maximum during peak hours, typically when refrigeration energy consumption is also high. Using energy storage to shift the refrigeration load from peak to off-peak hours can greatly reduce the operational costs in the supermarket. This study demonstrates a phase change material-based thermal energy storage (TES) system, specifically designed in stackable units, that can be integrated with an open vertical refrigerated display case. We perform numerical and experimental characterization that includes finite-difference modeling for the TES, prototype fabrication, and laboratory evaluation, followed by a preliminary system-level analysis to predict the impact of TES on the refrigerated case performance, energy use, and energy cost. The results show that the dedicated latent TES for refrigerated cases can provide a specific energy of 50.4 Wh/kg and a specific power of 15.5 W/kg. The TES can be charged during 12 h of the off-peak period and discharged at various rates during 4 to 6 h of the peak period, thereby shifting the refrigeration load from the peak to the off-peak period. Consequently, annual cost savings up to 19% can be achieved, depending on the thermal load, the summer/winter peak electricity pricing, and the transition temperature of the phase change material used.

25 ENERGY STORAGE↗

DOE OSTI2026

JetGP: A derivative enhanced Gaussian process library

Derivative enhanced Gaussian Processes (DEGPs) can significantly improve surrogate model accuracy over standard Gaussian Process (GP) formulations by incorporating derivative information. However, standard implementations scale poorly with dimension, limiting their use in high dimensional engineering problems. JetGP is a Python framework that unifies existing derivative enhanced GP methodologies into a single library and extends them to support arbitrary order derivative information. The library implements four complementary formulations: standard derivative enhanced Gaussian Processes (DEGP), directional DEGP (DDEGP), generalized directional DEGP (GDDEGP), and weighted DEGP (WDEGP). By unifying these approaches in a consistent interface with robust numerical implementations, JetGP enables practitioners to balance predictive accuracy and computational efficiency for high dimensional optimization, uncertainty quantification, and sensitivity analysis in engineering design.

Derivative enhanced Gaussian process↗

DOE OSTI2025

Intrinsic Thermal Hall Effect in Mott Insulators

Here, in light of recent experimental data indicating a substantial thermal Hall effect in square lattice antiferromagnetic Mott insulators, we investigate whether a simple Mott insulator can sustain a finite thermal Hall effect. We verify that the answer is “no” if one performs calculations within a spin-only low-energy effective spin model with noninteracting magnons. However, by performing determinant quantum Monte Carlo simulations, we show the single-band 𝑡−𝑡′−𝑈 Hubbard model coupled to an orbital magnetic field does support a finite thermal Hall effect when 𝑡′ ≠0 and 𝐵 ≠0 in the Mott insulating phase. We argue that the (carrier agnostic) necessary conditions for observing a finite thermal Hall effect are time-reversal and particle-hole symmetry breaking. By considering magnon-magnon scattering using a semiclassical Boltzmann analysis, we illustrate a physical mechanism by which finite transverse thermal conductivity may arise, consistent with our symmetry argument and numerical results. Our results indicate that square and triangular lattices with SU(2) symmetry can support a finite thermal Hall effect and call for a critical reexamination of thermal Hall effect data in insulating magnets, as the magnon contribution should not be excluded a priori .

Mott insulators↗