Search NASA⌕ Search

SEARCH · Search NASA

Results for “Two point model”

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 775 records · Page 43
All sourcesNASA NTRSDOE OSTI
NASA NTRS2006

Viking Afterbody Heating Computations and Comparisons to Flight Data

Computational fluid dynamics predictions of Viking Lander 1 entry vehicle afterbody heating are compared to flight data. The analysis includes a derivation of heat flux from temperature data at two base cover locations, as well as a discussion of available reconstructed entry trajectories. Based on the raw temperature-time history data, convective heat flux is derived to be 0.63-1.10 W/sq cm for the aluminum base cover at the time of thermocouple failure. Peak heat flux at the fiberglass base cover thermocouple is estimated to be 0.54-0.76 W/sq cm, occurring 16 seconds after peak stagnation point heat flux. Navier-Stokes computational solutions are obtained with two separate codes using an 8-species Mars gas model in chemical and thermal non-equilibrium. Flowfield solutions using local time-stepping did not result in converged heating at either thermocouple location. A global time-stepping approach improved the computational stability, but steady state heat flux was not reached for either base cover location. Both thermocouple locations lie within a separated flow region of the base cover that is likely unsteady. Heat flux computations averaged over the solution history are generally below the flight data and do not vary smoothly over time for both base cover locations. Possible reasons for the mismatch between flight data and flowfield solutions include underestimated conduction effects and limitations of the computational methods.

Edquist, Karl T.↗
NASA NTRS2006

Viking Afterbody Heating Computations and Comparisons to Flight Data

Computational fluid dynamics predictions of Viking Lander 1 entry vehicle afterbody heating are compared to flight data. The analysis includes a derivation of heat flux from temperature data at two base cover locations, as well as a discussion of available reconstructed entry trajectories. Based on the raw temperature-time history data, convective heat flux is derived to be 0.63-1.10 W/cm2 for the aluminum base cover at the time of thermocouple failure. Peak heat flux at the fiberglass base cover thermocouple is estimated to be 0.54-0.76 W/cm2, occurring 16 seconds after peak stagnation point heat flux. Navier-Stokes computational solutions are obtained with two separate codes using an 8- species Mars gas model in chemical and thermal non-equilibrium. Flowfield solutions using local time-stepping did not result in converged heating at either thermocouple location. A global time-stepping approach improved the computational stability, but steady state heat flux was not reached for either base cover location. Both thermocouple locations lie within a separated flow region of the base cover that is likely unsteady. Heat flux computations averaged over the solution history are generally below the flight data and do not vary smoothly over time for both base cover locations. Possible reasons for the mismatch between flight data and flowfield solutions include underestimated conduction effects and limitations of the computational methods.

Edquist, Karl T.↗
NASA NTRS2007

Analysis of Performance of Stereoscopic-Vision Software

A team of JPL researchers has analyzed stereoscopic vision software and produced a document describing its performance. This software is of the type used in maneuvering exploratory robotic vehicles on Martian terrain. The software in question utilizes correlations between portions of the images recorded by two electronic cameras to compute stereoscopic disparities, which, in conjunction with camera models, are used in computing distances to terrain points to be included in constructing a three-dimensional model of the terrain. The analysis included effects of correlation- window size, a pyramidal image down-sampling scheme, vertical misalignment, focus, maximum disparity, stereo baseline, and range ripples. Contributions of sub-pixel interpolation, vertical misalignment, and foreshortening to stereo correlation error were examined theoretically and experimentally. It was found that camera-calibration inaccuracy contributes to both down-range and cross-range error but stereo correlation error affects only the down-range error. Experimental data for quantifying the stereo disparity error were obtained by use of reflective metrological targets taped to corners of bricks placed at known positions relative to the cameras. For the particular 1,024-by-768-pixel cameras of the system analyzed, the standard deviation of the down-range disparity error was found to be 0.32 pixel.

Kim, Won↗
NASA NTRS2010

Composite Overwrapped Pressure Vessels (COPV) Stress Rupture Test: Part 2

