Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “mesh optimization”

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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Deploy threading in Nalu solver stack

The goal of the ExaWind project is to enable predictive simulations of wind farms composed of many MW-scale turbines situated in complex terrain. Predictive simulations will require computational fluid dynamics (CFD) simulations for which the mesh resolves the geometry of the turbines, and captures the rotation and large deflections of blades. Whereas such simulations for a single turbine are arguably petascale class, multi-turbine wind farm simulations will require exascale-class resources. The primary code in the ExaWind project is Nalu, which is an unstructured-grid solver for the acousticallyincompressible Navier-Stokes equations, and mass continuity is maintained through pressure projection. The model consists of the mass-continuity Poisson-type equation for pressure and a momentum equation for the velocity. For such modeling approaches, simulation times are dominated by linear-system setup and solution for the continuity and momentum systems. For the ExaWind challenge problem, the moving meshes greatly affect overall solver costs as re-initialization of matrices and re-computation of preconditioners is required at every time step In this Milestone, we examine the effect of threading on the solver stack performance against flat-MPI results obtained from previous milestones using Haswell performance data full-turbine simulations. Whereas the momentum equations are solved only with the Trilinos solvers, we investigate two algebraic-multigrid preconditioners for the continuity equations: Trilinos/Muelu and HYPRE/BoomerAMG. These two packages embody smoothed-aggregation and classical Ruge-Stiiben AMG methods, respectively. In our FY18 Q2 report, we described our efforts to improve setup and solve of the continuity equations under flat-MPI parallelism. While significant improvement was demonstrated in the solve phase, setup times remained larger than expected. Starting with the optimized settings described in the Q2 report, we explore here simulation performance where OpenMP threading is employed in the solver stack. For Trilinos, threading is acheived through the Kokkos abstraction where, whereas HYPRE/BoomerAMG employs straight OpenMP. We examined results for our mid-resolution baseline turbine simulation configuration (229M DOF). Simulations on 2048 Haswell cores explored the effect of decreasing the number of MPI ranks while increasing the number of threads. Both HYPRE and Trilinos exhibited similar overal solution times, and both showed dramatic increases in simulation time in the shift from MPI ranks to OpenMP threads. This increase is attributed to the large amount of work per MPI rank starting at the single-thread configuration. Decreasing MPI ranks, while increasing threads, may be increasing simulation time due to thread synchronization and start-up overhead contributing to the latency and serial time in the model. These result showed that an MPI+OpenMP parallel decomposition will be more effective as the amount per MPI rank computation per MPI rank decreases and the communication latency increases. This idea was demonstrated in a strong scaling study of our low-resolution baseline model (29M DOF) with the Trilinos-HYPRE configuration. While MPI-only results showed scaling improvement out to about 1536 cores, engaging threading carried scaling improvements out to 4128 cores — roughly 7000 DOF per core. This is an important result as improved strong scaling is needed for simulations to be executed over sufficiently long simulated durations (i.e., for many timesteps). In addition to threading work described above, the team examined solver-performance improvements by exploring communication-overhead in the HYPRE-GMRES implementation through a communicationoptimal- GMRE algorithm (CO-GMRES), and offloading compute-intensive solver actions to GPUs. To those ends, a HYPRE mini-app was allow us to easily test different solver approaches and HYPRE parameter settings without running the entire Nalu code. With GPU acceleration on the Summitdev supercomputer, a 20x speedup was achieved for the overall preconditioner and solver execution time for the mini-app. A study on Haswell processors showed that CO-GMRES provides benefits as one increases MPI ranks.

17 WIND ENERGY↗
NASA NTRSDate not supplied

TPSAS-NF1676L-33512-DND

