Search NASA⌕ Search

SEARCH · Search NASA

Results for “formulation”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Development of Property Composition Models For RPP-WTP LAW Glasses (Final Report)

This report describes the development of property-composition models for low-activity waste (LAW) glasses for the River Protection Project Waste Treatment Plant (RPP-WTP) at the Hanford site. The RPP-WTP will separate Hanford tank wastes into LAW and high-level waste (HLW) streams and each stream will be vitrified separately. The development of LAW and HLW glass formulations for the RPP-WTP has been reported previously and is an ongoing activity. Acceptable formulations must meet a variety of processability, product quality, and waste loading requirements that are dictated either by the RPP-WTP contract or by the characteristics of the particular treatment processes that have been selected. These requirements amount to constraints on the acceptable ranges of certain glass properties. These properties are determined first and foremost by the composition of the glass or glass melt. Thus, while there is no direct way of controlling the glass properties of interest during production, there are simple and extremely effective methods of achieving the same result by instead controlling the glass composition. This basic principle is no different from that used to produce enormous volumes of commercial glass to meet exacting product specifications. An essential difference in waste vitrification, however, is that one of the raw materials (the waste itself) can be subject to considerable compositional uncertainty and variability. Thus, waste vitrification facilities and associated process control systems (of which, the operating envelope in glass formulation space is a key part) must, of course, be designed to be robust with respect to such variations. The determination of quantitative relationships between the glass properties that must be controlled and the glass composition can play an important role in the development of such facilities.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

DuraMelter 100 Sub-Envelope Changeover Testing Using LAW SubEnvelopes A1 and C1 Feeds in Support of the LAW Pilot Melter (Final Report)

The primary goal of the testing described in this report was to develop and recommend a compliant HLW glass formulation to support the actual waste testing of AZ-101 Envelope D waste (blended with actual pretreatment products including Cs- and Tc-eluates from pretreatment of AP-101 and AZ-101 LAW). Testing of actual waste will be performed at Battelle, Pacific Northwest Division. The test objective was met by the development and recommendation of the glass formulation HLW98-95; the formulation has been transmitted to the WTP to support vitrification of HLW AZ-101 actual waste.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2025

Multi-scale, Multi-disciplinary, and Multi-agent Explainable AI with Koopman-Undergirded Learning, Prediction, and Analysis (M3EA KULPA) (Project Closeout Report)

The goal of this project was to develop and use domain-aware machine learning formulations, based on the Koopman Operator (KO), for modelling multi-scale, multi-disciplinary (e.g., multi-physics), and/or multi-agent systems. The project developed these formulations for the following cases: • Systems with dynamics at two separate time scales, • Systems with a bi-level hierarchical control structure, • Systems with bi-level hierarchical control and dynamics at two separate time scales (the lower level controls operating at the faster time scale), and • Systems with n separate but interacting agents/disciplines (with/without control, respectively); the controls for each agent could include bi-level hierarchical control and dynamics at two separate time scales as described above. The project then defined a set of dynamical systems consisting of different nonlinear oscillators that could be used to test these different formulations and then subsequently learned the KO models for those systems. With the KO models, we were able to do the following: • Quantify system stability, including both long-term and transient behavior, • Quantify the effects of feedbacks between the different time scales and agents/disciplines in terms of those feedbacks’ effects on system stability, • Replace a standard Proportional-Integral (PI) control in the hierarchical control structure with a KO-based Linear-Quadratic Regular (LQR), a form of optimal control, • Calculate optimal supervisory control policies a) with and without time scale separated dynamics at the lower level control levels and b) with both PI and KO-based LQR lower level control policies, and • Calculate dynamic Nash equilibria for multi-agent systems where each agent makes its own control decisions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Foam and Unreacted Material Reduction in Magnesium Oxysulfate for SRPPF Waste Solidification

