Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “geometric modelling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 163 records · Page 9

Finite elements for Matérn-type random fields: Uncertainty in computational mechanics and design optimization

This work highlights an approach for incorporating realistic uncertainties into scientific computing workflows based on finite elements, focusing on prevalent applications in computational mechanics and design optimization. We leverage Matérn-type Gaussian random fields (GRFs) generated using the SPDE method to model aleatoric uncertainties, including environmental influences, variating material properties, and geometric ambiguities. Our focus lies on delivering practical GRF realizations that accurately capture imperfections and variations and understanding how they impact the predictions of computational models as well as the shape and topology of optimized designs. Here we describe a numerical algorithm based on solving a generalized SPDE to sample GRFs on arbitrary meshed domains. The algorithm leverages established techniques and integrates seamlessly with the open-source finite element library MFEM and associated scientific computing workflows, like those found in industrial and national laboratory settings. Our solver scales efficiently for large-scale problems and supports various domain types, including surfaces and embedded manifolds. We showcase its versatility through biomechanics and topology optimization applications, emphasizing the potential to influence these domains. The flexibility and efficiency of SPDE-based GRF generation empowers us to run large-scale optimization problems on 2D and 3D domains, including finding optimized designs on embedded surfaces, and to generate design features and topologies beyond the reach of conventional techniques. Moreover, these capabilities allow us to model and quantify geometric uncertainties on reconstructed submanifolds, such as the interpolated surfaces of cerebral aneurysms provided by postprocessing CT scans. In addition to offering benefits in these specific domains, the proposed techniques transcend specific applications and generalize to arbitrary forward and backward problems in uncertainty quantification involving finite elements.

97 MATHEMATICS AND COMPUTING↗

Analysis of Postulated Accident Scenarios for the Transformational Challenge Reactor

The Transformational Challenge Reactor (TCR) is a helium-cooled, yttrium-hydride-moderated reactor that was designed for the U.S. Department of Energy Office of Nuclear Energy. A key objective of the TCR was to employ advanced manufacturing techniques in a nuclear system and demonstrate their potential for revolutionizing the nuclear reactor design process. One purpose of the present work is to demonstrate the safety of the TCR under postulated accidents. Based on RELAP5-3D and COMSOL analyses, the TCR remained below all current safety limits and far below the expected failure limits for the core materials. Another purpose of this work is to provide useful insights and recommendations regarding the application of RELAP5-3D to gas-cooled or other advanced reactors. A novel approach was implemented for simultaneously modeling conduction and radiation in RELAP5-3D, which was found to provide reasonable predictions of radial core, vessel, and ex-vessel heat transfer during postulated events. A multicode approach was also applied, in which high-fidelity COMSOL calculations were used to tune the radial heat transfer parameters in RELAP5-3D. The tuned RELAP5-3D model demonstrated comparable peak temperature predictions as COMSOL, despite a coarse treatment of the core in RELAP5-3D consisting of only two lumped heat structures. This high-fidelity tuning approach enabled enhanced accuracy as well as minimal complexity within the RELAP5-3D model, even for complex fuel geometric designs as in the TCR. Finally, investigations were made into the potential for flow reversal during a pressurized loss-of-forced-flow event in the TCR. The TCR is designed with downward helium flow through the core during normal operation. The RELAP5-3D model predicted that this downward flow would persist, without flow reversal, up to several days after the circulator trip. This was attributed to natural circulation hysteresis effects as have been noted in similar thermofluidic systems. Although flow stagnation and eventual reversal did not lead to unsafe TCR conditions, interesting spatial effects were observed which may have safety relevance for other reactor system designs and coolant types that are designed for downward core flow during normal operation, warranting closer investigation of the flow reversal phenomenon.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Static characterization of a highly optimized streak tube design incorporating a steering slot anode and an aberration-corrected Einzel lens

The BHx streak tube, under development at the Laboratory for Laser Energetics, incorporates a series of novel electron-optics elements to enable high fidelity measurements for high-energy-density physics experiments. The system is engineered to support a 25-mm active photocathode region compatible with sub-picosecond temporal resolution and with 70% internal photoelectron throughput. It exhibits negligible geometric distortion on a flat output screen, making it well suited for variety of measurements, such as ultrafast x-ray spectroscopy. Here, this paper presents characterization data from a prototype unit tested with an ultraviolet laser in a static deflection (non-swept) mode and shows good agreement with the predictions from numerical modeling, including focusing performance and geometric distortion measurements. Key design elements have been demonstrated and de-risked, laying the foundation for dynamic deflection testing. The anticipated improvements in data fidelity are expected to impact the fields as diverse as inertial confinement fusion, laboratory astrophysics, and materials science.

