Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “sub scale”

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 559 records · Page 31

Distributed Disdrometer and Rain Gauge Measurement Infrastructure Developed for GPM Ground Validation

Global Precipitation Mission (GPM)retrieval algorithm validation requires datasets characterizing the 4-D structure, variability, and correlation properties of hydrometeor particle size distributions (PSD) and accumulations over satellite fields of view (FOV;<10 km). Collection of this data provides a means to assess retrieval errors related to beam filling and algorithm PSD assumptions. Hence, GPM Ground Validation is developing a deployable network of precipitation gauges and disdrometers to provide fine-scale measurements of PSD and precipitation accumulation variability. These observations will be combined with dual-frequency, polarimetric, and profiling radar data in a bootstrapping fashion to extend validated PSD measurements to a large coverage domain. Accordingly, a total of 24 Parsivel disdrometers(PD), 5 3rd-generation 2D Video Disdrometers (2DVD), 70 tipping bucket rain gauges (TBRG),9 weighing gauges, 7 Hot-Plate precipitation sensors (HP), and 3 Micro Rain Radars (MRR) have been procured. In liquid precipitation the suite of TBRG, PD and 2DVD instruments will quantify a broad spectrum of rain rate and PSD variability at sub-kilometer scales. In the envisioned network configuration 5 2DVDs will act as reference points for 16 collocated PD and TBRG measurements. We find that PD measurements provide similar measures of the rain PSD as observed with collocated 2DVDs (e.g., D0, Nw) for rain rates less than 15 mm/hr. For heavier rain rates we will rely on 2DVDs for PSD information. For snowfall we will combine point-redundant observations of SWER distributed over three or more locations within a FOV. Each location will contain at least one fenced weighing gauge, one HP, two PDs, and a 2DVD. MRRs will also be located at each site to extend the measurement to the column. By collecting SWER measurements using different instrument types that employ different measurement techniques our objective is to separate measurement uncertainty from natural variability in SWER and PSD. As demonstrated using C3VP polarimetric radar, gauge, and 2DVD/PD datasets these measurements can be combined to bootstrap an area wide SWER estimate via constrained modification of density-diameter and radar reflectivity-snowfall relationships. These data will be combined with snowpack, airborne microphysics, radar, radiometer, and tropospheric sounding data to refine GPM snowfall retrievals. The gauge and disdrometer instruments are being developed to operate autonomously when necessary using solar power and wireless communications. These systems will be deployed in numerous field campaigns through 2016. Planned deployment of these systems include field campaigns in Finland (2010), Oklahoma (2011), Canada (2012) and North Carolina (2013). GPM will also deploy 20 pairs of TBRGs within a 25 km2 region along the Virginia coast under NASA NPOL radar coverage in order to quantify errors in point-area rainfall measurements.

Petersen, Walter A.↗

Modeling fluvial erosion on regional to continental scales

The fluvial system is a major concern in modeling landform evolution in response to tectonic deformation. Three stream bed types (bedrock, coarse-bed alluvial, and fine-bed alluvial) differ in factors controlling their occurrence and evolution and in appropriate modeling approaches. Spatial and temporal transitions among bed types occur in response to changes in sediment characteristics and tectonic deformation. Erosion in bedrock channels depends upon the ability to scour or pluck bed material; this detachment capacity is often a power function of drainage area and gradient. Exposure of bedrock in channel beds, due to rapid downcutting or resistant rock, slows the response of headwater catchments to downstream baselevel changes. Sediment routing through alluvial channels must account for supply from slope erosion, transport rates, abrasion, and sorting. In regional landform modeling, implicit rate laws must be developed for sediment production from erosion of sub-grid-scale slopes and small channels.

Howard, Alan D.↗

Resolving Mesoscale Convective Systems: Grid Spacing Sensitivity in the Tropics and Midlatitudes

