Engineering Papers⌕ Search

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.

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

At least 937 records · Page 52
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Using Thermodynamic Availability to Predict the Transitional Film Reynolds Number between the Droplet and Jet Modes in Falling Liquid between Horizontal Tubes

The objective of this paper is to propose and validate a model for predicting the transitional film Reynolds number between the droplet and jet modes of falling film flows. The proposed model relies on thermodynamic availability with no fitting to mode transition data. The model directly predicts the numerical values of the transitional film Reynolds numbers at different operating conditions. The results of the model are validated against 43 experimental data points from the literature, which cover 10 working fluids and a wide range of tube geometries. The accuracy of the proposed model in predicting the validation data is -on average- more than 20 % better than that of the widely used Re = a (Ga) b and Re = a (Ga) b (S/ξ) c models. The model also achieves low computation time (less than 30 seconds for the 43 validation data points). The model presented in this paper has the potential of integration with heat transfer, mass transfer, and pressure drop models for optimizing the design and operation of falling film heat and mass exchangers.

42 ENGINEERING↗
DOE OSTI2025

Market Driven Residential Energy Codes: Comparing Performance in a Changing Technological Environment

The research project is undertaken to better understand the changing relationship between the two basic methods of building energy code compliance – prescriptive and performance – and how those methods relate to each other with respect to advancements in building energy computer simulation standards and capabilities. The International Energy Efficiency Code (IECC) is a model code adopted by many jurisdictions across the United States. Historically, the prescriptive compliance methodology has been preferred in most jurisdictions. The prescriptive methodology requires meeting or exceeding specific efficiency minimums for each envelope component. This tends to be a simple method to teach and verify. A more involved prescriptive alternative called the Total UA alternative is sometimes used. This method requires some multiplication, summing, and comparison to compute, so it is done with a fairly simple computer program. However, advances in computer and building energy simulation technology have resulted in increased use of more detailed performance compliance methods. The performance compliance method establishes the annual energy cost threshold via hourly simulation models. The compliance threshold is determined with a comparison building model simulation with geometry similar to the proposed home and with energy feature parameters and efficiencies as specified in the IECC. This project examines relationships between the two methods of building energy code compliance, including: • Overall annual energy use based on utility bill analysis by compliance method • Code official work processes with respect to compliance methods • Gaps and issues associated with building code compliance methods • Simulated energy use difference between compliance methods • Code compliance cost as a function of compliance method • Code compliance labeling effectiveness for high performance residences • Getting to net zero energy use and net zero greenhouse gas emissions through high performance code alternatives • Electronic code permitting and compliance alternatives

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2025

Doppler Broadening and Other Temperature Effects

In this paper I attempt to document what I have learned and still remember about temperature effects, particularly Doppler broadening, on nuclear data, as it applies to ENDF formatted data. My focus is on the SIGMA1 method of Doppler broadening that I developed over 50 years ago. However, I do realize that there are currently many different methods used in computer codes, and I have tried to keep my discussion general as it applies to ALL of these methods. In particular, ALL of the MYTHS I describe below apply to all Doppler broadening methods: they all solve the parabolic diffusion equation in spherical geometry in (speed, reaction rate) versus temperature. These methods differ only in how the cross section is represented. So, at 85 years old here is what I still remember.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Progress Towards NQA-1 for Cardinal in FY25

