SEARCH · Search NASA

Results for “computational geometry”

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 181 records · Page 10

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Performant automatic differentiation of local coupled cluster theories: Response properties and ab initio molecular dynamics

In this work, we introduce a differentiable implementation of the local natural orbital coupled cluster (LNO-CC) method within the automatic differentiation framework of the PySCFAD package. The implementation is comprehensively tuned for enhanced performance, which enables the calculation of first-order static response properties on medium-sized molecular systems using coupled cluster theory with single, double, and perturbative triple excitations [CCSD(T)]. We evaluate the accuracy of our method by benchmarking it against the canonical CCSD(T) reference for nuclear gradients, dipole moments, and geometry optimizations. In addition, we demonstrate the possibility of property calculations for chemically interesting systems through the computation of bond orders and Mössbauer spectroscopy parameters for a [NiFe]-hydrogenase active site model, along with the simulation of infrared spectra via ab initio LNO-CC molecular dynamics for a protonated water hexamer.

Chemistry↗

DOE OSTI2024

Cold Plasma Measurements

We have continued the simulation campaign in support of our ongoing magnetospheric cold plasma research project. This project aims to develop the next-generation particle instruments to measure the properties of the cold particle populations in the Earth’s magnetosphere. For this purpose, simulations have been performed with a Particle-In-Cell (PIC) code called the Curvilinear PIC (CPIC). The code is formulated in curvilinear geometry and couples the standard PIC algorithm with algorithms for the generation and adaptation of the underlaying computational mesh. It conforms to complex objects like spacecraft and it can place more grid points in regions where higher resolution is needed. The code also features a scalable solver based on the multigrid algorithm and it is fully parallelized via domain decomposition and MPI.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

BlueCRAB Domain Overlapping Coupling: Theory & Verification

Coupling low- and high-fidelity codes is a useful way to model complex engineering systems. Thanks to the high-fidelity code, complex phenomena can be resolved in areas of the system where this is required, yet the efficiency of the lower-fidelity code is still retained in modeling the rest of the system. This document details the theory and implementation of coupling two different thermal-hydraulic codes: the system thermal-hydraulics (STH) code System Analysis Module (SAM) and the coarse-mesh computational fluid dynamics (CFD) code Pronghorn. Both applications are included in the Comprehensive Reactor Analysis Bundle (BlueCRAB) code suite and are based on the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework. The domain overlapping (DO) coupling approach was adopted, as it offers proven advantages over more conventional domain decomposition methods. In the DO coupling, SAM provides Pronghorn with boundary conditions that depend on the system-level simulation of the entire plant. In return, the overlapping coupled SAM components–termed “surrogate components”–receive friction factors and source terms computed online based on the Pronghorn simulation. The framework is designed to be generic and enable coupling regardless of geometry and the number of inlet/outlet boundaries in the DO coupled domain. The developed method leads to consistent pressure drops, enthalpies, and scalar concentrations when comparing the coupled SAM and Pronghorn simulations. This document presents the DO coupling approach, along with two verification and two demonstration cases. The proposed problems explore different physical aspects relevant to nuclear reactor analysis, including buoyancy-driven flows, complex flow patterns, and multiple inlets and outlets. Periodically, new versions of this “BlueCRAB Domain Overlapping Coupling: Theory & Verification” report will be issued to reflect future developments and verification tests.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2024

Hadamard products and BPS networks

We study examples of fourth-order Picard-Fuchs operators that are Hadamard products of two second-order Picard-Fuchs operators. Each second-order Picard-Fuchs operator is associated with a family of elliptic curves, and the Hadamard product computes period integrals on the fibred product of the two elliptic surfaces. We construct 3-cycles on this geometry as the union of 2-cycles in the fibre over contours on the base. We then use the special Lagrangian condition to constrain the contours on the base. This leads to a construction that is reminiscent of spectral networks and exponential networks that have previously appeared in string theory literature.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2026

AutoBEM: A scalable framework for nationwide building energy simulation and retrofit evaluation in the United States

