Search NASA⌕ Search

SEARCH · Search NASA

Results for “Numerical Algorithms”

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.

226 records · Page 9
All sourcesNASA NTRSDOE OSTI
NASA NTRS2026

Moon Base Transportation - Deliveries to the Lunar Surface

Development of the Moon Base will enable a home away from home for astronauts who will live and work at humanity’s first lunar outpost. In this effort, NASA’s Moon Transportation Office is responsible for enabling the transformational missions required to deliver habitats, supplies, science payloads, and all other elements needed to cultivate a permanent presence on the Lunar Surface. The Mission Concept (MC) is characterized through evaluation of an end-to-end architecture that can successfully deliver a generalized heavy large-volume payload, in excess of 4000 kg, to a precision landing and touchdown on the lunar surface. The mission architecture utilizes a single launch configuration of a Lunar Lander (LL) with a unique propellant system. The LL has an integral orbital transfer capability and features jettisonable elements. The design circumvents the need for prop transfer on orbit and multiple launch configurations. The launch vehicle (LV) for this work will assume the capability to deliver a payload in excess of 40,000 kg to orbit, affording multiple LV solutions. Considerations for the LL and payload deployment from the fairing are assumed to be handled through compliance with a launch providers’ Interface Requirements Document (IRD). The MC will span from launch at Kennedy Space Center (KSC) to terminal descent and touchdown on the lunar surface, requiring a total ΔV on the order of 6 km/s beyond what is required to get the vehicle stack to a 200 km circular Low Earth Orbit (LEO). Major mission phases include: launch and launch vehicle separation, transfer operations, pre-landing navigation, and lunar descent and touchdown. A Concept of Operations (ConOps) is used as the primary design driver for defining architecture of the vehicles necessary to achieve final payload delivery. Numerous ground rules and assumptions will be provided for each phase of the mission. Concept designs for the LL is presented. An emphasis of the design maximizes a feasible path for maturation, manufacturing, and operation. A self-imposed practical consideration for this effort is the incorporation of legacy designed hardware to minimize expensive, time-intensive, and high-risk hardware development cycles. The propulsion system of the LL adopts a conventional storable bipropellant configuration of monomethyl hydrazine (MMH) and mixed oxides of nitrogen (MON3). This effort will showcase a unique propellant delivery system to minimize the reliance on propellant management devices (PMDs) during descent. Numerous key constraints have been considered, across the multiple segments of the mission. These include the unique aspects of center of gravity (CG) management, thruster plume effects including self-impingement, propulsion system hardware limitations, navigation during multiple mission phases, and landing gear geometry for uneven terrain. These constraints shape the trades necessary for precision landing of heavy cargo and ensure compatibility with broader Moon Base Transportation concepts. The resulting insights inform future transportation strategies for the Moon and beyond; directly contributing to the development of cargo‑delivery standards that will support the long‑term buildup of a sustainable, continuously inhabited Moon Base.

Lunar Habitat↗
DOE OSTI2024

Impact of hydrogenation on the stability and mechanical properties of amorphous boron nitride

Abstract Interconnect materials with ultralow dielectric constant, and good thermal and mechanical properties are crucial for the further miniaturization of electronic devices. Recently, it has been demonstrated that ultrathin amorphous boron nitride (aBN) films have a very low dielectric constant, high density (above 2.1 g cm −3 ), high thermal stability, and mechanical properties. The excellent properties of aBN derive from the nature and degree of disorder, which can be controlled at fabrication, allowing tuning of the physical properties for desired applications. Here, we report an improvement in the stability and mechanical properties of aBN upon hydrogen doping. With the introduction of a Gaussian approximation potential for atomistic simulations, we investigate the changing morphology of aBN with varying H doping concentrations. We found that for 8 at% of H doping, the concentration ofsp 3 -hybridized atoms reaches to a maximum which leads to an improvement of thermal stability and mechanical properties by 20%. These results will be a guideline for experimentalists and process engineers to tune the growth conditions of aBN films for numerous applications.

Materials Science↗
DOE OSTI2024

Stellar Population Properties in the Stellar Streams around SPRC047