High energy density physics↗

Bubbling Water–Treating DBD Plasma Device Optimization Using Experimental and Computational Methods

A dry air atmospheric pressure volume dielectric barrier discharge is employed to fix nitrogen in water. Producing nitrate for use as nitrogen fertilizer is the primary motivation. A 0D chemistry model is developed and informed by the electrical, and geometric characteristics of the device and the plasma gas temperature. Modeled ozone and nitrate densities are compared to those measured experimentally in the plasma effluent and treated liquid for a range of gas temperatures. Modeled and measured ozone densities are in good agreement; however, the model lacks the liquid chemistry to properly represent the measured nitrate density. A gas temperature-based shift from ozone to NO x producing regimes is observed in both experiment and model, and the reactions responsible are evaluated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A sensitivity analysis to predict the neutronics behavior of samples irradiated in the VTR rabbit system

We report a low-order neutronics model is developed to carry out hundreds of simulations efficiently and investigate the neutronics behavior of samples being irradiated in a test reactor setting under different geometrical constraints. The low-order model allowed for simulations that yield the expected neutronics behavior of any irradiated sample in any environment and allows for the calculation of highly accurate spatially averaged statistics and idealized spatial distributions in the neutron flux. Several benchmarks are performed to evaluate the performance and limitations of the low-order model revealing many important findings. The low-order model predicted the LHGR in the EBR-II driver fuel to within 2.34% by only simulating the fuel rod by itself, which served as a validation for the model. Sensitivity studies investigated 3% enriched UO 2 and U-10Zr being irradiated in the Versatile Test Reactor rabbit system. The analyses investigated a range of combinations of 15 radii and 5 heights for each sample in the rabbit system. Similar data sets are also provided for irradiations in the Advanced Test Reactor’s B-10 irradiation position, which is a thermal neutron spectrum environment. Generalized fits and fit coefficients are obtained for sample heating, reaction rate densities, and local multiplication rate characteristics, allowing the predictions of the neutronics behavior of the samples based on their geometrical constraints. The analyses and fits laid the groundwork for developing a user-end Multiphysics analysis framework to assist and accelerate irradiation experiment design and optimization.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Modeling the Yaw Behavior of Tail Fins for Small Wind Turbines: November 22, 2021 - May 21, 2024

This report describes the aerodynamic modelling of the yaw behavior of tail fins for small wind turbines (SWTs). The analysis is based on unsteady slender body theory (USBT) and is formulated to be used with the proposed tail fin module to be added to OpenFAST. Linearized USBT, which has been applied in the past to SWT tail fins, is developed in three main ways. First, USBT is extended to high yaw angles and aspect ratios by modelling the associated nonlinear vortex dynamics and the chordwise load distribution respectively. Second, we consider the effect of time-varying wind speed. The extended theory is compared to recent unpublished measurements at the University of Perugia, Italy, for the yaw behavior of delta, elliptical, and rectangular tail fins without a rotor and nacelle. The fins were released from initial yaw angles of -40 degrees and -80 degrees ; the latter is of sufficient magnitude to show the importance of the nonlinear yaw dynamics. Further, the friction in the model tail fin bearings was measured and modelled. The constants in the aerodynamic models were taken from the literature but significant improvement in accuracy is achieved by using system identification techniques to optimize the values. This leads to the third main finding: extended USBT is accurate for all generic tail fin shapes that we considered and is, therefore, very suitable for inclusion in aeroelastic codes for SWTs. The final development is of a general nonlinear equation for yaw response that is simpler than the extended USBT and is also suitable for inclusion in OpenFAST. The two nonlinear models are used to highlight the geometric requirements for good tail fin design.

distributed wind↗

Circuit Quantum Electrodynamics in Hyperbolic Space: From Photon Bound States to Frustrated Spin Models

Circuit quantum electrodynamics is one of the most promising platforms for efficient quantum simulation and computation. In recent groundbreaking experiments, the immense flexibility of superconducting microwave resonators was utilized to realize hyperbolic lattices that emulate quantum physics in negatively curved space. Here we investigate experimentally feasible settings in which a few superconducting qubits are coupled to a bath of photons evolving on the hyperbolic lattice. We compare our numerical results for finite lattices with analytical results for continuous hyperbolic space on the Poincaré disk. We find good agreement between the two descriptions in the long-wavelength regime. We show that photon-qubit bound states have a curvature-limited size. We propose to use a qubit as a local probe of the hyperbolic bath, for example, by measuring the relaxation dynamics of the qubit. We find that, although the boundary effects strongly impact the photonic density of states, the spectral density is well described by the continuum theory. We show that interactions between qubits are mediated by photons propagating along geodesics. We demonstrate that the photonic bath can give rise to geometrically frustrated hyperbolic quantum spin models with finite-range or exponentially decaying interaction.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

