Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Sensitivity Simulation”

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 1,693 records · Page 94
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Novel usage of deep learning and high-performance computing in long-baseline neutrino oscillation experiments

Mención Internacional en el título de doctorDeep-learning methods are playing a crucial role in numerous scientific and industrialapplications. Over the past two decades, these techniques have helped in the collection,reconstruction, and analysis of large data samples in particle physics experiments. Themain topic of this PhD research is the study of deep-learning techniques in long-baselineneutrino oscillation experiments. Neutrinos are mysterious light elementary particles,and their investigation is essential to shed light on some of the remaining open questionsin physics. The work presented here describes an algorithm based on a convolutionalneural network developed to provide highly accurate and efficient selections of electronneutrino and muon neutrino interactions in the Deep Underground Neutrino Experiment(DUNE). With this algorithm, the electron neutrino (antineutrino) selection efficiencypeaks at 90% (94%) and exceeds 85% (90%) for reconstructed neutrino energies between2-5 GeV. The selection efficiency for muon neutrino (antineutrino) interactions is foundto have a maximum of 96% (97%) and exceeds 90% (95%) efficiency for reconstructedneutrino energies above 2 GeV. When considering all electron neutrino and antineutrinointeractions as signal (both those appearing from oscillations and those intrinsic tothe beam), a selection purity of 90% is achieved. These event selections are criticalto maximise the sensitivity of the experiment to CP-violating effects, key to furtherunderstand the matter-antimatter asymmetry of the Universe.In high-energy physics experiments, deep learning has also been explored for producingfast simulations and physically-motivated manipulations of simulated images. Some ofthose simulations, such as the light production and detection, are very computationallyexpensive and require novel methods to produce the necessary samples while controllingthe varied underlying physics model parameters. To do so, we invented the model-assistedgenerative adversarial network (MAGAN), first validated on simple generic case studiesand then successfully applied to the DUNE photon-detector simulation.Moreover, we also developed graph neural networks for 3D-voxel classification ofambiguities and optical crosstalk for a different particle physics experiment, most preciselyfor the proposed SuperFGD. This novel 3D-granular plastic-scintillator neutrino detectorwill be used to upgrade the near detector of the T2K neutrino oscillation experiment, and our method reports efficiencies and purities of 94-96% per event in the classificationof particle track voxels.Due to the growth and complexity of deep neural networks, researchers have beeninvestigating techniques to train those networks in a more computationally-efficient way.Many efforts have been made by the community to optimise deep-learning models byparallelising or distributing their training computation across multiple devices. In thisthesis, we study an approach based on data locality for those neural networks that cannotbenefit from scaling their computation due to a significant bottleneck in the data I/O.The research also includes a detailed study on the performance of deep neural networkson hardware accelerator boards.Los métodos de aprendizaje profundo son cada vez más utilizados en numerosas aplicacionescientíficas e industriales hoy en día. Durante las dos últimas décadas, estastécnicas se han empleado en la recolección, reconstrucción y análisis de la gran cantidadde datos generados por experimentos de física de partículas. El tema principal de estatesis doctoral es el uso de estos modelos de aprendizaje profundo en experimentos defísica de neutrinos, en concreto en los experimentos de larga distancia DUNE y T2K. Losneutrinos, partículas fundamentales neutras, de las más ligeras del Universo, pueden serclave para explicar algunas de las cuestiones todavía sin resolver en física fundamental.Entre las diferentes contribuciones que esta tesis ha hecho a su estudio, cabe destacar eldesarrollo de un algoritmo basado en una red de neuronas convolucional para seleccionarcon gran eficiencia y precisión las interacciones de neutrinos electrónicos y muónicos enel Deep Underground Neutrino Experiment (DUNE). La eficiencia de selección obtenidapara neutrinos (antineutrinos) electrónicos alcanza un máximo del 90% (94%) y supera el85% (90%) para neutrinos con energías reconstruidas en el rango 2-5 GeV. La selección deneutrinos (antineutrinos) muónicos tiene una eficiencia máxima del 96% (97%) y excedeel 90% (95%) para neutrinos con energías reconstruidas de más de 2 GeV. Considerandocomo señal todas las interacciones de neutrinos y antineutrinos electrónicos (procedentestanto de oscilaciones como intrínsecos en el haz inicial), se logra una pureza en la seleccióndel 90%. Dichas selecciones de eventos son fundamentales para maximizar la sensibilidaddel experimento a los efectos de violació...

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Charged particle transport coefficient challenges in high energy density plasmas

