Search NASA⌕ Search

SEARCH · Search NASA

Results for “Spatial variation”

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 613 records · Page 34

MO-FG-CAMPUS-TeP3-02: Benchmarks of a Proton Relative Biological Effectiveness (RBE) Model for DNA Double Strand Break (DSB) Induction in the FLUKA, MCNP, TOPAS, and RayStation™ Treatment Planning System

Purpose: Validate implementation of a published RBE model for DSB induction (RBEDSB) in several general purpose Monte Carlo (MC) code systems and the RayStation™ treatment planning system (TPS). For protons and other light ions, DSB induction is a critical initiating molecular event that correlates well with the RBE for cell survival. Methods: An efficient algorithm to incorporate information on proton and light ion RBEDSB from the independently tested Monte Carlo Damage Simulation (MCDS) has now been integrated into MCNP (Stewart et al. PMB 60, 8249–8274, 2015), FLUKA, TOPAS and a research build of the RayStation™ TPS. To cross-validate the RBEDSB model implementation LET distributions, depth-dose and lateral (dose and RBEDSB) profiles for monodirectional monoenergetic (100 to 200 MeV) protons incident on a water phantom are compared. The effects of recoil and secondary ion production ({sub 2}H{sub +}, {sub 3}H{sub +}, {sub 3}He{sub 2+}, {sub 4}He{sub 2+}), spot size (3 and 10 mm), and transport physics on beam profiles and RBEDSB are examined. Results: Depth-dose and RBEDSB profiles among all of the MC models are in excellent agreement using a 1 mm distance criterion (width of a voxel). For a 100 MeV proton beam (10 mm spot), RBEDSB = 1.2 ± 0.03 (− 2–3%) at the tip of the Bragg peak and increases to 1.59 ± 0.3 two mm distal to the Bragg peak. RBEDSB tends to decrease as the kinetic energy of the incident proton increases. Conclusion: The model for proton RBEDSB has been accurately implemented into FLUKA, MCNP, TOPAS and the RayStation™TPS. The transport of secondary light ions (Z > 1) has a significant impact on RBEDSB, especially distal to the Bragg peak, although light ions have a small effect on (dosexRBEDSB) profiles. The ability to incorporate spatial variations in proton RBE within a TPS creates new opportunities to individualize treatment plans and increase the therapeutic ratio. Dr. Erik Traneus is employed full-time as a Research Scientist at RaySearch Laboratories. The research build of the RayStation used in the study was made available to the University of Washington free of charge. RaySearch Laboratories did not provide any monetary support for the reported studies.

60 APPLIED LIFE SCIENCES↗

MO-FG-CAMPUS-TeP3-05: Limitations of the Dose Weighted LET Concept for Intensity Modulated Proton Therapy in the Distal Falloff Region and Beyond

Purpose: Dose-weighted linear energy transfer (dLET) has been shown to be useful for the analysis of late effects in proton therapy. This study presents the results of the testing of the dLET concept for intensity modulated proton therapy (IMPT) with a discrete spot scanning beam system without use of an aperture or compensator (AC). Methods: IMPT (no AC) and broad beams (BB) with (AC) were simulated in the TOPAS and FLUKA code systems. Information from the independently tested Monte Carlo Damage Simulation (MCDS) was integrated into the FLUKA code systems to account for spatial variations in the RBE for protons and other light ions using an endpoint of DNA double strand break (DSB) induction. Results: The proton spectra for IMPT beams at the depths beyond the distal edge contain a tail of high energy protons up to 100 MeV. The integral from the tail is compatible with the number of 5–8 MeV protons at the tip of the Bragg peak (BP). The dose averaged energy (dEav) decreases to 7 MeV at the tip of (BP) and then increases to about 15 MeV beyond the distal edge. Neutrons produced in the nozzle are two orders of magnitude higher for BB with AC than for IMPT in low energy part of the spectra. The dLET values beyond of the distal edge of the BP are 5 times larger for the IMPT than for BB with the AC. Contrarily, negligible differences are seen in the RBE estimates for IMPT and BB with AC beyond the distal edge of the BP. Conclusion: The analysis of late effects in IMPT with a spot scanning and double scattering or scanning techniques with AC may requires both dLET and RBE as quantitative parameters to characterize effects beyond the distal edge of the BP.

60 APPLIED LIFE SCIENCES↗

Topology control of human fibroblast cells monolayer by liquid crystal elastomer