Cardinal is a wrapping of the GPU-oriented spectral element Computational Fluid Dynamics (CFD) code NekRS and the Monte Carlo particle transport code OpenMC within the Multiphysics Object-Oriented Simulation Environment (MOOSE). Cardinal provides high-resolution thermal-hydraulics and/or radiation transport feedback to MOOSE multiphysics simulations. Multiphysics feedback is implemented in a geometry-agnostic manner which eliminates the need for rigid one-to-one mappings. A generic data transfer implementation also allows NekRS and OpenMC to couple to any MOOSE application, enabling a broad set of multiphysics capabilities. Cardinal simulations can also leverage combinations of MPI, OpenMP, and GPU resources. Cardinal continuous development and improvement efforts have led to the software being considered as a high-fidelity design and licensing tool for key areas of nuclear reactor relevant physics, including neutron transport, fluid flow, heat transfer, and mechanical processes. The fast development and expansion of the software from a pure R&D framework towards its application in the nuclear industry and regulation require a focus on developing, enhancing,and maintaining Cardinal’s software quality through strict adherence to a Software Quality Assurance (SQA) framework and SQA program. To facilitate compliance with SQA standards, the Cardinal SQA Program was initiated during Fiscal Year 2023 (FY23). During the development of the Cardinal SQA Program, multiple gaps have been identified. These gaps are primarily related to model verification and code pedigree as they relate to the use of Cardinal as an analysis tool. These gaps were captured in a report published in 2023. A second report highlighted the progress made during Fiscal Year 2024 (FY24) and described Argonne’s effort to document and integrate software verification within Cardinal’s software development process. This report documents the progress made towards NQA-1 for Cardinal in the Fiscal Year 2025 (FY25). All cases in the expanded Continuous Integration (CI) suite of NekRS are included in this report which test the solvers and modules available in NekRS exhaustively. The NekRS tests are integrated with the Cardinal CI suite and made available in publicly accessible Github documentation. Following the CI practice permits integrating of source code changes frequently and ensuring that the integrated codebase clears the verification testing for the software. Also in this report is a brief overview of the development of the Cardinal Software Quality Assurance Plan (SQAP) that was done in FY25, though it should be noted that the rest of the documentation for the SQA program needs to be developed in a future step of this task.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Design and Flow Considerations of Additively Manufactured, Internal Cooling Geometries for Small Industrial Gas Turbines

Additive manufacturing is now a mainstream technology and can be utilized to rapidly develop and test turbine airfoil cooling networks. This paper reports on an ongoing effort to integrate advanced internal cooling architectures in a realistic blade profile for test in a high-speed cascade. Airfoil cooling schemes were developed using reduced order modeling and computer aided design. However, the additive manufacturing impacts on cooling channel flow performance were unknown. Test articles consisting of typical cooling features and networks were derived from the designs and flow proved to identify additive manufacturing impacts on performance and develop guidelines to mitigate these impacts.

additive manufacturing↗
DOE OSTI2020

Pattern formation in a coupled membrane-bulk reaction-diffusion model for intracellular polarization and oscillations

Reaction-diffusion systems have been widely used to study spatio-temporal phenomena in cell biology, such as cell polarization. Coupled bulk-surface models naturally include compartmentalization of cytosolic and membrane-bound polarity molecules. Here we study the distribution of the polarity protein Cdc42 in a mass-conserved membrane-bulk model, and explore the effects of diffusion and spatial dimensionality on spatio-temporal pattern formation. We first analyze a one-dimensional (1-D) model for Cdc42 oscillations in fission yeast, consisting of two diffusion equations in the bulk domain coupled to nonlinear ODEs for binding kinetics at each end of the cell. In 1-D, our analysis reveals the existence of symmetric and asymmetric steady states, as well as anti-phase relaxation oscillations typical of slow-fast systems. We then extend our analysis to a two-dimensional (2-D) model with circular bulk geometry, for which species can either diffuse inside the cell or become bound to the membrane and undergo a nonlinear reaction-diffusion process. We also consider a nonlocal system of PDEs approximating the dynamics of the 2-D membrane-bulk model in the limit of fast bulk diffusion. In all three model variants we find that mass conservation selects perturbations of spatial modes that simply redistribute mass. In 1-D, only anti-phase oscillations between the two ends of the cell can occur, and in-phase oscillations are excluded. In higher dimensions, no radially symmetric oscillations are observed. Instead, the only instabilities are symmetry-breaking, either corresponding to stationary Turing instabilities, leading to the formation of stationary patterns, or to oscillatory Turing instabilities, leading to traveling and standing waves. Codimension-two Bogdanov—Takens bifurcations occur when the two distinct instabilities coincide, causing traveling waves to slow down and to eventually become stationary patterns. Our work clarifies the effect of geometry and dimensionality on behaviors observed in ma.ss-conserved cell polarity models.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Data Projection of the High Temperature Electrolysis System in the Dynamic Energy Transport and Integration Laboratory using Dynamic System Scaling