INTRODUCTION: In an effort to quantify the predictive capability of the Livermore Software Technology Corporation (LSTC)-provided Hybrid III (small female, mid-size male, and large male) anthropomorphic test device (ATD) finite element models (FEMs) within aerospace loading environments, the National Aeronautics and Space Administration (NASA) performed a series of sled tests with the isolated head-neck complex of each ATD. These tests and subsequent modeling efforts revealed significant deficiencies in the large male head neck FEM. Differences in both the shape and material properties between the FEM and the physical ATD were identified. In order to improve the predictive response of the FEM, a model development effort was undertaken to update the head-neck complex to better represent the geometric and material properties of the physical ATD.METHODS: To better match the physical geometry, a new neck component model was meshed and implemented into the Hybrid III large male head-neck complex. The updated geometry and weight were verified against measurements of the physical ATD to insure accuracy. The primary source of material discrepancy between the physical ATD and FEM was found to be in the material definition of the rubber neck pucks. Using LS-Opt the material properties of this part were calibrated to match physical response of the head-neck complex in the isolated head-neck testing performed. RESULTS AND DISCUSSION: The updated geometry and calibrated material properties were shown to significantly improve the FEM predictive accuracy under lateral, rear, and frontal impacts with combined horizontal and vertical loading. The results of this study demonstrate an effective means for improving ATD FEM response through isolated component level testing and material parameter optimization. The improved and quantified accuracy of the Hybrid III large male FEM head-neck complex stemming from with work lends to confidence in its future use for occupant protection evaluation of large male occupants

J Putnam↗
DOE OSTI2022

Heating, DPA, and He production in the Proton Beam Window Assembly for the Second Target Station

This report details the analysis of the heating rates, dpa values, and He production in the Proton Beam Window (PBW), the assembly around the PBW, and the shielding above it. Neutronics calculations use a detailed Unstructured Mesh (UM) converted from the solid CAD model of the PBW design of June 2022. The rest of the geometry is taken from the CSG master model (as of June 2022). For this analysis, a proton beam profile with a footprint of ≈90 cm 2 has been used. This proton profile was created by a Ring to Second Target magnet array configuration using octupole magnets. Maximum values in the PBW (aluminum) are 5.4 J/cc/pulse heating rate, 1.1 dpa per year, and 660 appm He-production per year. These results within ~ 10% with the previous calculations from June 2020. Average heating and dpa rates have been calculated for each component in the PBW assembly and will be used to evaluate and design the cooling needs and make decisions on removable components for the PBW. It is desirable to optimize and simplify the PBW design by only replacing components that are at the end of their lifetime. Additional analysis for the energy deposition in the PBW is performed for two super-Gaussian beam profiles with a footprint of 60 cm 2 and 30 cm 2 respectively. Also, an estimate of the PBW lifetime is performed for two realistic quadrupole beam profiles.

43 PARTICLE ACCELERATORS↗
DOE OSTI2026

Actinium–DOTA coordination in water from hybrid ML/MM: Structure, free energies, and water-exchange pathways

Quantitative simulation of trivalent ƒ-block chelates in water remains challenging because bonded and non-bonded force-field models make different approximations for coordination structure, exchange dynamics, and ion–ligand interactions in highly charged systems. Here, we develop a hybrid machine-learning/molecular-mechanics (ML/MM) framework for Ac 3+ –DOTA in explicit solvent by training an E(3)-equivariant neural network potential (MACELES) on mechanically embedded QM/MM data for Ac aquo and Ac–DOTA species and coupling it to NAMD 2.14 with particle-mesh Ewald electrostatics. Nanosecond ML/MM trajectories remain numerically stable and preserve chelate integrity, yielding a compact DOTA inner shell with an inner-sphere water coordination number of CN Ac,O w ≈ 1.7 arising from a dynamic equilibrium between one- and two-water states (37.5% and 59.9% of frames; three waters 2.5%). A 5 ns potential of mean force shows two low-lying basins at CN Ac,O w ≈ 1 and CN Ac,O w ≈ 2. DFT end-state free energies are consistent with the ML/MM profile, and DFT minimum-energy paths provide a qualitative electronic-structure reference for the observed basin connectivity. State-resolved kinetics reveal picosecond water-exchange pathways that couple hydration changes to transient DOTA arm fluctuations, and training-set comparisons show that temperature-matched Ac–DOTA data optimize energy/force accuracy while more diverse solvated data improve charge prediction. Overall, the present hybrid ML/MM model provides a practical description of Ac 3+ –DOTA hydration thermodynamics and short-time exchange behavior in explicit water at MD-like cost.

Actinium↗
DOE OSTI2021

Mathematical modeling of novel porous transport layer architectures for proton exchange membrane electrolysis cells