A magnesium oxysulfate (MOS) grout formulation has been proposed to solidify the liquid effluent from the Savannah River Plutonium Processing Facility (SRPPF) Aqueous Recovery System (ARS). This formulation uses a combination of light-burnt MgO, anhydrous MgSO 4 , and dead-burnt MgO resulting in a grout with good mixability and acceptable density. The setting time for the grout is within the operational limits of the SRPPF facility and the leachate pH is ~9.4, which is within the assumed pH range of the brine from the Waste Isolation Pilot Plant (WIPP). When first tested at the 1-gallon small-scale and 55-gallon full-scale, there was a significant foam layer that raised concerns of a nonhomogeneous final form. A nonhomogeneous mixture could have variable density throughout the waste form and therefore allow for unequal shielding. To reduce the potential for unequal shielding in the final form, a series of foam reduction tests were carried out and a final formulation of 28 wt% light burnt MgO, 10 wt% MgSO 4 , and 62 wt% dead burnt MgO was used to remove the foam in under two hours while still maintaining a mix that has an acceptable density, mix time, peak temperature, and set time.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2025

Water Sorption Study of Various Polymers and Polymeric Composites

In the formulations of many composite materials, polymeric materials are often used as binder materials to enhance processability, functionality, and safety during manufacturing. For instance, poly(ester urethane), such as Estane® 5703, is used in the formulation of PBX 9501, and vinyl copolymer elastomer (VCE) is used in the formulation of VCE/filler compos ites. The stability of these binders directly impacts the overall stability and performance of their composites. Estane and VCE are prone to hydrolytic degradation, with water sorption behavior playing a significant role in their long-term stability and aging behavior. Hence, it is critical to study their water sorption behavior under various humidity and tempera ture conditions. To improve the accuracy and reduce labor-intensive work involved in water sorption experiments, an automated Dynamic Vapor Sorption (DVS) system was employed, which can generate large datasets. To efficiently process and analyze these datasets, Python scripts were developed. These scripts not only streamlined data processing, but also accu rately derived the Arrhenius parameters related to water diffusion and sorption properties for the tested materials. This approach offers a robust framework for evaluating water trans port and sorption processes in polymeric materials and adds valuable capabilities to future water sorption testing efforts.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Effects of SRPPF ARS Effluent Variations on Magnesium Oxysulfate Solidification

A magnesium oxysulfate (MOS) grout formulation has been identified to solidify the liquid effluent from the Savannah River Plutonium Processing Facility (SRPPF) Aqueous Recovery System (ARS). The formulation uses reactive, light burnt MgO (28 wt%), anhydrous MgSO 4 (10 wt%), and dead burnt MgO (62 wt%), resulting in a grout with good mixability, acceptable density, acceptable setting times, and a leachate pH around 9.4, within the assumed Waste Isolation Pilot Plant (WIPP) brine pH range. A series of tests were then designed to observe how this formula reacted to variations in the liquid effluent. Liquid effluent temperature, the caustic pH level, concentration of NaNO 3 , and the presence of neutralization products, trace metals, and/or sodium sulfate (all potential residuals found in the effluent from upstream processing), were evaluated to see how each impacts the grout mixing time, peak temperature, setting time, presence of bleed water, leachate pH, and density. The standard mix has an average mixing time of 18 minutes, peak temperature of 100.5 °C, and set time of 60 minutes. Increasing the temperature of the liquid effluent increased the reaction rate of the mix and reduced mixing and setting times. Between 30 and 40°C the change was relatively small, raising no more than 6 °C compared to the standard mix maximum temperature. Increasing the pH of the liquid effluent had no significant impact on the mix. The addition of neutralization products and trace metals had an overall impact of decreasing the reactivity of the mix. The addition of more ions in the liquid effluent, such as an increase in the concentration of NaNO 3 and the addition of Na 2 SO 4 , generally decreased the reactivity of the mix. None of the variations to the liquid effluent hindered solidification as indicated by the lack of bleed water found on all the samples. The leachate pH values for all mixes tested did not significantly vary from the expected pH of 9.4 and none of the density values fell below 1.8 g/cm 3 . Overall, changes to the liquid effluent were found to have a minimal impact on the formulation suggesting the grout's ability to properly form despite increased temperatures and the presence of residuals typically found in the liquid effluent from the ARS.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2026

Melter Tests to Define LAW Halide Concentrations, Phase 2

The present work represents the second of two phases to address the formation of secondary sulfate phases while processing high waste loading LAW formulations with elevated halide, chromium, and phosphate feed concentrations over a wide range of LAW waste types, including those with high potassium concentrations. The Test Plan for the present work outlines a series of tests designed to assess the formation of secondary phases over a range of sulfur and halide concentrations for two LAW streams and associated glass formulations as well as crucible testing to identify additives that have the potential to suppress the formation of secondary phases for wastes with high chlorine contents. The combined results from both phases of testing and from previous tests with these LAW waste streams will be used as inputs to future work to develop concentration limits for these components in the LAW glass property composition models and to support the modification and updating of the current WTP LAW glass formulation algorithm.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2026