One of the major concerns for the aging Space Shuttle fleet is the stress rupture life of composite overwrapped pressure vessels (COPVs). Stress rupture life of a COPY has been defined as the minimum time during which the composite maintains structural integrity considering the combined effects of stress levels and time. To assist in the evaluation of the aging COPVs in the Orbiter fleet an analytical reliability model was developed. The actual data used to construct this model was from testing of COPVs constructed of similar, but not exactly same materials and pressure cycles as used on Orbiter vessels. Since no actual Orbiter COPV stress rupture data exists the Space Shuttle Program decided to run a stress rupture test to compare to model predictions. Due to availability of spares, the testing was unfortunately limited to one 40" vessel. The stress rupture test was performed at maximum operating pressure at an elevated temperature to accelerate aging. The test was performed in two phases. The first phase, 130 F, a moderately accelerated test designed to achieve the midpoint of the model predicted point reliability. A more aggressive second phase, performed at 160 F, was designed to determine if the test article will exceed the 95% confidence interval ofthe model. In phase 3, the vessel pressure was increased to above maximum operating pressure while maintaining the phase 2 temperature. After reaching enough effectives hours to reach the 99.99% confidence level of the model phase 4 testing began when the temperature was increased to greater than 170 F. The vessel was maintained at phase 4 conditions until it failed after over 3 million effect hours. This paper will discuss the results of this test, it's implications and possible follow-on testing.

Russell, Richard↗
NASA NTRS2009

Composite Overwrapped Pressure Vessels (COPV) Stress Rupture Test

One of the major concerns for the aging Space Shuttle fleet is the stress rupture life of composite overwrapped pressure vessels (COPVs). Stress rupture life of a COPV has been defined as the minimum time during which the composite maintains structural integrity considering the combined effects of stress levels and time. To assist in the evaluation of the aging COPVs in the Orbiter fleet an analytical reliability model was developed. The actual data used to construct this model was from testing of COPVs constructed of similar, but not exactly same materials and pressure cycles as used on Orbiter vessels. Since no actual Orbiter COPV stress rupture data exists the Space Shuttle Program decided to run a stress rupture test to compare to model predictions. Due to availability of spares, the testing was unfortunately limited to one 40" vessel. The stress rupture test was performed at maximum operating pressure at an elevated temperature to accelerate aging. The test was performed in two phases. The first phase, 130 F, a moderately accelerated test designed to achieve the midpoint of the model predicted point reliability. The more aggressive second phase, performed at 160 F was designed to determine if the test article will exceed the 95% confidence interval of the model. This paper will discuss the results of this test, it's implications and possible follow-on testing.

Russell, Richard↗
NASA NTRS2011

Sloshing Gas in the Core of the Most Luminous Galaxy Cluster RXJ1347.5-1145

We present new constraints on the merger history of the most X-ray luminous cluster of galaxies, RXJ1347.5-1145, based its unique multiwavelength morphology. Our X-ray analysis confirms the core gas is undergoing "sloshing" resulting from a prior, large scale, gravitational perturbation. In combination with extensive multiwavelength observations, the sloshing gas points to the primary and secondary clusters having had at least two prior strong gravitational interactions. The evidence supports a model in which the secondary subcluster with mass M=4.8+/-2.4 x 10(exp 14) Stellar Mass has previously (> or approx.=0.6 Gyr ago) passed by the primary cluster, and has now returned for a subsequent crossing where the subcluster's gas has been completely stripped from its dark matter halo. RXJ1347 is a prime example of how core gas sloshing may be used to constrain the merger histories of galaxy clusters through multiwavelength analyses.

Markevitch, Maxim↗
NASA NTRS2011

Sloshing Gas in the Core of the Most Luminous Galaxy Cluster RXJ1347.5-1145

We present new constraints on the merger history of the most X-ray luminous cluster of galaxies, RXJ1347.5-1145, based on its unique multiwavelength morphology. Our X-ray analysis confirms the core gas is undergoing "sloshing" resulting from a prior, large scale, gravitational perturbation. In combination with extensive multiwavelength observations, the sloshing gas points to the primary and secondary clusters having had at least two prior strong gravitational interactions. The evidence supports a model in which the secondary subcluster with mass M=4.8+/-2.4x10(exp 14) solar Mass has previously (> or approx.0.6 Gyr ago) passed by the primary cluster, and has now returned for a subsequent crossing where the subcluster's gas has been completely stripped from its dark matter halo. RXJ1347 is a prime example of how core gas sloshing may be used to constrain the merger histories of galaxy clusters through multiwavelength analyses.