Thin foil based porous transport layers (PTLs) that contain highly structured pore arrays have shown promise as anode PTLs in proton exchange membrane electrolysis cells. These novel PTLs, fabricated with advanced manufacturing techniques, produce thin, tunable, multifunctional layers with reduced flow and interfacial resistances and high thermal and electric conductivities. To further optimize their design, it is important to understand their fundamental impact on the transport of protons, electrons, and liquid/vapor mixtures in the electrode. In this work, we develop a two-dimensional multiphysics model to simulate the coupled electrochemistry and multiphase transport in an electrolysis cell operated with the novel PTL architecture. The results show that larger pores improve access of water to the anode catalyst layer, which is beneficial for both the oxygen evolution reaction and membrane hydration. Larger pore sizes also improve oxygen gas transport from the catalyst layer, because generated oxygen gas is forced to travel in-plane through the anode catalyst layer until it reaches a pore opening that is connected to a channel. The discussed results confirm that the proposed thin foil based PTLs are fundamentally different from conventional PTLs, such as felts or layered meshes. The model developed in this work also provides generalizable insight into fundamental PEMEC phenomena, such as the competition between liquid and gas phase transport, membrane hydration and water management, and nonuniform electrochemical reactions, which are processes relevant to all PEMEC designs.

25 ENERGY STORAGE↗
DOE OSTI2023

A subdivision-stabilized B-spline mixed material point method

Subjected to external loadings, polymeric materials, e.g., biological tissues, hydrogels, and elastomers, may undergo extreme, nearly incompressible, (self-)contact deformations. For numerical modeling employing mesh-based techniques such as the finite element method (FEM), these deformations pose significant challenges due to large distortions in the deformed geometry, accuracy issues stemming from volumetric locking effects, and increased computational cost from complex contact searches. As an alternative to mesh-based methods, the material point method (MPM), a continuum-based particle technique, is gaining attention for its ability to handle extreme distortions and capture no-slip contact without added cost. For nearly incompressible material behaviors, while mixed formulations can address locking effects by treating displacements and pressure as independent fields, they can suffer from numerical instabilities close to the incompressibility limit due to the violation of the inf-sup condition, leading to inaccurate nodal pressure solutions. Here we propose an efficient and stable mixed B-spline material point method with highest achievable regularity for quasi-compressible polymeric materials. Using the two-scale relation of B-splines, we introduce a subdivision-stabilization for the two-field mixed MPM and obtain numerically stable, oscillation-free nodal solutions with equal-order interpolations with optimal regularity. Building on the Eulerian-Lagrangian nature of MPM, a previously-converged solution framework is adopted to mitigate issues related to cell-crossing and numerical fracture artifact present in standard MPM. We assess the stability and accuracy of the developed mixed MPM at large deformations for soft materials through the benchmark Cook’s membrane problem. Additionally, we test the robustness of the proposed MPM by modeling several examples, including the compression and indentation of a circular block into a quasi-compressible substrate and the twisting deformation of a rectangular block. The findings demonstrate the MPM’s capabilities for modeling practical soft material applications.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Adaptive mesh based combustion simulations of direct fuel injection effects in a supersonic cavity flame-holder

Here, we present high-fidelity reacting simulations of a supersonic cavity flame-holder configuration. The focus of this work is on flame stabilization brought about by varying the location of fuel injection in a cavity stabilized supersonic flow of air. Central to our approach is a compressible multi-species reacting flow solver that uses adaptive-mesh-refinement (AMR), enabling the resolution of flame, shock-waves, boundary-layers, and small-scale structures in the computational domain. Our analysis indicates that fuel injection closer to the ramp at the aft end of the cavity allows for greater mixing and lower peak temperatures compared to fuel injection upstream that is closer to the backward facing step of the cavity. This difference is mainly due to greater turbulent fluctuations generated from the shear-layer towards the cavity ramp, thereby enhancing the mixing of fuel and air. A low frequency oscillatory behaviour in heat-release and pressure was also observed for the upstream injection case while a much higher-frequency phenomena was observed in the near-ramp injection case. By identifying the important physical determinants of the combustion processes, this study illustrates a promising pathway to design and optimize direct fuel injection strategies in supersonic cavity flame-holders that can improve flame stability, combustion efficiency, and reduce emissions.

42 ENGINEERING↗
DOE OSTI2019

