Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “potential inversion”

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 235 records · Page 13

Modulating the dynamics of Brønsted acid sites on PtWOx inverse catalyst

Here metal-metal oxide (M-MO) inverse catalysts are broadly applied. Brønsted acid sites on the oxide overlayers are often hypothesized to drive selective C-O bond activation. However, the Brønsted acid site nature and dynamics under working conditions remain poorly understood due to multiple materials functionalities. Here, we investigate the formation and the dynamics of Brønsted acid and redox sites on PtWO x /C under working conditions. DFT-based thermodynamic calculations and microkinetic modeling reveal a complex interplay between Brønsted acid and redox sites and potentially fast catalyst dynamics at comparable time scales to the chemistry. Combining in situ characterization and probe chemistry, we demonstrate that the density of Brønsted acid sites on the PtWO x /C inverse catalyst could be modulated by up to two orders of magnitude by altering the reaction parameters and by the chemistry itself. We elicit an order of magnitude increase in the acid-catalyzed dehydration average reaction rate by periodic hydrogen pulsing.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Inversion of Many-Beam Bragg Intensities for Phasing by Iterated Projections: Removal of Multiple Scattering Artifacts from Diffraction Data

An iterated projection algorithm (N-Phaser) is developed that reconstructs a scattering potential from N-beam multiple Bragg scattered intensities. The method may be used to eliminate multiple scattering artifacts from electron diffraction data, solving the phase problem and increasing the thicknesses of samples used in materials science, solid-state chemistry, and small molecule crystallography. For high-energy transmission electron diffraction, we show that the algorithm recovers accurate complex structure factors from a wide range of thicknesses, orientations, and relativistic beam energies, and does not require known thickness or atomic-resolution data if sufficient multiple scattering occurs. Furthermore, extensions to Cryo-electron microscopy and Micro-electron diffraction are suggested.

36 MATERIALS SCIENCE↗

Inversion of Time-Lapse Seismic Reservoir Monitoring Data Using CycleGAN: A Deep Learning-Based Approach for Estimating Dynamic Reservoir Property Changes

Carbon capture and storage is being pursued globally as a geoengineering measure for reducing the emission of anthropogenic CO 2 the atmosphere. Comprehensive monitoring, verification, and accounting programs must be established for demonstrating the safe storage of injected CO 2 . One of the most commonly deployed monitoring techniques is time-lapse seismic reservoir monitoring (also known as 4-D seismic), which involves comparing 3-D seismic survey data taken at the same study site but over different times. Analyses of 4-D seismic data volumes can help improve the quality of storage reservoir characterization, track the movement of injected CO 2 plume, and identify potential CO 2 spillover/leakage from the storage reservoirblue. However, the derivation of high-resolution CO 2 saturation maps from 4-D seismic data is a highly nonlinear and ill-posed inverse problem, often requiring significant computational effort. In this research, we apply a physics-based deep learning method to facilitate the solution of both the forward and inverse problems in seismic inversion while honoring physical constraints. A cycle generative adversarial neural network (CycleGAN) model is trained to learn the bidirectional functional mappings between the reservoir dynamic property changes and seismic attribute changes, such that both forward and inverse solutions can be obtained efficiently from the trained model. We show that our CycleGAN-based approach not only improves the reliability of 4-D seismic inversion but also expedites the quantitative interpretation. Our deep learning-based workflow is generic and can be readily used for reservoir characterization and reservoir model updates involving the use of 4-D seismic data.

58 GEOSCIENCES↗

SQMS Quantum R&D in Machine Learning, Optimization and Sensing beyond Fundamental Physics Applications