For nuclear power to be flexible in a functioning Integrated Energy System (IES), excess produced heat must be stored or utilized during times of low power demand to ensure a load factor of 1 while load balancing. The Dynamic Energy Transport and Integration Laboratory (DETAIL) is one facility that is under development to emulate IES conditions on the engineering-scale, planned to conduct virtual real time operations with industry-scale facilities, and is currently testing thermal storage and high temperature electrolysis. As part of the study to develop a method to preprocess input signals or postprocess output signals between systems of different scales via Dynamical System Scaling (DSS), the current research is one of the continued efforts branching from the data projection activity conducted for the Thermal Energy Distribution System and currently engages the High Temperature Electrolysis (HTE) System in DETAIL. The HTE SOEC electrical, fluid, and thermal dynamics Figure of Merits (FOM) were identified, governing equations and closure relations were successfully scaled, and relations between FOM scaling ratios were determined. Setting the scaling objectives to reform existing data to project a data set that doubly accelerated the electrolysis process while preserving the produced amount of hydrogen was generated for the full transient. The calculated boundary conditions were inlet temperature, stack current, and inlet steam mass flow rate at 1470 K, 121.1 A, and 1.886 g/s, respectively. The research outcomes demonstrated an output signal postprocessing case accelerating the hydrogen production without changing geometry, number of cells, and partial pressures.

08 HYDROGEN↗
DOE OSTI2021

Validation of Oregon State University High Temperature Test Facility Experiments Using Pronghorn

The OSU High Temperature Test Facility is a quarter-scale diameter, 1/64 scale volume test facility meant to replicate thermophysical phenomena in the prototypical General Atomics Modular High Temperature Gas Reactor. Tests pertaining to conduction cooldown events were performed from 2016-2019, providing a large database by which computational methods that are applicable to different length scales can be validated. One of these codes is Pronghorn, which is a coarse-meshed, porous-based subchannel thermal hydraulics code based on the MOOSE application, with the intention of better capturing the physics of both conduction and convection heat transfer within the OSU HTTF core. The goal of this summer project is to develop the framework by which Pronghorn can perform validation exercises of the HTTF core for benchmarking, by generating a mesh appropriate to the geometry of the HTTF core, developing input decks that accurately capture the initial and boundary conditions, materials, and relevant equations to the physics seen in the HTTF core, and using a postprocessor to compare simulation results to various experimental data. While validation of codes is a multi-year project, a mesh has been generated and tested in Pronghorn that meets mass conservation and basic heat transfer principles. The next step is to accurate depict the fluid inlet and outlet boundary conditions, which will be performed using computational fluid dynamics software.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2026

Identifying the Accuracy of Surface Roughness for Metal Additive Manufacturing Parts Captured with Computed Tomography Scanning

Conventional surface roughness measurement techniques require direct access to internal surfaces, often necessitating destructive sectioning of test articles. Computed tomography (CT) offers a non-destructive alternative for internal surface characterization, but its application in metrology remains unstandardized and sensitive to machine resolution and operator technique. This study investigates the feasibility of using CT scanning to quantify areal surface roughness in metal AM heat exchanger tubes with three distinct internal geometries: ribbed, discrete W, and featureless. CT-derived surface parameters were extracted using custom Python scripts and compared to measurements obtained from a focus variation microscope calibrated against a known standard. Results show that CT-based roughness measurements closely matched microscope values for the discrete W specimen, deviations between CT and microscope measurements were minimal—less than 1 µm—indicating reliable reconstruction. In contrast, the ribbed and featureless specimens, with lower roughness values showed greater discrepancies. The findings suggest that CT scanning can be a viable non-destructive metrology tool for AM parts with surface roughness above approximately 8 µm. For smoother surfaces, current CT capabilities may not provide sufficient accuracy, highlighting the need for resolution-aware workflows and further standardization in CT-based surface metrology.