This paper presents AutoBEM, an integrated, automated framework for nationwide building energy modeling and retrofit evaluation in the United States. Unlike prior UBEM platforms that either rely primarily on representative stock sampling or operate at city scale, AutoBEM automates the generation of building-resolved, physics-based EnergyPlus/OpenStudio simulation models at national scale using GIS-derived geometry, prototype-based assumptions, and standardized scalable workflows. Leveraging the Model America dataset and high-performance computing, AutoBEM generates and simulates energy models for 122.9 million buildings, representing 97.8% of the U.S. building stock. These models are being made publicly and freely available as the Model America v1.0 (MAv1) dataset. AutoBEM supports detailed, building-level assessments of energy consumption, CO2 emissions, and post-processed anthropogenic heat emissions (AHE), and evaluates 151 energy conservation measures (ECMs) using localized utility pricing and building characteristics. In addition, AutoBEM incorporates both typical and future climate conditions through integration with Typical Meteorological Year (TMY) and Future TMY (fTMY) weather data derived from IPCC scenarios. In a case study of Phoenix, Arizona, AutoBEM identified several high-efficiency HVAC upgrades and selected envelope measures with short modeled payback periods (1.5 years) for certain building types and standards. Simulations under future climate scenarios (SSP5–RCP8.5) project an 11.3% increase in electricity use and a 32% reduction in natural gas demand by 2100, underscoring the need for climate-adaptive retrofit planning. By enabling reproducible, bottom-up, and location-specific analysis at scale, AutoBEM provides a step toward a national digital twin of the built environment and supports data-driven screening and planning for decarbonization, resilience, and energy equity.

Li, Hang [ORNL] (ORCID:0000000306001920)↗

DOE OSTI2025

Multireference Equation-of-Motion Driven Similarity Renormalization Group: Theoretical Foundations and Applications to Ionized States

We present a formulation and implementation of an equation-of-motion (EOM) extension of the multireference driven similarity renormalization group (MR-DSRG) formalism for ionization potentials (IP-EOM-DSRG). The IP-EOM-DSRG formalism results in a Hermitian generalized eigenvalue problem, delivering accurate ionization potentials for strongly correlated systems. The EOM step scales as O(N 5 ) with the basis set size N, allowing for efficient calculation of spectroscopic properties, such as transition energies and intensities. The IP-EOM-DSRG formalism is combined with three truncation schemes of the parent MR-DSRG theory: an iterative nonperturbative method with up to two-body excitations [MR-LDSRG(2)] and second- and third-order perturbative approximations [DSRG-MRPT2/3]. We benchmark these variants by computing (1) the vertical valence ionization potentials of a series of small molecules at both equilibrium and stretched geometries; (2) the spectroscopic constants of several low-lying electronic states of the OH, CN, N 2 + , and CO + radicals; and (3) the binding curves of low-lying electronic states of the CN radical. A comparison with experimental data and theoretical results shows that all three IP-EOM-DSRG methods accurately reproduce the vertical ionization potentials and spectroscopic constants of these systems. Notably, the DSRG-MRPT3 and MR-LDSRG(2) versions outperform several state-of-the-art multireference methods of comparable or higher cost.

Hamiltonians↗

DOE OSTI2024

Measuring the Loschmidt Amplitude for Finite-Energy Properties of the Fermi-Hubbard Model on an Ion-Trap Quantum Computer

Calculating the equilibrium properties of condensed-matter systems is one of the promising applications of near-term quantum computing. Recently, hybrid quantum-classical time-series algorithms have been proposed to efficiently extract these properties from a measurement of the Loschmidt amplitude ⟨ ψ | e − i H ^ t | ψ ⟩ from initial states | ψ ⟩ and a time evolution under the Hamiltonian H ^ up to short times t . In this work, we study the operation of this algorithm on a present-day quantum computer. Specifically, we measure the Loschmidt amplitude for the Fermi-Hubbard model on a 16 -site ladder geometry (32 orbitals) on the Quantinuum H2-1 trapped-ion device. We assess the effect of noise on the Loschmidt amplitude and implement algorithm-specific error-mitigation techniques. By using a thus-motivated error model, we numerically analyze the influence of noise on the full operation of the quantum-classical algorithm by measuring expectation values of local observables at finite energies. Finally, we estimate the resources needed for scaling up the algorithm. Published by the American Physical Society 2024

Physics↗

DOE OSTI2025

Fermilab Automation of Coil Winding System