This newly formed team at SQMS under the Ecosystem Thrust is looking to develop capabilities impacting societal advances outside the core domain of HEP and condensed matter physics. We explicitly leverage the experimental and algorithmic innovations developed across all groups as well as connect to broad-scope external projects of the diverse team of PIs. As the inaugural set of projects, we are studying numerically quantum machine learning models inspired by efficiently trainable echo-state and orthogonal neural networks and developing designs for related experiments to be performed on quantum processors based on SQMS SRF cQED technology and Rigetti s transmon arrays. Investigated models exploit ideas and lessons learned from multiple prior work by SQMS team members in a variety of internal and external activities [R1]. Target initial applications include noisy signal processing, potentially captured by quantum sensors or noisy QPUs, as well as simulation and classification of healthcare data. For instance, image reconstruction of the brain s electrical properties by solving the inverse Maxwell equation problem with uncertainty [R2] through a hybrid quantum-classical physics-informed architecture for time-dependent processes [R3]. The group is also investigating the application and development of novel quantum sensors based on magnetic levitation of a superconducting sphere coupled to a superconducting qubit. This coupling enables high-precision measurements of the position of the sphere, which can be used for sensitive detection of forces, enabling practical applications such as gravimetry for geophysics analysis, or accelerometry for GPS-denied navigation [R4] [R1] Rieffel, Eleanor G., Ata Akbari Asanjan, M. Sohaib Alam, Namit Anand, David E. Bernal Neira, Sophie Block, Lucas T. Brady et al. "Assessing and advancing the potential of quantum computing: A NASA case study." Future Generation Computer Systems (2024). [R2] Yu, X., Serrall s, J.E., Giannakopoulos, I.I., Liu, Z., Daniel, L., Lattanzi, R. and Zhang, Z., 2023. Pifon-ept: Mr-based electrical property tomography using physics-informed fourier networks. IEEE Journal on Multiscale and Multiphysics Computational Techniques. [R3] Wudarski, Filip, Daniel OConnor, Shaun Geaney, Ata Akbari Asanjan, Max Wilson, Elena Strbac, P. Aaron Lott, and Davide Venturelli. "Hybrid quantum-classical reservoir computing for simulating chaotic systems." arXiv preprint arXiv:2311.14105 (2023). [R4] Higgins, Gerard, Saarik Kalia, and Zhen Liu. "Maglev for dark matter: Dark-photon and axion dark matter sensing with levitated superconductors." Physical Review D 109.5 (2024): 055024.

Venturelli, Davide↗

Open-source library for performance-portable neutrino reaction rates: Application to neutron star mergers

A realistic and detailed description of neutrinos in binary neutron star (BNS) mergers is essential to build reliable models of such systems. To this end, we present bns_nurates, a novel open-source numerical library designed for the efficient on-the-fly computation of neutrino interactions, with particular focus on regimes relevant to BNS mergers. bns_nurates targets a higher level of accuracy and realism in the implementation of commonly employed reactions by accounting for relevant microphysics effects on the interactions, such as weak magnetism and mean field effects. It also includes the contributions of inelastic neutrino scattering off electrons and positrons and (inverse) nucleon decays. Finally, it offers a way to reconstruct the neutrino distribution function in the framework of moment-based transport schemes. As a first application, we compute both energy-dependent and energy-integrated neutrino emissivities and opacities for conditions extracted from a BNS merger simulation with m1 transport scheme. We find some qualitative differences in the results when considering the impact of the additional relevant reactions and of microphysics effects. For example, neutrino-electron/positron scattering reactions are important for the energy exchange of heavy-type neutrinos as they do not undergo semileptonic charged-current processes, when μ± are not accounted for. Moreover, weak magnetism and mean field effects can significantly modify the contribution of β processes for electron-type (anti)neutrinos, increasing at the same time the importance of (inverse) neutron decays. Here, the improved treatment for the reaction rates also modifies the conditions at which neutrinos decouple from matter in the system, potentially affecting their emission spectra.

79 ASTRONOMY AND ASTROPHYSICS↗

A Scalable Gaussian Process Approach to Shear Mapping with MuyGPs

Analysis of cosmic shear is an integral part of understanding structure growth across cosmic time, which in turn provides us with information about the nature of dark energy. Conventional methods generate shear maps from which we can infer the matter distribution in the universe. Current methods (e.g., Kaiser–Squires inversion) for generating these maps, however, are tricky to implement and can introduce bias. Recent alternatives construct a spatial process prior for the lensing potential, which allows for inference of the convergence and shear parameters given lensing shear measurements. Realizing these spatial processes, however, scales cubically in the number of observations—an unacceptable expense as near-term surveys expect billions of correlated measurements. Therefore, we present a linearly scaling shear map construction alternative using a scalable Gaussian process prior called MuyGPs. MuyGPs avoids cubic scaling by conditioning interpolation on only nearest neighbors and fits hyperparameters using batched leave-one-out cross-validation. This work is the first step toward a full, scalable mass mapping method. We work in a simplified regime where we validate our method by interpolating and analyzing maps given noisy point-estimate data from all three shear fields, taken from a suite of N -body ray-tracing simulations. We also show that we can perform these operations at the scale of billions of galaxies on high-performance computing platforms.

