Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “gravity data”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Basin and Range Investigation for Developing Geothermal Energy: Exploration Data

This data package includes exploration material from the Basin & Range Investigation for Developing Geothermal Energy [in Hidden Systems] project (BRIDGE), which is part of a broader initiative to advance the exploration of hidden geothermal resources in the Basin & Range Province of the western U.S. Data modalities include a helicopter-borne time-domain electromagnetic survey, magnetotellurics, 2-meter temperature measurements, ground-based gravity and legacy aeromagnetic surveys, geochemistry, geologic mapping, LiDAR analysis, 3D models, associated geospatial data, and a bibliography of existing data and references utilized in prospect characterization and conceptual modeling. Key files are in CSV, Geosoft, and Geotools formats. Please refer to READMEs for dataset-specific information. Where applicable, acquisition data and inversion models for a particular prospect or area of interest are organized separately. This BRIDGE data package is the product of a collaboration led by Sandia National Laboratories with partners from Geologica Geothermal Group, Inc., the U.S. Navy Geothermal Program Office, and consultants Steven Sewell (Australis Geoscience Ltd) and William Cumming (Cumming Geoscience). The project's areas of interest (AOIs) are based off priority areas of interest in the southwestern portion of the Nevada Play Fairway map, distribution across tectonic provinces, accessibility, and the project team's extensive experience in the region. AOIs cover about a dozen basins that include unexplored prospects, partially explored prospects, and some developed analogue resources that provide validation cases. Many unexplored and partially explored prospects are on U.S. Department of Defense (DoD) land, though adjacent lands are included as well.

15 GEOTHERMAL ENERGY↗
DOE OSTI2024

Congo Basin Water Balance and Terrestrial Fluxes Inferred From Satellite Observations of the Isotopic Composition of Water Vapor

Large spatio-temporal gradients in the Congo basin vegetation and rainfall are observed. However, its water-balance (evapotranspiration minus precipitation, or ET - P) is typically measured at basin-scales, limited primarily by river-discharge data, spatial resolution of terrestrial water storage measurements, and poorly constrained ET. We use observations of the isotopic composition of water vapor to quantify the spatio-temporal variability of net surface water fluxes across the Congo Basin between 2003 and 2018. These data are calibrated at basin scale using satellite gravity and total Congo river discharge measurements and then used to estimate time-varying ET - P over four quadrants representing the Congo Basin, providing first estimates of this kind for the region. We find that the multi-year record, seasonality, and interannual variability of ET - P from both the isotopes and the gravity/river discharge based estimates are consistent. Additionally, we use precipitation and gravity-based estimates with our water vapor isotope-based ET - P to calculate time and space averaged ET and net river discharge within the Congo Basin. These quadrant-scale moisture flux estimates indicate (a) substantial recycling of moisture in the Congo Basin (temporally and spatially averaged ET/P > 70%), consistent with models and visible light-based ET estimates, and (b) net river outflow is largest in the Western Congo where there are more rivers and higher flow rates. Our results confirm the importance of ET in modulating the Congo water cycle relative to other water sources.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Real-Time Bayesian Inference at Extreme Scale: A Digital Twin for Tsunami Early Warning Applied to the Cascadia Subduction Zone

We present a Bayesian inversion-based digital twin that employs acoustic pressure data from seafloor sensors, along with 3D coupled acoustic–gravity wave equations, to infer earthquake-induced spatiotemporal seafloor motion in real time and forecast tsunami propagation toward coastlines for early warning with quantified uncertainties. Our target is the Cascadia subduction zone, with one billion parameters. Computing the posterior mean alone would require 50 years on a 512 GPU machine. Instead, exploiting the shift invariance of the parameter-to-observable map and devising novel parallel algorithms, we induce a fast offline–online decomposition. The offline component requires just one adjoint wave propagation per sensor; using MFEM, we scale this part of the computation to the full El Capitan system (43,520 GPUs) with 92% weak parallel efficiency. Moreover, given real-time data, the online component exactly solves the Bayesian inverse and forecasting problems in 0.2 seconds on a modest GPU system, a ten-billion-fold speedup.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Snowmass2021 Cosmic Frontier White Paper: Enabling Flagship Dark Energy Experiments to Reach their Full Potential