MULTIDIMENSIONAL NUMERICAL SIMULATIONS OF REACTING FLOW IN A NON-PREMIXED ROTATING DETONATION ENGINE

Over the last two decades, detonation based propulsion has received a great deal of attention as a potential means to achieve significant improvement in the performance of air-breathing and rocket engines. Detonative combustion mode is particularly interesting due to the resulting pressure gain from reactants to products, faster heat release, decreased entropy generation, more available work and higher thrust compared to conventional deflagrative combustion. Rotating detonation engine (RDE) is one such novel combustor concept. Realistic RDE configurations utilize separate fuel and air injection schemes, hence are not perfectly premixed. Moreover, RDE performance is governed by a large number of design parameters and operating conditions. In this context, computational fluid dynamics (CFD) has the potential to enhance the understanding of RDE combustion and aid future development/optimization of this technology. In the present work, a CFD model was developed to simulate a representative non-premixed RDE combustor. Unsteady Reynolds-Averaged Navier-Stokes (RANS) simulations were performed for the full combustor geometry (including the separate fuel and air injection ports), with hydrogen as fuel and air as the oxidizer. Adaptive mesh refinement (AMR) was incorporated to achieve a trade-off between model accuracy and computational expense. A finite-rate chemistry model along with a 10-species detailed kinetic mechanism was employed to describe the H2-Air combustion chemistry. Two operating conditions were simulated, corresponding to the same global equivalence ratio of unity but different fuel and air mass flow rates. For both conditions, the capability of the model to capture the essential detonation wave dynamics was assessed. A validation study was performed against experimental data available on detonation wave frequency/height, reactant fill height, oblique shock angle, axial pressure distribution in the channel, and fuel/air plenum pressure. The CFD model predicted the sensitivity of these wave characteristics to the operating conditions with good accuracy, both qualitatively and quantitatively. The present CFD model offers a potential capability to perform rapid design space exploration and/or performance optimization studies for realistic full-scale RDE configurations.

Pal, Pinaki↗
DOE OSTI2018

Efficient modeling of optically-complex, non-coplanar exterior shading: Validation of matrix algebraic methods

It has long been established that shading windows with overhangs, fins, and other types of non-coplanar systems (NCS) is one of the most effective ways of controlling solar heat gains in buildings because they intercept solar radiation prior to entry into the building. Designers however often specify non-opaque materials (e.g., louvers, fritted glass, expanded metal mesh) for these systems in order to admit daylight, reduce lighting energy use, and improve indoor environmental quality. Most simulation tools rely on geometric calculations and radiosity methods to model the solar heat gain impacts of NCS and cannot model optically-complex materials or geometries. For daylighting analysis, optically-complex NCS can be modeled using matrix algebraic methods, although time-efficient parametric analysis has not yet been implemented. Determining the best design and/or material for static or operable NCS that minimize cooling, heating, and lighting energy use and peak demand requires an iterative process. This study describes and validates a matrix algebraic method that enables parametric energy analysis of NCS. Such capabilities would be useful not only for design but also for development of prescriptive energy-efficiency standards, rating and labeling systems for commercial products, development of design guidelines, and development of more optimal NCS technologies. A facade or “F” matrix, which maps the transfer of flux from the NCS to the surface of the window, is introduced and its use is explained. A field study was conducted in a full-scale outdoor testbed to measure the daylight performance of an operable drop-arm awning. Simulated data were compared to measured data in order to validate the models. Results demonstrated model accuracy: simulated workplane illuminance was within 11–13%, surface luminance was within 16–18%, and the daylight glare probability was within 6–9% of measured results. Methods used to achieve accurate results are discussed. Results of the validation of daylighting performance are applicable to solar heat gain performance. Since exterior shading can also significantly reduce peak demand, these models enable stakeholders to more accurately assess HVAC and lighting impacts in support of grid management and resiliency goals.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

Differentiable Multiphysics Codes: A Breakthrough Technology for Simulation and Computing