Quantum Computing in Next-Generation Transportation Optimization

We explore how quantum computing (QC) can advance transportation optimization, with a focus on two high-impact areas: traffic signal control and vehicle electrification with grid integration. As transportation systems grow in complexity, classical optimization methods increasingly struggle to deliver scalable and efficient solutions, particularly for real-time, data-rich environments. This work identifies key challenges within these two domains where QC may offer advantages, particularly in handling combinatorial decision spaces and dynamic constraints. We begin by outlining the limitations of classical approaches for traffic signal control optimization and electric vehicle charging coordination, highlighting where computational limitations arise. Previous quantum formulations are presented and new formulations are proposed to demonstrate how emerging quantum algorithms, including quantum annealing and the Quantum Approximation Optimization Algorithm, could be leveraged to reformulate and address these problems. We also evaluate the suitability of current quantum hardware and discuss recent trends that indicate when QC may become a viable tool for transportation applications. While acknowledging the present limitations of QC technologies, this poster emphasizes the importance of preparing quantum-compatible models today. By reviewing and establishing formulations that align with the strengths of quantum algorithms, researchers and practitioners can better position themselves to take advantage of QC advancements as they occur. This work aims to provide a practical, forward-looking perspective on the near-term potential of quantum computing in transportation optimization.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2025

Direct Feed High-Level Waste APPS Model Glass Testing (DFHLW APPS) Matrix, Phase 2

This report summarizes the data collected during the batching and melting of a second matrix of Direct Feed High-Level Waste (DFHLW) glasses generated using the preliminary enhanced waste glass models (EWG2.5) and the Britton and Anderson (2024) preliminary DFHLW feed vector. The purpose of these glasses is two-fold: 1. Validate EWG2.5 glass calculations being used in the Aspen Process Performance Simulation (APPS) model. 2. Evaluate and ultimately improve the glass property models and formulation methods used for design of DFHLW glasses as part of an iterative process of data collection and model refinement. Some of the 16 APPS2 glasses tested did not satisfy all target property constraints due to the limited data on DFHLW glass supporting the EWG2.5 models. • One glass, APPS2-10, formed nepheline on canister centerline cooling (CCC) heat-treatment and failed the product consistency test (PCT) response limits. This glass also had high B and Cr release rates for the toxicity characteristic leaching procedure (TCLP). All other glasses were found to satisfy the PCT and TCLP constraints for both quenched and CCC samples. • One glass, APPS2-08, had higher than acceptable viscosity due to magnetite crystallization. • One glass, APPS2-09, formed greater than 2 vol% crystals at 950 °C. As the glass design criterion was that the temperature at 2 vol% crystal (T 2% ) be less than 950 °C, only one glass failed the criteria. However, this criterion is being reevaluated. Four additional glasses formed crystal fractions between 1 and 2 vol% at 950 °C (APPS2-03, -08, -12, and -14). • Four glasses – APPS2-01, -02, -04, and -16 – failed the Monofrax K-3 refractory neck corrosion (k neck ) design limit of 0.04 in. at 1208 °C for 6 d. This is another criterion being reevaluated. Four additional glasses (APPS2-05, -06, -11, and -13) exhibited 0.025 = k neck = 0.04 in. • All 16 glasses passed the sulfur solubility and TCLP constraints. The measured property values were compared to predicted values using EWG2.5 and a selection of other existing models. A few models (e.g., electrical conductivity, TCLP) were found to be adequate for designing DFHLW glasses in the near future, while others require refits or offsets. It is recommended that new property models be developed for EWG3.0, as a large amount of DFHLW glass property data (> 14 × existing data) is expected to be collected in the compositional spaces where no data was previously available. To enable near-term calculations and formulations for designing DFHLW glasses and processing rate estimations, a formulation algorithm with minor modifications will be developed, EWG2.6.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2025

Block-Structured Operator Inference for Coupled Multiphysics Model Reduction