The automated coil winder project creates a repeatable, scalable method for producing superconducting magnet coils with increased accuracy and efficiency. The design uses an ortho cyclic winding method to achieve high fill factors by combining computer numerical control (CNC), microcontroller-based controls, and sophisticated optimization techniques. The method optimizes coil turns, core geometry, and current consumption by utilizing analytical models based on magnetic moment, Ohm’s law, and power constraints. For particle accelerator applications, prototype testing using a 3D-printed coil winder has shown enhanced winding repeatability, thermal management, and manufacturability.

Dzida, Rafal [Northern Illinois U.]↗

DOE OSTI2024

Directional response of several geometries for reactor-neutrino detectors

Here we have modeled six abstracted detector designs, with the goal of determining their ability to resolve the direction to an antineutrino source, including two for which we have operational data for validating our computer modeling and analytical processes. We have found that the most promising options, regardless of scale and range, have angular resolutions on the order of a few degrees, which is better than any yet achieved in practice by a factor of at least two. We examine and compare several approaches to detector geometry for their ability not only to detect the inverse beta decay (IBD) reaction, but also to determine the source direction of incident antineutrinos. The information from these detectors provides insight into reactor power and burning profile, which is especially useful in constraining the clandestine production of weapons material. In a live deployment, a nonproliferation detector must be able to isolate the subject reactor, possibly from a field of much-larger power reactors; directional sensitivity can help greatly with this task. We also discuss implications for using such detectors in longer-distance observation of reactors, from a few kilometers to hundreds of kilometers.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2024

Cu site differentiation in tetracopper(I) sulfide clusters enables biomimetic N 2 O reduction

Copper clusters feature prominently in both metalloenzymes and synthetic nanoclusters that mediate catalytic redox transformations of gaseous small molecules. Such reactions are critical to biological energy conversion and are expected to be crucial parts of renewable energy economies. However, the precise roles of individual metal atoms within clusters are difficult to elucidate, particularly for cluster systems that are dynamic under operating conditions. Here, we present a metal site-specific analysis of synthetic Cu 4 (μ 4 -S) clusters that mimic the Cu Z active site of the nitrous oxide reductase enzyme. Leveraging the ability to obtain structural snapshots of both inactive and active forms of the synthetic model system, we analyzed both states using resonant X-ray diffraction anomalous fine structure (DAFS), a technique that enables X-ray absorption profiles of individual metal sites within a cluster to be extracted independently. Using DAFS, we found that a change in cluster geometry between the inactive and active states is correlated to Cu site differentiation that is presumably required for efficient activation of N 2 O gas. More precisely, we hypothesize that the Cu δ+ ∙∙∙Cu δ- pairs produced upon site differentiation are poised for N 2 O activation, as supported by computational modeling. These results provide an unprecedented level of detail on the roles of individual metal sites within the synthetic cluster system and how those roles interplay with cluster geometry to impact the reactivity function. We expect this fundamental knowledge to inform understanding of metal clusters in settings ranging from (bio)molecular to nanocluster to extended solid systems involved in energy conversion.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Experimental data for damage mechanics simulation challenge

While there are many computational approaches for simulating damage in rock and other materials, few have been ground truth tested with either known experimental data or with blind data sets. Here, in this work, we present a bench-mark laboratory data set for a damage mechanics challenge to compare computational approaches on damage evolution in brittle-ductile materials. The samples were fabricated through additive manufacturing to produce repeatable specimens designed to fail in controlled ways. The failure was induced in the samples using a 3-point bending test to produce different Modes such as Mode I and mixed Modes including I-II, I-III and I-II-III Modes to generate a calibration data set and a blind challenge data set. Data collected included spatial and temporal measurements from traditional digital load–displacement sensors, 2D digital image correlation measurement to map surface deformations, 3D X-ray microscopy to ground-truth the crack-failure geometry, and laser profilometry to capture surface roughness. The data sets are available, on a data repository, to the community to advance computational models to improve our ability to predict damage in brittle-ductile materials.

3-point bending↗

DOE OSTI2025

GPU-friendly surface model for Monte-Carlo detector simulations