Eukaryotic cells in living tissues form dynamic patterns with spatially varying orientational order that affects important physiological processes such as apoptosis and cell migration. The challenge is how to impart a predesigned map of orientational order onto a growing tissue. Here, we demonstrate an approach to produce cell monolayers of human dermal fibroblasts with predesigned orientational patterns and topological defects using a photoaligned liquid crystal elastomer (LCE) that swells anisotropically in an aqueous medium. The patterns inscribed into the LCE are replicated by the tissue monolayer and cause a strong spatial variation of cells phenotype, their surface density, and number density fluctuations. Unbinding dynamics of defect pairs intrinsic to active matter is suppressed by anisotropic surface anchoring allowing the estimation of the elastic characteristics of the tissues. The demonstrated patterned LCE approach has potential to control the collective behavior of cells in living tissues, cell differentiation, and tissue morphogenesis.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Distinct Biogeographic Patterns for Archaea, Bacteria, and Fungi along the Vegetation Gradient at the Continental Scale in Eastern China

Understanding biogeographic patterns is a precursor to improving our knowledge of the function of microbiomes and to predicting ecosystem responses to environmental change. Using natural forest soil samples from 110 locations, this study is one of the largest attempts to comprehensively understand the different patterns of soil archaeal, bacterial, and fungal biogeography at the continental scale in eastern China. These patterns in natural forest sites could ascertain reliable soil microbial biogeographic patterns by eliminating anthropogenic influences. This information provides guidelines for monitoring the belowground ecosystem’s decline and restoration. Meanwhile, the deviations in the soil microbial communities from corresponding natural forest states indicate the extent of degradation of the soil ecosystem. Moreover, given the association between vegetation type and the microbial community, this information could be used to predict the long-term response of the underground ecosystem to the vegetation distribution caused by global climate change. The natural forest ecosystem in Eastern China, from tropical forest to boreal forest, has declined due to cropland development during the last 300 years, yet little is known about the historical biogeographic patterns and driving processes for the major domains of microorganisms along this continental-scale natural vegetation gradient. We predicted the biogeographic patterns of soil archaeal, bacterial, and fungal communities across 110 natural forest sites along a transect across four vegetation zones in Eastern China. The distance decay relationships demonstrated the distinct biogeographic patterns of archaeal, bacterial, and fungal communities. While historical processes mainly influenced bacterial community variations, spatially autocorrelated environmental variables mainly influenced the fungal community. Archaea did not display a distance decay pattern along the vegetation gradient. Bacterial community diversity and structure were correlated with the ratio of acid oxalate-soluble Fe to free Fe oxides (Feo/Fed ratio). Fungal community diversity and structure were influenced by dissolved organic carbon (DOC) and free aluminum (Ald), respectively. The role of these environmental variables was confirmed by the correlations between dominant operational taxonomic units (OTUs) and edaphic variables. However, most of the dominant OTUs were not correlated with the major driving variables for the entire communities. These results demonstrate that soil archaea, bacteria, and fungi have different biogeographic patterns and driving processes along this continental-scale natural vegetation gradient, implying different community assembly mechanisms and ecological functions for archaea, bacteria, and fungi in soil ecosystems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Single-colony MALDI mass spectrometry imaging reveals spatial differences in metabolite abundance between natural and cultured Trichodesmium morphotypes

Trichodesmium, a globally significant N 2 -fixing marine cyanobacterium, forms extensive surface blooms in nutrient-poor ocean regions. These blooms consist of a dynamic assemblage of Trichodesmium species that form distinct colony morphotypes and are inhabited by diverse microorganisms. Trichodesmium colony morphotypes vary in ecological niche, nutrient uptake, and organic molecule release, differentially impacting ocean carbon and nitrogen biogeochemical cycles. Here, we assessed the poorly studied spatial abundance of metabolites within and between three morphologically distinct Trichodesmium colonies collected from the Red Sea. We also compared these results with two morphotypes of the cultivable Trichodesmium strain IMS101. Using matrix-assisted laser desorption/ionization (MALDI) mass spectrometry imaging (MSI) coupled with liquid extraction surface analysis (LESA) tandem mass spectrometry (MS2), we identified and localized a wide range of small metabolites associated with single-colony Trichodesmium morphotypes. Our untargeted MALDI-MSI approach revealed 80 unique features (metabolites) shared between Trichodesmium morphotypes. Discrimination analysis showed spatial variations in 57 shared metabolites, accounting for 62% of the observed variation between morphotypes. The greatest variations in metabolite abundance were observed between the cultured morphotypes compared to the natural colony morphotypes, suggesting substantial differences in metabolite production between the cultivable strain IMS101 and the naturally occurring colony morphotypes that the cultivable strain is meant to represent. This study highlights the variations in metabolite abundance between natural and cultured Trichodesmium morphotypes and provides valuable insights into metabolites common to morphologically distinct Trichodesmium colonies, offering a foundation for future targeted metabolomic investigations.