High energy density physics (HEDP) and inertial confinement fusion (ICF) research typically relies on computational modeling using radiation-hydrodynamics codes in order to design experiments and understand their results. These tools, in turn, rely on numerous charged particle transport and relaxation coefficients to account for laser energy absorption, viscous dissipation, mass transport, thermal conduction, electrical conduction, non-local ion (including charged fusion product) transport, non-local electron transport, magnetohydrodynamics, multi-ion-species thermalization, and electron-ion equilibration. In many situations, these coefficients couple to other physics, such as imposed or self-generated magnetic fields. Furthermore, how these coefficients combine are sensitive to plasma conditions as well as how materials are distributed within a computational cell. Uncertainties in these coefficients and how they couple to other physics could explain many of the discrepancies between simulation predictions and experimental results that persist in even the most detailed calculations. This paper reviews the challenges faced by radiation-hydrodynamics in predicting the results of HEDP and ICF experiments with regard to these and other physics models typically included in simulation codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Solving the Hele–Shaw flow using the Harrow–Hassidim–Lloyd algorithm on superconducting devices: A study of efficiency and challenges

The development of quantum processors for practical fluid flow problems is a promising yet distant goal. Recent advances in quantum linear solvers have highlighted their potential for classical fluid dynamics. In this study, we evaluate the Harrow–Hassidim–Lloyd (HHL) quantum linear systems algorithm (QLSA) for solving the idealized Hele–Shaw flow. Our focus is on the accuracy and computational cost of the HHL solver, which we find to be sensitive to the condition number, scaling exponentially with problem size. This emphasizes the need for preconditioning to enhance the practical use of QLSAs in fluid flow applications. Moreover, we perform shots-based simulations on quantum simulators and test the HHL solver on superconducting quantum devices, where noise, large circuit depths, and gate errors limit performance. Error suppression and mitigation techniques improve accuracy, suggesting that such fluid flow problems can benchmark noise mitigation efforts. Finally, our findings provide a foundation for future, more complex application of QLSAs in fluid flow simulations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

High yield polar direct drive fusion neutron sources at the National Ignition Facility

Polar direct drive neutron source experiments were performed at the National Ignition Facility showing substantial improvement in total neutron yield and efficiency of conversion of laser energy to fusion output. Plastic capsules 3–4 mm in diameter were filled with 1.5 mg/cc of deuterium–tritium (DT) fuel and imploded with laser beam pointing and defocus designed to compensate for polar asymmetry introduced by the facility beam entrance angles. Radiation-hydrodynamics simulations were employed to optimize the multi-dimensional laser and target parameter space, within facility and target fabrication constraints. Ensembles of 1D simulations tuned to match the outputs of early shots in the series were used to design subsequent shots in the series. This allowed the later shots to be designed based on empirically motivated sensitivities to laser and target input parameters, while eliminating the need to explicitly model phenomena such as hydrodynamic instabilities and nonlinear laser–plasma interactions. Additionally, one experiment with a 3.0 mm diameter CH capsule produced 13.6 kJ (4.81 × 10 15 DT neutrons) from a laser input below the NIF optics damage threshold at 585 kJ, 328 TW. Two experiments with 4.0 mm capsules produced 31.3 and 33.6 kJ of fusion output (1.11 × 10 16 and 1.19 × 10 16 DT neutrons) with 1.10 MJ, 390 TW and 1.26 MJ, 425 TW of laser input, respectively.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Formation of active sites on transition metals through reaction-driven migration of surface atoms

Adopting low-index single-crystal surfaces as models for metal nanoparticle catalysts has been questioned by the experimental findings of adsorbate-induced formation of subnanometer clusters on several single-crystal surfaces. Here we used density functional theory calculations to elucidate the conditions that lead to cluster formation and show how adatom formation energies enable efficient screening of the conditions required for adsorbate-induced cluster formation. We studied a combination of eight face-centered cubic transition metals and 18 common surface intermediates and identified systems relevant to catalytic reactions, such as carbon monoxide (CO) oxidation and ammonia (NH 3 ) oxidation. We used kinetic Monte Carlo simulations to elucidate the CO-induced cluster formation process on a copper surface. Scanning tunneling microscopy of CO on a nickel (111) surface that contains steps and dislocations points to the structure sensitivity of this phenomenon. Metal-metal bond breaking that leads to the evolution of catalyst structures under realistic reaction conditions occurs much more broadly than previously thought.

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