This document summarizes the findings of a strategic planning exercise commissioned by the Weapons Simulation and Computing, Computational Physics (WSC/CP) program at the Lawrence Livermore National Laboratory (LLNL) in FY24. During the year, the committee met with multiple stakeholder communities to gather input, opinions, suggestions and concerns which have been incorporated throughout this document. The key findings from this exercise are summarized: • The development of multiphysics modelling and simulation (mod/sim) codes and software technologies, their deployment on exascale compute platforms, and their broad adoption across the NNSA is a major success of the Advanced Simulation and Computing (ASC) program and the Exascale Computing Project (ECP). Sustained investment in these core technologies is essential. • Today’s state of the art involves running ensembles of O(100K) simulations to perform uncertainty quantification (UQ) and design studies using multiple statistical methods such as Bayesian optimization to understand sensitivities of our models and explore parameterized design spaces. Even with exascale computing, we are practically limited to O(10) parameters in these studies since the number of simulations required to sample the space scales exponentially with the number of design parameters. • The data from these simulation ensembles is increasingly being used to train machine learned (ML) surrogates (or reduced order models, ROMs) which can then be used for optimization or real time design exploration. However, the trained surrogates are still limited in the number of parameters they can represent due to the sampling limitations previously noted. • Augmenting our suite of integrated multiphysics simulation codes, both current and emerging, with the ability to compute gradients (solution derivatives) of arbitrary simulation outputs with respect to (some or all) simulation inputs would be a breakthrough technology, opening the door to a new era of efficient and automated inverse design based on verified and validated mod/sim capabilities. • This capability, which we refer to as differentiable multiphysics codes (DMCs), would revolutionize both UQ and optimization studies by breaking the curse of dimensionality that presently limits our “gradient-free” ensemble based computing approach. A similar breakthrough occurred in the AI/ML community once the ability to compute gradients of arbitrary loss functions using back-propagation became commonplace. Gradient information from the multiphysics codes can also be used to dramatically improve the efficiency and scale of training of ML/ROM surrogates for rapid assessments. • Achieving this in our suite of codes will be a grand challenge, similar to the amount of effort that was required to transition from CPU to GPU computing. It will require buy-in from the entire WSC/CP program and beyond, including all integrated codes, physics and engineering models, third-party library dependencies and performance portability abstractions. It will also require investment in research and development of numerical methods for computing adjoints of coupled physics across multiple adaptively refined moving meshes and of stochastic (Monte Carlo) and mesh free (SPH) methods. • New software and numerical techniques, largely pioneered by the AI/ML community, make this feasible. Chief among these is automatic differentiation (AD), the ability to employ AD at point-wise locations in a physics calculation (instead of traditional black-box approaches) and the ability to perform “back-propagation in time” (or reverse mode AD) for non-linear partial differential equations (PDEs). Fundamentally, the conclusion of this strategic planning exercise is that the time is right to undertake a large scale effort in WSC, centered on the existing integrated codes, to continue the natural evolution of mod/sim in the age of AI/ML. Instead of attempting to replace mod/sim with purely data driven AI/ML models, we believe the key to success is to integrate AI/ML by building on top of the decades of hard-won knowledge and the verified/validated multiphysics modelling capability that is the hallmark of the ASC program.

97 MATHEMATICS AND COMPUTING↗
NASA NTRSDate not supplied

Toward Verification of the γ-Reθt Transition Model in OVERFLOW and FUN3D

The results of an ongoing assessment of the transition modeling capability in NASA's OVERFLOW and FUN3D programs are presented, with a focus on the Langtry-Menter γ-Reθt transition model in combination with the shear-stress transport (SST) turbulence model. While the effect of numerics and boundary conditions on the accuracy and iterative convergence of RANS solutions for fully turbulent flows has been well documented, especially in the context of canonical flow configurations, the same cannot be said for transport-equation-based transition models coupled with RANS-based turbulence models. Given the criticality of transition modeling for new aircraft design and optimization, there has been a renewed focus on the accuracy of such transition models and their inconsistent implementation across different flow solvers as seen from the AIAA and NATO-AVT workshops. In this work, we aim to establish the verification of such models via these two well established CFD codes with different numerics. The goal is to produce high-quality data, such as grids, solutions, and other auxiliary data, that may be utilized for code verification by others in the CFD community. As a first step, the work reported here is focused on the Langtry-Menter transition model as applied to two simple 2D configurations, namely the flat plate, and the NLF-0416, respectively. The paper also highlights how the boundary conditions and baseline turbulence model can affect the solutions from the SST-based Langtry-Menter transition model. A preliminary evaluation of the automatic mesh adaption capabilities of these solvers and potential benefits for flow configurations involving a mix of laminar, transitional, and fully turbulent flows is also reported.

