Search NASA⌕ Search

SEARCH · Search NASA

Results for “data ordering”

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 847 records · Page 47

MARSAME Radiological Release Report for Replacement Water Lines from TA-48 to TA-55

Environmental Protection and Compliance, Environmental Stewardship (EPC-ES) has determined that only a portion of the soil associated with Replacement Water Lines from Technical Area (TA)-48 to TA55 (Figure 1) meets the criteria for unrestricted release to the public under Department of Energy (DOE) Order 458.1, Radiation Protection for the Public and the Environment (DOE, 2020). The remaining soil did not meet the criteria for unrestricted release and must be managed appropriately. These conclusions are based on the known history of the area combined with soil sample data collected in 2022; the findings are consistent with DOE Order 458.1 and Los Alamos National Laboratory (LANL) Functional Series Document EPC-ES-FSD-004, Environmental Radiation Protection (LANL, 2021). Sampling and data analysis, as described in this report, were sufficient to meet measurement quality objectives (MQOs) under the Multi-Agency Radiation Survey and Assessment of Materials and Equipment (MARSAME) manual (NUREG, 2009) and LANL procedures (LANL, 2020). Final approvals for waste disposition will come from LANL’s Waste Management Program.

54 ENVIRONMENTAL SCIENCES↗

MARSAME Release Report for Replacement Water Lines From TA-48 to TA-55, Revision 2

Environmental Protection and Compliance, Environmental Stewardship (EPC-ES) has determined that only a portion of the soil associated with Replacement Water Lines from Technical Area (TA)-48 to TA 55 (Figure 1) meets the criteria for unrestricted release to the public under Department of Energy (DOE) Order 458.1, Radiation Protection for the Public and the Environment (DOE, 2020). The remaining soil did not meet the criteria for unrestricted release and must be managed appropriately. These conclusions are based on the known history of the area combined with soil sample data collected in 2022; the findings are consistent with DOE Order 458.1 and Los Alamos National Laboratory (LANL) Functional Series Document EPC-ES-FSD-004, Environmental Radiation Protection (LANL, 2021). Sampling and data analysis, as described in this report, were sufficient to meet measurement quality objectives (MQOs) under the Multi-Agency Radiation Survey and Assessment of Materials and Equipment (MARSAME) manual (NUREG, 2009) and LANL procedures (LANL, 2020). Final approvals for waste disposition will come from LANL’s Waste Management Program.

54 ENVIRONMENTAL SCIENCES↗

First-order structural phase transition at low temperature in GaPt 5 P and its rapid enhancement with pressure

Single crystals of X Pt 5 ⁢P (X = Al, Ga, and In), belonging to the 1-5-1 family of compounds, were grown from a Pt-P solution at high temperatures, and measurements of the ambient pressure, temperature-dependent magnetization, resistivity, and x-ray diffraction were made. Additionally, the ambient-pressure Hall resistivity and temperature-dependent resistance under pressure were measured on GaPt 5 ⁢P. All three compounds have a tetragonal P4/mmm crystal structure at room temperature with metallic transport and weak diamagnetism over the 2–300 K temperature range. Surprisingly, at ambient pressure, both the transport and magnetization measurements on GaPt 5 ⁢P show a steplike feature in the 70–90 K region, suggesting a possible structural phase transition. Neither AlPt 5 ⁢P nor InPt 5 ⁢P have any signatures of a phase transition in their temperature-dependent electrical resistance and magnetization data. Both the hysteretic nature and sharpness of the features in the GaPt 5 ⁢P data suggest that the transition is first-order. Further, single-crystal x-ray diffraction measurements provided further details of the structural transition with a possibility of a crystal symmetry different from P⁢4/mmm below the transition temperature. The transition is characterized by anisotropic changes in the lattice parameters and a volume collapse with respect to the high-temperature tetragonal crystal structure. Furthermore, satellite peaks are observed at two distinct and nonequivalent wave vectors (0, 0, 0.5) and (0.5, 0.5, 0.5), and density functional theory calculations present phonon softening, especially at (0.5, 0.5, 0.5), as a possible driving mechanism. Additionally, we find that the structural transition temperature increases rapidly with increasing pressure, reaching room temperature by ~2.2 GPa, highlighting the high degree of pressure sensitivity of GaPt 5 ⁢P and fragile nature of its room-temperature structure. Even though the volume collapse and extreme pressure sensitivity suggest chemical pressure should drive a similar structural change in AlPt 5⁢ P, where both unit-cell dimensions and volume are smaller, its structure is found to be the same as that of the room-temperature GaPt 5 ⁢P. Overall, GaPt 5 ⁢P stands out as a sole member of the 1-5-1 family of compounds for which a temperature-driven structural change has been observed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Measurements of Zγ +jets differential cross sections in $pp$ collisions at $\sqrt{s}$ = 13 TeV with the ATLAS detector