A study of nuclear structure of light nuclei at the electron–ion collider

Understanding the substructure of atomic nuclei, particularly the clustering of nucleons inside them, is essential for comprehending nuclear dynamics. Various cluster configurations can emerge depending on excitation energy, the number and types of core clusters, and the presence of excess neutrons. Despite the prevalence of tightly bound cluster formations in low-lying states, understanding the correlation between clusters and their formation mechanisms remains incomplete. This exploring study investigates nuclear clustering at the electron–ion collider (EIC) using simulations based on the modified BeAGLE model. By simulating collisions involving e+ 9 Be, e+ 12 C, and e+ 16 O nuclei, we find that the average energy of particles and the system size ratios of particles at forward rapidity exhibit sensitivity to alpha clustering and its various configurations. These findings offer valuable insights into the dynamics of nuclear clustering and its implications for future studies at the EIC.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Understanding the Roles of Convective Trigger Functions in the Diurnal Cycle of Precipitation in the NCAR CAM5

The wrong diurnal cycle of precipitation is a common weakness of current global climate models (GCMs). To improve the simulation of the diurnal cycle of precipitation and understand what physical processes control it, we test a convective trigger function described in Xie et al. (2019) with additional optimizations in the NCAR Community Atmosphere Model version 5 (CAM5). The revised trigger function consists of three modifications: 1) replacing the Convective Available Potential Energy (CAPE) trigger with a dynamic CAPE (dCAPE) trigger, 2) allowing convection to originate above the top of planetary boundary layer (i.e., the unrestricted air parcel launch level - ULL), and 3) optimizing the entrainment rate and threshold value of the dynamic CAPE generation rate for convection onset based on observations. Results from 1°-resolution simulations show that the revised trigger can alleviate the long-standing GCM problem of too early maximum precipitation during the day and missing the nocturnal precipitation peak that is observed in many regions, including the U.S. southern Great Plains (SGP). The revised trigger also improves the simulation of the propagation of precipitation systems downstream of the Rockies and the Amazon region. A further composite analysis over the SGP unravels the mechanisms through which the revised trigger affects convection. Additional sensitivity tests show that both the peak time and the amplitude of the diurnal cycle of precipitation are sensitive to the entrainment rate and dCAPE threshold values.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

The Impact of Constrained Data Assimilation on the Forecasts of Three Convection Systems During the ARM MC3E Field Campaign

A constrained data assimilation (CDA) system based on the ensemble variational (EnVar) method and physical constraints of mass and water conservations is evaluated through three convective cases during the Midlatitude Continental Convective Clouds Experiment (MC3E) of the Atmospheric Radiation Measurement (ARM) program. Compared to the original data assimilation (ODA), the CDA is shown to perform better in the forecasted state variables and simulated precipitation. The CDA is also shown to greatly mitigate the loss of forecast skills in observation denial experiments when radar radial winds are withheld in the assimilation. Modifications to the algorithm and sensitivities of the CDA to the calculation of the time tendencies in the constraints are described.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Radiation-Hard Smart-Pixel Detector ASIC ReadOut with Digital AI in 28nm

Detectors at future high energy colliders will face enormous technical challenges. Disentangling the unprecedented numbers of particles expected in each event will require highly granular silicon pixel detectors with billions of readout channels. With event rates as high as 40 MHz, these detectors will generate petabytes of data per second. To enable discovery within strict bandwidth and latency constraints, future trackers must be capable of fast, power efficient, and radiation hard data-reduction at the source. This effort is pursuing the co-design development of high-performance readout smart pixel ASICs for a future Phase III High Luminosity upgrade of the Large Hadron Collider. A 1.6mm2 ASIC prototype was designed by Fermilab in CMOS 28 nm bulk process and submitted for manufacturing in February 2024. It leverages the analog front-end pixel design of a previous prototype fabricated and tested in 2023, which achieved a simulated detection level of ~400e- with 30fF input capacitance. The ROIC consists of two matrices of 16×16 smart pixels, each 25×25 μm2 in size. Each smart pixel contains a charge-sensitive preamplifier with leakage current compensation and three auto-zero comparators for a 2-bit flash-type ADC. There is digital space for the integration of our fully combinatorial AI that performs momentum classification at the bunch crossing rate. The total power consumption is ∼6μW per pixel, which corresponds to ~1mW/cm2. The ASIC incorporates programmable front-end charge injection circuitry to generate pixel cluster charges during characterization. The cluster profile will be generated to duplicate hit characteristics of various momentum (pT). We will present early results from chip testing.