“Understanding Robustness Lottery”: A Geometric Visual Comparative Analysis of Neural Network Pruning Approaches

Deep learning approaches have provided state-of-the-art performance in many applications by relying on large and overparameterized neural networks. However, such networks are very brittle and are difficult to deploy on resource-limited platforms. Model pruning, i.e., reducing the size of the network, is a widely adopted strategy that can lead to a more robust and compact model. Many heuristics exist for model pruning, but our understanding of the pruning process remains limited due to the black-box nature of a neural network model. Empirical studies show that some heuristics improve performance whereas others can make models more brittle. Here, this work aims to shed light on how different pruning methods alter the network’s internal feature representation and the corresponding impact on model performance. To facilitate a comprehensive comparison and characterization of the high-dimensional model feature space, we introduce a visual geometric analysis of feature representations. We evaluated a set of critical geometric concepts decomposed from the commonly adopted classification loss and used them to design a visualization system to compare and highlight the impact of pruning on model performance and feature representation. The proposed tool provides an environment for an in-depth comparison of pruning methods and a comprehensive understanding of how the model responds to common data corruption. By leveraging the proposed visualization, machine learning researchers can reveal the similarities between pruning methods and redundancy in robustness evaluation benchmarks, obtain geometric insights about the differences between pruned models that achieve superior robustness performance, and identify samples that are robust or fragile to model pruning and common data corruption.