Abstract We have investigated the properties (e.g., age, metallicity) of the stellar populations of a ringlike tidal stellar stream (or streams) around the edge-on galaxy SPRC047 (z= 0.031) using spectral energy distribution (SED) fits to integrated broadband aperture flux densities. We used visual images in six different bands and Spitzer/IRAC 3.6μm data. We have attempted to derive best-fit stellar population parameters (metallicity, age) in three noncontiguous segments of the stream. Due to the very low surface brightness of the stream, we have performed a deconvolution with a Richardson–Lucy–type algorithm of the low spatial resolution 3.6μm IRAC image, thereby reducing the effect of the point-spread function aliasedemissionfrom the bright edge-on central galaxy at the locations of our three stream segments. Our SED fits that used several different star formation (SF) history priors, from an exponentially decaying SF burst to continuous SF, indicate that the age–metallicity–dust degeneracy is not resolved, most likely because of inadequate wavelength coverage and low signal-to-noise ratios of the low surface brightness features. We also discuss how future deep visual–near-infrared observations, combined with absolute flux calibration uncertainties at or below the 1% level, complemented by equally well absolute flux-calibrated observations in ultraviolet and mid-infrared bands, would improve the accuracy of broadband SED fitting results for low surface brightness targets, such as stellar streams around nearby galaxies that are not resolved into stars.

Astronomy & Astrophysics↗
DOE OSTI2025

Terahertz conductivity of two-dimensional materials: a review

Two-dimensional (2D) van der Waals materials are shaping the landscape of next-generation devices, offering significant technological value thanks to their unique, tunable, and layer-dependent electronic and optoelectronic properties. Time-domain spectroscopic techniques at terahertz (THz) frequencies offer noninvasive, contact-free methods for characterizing the dynamics of carriers in 2D materials. They also pave the path toward the applications of 2D materials in detection, imaging, manufacturing, and communication within the increasingly important THz frequency range. In this paper, we overview the synthesis of 2D materials and the prominent THz spectroscopy techniques: THz time-domain spectroscopy, optical-pump THz-probe technique, and optical pump–probe THz spectroscopy. Through a confluence of experimental findings, numerical simulation, and theoretical analysis, we present the current understanding of the rich ultrafast physics of technologically significant 2D materials: graphene, transition metal dichalcogenides, MXenes, perovskites, topological 2D materials, and 2D heterostructures. Finally, we offer a perspective on the role of THz characterization in guiding future research and in the quest for ideal 2D materials for new applications.

2D materials↗
DOE OSTI2026

Time Correlations from Steady-State Expectation Values

Recovering properties of correlation functions is typically challenging. On the one hand, experimentally, it requires measurements with a temporal resolution finer than the system’s dynamics. On the other hand, analytical or numerical analysis requires solving the system evolution. Here, we use recent results of quantum metrology with continuous measurements to derive general lower bounds on the relaxation and second-order correlation times that are both easy to calculate and measure. These bounds are based solely on steady-state expectation values and their derivatives with respect to a system parameter, and can be readily extended to the autocorrelation of arbitrary observables. We validate our method on two examples of critical quantum systems: a critical driven-dissipative resonator, where the bound matches analytical results for the dynamics, and the infinite-range Ising model, where only the steady state is solvable, and thus the bound provides information beyond the reach of existing analytical approaches. Our results can be applied to the experimental characterization of ultrafast systems and to the theoretical analysis of many-body models whose dynamics are hard to compute.