36 MATERIALS SCIENCE↗
DOE OSTI2022

User’s Manual for RESRAD-BUILD Code V.4: Vol. 1 – Methodology and Models Used in RESRAD-BUILD Code

The RESRAD-BUILD computer code models radionuclide release and transport in indoor environments and performs pathway analyses to evaluate the potential radiological dose and risk incurred by an individual who works or lives in a building contaminated with radioactive material or housing radioactively contaminated furniture or equipment. The code provides four geometries to characterize a radiation source: point, line, area, and volume, in which radionuclides are homogeneously distributed. Radionuclides contained in a source are considered to be released to the indoor air due to various processes including erosion (mechanically or weathering), diffusion (for tritium and radon in a volume source), or emanation (radon in a point, line, or area source). The release can proceed through different time phases with different rates. In RESRAD-BUILD Version 4.0, a dynamic ventilation model is implemented to simulate the fate and transport of source material particles and radionuclides after their releases. This dynamic ventilation model considers (1) air exchange between rooms in the building and between the rooms and the outdoor environment, (2) deposition from air to floor, (3) resuspension from the floor to the air, and (4) periodical vacuuming that reduces the floor deposition. The fate and transport modeling provides estimates of radionuclide concentrations in the source, in the air, and on the floor at different times, which are then integrated over the exposure duration for the calculation of radiation doses and cancer risks. A single run of the RESRAD-BUILD code can model a building with up to 9 rooms, 10 sources, and 10 receptors. The potential radiation dose and cancer risk incurred by each receptor are calculated for seven exposure pathways: (1) external radiation directly from the sources (accounting for shielding), (2) external radiation from radioactive particles deposited on the floors, (3) external radiation from airborne radionuclides, (4) inhalation of airborne radionuclides, (5) inhalation of radon and radon progenies, (6) inadvertent ingestion of radioactive particles directly from the source, and (7) ingestion of radioactive particles deposited on the floors. Various exposure scenarios can be modeled with RESRAD-BUILD, including but are not limited to, office worker, renovation worker, decontamination worker, building visitor, and resident. Both deterministic and probabilistic analyses can be performed to obtain results in both text reports and graphic displays.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2020

Chaos and complementarity in de Sitter space

We consider small perturbations to a static three-dimensional de Sitter geometry. For early enough perturbations that satisfy the null energy condition, the result is a shockwave geometry that leads to a time advance in the trajectory of geodesics crossing it. This brings the opposite poles of de Sitter space into causal contact with each other, much like a traversable wormhole in Anti-de Sitter space. In this background, we compute out-of-time-order correlators (OTOCs) to asses the chaotic nature of the de Sitter horizon and find that it is maximally chaotic: one of the OTOCs we study decays exponentially with a Lyapunov exponent that saturates the chaos bound. We discuss the consequences of our results for de Sitter complementarity and inflation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Peptide nucleic acid Hoogsteen strand linker design for major groove recognition of DNA thymine bases

Sequence-specific targeting of double-stranded DNA and non-coding RNA via triple-helix-forming peptide nucleic acids (PNAs) has attracted considerable attention in therapeutic, diagnostic and nanotechnological fields. An E-base (3-oxo-2,3-dihydropyridazine), attached to the polyamide backbone of a PNA Hoogsteen strand by a side-chain linker molecule, is typically used in the hydrogen bond recognition of the 4-oxo group of thymine and uracil nucleic acid bases in the major groove. We report on the application of quantum chemical computational methods, in conjunction with spatial constraints derived from the experimental structure of a homopyrimidine PNA·DNA-PNA hetero-triplex, to investigate the influence of linker flexibility on binding interactions of the E-base with thymine and uracil bases in geometry-optimised model systems. Hydrogen bond formation between the N2 E-base atom and target pyrimidine base 4-oxo groups in model systems containing a β-alanine linker (J Am Chem Soc 119:11116, 1997) was found to incur significant internal strain energy and the potential disruption of intra-stand aromatic base stacking interactions in an oligomeric context. Additionally, in geometry-optimised model systems containing a 3-trans olefin linker (Bioorg Med Chem Lett 14:1551, 2004) the E-base swung out away from the target pyrimidine bases into the solvent. These findings are in qualitative agreement with calorimetric measurements in hybridisation experiments at T–A and U–A inversion sites. In contrast, calculations on a novel 2-cis olefin linker design indicate that it could permit simultaneous E-base hydrogen bonding with the thymine 4-oxo group, circumvention and solvent screening of the thymine 5-methyl group, and maintenance of triplex intra-stand base stacking interactions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