Abstract Mesoscale convective systems (MCSs) are a critical global water cycle component and drive extreme precipitation events in tropical and midlatitude regions. However, simulating deep convection remains challenging for modern numerical weather and climate models due to the complex interactions of processes from microscales to synoptic scales. Recent models with kilometer‐scale horizontal grid spacings offer notable improvements in simulating deep convection compared to coarser‐resolution models. Still, deficiencies in representing key physical processes, such as entrainment, lead to systematic biases. Additionally, evaluating model outputs using process‐oriented observational data remain difficult. This study presents an ensemble of MCS simulations with spanning the deep convective gray zone ( from 12 km to 125 m) in the Southern Great Plains of the U.S. and the Amazon Basin. Comparing these simulations with Atmospheric Radiation Measurement (ARM) wind profiler observations, we find greater sensitivity in the Amazon Basin compared to the Great Plains. Convective drafts converge structurally at sub‐kilometer scales, but some deficiencies remain. In both regions, simulated up and downdrafts are too deep and extreme downdrafts are not strong enough. Furthermore, Amazonian updrafts are too strong. Overall, we observe higher sensitivity in the tropics, including an artificial buildup in vertical kinetic energy at scales of , suggesting a need for 250 m in this region. Nevertheless, bulk convergence—agreement of storm‐average statistics—is achievable with kilometer‐scale simulations within a 10% error margin with 1 km providing a good balance between accuracy and computational cost.

54 ENVIRONMENTAL SCIENCES↗

Formation of Large Regular Satellites of Giant Planets in an Extended Gaseous Nebula: Subnebula Model and Accretion of Satellites

We model the subnebulae of Jupiter and Saturn wherein satellite accretion took place. We expect a giant planet subnebula to be composed of an optically thick (given gaseous opacity) inner region inside of the planet's centrifugal radius (located at r(sub c, sup J) = l5R(sub J) for Jupiter and r(sub c, sup S) = 22R(sub S) for Saturn), and an optically thin, extended outer disk out to a fraction of the planet's Roche lobe, which we choose to be R(sub roche)/5 (located at approximately 150R(sub J) near the inner irregular satellites for Jupiter, and approximately 200R(sub S) near Phoebe for Saturn). This places Titan and Ganymede in the inner disk, Callisto and Iapetus in the outer disk, and Hyperion in the transition region. The inner disk is the leftover of the gas accreted by the protoplanet. The outer disk results from the solar torque on nebula gas flowing into the protoplanet during the time of giant planet gap opening. For the sake of specificity, we use a cosmic mixture 'minimum mass' model to constrain the gas densities of the inner disks of Jupiter and Saturn (and also Uranus). For the total mass of the outer disk we use the simple scaling M(sub disk) = M(sub P)tau(sub gap)/tau(sub acc), where M(sub P) is the mass of the giant planet, tau(sub gap) is the gap opening timescale, and tau(sub acc) is the giant planet accretion time. This gives a total outer disk mass of approximately 100M(sub Callisto) for Jupiter and possibly approximately 200M(sub Iapetus) for Saturn (which contain enough condensables to form Callisto and Iapetus respectively). Our model has Ganymede at a subnebula temperature of approximately 250 K and Titan at approximately 100 K. The outer disks of Jupiter and Saturn have constant temperatures of 130 K and 90 K respectively.

Mosqueira, I.↗

Transmission-distribution long-term volt-var planning considering reactive power support capability of distributed PV

High penetration of grid-edge, inverter-based photovoltaic (PV) can cause significant voltage fluctuations not only at the distribution but also at the sub-transmission levels due to PV output intermittency. Traditional reactive power planning approaches do not consider intermittency, nor the possibility of coordinating the control of existing and future volt-ampere reactive resources. This paper proposes a reactive power planning tool for sub-transmission systems to mitigate voltage violations and fluctuations caused by high PV penetration and intermittency with a minimum investment cost. The planning tool coordinates with an optimization-based volt-var operational tool for: a) modeling the coordination of all existing var assets in both sub-transmission and distribution systems to reduce the need of new equipment, and b)selecting a set of scenarios with voltage violations, derived from PV intermittency c) testing the final investment decision. The tool obtains an investment need for each intermittency scenario with a proposed optimal power-flow framework with efficient techniques to handle a high number of discrete variables. Two options are provided for final planning decision: i) a conservative direct combination of investment need solutions and ii) a machine learning-based selection of representative investment needs at most time steps. The final investment decision options are verified using a realistic large-scale sub-transmission system and 5-minute PV and load data. The results show a significant voltage performance improvement with a lower investment cost for additional var equipment compared to conventional approaches.