Górecki, Wojciech [INFN, Pavia] (ORCID:00000001991↗
NASA NTRS2025

Improving Adhesive Bondline Time of Flight Predictions During Autoclave Cure Utilizing Machine Learning

Composite materials are increasingly being used in aerospace applications due to their superior strength-to-weight ratio compared to commonly used metals. A current limitation to widespread adoption is the certification of adhesively bonded joints. One approach to improving adhesive bonding in composites is accurately measuring the thickness of adhesive bondlines in composite laminates. Precise bondline thickness control is essential for aerospace applications where adhesive layer thickness directly affects joint fracture properties and structural performance. This study focused on implementing machine learning techniques to determine the ultrasonic time of flight (directly correlated to thickness) in adhesive bondlines throughout autoclave cure cycles. A high-temperature (use up to 180°C) ultrasonic scanning system was deployed in an autoclave to provide time of flight data through composite panels. Three experiments were conducted on the curing of 305 mm × 305 mm unidirectional composite panels. In the first experiment, a piecewise function was fit for the temperature correction factor to account for changing autoclave temperatures. Due to deficiencies in the first calibration experiment, a second experiment was run, and the results were used to train a machine learning model. The revised experiment, in combination with the machine learning model, significantly increased the accuracy of the bondline time of flight predictions (~14% error reduced to <1%). Data was processed using the Regression Learner Application in MATLAB®, with a Support Vector Machine selected for the model. The result was a machine learning algorithm capable of reliably quantifying ultrasonic time of flight through adhesive bondlines. The third experiment provided independent test data for the machine learning model, demonstrating that the model produces accurate predictions from data beyond its training set.

Machine Learning↗
DOE OSTI2024

User-Centric Communication With Aerial Network for 6G: A Reinforcement Learning Approach

Meeting the diverse needs of user verticals requires innovative cellular architectures that can offer additional degrees of freedom to provide on-demand services. The terrestrial user-centric radio access network (UC-RAN) stands out as an excellent choice for this purpose. However, a drawback of UC-RAN is its tendency to prioritize high-priority verticals, often resulting in a subpar quality of experience for low-priority verticals. This issue is particularly exacerbated in hotspot areas. Here, to address this problem, we introduce an aerial network integrated with terrestrial UC-RAN to provide coverage to users which are not served by the terrestrial network. Furthermore, we analyze the impact of key configuration and optimization parameters (COPs), such as location, transmit power, altitude, and beamwidth of aerial base stations (ABSs) on system key performance indicators (KPIs), such as coverage, latency satisfaction, average spectral efficiency, and energy efficiency. We formulate a robust multiobjective function to maximize these KPIs without biasing toward any specific KPI(s). Finally, we propose a deep reinforcement learning optimization framework based on the state-of-the-art soft actor-critic algorithm to control ABS COPs and optimize system KPIs. Experimental evaluations demonstrate that the proposed optimization framework can converge to near-optimal solutions derived from the pseudo brute force in a few thousand epochs.

6G↗
NASA NTRSDate not supplied

Preliminary Findings of the Experimental Development Unit Cold Flow Test for a Generation Zero Nuclear Propulsion Engine

Nuclear thermal propulsion (NTP) technology will greatly benefit human travel to Mars by significantly shortening transit times, improving crew safety, and providing more mission flexibility than traditional chemical rockets. As part of DRACO follow-on work to develop, build, and fly a generation zero NTP engine, a full scale flight-like experimental design unit (EDU) reactor was constructed to collect sufficient on-ground performance data to characterize flow induced vibrations (FIV) of critical reactor structures/components, inform development of the engine and reactor control algorithm, and collect pressure drop and flow distribution data across the reactor. The fluid conditions for the test program were designed to achieve system responses equivalent to that of an operational engine through all phases of engine operation including reactor startup, mainstage operation, reactor shutdown, and reactor cooldown. Over 100 tests were executed, flowing either GN2 or GHe through the EDU at varying flow rates and pressures. This experiment provided early validation of flow behavior and vibration risks before nuclear testing, boosted critical subsystem TRLs, informed design iterations, and reduced future test costs. The steady-state flow parameters for the experiment were modeled in Ansys Thermal Desktop, allowing rapid tuning and experiment-informed updates to a flight-like test matrix. The EDU dynamic environment was characterized with accelerometers, strain gauges, and high-frequency pressure transducers all sampled at 20 kHz. While many narrow-band oscillations were identified, no significant FIV occurred; the reactor structural responses tend to be enveloped by typical launch vehicle ascent vibration environments (defined up to 2 kHz), although significant energy is also present at higher frequencies.

Flow Induced Vibration↗
DOE OSTI2025

Polaron catastrophe within quantum acoustics

The quantum acoustic framework has recently emerged as a nonperturbative, coherent approach to electron–lattice interactions, uncovering rich physics often obscured by perturbative methods with incoherent scattering events. Here, we model the strongly coupled dynamics of electrons and acoustic lattice vibrations within this framework, representing lattice vibrations as coherent states and electrons as quantum wave packets, in a manner distinctively different from tight-binding or discrete hopping-based approaches. We derive and numerically implement electron backaction on the lattice, providing both visual and quantitative insights into electron wave packet evolution and the formation of acoustic polarons. We investigate polaron binding energies across varying material parameters and compute key observables—including mean square displacement, kinetic energy, potential energy, and vibrational energy—over time. Our findings reveal the conditions that favor polaron formation, which is enhanced by low temperatures, high deformation potential constants, slow sound velocities, and high effective masses. Additionally, we explore the impact of external electric and magnetic fields, showing that while polaron formation remains robust under moderate fields, it is weakly suppressed at higher field strengths. These results deepen our understanding of polaron dynamics and pave the way for future studies into nontrivial transport behavior in quantum materials.

Science & Technology - Other Topics↗
NASA NTRS2026

Impact Testing of 410 Stainless Steel for Material Impact Model Development

A project is underway to develop a consistent set of material properties, impact test data, and failure analysis for a variety of metallic aircraft materials that can be used to develop improved impact failure and deformation models. This project is jointly funded by the NASA Glenn Research Center and the Federal Aviation Administration William J. Hughes Technical Center. Particular features of this set of data are that all material property and impact test data are obtained using traceable material, the test methods and procedures are extensively documented, and all the raw data are available. Four parallel efforts are currently underway: measurement of material deformation and failure response over a wide range of strain rates and temperatures, failure analysis of material property specimens and impact test articles, development of improved numerical modeling techniques for deformation and failure, and impact testing of flat panels and substructures for model validation. This report describes impact testing performed on 410 stainless steel sheet and plate samples of different thicknesses with two different types of projectiles, one a regular cylinder and one with a more complex geometry incorporating features representative of a jet engine fan blade. Data from this testing will be used in validating material models developed under this program. The material tests and the material models developed in this program will be published in separate reports.

410 Stainless Steel↗
DOE OSTI2026

Sparse non-Markovian Noise Modeling of Transmon-Based Multi-Qubit Operations

The influence of noise on quantum dynamics is one of the main factors preventing current quantum processors from performing accurate quantum computations. Sufficient noise characterization and modeling can provide key insights into the effect of noise on quantum algorithms and inform the design of targeted error protection protocols. However, constructing effective noise models that are sparse in model parameters, yet predictive can be challenging. In this work, we present an approach for effective noise modeling of multi-qubit operations on transmon-based devices. Through a comprehensive characterization of seven devices offered by the IBM Quantum Platform, we show that the model can capture and predict a wide range of single- and two-qubit behaviors, including non-Markovian effects resulting from spatiotemporally correlated noise sources. The model’s predictive power is further highlighted through multi-qubit dynamical decoupling demonstrations and an implementation of the variational quantum eigensolver. As a training proxy for the hardware, we show that the model can predict expectation values within a relative error of 0.5%; this is a sevenfold improvement over default hardware noise models. Through these demonstrations, we highlight key error sources in superconducting qubits and illustrate the utility of reduced noise models for predicting hardware dynamics.

open quantum systems & decoherence↗
NASA NTRS2025

Erosion Results of the MISSE 9-15 Polymers and Composites Experiment 1-4 (PCE 1-4)

Polymers and other oxidizable materials on the exterior of spacecraft in the low Earth orbit (LEO) space environment can be eroded from reaction with atomic oxygen (AO). Therefore, in order to design durable spacecraft it is important to know the extent of erosion that will occur during a mission. This can be determined by knowing the LEO AO erosion yield, E y (volume loss per incident oxygen atom), of materials susceptible to AO reaction. In addition, recent flight experiments have shown that the AO E y can vary with the AO fluence and/or solar exposure. Therefore, obtaining AO E y data for materials flown on various spaceflight missions is important. NASA Glenn Research Center has flown numerous experiments as part of the Materials International Space Station Experiment (MISSE) missions on the exterior of the International Space Station (ISS) to characterize the LEO E y of polymers, composites, protective coatings, and other spacecraft materials. Recently, four Glenn experiments with 365 flight (F) samples were flown on ISS’s MISSE-Flight Facility (MISSE-FF). These experiments are the Polymers and Composites Experiment-1 (PCE-1) flown as part of the MISSE-9 mission, the PCE-2 flown as part of the MISSE-10 mission, the PCE-3 flown as part of the MISSE-12 and MISSE-15 missions, and the PCE-4 flown as part of the MISSE-13 mission. Although each experiment had numerous sample objectives, the primary objective was to determine the LEO AO E y of various spacecraft materials as a function of solar irradiation and AO fluence. This paper provides a summary of the erosion data for the PCE 1-4 AO E y samples. The AO E y for 150 samples flown in either the LEO ram, wake, zenith or nadir directions are provided. The AO ram fluence varied from 2.97×10 20 atoms/cm 2 after 0.89 years of direct space exposure (with relatively high levels of Si contamination) on MISSE 12 to 3.93×10 20 atoms/cm 2 after 1.17 years of direct space exposure on MISSE-10. The ram AO E y values for uncoated polymers range from 3.81×10 –25 cm 3 /atom for polytetrafluoroethylene (M9R-C20 F) exposed to an AO fluence of 3.44×10 20 atoms/cm 2 on MISSE-9 to 4.43×10 –23 cm 3 /atom for AO etched low density polyimide aerogel (M12R-C21 F) exposed to an AO fluence of 2.97×10 21 atoms/cm 2 on MISSE-12. Because of the low AO fluence and relatively high Si contamination, a number of PCE-3 wake samples experienced mass gain. Thus, AO E y values are not provided for these samples. Although there are calculated AO E y values for the zenith, wake and nadir samples, the ram AO E y for a particular material is a more reliable value in terms of AO exposure because the zenith, wake and nadir directions were exposed to either no or very little AO fluence and thus other space environmental factors (i.e. vacuum, thermal extremes and thermal cycling, and/or various types of radiation) are responsible for the mass loss.

Erosion yield↗
NASA NTRS2006

2005 NASA Seal/Secondary Air System Workshop

The 2005 NASA Seal/Secondary Air System workshop covered the following topics: (i) Overview of NASA’s new Exploration Initiative program aimed at exploring the Moon, Mars, and beyond; (ii) Overview of the NASA-sponsored Propulsion 21 Project; (iii) Overview of NASA Glenn’s seal project aimed at developing advanced seals for NASA’s turbomachinery, space, and reentry vehicle needs; (iv) Reviews of NASA prime contractor, vendor, and university advanced sealing concepts including tip clearance control, test results, experimental facilities, and numerical predictions; and (v) Reviews of material development programs relevant to advanced seals development. Turbine engine studies have shown that reducing high-pressure turbine (HPT) blade tip clearances will reduce fuel burn, lower emissions, retain exhaust gas temperature margin, and increase range. Several organizations presented development efforts aimed at developing faster clearance control systems and associated technology to meet future engine needs. The workshop also covered several programs NASA is funding to develop technologies for the Exploration Initiative and advanced reusable space vehicle technologies. NASA plans on developing an advanced docking and berthing system that would permit any vehicle to dock to any on-orbit station or vehicle. Seal technical challenges (including space environments, temperature variation, and seal-on-seal operation) as well as plans to develop the necessary “androgynous” seal technologies were reviewed. Researchers also reviewed tests completed for the shuttle main landing gear door seals.

Hypersonic↗
NASA NTRS2026

PERSIANN-Unet: A Global Deep Learning Framework for Near-Real-Time Precipitation Estimation Using Infrared Data

Access to high-quality, high-resolution, near-real-time precipitation data is essential for hydrological and meteorological research and disaster mitigation. Traditional tools such as rain gauges and radar networks, though effective, have limitations, including sparse coverage in remote areas and high operational costs. Satellite data, with its global coverage and high spatial and temporal resolutions, mitigates limitations in coverage. Satellite precipitation products like Hydro Estimator (HE), Integrated Multi-satellitE Retrievals for Global Precipitation Measurement (IMERG), and Precipitation Estimation from Remotely Sensed Information using Artificial Neural Networks (PERSIANN) utilize both geosynchronous thermal infrared (IR) and passive microwave (PMW) data in their operation. PMW sensors offer detailed atmospheric profiles but suffer from higher latency, whereas IR sensors provide lower latency but only capture cloud-top information. Despite this constraint, IR data remains attractive for low-latency precipitation estimation. Recent advances in deep learning, particularly convolutional neural networks (CNNs), have further improved satellite precipitation retrievals. This study introduces PERSIANN-Unet (PUnet or PERSIANN V3), a quasi-global algorithm covering 60°N–60°S that combines IR data, monthly climatology, and the UNet architecture to produce half-hourly precipitation estimates at 0.04° resolution. The product is evaluated against HE, IMERG, and PDIR-Now for 2022–2023. Results show that PUnet closely matches its training target, IMERG V07 Final, at the global scale, and performance is further evaluated against Stage IV as a reference over CONUS. Training PUnet on IMERG (2016–2021) leverages a high-quality, integrated PMW IR-gauge precipitation product while developing an IR-based framework not reliant on PMW availability. By operating on a single global image, PUnet avoids tile partitioning and blending steps, reducing edge discontinuities, and produces more spatially consistent precipitation fields across hemispheres.

Phu Nguyen↗
DOE OSTI2026

Active learning using hybrid surrogate tool life modeling for machining process optimization

Here, this paper describes an active learning approach for part-to-part iterative machining process optimization using a hybrid surrogate tool life model. A probabilistic interpolating tool life model is developed by combining the empirical Taylor-type tool life equation and the model fit error. The probabilistic tool life model is then used to calculate the machining cost per part distribution. The optimal machining parameters are selected using an expected improvement in machining cost per part criterion. The method is validated numerically using experimental results; the results show a median convergence error of 2.2% after three tests over 400 simulations. The method is validated experimentally on two industrial applications for Ti-6Al-4V roughing resulting in a cost per part reduction greater than 23% after two tests. The described method is a robust solution for rapid convergence to optimal machining parameters in an industrial production environment.

Active learning↗
DOE OSTI2026

ClassNMSW- a real-time classification approach for non-recycled municipal solid waste using hyperspectral imaging

Real-time classification of non-recycled municipal solid waste (NMSW) is essential for efficient valorization. This study introduces ClassNMSW, a comprehensive framework for classifying 22 NMSW subclasses under industrial constraints by using hyperspectral imaging (HSI). A primary innovation of this work is the development of a variance-controlled spectral extraction algorithm. Unlike traditional methods that rely on simple averaging, this approach systematically investigates the extent of pixel extraction to minimize the loss of critical chemical information while maximizing data reduction thus ensuring high spectral fidelity with low computational cost. The approach developed in this work integrates automated, computer-vision-based background removal, eliminating the need for the manual thresholding common in current literature. To resolve ambiguities among chemically similar subclasses, a tiered classification and multi-camera fusion strategy (NIR17 and NIR22) is implemented. Results demonstrate that ClassNMSW achieves an object-wise weighted accuracy of 98.70% for single-sensor configurations and 100% under sensor fusion. A novel rolling-window strategy satisfies desired end-to-end latency of <2 s, satisfying the strict deterministic requirements of high-speed industrial sorting environments. The ClassNMSW framework provides a scalable foundation for advancing circularity and resource recovery in large-scale waste valorization operations.

99 - GENERAL AND MISCELLANEOUS↗
DOE OSTI2024

Asymmetric fluctuations and self-folding of active interfaces

We study the structure and dynamics of the interface separating a passive fluid from a microtubule-based active fluid. Turbulent-like active flows power giant interfacial fluctuations, which exhibit pronounced asymmetry between regions of positive and negative curvature. Experiments, numerical simulations, and theoretical arguments reveal how the interface breaks up the spatial symmetry of the fundamental bend instability to generate local vortical flows that lead to asymmetric interface fluctuations. The magnitude of interface deformations increases with activity: In the high activity limit, the interface self-folds invaginating passive droplets and generating a foam-like phase, where active fluid is perforated with passive droplets. These results demonstrate how active stresses control the structure, dynamics, and break-up of soft, deformable, and reconfigurable liquid–liquid interfaces.

active fluid↗
DOE OSTI2026

SAM Finite Volume Method Development Status Update: GCR Application, Restart, and MultiApp

The System Analysis Module (SAM) is being developed as a modern system analysis code for advanced non-light-water-reactor safety analysis under the U.S. DOE NEAMS program. Previous feasibility studies have demonstrated that a staggered-grid finite volume method (SG-FVM), implemented under the MOOSE framework, can deliver more than an order of magnitude speedup over the existing continuous Galerkin finite element method (CG-FEM) solver for liquid-cooled, incompressible but thermally expandable flow systems. This work extends the previous effort to compressible, gas-cooled reactor applications, where pressure couples directly into the mass equation adding additional nonlinearity into the equation system. New code capabilities are implemented for pebble bed high-temperature gas-cooled reactor (PB-HTGR) analysis, including a pebble bed CoreChannel component, built-in pebble bed effective thermal conductivity model and channel-to-channel crossflow model. The capabilities are tested, benchmarked, and demonstrated for problems with increased level of model and physical complexities, including the HTTU effective thermal conductivity test, the SANA passive cooling test, and a demonstration case using the GPBR200 reactor design covering steady-state operation, DLOFC and PLOFC transients. Across all cases, the SG-FVM solver demonstrated strong robustness and efficiency, and the solutions agree well with reference results and data. The finding of this work proves that SG-FVM is a viable and efficient solver pathway for compressible, gas-cooled reactor system analysis in SAM. In addition, work has been done to successfully support SAM-FVM recover/restart code feature that is essential to reactor safety analysis applications, and MultiApp code feature that is essential to multi-scale and multi-physics simulations. In summary, this work continued from previous feasibility studies, and further demonstrated that the SG-FVM will serve as a strong foundation for SAM’s advanced solver algorithm for future deployment.

Zou, Ling↗