A new generation of powerful dark energy experiments will open new vistas for cosmology in the next decade. However, these projects cannot reach their utmost potential without data from other telescopes. This white paper focuses in particular on the compelling benefits of ground-based spectroscopic and photometric observations to complement the Vera C. Rubin Observatory, as well as smaller programs in aid of a DESI-2 experiment and CMB-S4. These additional data sets will both improve dark energy constraints from these flagship projects beyond what would possible on their own and open completely new windows into fundamental physics. For example, additional photometry and single-object spectroscopy will provide necessary follow-up information for supernova and strong lensing cosmology, while highly-multiplexed spectroscopy both from smaller facilities over wide fields and from larger facilities over narrower regions of sky will yield more accurate photometric redshift estimates for weak lensing and galaxy clustering measurements from the Rubin Observatory, provide critical spectroscopic host galaxy redshifts for supernova Hubble diagrams, provide improved understanding of limiting astrophysical systematic effects, and enable new measurements that probe the nature of gravity. A common thread is that access to complementary data from a range of telescopes/instruments would have a substantial impact on the rate of advance of dark energy science in the coming years.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Integration of Dynamical System Scaling to RAVEN and Facility Application

As part of the design optimization and model validation effort within the hybrid energy systems program, the research to implement dynamical system scaling (DSS) code to RAVEN is funded under the Integrated Energy System program in collaboration with the Digital Reactor Technology & Development department. The DSS data processing algorithm has been coded within the RAVEN framework along with other metrics and postprocessing models. The implemented code was tested using a gravity-driven draining tank draining case modeled in RELAP5-3D and generated data was successfully postprocessed based on the scaling analyses defined by DSS. Postprocessed data concluded less agreement for the tank exit velocity and indicated that the modeled input deck for RELAP5-3D requires modification to enforce minimal pressure differential effects to the draining process.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Scattering Amplitudes and QFT Insight (Final Scientific Report)

The calculation of scattering amplitudes provides an invariant window into the physical dynamic content of relativistic quantum field theory. Yet with traditional methods, for empirically relevant theories, these calculations scale with a factorial complexity in external particles and precision. If we aspire to collapse the theoretical uncertainty obscuring new physics hidden at all scales from the microscopic probed at high energy colliders to the largest effective field theory in the universe governing the evolution of large scale structure, this complexity challenge necessitates new ideas and methods in calculation. Novel approaches like the color-kinematics duality and the associated double-copy construction have drastically simplified the situation, relating both gauge and gravity theory predictions to a much smaller kernel of invariant kinematic data. This project supporting research by the PI's Amplitudes and Insights group at Northwestern University looked to push insight deep into both the IR and the UV by establishing how novel structures must constrain the predictions of counterterms in both gauge and gravity theories, setting the groundwork to exploring the high energy behavior of particular gravity theories via constituent gauge-theory calculations, and confronting fundamental challenges at the interface between QFT amplitudes analysis and next generation gravitational wave science as well as well as inflationary and large scale structure cosmology.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Planck constraints and gravitational wave forecasts for primordial black hole dark matter seeded by multifield inflation

Here, we perform a Markov Chain Monte Carlo analysis of a simple yet generic multifield inflation model characterized by two scalar fields coupled to each other and nonminimally coupled to gravity, fit to Planck 2018 cosmic microwave background (CMB) data. In particular, model parameters are constrained by data on the amplitude of the primordial power spectrum of scalar curvature perturbations on CMB scales A s , the spectral index n s , and the ratio of power in tensor to scalar modes r, with a prior that the primordial power spectrum should also lead to primordial black hole (PBH) production sufficient to account for the observed dark matter abundance. We find that n s in particular controls the constraints on our model. Whereas previous studies of PBH formation from an ultraslow-roll phase of inflation have highlighted the need for at least one model parameter to be highly fine-tuned, we identify a degeneracy direction in parameter space such that shifts by ~10% of one parameter can be compensated by comparable shifts in other parameters while preserving a close fit between model predictions and observations. Furthermore, we find this allowed parameter region produces observable gravitational wave signals in the frequency ranges to which upcoming experiments are projected to be sensitive, including Advanced LIGO and Virgo, the Einstein Telescope, Cosmic Explorer, DECIGO, and LISA.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Searches for beyond-Riemann gravity