14 SOLAR ENERGY↗

Multiscale petrophysical characterization and flow unit classification of the Minnelusa eolian sandstones

Integration of petrophysical and geological information is critical to simulation of subsurface carbon storage (GCS). In this sense, two depositional facies were identified from the core description and well-log interpretation, namely massive (MS) and cross-bedded (CB) facies groups. Additionally, pore-scale characteristics were studied by a combination of techniques, e.g. Nuclear Magnetic Resonance (NMR) and mercury intrusion capillary pressure (MICP). Scanning electron microscope (SEM) and petrographic analyses show that the pore structure is dominantly controlled by the depositional environment and dolomite cementation. NMR-T 2 distributions of MS and CB facies show triple and quadruple modes, respectively. In addition, MICP of high- and low-permeability MS facies samples, and their CB facies group mixtures were collected. The MS sample pore-throat size distribution is uni-modal, while the triple-modal characteristic of the mixtures indicates heterogeneous pore structures at the sub-core scale for CB facies. The reliably estimates of porosity and permeability for both facies groups via NMR techniques and the MLR (Multiple Linear Regression) approach demonstrate the applicability of these techniques to eolian sandstone. Moreover, irreducible water saturation via the T 2 -cutoff method correlates strongly with T 2LM instead of porosity. Finally, the rock quality index and flow zone indicator were calculated based on Combinable Magnetic Resonance (CMR) log interpretations. This provides direct connection to properties measured in the well. Four flow units were classified for both facies groups. Finally, results show that better reservoir quality with significant heterogeneities is observed in the CB facies. This study highlights the importance integrating multiscale petrophysical properties including facies, pore architecture and diagenesis analysis with core- to log-scale property characterization. The results herein validate our reservoir characterization and flow unit classification in eolian reservoirs.

58 GEOSCIENCES↗

Integrating NDVI-Based Within-Wetland Vegetation Classification in a Land Surface Model Improves Methane Emission Estimations

Earth system models (ESMs) are a common tool for estimating local and global greenhouse gas emissions under current and projected future conditions. Efforts are underway to expand the representation of wetlands in the Energy Exascale Earth System Model (E3SM) Land Model (ELM) by resolving the simultaneous contributions to greenhouse gas fluxes from multiple, different, sub-grid-scale patch-types, representing different eco-hydrological patches within a wetland. However, for this effort to be effective, it should be coupled with the detection and mapping of within-wetland eco-hydrological patches in real-world wetlands, providing models with corresponding information about vegetation cover. In this short communication, we describe the application of a recently developed NDVI-based method for within-wetland vegetation classification on a coastal wetland in Louisiana and the use of the resulting yearly vegetation cover as input for ELM simulations. Processed Harmonized Landsat and Sentinel-2 (HLS) datasets were used to drive the sub-grid composition of simulated wetland vegetation each year, thus tracking the spatial heterogeneity of wetlands at sufficient spatial and temporal resolutions and providing necessary input for improving the estimation of methane emissions from wetlands. Our results show that including NDVI-based classification in an ELM reduced the uncertainty in predicted methane flux by decreasing the model’s RMSE when compared to Eddy Covariance measurements, while a minimal bias was introduced due to the resampling technique involved in processing HLS data. Our study shows promising results in integrating the remote sensing-based classification of within-wetland vegetation cover into earth system models, while improving their performances toward more accurate predictions of important greenhouse gas emissions.

54 ENVIRONMENTAL SCIENCES↗

Spatial and Temporal Uncertainty of Crop Yield Aggregations