A high order Cartesian grid, finite volume method for elliptic interface problems

We present a higher-order finite volume method for solving elliptic PDEs with jump conditions on interfaces embedded in a 2D Cartesian grid. Second, fourth, and sixth order accuracy is demonstrated on a variety of tests including problems with high-contrast and spatially varying coefficients, large discontinuities in the source term, and complex interface geometries. We include a generalized truncation error analysis based on cell-centered Taylor series expansions, which then define stencils in terms of local discrete solution data and geometric information. In the process, we develop a simple method based on Green's theorem for computing exact geometric moments directly from an implicit function definition of the embedded interface. This approach produces stencils with a simple bilinear representation, where spatially-varying coefficients and jump conditions can be easily included and finite volume conservation can be enforced.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Local prediction of Laser Powder Bed Fusion porosity by short-wave infrared imaging thermal feature porosity probability maps

We report that local thermal history can significantly vary in parts during metal Additive Manufacturing (AM), leading to local defects. However, the sequential layer-by-layer nature of AM facilitates in-situ part voxelmetric observations that can be used to detect and correct these defects for part qualification and quality control. The challenge is to relate this local radiometric data with local defect information to estimate process error likelihood in future builds. This paper uses a Short-Wave Infrared (SWIR) camera to record the temperature history for parts manufactured with Laser Powder Bed Fusion (LPBF) processes. The porosity from a cylindrical specimen is measured by ex-situ micro-computed tomography (μCT). Specimen data from the SWIR camera, combined with the μCT data, are used to generate thermal feature-based porosity probability maps. The porosity predictions made by various SWIR thermal feature-porosity probability maps of a specimen with a complex geometry are scored against the true porosity obtained via μCT. The receiver operating characteristic curves constructed from the predictions for the complex sample demonstrate the porosity probability mapping methodology’s potential for in-situ based porosity detection.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Modelling the nonlinear plasma response to externally applied three-dimensional fields with the Stepped Pressure Equilibrium Code

Small-amplitude, symmetry-breaking magnetic field perturbations, including resonant magnetic perturbations (RMPs) and error fields, can profoundly impact plasma properties in both tokamaks and stellarators. In this work, we perform the first comparison between the Stepped Pressure Equilibrium Code (SPEC) (a comparatively fast and efficient equilibrium code based on energy-minimisation principles) and M3D-C 1 (a high-fidelity albeit computationally expensive initial-value extended-magnetohydro-dynamic (MHD) code) to assess the conditions under which SPEC can be used to model the nonlinear, non-ideal plasma response to an externally applied (m = 2,n = 1) RMP field in an experimentally relevant geometry. We find that SPEC is able to capture the plasma response in the weakly nonlinear regime – meaning perturbation amplitudes below the threshold for break up of the separatrix and onset of secondary magnetic island formation – when around half of the total toroidal flux is enclosed in the volume containing the q = 2 resonant surface. Further, the observed dependence of SPEC solutions on input parameters, including toroidal flux and the number of volumes into which the plasma is partitioned, indicates that additional exploration of the underlying Multi-Region Relaxed MHD physics model is needed to constrain the choice of parameters. Nonetheless, this work suggests promising applications of SPEC to optimisation and fusion plasma design.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Accelerating FEM-Based Corrosion Predictions Using Machine Learning