Many effective field theories describing gravity cannot arise from an underlying theory based on Riemann geometry or its extensions to include torsion and nonmetricity but may instead emerge from another geometry or may have a nongeometric basis. The Lagrange density for a broad class of such theories is investigated. The action for fermions coupled to gravity is linearized about a Minkowski background and is found to include terms describing small deviations from Lorentz invariance and gravitational gauge invariance. The corresponding nonrelativistic hamiltonian is derived at second order in the fermion momentum. The implications for laboratory experiments and astrophysical observations with fermions are studied, with primary focus on anomalous spin-gravity couplings. First constraints on some coefficients are extracted from existing data obtained via measurements at different potentials, comparisons of gravitational accelerations, interferometric methods, and investigations of gravitational bound states. As a result, some prospects for future experimental studies are discussed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

New Tests of Low-Scale Quantum Gravity with Cosmic-Ray Collisions

Cosmic ray collisions at high center of mass energy could enable graviton and black hole production as expected in theories of low-scale quantum gravity, such as extra-dimensions, many species, or some versions of string theory. Here we propose three novel phenomenological tests of these theories. We first consider the collision of cosmic rays with ambient protons, electrons and photons in Active Galactic Nuclei (AGN), finding that high-energy neutrino data from the blazar TXS 0506+056 places a constraint on the fundamental scale of gravity of $M_f \gtrsim 0.3$ TeV, and future high-energy neutrino data could raise this bound to $M_f \gtrsim 200$ TeV. We then point out that collisions of pairs of cosmic rays could occur at a sizable rate in AGN where the accelerated cosmic rays are not collimated, or on supermassive black hole binaries. This consideration could potentially let us test unprecedented large fundamental scales of $M_f \gtrsim 2$ PeV. We further compute the corresponding thermal neutrino emission arising from the Hawking evaporation of black holes produced in cosmic ray collisions, finding a spectrum that clearly differs from that expected in meson decays. Finally, we speculate with an scenario which would produce high-energy neutrino and gamma-ray emission from regions in the sky where no multi-wavelength counterparts would be expected, via graviton propagation from a different brane, which then decays in our Universe.

Ettengruber, Manuel [IPhT, Saclay]↗
DOE OSTI2025

Modified gravity constraints from the full shape modeling of clustering measurements from DESI 2024

We present cosmological constraints on deviations from general relativity (GR) from the first-year of clustering observations from the Dark Energy Spectroscopic Instrument (DESI) in combination with other available datasets including the CMB data from Planck with CMB-lensing from Planck and ACT, BBN constraints on the physical baryon density, the galaxy weak lensing and clustering from DESY3 and supernova data from DESY5. We first consider the μ(a,k)–Σ(a,k) modified gravity (MG) parameterization (as well as η(a,k)) in a ΛCDM and a w 0 w a CDM cosmological backgrounds. Using a functional form for time-only evolution gives μ 0 = 0.11 +0.44 -0.54 from DESI(FS+BAO)+BBN and a wide prior on n s . Using DESI(FS+BAO)+CMB+DESY3+DESY5-SN, we obtain μ 0 = 0.05 ± 0.22 and Σ 0 = 0.008 ± 0.045 and similarly μ 0 = 0.02 +0.19 -0.24 and η 0 = 0.09 +0.36 -0.60 , in an ΛCDM background. In w 0 w a CDM we obtain μ 0 = -0.24 +0.32 -0.28 and Σ 0 = 0.006 ± 0.043, consistent with GR, and we still find a preference of the data for a dynamical dark energy with w 0 > -1 and w a < 0. Using functional dependencies in both time and scale gives μ 0 and Σ 0 with a same level of precision as above but other scale MG parameters remain hard to constrain. We then move to binned parameterizations in a ΛCDM background starting with two bins in redshift and obtain, μ 1 = 1.02 ± 0.13, μ 2 = 1.04 ± 0.11, Σ 1 = 1.021 ± 0.029 and Σ 2 = 1.022 +0.027 -0.023 , all consistent with the unity value of GR in the binning formalism. We then extend the analysis to combine two bins in redshift and two in scale giving 8 MG parameters that we find all consistent with GR. We note that we find here that the tension reported in previous studies about Σ 0 being inconsistent with GR when using Planck PR3 data goes away when we use the recent LoLLiPoP+HiLLiPoP likelihoods. As noted in previous studies, this seems to indicate that the tension is indeed related to the CMB lensing anomaly in PR3 which is also resolved when using the recent likelihoods. We then constrain the class of Horndeski theory in the effective field theory of dark energy approach. We consider both EFT-basis and α-basis in the analysis. Assuming a power law parameterization for the EFT function Ω, which controls non-minimal coupling, we obtain Ω 0 = 0.012 +0.001 -0.012 and s 0 = 0.996 +0.54 -0.20 from the combination of DESI(FS+BAO)+DESY5SN+CMB in a ΛCDM background, which are consistent with GR. Similar results are obtained when using the α-basis and assuming no-braiding (α B = 0) giving c M < 1.14 at 95% CL in a ΛCDM background, also in agreement with GR. However, we see a mild yet consistent indication for c B > 0 when α B is allowed to vary which will require further study to determine whether this is due to systematics or new physics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Modelling the Lyman-α forest with Eulerian and SPH hydrodynamical methods