This work presents a block-structured formulation of Operator Inference as a way to learn structured reduced-order models for multiphysics systems. The approach specifies the governing equation structure for each physics component and the structure of the coupling terms. Once the multiphysics structure is specified, the reduced-order model is learned from snapshot data following the nonintrusive Operator Inference methodology. In addition to preserving physical system structure, which in turn permits preservation of system properties such as stability and second-order structure, the block-structured approach has the advantages of reducing the overall dimensionality of the learning problem and admitting tailored regularization for each physics component. The numerical advantages of the block-structured formulation over a monolithic Operator Inference formulation are demonstrated for aeroelastic analysis, which couples aerodynamic and structural models. For the benchmark test case of the AGARD 445.6 wing, block-structured Operator Inference provides an average 20% online prediction speedup over monolithic Operator Inference across subsonic and supersonic flow conditions in both the stable and fluttering parameter regimes while preserving the accuracy achieved with monolithic Operator Inference.

42 ENGINEERING↗
DOE OSTI2025

Nanolipoprotein particle (NLP) vaccine confers protection against Yersinia pestis aerosol challenge in a BALB/c mouse model

Introduction: Yersinia pestis is the etiological agent of plague, a disease that remains a concern as demonstrated by recent outbreaks in Madagascar. Infection with Y. pestis results in a rapidly progressing illness that can only be successfully treated with antibiotics given shortly after symptom onset. Live attenuated or whole cell inactivated vaccines confer protection against bubonic plague, but pneumonic plague has been more difficult to prevent. Novel effective subunit vaccine formulations may circumvent some of these shortfalls. Here, we compare the immunogenicity generated by an advanced subunit vaccine (F1V fusion protein) and a nanolipoprotein particle (NLP)-based vaccine. Methods: The NLP, a high-density lipoprotein mimetic, provides a nanoscale delivery platform for recombinant Y. pestis antigens LcrV (V) and F1. BALB/c mice were immunized via subcutaneous injection twice, three or four weeks apart. Four weeks later, splenocytes and sera were collected for immune profiling, and mice were challenged with aerosolized Y. pestis CO92. Results: Both formulations induced a strong IgG response against the F1 and V proteins, along with a robust memory B cell response and a balanced cell-mediated immune response as evidenced by both Th1- and Th2-related cytokines. The NLP-based vaccine induced a stronger cytokine response against F1, V, and F1V proteins relative to the F1V vaccine. As with F1V, the inclusion of Alhydrogel (Alu) in NLP vaccine formulations was critical for enhanced immunogenicity and protective efficacy. Mice that received two doses of F1:V:NLP + Alu and CpG were completely protected from a challenge with approximately eight median lethal doses of aerosolized Y. pestis CO92 and this protection confirmed the well-documented synergy between the F1 and V antigens in context of pneumonic plague. The NLPs have defined regions of polarity that facilitates the incorporation of a wide range of adjuvants and antigens with distinct physicochemical properties and are an excellent candidate platform for the development of multi-antigen vaccines.

F1↗
DOE OSTI2025

A Finite Element Method for Compressible and Turbulent Multiphase Flow Instabilities with Heat Transfer

We present a new finite element framework for modeling compressible, turbulent multiphase flows with heat transfer. For two-fluid systems with a free surface, the Volume of Fluid (VOF) method is implemented without the need for interface reconstruction, while turbulence is resolved using a dynamic Vreman large eddy simulation (LES) model. Unlike most two-phase VOF studies, which neglect heat transfer, the present approach incorporates energy transport equations within the VOF formulation to account for heat exchange, an effect particularly important in turbulent flows. Conjugate heat transfer is often challenging in finite volume methods, which require explicit specification of heat fluxes at the solid–fluid interface, limiting accuracy and predictive capability. By contrast, the finite element formulation does not require heat flux inputs, allowing more accurate and robust simulation of heat transfer between solids and fluids. The method is demonstrated through three representative cases. First, a two-fluid instability with a single-mode perturbation is simulated and validated against analytical growth rates. Second, conjugate heat transfer is examined in a high-temperature flow over a cold metal cylinder, with validation performed both quantitatively—via pressure coefficient comparisons with experimental data—and qualitatively using vector field topology. Finally, compressible spray injection and breakup are modeled, demonstrating the ability of the framework to capture interfacial dynamics and atomization under turbulent, high-speed conditions. In the compressible spray injection and breakup case, the results indicate that the finite element formulation achieved higher predictive accuracy and robustness than the finite-volume method. With the same mesh resolution, the FEM reduced the root mean square error (RMSE) and mean absolute percentage error (MAPE) from 6.96 mm and 26.0% (for the FVM) to 4.85 mm and 12.7%, respectively, demonstrating improved accuracy and robustness in capturing interfacial dynamics and heat transfer. The study also introduced vector field topology to visualize and interpret coherent flow structures and instabilities, offering insights beyond conventional scalar-field analyses.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Structural Aspects of Neutron Survival Probabilities