Differential cross-section measurements of $Zγ$ production in association with hadronic jets are presented, using the full 139 fb -1 dataset of $\sqrt{s}$ = 13 TeV proton–proton collisions collected by the ATLAS detector during Run 2 of the LHC. Distributions are measured using events in which the $Z$ boson decays leptonically and the photon is usually radiated from an initial-state quark. Measurements are made in both one and two observables, including those sensitive to the hard scattering in the event and others which probe additional soft and collinear radiation. Different Standard Model predictions, from both parton-shower Monte Carlo simulation and fixed-order QCD calculations, are compared with the measurements. In general, good agreement is observed between data and predictions from MATRIX and MiNNLO PS , as well as next-to-leading-order predictions from MADGRAPH5_AMC@NLO and SHERPA.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Powder X-ray diffraction of fluorometholone, C 22 H 29 FO 4

Commercial fluorometholone, CAS #426-13-1, crystallizes in the monoclinic space group P 2 1 (#4) with a = 6.40648(2), b = 13.43260(5), c = 11.00060(8) Å, β = 92.8203(5)°, V = 945.517(5) Å 3 , and Z = 2. A reduced cell search in the Cambridge Structural Database yielded one previous structure determination, using single-crystal data at 292 K. In this work, the sample was ordered from the United States Pharmacopeial Convention (Lot # R032K0) and analyzed as-received. The room temperature (295 K) crystal structure was refined using synchrotron ( λ = 0.412826 Å) powder diffraction data and optimized using density functional theory (DFT) techniques. Hydrogen positions were included as a part of the structure and were re-calculated during the refinement. The diffraction data were collected on beamline 11-BM at the Advanced Photon Source, Argonne National Laboratory, and the powder X-ray diffraction pattern of the compound has been submitted to ICDD® for inclusion in the Powder Diffraction File™. The agreement of the Rietveld-refined and DFT-optimized structures is excellent; the root-mean-square Cartesian displacement is 0.060 Å. In addition to the O–H⋯O hydrogen bonds observed by Park et al. (Park, Y. J., Lee, M. Y., and Cho, S. I. ( 1992 ). “Fluorometholone,” J. Korean Chem. Soc. 36 , 812–817), C–H⋯O hydrogen bonds contribute to the crystal energy.

Gonzalez, Diana↗

Acceleration of Power System Dynamic Simulations Using a Deep Equilibrium Layer and Neural ODE Surrogate

The dominant paradigm for power system dynamic simulation is to build system-level simulations by combining physics-based models of individual components. The sheer size of the system along with the rapid integration of inverter-based resources exacerbates the computational burden of running time domain simulations. Here, in this paper, we propose a data-driven surrogate model based on implicit machine learningspecifically deep equilibrium layers and neural ordinary differential equationsto learn a reduced order model of a portion of the full underlying system. The data-driven surrogate achieves similar accuracy and reduction in simulation time compared to a physics-based surrogate, without the constraint of requiring detailed knowledge of the underlying dynamic models. This work also establishes key requirements needed to integrate the surrogate into existing simulation workflows; the proposed surrogate is initialized to a steady state operating point that matches the power flow solution by design.

Neural ordinary differential equations↗