ABSTRACT We compare two state-of-the-art numerical codes to study the overall accuracy in modelling the intergalactic medium and reproducing Lyman-α forest observables for DESI and high-resolution data sets. The codes employ different approaches to solving both gravity and modelling the gas hydrodynamics. The first code, Nyx, solves the Poisson equation using the Particle-Mesh (PM) method and the Euler equations using a finite-volume method. The second code, CRK-HACC , uses a Tree-PM method to solve for gravity, and an improved Lagrangian smoothed particle hydrodynamics (SPH) technique, where fluid elements are modelled with particles, to treat the intergalactic gas. We compare the convergence behaviour of the codes in flux statistics as well as the degree to which the codes agree in the converged limit. We find good agreement overall with differences being less than observational uncertainties, and a particularly notable ≲1 per cent agreement in the 1D flux power spectrum. This agreement was achieved by applying a tessellation methodology for reconstructing the density in CRK-HACC instead of using an SPH kernel as is standard practice. We show that use of the SPH kernel can lead to significant and unnecessary biases in flux statistics; this is especially prominent at high redshifts, z ∼ 5, as the Lyman-α forest mostly comes from lower-density regions that are intrinsically poorly sampled by SPH particles.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Stellar Characterization and Radius Inflation of Hyades M-dwarf Stars from the APOGEE Survey

We present a spectroscopic analysis of a sample of 48 M-dwarf stars (0.2 M ⊙ < M < 0.6 M ⊙ ) from the Hyades open cluster using high-resolution H-band spectra from the Sloan Digital Sky Survey/Apache Point Observatory Galactic Evolution Experiment (APOGEE) survey. Our methodology adopts spectrum synthesis with LTE MARCS model atmospheres, along with the APOGEE Data Release 17 line list, to determine effective temperatures, surface gravities, metallicities, and projected rotational velocities. The median metallicity obtained for the Hyades M dwarfs is [M/H] = 0.09 ± 0.03 dex, indicating a small internal uncertainty and good agreement with optical results for Hyades red giants. Overall, the median radii are larger than predicted by stellar models by 1.6% ± 2.3% and 2.4% ± 2.3%, relative to a MIST and DARTMOUTH isochrone, respectively. We emphasize, however, that these isochrones are different, and the fractional radius inflation for the fully and partially convective regimes have distinct behaviors depending on the isochrone. Using a MIST isochrone there is no evidence of radius inflation for the fully convective stars, while for the partially convective M dwarfs the radii are inflated by 2.7% ± 2.1%, which is in agreement with predictions from models that include magnetic fields. For the partially convective stars, rapid rotators present on average higher inflation levels than slow rotators. The comparison with SPOTS isochrone models indicates that the derived M-dwarf radii can be explained by accounting for stellar spots in the photosphere of the stars, with 76% of the studied M dwarfs having up to 20% spot coverage, and the most inflated stars with ~20%–40% spot coverage.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Monsoonal MCS Initiation, Rainfall, and Diurnal Gravity Waves over the Bay of Bengal: Observation and a Linear Model