59 BASIC BIOLOGICAL SCIENCES↗

Stimulated Brillouin scattering in the field of a two-dimensionally localized pumping wave

Stimulated Brillouin scattering of electromagnetic waves in the field of a two-dimensionally localized pump wave at arbitrary scattering angles in the regime of forward scattering is analyzed. Spatial variations in the amplitudes of interacting waves are studied for different values of the pump field and different dimensions of the pump wave localization region. The intensity of scattered radiation is determined as a function of the scattering angle and the dimensions of the pump wave localization region. It is shown that the intensity increases with increasing scattering angle.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Influence of grain charge gradients on the dynamics of macroparticles in an electrostatic trap

An analytical model of anomalous heating of charged dust grains (macroparticles) caused by their stochastic motion in a bounded plasma volume is proposed. Analytical expressions allowing one to describe the pumping (heating) of interacting grains with additional stochastic energy due to grain charge gradients are derived. The analytical results are verified by numerical simulation of the problem. It is shown that spatial variations in the charges of dust grains can lead to their anomalous heating in laboratory plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Do solar cycles influence giant cell arteritis and rheumatoid arthritis incidence?

Our objective was to examine the influence of solar cycle and geomagnetic effects on the incidence of giant cell arteritis (GCA) and rheumatoid arthritis (RA). Methods: We used data from patients with GCA (1950-2004) and RA (1955-2007) obtained from population-based cohorts. Yearly trends in age-adjusted and sex-adjusted incidence were correlated with the F10.7 index (solar radiation at 10.7 cm wavelength, a proxy for the solar extreme ultraviolet radiation) and AL index (a proxy for the westward auroral electrojet and a measure of geomagnetic activity). Fourier analysis was performed on AL, F10.7, and GCA and RA incidence rates. Results: The correlation of GCA incidence with AL is highly significant: GCA incidence peaks 0-1 year after the AL reaches its minimum (ie, auroral electrojet reaches a maximum). The correlation of RA incidence with AL is also highly significant. RA incidence rates are lowest 5-7 years after AL reaches maximum. AL, GCA and RA incidence power spectra are similar: they have a main peak (periodicity) at about 10 years and a minor peak at 4-5 years. However, the RA incidence power spectrum main peak is broader (8-11 years), which partly explains the lower correlation between RA onset and AL. The auroral electrojets may be linked to the decline of RA incidence more strongly than the onset of RA. The incidences of RA and GCA are aligned in geomagnetic latitude. Conclusions: AL and the incidences of GCA and RA all have a major periodicity of about 10 years and a secondary periodicity at 4-5 years. Geomagnetic activity may explain the temporal and spatial variations, including east-west skewness in geographic coordinates, in GCA and RA incidence, although the mechanism is unknown. Lastly, the link with solar, geospace and atmospheric parameters need to be investigated. These novel findings warrant examination in other populations and with other autoimmune diseases.

58 GEOSCIENCES↗

Polarization phenomenon in heavy-ion collisions

The strongly interacting system created in ultrarelativistic nuclear collisions behaves almost as an ideal fluid with rich patterns of the velocity field exhibiting strong vortical structure. Vorticity of the fluid, via spin-orbit coupling, leads to particle spin polarization. Due to the finite orbital momentum of the system, the polarization on average is not zero. It depends on the particle momenta reflecting the spatial variation of the local vorticity. In the last few years, this field experienced a rapid growth due to experimental discoveries of the global and local polarizations. Recent measurements triggered further development of the theoretical description of the spin dynamics and suggestions of several new mechanisms for particle polarization. In this review, we focus mostly on the experimental results. We compare the measurements with the existing theoretical calculations but try to keep the discussion of possible underlying physics at the qualitative level. Future measurements and how they can help to answer open theoretical questions are also discussed. We pay a special attention to the employed experimental methods, as well as to the detector effects and associated corrections to the measurements.

