SEARCH · Engineering Papers

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 1,171 records · Page 65

All sources NASA NTRS DOE OSTI

DOE OSTI2022

LivermorE Al Projector for Computed Tomography Tasks

With recent computed tomography (CT) efforts using Artificial Intelligence (AI) and Deep Learning (DL)techniques, there is a strong need for differentiable forward projection models that can be integrated into existing DL frameworks. We developed a pytorch-based package library providing differentiable forward and back projection functions and classes to facilitate forward and back propagation of CT operations in the training procedure. This forward projectors support three CT projection geometries: cone, parallel and modular beams. This package can be used with both CPU and GPU with CUDA.

Kim, Hyojin↗

DOE OSTI2025

Effective Polarizability in Near-Field Microscopy of Phonon-Polariton Resonances

We investigate the resonant characteristics of planar surfaces and distinct edges of structures with the excitation of phonon-polaritons. We analyze two materials supporting phonon-polariton excitations in the mid-infrared spectrum: silicon carbide, characterized by an almost isotropic dielectric constant, and hexagonal boron nitride, notable for its pronounced anisotropy in a spectral region exhibiting hyperbolic dispersion. We formulate a theoretical framework that accurately captures the excitations of the structure involving phonon-polaritons, predicts the response in scattering-type near-field optical microscopy, and is effective for complex resonant geometries where the locations of hot spots are uncertain. We account for the tapping motion of the probe, perform analysis for different heights of the probe, and demodulate the signal using a fast Fourier transform. Using this Fourier demodulation analysis, we show that light enhancement across the entire apex is the most accurate characteristic for describing the response of all resonant excitations and hot spots. We demonstrate that computing the demodulation orders of light enhancement in the microscope probe accurately predicts its imaging.

Chemistry↗

DOE OSTI2021

Simulation and Modeling of Time-Resolved X-Ray Detector for the Saturn Accelerator

In this work, we present the technology-aided computer design (TCAD) device simulation and modeling of a silicon p-i-n diode for detecting time-dependent X-ray radiation. We show that the simulated forward and reverse breakdown current–voltage characteristics agree well with the measured data under nonradiation environment by only calibrating carrier lifetimes for the forward bias case and avalanche model critical fields for the reverse bias condition. Using the calibrated parameters and other nominal material properties, we simulated the radiation responses of the p-i-n diode and compared with experimental data when the diode was exposed to X-ray radiation at Sandia’s Saturn facility and the Idaho State University (ISU) TriMeV facility. For Saturn’s Gaussian dose-rate pulses, we show three findings from TCAD simulations. First, the simulated photocurrents are in excellent agreement with the measured data for two dose-rate pulses with peak values of 1.16×10 10 and 1.88×10 10 rad(Si)/s. Second, the simulation results of high dose-rate pulses predict increased delayed photocurrents with longer time tails in the diode electrical responses due to excess carrier generation. Third, simulated peak values of diode radiation responses versus peak dose rates at different bias conditions provide useful guidance to determine the dose-rate range that the p-i-n diode can reliably detect in experiment. For TriMeV’s non-Gaussian dose-rate pulse, our simulated diode response is in decent agreement with the measured data without further calibration. We also studied the effects of device geometry, recombination process, and dose-rate enhancement via TCAD simulations to understand the higher measured response in the time after the peak dose-rate radiation for the p-i-n diode exposed to TriMeV irradiation.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2025

Generalized entropy of gravitational fluctuations

The corrections to holographic entanglement entropy from bulk quantum fields in a classical gravitational background are now well understood. They lead, in particular, to unitary Page curves for evaporating black holes. However, the correct treatment of quantum fluctuations of the metric, including graviton excitations, is a longstanding problem. We provide a gauge-invariant prescription for the generalized entropy of gravitons in anti-de Sitter space in terms of areas and bulk entanglement entropy, generalizing the quantum extremal surface prescription to accommodate fluctuations in the semiclassical spacetime geometry. This task requires a careful treatment of the area operator on the graviton Hilbert space and the definition of a “quantum extremal gauge” in which the extremal surface is unperturbed. It also requires us to determine the correct vacuum modular Hamiltonian for the graviton field, which we fix by requiring that it doesn’t contain a boundary term in extremal gauge. We check our prescription with an explicit computation of the vacuum-subtracted generalized entropy of states containing a graviton in an AdS-Rindler background. Our results exactly match vacuum-subtracted von Neumann entropies for stress-tensor excited states in holographic conformal field theory with d > 2 dimensions. We also use covariant phase space techniques to give a partial proof of our prescription when the entanglement wedge for the background spacetime has a bifurcate Killing horizon. Along the way, we identify a class of perturbative graviton states that have parametrically larger generalized entropy, in the small G N expansion, than any low-energy excitations of an ordinary quantum field.