Abstract Previous observational studies have indicated that mesoscale convective systems (MCSs) contribute the majority of precipitation over the Bay of Bengal (BoB) during the summer monsoon season, yet their initiation and propagation remain incompletely understood. To fill this knowledge gap, we conducted a comprehensive study using a combination of 20-yr satellite observations, MCS tracking, reanalysis data, and a theoretical linear model. Satellite observations reveal clear diurnal propagation signals of MCS initiation frequency and rainfall from the west coast of the BoB toward the central BoB, with the MCS rainfall propagating slightly slower than the MCS initiation frequency. Global reanalysis data indicate a strong association between the offshore-propagating MCS initiation frequency/rainfall and diurnal low-level wind perturbations, implying the potential role of gravity waves. To verify the hypothesis, we developed a 2D linear model that can be driven by realistic meteorological fields from reanalysis. The linear model realistically reproduces the characteristics of offshore-propagating diurnal wind perturbations. The wind perturbations, as well as the offshore propagation signals of MCS initiation frequency and rainfall, are associated with diurnal gravity waves emitted from the coastal regions, which in turn are caused by the diurnal land–sea thermal contrast. The ambient wind speed and vertical wind shear play crucial roles in modulating the timing, propagation, and amplitude of diurnal gravity waves. Using the linear model and satellite observations, we further show that the stronger monsoonal flows lead to faster offshore propagation of diurnal gravity waves, which subsequently control the offshore propagation signals of MCS initiation and rainfall. Significance Statement Rainfall over the Bay of Bengal (BoB) is primarily contributed by large and organized rainfall systems in the summer monsoon season. During this season, these systems are commonly observed over the east coast of India around midnight, the western BoB in the morning, and the central BoB in the afternoon. This eastward rainfall propagation is confirmed by observations, reanalysis data, and a theoretical model to have a strong association with atmospheric diurnal gravity waves. These waves are caused by land–sea thermal contrast and can trigger rainfall systems over the offshore regions. We also found that the diurnal gravity waves, as well as their triggered rainfall systems, can be greatly modulated by the large-scale monsoonal flows. All the above findings improve our understanding of diurnal rainfall cycle over the BoB and may contribute to the future improvement of rainfall forecast over the region.

Meteorology & Atmospheric Sciences↗
DOE OSTI2022

Detailed Chemical Abundances for a Benchmark Sample of M Dwarfs from the APOGEE Survey

Individual chemical abundances for 14 elements (C, O, Na, Mg, Al, Si, K, Ca, Ti, V, Cr, Mn, Fe, and Ni) are derived for a sample of M dwarfs using high-resolution, near-infrared H-band spectra from the Sloan Digital Sky Survey-IV/Apache Point Observatory Galactic Evolution Experiment (APOGEE) survey. The quantitative analysis included synthetic spectra computed with 1D LTE plane-parallel MARCS models using the APOGEE Data Release 17 line list to determine chemical abundances. The sample consists of 11 M dwarfs in binary systems with warmer FGK dwarf primaries and 10 measured interferometric angular diameters. To minimize atomic diffusion effects, [X/Fe] ratios are used to compare M dwarfs in binary systems and literature results for their warmer primary stars, indicating good agreement (<0.08 dex) for all studied elements. The mean abundance difference in primaries minus this work's M dwarfs is –0.05 ± 0.03 dex. It indicates that M dwarfs in binary systems are a reliable way to calibrate empirical relationships. A comparison with abundance, effective temperature, and surface gravity results from the APOGEE Stellar Parameter and Chemical Abundances Pipeline (ASPCAP) Data Release 16 finds a systematic offset of [M/H], T eff , log g = +0.21 dex, –50 K, and 0.30 dex, respectively, although ASPCAP [X/Fe] ratios are generally consistent with this study. The metallicities of the M dwarfs cover the range of [Fe/H] = –0.9 to +0.4 and are used to investigate Galactic chemical evolution via trends of [X/Fe] as a function of [Fe/H]. The behavior of the various elemental abundances [X/Fe] versus [Fe/H] agrees well with the corresponding trends derived from warmer FGK dwarfs, demonstrating that the APOGEE spectra can be used to examine Galactic chemical evolution using large samples of selected M dwarfs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Deciphering the Role of Total Water Storage Anomalies in Mediating Regional Flooding