The aggregation of simulated gridded crop yields to national or regional scale requires information on temporal and spatial patterns of crop-specific harvested areas. This analysis estimates the uncertainty of simulated gridded yield time series related to the aggregation with four different harvested area data sets. We compare aggregated yield time series from the Global Gridded Crop Model Inter-comparison project for four crop types from 14 models at global, national, and regional scale to determine aggregation-driven differences in mean yields and temporal patterns as measures of uncertainty. The quantity and spatial patterns of harvested areas differ for individual crops among the four datasets applied for the aggregation. Also simulated spatial yield patterns differ among the 14 models. These differences in harvested areas and simulated yield patterns lead to differences in aggregated productivity estimates, both in mean yield and in the temporal dynamics. Among the four investigated crops, wheat yield (17% relative difference) is most affected by the uncertainty introduced by the aggregation at the global scale. The correlation of temporal patterns of global aggregated yield time series can be as low as for soybean (r = 0.28).For the majority of countries, mean relative differences of nationally aggregated yields account for10% or less. The spatial and temporal difference can be substantial higher for individual countries. Of the top-10 crop producers, aggregated national multi-annual mean relative difference of yields can be up to 67% (maize, South Africa), 43% (wheat, Pakistan), 51% (rice, Japan), and 427% (soybean, Bolivia).Correlations of differently aggregated yield time series can be as low as r = 0.56 (maize, India), r = 0.05∗Corresponding (wheat, Russia), r = 0.13 (rice, Vietnam), and r = −0.01 (soybean, Uruguay). The aggregation to sub-national scale in comparison to country scale shows that spatial uncertainties can cancel out in countries with large harvested areas per crop type. We conclude that the aggregation uncertainty can be substantial for crop productivity and production estimations in the context of food security, impact assessment, and model evaluation exercises.

Aggregation uncertainty↗

High-Fidelity Simulations of Gravity Currents Using a High-Order Finite-Difference Spectral Vanishing Viscosity Approach

This numerical work investigates the potential of a high-order finite-difference spectral vanishing viscosity approach to simulate gravity currents at high Reynolds numbers. The method introduces targeted numerical dissipation at small scales through altering the discretisation of the second derivatives of the viscous terms in the incompressible Navier-Stokes equations to mimic the spectral vanishing viscosity (SVV) operator, originally designed for the regularisation of spectral element method (SEM) solutions of pure advection problems. Using a sixth-order accurate finite-difference scheme, the adoption of the SVV method is straightforward and comes with a negligible additional computational cost. In order to assess the ability of this high-order finite-difference spectral vanishing viscosity approach, we performed large-eddy simulations (LES) of a gravity current in a channelised lock-exchange set-up with our SVV model and with the well-known explicit static and dynamic Smagorinsky sub-grid scale (SGS) models. The obtained data are compared with a direct numerical simulation (DNS) based on more than 800 million mesh nodes, and with experimental measurements. A framework for the energy budget is introduced to investigate the behaviour of the gravity current. First, it is found that the DNS is in good agreement with the experimental data for the evolution of the front location and velocity field as well as for the stirring and mixing inside the gravity current. Secondly, the LES performed with less than 0.4% of the total number of mesh nodes compared to the DNS, can reproduce the main features of the gravity currents, with the SVV model yielding slightly more accurate results. It is also found that the dynamic Smagorinsky model performs better than its static version. For the present study, the static and dynamic Smagorinsky models are 1.8 and 2.5 times more expensive than the SVV model, because the latter does not require the calculation of explicit SGS terms in the Navier-Stokes equations nor spatial filtering operations.

direct numerical simulation↗

Use of In Situ and Airborne Multiangle Data to Assess MODIS and Landsat-Based Estimates of Directional Reflectance and Albedo