Li, Zhimin [Univ. of Utah, Salt Lake City, UT (Uni↗

FAIR Data and Interpretable AI Framework for Architectured Metamaterials (Final Report)

This research program established a transformative framework for the discovery and design of mechanical metamaterials, which are architected structures engineered to control physical phenomena like sound and vibration in ways natural materials cannot. To overcome the traditional reliance on trial-and-error, the project developed an interpretable Artificial Intelligence (AI) framework that moves beyond "black box" models to reveal the specific geometric patterns—such as "unit-cell templates"—that govern a material’s performance. A major breakthrough was the development of a hierarchical design method, which allows a single material to block vibrations across multiple frequency ranges simultaneously by layering patterns at different scales without them interfering with one another. This was further expanded to include irregular, graph-based designs that use spanning tree algorithms to ensure structural connectivity while allowing for customized, direction-dependent properties like stiffness and acoustic impedance. Beyond design, the project addressed the practicalities of real-world production by developing uncertainty quantification techniques that account for manufacturing defects and material variability, reducing the need for expensive physical testing by orders of magnitude. To speed up the discovery process, the team implemented Gaussian Process Regression and other surrogate models that provide accurate performance predictions at a fraction of the traditional computational cost. The AI-generated designs were successfully validated through fabrication of physical samples and wave propagation experiments, confirming their ability to accurately guide or reflect waves as predicted. By contributing these tools and high-quality FAIR benchmark datasets to the wider scientific community, this work provides a scalable foundation for advancing technologies in aerospace vibration control, medical imaging, and noise reduction.

36 MATERIALS SCIENCE↗

Experiment and non-local crystal plasticity finite element study of nanoindentation on Al-8Ce-10Mg alloy

Cerium and magnesium strengthened aluminum alloys, or Al-Ce-Mg, is a recently developed alloy family that exhibits good mechanical properties at elevated temperatures (~300 °C). To examine the single-crystal properties of Al-Ce-Mg alloys, nanoindentation experiments are conducted in this study. A crystal plasticity finite element model (CPFEM) with the evolution of geometrically necessary dislocations (GNDs) is applied to simulate the indentation in individual grains. A parametric study is carried out to investigate the sensitivity of each crystal plasticity model parameter to the indentation behavior. The highly sensitive parameters are calibrated by matching the indentation load-depth curves, while the rest parameters are obtained from bulk polycrystal uniaxial tension tests. Overall, satisfactory matching between experiment and simulation is obtained for each individual grain. The calculated hardness, as determined from the experiment, shows the dependence on indent depth, which is captured by the GND model. Furthermore, the effect of grain orientation and neighboring grains to nanoindentation behavior have been discussed with the comparison between the simulation and experiments.

36 MATERIALS SCIENCE↗

Geometric Scale-up Experiments on Fluidization of Geldart B Glass Beads

The objective of this work is to provide a valuable database from controlled experiments for validating computational models. Recently, coarse-grained techniques such as particle-in-cell (PIC) or coarse-grained discrete element modeling (DEM) have gained popularity due to their computational efficiency while modeling large-scale systems; however, the influence of model parameters and their sensitivities at different geometric scales and flow conditions remain to be analyzed. These datasets are critical for the multiphase flow research community to assess predictive capability of modeling techniques as well as elucidate the hydrodynamic behavior in these systems. This study performed fluidization experiments using three different test sections with internal diameters of 2.5, 4, and 6 in. The operating conditions, bed material, and range of flow velocities at the inlet were constant in all the units, which were not hydrodynamically scaled. Glass beads having a Sauter Mean Diameter of 332 μm were used. Superficial velocity was varied from 2.97 to 5.35 times the minimum fluidization velocity and the initial static bed height was 0.1524 m. The order in which the experiments were performed was randomized and replicates were included to provide uncertainty in measurements. Statistics of differential pressure and bed height from these tests were reported. Future plans include validating PIC methodology in the open-source software, MFiX (Multiphase Flow with Interphase Exchanges) using results from this study. This could further be extended to determine optimal model parameters using inverse techniques such as deterministic calibration or Bayesian inference.

20 FOSSIL-FUELED POWER PLANTS↗

Geometric Scale-up Experiments on Fluidization of Geldart B Glass Beads

The objective of this work is to provide a valuable database from controlled experiments for validating computational models. Recently, coarse-grained techniques such as particle-in-cell (PIC) or coarse-grained discrete element modeling (DEM) have gained popularity due to their computational efficiency while modeling large-scale systems; however, the influence of model parameters and their sensitivities at different geometric scales and flow conditions remain to be analyzed. These datasets are critical for the multiphase flow research community to assess predictive capability of modeling techniques as well as elucidate the hydrodynamic behavior in these systems. This study performed fluidization experiments using three different test sections with internal diameters of 2.5, 4, and 6 in. The operating conditions, bed material, and range of flow velocities at the inlet were constant in all the units, which were not hydrodynamically scaled. Glass beads having a Sauter Mean Diameter of 332 μm were used. Superficial velocity was varied from 2.97 to 5.35 times the minimum fluidization velocity and the initial static bed height was 0.1524 m. The order in which the experiments were performed was randomized and replicates were included to provide uncertainty in measurements. Statistics of differential pressure and bed height from these tests were reported. Future plans include validating PIC methodology in the open-source software, MFiX (Multiphase Flow with Interphase Exchanges) using results from this study. This could further be extended to determine optimal model parameters using inverse techniques such as deterministic calibration or Bayesian inference.

42 ENGINEERING↗

Parametric reduced order models for graded lattice structures

Graded lattice structures, characterized by smoothly varying mechanical properties, hold significant promise for optimizing material distribution in advanced engineering applications. However, accurately modeling these structures poses substantial computational challenges due to the continuous geometric variations within their unit cells. Here, to address these challenges, this paper introduces a novel Efficient Reduced Order Model (EROM) that integrates the Matrix Discrete Empirical Interpolation Method (MDEIM) and Discrete Empirical Interpolation Method (DEIM) with polynomial regression to efficiently manage geometric parametrization in lattice structures. Unlike traditional reduced order models (ROMs) that require extensive precomputed libraries for each geometric configuration, our approach enables continuous geometric variations through a flexible algebraic formulation, significantly reducing computational costs while preserving high accuracy. The method constructs projection matrices for individual unit cells that can be efficiently assembled into global systems, leveraging the repetitive nature of lattice structures. Numerical studies demonstrate that our EROM achieves displacement errors below 1% and von Mises stress prediction errors below 4%, coupled with computational speedups exceeding two orders of magnitude compared to full-order simulations. The proposed method's modularity and scalability make it particularly suitable for design optimization and real-time simulation of functionally graded lattice structures, with applications spanning aerospace to biomedical engineering.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

EWPD in the SMEFT to dimension eight

We calculate the O H † H 2 / Λ 4 corrections to LEP electroweak precision data using the geometric formulation of the Standard Model Effective Field Theory (SMEFT). We report our results in simple-to-use interpolation tables that allow the interpretation of this data set to dimension eight for the first time. We demonstrate the impact of these previously unknown terms in the case of a general analysis in the SMEFT, and also in the cases of two distinct models matched to dimension eight. Neglecting such dimension-eight corrections to LEP observables introduces a theoretical error in SMEFT studies. We report some preliminary studies defining such a theory error, explicitly demonstrating the effect of previously unknown dimension-eight SMEFT corrections on LEP observables.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

PQ axiverse

We show that the strong CP problem is solved in a large class of compactifications of string theory. The Peccei-Quinn mechanism solves the strong CP problem if the CP-breaking effects of the ultraviolet completion of gravity and of QCD are small compared to the CP-preserving axion potential generated by low-energy QCD instantons. We characterize both classes of effects. To understand quantum gravitational effects, we consider an ensemble of flux compactifications of type IIB string theory on orientifolds of Calabi-Yau hypersurfaces in the geometric regime, taking a simple model of QCD on D7-branes. We show that the D-brane instanton contribution to the neutron electric dipole moment falls exponentially in N 4 , with N the number of axions. In particular, this contribution is negligible in all models in our ensemble with N > 17. We interpret this result as a consequence of large N effects in the geometry that create hierarchies in instanton actions and also suppress the ultraviolet cutoff. We also compute the CP breaking due to high-energy instantons in QCD. In the absence of vectorlike pairs, we find contributions to the neutron electric dipole moment that are not excluded, but that could be accessible to future experiments if the scale of supersymmetry breaking is sufficiently low. The existence of vectorlike pairs can lead to a larger dipole moment. Finally, we show that a significant fraction of models are allowed by standard cosmological and astrophysical constraints.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Improved machine learning algorithm for predicting ground state properties

Finding the ground state of a quantum many-body system is a fundamental problem in quantum physics. In this work, we give a classical machine learning (ML) algorithm for predicting ground state properties with an inductive bias encoding geometric locality. The proposed ML model can efficiently predict ground state properties of an n-qubit gapped local Hamiltonian after learning from only $\mathcal{O}$(log(n)) data about other Hamiltonians in the same quantum phase of matter. This improves substantially upon previous results that require $\mathcal{O}$(n c ) data for a large constant c. Furthermore, the training and prediction time of the proposed ML model scale as $\mathcal{O}$(n log n) in the number of qubits n. Numerical experiments on physical systems with up to 45 qubits confirm the favorable scaling in predicting ground state properties using a small training dataset.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A generalized approach to x-ray data modeling for high-energy-density plasma experiments

Accurate understanding of x-ray diagnostics is crucial for both interpreting high-energy-density experiments and testing simulations through quantitative comparisons. X-ray diagnostic models are complex. Past treatments of individual x-ray diagnostics on a case-by-case basis have hindered universal diagnostic understanding. Here, in this study, we derive a general formula for modeling the absolute response of non-focusing x-ray diagnostics, such as x-ray imagers, one-dimensional space-resolved spectrometers, and x-ray power diagnostics. The present model is useful for both data modeling and data processing. It naturally accounts for the x-ray crystal broadening. The new model verifies that standard approaches for a crystal response can be good approximations, but they can underestimate the total reflectivity and overestimate spectral resolving power by more than a factor of 2 in some cases near reflectivity edge features. We also find that a frequently used, simplified-crystal-response approximation for processing spectral data can introduce an absolute error of more than an order of magnitude and the relative spectral radiance error of a factor of 3. The present model is derived with straightforward geometric arguments. It is more general and is recommended for developing a unified picture and providing consistent treatment over multiple x-ray diagnostics. Such consistency is crucial for reliable multi-objective data analyses.

Nagayama, Taisuke↗

Oak Ridge National Laboratory Building Envelope Library (ORNOBEL)

The Oak Ridge National Laboratory Building Envelope Library (ORNOBEL) is a collection of dense exterior building-facade point clouds acquired using a survey-grade terrestrial laser scanner. Each file represents an individual facade from a building on the Oak Ridge National Laboratory (ORNL) campus or in Knoxville, Tennessee, with an average point-cloud resolution of approximately 3 mm. The points in each facade are semantically labeled into three classes: (1) window/door, representing openings in the building envelope; (2) wall, representing planar opaque envelope surfaces; and (3) other, representing the remaining facade-adjacent elements, architectural features, and protrusions. ORNOBEL supports the development, training, and evaluation of advanced deep-learning methods for automated building-envelope segmentation, geometric reconstruction, and building information modeling (BIM).

Maldonado Puente, Bryan [ORNL] (ORCID:000000033880↗