The demands for Monte-Carlo simulation are drastically increasing with the Large Hadron Collider’s high-luminosity upgrade, and are expected to exceed the currently available compute resources. At the same time, modern high-performance computing has adopted powerful hardware accelerators, particularly GPUs. The AdePT and Celeritas projects aim to address the demanding computational needs by leveraging these heterogeneous computing architectures. While both have successfully ported realistic detector simulations to GPUs using the VecGeom library, the complexity of geometry modeling emerged as a bottleneck. Thread divergence and high register usage were degrading the GPU performance. Therefore, a new, GPU-friendly surface-based model has been introduced in the VecGeom library that decomposes the divergent code of the 3D primitive solids into simpler and more balanced surface algorithms. In this work, we present the latest developments, focusing on the additions required to efficiently model complex setups like the CMS Phase-2 geometry. This includes memory reduction techniques, and adding accelerating structures for faster traversal.

Diederichs, Severin [CERN]↗

DOE OSTI2025

Design and Analysis of an Integrated Additively Manufactured Test Article for Plasma-Facing Components

Current plasma-facing components (PFCs) used in helium-cooled divertor modules are complex structures with tungsten tile, steel sleeve components, and cartridges, all assembled in a helium-cooled multiple jet (HEMJ) structure. The goal of this project is to simplify the complex PFC design using additive manufacturing techniques to create a single integrated tungsten test article. Apart from the flexibility this opens up in exploring a wide array of geometries for the article, having a single integrated article significantly reduces the number of joints and parts in the article, thus reducing chances of leaks. A process called electron beam melting has shown to produce very high-density samples and unique geometries, enabling HEMJ or similar designs. To validate and optimize this novel design, the model underwent a series of computational fluid dynamics and finite element analysis simulations to replicate steady-state heat flux in the divertors. The simulations presented in this study consider a steady-state base heat flux of 5 MW/m 2 , with water serving as the coolant. Future research will explore the use of helium as a coolant, simulate edge-localized-mode conditions, and include experimental validation. Since 3D-printed tungsten is anisotropic, the build direction versus build plane of the article are taken into consideration for the test article strength. Because of the high operating temperatures and low ductility of tungsten, thermal creep and brittle fracture are important failure mechanisms to consider. In conclusion, the cap is evaluated with various flow velocities and nozzle diameters, and an optimal design choice is made for which this cap will survive the divertor conditions with a conservative safety margin.

Divertor↗

DOE OSTI2025

A 3D helical filament surrogate model for 3D tokamak equilibria

A novel approach for efficient representation of three-dimensional (3D) tokamak equilibria is investigated, where a set of helical current filaments occupying the plasma region are employed to resolve deviations from the two-dimensional (2D) axi-symmetric state. A discrete set of 3D filaments, located at rational surfaces for a given toroidal mode number n and following the 2D equilibrium field lines (thus forming closed current loops), are found to provide a surrogate model of 3D equilibria with reasonable accuracy. Specifically, application of the filament model to 3D perturbed equilibria, due to the resonant magnetic perturbation (RMP) in DIII-D and MAST-U discharges, reveals that (1) a single helical filament per rational surface is sufficient; (2) 21 such helical filaments are capable of representing the n = 2 3D response field in MAST-U with less than 10% relative error as compared to that computed by a full magnetohydrodynamic code; (3) optimizing currents (both amplitude and phase) flowing in 3D filaments with fixed geometry, the highest accuracy fitting is found to depend on the characteristics of the 3D equilibria such as the coil current phasing of the RMP coils in our case studies. Here, whis filament approach is also applicable for generating surrogate models of other type of 3D tokamak equilibria, including those during the initial phase of the plasma disruption.

MARS-F↗

DOE OSTI2025

Electromagnetic Modeling of Lattice Structures in Additively Manufactured Electric Machines

Recent advancements in metal additive manufacturing (AM) show great potential to revolutionize the design and manufacturing of electromagnetic components used in the field of electrical engineering. Lattice structures directly printed by AM processes typically offer better structural performance with reduced weight, such as high stiffness, surface area, elongation, energy absorption, and porosity, than the solid counterpart. This paper aims to study the electromagnetic modeling of lattice structures used in additively manufactured magnetic cores or windings for electric machines. Three dimensional (3D) electromagnetic finite element (FE) analysis with high performance computing (HPC) shows the highest fidelity in predicting the electromagnetic performance of designs with lattice structures by preserving complex geometry details. FE-based homogenization methods have also been explored to potentially speed up concept design. A case study based on an additively manufactured axial-flux permanent magnet machine with a Hilbert pattern stator validates the discussed electromagnetic modeling approaches.