Parpillon, Benjamin↗
DOE OSTI2026

Electrochemical testing of metals for petroleum wastewater monitoring

The purpose of this project was to determine the feasibility of ion-selective electrode (ISE) as compared to x-ray fluorescence (XRF) for determining the concentrations of metal ions in petroleum wastewater. ISEs for silver and cadmium were tested in a laboratory under increasingly realistic conditions, then deployed in the sewer wastewater monitoring complex to simulate a petroleum wastewater system. By verifying the ISE calibration after wastewater exposure and comparing its measurements to the XRF systems within the wastewater monitoring system, we determined that the ISE is more sensitive but less feasible for use in live wastewater monitoring due to fouling concerns.

Energy - Petroleum↗
DOE OSTI2023

High-energy Radiation and Ion Acceleration in Three-dimensional Relativistic Magnetic Reconnection with Strong Synchrotron Cooling

Abstract We present the results of 3D particle-in-cell simulations that explore relativistic magnetic reconnection in pair plasma with strong synchrotron cooling and a small mass fraction of nonradiating ions. Our results demonstrate that the structure of the current sheet is highly sensitive to the dynamic efficiency of radiative cooling. Specifically, stronger cooling leads to more significant compression of the plasma and magnetic field within the plasmoids. We demonstrate that ions can be efficiently accelerated to energies exceeding the plasma magnetization parameter, ≫ σ , and form a hard power-law energy distribution, f i ∝ γ −1 . This conclusion implies a highly efficient proton acceleration in the magnetospheres of young pulsars. Conversely, the energies of pairs are limited to either σ in the strong cooling regime or the radiation burnoff limit, γ syn , when cooling is weak. We find that the high-energy radiation from pairs above the synchrotron burnoff limit, ε c ≈ 16 MeV, is only efficiently produced in the strong cooling regime, γ syn < σ . In this regime, we find that the spectral cutoff scales as ε cut ≈ ε c ( σ / γ syn ) and the highest energy photons are beamed along the direction of the upstream magnetic field, consistent with the phenomenological models of gamma-ray emission from young pulsars. Furthermore, our results place constraints on the reconnection-driven models of gamma-ray flares in the Crab Nebula.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Kinetics of Carbon Condensation in Detonation of High Explosives: First-Order Phase Transition Theory Perspective

We report the kinetics of carbon condensation, or carbon clustering, in detonation of carbon-rich high explosives is modeled by solving a system of rate equations for concentrations of carbon particles. Unlike previous efforts, the rate equations account not only for the aggregation of particles but also for their fragmentation in a thermodynamically consistent manner. Numerical simulations are performed, yielding the distribution of particle concentrations as a function of time. In addition to that, analytical expressions are obtained for all the distinct steps and regimes of the condensation kinetics, which facilitates the analysis of the numerical results and allows one to study the sensitivity of the kinetic behavior to the variation of system parameters. The latter is important because the numerical values of many parameters are not reliably known at present. The theory of the kinetics of first-order phase transitions is found adequate to describe the general kinetic trends of carbon condensation, as described by the rate equations. Such physical phenomena and processes as the coagulation, nucleation, growth, and Ostwald ripening are observed, and their dependence on various system parameters is studied and reported. It is believed that the present work will become useful when analyzing the present and future results for the kinetics of carbon condensation, obtained from experiments or atomistic simulations.

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

Real-time signal detection for Cyclotron Radiation Emission Spectroscopy measurements using antenna arrays

Cyclotron Radiation Emission Spectroscopy (CRES) is a technique for precision measurement of the energies of charged particles, which is being developed by the Project 8 Collaboration to measure the neutrino mass using tritium beta-decay spectroscopy. Project 8 seeks to use the CRES technique to measure the neutrino mass with a sensitivity of 40 meV, requiring a large supply of tritium atoms stored in a multi-cubic meter detector volume. Antenna arrays are one potential technology compatible with an experiment of this scale, but the capability of an antenna-based CRES experiment to measure the neutrino mass depends on the efficiency of the signal detection algorithms. Here, in this paper, we develop efficiency models for three signal detection algorithms and compare them using simulations from a prototype antenna-based CRES experiment as a case-study. The algorithms include a power threshold, a matched filter template bank, and a neural network based machine learning approach, which are analyzed in terms of their average detection efficiency and relative computational cost. It is found that significant improvements in detection efficiency and, therefore, neutrino mass sensitivity are achievable, with only a moderate increase in computation cost, by utilizing either the matched filter or machine learning approach in place of a power threshold, which is the baseline signal detection algorithm used in previous CRES experiments by Project 8.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