1/N expansion↗

DOE OSTI2020

An experimental investigation of the effectiveness of Ar-CO 2 shielding gas mixture for the wire arc additive process

Wire arc additive manufacturing (AM) is the process by which a large, metallic structure is built layer-by-layer using a welding arc to melt a wire feedstock. A novel opportunity exists to alter the shielding gas composition to fabricate distinct geometrical features without altering the other AM parameters. In this paper, shielding gases with varied concentrations of Argon (Ar) and CO 2 was used to deposit three distinct geometric shapes (walls, infill, and overhang) using a wire-based additive manufacturing system utilizing the gas metal arc welding (GMAW) surface tension transfer (STT) process. Computer-aided design (CAD) models were sliced with a custom-built slicer, and the sliced algorithm was converted into optimal robotic toolpaths. A custom virtual instrument (VI) was built in LabVIEW to compare the temperature profiles on the surface during each deposition process. After each deposition, the geometric features were scanned, and the surface waviness value was evaluated. Tensile and Charpy impact coupons were extracted from the wall geometries in the longitudinal and transverse directions and tested. The results indicated that a higher CO 2 content produced higher melt pool temperatures to an extent, while lower contents of CO 2 resulted in a dimensionally accurate geometry. The data also indicated that the 2%/98% CO 2 /Ar blend produced scatter in tensile strength and the analysis of variance (ANOVA) shows significant difference. However, the intermediate range of CO 2 (5–10%) resulted in uniform tensile properties. Altogether, these results indicate that a 5%/95% CO 2 /Ar blend is the ideal shielding gas for lowering process temperatures and improving mechanical properties in wire arc additive manufacturing using the gas metal arc welding surface tension transfer process. Additionally, varying concentrations of Ar/CO 2 can be used within the same part in order to modify the local properties or process parameters such as strength, toughness, temperature, or dimensional features. This may improve overall manufacturing quality without sacrificing specific properties.

42 ENGINEERING↗

DOE OSTI2025

Photon counting with intensified charge coupled device (ICCD) – I. In-depth detector characterization

While the adsorption properties of transition metal catalysts have been widely studied, leading to the discovery of various scaling relations, descriptors of catalytic activity, and well-established computational models, a similar understanding of semiconductor catalysts has not yet been achieved. In this work, we present a high-throughput density functional theory investigation into the adsorption properties of 5 oxides of interest to the photocatalytic CO 2 reduction reaction: TiO 2 (rutile and anatase), SrTiO 3 , NaTaO 3 , and CeO 2 . Using a systematic approach, we exhaustively identify unique surfaces and construct adsorption structures to undergo geometry optimizations. We then perform a data-driven analysis, which reveals the presence of weak adsorption energy scaling relations, the propensity of adsorbates of interest to interact with oxygen surface sites, and the importance of slab deformation upon adsorption. Our findings are presented in the context of experimental observations and in comparison to previously studied classes of catalysts, such as pure metals and tellurium-containing semiconductors, and reinforce the need for a comprehensive approach to the study of site-specific surface phenomena on semiconductors.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Systematic computational study of oxide adsorption properties for applications in photocatalytic CO 2 reduction

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Vulnerable window of yield strength for swelling-driven fracture of phase-transforming battery materials