Physics↗

Design Principles to Govern Electrode Fabrication for the Lithium Trivanadate Cathode

A full depth of discharge mathematical model for the lithium trivanadate cathode, considering lithiation of the layered α-phase, phase change, and lithiation of the rock-salt like β-phase at lower potentials, is developed. The coupled electrode-scale and crystal-scale model is fit to electrochemical data, and additionally validated with operando EDXRD. There is good agreement between the simulated and measured spatial variation of the volume fraction of the β-phase. This mathematical model is used to guide electrode fabrication, accounting for both ionic and electronic transport effects. Values of three design parameters—electrode thickness, porosity, and volume fraction of conductor—are identified, and the sensitivity of the energy density to these design parameters is quantified. The model is also used to investigate electrode design to create electrodes that deliver the maximum achievable energy density under the constraint that the α to β-phase transition is avoided, since phase change has been demonstrated to reduce cycle life. The energy density sacrificed to avoid phase change decreases at higher discharge rates, but the target values for electrode fabrication remain the same as those when optimizing the electrode for the full depth of discharge.

25 ENERGY STORAGE↗

Plating and Stripping of Lithium Metal Stabilized by a Block Copolymer Electrolyte: Local Current Density Measurement and Modeling

Here, plating and stripping of lithium protrusions in lithium metal symmetric cells containing a solid block copolymer electrolyte was studied as a function of time in 3D using time-resolved X-ray tomography. These measurements enabled determination of the spatial variation in current densities at the plating and stripping electrodes. The initial interelectrode distance was 27 μm. Correlation functions were calculated to reveal the relationships between current densities at the two electrodes and local electrolyte thickness. Current densities at opposing electrode locations during protrusion growth is uncorrelated until the local interelectrode distance decreases to less than 6 μm, just before the cell shorts. Mass balance was used to determine the area from which lithium ions that form a protrusion were stripped. Computational modeling of the plating and stripping process reveals the interplay between electrochemical and mechanical driving forces and their effect on nonuniform current distribution. Model predictions were compared with experiments without resorting to any adjustable parameters. The computed correlation functions were in qualitative agreement with experiments. Finally, the model was used to calculate contour plots of electrochemical potential within the electrolyte, shedding light on how geometry, salt concentration, interelectrode distance, and mechanical stress influence local rates of electrochemical reaction.

36 MATERIALS SCIENCE↗

Exploring the Limits of Bottom-Up Gold Filling to Fabricate Diffraction Gratings

Gold deposition on rotating disk electrodes, Bi 3+ adsorption on planar Au films and superconformal Au filling of trenches up to 45 μm deep are examined in Bi 3+ -containing Na 3 Au(SO 3 ) 2 electrolytes with pH between 9.5 and 11.5. Higher pH is found to increase the potential-dependent rate of Bi 3+ adsorption on planar Au surfaces, shortening the incubation period that precedes active Au deposition on planar surfaces and bottom-up filling in patterned features. Decreased contact angles between the Au seeded sidewalls and bottom-up growth front also suggest improved wetting. The bottom-up filling dynamic in trenches is, however, lost at pH 11.5. The impact of Au concentration, 80 mmol/L versus 160 mmol/L Na 3 Au(SO 3 ) 2 , on bottom-up filling is examined in trenches up to ≈ 210 μm deep with aspect ratio of depth/width ≈ 30. The microstructures of void-free, bottom-up filled trench arrays used as X-ray diffraction gratings are characterized by scanning electron microscopy (SEM) and Electron Backscatter Diffraction (EBSD), revealing marked spatial variation of the grain size and orientation within the filled features.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Geometry-Driven Longitudal Topic Model

A simple and scalable framework for longitudinal analysis of Twitter data is developed that combines latent topic models with computational geometric methods. Dimensionality reduction tools from computational geometry are applied to learn the intrinsic manifold on which the latent, temporal topics reside. Then shortest path distances on the manifold are used to link together these topics. The proposed framework permits visualization of the low-dimensional embedding which provides clear interpretation of the complex, high-dimensional trajectories that may exist among latent topics. Practical application of the proposed framework is demonstrated through its ability to capture and effectively visualize natural progression of latent COVID-19 related topics learned from Twitter data. Interpretability of the trajectories is achieved by comparing to real-world events. In addition, the framework permits study of spatial variation in Twitter behavior for learned topics. The analysis demonstrates that the proposed framework is able to capture granular-level impact of COVID-19 on public discussions. We end by arguing that Twitter data, when analyzed within the proposed framework, can serve as a valuable supplementary data stream for COVID-related studies.