Johnson, Ryan E.↗
DOE OSTI2017

Delineating Facies Spatial Distribution by Integrating Ensemble Data Assimilation and Indicator Geostatistics with Level Set Transformation.

A new approach is developed to delineate the spatial distribution of discrete facies (geological units that have unique distributions of hydraulic, physical, and/or chemical properties) conditioned not only on direct data (measurements directly related to facies properties, e.g., grain size distribution obtained from borehole samples) but also on indirect data (observations indirectly related to facies distribution, e.g., hydraulic head and tracer concentration). Our method integrates for the first time ensemble data assimilation with traditional transition probability-based geostatistics. The concept of level set is introduced to build shape parameterization that allows transformation between discrete facies indicators and continuous random variables. The spatial structure of different facies is simulated by indicator models using conditioning points selected adaptively during the iterative process of data assimilation. To evaluate the new method, a two-dimensional semi-synthetic example is designed to estimate the spatial distribution and permeability of two distinct facies from transient head data induced by pumping tests. The example demonstrates that our new method adequately captures the spatial pattern of facies distribution by imposing spatial continuity through conditioning points. The new method also reproduces the overall response in hydraulic head field with better accuracy compared to data assimilation with no constraints on spatial continuity on facies.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2001

Modeling the Gas Dynamics Environment in a Subscale Solid Rocket Test Motor

Subscale test motors are often used for the evaluation of solid rocket motor component materials such as internal insulation. These motors are useful for characterizing insulation performance behavior, screening insulation material candidates and obtaining material thermal and ablative property design data. One of the primary challenges associated with using subscale motors however, is the uncertainty involved when extrapolating the results to full-scale motor conditions. These uncertainties are related to differences in such phenomena as turbulent flow behavior and boundary layer development, propellant particle interactions with the wall, insulation off-gas mixing and thermochemical reactions with the bulk flow, radiation levels, material response to the local environment, and other anomalous flow conditions. In addition to the need for better understanding of physical mechanisms, there is also a need to better understand how to best simulate these phenomena using numerical modeling approaches such as computational fluid dynamics (CFD). To better understand and model interactions between major phenomena in a subscale test motor, a numerical study of the internal flow environment of a representative motor was performed. Simulation of the environment included not only gas dynamics, but two-phase flow modeling of entrained alumina particles like those found in an aluminized propellant, and offgassing from wall surfaces similar to an ablating insulation material. This work represents a starting point for establishing the internal environment of a subscale test motor using comprehensive modeling techniques, and lays the groundwork for improving the understanding of the applicability of subscale test data to full-scale motors. It was found that grid resolution, and inclusion of phenomena in addition to gas dynamics, such as two-phase and multi-component gas composition are all important factors that can effect the overall flow field predictions.

Eaton, Andrew M.↗
NASA NTRS2024

Wind Tunnel Test of A Scale Model of A Venus Probe to Determine Aeroacoustics Environment

The Zephyr probe for NASA’s DAVINCI project is being designed to fly to Venus and collect measurements to characterize the Venusian atmosphere while it descends towards the surface. Despite a slow descent speed as the probe approaches the surface, the very dense atmosphere means the vibro-acoustic environment that the probe must withstand could be very harsh. A wind tunnel test campaign was conducted to characterize the surface-pressure fluctuations (acoustics) on the external surface of the probe during the descent. The goal was to provide forcing functions for vibro-acoustic analysis of all instruments inside the probe. A 25% scale-model of the Zephyr probe was manufactured, instrumented with 18 microphones, and tested in two wind tunnels located in the Fluid Mechanics Lab at NASA Ames Research Center at velocities ranging from 13 m/s to 48 m/s. Trip strips were added to the model to ensure a turbulent boundary layer on the model. Microphones were located appropriately to provide spectra of pressure fluctuations in regions with different flow characteristics as well as for calculating various two-point statistics. The acoustic spectra measured in the wind tunnel tests were scaled up to flight-scale using the estimated flight trajectory information. Besides fluctuating surface-pressure measurements, hot-film velocity measurements as well as smoke-laser flow visualization were conducted to better understand the flow around the probe. Lastly, several different geometry configurations were tested to see whether any reduction in acoustic levels could be achieved, but none of the configurations that were tested offered significant improvements. The flow visualization showed 3 primary flow regimes on the model: an attached turbulent boundary layer on the front of the model, a region of separated flow upstream of the flared lip, and a large, separated wake at the rear of the model. The acoustic spectra generally show three different shapes corresponding to these flow regions, with the highest levels seen at the rear of the model, and directly in front of the drag plate. The scaled overall levels of pressure fluctuations were found to be high: in the range of 130dB to 155dB.