CFD Modeling↗
NASA NTRS1993

Towards a generalized computational fluid dynamics technique for all Mach numbers

Currently there exists no single unified approach for efficiently and accurately solving computational fluid dynamics (CFD) problems across the Mach number regime, from truly low speed incompressible flows to hypersonic speeds. There are several CFD codes that have evolved into sophisticated prediction tools with a wide variety of features including multiblock capabilities, generalized chemistry and thermodynamics models among other features. However, as these codes evolve, the demand placed on the end user also increases simply because of the myriad of features that are incorporated into these codes. In order for a user to be able to solve a wide range of problems, several codes may be needed requiring the user to be familiar with the intricacies of each code and their rather complicated input files. Moreover, the cost of training users and maintaining several codes becomes prohibitive. The objective of the current work is to extend the compressible, characteristic-based, thermochemical nonequilibrium Navier-Stokes code GASP to very low speed flows and simultaneously improve convergence at all speeds. Before this work began, the practical speed range of GASP was Mach numbers on the order of 0.1 and higher. In addition, a number of new techniques have been developed for more accurate physical and numerical modeling. The primary focus has been on the development of optimal preconditioning techniques for the Euler and the Navier-Stokes equations with general finite-rate chemistry models and both equilibrium and nonequilibrium thermodynamics models. We began with the work of Van Leer, Lee, and Roe for inviscid, one-dimensional perfect gases and extended their approach to include three-dimensional reacting flows. The basic steps required to accomplish this task were a transformation to stream-aligned coordinates, the formulation of the preconditioning matrix, incorporation into both explicit and implicit temporal integration schemes, and modification of the numerical flux formulae. In addition, we improved the convergence rate of the implicit time integration schemes in GASP through the use of inner iteration strategies and the use of the GMRES (General Minimized Resisual) which belongs to the class of algorithms referred to as Krylov subspace iteration. Finally, we significantly improved the practical utility of GASP through the addition of mesh sequencing, a technique in which computations begin on a coarse grid and get interpolated onto successively finer grids. The fluid dynamic problems of interest to the propulsion community involve complex flow physics spanning different velocity regimes and possibly involving chemical reactions. This class of problems results in widely disparate time scales causing numerical stiffness. Even in the absence of chemical reactions, eigenvalue stiffness manifests itself at transonic and very low speed flows which can be quantified by the large condition number of the system and evidenced by slow convergence rates. This results in the need for thorough numerical analysis and subsequent implementation of sophisticated numerical techniques for these difficult yet practical problems. As a result of this work, we have been able to extend the range of applicability of compressible codes to very low speed inviscid flows (M = .001) and reacting flows.

Walters, R. W.↗
DOE OSTI2026

Integrated Neutronics Modeling for Inertial Fusion Energy Systems: Development and Application to LD-FIRST

Lawrence Livermore National Laboratory (LLNL) is proposing a new Laser Driven Fusion Integration Research and Science Test Facility (LD-FIRST) with the goal of providing an experimental testbed for future Inertial Fusion Energy (IFE) systems. However, IFE systems require detailed and accurate multiphysics modeling to quantify material damage, thermal loading, and tritium breeding within complex chamber environments. This article presents the first step in an integrated multiphysics framework that couples meshed CAD-based geometry within Monte Carlo neutronic simulations to enable high-fidelity analysis of IFE chamber concepts, with future coupling to external codes. The neutronics workflow utilizes OpenMC and its third-party capability to use CAD-based geometries through DAGMC and tally on unstructured meshes with Libmesh to evaluate neutron transport behavior, geometric fidelity, and material performance under reactor-relevant conditions. The use of tailored tallies on unstructured meshes in this framework allows direct transfer without interpolating to CFD simulation tools. Two IFE chambers were evaluated, both conceived by LLNL: HYLIFE-II and Laser IFE (LIFE). This work produced high-fidelity conformal surface and volumetric meshes of the HYLIFE-II and LIFE chambers with mapped spatial insight into material damage, thermal loading, and tritium breeding. The HYLIFE-II model was built utilizing available resources and used as a test case to verify that the neutronics framework can handle complex geometries. The LIFE chamber CAD was provided by LLNL and was the main focus of this work. This work analyzes multiple ternary alloy breeding materials for the LIFE chamber, across different 6 Li enrichments to produce data relevant to the LD-FIRST project. This work also investigates the level of model fidelity for the LIFE chamber, and results show that inclusion of detailed first wall and coolant structures increased the predicted tritium breeding ratio (TBR) by ~30%, highlighting the sensitivity of tritium breeding and the need for a high-fidelity simulation framework for IFE chambers. These developments provide a scalable toolset for the design and optimization of next-generation IFE chambers, forming a solid foundation for future coupled multiphysics analysis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2019