Chung, Shi-Uk [ANSYS, Inc.]↗

DOE OSTI2024

Interlaced Characterization and Calibration (ICC) for Improved Computational Simulation Credibility

Accurate material characterization and model calibration are pivotal for simulations used for high-consequence engineering decisions. Current characterization and calibration methods (1) use simplified test specimen geometries and global data, (2) cannot guarantee that sufficient characterization data is collected for a specific model of interest, (3) provide only mean parameter values with no uncertainty quantification, and (4) are sequential, inflexible, and time-consuming. This work developed a new paradigm—coined Interlaced Characterization and Calibration (ICC)—which drives forward the state-of-the-art in model calibration by bringing together recent advancements into one improved workflow. The ICC paradigm (1) employs tools to efficiently use full-field data to calibrate high-fidelity material models, (2) aligns the data needed with the data collected by adopting an optimal experimental design protocol, (3) provides uncertainty metrics on the calibrated model parameters, and (4) incorporates these advances into a quasi real-time feedback loop. The ICC framework was validated synthetically with both low-fidelity and high-fidelity simulations paired with several different elastoplastic material models, and was also demonstrated experimentally with an aluminum 6061 cruciform exemplar specimen. Results showed that the ICC framework—in which Bayesian optimal experimental design actively guided the experiment— resulted in calibrations with similar or better accuracy than predetermined experiments based on subject matter expertise. Moreover, the ICC framework produced a complete model calibration— with quantified uncertainties on model parameters—in 1 week, a 5 - 10× increase in efficiency over traditional approaches. Thus, the ICC paradigm improves both the calibration process and quality, by (1) improving efficiency, which increases agility of solid mechanics modeling and enables utilization of computational simulation (CompSim) at earlier stages of the design cycle and (2) providing quantified, and in some cases reduced, parameter uncertainties, which increases confidence in model predictions and supports credible decision making.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

CFD modeling of turbulent air flow in self-heated gyroid TPMS structures: Thermal-hydraulic performance and validation

The application of mathematically derived geometries, such as triply periodic minimal surface (TPMS) lattices, has garnered significant interest across various fields, including the nuclear sector, due to their superior thermal-hydraulic characteristics for heat transfer compared to traditional plain or finned tubes. Here, this study validates a computational fluid dynamics (CFD) model, evaluates different turbulence models and CFD model settings, and performs uncertainty quantification to provide a comprehensive analysis. Despite extensive research on CFD modeling of TPMS lattices, such as gyroid and diamond geometries, there is a notable lack of publicly available literature providing comprehensive details on numerical analysis aspects, including convergence and methodological best practices. This study embarks on a benchmark analysis of a gyroid geometry to evaluate its thermal-hydraulic performance under turbulent flow conditions and scrutinize various CFD model configurations. The main contributions of this work include validating the CFD model, assessing and comparing different turbulence models, and enhancing pressure drop and temperature prediction capabilities. The results aim to support the development of methodologies needed to benchmark and enhance numerical analysis techniques for TPMS lattices. This work seeks to complement the existing body of knowledge, support the development of TPMS reactor concepts, and improve best practices for CFD modeling of TPMS lattices, ultimately advancing methodologies to support future applications in this domain.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2025

Additional considerations in analytical solution for time-dependent heat conduction in a three-dimensional multilayer sphere

This work presents an analytical method to solve the heat conduction equation in three dimensions for problems consisting of multilayer concentric spheres. The method can be used to treat time-varying heat conduction problems where the heat source that drives the transient is time-invariant. Equally applicable to all Poisson-type problems with concentric spherical geometry, the method consists of representing the solution as a summation of weighted eigenfunctions. The weights for each eigenfunction are computed algebraically. Previous work has already established the core constituents of the methodology. The current work augments the existing methods by including consideration of nonzero interface resistance between layers and explicit discussion on the boundary condition homogenization required to treat inhomogeneous problems. Also, two demonstration problems are presented. One demonstration problem is based on the method of manufactured solutions and therefore allows for comparison with exact expressions for the solution temperature distribution. The second, more complex, demonstration problem relies on the finite element method for comparisons. The expected convergence behavior is observed for both demonstration problems.

97 - MATHEMATICS AND COMPUTING↗