97 MATHEMATICS AND COMPUTING↗

Implementing and Evaluating Variable Soil Thickness in the Community Land Model, Version 4.5 (CLM4.5)

One of the recognized weaknesses of land surface models as used in weather and climate models is the assumption of constant soil thickness because of the lack of global estimates of bedrock depth. As such, using a 30-arc-s global dataset for the thickness of relatively porous, unconsolidated sediments over bedrock, spatial variation in soil thickness is included here in version 4.5 of the Community Land Model (CLM4.5). The number of soil layers for each grid cell is determined from the average soil depth for each 0.9° latitude × 1.25° longitude grid cell. The greatest changes in the simulation with variable soil thickness are to baseflow, with the annual minimum generally occurring earlier. Smaller changes are seen in latent heat flux and surface runoff primarily as a result of an increase in the annual cycle amplitude. These changes are related to soil moisture changes that are most substantial in locations with shallow bedrock. Total water storage (TWS) anomalies are not strongly affected over most river basins since most basins contain mostly deep soils, but TWS anomalies are substantially different for a river basin with more mountainous terrain. Additionally, the annual cycle in soil temperature is partially affected by including realistic soil thicknesses resulting from changes in the vertical profile of heat capacity and thermal conductivity. However, the largest changes to soil temperature are introduced by the soil moisture changes in the variable soil thickness simulation. This implementation of variable soil thickness represents a step forward in land surface model development.

54 ENVIRONMENTAL SCIENCES↗

Uncertainties in Projecting Future Changes in Atmospheric Rivers and Their Impacts on Heavy Precipitation over Europe

This study investigates the North Atlantic atmospheric rivers (ARs) making landfall over western Europe in the present and future climate from the multi-model ensemble of the Coupled Model Intercomparison Project Phase 5 (CMIP5). Overall, CMIP5 captures the seasonal and spatial variations of historical landfalling AR days, with the large inter-model variability strongly correlated with the inter-model spread of historical jet position. Under RCP 8.5, AR frequency is projected to increase a few times by the end of this century. While thermodynamics plays a dominate role in the future increase of ARs, wind changes associated with the midlatitude jet shifts also significantly contribute to AR changes, resulting in dipole change patterns in all seasons. In the North Atlantic, the model projected jet shifts are strongly correlated with the simulated historical jet position. As models exhibit predominantly equatorward biases in the historical jet position, the large poleward jet shifts reduce AR days south of the historical mean jet position through the dynamical connections between the jet positions and AR days. Using the observed historical jet position as an emergent constraint, dynamical effects further increase AR days in the future above the large increases due to thermodynamical effects. In the future, both total and extreme precipitation induced by AR contribute more to the seasonal mean and extreme precipitation compared to present primarily because of the increase in AR frequency. While AR precipitation intensity generally increases more relative to the increase in integrated vapor transport, AR extreme precipitation intensity increases much less.

atmospheric rivers↗

Modeling and Calibration of Supplier Selection Problem in Freight Agent-Based Simulations

Freight transportation modeling often struggles with data limitations, especially in accurately representing complex supplier selection processes and their impact on network flows. This research addresses this critical gap by developing a large-scale, calibrated agent-based model for supplier selection, complemented by a probabilistic heuristic for international shipments. Our approach integrates trade relationships between industry sectors, transportation costs, and a supplier-rating model adapted from existing literature. The model’s core objective is to minimize the discrepancy between modeled and observed commodity flows while ensuring a close match to regional shipping distance distributions. Implemented and tested across four major U.S. metropolitan areas—Atlanta, Chicago, Dallas–Fort Worth, and Los Angeles—the model demonstrates high fidelity in replicating observed freight patterns. Key findings reveal consistent alignment with national shipping distance trends and highlight significant spatial variations in commodity trade assignments and demand across the study regions. This behaviorally informed and transport-sensitive framework is designed to approximate real-world decision making, providing a robust tool for policymakers and planners to evaluate targeted interventions, assess infrastructure investments, and enhance supply chain resilience in the face of disruptions.

Ismael, Abdelrahman (ORCID:0000000303712110)↗