Regional floods result from various flood generation mechanisms. Traditional analyses mainly link flooding to extreme rainfall, with limited input from soil moisture. Total water storage (TWS) is a holistic measure of basin wetness, including additional storage components from surface water, snow, and groundwater. Utilizing a new 5-day Gravity Recovery and Climate Experiment and its Follow On (GRACE(-FO)) data set, we investigated the linkage between short-term TWS anomaly (TWSA) and regional flooding. The 5-day TWSA solutions revealed flood signals missed by monthly TWSA solutions. Global basins exhibit distinct storage-discharge co-evolution patterns, offering new insights into flood mechanisms and propensity. Our bivariate event analyses show the annual maximum river discharges co-occur more often with the TWSA maxima than with precipitation in many basins. Further analyses revealed TWSA's time-lagged effect on river discharge, particularly in basins susceptible to floods triggered by saturation-excess runoff. The 5-day TWSA provides a new source of information for enhancing global flood preparedness.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Measuring cosmological distances using cluster edges as a standard ruler

ABSTRACT The line-of-sight velocity dispersion profile of galaxy clusters exhibits a ‘kink’ corresponding to the spatial extent of orbiting galaxies. Because the spatial extent of a cluster is correlated with the amplitude of the velocity dispersion profile, we can utilize this feature as a gravity-calibrated standard ruler. Specifically, the amplitude of the velocity dispersion data allows us to infer the physical cluster size. Consequently, observations of the angular scale of the ‘kink’ in the profile can be translated into a distance measurement to the cluster. Assuming the relation between cluster radius and cluster velocity dispersion can be calibrated from simulations, we forecast that with existing data from the Sloan Digital Sky Survey we will be able to measure the Hubble constant with 3.0 per cent precision. Implementing our method with data from the Dark Energy Spectroscopic Instrument (DESI) will result in a 1.3 per cent measurement of the Hubble constant. Adding cosmological supernova data improves the uncertainty of the DESI measurement to 0.7 per cent. While these error estimates are statistical only, they provide strong motivation for pursuing the necessary simulation program required to characterize and calibrate the systematic uncertainties impacting our proposed measurement. Whether or not our proposed measurement can in fact result in competitive H0 constraints will depend on what the eventual systematics floor for this method is.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

DEEPEN Leapfrog Geodata Model Cleaned and Reformatted Exploration Datasets from Newberry Volcano

DEEPEN stands for DE-risking Exploration of geothermal Plays in magmatic ENvironments. As part of the DEEPEN 3D play fairway analysis (PFA) conducted at Newberry Volcano for multiple play types (conventional hydrothermal, superhot EGS, and supercritical), existing geoscientific exploration datasets needed to be acquired, cleaned, reformatted, and assembled in Leapfrog Geothermal. This GDR submission includes all of the cleaned and reformatted (X (m), Y (m), elevation (m), processed data values) datasets used to build the Leapfrog Geodata model. Existing datasets were acquired from the GDR, from AltaRock, and from other sources. This yielded the following datasets: - Digital elevation model produced from LiDAR data by Ramsey and Bard, 2016 - MT surveys from 2006, 2011, 2014, and 2017 (including single inversions) - Gravity surveys from 2006, 2007, and 2011 (including single) - Earthquake catalogs from PNSN, LLNL, and the Newberry EGS Demonstration project - Seismic velocity model from Templeton et al., 2014 - The Frone, 2015 temperature model and a new one produced through extrapolating downhole temperature measurements and the SMU temperature at depth maps. Two versions of the new model are provided: 250 m spacing and 500 m spacing - EarthVision geologic model with alteration from Moser et al., 2016 - Well data from EGS well 55-29, deep geothermal wells, coreholes (GEO N-2 through 5) and several thermal gradient holes - "Newberry Well Data:" Location, simple lithology, directional survey data, and temperature data for the 34 wells and coreholes used in the Newberry PFA Although there are additional 2D datasets available in the area, such as aeromagnetic surveys, these were not included in the analysis. While it may be possible to project these datasets into three dimensions by assuming the surface measurements do not vary with depth, this method is associated with high uncertainty. Preexisting inversions of these data were unavailable, and inverting additional geophysical datasets is outside the scope of this project.

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

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↗