79 ASTRONOMY AND ASTROPHYSICS↗

Efficient Simulation of Open Quantum Systems on NISQ Trapped‐Ion Hardware

Abstract Simulating open quantum systems, which interact with external environments, presents significant challenges on noisy intermediate‐scale quantum (NISQ) devices due to limited qubit resources and noise. In this study, an efficient framework is proposed for simulating open quantum systems on NISQ hardware by leveraging a time‐perturbative Kraus operator representation of the system's dynamics. This approach avoids the computationally expensive Trotterization method and exploits the Lindblad master equation to represent time evolution in a compact form, particularly for systems satisfying specific commutation relations. The efficiency of this method is demonstrated by simulating quantum channels, such as the continuous‐time Pauli channel and damped harmonic oscillators, on NISQ trapped‐ion hardware, including IonQ Harmony and Quantinuum H1‐1. Additionally, hardware‐agnostic error mitigation techniques are introduced, including Pauli channel fitting and quantum depolarizing channel inversion, to enhance the fidelity of quantum simulations. These results show strong agreement between the simulations on real quantum hardware and exact solutions, highlighting the potential of Kraus‐based methods for scalable and accurate simulation of open quantum systems on NISQ devices. This framework opens pathways for simulating more complex systems under realistic conditions in the near term.

