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 865 records · Page 48

A computer program for modeling non-spherical eclipsing binary star systems

The accurate analysis of eclipsing binary light curves is fundamental to obtaining information on the physical properties of stars. The model described accounts for the important geometric and photometric distortions such as rotational and tidal distortion, gravity brightening, and reflection effect. This permits a more accurate analysis of interacting eclipsing star systems. The model is designed to be useful to anyone with moderate computing resources. The programs, written in FORTRAN 4 for the IBM 360, consume about 80k bytes of core. The FORTRAN program listings are provided, and the computational aspects are described in some detail.

Wood, D. B.↗

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↗

Influence of Coronal Abundance Variations

During the past year we have developed the multispecies code to include large scale flows and shown how those solutions compare with the observations of long lived loops. The time-dependent aspects of this work are under development. We will extend the steady flow investigations to study the effect these flows have on coronal structure as observed with TRACE. Coronal observations derive from heavy-ion emission; thus, we focus on the extent to which flow may modify coronal abundances by examining the heavy-ion abundance stratification within long-lived loops. We discuss the magnitudes of the physical effects modeled and compare simulated results with TRACE observations. These results can have a profound effect on the interpretation of TRACE observations. We will also apply the potential field matching programs developed by Aad Van Ballegooijen and Harvard undergraduate, Ellen Lee, to loops modeled for this project. Having the best possible understanding of the geometric properties of the loops is important for accurate modeling.

Gurman, Joseph↗

DNA Break Clustering as a Predictor of Cell Death across Various Radiation Qualities: Influence of Cell Size, Cell Asymmetry, and Beam Orientation

Cosmic radiation, composed of high charge and energy (HZE) particles, causes cellular DNA damage that can result in cell death or mutation that can evolve into cancer. In this work, a cell death model is applied to several cell lines exposed to HZE ions spanning a broad range of linear energy transfer (LET) values. We hypothesize that chromatin movement leads to the clustering of multiple double strand breaks (DSB) within one radiation-induced foci (RIF). The survival probability of a cell population is determined by averaging the survival probabilities of individual cells, which is function of the number of pairwise DSB interactions within RIF. The simulation code RITCARD was used to compute DSB. Two clustering approaches were applied to determine the number of RIF per cell. RITCARD outputs were combined with experimental data from four normal human cell lines to derive the model parameters and expand its predictions in response to ions with LET ranging from ∼0.2keV/μmto∼3000keV/μm. Spherical and ellipsoidal nuclear shapes and two ion beam orientations were modeled to assess the impact of geometrical properties on cell death. The calculated average number of RIF per cell reproduces the saturation trend for high doses and high-LET values that is usually experimentally observed. The cell survival model generates the recognizable bell shape of LET dependence for the relative biological effectiveness (RBE). At low LET, smaller nuclei have lower survival due to increased DNA density and DSB clustering. At high LET, nuclei with a smaller irradiation area either because of a smaller size or a change in beam orientation have a higher survival rate due to a change in the distribution of DSB/RIF per cell. If confirmed experimentally, the geometric characteristics of cells would become a significant factor in predicting radiation-induced biological effects.

cell survival↗

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↗

New capabilities of the MORET 6 Monte Carlo neutron transport code

The MORET code is a simulation tool that solves the transport equation for neutrons using the Monte Carlo method. It allows users to model complex three-dimensional geometrical configurations in a user-friendly way. New features have been introduced to extend the application field of MORET beyond the usual criticality calculations for which it has been initially designed. The most important change is the addition of an analog fixed source mode which allows studies of systems of any reactivity combined with very flexible outputs. Other useful improvements have been added concerning the geometric part, the fission matrix, the multigroup sensitivity coefficients and the outputs. This paper presents an overview of these new features. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Model Development and Analysis of a High-Fidelity Neutron Transport Sensor: The Quadrupole Detector Concept for Measurement of the Neutron Flux Gradient