planetary probe↗
DOE OSTI2022

Constraining IGM enrichment and metallicity with the C iv forest correlation function

ABSTRACT The distribution and abundance of metals in the diffuse intergalactic medium (IGM) have implications for galaxy formation and evolution models, and has been argued to be sensitive to the Universe’s reionization history. However, reduced sensitivity in the near-IR implies that probing IGM metals at z > 4 is currently out of reach with the traditional method of detecting individual absorbers. We present a new technique based on clustering analysis that enables the detection of these weak IGM absorbers. We investigate the two-point correlation function (2PCF) of the ${\rm C\, {\small IV}}$ forest as a probe of IGM metallicity and enrichment topology by simulating the z = 4.5 IGM with models of inhomogeneous metal distributions. The 2PCF of the ${\rm C\, {\small IV}}$ forest demonstrates a clear peak at a characteristic separation corresponding to the doublet separation of the ${\rm C\, {\small IV}}$ line.The peak amplitude scales quadratically with metallicity, while enrichment topology affects both the shape and amplitude of the 2PCF. For models consistent with the distribution of metals at z ∼ 3, we find that we can constrain [C/H] to within 0.2 dex, log$\, M_{\rm {min}}$ to within 0.4 dex, and R to within 15 per cent. We show that CGM absorbers can be reliably identified and masked, thus recovering the underlying IGM signal. The auto-correlation of the metal-line forest presents a compelling avenue to constrain the IGM metallicity and enrichment topology with high precision at z > 4, thereby pushing such measurements into the Epoch of Reionization.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2021

Magnetic Structure and Propagation of Two Interacting CMEs from the Sun to Saturn