Abstract Despite numerous experimental and theoretical investigations of the mechanical behavior of high-capacity Si and Ge Li-ion battery anodes, our basic understanding of swelling-driven fracture in these materials remains limited. Existing theoretical studies have provided insights into elasto-plastic deformations caused by large volume change phase transformations, but have not modeled fracture explicitly beyond Griffith’s criterion. Here, we use a multi-physics phase-field approach to model self-consistently anisotropic phase transformation, elasto-plastic deformation, and crack initiation and propagation during lithiation of Si nanopillars. Our computational results reveal that fracture occurs within a “vulnerable window” inside the two-dimensional parameter space of yield strength and fracture energy and highlight the importance of taking into account the surface localization of plastic deformation to accurately predict the magnitude of tensile stresses at the onset of fracture. They further demonstrate how the increased robustness of hollow nanopillars can be understood as a direct effect of anode geometry on the size of this vulnerable window. Those insights provide an improved theoretical basis for designing next-generation mechanically stable phase-transforming battery materials undergoing large volume changes.

25 ENERGY STORAGE↗

DOE OSTI2022

Benchmark evaluation of one dimensional array of HEU moderated and reflected by Lucite

The evaluation of two HEU-Lucite experiments is presented. The critical experiments were performed using the Planet Universal Critical Assembly Machine at Los Alamos National Laboratory in 2019. In these experiments, HEU foils were interleaved with Lucite in a column stack and were moderated and reflected from the square Lucite plates. The neutron data produced were used to approximate the multiplication of each measured configuration. The evaluation of the experiments using calculational techniques is presented in the paper. A Monte-Carlo computational model to be used with MCNP6 was developed for the determination of uncertainties and the effects of various uncertainties on the multiplication factor were studied. The evaluation of the benchmark was performed for comparison of calculation results with the experimental data and validation of calculation techniques in particular the newly available thermal scattering data for the Lucite material. The sensitivity studies include parametric variations of material composition and geometry. Experimental uncertainties were low, with the overall height of the HEU/Lucite stack being the largest contributor to the overall uncertainty in k{sub eff}. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2020

BISON Robustness and Performance Improvements

BISON is a modern finite-element based nuclear fuel performance code that has been under development at the Idaho National Laboratory (USA) since 2009 [1]. The code is applicable to both steady and transient fuel behavior and can be used to analyze 1D (spherically symmetric), 2D (axisymmetric and generalized plane strain) or 3D geometries. BISON is the fuel performance code used within CASL for LWR fuel under both normal operating and accident conditions. BISON is built using the INL Multiphysics ObjectOriented Simulation Environment, or MOOSE [2, 3]. MOOSE is a massively parallel, finite element-based framework to solve systems of coupled non-linear partial differential equations using the Jacobian-Free Newton Krylov (JFNK) method [4]. This enables investigation of computationally large problems, for example a full stack of discrete pellets in a LWR fuel rod, or every rod in a full reactor core. MOOSE supports the use of complex two and three-dimensional meshes and uses implicit time integration, important for the widely varied time scale in nuclear fuel simulation. An object-oriented architecture is employed which greatly minimizes the programming effort required to add new material and behavioral models. The flexibility of the implicit and fully coupled multiphysics approach comes with a need for constructing suitable approximations for the Jacobian matrix of the coupled system used for either preconditioning a Krylov solve or in a direct Newton solve. Preconditioning options for Bison problems need to be revisited with new preconditioning methods becoming available.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2025

A geometric framework for momentum-based optimizers for low-rank training

Low-rank pre-training and fine-tuning have recently emerged as promising techniques for reducing the computational and storage costs of large neural networks. Training low-rank parameterizations typically relies on conventional optimizers such as heavy ball momentum methods or Adam. In this work, we identify and analyze potential difficulties that these training methods encounter when used to train low-rank parameterizations of weights. In particular, we show that classical momentum methods can struggle to converge to a local optimum due to the geometry of the underlying optimization landscape. To address this, we introduce novel training strategies derived from dynamical low-rank approximation, which explicitly account for the underlying geometric structure. Our approach leverages and combines tools from dynamical low-rank approximation and momentum-based optimization to design optimizers that respect the intrinsic geometry of the parameter space. We validate our methods through numerical experiments, demonstrating faster convergence, and stronger validation metrics at given parameter budgets.