Accurate reconstruction of the neutron flux distribution within a reactor core is essential for safe and efficient reactor operation. Traditional power shape synthesis in Light Water Reactors relies on hundreds of in-core detectors. However, this approach becomes impractical for Advanced Reactors and Microreactors due to limited space and harsh environments. To address this challenge, we propose a data-driven methodology that combines high-fidelity modeling with real-time ex-core sensor measurements, enabling the reconstruction of core power distribution while minimizing the reliance on intrusive in-core instrumentation. This project began in FY24 and achieved two initial milestones: (1) the definition of a three-year development plan for a Digital Twin framework and (2) the development of high-fidelity neutronics models of the Purdue University Reactor One (PUR-1) using both MCNP6 and OpenMC. The PUR-1 reactor, a zero-power facility, was selected due to its suitability for neutronics-focused modeling and the availability of experimental data for validation. Both models were benchmarked using neutron flux measurements obtained from irradiated gold foils, which were strategically placed within the core during a dedicated campaign in July 2024. This report marks the continuation and completion of those foundational tasks. The OpenMC model has been refined (improved geometric accuracy, expanded cross-section libraries, and refined sampling) and validated using additional experimental data. An updated sensor design—based on quadrupole configuration—was designed to measure both ex-core flux and its spatial gradient. These measurements will serve as inputs to a neural network-based reconstruction algorithm. Finally, the methodology was demonstrated on a two-dimensional test case representative of the heterogeneous material composition of the PUR-1 reactor core. A neural network implementation of the Kirchhoff-Helmholtz integral equation was employed to solve the boundary value problem using peripheral sensor measurements. The preliminary results confirm the strong potential of the proposed approach for accurate and minimally invasive neutron flux reconstruction.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Performance Comparison at Mach Numbers 1.8 and 2.0 of Full Scale and Quarter Scale Translating-Spike Inlets

The performance of a full-scale translating-spike inlet was obtained at Mach numbers of 1.8 and 2.0 and at angles of attach from 0 deg to 6 deg. Comparisons were made between the full-scale production inlet configuration and a geometrically similar quarter-scale model. The inlet pressure-recovery, cowl pressure-distribution, and compressor-face distortion characteristics of the full-scale inlet agreed fairly well with the quarter-scale results. In addition, the results indicated that bleeding around the periphery ahead of the compressor-face station improved pressure recovery and compressor-face distortion, especially at angle of attack.

Anderson, B. H.↗

Conjugate point determination for multitemporal data overlay

The machine processing is discussed of spatially variant multitemporal data such as imagery obtained at different times which requires that these data be in geometrical registration so that the analysis processor may obtain the datum for a specified ground resolution element in each of the sets of imagery being utilized for analysis. Misregistration between corresponding subsets of imagery contains both a displacement and a geometrical distortion component. Search techniques utilizing the moduli of the Fourier Transforms are developed for estimating the coefficients of geometrical distortion components of this model. Following the correction of these distortion components, the displacement is located by the crosscorrelation of a template obtained from one set of data, termed the reference, with the second, or background data. This template, derived for the optimum discrimination of the reference data embedded in the background, is determined by the solution of a system of equations involving the reference data and the covariance matrix of these data.

Emmert, R. A.↗

Transonic flutter study of a 50.5 deg cropped-delta wing with two rearward-mounted nacelles

Transonic flutter characteristics of three geometrically similar delta-wing models were experimentally determined in the Langley transonic dynamics tunnel at Mach numbers from about 0.6 to 1.2. The models were designed to be simplified versions of an early supersonic transport wing design. The model was an aspect-ratio-1.28 cropped-delta wing with a leadingedge sweep of 50.5 deg. The flutter characteristics obtained for this wing configuration indicated a minimum flutter-speed index near a Mach number of 0.92 and a transonic compressibility dip amounting to about a 27-percent decrease in the flutter-speed index relative to the value at a Mach number of 0.6. Analytical studies were performed for one wing model at Mach numbers of 0.6, 0.7, 0.8, and 0.9 by using both doublet-lattice and lifting-surface (kernel-function) unsteady aerodynamic theory. A comparison of the analytical and experimental flutter results showed good agreement at all Mach numbers investigated.

Sandford, M. C.↗

Application of the CINGEN program a thermal network data generator

The application of the CINGEN computer program and two of its supporting programs for the evaluation of structural and thermal performance of physical systems was described. The CINGEN program was written and implemented to avoid the duplication effort of performing a finite element approach for structural analysis and a finite differencing technique for thermal analysis, as well as the desire for a geometrical representation of the thermal model to reduce modeling errors. The program simplifies the thermal modeling process by performing all of the capacitance and conductance calculations normally done by the analyst. Each solid element is divided into five tetrahedrons, allowing the total volume to be calculated precisely. A sample problem was illustrated.

Shultz, W. E.↗

Vertical velocity structure and geometry of clear air convective elements

The paper discusses observations of individual convective elements with a high-power narrow-beam scanning radar, an FM-CW radar, and an acoustic sounder, including the determination of the vertical air velocity patterns of convective structures with the FM-CW radar and acoustic sounder. Data are presented which link the observed velocity structure and geometrical patterns to previously proposed models of boundary layer convection. It is shown that the high-power radar provides a clear three-dimensional picture of convective cells and fields over a large area with a resolution of 150 m, where the convective cells are roughly spherical. Analysis of time-height records of the FM-CW radar and acoustic sounder confirms the downdraft-entrainment mechanism of the convective cell. The Doppler return of the acoustic sounder and the insect-trail slopes on FM-CW radar records are independent but redundant methods for obtaining the vertical velocity patterns of convective structures.

Rowland, J. R.↗