The quantification of uncertainty in satellite-derived global surface albedo products is a critical aspect in producing complete, physically consistent, and decadal land property data records for studying ecosystem change. A challenge in validating albedo measurements acquired from space is the ability to overcome the spatial scaling errors that can produce disagreements between satellite and field-measured values. Here, we present the results from an accuracy assessment of MODIS and Landsat-TM albedo retrievals, based on collocated comparisons with tower and airborne Cloud Absorption Radiometer (CAR) measurements collected during the 2007 Cloud and Land Surface Interaction Campaign (CLASIC). The initial focus was on evaluating intersensor consistency through comparisons of intrinsic bidirectional reflectance estimates. Local and regional assessments were then performed to obtain estimates of the resulting scaling uncertainties, and to establish the accuracy of albedo reconstructions during extended periods of precipitation. In general, the satellite derived estimates met the accuracy requirements established for the high-quality MODIS operational albedos at 500 m (the greater of 0.02 units or +/-10% of surface measured values). However, results reveal a high degree of variability in the root-mean-square error (RMSE) and bias of MODIS visible (0.3-0.7 μm) and Landsat-TM shortwave (0.3-5.0 μm) albedos; where, in some cases, retrieval uncertainties were found to be in excess of 15%. Results suggest that an overall improvement in MODIS shortwave albedo retrieval accuracy of 7.8%, based on comparisons between MODIS and CAR albedos, resulted from the removal of sub-grid scale mismatch errors when directly scaling-up the tower measurements to the MODIS satellite footprint.

biosphere↗

Few-cycle atomic layer deposition to nanoengineer polybenzimidazole for H 2 /CO 2 separation

Atomic layer deposition (ALD) creates uniform sub-nanometer films on a variety of surfaces and nanopore walls and has been used to modify polymers to improve surface affinity towards specific molecules, solvent resistance, and barrier properties to gases and vapors. Here, for the first time, we demonstrate that few-cycle ALD can be used to engineer functional polymers at a sub-nanometer scale to improve both molecular size-sieving ability and counterintuitively, gas permeability. Particularly, 1-cycle ALD treatment of polybenzimidazole (PBI) by sequential exposure to trimethylaluminum (TMA) and water vapor remarkably increases H 2 permeability by 120% - 270% and H 2 /CO 2 selectivity by 30% at 35–200 °C. The ALD not only deposits an AlOx layer on the surface but also enables the TMA to infiltrate and react with the bulk PBI to form an AlO x network, disrupting polymer chain packing and increasing chain rigidity. The membrane exhibits excellent stability when challenged with simulated syngas, overcoming the permeability/selectivity tradeoff for H 2 /CO 2 separation. In conclusion, this study showcases a facile and scalable way of engineering polymeric membranes at a sub-nanometer level to improve molecular separation performance.

36 MATERIALS SCIENCE↗

New numerical solutions of three-dimensional compressible hydrodynamic convection

Numerical solutions of three-dimensional compressible hydrodynamics (including sound waves) in a stratified medium with open boundaries are presented. Convergent/divergent points play a controlling role in the flows, which are dominated by a single frequency related to the mean sound crossing time. Superposed on these rapid compressive flows, slower eddy-like flows eventually create convective transport. The solutions contain small structures stacked on top of larger ones, with vertical scales equal to the local pressure scale heights, H sub p. Although convective transport starts later in the evolution, vertical scales of H sub p are apparently selected at much earlier times by nonlinear compressive effects.

Hossain, Murshed↗

The Global Atmosphere‐aerosol Model ICON‐A‐HAM2.3–Initial Model Evaluation and Effects of Radiation Balance Tuning on Aerosol Optical Thickness

Abstract The Hamburg Aerosol Module version 2.3 (HAM2.3) from the ECHAM6.3‐HAM2.3 global atmosphere‐aerosol model is coupled to the recently developed icosahedral nonhydrostatic ICON‐A (icon‐aes‐1.3.00) global atmosphere model to yield the new ICON‐A‐HAM2.3 atmosphere‐aerosol model. The ICON‐A and ECHAM6.3 host models use different dynamical cores, parameterizations of vertical mixing due to sub‐grid scale turbulence, and parameter settings for radiation balance tuning. Here, we study the role of the different host models for simulated aerosol optical thickness (AOT) and evaluate impacts of using HAM2.3 and the ECHAM6‐HAM2.3 two‐moment cloud microphysics scheme on several meteorological variables. Sensitivity runs show that a positive AOT bias over the subtropical oceans is remedied in ICON‐A‐HAM2.3 because of a different default setting of a parameter in the moist convection parameterization of the host models. The global mean AOT is biased low compared to MODIS satellite instrument retrievals in ICON‐A‐HAM2.3 and ECHAM6.3‐HAM2.3, but the bias is larger in ICON‐A‐HAM2.3 because negative AOT biases over the Amazon, the African rain forest, and the northern Indian Ocean are no longer compensated by high biases over the sub‐tropical oceans. ICON‐A‐HAM2.3 shows a moderate improvement with respect to AOT observations at AERONET sites. A multivariable bias score combining biases of several meteorological variables into a single number is larger in ICON‐A‐HAM2.3 compared to standard ICON‐A and standard ECHAM6.3. In the tropics, this multivariable bias is of similar magnitude in ICON‐A‐HAM2.3 and in ECHAM6.3‐HAM2.3. In the extra‐tropics, a smaller multivariable bias is found for ICON‐A‐HAM2.3 than for ECHAM6.3‐HAM2.3.