Atmospheric corrosion of metallic parts is a widespread materials degradation phenomena that is challenging to predict given its dependence on many factors (e.g. environmental, physiochemical, and part geometry). For materials with long expected service lives, accurately predicting the degree to which corrosion will degrade part performance is especially difficult due to the stochastic nature of corrosion damage spread across years or decades of service. The Finite Element Method (FEM) is a computational technique capable of providing accurate estimates of corrosion rate by numerically solving complex differential Eqs. characterizing this phenomena. Nevertheless, given the iterative nature of FEM and the computational expense required to solve these complex equations, FEM is ill-equipped for an efficient exploration of the design space to identify factors that accelerate or deter corrosion, despite its accuracy. In this work, a machine learning based surrogate model capable of providing accurate predictions of corrosion with significant computational savings is introduced. Specifically, this work leverages AdaBoosted Decision trees to provide an accurate estimate of corrosion current per width given different values of temperature, water layer thickness, molarity of the solution, and the length of the cathode for a galvanic couple of aluminum and stainless steel.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Evaluation of DED and LPBF Fe-based Alloys Process Application Envelopes based on Performance, Process Economics, Supply Chain Risks, and Reactor-specific Targeted Components

The U.S. Department of Energy (DOE), Office of Nuclear Energy (NE), Advanced Materials and Manufacturing Technologies (AMMT) program aims to develop extreme-environment materials solutions for use in the deployment of advanced nuclear reactors and the sustainment of the current fleet. To achieve this objective, a combination of experiment, a computational tool, and machine learning (ML) for the design of materials is adopted for the maturation of materials for nuclear technology. Through advanced manufacturing techniques such as laser powder bed fusion (LPBF) and laser powder direct energy deposition (LP-DED), components with complex geometries can be fabricated with reduced time and effort. Such advanced manufacturing methods can also provide the opportunity to improve materials performance through optimized microstructures and mechanical properties. However, existing engineering alloys are not always well suited for fabrication with additive manufacturing (AM), as their compositions have been tuned to optimize fabrication via conventional methods. Thus, similar alloys with modified compositions that are better suited for AM can be studied for improved performance. Over the past three years, the AMMT teams from Argonne National Laboratory (ANL) and Pacific Northwest National Laboratory (PNNL) studied various known Fe-based alloys by evaluating their initial printability using LPBF, and an AMMT-developed down-selection and decision matrix reduced the number of alloys to be studied from six to three in fiscal year (FY) 2024. Additionally, in FY 2024, for parallel evaluation, these three alloys were studied using LPDED. While LPBF is better for small- to medium-sized components with high detail and internal features, LP-DED combines a material feed system to place the powder onto the exact spot where the laser will melt the material. This AM method can be easily scaled to extremely large components and provides high build rate speeds compared to those of conventional LPBF systems. Additionally, DED is a better choice for complex geometries and compositional gradients.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Towards Exascale Astrophysics of Mergers and Supernovae (TEAMS)

The TEAMS project brought together cutting-edge simulations, theoretical insights, and collaborative efforts to deepen our understanding of some of the universe’s most extreme phenomena—supernovae, neutron star mergers, and the powerful signals they emit. Using one of the largest suites of 3D supernova simulations ever conducted, researchers uncovered new insights into how massive stars explode, how those explosions vary by stellar mass, and what conditions lead to the birth of neutron stars or black holes. They also studied the radiation and gravitational wave signals emitted during these events, revealing how future observations can be used to uncover what happens deep inside collapsing stars. The team developed improved tools for modeling how light and neutrinos behave in such explosive environments, enabling more accurate predictions of what astronomers might observe. Work also explored how the chemical composition and geometry of kilonovae—the visible explosions that follow neutron star mergers—influence their signals and can reveal the origins of heavy elements like gold. These efforts not only advanced scientific knowledge, but also trained a new generation of researchers at the intersection of astrophysics, computational science, and nuclear theory.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