High-Performance Simulation of Dynamic Hydrologic Exchange and Implications for Surrogate Flow and Reactive Transport Modeling in a Large River Corridor

Hydrologic exchange flows (HEFs) have environmental significance in riverine ecosystems. Key river channel factors that influence the spatial and temporal variations of HEFs include river stage, riverbed morphology, and riverbed hydraulic conductivity. However, their impacts on HEFs were often evaluated independently or on small scales. In this study, we numerically evaluated the combined interactions of these factors on HEFs using a high-performance simulator, PFLOTRAN, for subsurface flow and transport. The model covers 51 square kilometers of a selected river corridor with large sinuosity along the Hanford Reach of the Columbia River in Washington, US. Three years of spatially distributed hourly river stages were applied to the riverbed. Compared to the simulation when riverbed heterogeneity is not ignored, the simulation using homogeneous riverbed conductivity underestimated HEFs, especially upwelling from lateral features, and overestimated the mean residence times derived from particle tracking. To derive a surrogate model for the river corridor, we amended the widely used transient storage model (TSM) for riverine solute study at reach scale with reactions. By treating the whole river corridor as a batch reactor, the temporal changes in the exchange rate coefficient for the TSM were derived from the dynamic residence time estimated from the hourly PFLOTRAN results. The TSM results were evaluated against the effective concentrations in the hyporheic zone calculated from the PFLOTRAN simulations. Our results show that there is potential to parameterize surrogate models such as TSM amended with biogeochemical reactions while incorporating small-scale process understandings and the signature of time-varying streamflow to advance the mechanistic understanding of river corridor processes at reach to watershed scales. However, the assumption of a well-mixed storage zone for TSM should be revisited when redox-sensitive reactions in the storage zones play important roles in river corridor functioning.

Fang, Yilin↗
DOE OSTI2020

Technology Summary DNATRAX

DNATrax (DNA Tagged Reagents for Aerosol eXperiments) was developed to provide a safe simulant for understanding the transport and dispersion of pathogens. Specifically, it was developed as a safe surrogate for the Bacillus anthrasis spore bioagent, since most naturally occurring spores are not safe for public release. As DNA-tagged sugar particles, DNATrax aerosols are safe for public release. Scientists can add unique DNA barcode sequences to the sugar particles to produce a wide variety of test particles that can be easily differentiated, making it possible to release and test them simultaneously. The unique barcodes also eliminate the need to conduct decontamination between tests, dramatically reducing the cost and time to run experiments. Additionally, the surrogate offers a high-sensitivity detection method to understand key public health issues. We conducted both lab-based and largescale field studies using DNATrax to validate its capabilities for use as a pathogen surrogate. We have demonstrated the ability to produce DNATrax in a range of sizes, adjusting the microparticles to simulate different sized pathogens with ability to create unique sequences of DNA bases, enabling them to create a nearly unlimited variety of test particles.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

PPPL Report on Reduced Modeling of Fusion Alpha Transport in ARC Burning Plasmas