Cognitive Communications for NASA Space Systems

The growing complexity of spacecraft constellations, communication relay offerings, and mission architectures drives the need for the development of autonomous communication systems. NASA has traditionally launched single spacecraft missions that are served by the Space Communication and Navigation (SCaN) program. Operations on SCaN networks are typically scheduled weeks in advance, and often each asset serves a single user spacecraft at a time. Recent movement towards swarm missions could make the current approach unsustainable. Additionally, the integration of commercial communication service providers will substantially increase the data transfer options available to new missions. NASA science missions have found benefit in launching swarms of spacecraft, allowing coordinated simultaneous observations from different perspectives. Inter-spacecraft communication (mesh networking) is an enabler for this architecture, as are CubeSats that allow cost-effective provisioning of distributed mission assets. As more complex swarm missions launch, one challenge is coordinating communication within the swarm and choosing the appropriate mechanism for telemetry, tracking, control, and data services to and from Earth. Cognitive communications research conducted by SCaN aims to mitigate the increasing communication complexity for mission users by increasing the autonomy of links, networks, and service scheduling. By considering automation techniques including recent advances in artificial intelligence and machine learning, cognitive algorithms and related approaches enable increased mission science return, improved resource utilization for service provider networks, and resiliency in unpredictable or unplanned environments. The Cognitive Communications Project at the NASA Glenn Research Center develops applications of data-driven, non-deterministic methods to improve the autonomy of space communication. The project emphasizes development of decentralized space networks with artificial intelligence agents optimizing communication link throughput, data routing, and system-wide asset management. This paper discusses the objectives, approaches, and opportunities of the research to address growing needs of the space communications community.

Chelmins, David↗
NASA NTRS1999

Acoustic Treatment Design Scaling Methods: Test Plans, Hardware, Results, and Evaluation - Volume 3

The ability to design, build, and test miniaturized acoustic treatment panels on scale-model fan rigs representative of the full-scale engine provides not only a cost-savings, but an opportunity to optimize the treatment by allowing tests of different designs. To be able to use scale model treatment as a full-scale design tool, it is necessary that the designer be able to reliably translate the scale model design and performance to an equivalent full-scale design. The primary objective of the study presented in this volume of the final report was to conduct laboratory tests to evaluate liner acoustic properties and validate advanced treatment impedance models. These laboratory tests include DC flow resistance measurements, normal incidence impedance measurements, DC flow and impedance measurements in the presence of grazing flow, and in-duct liner attenuation as well as modal measurements. Test panels were fabricated at three different scale factors (i.e., full-scale, half-scale, and one-fifth scale) to support laboratory acoustic testing. The panel configurations include single-degree-of-freedom (SDOF) perforated sandwich panels, SDOF linear (wire mesh) liners, and double-degree-of-freedom (DDOF) linear acoustic panels.

Yu, J.↗
DOE OSTI2026

Reduced-Order CFD Modeling to Support Waste Loading Optimization in Hanford WTP Vitrification