The neutron survival probability (and related quantities including probabilities of extinction and initiation) is a central element of the broader stochastic theory of neutron populations and finds application in fields including reactor start-up, analysis of reactor power bursts and criticality accidents, and safeguards. In a full neutron transport formulation, the equation governing the single-neutron survival probability is a backward or adjoint-like integro-partial differential equation with the added complexity of being highly nonlinear. Analogous formulations of this equation exist in the context of many approximate theories of neutron transport, with the point kinetics formulation having received significant theoretical attention since the 1940s. This work continues this tradition by providing a novel analysis of the single-neutron survival probability equation using the tools of boundary layer theory. The analysis reveals that the “fully dynamic” solution of the single-neutron survival probability equation—and some key probability distributions derived from it—may be cast as a singular perturbation around the underlying quasi-static single-neutron probability of initiation. In this perturbation solution, the expansion parameter is the ratio of the neutron generation time to a macroscopic time scale characterizing the overall system evolution; this interpretation illuminates some of the fundamental structural aspects of neutron survival phenomena.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Influence of Carboxymethyl Cellulose as a Thickening Agent for Glauber’s Salt-Based Low Temperature PCM

This work is focused on a novel, promising low temperature phase change material (PCM), based on the eutectic Glauber’s salt composition. To allow phase transition within the refrigeration range of temperatures of +5 °C to +12 °C, combined with a high repeatability of melting–freezing processes, and minimized subcooling, the application of three variants of sodium carboxymethyl cellulose (Na-CMC) with distinct molecular weights (700,000, 250,000, and 90,000) is considered. The primary objective is to optimize the stabilization of this eutectic PCM formulation, while maintaining the desired enthalpy level. Preparation methods are refined to ensure repeatability in mixing components, thereby optimizing performance and stability. Additionally, the influence of Na-CMC molecular weight on stabilization is examined through differential scanning calorimetry (DSC), T-history, and rheology tests. The PCM formulation of interest builds upon prior research in which borax, ammonium chloride, and potassium chloride were used as additives to sodium sulfate decahydrate (Glauber’s salt), prioritizing environmentally responsible materials. The results reveal that CMC with molecular weights of 250 kg/mol and 90 kg/mol effectively stabilize the PCM without phase separation issues, slowing crystallization kinetics. Conversely, CMC of 700 kg/mol proved ineffective due to the disruption of gel formation at its low gel point, hindering higher concentrations. Calculations of ionic concentration indicate higher Na ion content in PCM stabilized with 90 kg/mol CMC, suggesting increased ionic interactions and gel strength. A tradeoff is discovered between the faster crystallization in lower molecular weight CMC and the higher concentration required, which increases the amount of inert material that does not participate in the phase transition. After thermal cycling, the best formulation had a latent heat of 130 J/g with no supercooling, demonstrating excellent performance. This work advances PCM’s reliability as a thermal energy storage solution for diverse applications and highlights the complex relationship between Na-CMC molecular weight and PCM stabilization.

25 ENERGY STORAGE↗
DOE OSTI2025

Relationship Between Filler Type, Thermomechanical Properties, and Aging of RTV Silicone Foams

Room-temperature vulcanizing (RTV) silicone foams are used in many industrial applications that require the material to perform over long time periods. However, mechanical properties tend to deteriorate when these foams age under a compressive load. The chemical aging is attributed to the presence of unreacted functional groups of the prepolymers, residues from acid, and catalytically active tin (II) species. Here, an optimized thermal treatment of an RTV foam that achieves completion of curing reactions and deactivation of reactive species is proposed. Foams that were thermally aged for three months under compressive load showed no signs of compression set, indicative of the effectiveness of the implemented post-curing approach. In addition, the effects of fillers (diatomaceous earth, fumed silica, and carbon nanofibers) on thermomechanical properties were investigated. Tensile strength, tear strength, and thermal conductivity increased when these fillers were added to the unfilled RTV formulation, with carbon nanofibers (CNFs) being the most effective filler. Rheological studies of RTV formulations indicated that 2.5 wt.% of CNFs is the upper limit that can be added to the RTV formulation.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Cell-Free Screening, Production and Animal Testing of a STI-Related Chlamydial Major Outer Membrane Protein Supported in Nanolipoproteins