We are reporting on the modeling of fusion alpha particle transport in the planned ARC fusion device being designed by the CFS (Commonwealth Fusion Systems: https://cfs.energy). The ARC tokamak is designed to operate in a burning-plasma regime characterized by a substantial population of fusion-born alpha particles. Alfvén eigenmode (AE) stability is assessed both analytically and numerically, incorporating alpha-particle drive, ion Landau, and radiative damping from thermal species and collisional damping from trapped electrons. Regions of unstable and near-threshold AE activity are mapped across ARC’s operational parameter space. Linear stability analysis with NOVA indicates multiple, often marginally unstable AEs, extending to toroidal mode numbers up to n= 30. The present report focuses on the ARC flat-top operating point prior to the sawtooth event. Alpha-particle transport on timescales exceeding the neoclassical slowing-down time is assessed using the NUBEAM module [1][2] of the TRANSP code [3], employing transport coefficients derived from the RBQ quasilinear modeling (cf. Appendix B). These global simulations identify favorable and unfavorable operating regimes with respect to alpha confinement, pressure redistribution, and overall alpha-heating efficiency. We also evaluate additional transport mechanisms—including neoclassical tearing mode (TM)–induced stochasticity, sawtooth-driven redistribution, and toroidal-field ripple using the kick model (cf. Appendix C) which makes use of the guiding-center code ORBIT, see Section 5. The kick model is integrated into TRANSP to enable self-consistent predictions of alpha-driven current formation and sustainment within the ARC scenario. Sensitivity scans are performed over the mode frequency, rational-surface alignment, island width, mode amplitude, and proximity of the limiter to the plasma. Our study provides an initial, physics-based guidance for machine design, operational planning, and equilibrium control, ensuring adequate alpha confinement and robust self-heating performance in ARC. Our simulations mostly targeted worst case scenarios, e.g. for TMs and sawteeth. Overall, we expect benign effects for the ARC scenario investigated in this work on fusion alpha confinement and losses in the presence of AEs, tearing modes and sawteeth. This report addresses three thrusts identified at the outset. The first thrust focuses on analytic estimates of the parametric dependencies of EP relaxation based on local AE stability simulations (Section 3). The second thrust involves global evaluations of AE stability using the NOVA, RBQ, and NUBEAM codes (Section 4). Finally, we investigate alpha-particle transport driven by low-frequency instabilities associated with sawteeth and tearing modes (Section 5).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Investigating permafrost carbon dynamics in Alaska with artificial intelligence

Abstract Positive feedbacks between permafrost degradation and the release of soil carbon into the atmosphere impact land–atmosphere interactions, disrupt the global carbon cycle, and accelerate climate change. The widespread distribution of thawing permafrost is causing a cascade of geophysical and biochemical disturbances with global impacts. Currently, few earth system models account for permafrost carbon feedback (PCF) mechanisms. This research study integrates artificial intelligence (AI) tools and information derived from field-scale surveys across the tundra and boreal landscapes in Alaska. We identify and interpret the permafrost carbon cycling links and feedback sensitivities with GeoCryoAI, a hybridized multimodal deep learning (DL) architecture of stacked convolutionally layered, memory-encoded recurrent neural networks (NN). This framework integratesin-situmeasurements and flux tower observations for teacher forcing and model training. Preliminary experiments to quantify, validate, and forecast permafrost degradation and carbon efflux across Alaska demonstrate the fidelity of this data-driven architecture. More specifically, GeoCryoAI logs the ecological memory and effectively learns covariate dynamics while demonstrating an aptitude to simulate and forecast PCF dynamics—active layer thickness (ALT), carbon dioxide flux (CO 2 ), and methane flux (CH 4 )—with high precision and minimal loss (i.e. ALT RMSE : 1.327 cm [1969–2022]; CO 2 RMSE : 0.697µmolCO 2 m −2 s −1 [2003–2021]; CH 4 RMSE : 0.715 nmolCH 4 m −2 s −1 [2011–2022]). ALT variability is a sensitive harbinger of change, a unique signal characterizing the PCF, and our model is the first characterization of these dynamics across space and time.

Environmental Sciences & Ecology↗
DOE OSTI2020

Comparison of momentum transport models for numerical relativity

The main problems of nonvacuum numerical relativity, compact binary mergers and stellar collapse, involve hydromagnetic instabilities and turbulent flows, so that kinetic energy at small scales leads to mean effects at large scale that drive the secular evolution. Notable among these effects is momentum transport. We investigate two models of this transport effect, a relativistic Navier-Stokes system and a turbulent mean stress model, that are similar to all of the prescriptions that have been attempted to date for treating subgrid effects on binary neutron star mergers and their aftermath. Our investigation involves both stability analysis and numerical experimentation on star and disk systems. We also begin the investigation of the effects of particle and heat transport on postmerger simulations. Here, we find that correct handling of turbulent heating is crucial for avoiding unphysical instabilities. Given such appropriate handling, the evolution of a differentially rotating star and the accretion rate of a disk are reassuringly insensitive to the choice of prescription. However, disk outflows can be sensitive to the choice of method, even for the same effective viscous strength. We also consider the effects of eddy diffusion in the evolution of an accretion disk and show that it can interestingly affect the composition of outflows.

79 ASTRONOMY AND ASTROPHYSICS↗