The U.S. DOE Hanford Site stores over 56 million gallons of radioactive liquid tank waste that must be treated and immobilized for long-term disposal The Waste Treatment and Immobilization Plant (WTP) will vitrify this waste by feeding it into Joule-heated melters, where it is incorporated into a stable borosilicate glass Computational fluid dynamics (CFD) simulations of glass melters can provide insight into the maximum achievable waste loading under varying melter operating conditions Fully resolved VOF multiphase simulations were used as the reference model to capture bubble-driven convection in the melter, including bubble formation, rise behavior, and induced glass melt circulation Effective bubble column diameter and rise velocity were extracted from the resolved simulations, compared with empirical correlations, and refit across relevant viscosity and gas flow rate conditions Explicit gas–liquid interface tracking was replaced with a single-phase momentum source term model, enabling faster steady-state CFD simulations while preserving the dominant hydrodynamic effects of bubbling New empirical correlations were developed for effective bubble column diameter and bubble rise velocity by fitting resolved simulation data across expected melter viscosity and gas flow rate ranges, providing improved inputs for the momentum source term model compared with existing literature correlations The momentum source term model reduced fluid-domain mesh size by 89% and achieved an 8.4× computational speedup relative to resolved bubbling simulations The validated momentum source term approach enables prediction of process-relevant heat transfer behavior in the integrated melter model, including heat transfer from the molten glass to the cold cap, plenum, refractory walls, and surrounding structural regions under varying melter operating conditions

12 - MGMT OF RADIOACTIVE AND NON-RADIOACTIVE WASTE↗
DOE OSTI2017

A novel vibration-based fault diagnostic algorithm for gearboxes under speed fluctuations without rotational speed measurement

The localized failures of gears introduce cyclic-transient impulses in the measured gearbox vibration signals. These impulses are usually identified from the sidebands around gear-mesh harmonics through the spectral analysis of cyclo-stationary signals. However, in practice, several high-powered applications of gearboxes like wind turbines are intrinsically characterized by nonstationary processes that blur the measured vibration spectra of a gearbox and deteriorate the efficacy of spectral diagnostic methods. Although order-tracking techniques have been proposed to improve the performance of spectral diagnosis for nonstationary signals measured in such applications, the required hardware for the measurement of rotational speed of these machines is often unavailable in industrial settings. Moreover, existing tacho-less order-tracking approaches are usually limited by the high time-frequency resolution requirement, which is a prerequisite for the precise estimation of the instantaneous frequency. To address such issues, a novel fault-signature enhancement algorithm is proposed that can alleviate the spectral smearing without the need of rotational speed measurement. This proposed tacho-less diagnostic technique resamples the measured acceleration signal of the gearbox based on the optimal warping path evaluated from the fast dynamic time-warping algorithm, which aligns a filtered shaft rotational harmonic signal with respect to a reference signal assuming a constant shaft rotational speed estimated from the approximation of operational speed. The effectiveness of this method is validated using both simulated signals from a fixed-axis gear pair under nonstationary conditions and experimental measurements from a 750-kW planetary wind turbine gearbox on a dynamometer test rig. Lastly, the results demonstrate that the proposed algorithm can identify fault information from typical gearbox vibration measurements carried out in a resource-constrained industrial environment.

17 WIND ENERGY↗
DOE OSTI2022

An OpenMP GPU-offload implementation of a non-equilibrium solidification cellular automata model for additive manufacturing

Here, in this paper, performance strategies on GPU-based HPC platforms of a cellular automata (CA) simulation code for non-equilibrium solidification, including nucleation, grain growth, solute partitioning and transport for the metal additive manufacturing (AM) process are investigated using OpenMP 4.5. To accurately report the speed-up for multicore CPUs and GPUs, a rigorous performance analysis employed optimizations appropriate for both CPU-only code (baseline) and GPU offload codes for an isothermal test problem. The performance results on Summit at the Oak Ridge Leadership Computing Facility indicate that using a precomputed list of interface cells significantly decreased the wall-clock time on GPUs. The speedup due to GPU acceleration was evaluated for a full Summit node and measured to be 1.8X when comparing a 6 MPI tasks run with 6 GPUs versus 36 MPI tasks on the CPU only. That speed-up was found to be 7.9X when comparing 6 MPI tasks with 6 GPUs versus the 6 MPI tasks running on the CPU only. Performance measurements showed that system total time is almost constant for runs with more than 96 MPI tasks (or GPUs), indicating that the GPU-accelerated code showed an excellent weak scaling performance. Finally, a rapid directional solidification problem was considered to demonstrate the CA code capability on Summit. It was found that a mesh size of at least 0.05 μm is recommended for the AM-like simulations in order to obtain accurate elongated grain microstructure and elongated subgrain features, which are in qualitative good agreement with experimental data. The results presented in this study indicate that the performance strategies on GPU-based HPC platforms for the CA code are appropriate for novel HPC exascale platforms.

36 MATERIALS SCIENCE↗