54 ENVIRONMENTAL SCIENCES↗

Maximum entropy modeling to identify physical drivers of shallow snowpack heterogeneity using unpiloted aerial system (UAS) lidar

Understanding the spatial variability of the snowpack is valuable for hydrologists and ecologists seeking to predict hydrological processes in a cold region. Snow distribution is a function of interactions among static variables, such as terrain, vegetation, and soil properties, and dynamic meteorological variables, such as solar radiation, wind speed and direction, and soil moisture. However, identifying the dominant physical drivers responsible for spatial patterns of the snowpack, particularly for ephemeral, shallow snowpacks, has been challenging due to the lack of the high-resolution snowpack and physical variables with high vertical accuracy as well as inherent limitations in traditional approaches. This study uses an Unpiloted Aerial System (UAS) lidar-based snow depth and static variables (1-m spatial resolution) to analyze field-scale spatial structures of snow depth and apply the Maximum Entropy (MaxEnt) model to identify primary controls over open terrain and forests at the University of New Hampshire Thompson Farm Research Observatory, New Hampshire, United States. We found that, among nine topographic and soil variables, plant functional type and terrain roughness contribute up to 80% and 76% of relative importance in the MaxEnt framework to predict locations of deeper or shallower snowpacks, respectively, across a mixed temperate forested and field landscape. Soil variables, such as organic matter and saturated hydraulic conductivity, were also important controls (up to 70% and 81%) on snow depth spatial variations for both open and forested landscapes suggesting spatial variations in soil variables under snow can control thermal transfer among soil, snowpack, and surface-atmosphere. This work contributes to improving land surface and snow models by informing parameterization of the sub-grid scale snow depths, down-scaling remotely sensed snow products, and understanding field scale snow states.

Eunsang Cho↗

NASA's High-Resolution GEOS Forecasting and Reanalysis Products: A Unified Tool from Local to Global Scales

NASA's GMAO produces high-resolution global forecasts for weather, aerosols, and air quality. The NASA Global Earth Observing System (GEOS) model has been expanded to provide global near-real-time 5-day forecasts of atmospheric chemical composition at unprecedented horizontal resolution of 0.25 degrees (~25 km), freely available to the public. This composition forecast system (GEOS-CF) combines the operational GEOS weather forecasting model with the state-of-the-science GEOS-Chem chemistry module to provide detailed analysis of a wide range of air pollutants such as ozone, carbon monoxide, nitrogen oxides, and fine particulate matter (PM2.5). GEOS-CF also assimilated satellite observations into the system for improved representation of weather and smoke. The assimilation system is currently being expanded to include chemically reactive trace gases. While the main focus of this new product is on tropospheric air quality information, the GEOS-Chem chemistry model used in this system includes the unified tropospheric stratospheric chemistry mechanism for improved forecasts of total column ozone during anomalous dynamical and chemical events. I will discuss current capabilities of the GEOS Constituent Data Assimilation System (CoDAS) to improve atmospheric composition modeling and possible future directions for GEOS-CF and reanalysis products. In addition, I will show how machine learning techniques can be used to correct for sub-grid-scale variability, which further improves model estimates at a given observation site.

Co-DAS↗

Hardware-accelerated inference for real-time gravitational-wave astronomy