One of the grand challenges in heliophysics is the characterization of coronal mass ejection (CME) magnetic structure and evolution from eruption at the Sun through heliospheric propagation. At present, the main difficulties are related to the lack of direct measurements of the coronal magnetic fi elds and the lack of 3D in-situ measurements of the CME body in interplanetary space. Nevertheless, the evolution of a CME magnetic structure can be followed using a combination of multi-point remote-sensing observations and multi-spacecraft in-situ measurements as well as modelling. Accordingly, we present in this work the analysis of two CMEs that erupted from the Sun on 28 April 2012. We follow their eruption and early evolution using remote-sensing data, fi nding indications of CME{CME interaction, and then analyse their interplanetary counterpart(s) using in-situ measurements at Venus, Earth, and Saturn. We observe a seemingly single flux rope at all locations, but find possible signatures of interaction at Earth, where high-cadence plasma data are available. Reconstructions of the in-situ flux ropes provide almost identical results at Venus and Earth but show greater discrepancies at Saturn, suggesting that the CME was highly distorted and/or that further interaction with nearby solar wind structures took place before 10 AU. This work highlights the difficulties in connecting structures from the Sun to the outer heliosphere and demonstrates the importance of multi-spacecraft studies to achieve a deeper understanding of the magnetic confi guration of CMEs.

Erika Palmerio↗
DOE OSTI2025

Towards a Deeper Fundamental Understanding of (Al,Sc)N Ferroelectric Nitrides

Density functional theory (DFT) calculations, within the virtual crystal alloy approximation, are performed, along with the development of a Landau-type model employing a symmetry-allowed analytical expression of the internal energy and having parameters determined from first principles, to investigate properties and energetics of Al1-xScxN ferroelectric nitrides in their hexagonal forms. These DFT computations and this model predict the existence of two different types of minima, namely, the fourfold-coordinated wurtzite (WZ) polar structure and a five-fold coordinated paraelectric hexagonal phase (denoted as H5), for any Sc composition up to 40%. The H5 minimum progressively becomes the lowest-energy state within hexagonal symmetry as the Sc concentration increases from 0 to 0.4. Furthermore, the model points to several key findings. Examples include the crucial role of the coupling between polarization and strains to create the WZ minimum, in addition to polar and elastic energies, and that the origin of the H5 state overcoming the WZ phase as the global minimum within hexagonal symmetry when increasing the Sc composition mostly lies in the compositional dependency of only two parameters-one linked to the polarization and another one being purely elastic in nature. Other examples are that forcing Al1-xScxN systems to have no or a weak change in lattice parameters when heating them allows us to reproduce their finite-temperature polar properties well and that a value of the axial ratio close to that of the ideal WZ structure implies a large polarization at low temperatures but not necessarily at high temperatures because of the ordered-disordered character of the temperature-induced formation of the WZ state. Such findings should allow for a better fundamental understanding of (Al,Sc)N ferroelectric nitrides, which may be used to design efficient devices having, e.g., low operating voltages.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2019

Kinetics of the γ–δ phase transition in energetic nitramine-octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine

The solid, secondary explosive nitramine-octahydro-1,3,5,7-tetranitro-1,3,5,7 or HMX has four different stable polymorphs which have different molecular conformations, crystalline structures, and densities, making structural phase transitions between these nontrivial. Previous studies of the kinetics of the β–δ HMX structural transition found this to happen by a nucleation and growth mechanism, where growth was governed by the heat of fusion, or melting, even though the phase transition temperature is more than 100 K below the melting point. A theory known as virtual melting could easily justify this since the large volume difference in the two phases creates a strain at their interface that can lower the melting point to the phase transition temperature through a relaxation of the elastic energy. To learn more about structural phase transitions in organic crystalline solids and virtual melting, here we use time-resolved X-ray diffraction to study another structural phase transition in HMX, γ–δ. Again, second order kinetics are observed which fit to the same nucleation and growth model associated with growth by melting even though the volume change in this transition is too small to lower the melting point by interfacial strain. To account for this, we present a more general model illustrating that melting over a very thin layer at the interface between the two phases reduces the total interfacial energy and is therefore thermodynamically favorable and can drive the structural phase transition in the absence of large volume changes. Finally, our work supports the idea that virtual melting may be a more generally applicable mechanism for structural phase transitions in organic crystalline solids.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Inference of Rock Flow and Mechanical Properties from Injection-Induced Microseismic Events During Geologic CO 2 Storage

Monitoring microseismic activities during CO 2 injection into geologic formations is important for ensuring the safety of the storage operations. The resulting data provide insight into the response of the storage formation to CO 2 injection and can be used to infer the underlying rock flow and mechanical properties. In this paper, assimilation of microseismic data is performed for dynamic characterization of the storage formation by using a stochastic simulation model to forecast the microseismic response of a geologic formation during CO 2 injection. Two modeling approaches are adopted to predict the space-time distribution of the injection-induced microseismicity. The first model is based on pore pressure relaxation assumption, while the second model uses coupled flow and geomechanics simulation to establish the complex physical relation between the storage formation properties and the corresponding microseismic responses during CO 2 injection. The stochastic predictive models in each case are used in ensemble data assimilation frameworks to estimate rock properties from the observed microseismic data. Two data assimilation methods are considered: (i) a new ensemble-based stochastic point process filter (EnPPF) that can directly integrate discrete microseismic events, and (ii) a variant of ensemble smoother, known as the ensemble smoother with multiple data assimilation (ES-MDA), which requires continuous representation of microseismic events for assimilation. The two methods are successfully applied to a geologically realistic model of the Farnsworth Field in Texas, with complex geologic flow units and interacting fault systems.

42 ENGINEERING↗
NASA NTRS1995

Two-micron Laser Atmospheric Wind Sounder (LAWS) pointing/tracking study

The objective of the study was to identify and model major sources of short-term pointing jitter for a free-flying, full performance 2 micron LAWS system and evaluate the impact of the short-term jitter on wind-measurement performance. A fast steering mirror controls system was designed for the short-term jitter compensation. The performance analysis showed that the short-term jitter performance of the controls system over the 5.2 msec round-trip time for a realistic spacecraft environment was = 0.3 micro rad, rms, within the specified value of less than 0.5 micro rad, rms, derived in a 2 micron LAWS System Study. Disturbance modes were defined for: (1) the Bearing and Power Transfer Assembly (BAPTA) scan bearing, (2) the spacecraft reaction wheel torques, and (3) the solar array drive torques. The scan bearing disturbance was found to be the greatest contributing noise source to the jitter performance. Disturbances from the fast steering mirror reaction torques and a boom-mounted cross-link antenna clocking were also considered but were judged to be small compared to the three principal disturbance sources above and were not included in the final controls analysis.

Manlief, Scott↗
NASA NTRS1993

Space transportation vehicle design evaluation using saturated designs

An important objective in the preliminary design and evaluation of space transportation vehicles is to find the best values of design variables that optimize the performance characteristic (e.g. dry weight). For a given configuration, the vehicle performance can be determined by the use of complex sizing and performance evaluation computer programs. These complex computer programs utilize iterative algorithms and they are generally too expensive and/or difficult to use directly in multidisciplinary design optimization. An alternative is to use response surface methodology (RSM) and obtain quadratic polynomial approximations to the functional relationships between performance characteristics and design variables. In RSM, these approximation models are then used to determine optimum design parameter values and for rapid sensitivity studies. Constructing a second-order model requires that 'n' design parameters be studied at least at 3 levels (values) so that the coefficients in the model can be estimated. There, 3(n) factorial experiments (point designs or observations) may be necessary. For small values of 'n' such as two or three, this design works well. However, when a large number of design parameters are under study, the number of design points required for a full-factorial design may become excessive. Fortunately, these quadratic polynomial approximations can be obtained by selecting an efficient design matrix using central composite designs (CCD) from design of experiments theory. Each unique point design from the CCD matrix is then conducted using computerized analysis tools (e.g. POST, CONSIZ, etc.). In the next step, least squares regression analysis is used to calculate the quadratic polynomial coefficients from the data. However, in some multidisciplinary applications involving a large number of design variables and several disciplines, the computerized performance synthesis programs may get too time consuming and expensive to run even with the use of efficient central composite designs. In such cases, it may be preferable to keep the number of design points to an absolute minimum and trade some model accuracy with cost. For this purpose, another class of experimental designs, called saturated D-optimal designs may be utilized for generating a matrix of vehicle designs. A design is called saturated when the number of design points is exactly equal to the number of terms in the model to be fitted. As a result, saturated designs require the absolute minimum number of design points ((n+1)(n+2)/s) to estimate the quadratic polynomial model coefficients. Saturated designs can be generated using the D-optimality criterion. A good saturated design should give rise to least squares estimates with minimum generalized variance.

Unal, Resit↗
DOE OSTI2026

MFC 5.0: An exascale many-physics flow solver

Many problems of interest in engineering, medicine, and the fundamental sciences rely on high-fidelity flow simulation, making performant computational fluid dynamics solvers a mainstay of the open-source software community. Previous work MFC 3.0 was made a published, documented, and open-source solver via Bryngelson et al. Comp. Phys. Comm. (2021) with numerous physical features, numerical methods, and scalable infrastructure. MFC 5.0 is a significant update to MFC 3.0, featuring a broad set of well-established and novel physical models and numerical methods, as well as the introduction of GPU and APU (or superchip) acceleration. Here, we exhibit state-of-the-art performance and ideal scaling on the first two exascale supercomputers, OLCF Frontier and LLNL El Capitan. Combined with MFC’s single-accelerator performance, MFC achieves exascale computation in practice, and achieved the largest-to-date public CFD simulation at 200 trillion grid points as a 2025 ACM Gordon Bell Prize finalist. New physical features include the immersed boundary method, N-fluid phase change, Euler–Euler and Euler–Lagrange sub-grid bubble models, fluid-structure interaction, hypo- and hyper-elastic materials, chemically reacting flow, two-material surface tension, magnetohydrodynamics (MHD), and more. Numerical techniques now represent the current state-of-the-art, including general relaxation characteristic boundary conditions, WENO variants, Strang splitting for stiff sub-grid flow features, and low Mach number treatments. Weak scaling to tens of thousands of GPUs on OLCF Summit and Frontier and LLNL El Capitan achieves efficiencies within 5% of ideal to over 90% of their respective system sizes. Strong scaling results for a 16-times increase in device count show parallel efficiencies over 90% on OLCF Frontier. MFC’s software stack has undergone further improvements, including continuous integration, which ensures code resilience and correctness through over 300 regression tests; metaprogramming, which reduces code length while maintaining performance portability; and code generation for computing chemical reactions

Computational fluid dynamics↗