Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “differential 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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
NASA NTRS2010

Oscillating Cell Culture Bioreactor

To better exploit the principles of gas transport and mass transport during the processes of cell seeding of 3D scaffolds and in vitro culture of 3D tissue engineered constructs, the oscillatory cell culture bioreactor provides a flow of cell suspensions and culture media directly through a porous 3D scaffold (during cell seeding) and a 3D construct (during subsequent cultivation) within a highly gas-permeable closed-loop tube. This design is simple, modular, and flexible, and its component parts are easy to assemble and operate, and are inexpensive. Chamber volume can be very low, but can be easily scaled up. This innovation is well suited to work with different biological specimens, particularly with cells having high oxygen requirements and/or shear sensitivity, and different scaffold structures and dimensions. The closed-loop changer is highly gas permeable to allow efficient gas exchange during the cell seeding/culturing process. A porous scaffold, which may be seeded with cells, is fixed by means of a scaffold holder to the chamber wall with scaffold/construct orientation with respect to the chamber determined by the geometry of the scaffold holder. A fluid, with/without biological specimens, is added to the chamber such that all, or most, of the air is displaced (i.e., with or without an enclosed air bubble). Motion is applied to the chamber within a controlled environment (e.g., oscillatory motion within a humidified 37 C incubator). Movement of the chamber induces relative motion of the scaffold/construct with respect to the fluid. In case the fluid is a cell suspension, cells will come into contact with the scaffold and eventually adhere to it. Alternatively, cells can be seeded on scaffolds by gel entrapment prior to bioreactor cultivation. Subsequently, the oscillatory cell culture bioreactor will provide efficient gas exchange (i.e., of oxygen and carbon dioxide, as required for viability of metabolically active cells) and controlled levels of fluid dynamic shear (i.e., as required for viability of shear-sensitive cells) to the developing engineered tissue construct. This bioreactor was recently utilized to show independent and interactive effects of a growth factor (IGF-I) and slow bidirectional perfusion on the survival, differentiation, and contractile performance of 3D tissue engineering cardiac constructs. The main application of this system is within the tissue engineering industry. The ideal final application is within the automated mass production of tissue- engineered constructs. Target industries could be both life sciences companies as well as bioreactor device producing companies.

Freed, Lisa E.↗
DOE OSTI2026

When Do Band Gap Calculations Agree with Experiments in Monolayer-Protected Cu 14 and Au 20 Atomically Precise Nanoclusters? A (TD)-DFT Comparison of HOMO–LUMO, Fundamental, Optical, and Electrochemical Energy Gaps

In view of the tremendous progress in atomically precise metal nanoclusters where electrochemical and optical energetics are routinely supported by computations to establish structure−function correlations, we explore the relationship between different protocols for measuring and computing band gaps of two distinct organic ligand-protected nanoclusters: [Cu 14 H 10 (MBN) 3 (PPh 3 ) 8 ] + and Au 20 (TBBT) 16 . Through UV/visible spectroscopy and differential pulse voltammetry, we measure optical and electrochemical band gaps in those systems. We then compare these experimentally determined gaps to HOMO−LUMO gaps, fundamental gaps, vertical excitation energies, and E o ox − E o red potentials computed using different density functional theory (DFT) or time-dependent DFT (TDDFT) methods. Specifically, in both copper and gold nanoclusters, we test the effect of truncating inert ligands from the model and compare density functionals with varying degrees of Hartree−Fock (HF) exchange from 0 to 50%, range-separated hybrids with a varying long-range tuning parameter, different correlation functionals, basis sets, and (equilibrium and nonequilibrium) continuum solvation models. Despite having different frontier orbital characters (the copper nanocluster has a metal-to-ligand charge transfer character while the gold nanocluster has metal-centered frontier orbitals), both nanoclusters display a similar sensitivity of the HOMO−LUMO gap to the HF exchange that is partially mitigated when computing the fundamental, optical, and electrochemical gaps. Other factors, such as the nature of the correlation functional, basis set, and geometry relaxation, have a considerably smaller effect on computed band gaps in these systems. Overall, this work provides guidelines for factors of varied importance for correlating computed and experimental band gap values.