Did You Win the GPU Cloud Lottery? Benchmarking from TFLOPS to Tokens/$

Cloud GPUs are commonly assumed to deliver consistent performance for a given GPU model. This assumption does not always hold: cloud providers employ diverse system configurations and virtualization mechanisms, and GPUs themselves exhibit non-negligible manufacturing variability (the silicon lottery). In this work, we present a large-scale measurement study of GPU performance variability across 11 cloud providers, covering over 3,500 physical GPUs and 6,800 benchmark runs. Our hierarchical analysis shows that while execution-level variation stays below 9%, performance varies by up to 38% across devices and providers for the same GPU model. Regression analysis indicates that driver- and OS-related software factors contribute less than 1% of the variance; instead, silicon lottery effects dominate observed performance variation, and cloud providers further amplify them through persistent, systematic second-order effects.

Slynko, Platon [Silicon Data, New York, USA] (ORCI↗

Whence the Interstellar Magnetic Field Shaping the Heliosphere?

Measurements of starlight polarized by aligned interstellar dust grains are used to probe the relation between the orientation of the ambient interstellar magnetic field (ISMF) and the ISMF traced by the ribbons of energetic neutral atoms discovered by the Interstellar Boundary Explorer spacecraft. We utilize polarization data, many acquired specifically for this study, to trace the configuration of the ISMF within 40 pc. A statistical analysis yields a best-fit ISMF orientation, B magpol , aligned with Galactic coordinates ℓ = 42°, b = 49°. Further analysis shows the ISMF is more orderly for "downfield" stars located over 90° from Bmagpol. The data subset of downfield stars yields an orientation for the nearby ISMF at ecliptic coordinates λ, β ≈ 219° ± 15°, 43° ± 9° (Galactic coordinates l, b ≈ 40°, 56°, ±17°). This best-fit ISMF orientation from polarization data is close to the field direction obtained from ribbon models. This agreement suggests that the ISMF shaping the heliosphere belongs to an extended ordered magnetic field. Extended filamentary structures are found throughout the sky. A previously discovered filament traversing the heliosphere nose region, "Filament A," extends over 300° of the sky, and crosses the upwind direction of interstellar dust flowing into the heliosphere. Filament A overlaps the locations of the Voyager kilohertz emissions, three quasar intraday variables, cosmic microwave background (CMB) components, and the inflow direction of interstellar grains sampled by Ulysses and Galileo. These features are likely located in the upstream outer heliosheath where ISMF drapes over the heliosphere, suggesting Filament A coincides with a dusty magnetized plasma. A filament 55° long is aligned with a possible shock interface between local interstellar clouds. A dark spot in the CMB is seen within 5° of the filament and within 10° of the downfield ISMF direction. Two large magnetic arcs are centered on the directions of the heliotail. The overlap between CMB components and the aligned dust grains forming Filament A indicates the configuration of dust entrained in the ISMF interacting with the heliosphere provides a measurable foreground to the CMB.

79 ASTRONOMY AND ASTROPHYSICS↗

Modelling nearest neighbour distributions of biased tracers using hybrid effective field theory

ABSTRACT We investigate the application of hybrid effective field theory (HEFT) – which combines a Lagrangian bias expansion with subsequent particle dynamics from N-body simulations – to the modelling of k-nearest neighbour cumulative distribution functions (kNN-CDFs) of biased tracers of the cosmological matter field. The kNN-CDFs are sensitive to all higher order connected N-point functions in the data, but are computationally cheap to compute. We develop the formalism to predict the kNN-CDFs of discrete tracers of a continuous field from the statistics of the continuous field itself. Using this formalism, we demonstrate how kNN-CDF statistics of a set of biased tracers, such as haloes or galaxies, of the cosmological matter field can be modelled given a set of low-redshift HEFT component fields and bias parameter values. These are the same ingredients needed to predict the two-point clustering. For a specific sample of haloes, we show that both the two-point clustering and the kNN-CDFs can be well-fit on quasi-linear scales (≳ 20h−1Mpc) by the second-order HEFT formalism with the same values of the bias parameters, implying that joint modelling of the two is possible. Finally, using a Fisher matrix analysis, we show that including kNN-CDF measurements over the range of allowed scales in the HEFT framework can improve the constraints on σ8 by roughly a factor of 3, compared to the case where only two-point measurements are considered. Combining the statistical power of kNN measurements with the modelling power of HEFT, therefore, represents an exciting prospect for extracting greater information from small-scale cosmological clustering.

79 ASTRONOMY AND ASTROPHYSICS↗

Ab initio leading order effective potential for elastic proton scattering based on the symmetry-adapted no-core shell model

Calculating microscopic optical potentials for elastic scattering at intermediate energies from light nuclei in an ab initio fashion within the Watson expansion has been established within the last few years. Based on the Watson expansion of the multiple scattering series, we employ a nonlocal translationally invariant nuclear density derived within the symmetry-adapted no-core shell model (SA-NCSM) framework from a chiral next-to-next-to-leading order (NNLO) nucleon-nucleon interaction and the very same interaction for a consistent full-folding calculation of the effective (optical) potential for nucleon-nucleus scattering for medium-heavy nuclei. The leading order effective (optical) folding potential is computed by integrating over a translationally invariant SA-NCSM one-body scalar density, spin-projected momentum distribution, and the Wolfenstein amplitudes 𝐴, 𝐶, and 𝑀. The resulting nonlocal potentials serve as input for a momentum space Lippmann-Schwinger equation. In the SA-NCSM, the model space is systematically up-selected using Sp⁡(3,ℝ) symmetry considerations. For the light nucleus of 6 He, we establish a systematic selection scheme in the SA-NCSM for scattering observables. Then, we apply this scheme to calculations of scattering observables, such as differential cross sections, analyzing powers, and spin rotation functions for elastic proton scattering from 20 Ne and 40 Ca in the energy regime between 65 and 200 MeV, and compare to available data. Furthermore, our calculations show that the leading order effective nucleon-nucleus potential in the Watson expansion of multiple scattering theory obtained from an up-selected SA-NCSM model space describes 40 Ca elastic scattering observables reasonably well to about 60 degrees in the center-of-mass frame, which coincides roughly with the validity of the NNLO chiral interaction used to calculate both the nucleon-nucleon amplitudes and the one-body scalar and spin nuclear densities.

Ab initio calculations↗

Performance evaluation of a Terry GS-2 steam impulse turbine with air-water mixtures

Terry steam turbines are widely used in various industries because of their robust design. Within the nuclear power generation industry, they are used in the Reactor Core Isolation Cooling System to remove decay heat during reactor isolation events. During the Fukushima Daiichi nuclear power station disaster in Japan in 2011, the Reactor Core Isolation Cooling System and associated Terry turbine operated for over 70 hours in Unit 2; this runtime is well beyond the expected operating duration. Theories suggest the turbine was subjected to a two-phase inlet flow, which could degrade the turbine performance. In this work, an experimental test rig was constructed to test a full-scale Terry model GS-2 steam turbine under two-phase air/water flows. Steady-state efficiency and torque performance maps of the turbine were developed over a range of turbine inlet pressures (1.38–4.83 bar or 20–70 psia), air mass fractions (0.05–1.0) and rotational speeds up to 4000 RPM. Furthermore, turbine performance followed expected trends with torque varying linearly and efficiency varying quadratically with rotational speed. In addition, high-speed images of the two-phase flow entering the turbine were also analyzed to understand how changes in inlet pressure and air mass fraction affect the flow regime and homogenization. The present tests with air–water two-phase mixtures are an important step towards providing an understanding of the full-scale Terry turbine’s behavior and performance curves under two-phase conditions. The results of this work will be combined with air/water and steam/water data gathered using a small-scale Terry ZS-1 steam turbine in order to understand the scaling relationship between large and small size Terry turbines and fluid pairs. The combined data set will enable further development of analytical models over a wide range of conditions and may be used to provide technical justification for expanded use of the Terry turbines in nuclear power plant safety systems and other systems.

42 ENGINEERING↗

Survey-wide asteroid discovery with a high-performance computing enabled non-linear digital tracking framework

Modern astronomical surveys detect asteroids by linking together their appearances across multiple images taken over time. This approach faces limitations in detecting faint asteroids and handling the computational complexity of trajectory linking. Here, we present a novel method that adapts “digital tracking” – traditionally used for short-term linear asteroid motion across images – to work with large-scale synoptic surveys such as the Vera Rubin Observatory Legacy Survey of Space and Time (Rubin/LSST). Our approach combines hundreds of sparse observations of individual asteroids across their non-linear orbital paths to enhance detection sensitivity by several magnitudes. To address the computational challenges of processing massive data sets and dense orbital phase spaces, we developed a specialized high-performance computing architecture. We demonstrate the effectiveness of our method through experiments that take advantage of the extensive computational resources at Lawrence Livermore National Laboratory. This work enables the detection of significantly fainter asteroids in existing and future survey data, potentially increasing the observable asteroid population by orders of magnitude across different orbital families, from near-Earth objects (NEOs) to Kuiper belt objects (KBOs).

Asteroid discovery↗

Analysis of HolosGen Sub-Scale Simulator with Plant Dynamics Code

The Subcritical Power Module Sub-scale Simulator (SPM-SS) has been designed and constructed by HolosGen LLC under the ARPA-E MEITNER program to simulate the thermal-hydraulic and heat transfer behavior of the full-scale Holos-Quad Subcritical Power Modules (SPMs). Four coupled SPMs, each rated at 5.5MW, form the Holos-Quad gas-cooled microreactor design. The SPM-SS represents a substantially scaled-down system with a power rating less than 40 kW, equipped with an electrically heated fuel cartridge heat exchanger, an electrically heated compressor heat exchanger, and a valve actuated turbine heat exchanger, in addition to a recuperator and a cooler heat exchanger. The fuel cartridge represents a portion of the full-scale SPM core, the compressor heat exchanger mimics the temperature changes resulting from the compressor’s turbomachinery inefficiencies, the turbine heat exchanger mimics the expansion process normally occurring through the turbine, while the recuperator and cooler heat exchangers complete the subscale simulator loop. The heaters equipping the fuel cartridge and the compressor heat exchangers are electronically controlled to simulate normal and off-normal SPM operating conditions. The full-scale Holos-Quad SPM design eliminates the traditional balance of plant and executes thermal-to-electric energy conversion by means of an intercooled Brayton cycle with decoupled compressor-turbine turbomachinery. The Holos-Quad full-scale design is equipped with a multi-stage axial Low- and High-Pressure compressor, and a multistage axial turbine. The SPM-SS is designed for testing and validation of selected components which are instead coupled by a traditional balance of plant. The SPMSS is not equipped with turbo-machinery (compressor and turbine) as the development of these components were excluded from the scope of work under the ARPA-E MEITNER funding program. The SPM-SS balance of plant enables modifications, replacement and testing of individual components with different working fluids and is designed to include the turbo-machinery components that will be developed in future research . The SPM-SS can be operated with different gases, variable mass-flow-rates, pressures, and temperatures to obtain test data for selected components, whose performance can be scaled to validate the computer model of the full-scale SPM at various conditions (e.g., start-up, transients conditions). The SPM-SS can operate at the maximum Holos-Quad design pressure of 7 MPa, and a maximum temperature limited to 650 °C by the electrical heaters. Several SPM-SS tests have been conducted and analyzed with the Plant Dynamics Code (PDC) developed at the Argonne National Laboratory (ANL). These tests aimed at validating the PDC modeled predictions of the full-scale Holos-Quad design with data from selected SPM-SS components. In order to address SPM-SS specific characteristics, such as components heat losses and absence of turbomachinery components, some modifications to the PDC have been implemented to factor the design differences from the full-scale Holos-Quad SPM to the SPM-SS. As the PDC offers capabilities to analyze systems with different working fluids, air, nitrogen, and helium were utilized as the SPM-SS working fluids. Air was utilized to fine-tune the SPM-SS Systems Structures and Components (SSCs), nitrogen was utilized to pressure test the SPM-SS loop at the SPMs design maximum pressure of 7MPa. Helium was utilized as the working fluid circulating through the SPM-SS SSCs for specific tests to validate the PDC predictions of the fuel cartridge heat exchanger. SPM-SS tests data were also analyzed with both the steady-state and transient analysis capabilities offered by the PDC. This report describes the PDC analysis of the SPM-SS tests data, including the necessary code modifications and comparison of the code results with the experimental data. Based on the results, a discussion is presented on how the analysis supports design and transient calculations of the full-scale Holos-Quad microreactor. Also based on the results of this work, recommendations are made for future optimizations of the SPM-SS components and PDC model development needs.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Labeling sequential data from noisy annotations

Crowdsourcing algorithms often work under the assumption that the data samples are independent. Recent work has shown that data dependence, such as temporal correlations in sequential data, can be leveraged to improve the label quality. Existing methods that exploit this special structure rely on third-order statistics of the annotator outputs to ensure the identifiability of key latent parameters, which are costly to acquire. This work proposes an approach for integrating crowdsourced annotations under the Dawid-Skene/Hidden Markov Model (DS-HMM) for sequential data based on second-order statistics, which naturally enjoys a lower sample complexity. An effective algorithm is proposed to tackle the challenging optimization problem associated with the proposed estimator. Numerical experiments showcase the effectiveness of the data labeling paradigm.

Marrinan, Timothy P.↗

The Installation of Direct Water-Cooling Systems to Reduce Cooling Energy Requirements for High-Performance Computing Centers

A large cluster of High-Performance Computing (HPC) equipment at the Lawrence Livermore National Laboratory in California was retrofitted with an Asetek cooling system. The Asetek system is a hybrid scheme with water-cooled cold plates on high-heat-producing components in the information technology (IT) equipment, and with the remainder of the heat being removed by conventional air-cooling systems. In order to determine energy savings of the Asetek system, data were gathered and analyzed two ways: using top-down statistical models, and bottom-up engineering models. The cluster, “Cabernet”, rejected its heat into a facilities cooling water loop which in turn rejected the heat into the same chilled water system serving the computer-room air handlers (CRAHs) that provided the air-based cooling for the room. Because the “before” and “after” cases both reject their heat into the chilled water system, the only savings is due to reduction in CRAH fan power. The top-down analysis showed a 4% overall energy savings for the data center (power usage effectiveness (PUE) —the ratio of total data center energy to IT energy— dropped from 1.60 to 1.53, lower is better); the bottom-up analysis showed a 3% overall energy savings (PUE from 1.70 to 1.66) and an 11% savings for the Cabernet system by itself (partial PUE of 1.51). Greater savings, on the order of 15-20%, would be possible if the chilled water system was not used for rejecting the heat from the Asetek system. About 37% of the heat from the Cab system was rejected to the cooling water, lower than at other installations.

Earni, Shankar↗

Intercomparison of Aerosol Instrumental Responses Using AEROTRAK, SMPS, and APS in a Simulated Building

In this work, three aerosol instruments were intercompared by using data collected from field campaigns in a simulated building environment during the last quarter of 2020 from October to December at Oak Ridge National Laboratory. The three aerosol instruments—AeroTrak Model 9110, Scanning Mobility Particle Spectrometer (SMPS) 3080L + 3025A, and Aerodynamic Particle Sizer (APS) 3320—were manufactured by TSI Inc. to perform in situ measurements of particle number concentration in the air. These instruments were co-located in one room in which an aerosol generator was used to produce a defined aerosol population for the instruments to sample and detect in each campaign. The SMPS measured aerosol particles with diameters from 0.01 to 0.7 µm, APS measured aerosol particles with diameters from 0.5 to 20 µm, and AeroTrak measured aerosol particles with diameters from 0.1 to 10 µm. The results suggest that although AeroTrak is fast and easy to use, it does not provide measurements comparable with those obtained by SMPS and APS. The differences in particle counts were orders of magnitude in size (i.e., the AeroTrak data were unreasonably small in the tested scenarios). The results clearly indicate that AeroTrak should not be used in ambient applications in which the level of particle concentration is expected to be high. This would cause difficulty to a campaign in which the dilution factor is unknow a priori. The high particle concentration could also lead to an erroneous aerosol size distribution (e.g., a bimodal size distribution was measured as opposed to a single-mode one). The number concentration limit that AeroTrak is expected to run reliably (i.e., with a counting coincidence error less than ±10%) is approximately three to four per cubic centimeter. This suggests that the instrument is best used in a clean room environment, not in an ambient or indoor application, such as the Flexible Research Platform campaigns.

47 OTHER INSTRUMENTATION↗

Macroscopic Traffic Modeling Using Probe Vehicle Data: A Machine Learning Approach

Abstract The macroscopic fundamental diagram (MFD) captures an orderly relationship among traffic flow, density, and speed at the network level. It is a simple yet powerful tool for modeling traffic dynamics in large urban networks with broad application in traffic control and management. However, empirically derived MFDs in urban regions require high-resolution traffic data from the network. Having the network flow and vehicular density estimated at the (granular) census tract level using vehicle probe data, we apply machine learning methods to predict the MFDs across U.S. urban areas and capture the impacts of location-specific input features on the network flow–density relationships at a large scale. The results show that, among the four tested machine learning approaches (Random Forest, XGBoost, Support Vector Machine, and Neural Network), XGBoost delivers the best performance in predicting network traffic flow based on vehicular density and location attributes. Using interaction Shapley Additive explanation (SHAP) values and partial correlation analysis, we examine the factors influencing MFD shapes across different locations. Our empirical findings reveal that across U.S. urban areas, network topology, transportation infrastructure, and land use are primary factors shaping MFD curves, while demand and trip-related factors play a lesser role. Specifically, higher ranking roads, centrality, and development levels correlate positively with network capacity and critical density, whereas negative associations are observed for network connectivity, mixed-use development, and road roughness levels.

Jin, Ling↗

Hutchinson Trace Estimation for high-dimensional and high-order Physics-Informed Neural Networks

Physics-Informed Neural Networks (PINNs) have proven effective in solving partial differential equations (PDEs), especially when some data are available by seamlessly blending data and physics. However, extending PINNs to high-dimensional and even high-order PDEs encounters significant challenges due to the computational cost associated with automatic differentiation in the residual loss function calculation. Herein, we address the limitations of PINNs in handling high-dimensional and high-order PDEs by introducing the Hutchinson Trace Estimation (HTE) method. Starting with the second-order high-dimensional PDEs, which are ubiquitous in scientific computing, HTE is applied to transform the calculation of the entire Hessian matrix into a Hessian vector product (HVP). This approach not only alleviates the computational bottleneck via Taylor-mode automatic differentiation but also significantly reduces memory consumption from the Hessian matrix to an HVP’s scalar output. We further showcase HTE’s convergence to the original PINN loss and its unbiased behavior under specific conditions. Comparisons with the Stochastic Dimension Gradient Descent (SDGD) highlight the distinct advantages of HTE, particularly in scenarios with significant variability and variance among dimensions. We further extend the application of HTE to higher-order and higher-dimensional PDEs, specifically addressing the biharmonic equation. By employing tensor-vector products (TVP), HTE efficiently computes the colossal tensor associated with the fourth-order high-dimensional biharmonic equation, saving memory and enabling rapid computation. The effectiveness of HTE is illustrated through experimental setups, demonstrating comparable convergence rates with SDGD under memory and speed constraints. Additionally, HTE proves valuable in accelerating the Gradient-Enhanced PINN (gPINN) version as well as the Biharmonic equation. Overall, HTE opens up a new capability in scientific machine learning for tackling high-order and high-dimensional PDEs.

Curse of dimensionality↗