The field of transient astronomy has seen a revolution with the first gravitational-wave detections and the arrival of multi-messenger observations they enabled. Transformed by the first detection of binary black hole and binary neutron star mergers, computational demands in gravitational-wave astronomy are expected to grow by at least a factor of two over the next five years as the global network of kilometer-scale interferometers are brought to design sensitivity. With the increase in detector sensitivity, real-time delivery of gravitational-wave alerts will become increasingly important as an enabler of multi-messenger followup. In this work, we report a novel implementation and deployment of deep learning inference for real-time gravitational-wave data denoising and astrophysical source identification. This is accomplished using a generic Inference-as-a-Service model that is capable of adapting to the future needs of gravitational-wave data analysis. Overall, our implementation allows seamless incorporation of hardware accelerators and also enables the use of commercial or private (dedicated) as-a-service computing. Based on our results, we propose a paradigm shift in low-latency and offline computing in gravitational-wave astronomy. Such a shift can address key challenges in peak-usage, scalability and reliability, and provide a data analysis platform particularly optimized for deep learning applications. The achieved sub-millisecond scale latency will also be relevant for any machine learning-based real-time control systems that may be invoked in the operation of near-future and next generation ground-based laser interferometers, as well as the front-end collection, distribution and processing of data from such instruments.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Improved Understanding of the Modeled QBO Using MLS Observations and MERRA Reanalysis

The Quasi-Biennial Oscillation (QBO) dominates the variability of the tropical stratosphere on interannual time scales. The QBO has been shown to extend its influence into the chemical composition of this region through dynamical mechanisms. We have started our analysis using the realistic QBO internally generated by the Goddard Earth Observing System Version 5 (GEOS-5) general circulation model coupled to a comprehensive stratospheric and tropospheric chemical mechanism forced with observed sea surface temperatures over the past 33 years. We will show targeted comparisons with observations from NASAs Aura satellite Microwave Limb Sounder (MLS) and the Modern Era Retrospective-Analysis for Research and Applications (MERRA) reanalysis to provide insight into the simulation of the primary and secondary circulations associated with the QBO. Using frequency spectrum analysis and multiple linear regression we can illuminate the resulting circulations and deduce the strengths and weaknesses in their modeled representation. Inclusion of the QBO in our simulation improves the representation of the subtropical barriers and overall tropical variability. The QBO impact on tropical upwelling is important to quantify when calculating trends in sub-decadal scale datasets.

MERRA↗

Machine Learning-Driven Conservative-to-Primitive Conversion in Hybrid Piecewise Polytropic and Tabulated Equations of State

We present a novel machine learning (ML)-based method to accelerate conservative-to-primitive inversion, focusing on hybrid piecewise polytropic and tabulated equations of state. Traditional root-finding techniques are computationally expensive, particularly for large-scale relativistic hydrodynamics simulations. To address this, we employ feedforward neural networks (NNC2PS and NNC2PL), trained in PyTorch (2.0+) and optimized for GPU inference using NVIDIA TensorRT (8.4.1), achieving significant speedups with minimal accuracy loss. The NNC2PS model achieves 𝐿 1 and 𝐿 ∞ errors of 4.54 × 10 −7 and 3.44 × 10−6, respectively, while the NNC2PL model exhibits even lower error values. TensorRT optimization with mixed-precision deployment substantially accelerates performance compared to traditional root-finding methods. Specifically, the mixed-precision TensorRT engine for NNC2PS achieves inference speeds approximately 400 times faster than a traditional single-threaded CPU implementation for a dataset size of 1,000,000 points. Ideal parallelization across an entire compute node in the Delta supercomputer (dual AMD 64-core 2.45 GHz Milan processors and 8 NVIDIA A100 GPUs with 40 GB HBM2 RAM and NVLink) predicts a 25-fold speedup for TensorRT over an optimally parallelized numerical method when processing 8 million data points. Moreover, the ML method exhibits sub-linear scaling with increasing dataset sizes. We release the scientific software developed, enabling further validation and extension of our findings. By exploiting the underlying symmetries within the equation of state, these findings highlight the potential of ML, combined with GPU optimization and model quantization, to accelerate conservative-to-primitive inversion in relativistic hydrodynamics simulations.

conservative-to-primitive conversion↗