Background: Vaccine development against Chlamydia, a prevalent sexually transmied infection (STI), is imperative due to its global public health impact. However, significant challenges arise in the production of effective subunit vaccines based on recombinant protein antigens, particularly with membrane proteins like the Major Outer Membrane Protein (MOMP). Methods: Cellfree protein synthesis (CFPS) technology is an aractive approach to address these challenges as a method of high-throughput membrane protein and protein complex production coupled with nanolipoprotein particles (NLPs). NLPs provide a supporting scaffold while allowing easy adjuvant addition during formulation. Over the last decade, we have been working toward the production and characterization of MOMP-NLP complexes for vaccine testing. Results: The work presented here highlights the expression and biophysical analyses, including transmission electron microscopy (TEM) and dynamic light scaering (DLS), which confirm the formation and functionality of MOMP-NLP complexes for use in animal studies. Moreover, immunization studies in preclinical models compare the past and present protective efficacy of MOMP-NLP formulations, particularly when co-adjuvanted with CpG and FSL1. Conclusion: Ex vivo assessments further highlight the immunomodulatory effects of MOMP-NLP vaccinations, emphasizing their potential to elicitrobust immune responses. However, further research is warranted to optimize vaccine formulations further, validate efficacy against Chlamydia trachomatis, and beer understand the underlying mechanisms of immune response.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Performance-portable Binary Neutron Star Mergers with AthenaK

We introduce an extension to the AthenaK code for general-relativistic magnetohydrodynamics (GRMHD) in dynamical spacetimes using a 3+1 conservative Eulerian formulation. Like the fixed-spacetime GRMHD solver, we use standard finite-volume methods to evolve the fluid and a constrained-transport scheme to preserve the divergence-free constraint for the magnetic field. We also utilize a first-order flux correction (FOFC) scheme to reduce the need for an artificial atmosphere and optionally enforce a maximum principle to improve robustness. We demonstrate the accuracy of AthenaK using a set of standard tests in flat and curved spacetimes. Using a SANE accretion disk around a Kerr black hole, we compare the new solver to the existing solver for stationary spacetimes using the so-called "HARM-like" formulation. We find that both formulations converge to similar results. We also include the first published binary neutron star (BNS) mergers performed on graphical processing units (GPUs). Thanks to the FOFC scheme, our BNS mergers maintain a relative error of $\mathcal{O}$(10 –11 ) or better in baryon mass conservation up to collapse. Finally, we perform scaling tests of AthenaK on OLCF Frontier, where we show excellent weak scaling of ≥80% efficiency up to 32,768 GPUs and 74% up to 65,536 GPUs for a GRMHD problem in dynamical spacetimes with six levels of mesh refinement. AthenaK achieves an order-of-magnitude speedup using GPUs compared to CPUs, demonstrating that it is suitable for performing numerical relativity problems on modern exascale resources.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

A novel transformation of the ice sheet Stokes equations and some of its properties and applications

We introduce a novel transformation of the Stokes equations into a form closely resembling the shallow Blatter–Pattyn equations. The two forms differ by only a few additional terms, while their variational formulations differ only by a single term in each horizontal direction. Specifically, the variational formulation of the Blatter–Pattyn model drops the vertical velocity in the second invariant of the strain rate tensor. Here we make use of the new transformation in two ways. First, we consider incorporating the transformed equations into a code that can be very easily converted from a Stokes to a Blatter–Pattyn model, and vice versa, by switching these terms on or off. This may be generalized so that the Stokes model is switched on adaptively only where the Blatter–Pattyn model loses accuracy. Second, the key role played by the vertical velocity in the Blatter–Pattyn approximation motivates new approximations. Two examples are presented. These require a mesh that enables the discrete continuity equation to be invertible for the vertical velocity in terms of the horizontal velocity components. Examples of such meshes, such as the first-order P1–E0 mesh and the second-order P2–E1 mesh, are given in both 2D and 3D. However, the transformed Stokes model has the same type of gravity forcing as the Blatter–Pattyn model, determined by the ice surface slope, thereby forgoing some of the mesh generality of the traditional formulation of the Stokes model.

58 GEOSCIENCES↗