Verification of diesel spray ignition phenomenon in dual-fuel diesel-piloted premixed natural gas engine

Dual-fuel (DF) engines, in which premixed natural gas and air in an open-type combustion chamber is ignited by diesel-fuel pilot sprays, have been more popular for marine use than pre-chamber spark ignition (PCSI) engines because of their superior durability. However, control of ignition and combustion in DF engines is more difficult than in PCSI engines. In this context, here we focus on the ignition stability of n-heptane pilot-fuel jets injected into a compressed premixed charge of natural gas and air at low-load conditions. To aid understanding of the experimental data, chemical-kinetics simulations were carried out in a simplified engine-environment that provided insight into the chemical effects of methane (CH4) on pilot-fuel ignition. The simulations reveal that CH4 has an effect on both stages of n-heptane autoignition: the small, first-stage, cool-flame-type, low-temperature ignition (LTI) and the larger, second-stage, high-temperature ignition (HTI). As the ratio of pilot-fuel to CH4 entrained into the spray decreases, the initial oxidization of CH4 consumes the OH radicals produced by pilot-fuel decomposition during LTI, thereby inhibiting its progression to HTI. Using imaging diagnostics, the spatial and temporal progression of LTI and HTI in DF combustion are measured in a heavy-duty optical engine, and the imaging data are analyzed to understand the cause of severe fluctuations in ignition timing and combustion completeness at low-load conditions. Images of cool-flame and hydroxyl radical (OH*) chemiluminescence serve as indicators of LTI and HTI, respectively. The cycle-to-cycle and spatial variation in ignition extracted from the imaging data are used as key metrics of comparison. The imaging data indicate that the local concentration of the pilot-fuel and the richness of the surrounding natural-gas air mixture are important for LTI and HTI, but in different ways. In particular, higher injection pressures and shorter injection durations increase the mixing rate, leading to lower concentrations of pilot-fuel more quickly, which can inhibit HTI even as LTI remains relatively robust. Decreasing the injection pressure from 80 MPa to 40 MPa and increasing the injection duration from 500 µs to 760 µs maintained constant pilot-fuel mass, while promoting robust transition from LTI to HTI by effectively slowing the mixing rate. This allows enough residence time for the OH radicals, produced by the two-stage ignition chemistry of the pilot-fuel, to accelerate the transition from LTI to HTI before being consumed by CH4 oxidation. Thus from a practical perspective, for a premixed natural gas fuel–air equivalence-ratio, it is possible to improve the “stability” of the combustion process by solely manipulating the pilot-fuel injection parameters while maintaining constant mass of injected pilot-fuel. This allows for tailoring mixing trajectories to offset changes in fuel ignition chemistry, so as to promote a robust transition from LTI to HTI by changing the balance between the local concentration of the pilot-fuel and richness of the premixed natural gas and air. This could prove to be a valuable tool for combustion design to improve fuel efficiency or reduce noise or perhaps even reduce heat-transfer losses by locating early combustion away from in-cylinder walls.

42 ENGINEERING↗

Influence of extreme fracture flow channels on the thermal performance of open-loop geothermal systems at commercial scale

Abstract Adequate stewardship of geothermal resources requires accurate forecasting of long-term thermal performance. In enhanced geothermal systems and other fracture-dominated reservoirs, predictive models commonly assume constant-aperture fractures, although spatial variations in aperture can greatly affect reservoir permeability, fluid flow distribution, and heat transport. Whereas previous authors have investigated the effects of theoretical random aperture distributions on thermal performance, here we further explore the influence of permeability heterogeneity considering field-constrained aperture distributions from a meso-scale field site in northern New York, USA. Using numerical models of coupled fluid flow and heat transport, we conduct thermal–hydraulic simulations for a hypothetical reservoir consisting of a relatively impervious porous matrix and a single, horizontal fracture. Our results indicate that in highly channelized fields, most well design configurations and operating conditions result in extreme rates of thermal drawdown (e.g., 50% drop in production well temperatures in under 2 years). However, some other scenarios that account for the risks of short-circuiting can potentially enhance heat extraction when mass flow rate is not excessively high, and the direction of geothermal extraction is not aligned with the most permeable features in the reservoir. Through a parametric approach, we illustrate that well separation distance and relative positioning play a major role in the long-term performance of highly channelized fields, and both can be used to help mitigate premature thermal breakthrough.

Rangel-Jurado, Nicolás↗