Chemical calculations↗
DOE OSTI2025

X-ray spectroscopy of multi-temperature plasmas using the differential emission measure formalism

We present a theoretical construct that nominally underlies spectroscopic data analysis of multi-temperature plasmas, known as the differential emission measure (DEM). From a data analytic perspective, the DEM formalism is used to derive temperature distributions from line spectra that are formed in the presence of temperature gradients and by time integrations of evolving plasmas. From a modeling perspective, DEMs are convenient intermediaries between radiation hydrodynamics simulations and spectroscopic measurements acquired in the laboratory. The DEM concept and its associated methodologies were originally developed by spectroscopists working with astrophysical data. We borrow from these earlier investigations. In this manuscript, intended primarily as a tutorial, we discuss the basic concepts, but also augment various aspects of the theory by the way of extension and example, including a detailed treatment of various weighting and averaging schemes, intended to mitigate ambiguities that often arise when reporting temperature information. We focus on high-temperature plasmas that are not in local thermodynamic equilibrium and the x-ray spectra that they produce, although the core ideas presented here are applicable to spectroscopy in other energy bands. A few examples involving the derivation and manipulation of model DEMs in simple geometries are provided.

Liedahl, Duane A. [Lawrence Livermore National Lab↗
NASA NTRS1993

Al-26 production profile and model comparisons in Canyon Diablo

The large preatmospheric size of the Canyon Diablo meteorite, a radius of about 15 m, makes it especially suitable for systematic studies of cosmogenic nuclide production rates of iron objects in a 2 pi geometry. To reconstruct the exposure history of the meteoroid, Heymann et al. investigated several fragments recovered from known geographic locations around the crater for their shock features and cosmogenic nobel gases. They applied the Signer-Nier noble gas production rate model to establish the preatmospheric depth of the specimens in the meteoroid. Cosmic ray exposure ages suggested a multi-episodic irradiation, with 170 or 540 Ma being inferred for most of the samples studied while two anomalous specimens indicated a possible third exposure age at 940 Ma. Be-10 and Cl-36 have been measured in a number of these same samples by accelerator mass spectrometry (AMS), with use being made of the preatmospheric depths determined in Heymann et al. to construct production profiles. The present study extends the cosmogenic radionuclide data to Al-26 and compares the results with both the production rate model of Reedy and Arnold and production rates determined from the cross sections used by the Reedy-Arnold model (for the major nuclear reactions making Al-26) in combination with differential fluxes calculated using the Los Alamos High Energy Transport (LAHET) Code System. Model calculations for Be-10 and Cl-36 have also been obtained.

Michlovich, E.↗
NASA NTRS2009

T-Slide Linear Actuators

T-slide linear actuators use gear bearing differential epicyclical transmissions (GBDETs) to directly drive a linear rack, which, in turn, performs the actuation. Conventional systems use a rotary power source in conjunction with a nut and screw to provide linear motion. Non-back-drive properties of GBDETs make the new actuator more direct and simpler. Versions of this approach will serve as a long-stroke, ultra-precision, position actuator for NASA science instruments, and as a rugged, linear actuator for NASA deployment duties. The T slide can operate effectively in the presence of side forces and torques. Versions of the actuator can perform ultra-precision positioning. A basic T-slide actuator is a long-stroke, rack-and-pinion linear actuator that, typically, consists of a T-slide, several idlers, a transmission to drive the slide (powered by an electric motor) and a housing that holds the entire assembly. The actuator is driven by gear action on its top surface, and is guided and constrained by gear-bearing idlers on its other two parallel surfaces. The geometry, implemented with gear-bearing technology, is particularly effective. An electronic motor operating through a GBDET can directly drive the T slide against large loads, as a rack and pinion linear actuator, with no break and no danger of back driving. The actuator drives the slide into position and stops. The slide holes position with power off and no brake, regardless of load. With the T slide configuration, this GBDET has an entire T-gear surface on which to operate. The GB idlers coupling the other two T slide parallel surfaces to their housing counterpart surfaces provide constraints in five degrees-of-freedom and rolling friction in the direction of actuation. Multiple GB idlers provide roller bearing strength sufficient to support efficient, rolling friction movement, even in the presence of large, resisting forces. T-slide actuators can be controlled using the combination of an off-the-shelf, electric servomotor, a motor angle resolution sensor (typically an encoder or resolver), and microprocessor-based intelligent software. In applications requiring precision positioning, it may be necessary to add strain gauges to the T-slide housing. Existing sensory- interactive motion control art will work for T slides. For open-loop positioning, a stepping motor emulation technique can be used.

Vranish, John↗
NASA NTRS2001

The Role of Convection and Growth Competition in Phase Selection in Microgravity: Controlled Convection in the Containerless Processing of Steel Alloys

Containerless processing using electromagnetic levitation (EML) is a powerful technique in the investigation of reactive molten metal systems. On ground, the power required to overcome the weight of the sample is sufficient to cause significant heating and induce substantial melt convection. In microgravity, the heating and positioning fields may be decoupled and the field strength may be varied to achieve the desired level of convection within the limits set by the geometry of the levitation coil and the sample size. From high-speed digital images of the double recalescence behavior of Fe-Cr-Ni alloys in ground-based testing and in reduced-gravity aboard the NASA KC-135 parabolic aircraft, we have shown that phase selection can be predicted based on a growth competition model. An important parameter in this model is the delay time between primary nucleation and subsequent nucleation of the stable solid within the liquid/metastable solid array. This delay time is a strong function of composition and a weak function of the undercooling of the melt below the metastable liquidus. From the results obtained during the first Microgravity Sciences Laboratory (MSL-1) mission, we also know that convection may significantly influence the delay time, especially at low undercoolings. Currently, it is unclear what mechanism controls the formation of a heterogeneous site that allows nucleation of the austenitic phase on the pre-existing ferrite skeleton. By examining the behavior of the delay time under different convective conditions, we hypothesize that we can differentiate between several of these mechanisms to gain an understanding of how to control microstructural. evolution. We will anchor these predictions by examining samples quenched at different times following primary recalescence in microgravity. A second important parameter in the growth competition model is the identification of the growth rate of the stable phase into the semi-solid array that formed during primary recalescence. Current dendritic growth theory is inadequate in predicting solidification behavior under these conditions as metallographic analyses show that stable phase growth proceeds along the interface between the metastable solid and residual liquid. Since growth velocity is independent of the initial undercooling relative to the metastable liquidus, we hypothesize that purely thermal effects can be separated from other important growth model parameters by careful selection of the liquid composition in a ternary system.

Matson, D. M.↗
DOE OSTI2020

Enhancing nuclear data validation analysis by using machine learning

In this work, we demonstrate how machine learning (ML) techniques can be used as an auxiliary tool for nuclear data validation analysis. The output of the ML analysis can inform evaluators and validators on the quality, or lack thereof, of specific nuclear data and benchmark experiments with respect to simulating these experimental benchmark values. To this end, measured and simulated effective neutron multiplication factors, k eff , of 875 selected ICSBEP critical assemblies and the latter's sensitivities with respect to nuclear data as well as benchmarks' features (e.g., material nuclide, core geometry) are used as input for a random forest (RF) regression model. The RF is used to encode the complex inter-dependencies between thousands of nuclear data features (i.e., sensitivity profiles and aspects describing the measurements) and 875 simulated and experimental k eff values in order to understand what nuclear data features are most informative for predicting bias. The complexity of relationships and high-dimensional space is difficult-to-impossible to search fully using simply expert judgment. As a first proof-of-concept—step, it is demonstrated that this technique is able to correctly trace large discrepancies between simulated and experimental k eff back to fabricated shortcomings in nuclear data that were perturbed to simulated k eff values. In a second, real-case scenario, step, the RF algorithm is used to validate the ENDF/B-VIII.0 library in comparison to ENDF/B-VII.1 nuclear data. One case is showcased where the chosen ML algorithms highlighted nuclear data (the 19 F(n,inl) cross section from 0.4–0.9 MeV) that are shown to be problematic by comparing them to associated differential experimental data and nuclear data from other libraries. In addition to that, it is shown that the RF results point towards poor benchmark experiments and associated underestimated uncertainties (e.g., the PU-SOL-THERM-028 series). However, using the RF algorithm for validating nuclear data with respect to k eff is currently limited to pinpointing groups of questionable nuclear data due to the inherent correlations between features introduced by the nuclear data themselves and how k eff is simulated. Due to this, we recommend that the ML methods presented be used to augment—rather than replace—the expert knowledge of evaluators and validators.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗
DOE OSTI2026

Heat Transfer Experiments of a 1st Stage Blade Cascade for Supercritical CO2 Oxy-Combustion Turbine Application

The results of internally cooled 1st stage blade (S1B) cascade testing in a supercritical CO2 environment is presented. The turbine blade design has been previously established for the end application of an oxy-combustion turbine operating in the Allam-Fetvedt cycle with turbine inlet conditions of 305 bar and 1150°C. The internally cooled blade features leading edge (LE) region impingement cooling, mid-section ribbed serpentine passages, and a pin-finned trailing edge (TE) region before cooling ejection holes. The geometry for the tested blade cascade has a cooled central blade with un-cooled blades on either side to match flowpath areas of the actual turbine. The flowpath reuses internal components previously employed for mid-section region ribbed serpentine passage experiments that established Nusselt number enhancement ratios over a range of Reynolds numbers from 100,000-400,000. New components include flow conditioning plates upstream and downstream of the blade cascade to adequately represent the flow field and blade external heat transfer coefficient profiles for the actual turbine. The cooled central blade utilizes uniform crystal temperature sensors (UCTS) with six sensors each on the blade pressure and suction surfaces distributed radially and from LE to TE. The post-processed UCTS quantified the maximum wall temperature seen at each installed sensor location. The test procedure consisted of establishing supercritical CO2 cooling flow temperature and flow rate and maintaining it throughout the test. The flow rate aims to match that for the actual in-service turbine blade design and is maintained through an orifice restriction to keep the pressure differential between internal cooling flow and external hot flow nearly constant. For the sCO2 flow path external to the blade, temperatures were ramped throughout the test via control of the test loop’s natural gas burner heater. The maximum temperature seen was 468°C and held constant for a duration of 10 minutes at which the blade metal temperature was predicted to be at its maximum before ramping down. For the turbine blade design for service inlet conditions, external flow path computational fluid dynamics (CFD) results and an internal cooling 1-D thermal and hydraulic flow network model using experimentally validated correlations served as thermal finite element (FE) boundary conditions to predict blade metal temperatures. These predicted temperatures were subsequently utilized in a structural FE model to predict blade life ratings dictated by Haynes 282 creep strength data, having a strong dependence on temperature. The boundary conditions experienced during testing are used in the same workflow and compared to the experimental results, with the goal of validating the analysis methodology and providing insight on the uncertainty in local metal temperature predictions.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2025

Ligand‐Directed Actinide Oxo‐Bond Manipulation in Actinyl Thiacalix[4]arene Complexes

Understanding the chemistry of the inert actinide oxo bond in actinyl ions AnO2 2+ is important for controlling actinide behavior in the environment, during separations, and in nuclear waste (An=U, Np, Pu). The thioether calixarene TC4A (4-tert-butyltetrathiacalix[4]arene) binds equatorially to the actinyl cation forming a conical pocket that differentiates the two trans-oxo groups. The 'ate' complexes, [A]2[UO2(TC4A)] (A=[Li(DME)2], HNEt3) and [HNEt3]2[AnO2(TC4A)] (An=U, Np, Pu), enable selective oxo chemistry. Silylation of the UVI oxo groups by bis(trimethylsilyl)pyrazine occurs first at only the unencapsulated exo oxo and only one silylation is needed to enable migration of the endo oxo out of the cone, whereupon a second silylation affords the stable UIV cis-bis(siloxide) [A]2[U(OSiMe3)2(TC4A)]. Calculations confirm that only one silylation event is needed to initiate oxo rearrangement, and that the putative cis dioxo isomer of [UO2(TC4A)]2- would be stable if it could be accessed synthetically, at only 23 kcal.mol-1 in energy above the classical trans dioxo. Calculations for the transuranic cis[AnO2(TC4A)]2- (An=Np, Pu) are at higher energies, 30-35 kcal.mol-1, retaining the U complexes as the more obvious target for a cis-dioxo actinyl ion. The aryloxide (OAr) groups of the macrocycle are essential in stabilizing this as-yet unseen uranyl geometry as further bonding in the TC4A U-OAr groups stabilizes the U=O 'yl' bonds, explaining the stability of the putative cis[UO2(TC4A)]2- in this ligand framework.

Pyrch, Mikaela M↗
NASA NTRS1996

Comparative Study of Advanced Turbulence Models for Turbomachinery

A computational study has been undertaken to study the performance of advanced phenomenological turbulence models coded in a modular form to describe incompressible turbulent flow behavior in two dimensional/axisymmetric and three dimensional complex geometry. The models include a variety of two equation models (single and multi-scale k-epsilon models with different near wall treatments) and second moment algebraic and full Reynolds stress closure models. These models were systematically assessed to evaluate their performance in complex flows with rotation, curvature and separation. The models are coded as self contained modules that can be interfaced with a number of flow solvers. These modules are stand alone satellite programs that come with their own formulation, finite-volume discretization scheme, solver and boundary condition implementation. They will take as input (from any generic Navier-Stokes solver) the velocity field, grid (structured H-type grid) and computational domain specification (boundary conditions), and will deliver, depending on the model used, turbulent viscosity, or the components of the Reynolds stress tensor. There are separate 2D/axisymmetric and/or 3D decks for each module considered. The modules are tested using Rocketdyn's proprietary code REACT. The code utilizes an efficient solution procedure to solve Navier-Stokes equations in a non-orthogonal body-fitted coordinate system. The differential equations are discretized over a finite-volume grid using a non-staggered variable arrangement and an efficient solution procedure based on the SIMPLE algorithm for the velocity-pressure coupling is used. The modules developed have been interfaced and tested using finite-volume, pressure-correction CFD solvers which are widely used in the CFD community. Other solvers can also be used to test these modules since they are independently structured with their own discretization scheme and solver methodology. Many of these modules have been independently tested by Professor C.P. Chen and his group at the University of Alabama at Huntsville (UAH) by interfacing them with own flow solver (MAST).

Hadid, Ali H.↗
NASA NTRS1987

Damage growth in composite laminates with interleaves

The influence of placing interleaves between fiber reinforced plies in multilayered composite laminates is investigated. The geometry of the composite is idealized as two dimensional, isotropic, linearly elastic media made of a damaged layer bonded between two half planes and separated by thin interleaves of low extensional and shear moduli. The damage in the layer is taken in the form of a symmetric crack perpendicular to the interface and may extend up to the interface. The case of an H-shaped crack in the form of a broken layer with delamination along the interface is also analyzed. The interleaves are modeled as distributed shear and tension springs. Fourier integral transform techniques are used to develop solutions in terms of singular integral equations. An asymptotic analysis of the integral equations based on Muskhelishvili's techniques reveals logarithmically singular axial stresses in the half plane at the crack tips for the broken layer. For the H shaped crack, similar singularities are found to exist in the axial stresses at the interface crack tips in the layer and the half plane. The solution of the equations is found numerically for the stresses and displacements by using the Hadamard's concept of direct differentiation of Cauchy integrals as well as Gaussian integration techniques.

Goree, James G.↗
DOE OSTI2024

A double copy from twisted (co)homology at genus one

We study the twisted (co)homology of a family of genus-one integrals — the so called Riemann-Wirtinger integrals. These integrals are closely related to one-loop string amplitudes in chiral splitting where one leaves the loop-momentum, modulus and all but one puncture un-integrated. While not actual one-loop string integrals, they share many properties and are simple enough that the associated twisted (co)homologies have been completely characterized [1]. Using intersection numbers — an inner product on the vector space of allowed differential forms — we derive the Gauss-Manin connection for two bases of the twisted cohomology providing an independent check of [2]. We also use the intersection index — an inner product on the vector space of allowed contours — to derive a double-copy formula for the closed-string analogues of Riemann-Wirtinger integrals (one-dimensional integrals over the torus). Similar to the celebrated KLT formula between open- and closed-string tree-level amplitudes, these intersection indices form a genus-one KLT-like kernel defining bilinears in meromorphic Riemann-Wirtinger integrals that are equal to their complex counterparts.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Ligand‐Directed Actinide Oxo‐Bond Manipulation in Actinyl Thiacalix[4]arene Complexes

Abstract Understanding the chemistry of the inert actinide oxo bond in actinyl ions AnO 2 2+ is important for controlling actinide behavior in the environment, during separations, and in nuclear waste (An=U, Np, Pu). The thioether calixarene TC4A (4‐ tert ‐butyltetrathiacalix[4]arene) binds equatorially to the actinyl cation forming a conical pocket that differentiates the two trans‐ oxo groups. The ‘ate’ complexes, [A] 2 [UO 2 (TC4A)] (A=[Li(DME) 2 ], HNEt 3 ) and [HNEt 3 ] 2 [AnO 2 (TC4A)] (An=U, Np, Pu), enable selective oxo chemistry. Silylation of the U VI oxo groups by bis(trimethylsilyl)pyrazine occurs first at only the unencapsulated exo oxo and only one silylation is needed to enable migration of the endo oxo out of the cone, whereupon a second silylation affords the stable U IV cis‐ bis(siloxide) [A] 2 [U(OSiMe 3 ) 2 (TC4A)]. Calculations confirm that only one silylation event is needed to initiate oxo rearrangement, and that the putative cis dioxo isomer of [UO 2 (TC4A)] 2− would be stable if it could be accessed synthetically, at only 23 kcal.mol −1 in energy above the classical trans dioxo. Calculations for the transuranic cis [AnO 2 (TC4A)] 2− (An=Np, Pu) are at higher energies, 30–35 kcal.mol −1 , retaining the U complexes as the more obvious target for a cis ‐dioxo actinyl ion. The aryloxide (OAr) groups of the macrocycle are essential in stabilizing this as‐yet unseen uranyl geometry as further bonding in the TC4A U‐O Ar groups stabilizes the U=O ‘yl’ bonds, explaining the stability of the putative cis [UO 2 (TC4A)] 2− in this ligand framework.

Pyrch, Mikaela M. [College of Chemistry University↗
DOE OSTI2021

A multi-dimensional parametric study of variability in multi-phase flow dynamics during geologic CO 2 sequestration accelerated with machine learning

Successful geologic CO 2 storage projects depend on numerical simulations to predict reservoir performance during site selection, injection verification, and post-injection monitoring phases of the project. These numerical simulations solve non-linear sets of coupled partial differential equations, while accounting for multi-phase fluid dynamics on the basis of constitutive equations that are embedded into the solution scheme. As a consequence, individual simulations often require tens to hundreds of hours to complete on high-performance computing clusters. Moreover, laboratory experiments reveal that parametric functions for capillary pressure and relative permeability exhibit substantial variability, even within the same rock type. This combination of computational expense and wide-ranging parametric variability means that there remains substantial uncertainty in the behavior of multi-phase CO 2 -water systems, particularly in the context of feedbacks between relative permeability and capillary pressure. To bridge this knowledge gap, here we develop a novel workflow that utilizes physics-based numerical simulation to train an artificial neural network (ANN) emulator for interrogating the multivariate parameter space that governs both capillary pressure and relative permeability. With this approach, the ANN is trained to emulate both fluid pressure distribution and CO 2 saturation, which are then interrogated quantitatively to generate parametric response surface mappings with high-fidelity resolution. Results from this study initially show that capillary entry pressure is the dominant control on both CO 2 plume geometry and fluid pressure propagation when considering the combined effects of capillary pressure and relative permeability, particularly when phase interference is low and residual CO 2 saturation is high. Moreover, the ANN emulator provides tremendous computational speed-up by computing 2691 individual simulations in several minutes; whereas, the same simulation ensemble would have required ~3 years of simulation time using only physics-based simulation methods (25,000 times speed up).

58 GEOSCIENCES↗
DOE OSTI2024

Sodium-Cooled Fast Reactor Reference Plant Model

This report details the progress of Idaho National Laboratory (INL) in creating a reference plant multiphysics model for the Advanced Burner Test Reactor (ABTR). This model was developed under Task 13 of the U.S. Nuclear Regulatory Commission project “Development and Modeling Support for Advanced Non-Light Water Reactors,” and is an extension of the reference plant model developed in Task 4b, which was improved upon in the following ways. (1) The discrete ordinates method was used in lieu of the super-homogenization (SPH)-corrected diffusion approximation in order to better capture the anisotropic scattering contribution and the neutron leakage change due to thermal expansion. (2) The novel neutronic spatial discretization approach, termed the ring-heterogeneous (RH) approximation, was conceptualized and introduced to capture the differential expansion of the materials in the core. This new technique proved capable of preserving fission rates and maintaining the eigenvalue within 2.5% and 266 pcm with 9 neutron energy groups, respectively. Separating the different materials in the core enables the differential expansion of materials to be explicitly accounted for, eliminating the need for problem-specific cross-section functionalization techniques. (3) The SAM model for the core and system thermal-hydraulics analysis was updated to include 61 channels instead of just four representative ones. This enables users to obtain improved spatial resolution for sodium temperature and density scalar fields. (4) All the mesh files were created via the Multiphysics Object-Oriented Simulation Environment (MOOSE) Reactor module, eliminating all reliance on external tools for mesh creation. (5) Finally, the fuel axial expansion now leverages the HT9 and UPuZr material properties that have been validated against experimental data. The reference plant model was used to perform a full-core unprotected loss of flow (ULOF) transient calculation, including neutronics, thermal and mechanical feedback mechanisms. Future work will be devoted to further enhancements of the model. Potential improvements to the model include the addition of the control rod driveline expansion feedback and the upgrading of the support plate model so as to explicitly include 3D effects. Additionally, a Nuclear Energy Advanced Modeling and Simulation funded parallel effort has completely automated the creation of the ring-heterogeneous (RH) mesh from the fully heterogeneous (FH) geometry, thus maximizing user friendliness for the sodium fast reactor sodium-cooled fast reactor (SFR) workflow and will be incorporated in future work.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Effects of Internal and External Heat Sources on Cladding Microstructure and Rupture Performance

In the event of a Loss of Coolant Accident (LOCA), the primary supply of cooling water for a nuclear reactor is lost, leading to a significant pressure differential across the cladding wall. Without adequate cooling, the fuel rods continue to heat as a result of fission reactions. Research at Oak Ridge National Laboratory’s (ORNL’s) Severe Accident Test Station (SATS) is currently focused on evaluating fuel cladding performance using an external infrared lamp as a heat source, whereas legacy testing primarily utilized an internal heating approach. While external heating may better simulate the effect of neighboring fuel rods heating a central rod, it may not accurately represent the internal heat absorption from the fuel during accident transients. The heating dynamics depend greatly on the fuel rod's position within the bundle and the reactor. To fully understand the implications of a LOCA event and assess the influence of internal heating on cladding performance, combined internal heating and pressurization capability was developed at ORNL. Tests were conducted to compare cladding segments heated internally (representing heat from the fuel within the rod) and externally (representing heat from adjacent fuel rods). The findings indicate that both internal and external heating result in comparable rupture temperatures during 5°C/s laboratory LOCA tests and also agree with legacy test data. Axial temperature gradients and internal heat source dispersal were found to significantly impact cladding deformation and rupture geometry. Clear modifications were outlined to further improve the capability with heating rates above 5°C/s.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Effects of Internal and External Heat Sources on Cladding Microstructure and Rupture Performance

In the event of a Loss of Coolant Accident (LOCA), the primary supply of cooling water for a nuclear reactor is lost, leading to a significant pressure differential across the cladding wall. Without adequate cooling, the fuel rods continue to heat as a result of fission reactions. Research at Oak Ridge National Laboratory’s (ORNL’s) Severe Accident Test Station (SATS) is currently focused on evaluating fuel cladding performance using an external infrared lamp as a heat source, whereas legacy testing primarily utilized an internal heating approach. While external heating may better simulate the effect of neighboring fuel rods heating a central rod, it may not accurately represent the internal heat absorption from the fuel during accident transients. The heating dynamics depend greatly on the fuel rod's position within the bundle and the reactor. To fully understand the implications of a LOCA event and assess the influence of internal heating on cladding performance, combined internal heating and pressurization capability was developed at ORNL. Tests were conducted to compare cladding segments heated internally (representing heat from the fuel within the rod) and externally (representing heat from adjacent fuel rods). The findings indicate that both internal and external heating result in comparable rupture temperatures during 5°C/s laboratory LOCA tests and also agree with legacy test data. Axial temperature gradients and internal heat source dispersal were found to significantly impact cladding deformation and rupture geometry. Clear modifications were outlined to further improve the capability with heating rates above 5°C/s.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

The Schwarz Alternating Method for the Seamless Coupling of Nonlinear Reduced Order Models and Full Order Models

Projection-based model order reduction allows for the parsimonious representation of full order models (FOMs), typically obtained through the discretization of a set of partial differential equations (PDEs) using conventional techniques (e.g., finite element, finite volume, finite difference methods) where the discretization may contain a very large number of degrees of freedom. As a result of this more compact representation, the resulting projection-based reduced order models (ROMs) can achieve considerable computational speedups, which are especially useful in real-time or multi-query analyses. One known deficiency of projection-based ROMs is that they can suffer from a lack of robustness, stability and accuracy, especially in the predictive regime, which ultimately limits their useful application. Another research gap that has prevented the widespread adoption of ROMs within the modeling and simulation community is the lack of theoretical and algorithmic foundations necessary for the “plug-and-play” integration of these models into existing multi-scale and multi-physics frameworks. This paper describes a new methodology that has the potential to address both of the aforementioned deficiencies by coupling projection-based ROMs with each other as well as with conventional FOMs by means of the Schwarz alternating method [41]. Leveraging recent work that adapted the Schwarz alternating method to enable consistent and concurrent multiscale coupling of finite element FOMs in solid mechanics [35, 36], we present a new extension of the Schwarz framework that enables FOM-ROM and ROM-ROM coupling, following a domain decomposition of the physical geometry on which a PDE is posed. In order to maintain efficiency and achieve computation speed-ups, we employ hyper-reduction via the Energy-Conserving Sampling and Weighting (ECSW) approach [13]. We evaluate the proposed coupling approach in the reproductive as well as in the predictive regime on a canonical test case that involves the dynamic propagation of a traveling wave in a nonlinear hyper-elastic material.

97 MATHEMATICS AND COMPUTING↗