Burdine, Colin [Department of Electrical and Compu↗

Input specific neural networks

Neural networks have emerged as powerful tools for mapping between inputs and outputs. However, their black-box nature limits the ability to encode or impose specific structural relationships between inputs and outputs. Many scientific and engineering problems, such as constitutive modeling in solid mechanics, require networks that can enforce convexity, monotonicity, or other structural constraints to ensure physical consistency. Here, we introduce the Input Specific Neural Network (ISNN), a new architecture that enables multiple, distinct constraints to be imposed on different input subsets for scalar-valued outputs. This framework unifies convex, monotone–convex, monotone, and arbitrary mappings within a single network for the first time. Two ISNN architectures with analytical first- and second-order derivatives are developed. We demonstrate the performance on synthetic toy problems, inverse problems in isotropic hyperelasticity, and finite element simulations. ISNNs achieve improved extrapolation behavior, require fewer invariant inputs than standard input convex networks for polyconvex potentials, and enable significant computational savings via manual differentiation. We also show how ISNNs can be used to learn structural relationships between inputs and outputs via a binary gating mechanism. Particularly, ISNNs are employed to model a homogenized anisotropic free energy potential in a decoupled multiscale setting. The network learns whether or not the potential should be modeled as polyconvex and retains only the relevant layers while using the minimum number of inputs. ISNNs provide a flexible foundation for embedding structural priors into neural networks, enhancing both interpretability and stability. They are broadly applicable across computational mechanics and other scientific domains requiring constrained functional relationships.

Jadoon, Asghar A. [Univ. of Texas, Austin, TX (Uni↗

Reactor neutrino applications and coherent elastic neutrino nucleus scattering

Potential applications of neutrino detection to nuclear security have been discussed since the 1970s. Recent years have seen great progress in detector technologies based on inverse beta decay, with the demonstration of ton-scale surface-level detectors capable of high quality neutrino spectrum measurements. At the same time coherent elastic neutrino nucleus scattering has been experimentally confirmed in 2017 with neutrinos from stopped pion decay and there is a number of experiments aimed at seeing this reaction with reactor neutrinos. Here, the large cross section and threshold-less nature of this reaction make it plausible to consider it for applications to nuclear security and here, we present a first direct comparison of the two reaction modes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Self-aligned heterogeneous quantum photonic integration

Integrated quantum photonics holds significant promise for scalable photonic quantum information processing, quantum repeaters, and quantum networks, but its development is hindered by the mismatch between materials hosting high-quality quantum emitters and those compatible with mature photonic technologies. Heterogeneous integration offers a potential solution to this challenge, yet practical implementations have been limited by inevitable insertion losses at material interfaces. Here, we present a self-aligned heterogeneous quantum photonic integration approach that enables near-unity coupling efficiency at the interface. To showcase our approach, we demonstrate Purcell enhancement of a silicon vacancy (SiV) center in diamond induced by a heterogeneous photonic crystal cavity defined by titanium dioxide (TiO 2 ), as well as optical spin control and readout via a TiO 2 photonic circuit. We further show that, when combined with inverse photonic design, our approach enables efficient and broadband collection of single photons from a color center into a heterogeneous waveguide. Our approach is not restricted to SiV centers or TiO 2 ; it has the potential to be broadly applied to integrate diverse solid-state quantum emitters with thin-film photonic devices where conformal deposition is possible. Together, these results establish a practical route to scalable quantum photonic integrated circuits that combine high-quality quantum emitters with technologically mature photonic platforms.

36 MATERIALS SCIENCE↗

Boundary Layer Height Variation over the Savannah River Site during a Heat Wave

With the number of heat waves expect to increase in both in their severity and frequency in the coming years, an insight into their impact on planetary boundary layer (PBL) heights is essential for one development of future emergency and safety plans. Using High-Resolution Rapid Refresh (HRRR) atmospheric model output, an analysis is done on the effects a heat wave has on PBL heights and temperature inversions across the Savannah River Site (SRS). This study showed that due to its adjacency to anthropogenic land use, SRS experiences a non-negligible gradient in its PBL height, with higher heights to the northwest. The singular location PBL heigh measurement and forecasts for SRS do not account for this gradient, which has potential implications on pollutant and smoke dispersion onsite. During a heat wave this gradient is present but PBL heights are less impacted than surrounding developed areas, possibly attributable to the woody wetlands that surround the Savannah River's banks. Additionally, the heat wave appears to strongly enhance the nocturnal temperature inversions that the site and surrounding areas experience. This is possibly due to increased stability from higher pressure conditions, which are often associated with heat waves. This influence on nocturnal inversions, however, is limited to their intensity, leaving their timing relatively unaltered. The impact a heat wave has on PBL heights and temperature inversions across SRS is seen to be more dependent on adjacent land use and features rather than the unique land use types onsite.

54 ENVIRONMENTAL SCIENCES↗

Enhancement of the spin–orbit torque efficiency in W/Cu/CoFeB heterostructures via interface engineering

Here, the spin-torque ferromagnetic resonance signal and the spin Hall magnetoresistance induced by the spin Hall effect of W/Cu/CoFeB heterostructures with different Cu layer thicknesses (t Cu ) have been systemically studied. The effective spin mixing conductance $g^{↑↓}_{eff}$, the interfacial spin transparency T, and the real spin–orbit torque efficiency ( J s J c ) real show a significant increase compared to the W/CoFeB heterostructure. ( J s J c ) real reaches its maximum of ~0.54, increased up to ~50% at the optimized tCu ~ 0.52 nm according to our theoretical prediction. More importantly, the intrinsic spin Hall angle of W, $θ^{int}_{SH}$ ~ 0.79 ± 0.20, has also been obtained after the correction of the inverse spin Hall effect and T. This suggests that the Cu insertion improves the interface quality and, therefore, assists the spin transport in the heterostructures, which potentially improves the performance of next-generation spintronic devices.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Enhanced Dirac node separation in the strained Cd 3 As 2 topological semimetal

In topological semimetals, nodes appear at symmetry points in the Brillouin zone as a result of band inversion, and yield quasirelativistic massless fermions at low energies. Cd 3 As 2 is a three-dimensional topological semimetal that hosts two Dirac cones responsible for a variety of quantum phenomena. In this work, we demonstrate the strain tuning of the Dirac nodes of Cd 3 As 2 through a combination of magneto-optical infrared spectroscopy and high-resolution x-ray-diffraction studies performed on epitaxial films. In these thin films, we observe a giant enhancement of the node separation in momentum space by close to a factor of 4. A combination of experimental measurements and theoretical modeling allows us to relate the origin of this enhancement to a strengthening of the topological band inversion driven by lattice strain. Finally, our results demonstrate how strain can be used as a knob to tune the topological properties of semimetals and to potentially enhance their performance and response for various applications.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

STILGAR End-of-Project Report

The Subsurface Tunnel Imaging LeveraGed by Analysis of Rayleigh wave ellipticity (STILGAR) project demonstrated an integrated geophysical approach for detecting, locating, and characterizing underground structural changes using dense seismic arrays and advanced inversion techniques. Field campaigns were conducted at two operational mines—the Redmond salt mine (Utah) and Graymont Pleasant Gap limestone mine (Pennsylvania)—providing real-world testbeds for monitoring anthropogenic subsurface activity. At the Redmond salt mine, seismic interferometry combined with back-projection inversion successfully identified continuous, low-amplitude signals from mining operations. The approach differentiated stationary from migrating anthropogenic sources, captured daily operational cycles, and validated the potential of passive seismic monitoring for remote detection of underground activity. At the Graymont Pleasant Gap mine, two dense seismic deployments in the spring and fall of 2023 generated over 4 TB of high-resolution data. Key outcomes included the relocation of 199 underground and 8 surface explosions with accuracies within tens of meters and the development of a 3D P-wave velocity model using the triple-difference tomography algorithm (tomoTD) that resolved major structural features such as the mine entrance, low-velocity tunnels, and roof-collapse areas. Ambient noise cross-correlation and back-projection analyses revealed persistent sources linked to ongoing mining activity, whereas horizontal-to-vertical spectral ratio (HVSR) and ellipticity studies confirmed stable site responses across seasons and identified soil thickness trends consistent with regional erosional and depositional processes. Checkerboard and sensitivity tests further validated the robustness of the tomographic results. Overall, the findings emphasize that although significant progress has been made in subsurface imaging, further work is needed to enhance the detection and localization of underground structures. Accurate imaging requires higher frequencies, yet anthropogenic sources tend to dominate the seismic record at those frequencies, and high-frequency surface waves are affected by higher modes that complicate interpretation. The improved detection and localization of human-induced signals enabled detailed temporal and spatial mapping of daily mine operations, demonstrating the feasibility of continuous anthropogenic source monitoring. Sensitivity to signals from nontraditional sources, such as fan operations, highlights the broader applicability of this approach to other industrial environments in which continuous and impulsive signals are present. The field campaigns produced a substantial volume of high-quality seismic data, supporting the development and testing of new methods for seismic source characterization and subsurface imaging. Future deployments should include sensors capable of recording lower frequencies to probe deeper structures, increase bandwidth to enhance resolution and sensitivity to both shallow and deep targets, and collect additional large-scale datasets to refine imaging and source characterization techniques. Moreover, conducting 3D modeling studies of seismic wavefields at higher frequencies will provide a better understanding of wave scattering and cavity–wavefield interactions in complex underground environments. In conclusion, the STILGAR project demonstrated that integrated seismic monitoring can effectively characterize underground operations, capturing both natural and anthropogenic signals. The approaches developed provide a foundation for improved detection, localization, and imaging of subsurface structures and are directly transferable to broader industrial monitoring applications.

58 GEOSCIENCES↗

Surface Charge Distribution Mapping Using an Electrostatic Probe

Undesired accumulation of charge on dielectric materials causing electrostatic discharges can be an issue in pulsed power systems and electrical systems in general. As such, an understanding of surface charge distribution, charge accumulation, and decay is required. An in-house no-contact electrostatic probe designed with a 2 mm resolution was fabricated to measure and map surface charge distribution. Designed as a contactless instrument, the probe consisting of a metal rod and low leakage amplifier circuit probes the electric potential near the surface. A scan of the surface then provides the raw data, a potential distribution. The actual probe response - i.e., the surface charge to voltage transfer function, is accounted for in post-processing. This is accomplished via an Inverse Wien Filter - a technique often applied in image processing - to deconvolve the probe response from the measured data. A commercially available electrostatic probe, the Trek 341B meter with a 3455ET probe, capable of measuring +/- 20 kV was compared to the in- house probe that is designed to cover a wider range from +/- 35 kV. A resolution better than the simple probe resolution is achieved through the distribution’s scanning voltage method and deconvolution. Applying repeated scans, surface charge decay was tracked on various polymer materials to determine the material and environment dependence; materials included were PA6, PTFE, and others. As an example of material dependence, samples charged to 20 kV at 65% humidity experience full charge decay in approximately 45 seconds for PA5 (152 mm dia.), while it took about 100 times longer for PTFE (51 mm dia.) to observe the same decay/redistribution of charge.

36 MATERIALS SCIENCE↗

Characterization of mechanical deformation in aluminum by optical second harmonic generation

Second harmonic generation (SHG), a laser-based surface technique, was used to characterize near-surface defect evolution associated with mechanical deformation in 2024-T3 aluminum due to tensile loads. Electron backscatter diffraction (EBSD) measurements were also performed to determine the geometrically necessary dislocation density before and after tensile deformation. Comparison of the EBSD and SHG results found that dislocation density and SHG signal are inversely proportional; tensile deformation in the plastic regime led to an increase in dislocation density and a decrease in SHG response. A similar trend was seen in the SHG results on as-received (unpolished) Al samples; SHG response decreases with extent of plastic deformation. Furthermore, this study suggests that SHG has potential as a nondestructive evaluation technique.

36 MATERIALS SCIENCE↗

Inverse prediction of PuO2 processing conditions using Bayesian seemingly unrelated regression with functional data

Over the past decade, a variety of innovative methodologies have been developed to better characterize the relationships between processing conditions and the physical, morphological, and chemical features of special nuclear material (SNM). Different processing conditions generate SNM products with different features, which are known as “signatures” because they are indicative of the processing conditions used to produce the material. These signatures can potentially allow a forensic analyst to determine which processes were used to produce the SNM and make inferences about where the material originated. This article investigates a statistical technique for relating processing conditions to the morphological features of PuO 2 particles. We develop a Bayesian implementation of seemingly unrelated regression (SUR) to inverse-predict unknown PuO 2 processing conditions from known PuO 2 features. Model results from simulated data demonstrate the usefulness of the technique. Applied to empirical data from a bench-scale experiment specifically designed with inverse prediction in mind, our model successfully predicts nitric acid concentration, while results for Pu concentration and precipitation temperature were equivalent to a simple mean model. Our technique compliments other recent methodologies developed for forensic analysis of nuclear material and can be generalized across the field of chemometrics for application to other materials.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Relocating microseismicity from downhole monitoring of the Decatur CCS site using a modified double-difference algorithm

SUMMARY The injection of CO2 at the Decatur carbon capture and storage site has generated significant microseismic activity, which occurs in distinct spatial clusters up to approximately 2.2 km from the primary injection well. Accurate and precise event locations are vital for the characterization of the microseismicity to help understand the reservoir response to the CO2 injection, whilst enabling the identification of minor faults and fractures below the resolution of conventional active seismic imaging. However, microseismic monitoring of fluid injection sites, such as Decatur, is often performed using a network of borehole sensors often from a single well. While these downhole sensors have excellent detection capabilities, their poor azimuthal coverage limits the ability to precisely determine event locations. We have developed a modified double-difference relocation algorithm suitable for both 1-D and 3-D velocity models, and which incorporates differential back azimuth observations to allow the benefits of the original double-difference algorithm to be applicable to a downhole microseismic monitoring setting. Applying the modified double-difference algorithm to the microseismicity at Decatur, we have successfully relocated 4293 events. The relocation included over 59 million observations for 757 285 event pairs, split across seven geographic regions. Despite the majority of observations being recorded in only two boreholes, with an almost identical azimuthal coverage, the results have shown to be reliable with significantly reduced residuals and low uncertainties associated with the final locations. We have analysed the residuals in terms of their association with each geographic region, data type, station and individual events, to fully appreciate their influence in the inversion and the fit of the data to the final set of event locations. For each region, the relocated seismicity has become less diffuse with improved clustering, and with newly visible linear features often orientated in a NE–SW direction. These results show the potential improvements that can be made to microseismic event locations recorded by a borehole network with a limited and variable azimuthal distribution.

58 GEOSCIENCES↗