Schotthoefer, Steffen [ORNL] (ORCID:00000002156965↗

DOE OSTI2023

Butterfly, Vinylidene-Like, Monobridged and Trans Structures of Si 2 H 2 + : Comparison to the Well-Characterized Neutral Si 2 H 2

This paper investigates four of the constitutional isomers of Si 2 H 2 + , namely the butterfly, vinylidene-like, monobridged, and trans structures. These isomer geometries were all studied using the CCSD(T) method with basis sets as large as cc-pV5Z. Higher level methods CCSDT and CCSDT(Q) were used for final energetics. It is found that the butterfly isomer has the lowest energy, followed by vinylidene-like and monobridged structures; the trans structure has the highest energy of the four. All structures were compared with their neutral counterparts. Partial charges are computed and it is found that all isomers share the positive charge between the silicon atoms. NBO analysis shows that the cation Si-Si bond order is reduced by 0.5 relative to the neutral versions of the same isomer in the case of the butterfly, vinylidene-like, and monobridged; and reduced by 0.6 for the trans isomer. Vibrational frequencies, infrared intensities, and dipole moments are predicted to encourage the spectroscopic identification of Si 2 H 2 + .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Reduced-order CFD modeling of cryogenic hydrogen isotope extrusion for pellet fueling

This study presents a reduced-order model (ROM) for computational fluid dynamics (CFD) simulations of cryogenic hydrogen isotope extrusions, focusing on protium (H₂) and deuterium (D₂) piston extruders. Using a 2D axisymmetric ROM in ANSYS-Polyflow, significant computational savings were achieved (runtime reduced from 9∼24 h to 3∼5 min), with extrusion force discrepancies between the 2D ROM and 3D models being on the order of 1%. Parametric studies identified optimal cutoff shear rates in the viscosity model (0.01/s for H₂ and 0.001/s for D₂), providing recommendations for future simulations. Finally, a comprehensive comparison of ROM results with experimental data was performed across varying geometries, cryogenic materials, temperatures, extrusion lengths, and piston velocities. Predictions at low extrusion temperatures met the objective of providing quick and efficient solutions with an acceptable extrusion force error of approximately 10% or less, validating the effectiveness of the 2D ROM approach. However, at high temperatures closer to the triple point, extrusion force error grows, which necessitates developing an improved model that accounts for temperature effects, e.g. melting. Nevertheless, the findings still represent a significant improvement in efficiency of CFD modeling of cryogenic hydrogenic extrusion. The ROM framework can also be extended to tritium (T2) and screw extruders, which will ultimately provide a fast and effective tool for optimizing pellet injector design for ITER and future reactor systems.

Fan, Joy [ORNL] (ORCID:0000000229751735)↗

DOE OSTI2020

Spin Splitting Energy of Transition Metals: A New, More Affordable Wave Function Benchmark Method and Its Use to Test Density Functional Theory

Accurately predicting the spin splitting energy of chemical species is important for understanding their reactivity and magnetic properties, but it is very challenging, especially for molecules containing transition metals. One impediment to progress is the scarcity of accurate benchmark data. Here we report a set of calculations designed to yield reliable benchmarks for simple transition-metal complexes that can be used to test density functional methods that are affordable for large systems of more practical interest. Various wave function methods are tested against experiment for Fe 2+ , Fe 3+ , and Co 3+ , including CASSCF, CASPT2, CASPT3, MRCISD, MRCISD+Q, ACPF, AQCC, CCSD(T), and CASPT2/CCSD(T) and also a new method called CASPT2.5, which is performed by taking the average of the CASPT2 and CASPT3 energies. Furthermore, we find that MRCISD+Q, ACPF, and AQCC require smaller active spaces for good accuracy than are required by CASPT2 and CASPT3, and this aspect may be important for calculations on larger molecules; here we find that CASPT2.5 extrapolated to a complete basis set is the most suitable method—in terms of computational cost and in terms of accuracy on monatomic systems—and therefore we chose this method for molecular benchmarks. Then Kohn–Sham density functional calculations with 60 exchange-correlation functionals are tested for FeF 2 , FeCl 2 , and CoF 2 . We find that MN15-L, M06-SX, and revM06 have very good agreement with CASPT2.5 benchmarks in terms of both the spin splitting energy and the optimized geometry for each spin state. In addition, we recommend def2-TZVP as the most suitable basis set to perform density functional calculations for molecular spin splitting energies; extra polarization functions in the basis set do not help to increase the accuracy of the spin splitting energy in KS calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Experimental and computational studies of the kinetics of the reaction of hydrogen peroxide with the amidogen radical

In this work, the pulsed-laser photolysis/laser-induced fluorescence method is used to study the kinetics of the reaction of NH 2 with H 2 O 2 to yield a second-order rate constant of (2.42 ± 0.55) × 10 –14 cm 3 molecule –1 s –1 at 412 K in 10–22 mbar of Ar bath gas. There are no prior measurements for comparison. To check this value and enable reliable extrapolation to other temperatures, we also compute thermal rate constants for this process over the temperature range 298–3000 K via multi-structural canonical variational transition-state theory with small-curvature multidimensional tunneling (MS-CVT/SCT). The CVT/SCT rate constants are derived using a dual-level direct dynamics approach utilizing single-point CCSD(T)-F12b/cc-pVQZ-F12 energies—corrected for core-valence and scalar relativistic effects—and M06-2X/MG3S geometries, gradients, and Hessians—for all stationary and non-stationary points along the reaction path. The multistructural method with torsional anharmonicity, based on a coupled torsional potential, is then employed to calculate correction factors for the rate constants, accounting for the comprehensive effects of torsional anharmonicity on the kinetics of this reaction system. The final MS-CVT/SCT rate constants are found to be in good agreement with our measurements and can be expressed in modified Arrhenius form as 2.13 × 10 –15 (T/298 K) 4.02 exp(–513 K/T) cm 3 molecule –1 s –1 over the temperature range of 298–3000 K.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Diagonal isometric form for tensor network states in two dimensions

Isometric tensor network states (isoTNS) generalize the isometric form of the one-dimensional matrix product states to tensor networks in two and higher dimensions. Here, we introduce an alternative isometric form for isoTNS by incorporating auxiliary tensors to represent the orthogonality hypersurface. We implement the time evolving block decimation algorithm on this new isometric form and benchmark the method by computing ground states and the real time evolution of the transverse field Ising model in two dimensions on large square lattices of up to 1250 sites. Our results demonstrate that isoTNS can efficiently capture the entanglement structure of two-dimensional area law states. The short-time dynamics is also accurately reproduced even at the critical point. Our isoTNS formulation further allows for a natural extension to different lattice geometries, such as the honeycomb or kagome lattice.

2-dimensional systems↗

DOE OSTI2020

Localization of fluorescent targets in deep tissue with expanded beam illumination for studies of cancer and the brain

Imaging fluorescence through millimeters or centimeters of tissue has important in vivo applications, such as guiding surgery and studying the brain. Often, the important information is the location of one of more optical reporters, rather than the specifics of the local geometry, motivating the need for a localization method that provides this information. We present an optimization approach based on a diffusion model for the fast localization of fluorescent inhomogeneities in deep tissue with expanded beam illumination that simplifies the experiment and the reconstruction. We show that the position of a fluorescent inhomogeneity can be estimated while assuming homogeneous tissue parameters and without having to model the excitation profile, reducing the computational burden and improving the utility of the method. We perform two experiments as a demonstration. First, a tumor in a mouse is localized using a near infrared folate-targeted fluorescent agent (OTL38). This result shows that localization can quickly provide tumor depth information, which could reduce damage to healthy tissue during fluorescence-guided surgery. Second, another near infrared fluorescent agent (ATTO647N) is injected into the brain of a rat, and localized through the intact skull and surface tissue. This result will enable studies of protein aggregation and neuron signaling.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2023

Topological aspects of matters and Langlands program

The Langlands program is a vast mathematical projection linking number theory and geometry. In high-energy physics, a connection with mirror symmetry has been suggested in string theory, but it has been little studied in low-energy physics. In the framework of the Langlands program, we present a unified description of the integer and fractional quantum Hall effect and the duality found in the fractal nature of the energy spectrum of two-dimensional Bloch electrons, statistical physics, and quantum computation. The new unified view of existing dualism presented in this paper raises the entirely new question of how each theory of physics is connected as a piece of the Langlands program.

Physics↗