Search NASA⌕ Search

SEARCH · Search NASA

Results for “Natural Fractures”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

High Thermoelectric Performance in 2D Sb 2 Te 3 and Bi 2 Te 3 Nanoplate Composites Enabled by Energy Carrier Filtering and Low Thermal Conductivity

Thermoelectrics are an important class of materials with great potential in alternative energy applications. In this study, two-dimensional (2D) nanoplates of the layered chalcogenides, Sb 2 Te 3 and Bi 2 Te 3 , are synthesized and composites of the two are investigated for their thermoelectric properties. The two materials, Sb 2 Te 3 and Bi 2 Te 3 , were synthesized as hexagonal, 2D nanoplates via a colloidal polyol route. The as-synthesized Sb 2 Te 3 and Bi 2 Te 3 vary drastically from one another in their lateral and vertical dimensions as revealed by scanning electron microscopy and atomic force microscopy. The single crystalline nanoplate nature is deduced by high-resolution transmission electron microscopy and selected area electron diffraction. Nanoplates have well-defined hexagonal facets as seen in the scanning and transmission electron microscopy images. The nanoplates were consolidated as an anisotropic nanostructured pellet via spark plasma sintering. Preferred orientation observed in the powder X-ray diffraction pattern and scanning electron microscopy images of the fractured pellets confirm the anisotropic structure of the nanoplates. Thermoelectric properties in the parallel and perpendicular directions were measured, revealing strong anisotropy with a significant reduction to thermal conductivity in the perpendicular direction due to increased phonon scattering at nanoplate interfaces. All compositions, except that of the 25% Bi 2 Te 3 nanoplate composite, behave as degenerate semiconductors with increasing electrical resistivity as the temperature increases. The Seebeck coefficient is also increased dramatically in the nanocomposites, the highest reaching 210 μV/K for 15% Bi 2 Te 3 . The increase in Seebeck is attributed to energy carrier filtering at the nanoplate interfaces. Overall, these enhanced thermoelectric properties lead to a drastic increase in the thermoelectric performance in the perpendicular direction, with zT ~ 1.26, for the 15% Bi 2 Te 3 nanoplate composite at 450 K.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Machine Learning Application to Assess Occurrence and Saturations of Methane Hydrate in Marine Deposits Offshore India

Artificial Neural Networks (ANN) were used to assess methane hydrate occurrence and saturation in marine sediments offshore India. The ANN analysis classifies the gas hydrate occurrence into three types: methane hydrate in pore space, methane hydrate in fractures, or no methane hydrate. Further, predicted saturation characterizes the volume of gas hydrate with respect to the available void volume. Log data collected at six wells, which were drilled during the India National Gas Hydrate Program Expedition 02 (NGHP-02), provided a combination of well log measurements that were used as input for machine learning (ML) models. Well log measurements included density, porosity, electrical resistivity, natural gamma radiation, and acoustic wave velocity. Combinations of well logs used in the ML models provide good overall balanced accuracy (0.79 to 0.86) for the prediction of the gas hydrate occurrence and good accuracy (0.68 to 0.92) for methane hydrate saturation prediction in the marine accumulations against reference data. The accuracy scores indicate that the ML models can successfully predict reservoir characteristics for marine methane hydrate deposits. In conclusion, the results indicate that the ML models can either augment physics-driven methods for assessing the occurrence and saturation of methane hydrate deposits or serve as an independent predictive tool for those characteristics.

58 GEOSCIENCES↗
DOE OSTI2025

Processes in Salt Repositories for Radioactive Waste Disposal

This document summarizes the key processes (thermal, hydrological, mechanical, and chemical; THMC) impacting the features of a deep geological repository for radioactive waste in salt. Some processes are natural and on-going whether the repository is there or not, and other processes are driven by the perturbation associated with the repository. The features considered here include both engineered and natural components of the repository system. The engineered barrier system (EBS) in a salt repository is quite different from those implemented for a repository in clay or crystalline rocks, because it is comprised mostly of granular salt and salt-compatible cements, rather than bentonite. When compared to other rocks (i.e., silicates), salt has unique properties that make it an excellent potential host rock. Openings and fractures in salt creep closed readily. Salt has high thermal conductivity, which can reduce peak temperatures. Additionally, far away from the excavations the porosity of salt is unconnected, which leads to essentially zero advective or diffusive transport. The small amount of hypersaline brine occurring in salt minimizes microbial activity, reduces colloid-assisted transport, and eliminates in-package criticality (i.e., chloride is a neutron poison). At the end of the report, we present a brief outline for a potential salt repository, including considerations avoided in previous repository disposal concepts. We propose considering higher-temperature processes in future disposal concepts, rather than trying to minimize the thermal perturbation of the repository. Since hot salt is drier, a dry repository would limit corrosion, gas generation, and solute transport. Openings and fractures creep shut faster in hot salt. Therefore, higher temperatures could be seen as beneficial, rather than something to minimize, through increased spacing between waste packages (increasing repository costs).

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

A review of spatial Markov models for predicting pre-asymptotic and anomalous transport in porous and fractured media

Heterogeneity across a broad range of scales in geologic porous media often manifests in observations of non-Fickian or anomalous transport. While traditional anomalous transport models can successfully make predictions in certain geological systems, increasing evidence suggests that assumptions relating to independent and identically distributed increments constrain where and when they can be reliably applied. A relatively novel model, the Spatial Markov model (SMM), relaxes the assumption of independence. The SMM belongs to the family of correlated continuous time random walks and has shown promise across a wide range of transport problems relevant to natural porous media. It has been successfully used to model conservative as well as more recently reactive transport in highly complex flows ranging from pore scales to much larger scales of interest in geology and subsurface hydrology. In this review paper we summarize its original development and provide a comprehensive review of its advances and applications as well as lay out a vision for its future development.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Toughening Brittle Bio‐P3HB with Synthetic P3HB of Engineered Stereomicrostructures

Abstract Poly(3‐hydroxybutyrate) (P3HB), a biologically produced, biodegradable natural polyester, exhibits excellent thermal and barrier properties but suffers from mechanical brittleness, largely limiting its applications. Here we report a mono‐material product design strategy to toughen stereoperfect, brittle bio or synthetic P3HB by blending it with stereomicrostructurally engineered P3HB. Through tacticity ([ mm ] from 0 to 100 %) and molecular weight ( M n to 788 kDa) tuning, high‐performance synthetic P3HB materials with tensile strength to ≈30 MPa, fracture strain to ≈800 %, and toughness to 126 MJ m −3 (>110× tougher than bio‐P3HB) have been produced. Physical blending of the brittle P3HB with such P3HB in 10 to 90 wt % dramatically enhances its ductility from ≈5 % to 95–450 % and optical clarity from 19 % to 85 % visible light transmittance while maintaining desirably high elastic modulus (>1 GPa), tensile strength (>35 MPa), and melting temperature (160–170 °C). This P3HB‐toughening‐P3HB methodology departs from the traditional approach of incorporating chemically distinct components to toughen P3HB, which hinders chemical or mechanical recycling, highlighting the potential of the mono‐material product design solely based on biodegradable P3HB to deliver P3HB materials with diverse performance properties.

Zhang, Zhen↗
DOE OSTI2023

Toughening Brittle Bio‐P3HB with Synthetic P3HB of Engineered Stereomicrostructures

Abstract Poly(3‐hydroxybutyrate) (P3HB), a biologically produced, biodegradable natural polyester, exhibits excellent thermal and barrier properties but suffers from mechanical brittleness, largely limiting its applications. Here we report a mono‐material product design strategy to toughen stereoperfect, brittle bio or synthetic P3HB by blending it with stereomicrostructurally engineered P3HB. Through tacticity ([ mm ] from 0 to 100 %) and molecular weight ( M n to 788 kDa) tuning, high‐performance synthetic P3HB materials with tensile strength to ≈30 MPa, fracture strain to ≈800 %, and toughness to 126 MJ m −3 (>110× tougher than bio‐P3HB) have been produced. Physical blending of the brittle P3HB with such P3HB in 10 to 90 wt % dramatically enhances its ductility from ≈5 % to 95–450 % and optical clarity from 19 % to 85 % visible light transmittance while maintaining desirably high elastic modulus (>1 GPa), tensile strength (>35 MPa), and melting temperature (160–170 °C). This P3HB‐toughening‐P3HB methodology departs from the traditional approach of incorporating chemically distinct components to toughen P3HB, which hinders chemical or mechanical recycling, highlighting the potential of the mono‐material product design solely based on biodegradable P3HB to deliver P3HB materials with diverse performance properties.

Chemistry↗
DOE OSTI2020

Cooperative Research and Development Agreement between National Energy Technology Laboratory and Ohio University [Abstract]

Ohio University and NETL will collaborate in the area of hydraulic fracturing and chemical treatment of gas and oil producing subsurface formations. The collaboration will provide Ohio University access to the J200 LIBS measurement machine, as well as training, facilities, and experts at NETL. NETL will gain access to Ohio University’s Point Pleasant formation shale samples and analysis results. The goal of this project is to determine the effects of chemical exposure on the face of the shale samples and how this exposure may affect hydrocarbon production. This work could help bring understanding to hydrocarbon movement within this formation and the effect hydraulic fracturing and other chemical treatments have on long-term hydrocarbon production. During this project the two parties will work together to generate and analyze experimental data from shale samples to develop a journal paper co-authored by both parties.

02 PETROLEUM↗
DOE OSTI2022

Stress-controlled fatigue of HfNbTaTiZr high-entropy alloy and associated deformation and fracture mechanisms

The stress-controlled fatigue tests are carried out at a stress ratio of 0.1 and a frequency of 10 Hz, and span both low-cycle and high-cycle regimes by varying the applied stress amplitudes. The high-cycle fatigue regime gives a fatigue strength of 497 MPa and a fatigue ratio of 0.44. At equivalent conditions, the alloy's fatigue strength is greater than all other high-entropy alloys (HEAs) with reported high-cycle fatigue data, dilute body-centered cubic alloys, and many structural alloys such as steels, titanium alloys, and aluminum alloys. Through in-depth analyses of crack-propagation trajectories, fracture-surface morphologies and deformation plasticity by means of various microstructural analysis techniques and theoretical frameworks, the alloy's remarkable fatigue resistance is attributed to delayed crack initiation in the high-cycle regime, which is achieved by retarding the formation of localized persistent slip bands, and its good resistance to crack propagation in the low-cycle regime, which is accomplished by intrinsic toughening backed up by extrinsic toughening. Additionally, the stochastic nature of the fatigue data is neatly captured with a 2-parameter Weibull model.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Mechanically-, Thermally, and Chemically-Robust High-Temperature Ceramic Composites

In an effort to lower the levelized cost of electricity generated by Concentrated Solar Power (CSP) plants, the heat-to-electricity conversion efficiency may be increased by operating such plants with higher turbine inlet temperatures using closed Brayton cycles with high-pressure, supercritical carbon dioxide (sCO 2 ) as the working fluid. For example, if the turbine inlet temperature is raised to ≥750°C, in combination with sCO2 closed Brayton cycles, the relative heat-to-electricity efficiency may be raised by more than 20% (compared to subcritical steam Rankine cycles operating with turbine inlet temperatures ≤550°C). The associated reduction in the cost of dispatchable electricity from such CSP plants (utilizing high-temperature thermal energy storage) would be an important breakthrough on the path towards direct competition with fossil-fuel-based plants. A key barrier to achieve such high-temperature efficient operation has been the limited thermomechanical performance of metal alloys used in compact, primary printed circuit- type heat exchangers (PCHEXs) for heat transfer to high-pressure sCO 2 . The maximum allowed stresses for the use of conventional stainless steels and nickel-based superalloys at high sCO 2 pressures (≥20 MPa) decline rapidly at temperatures ≥550°C. This project has been focused on demonstrating the desired high-temperature properties and attractive manufacturing characteristics of mechanically-, thermally-, and chemically-robust ceramic/metal composites (cermets), in order to allow for the use of such composites in primary PCHEXs for heat transfer to sCO 2 at ≥750°C and ≥20 MPa in CSP plants. Prior work with high-temperature, co-continuous carbide/refractory metal cermets, such as ZrC/W composites, has already demonstrated: i) the attractive combination of high-temperature properties exhibited by such materials (e.g., high values of stiffness, failure strength, and thermal conductivity at 800°C), and ii) cost-effective processes for manufacturing PCHEX plates comprised of such cermets with tailorable channel patterns. However, ZrC/W cermets were not found to be inherently corrosion resistant in sCO 2 at 750°C and 20 MPa. While the application of a Cu coating to ZrC/W surfaces, along with a modest addition (50 ppm) of CO to CO 2 , rendered such composites resistant to corrosion for 1000 h in such sCO 2 -based fluids at 750°C, such additional steps would increase the design complexity of such cermet-based components. An alternative approach (the focus of this project) is to develop mechanically-robust, cost- effective cermets that are inherently resistant to high-temperature oxidation in air and CO 2 . The overall goals of this project have been to demonstrate: i) the oxidation resistance of at least one cermet at 750°C in air and in CO 2 (i.e., with a projected annual corrosion ≤30 μm), ii) the mechanically robust nature of at least one cermet at 750°C (i.e., with an average failure strength ≥200 MPa), and iii) that the selected cermet can be manufactured via low-cost forming techniques.

14 SOLAR ENERGY↗
DOE OSTI2021

Oil and Gas Produced Water Reuse: Opportunities, Treatment Needs, and Challenges

Advances in water treatment technologies paired with potential restrictions on oil and gas (O&G) produced water disposal could incentivize the beneficial reuse of treated produced water in the O&G industry. However, the remote nature of O&G operations limits the applicability of many of these solutions, which may be spatially inefficient, require operator supervision, or are ill-suited for the complex nature of produced water. Furthermore, the responsible, sustainable reuse of produced water as an alternative water source requires standardized analytical techniques for characterizing and determining the toxicity of treated produced water and improving our understanding of the fate and transport of various constituents. In the past decade, we made little progress in economically treating produced water for beneficial reuse outside of oilfield operations; the sole major breakthrough has been in the development of salt-tolerant fracturing chemicals that allow for reuse of produced water for fracking operations. Guided research should assist in the development of fit-for-purpose solutions to maximize the reuse of treated produced water. This is exemplified by the case studies presented here that detail currently operating treatment facilities for reclamation and reuse of produced water.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2021

Fracture Sustainability in Enhanced Geothermal Systems: Experimental and Modeling Constraints

Enhanced geothermal systems (EGS) offer the potential for a much larger energy source than conventional hydrothermal systems. Hot, low-permeability rocks are prevalent at depth around the world, but the challenge of extracting thermal energy depends on the ability to create and sustain open fracture networks. Laboratory experiments were conducted using a suite of selected rock cores (granite, metasediment, rhyolite ash-flow tuff, and silicified rhyolitic tuff) at relevant pressures (uniaxial loading up to 20.7 MPa and fluid pressures up to 10.3 MPa) and temperatures (150–250 °C) to evaluate the potential impacts of circulating fluids through fractured rock by monitoring changes in fracture aperture, mineralogy, permeability, and fluid chemistry. Because a fluid in disequilibrium with the rocks (deionized water) was used for these experiments, there was net dissolution of the rock sample: this increased with increasing temperature and experiment duration. Thermal-hydrological-mechanical-chemical (THMC) modeling simulations were performed for the rhyolite ash-flow tuff experiment to test the ability to predict the observed changes. These simulations were performed in two steps: a thermal-hydrological-mechanical (THM) simulation to evaluate the effects of compression of the fracture, and a thermal-hydrological-chemical (THC) simulation to evaluate the effects of hydrothermal reactions on the fracture mineralogy, porosity, and permeability. Furthermore, these experiments and simulations point out how differences in rock mineralogy, fluid chemistry, and geomechanical properties influence how long asperity-propped fracture apertures may be sustained. Such core-scale experiments and simulations can be used to predict EGS reservoir behavior on the field scale.

15 GEOTHERMAL ENERGY↗
DOE OSTI2023

Compatibility of mesophase pitch and linear low-density polyethylene for low-cost carbon fiber

Mesophase pitch and linear low-density polyethylene (LLDPE) were blended by shear mixing and extruded producing precursor fibers which were then converted to carbon fibers (CFs). LLDPE additives agglomerate within the isotropic spheres rather than uniformly dispersing throughout the entire matrix. The blends possess different thermal properties and precursor tensile properties evident by two-phase tensile stress nature not present in precursor fibers produced from pure pitch. The resultant CFs were oxidized at 280 °C for 4 h and carbonized at 1000 °C for 15 min. This process created pitch CFs with a tensile strength of 1.23 GPa and a Young’s modulus of 99 GPa. All CFs produced from blends of pitch and LLDPE possessed a tensile strength of 0.63 GPa and less. These CFs showed a microstructure slightly different from the pure pitch CFs with fractured CF cross sections revealing voids present in all blended CFs due to burnt off LLDPE stemming from carbonization. Further study into compatibilizers and new conversion processes is required to successfully blend LLDPE with mesophase pitch for the purpose of producing a new low-cost CF.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Arctic Aerosol Sources and Mixing States Field Campaign Report

The Arctic is warming at a faster rate than anywhere else on Earth, with rapidly shrinking sea ice extent transforming the region. Depending on chemical and physical properties, atmospheric aerosols directly scatter and/or absorb radiation, serve as cloud droplet and/or ice crystal nuclei, and/or reduce the reflectiveness of the snow surface, thereby altering the atmospheric energy budget. There is a wide spread in the magnitude of simulated arctic aerosol radiative forcing, and significant differences in aerosol concentration levels and seasonal cycles often exist between models and observations. Increasing local natural and anthropogenic emissions are significant, with uncertain climate impacts due to complex feedbacks. Model evaluations, however, are limited by the dearth of arctic aerosol observations available and an inadequate understanding of arctic aerosol processes. The majority of previous arctic aerosol observations have been made through intensive spring/summer field campaigns, with few intensive measurement studies focused on the fall-winter transition, a period when freeze-up is occurring later and thinning sea ice is resulting in wintertime ice fracturing. Aerosol monitoring at arctic coastal stations has provided knowledge of long-term seasonal trends in aerosol concentrations. The completed ARM field campaign addresses observational and knowledge gaps through detailed aerosol size and chemical composition measurements during the fall-winter transition in the coastal Arctic and through the entire annual cycle in the central Arctic. The observations are improving our understanding of the sources and processes controlling the aerosol population in the rapidly changing Arctic.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Physics-embedded inverse analysis with algorithmic differentiation for the earth’s subsurface

Abstract Inverse analysis has been utilized to understand unknown underground geological properties by matching the observational data with simulators. To overcome the underconstrained nature of inverse problems and achieve good performance, an approach is presented with embedded physics and a technique known as algorithmic differentiation. We use a physics-embedded generative model, which takes statistically simple parameters as input and outputs subsurface properties (e.g., permeability or P-wave velocity), that embeds physical knowledge of the subsurface properties into inverse analysis and improves its performance. We tested the application of this approach on four geologic problems: two heterogeneous hydraulic conductivity fields, a hydraulic fracture network, and a seismic inversion for P-wave velocity. This physics-embedded inverse analysis approach consistently characterizes these geological problems accurately. Furthermore, the excellent performance in matching the observational data demonstrates the reliability of the proposed method. Moreover, the application of algorithmic differentiation makes this an easy and fast approach to inverse analysis when dealing with complicated geological structures.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Hydroxylpyromorphite, a mineral important to lead remediation: Modern description and characterization

Abstract Hydroxylpyromorphite, Pb5(PO4)3(OH), has been documented in the literature as a synthetic and naturally occurring phase for some time but has not previously been formally described as a mineral. It is fully described here for the first time using crystals collected underground in the Copps mine, Gogebic County, Michigan. Hydroxylpyromorphite occurs as aggregates of randomly oriented hexagonal prisms, primarily between about 20–35 μm in length and 6–10 μm in diameter. The mineral is colorless and translucent with vitreous luster and white streak. The Mohs hardness is ~3½–4; the tenacity is brittle, the fracture is irregular, and indistinct cleavage was observed on {001}. Electron microprobe analyses provided the empirical formula Pb4.97(PO4)3(OH0.69F0.33Cl0.06)Σ1.08. The calculated density using the measured composition is 7.32 g/cm3. Powder X-ray diffraction data for the type material is compared to data previously reported for hydroxylpyromorphite from the talc mine at Rabenwald, Austria, and from Whytes Cleuch, Wanlockhead, Scotland. Hydroxylpyromorphite is hexagonal, P63/m, at 100 K, a = 9.7872(14), c = 7.3070(10) Å, V = 606.16(19) Å3, and Z = 2. The structure [R1 = 0.0181 for 494 F>4σ(F) reflections] reveals that hydroxylpyromorphite adopts a column anion arrangement distinct from other members of the apatite supergroup due to the presence of fluorine and steric constraints imposed by stereoactive lone-pair electrons of Pb2+ cations. The F– anion sites are displaced slightly from hydroxyl oxygen anions, which allows for stronger hydrogen-bonding interactions that may in turn stabilize the observed column-anion arrangement and overall structure. Our modern characterization of hydroxylpyromorphite provides deeper understanding to a mineral useful for remediation of lead-contaminated water.

Geochemistry & Geophysics↗
DOE OSTI2021

Surfactants are Ineffective for Reducing Imbibition of Water-Based Fracturing Fluids in Deep Gas Reservoirs

Minimizing loss of injected hydraulic fracturing fluids into shale along fracture-matrix boundaries is desired because imbibed water restricts gas production and wastes valuable water resources. This problem has motivated the addition of surfactants into water-based hydraulic fracturing fluids in order to reduce the capillary driving force for imbibition. Here, we show that reduction in interfacial tension and wettability alteration has negligible ability to reduce imbibition in deep gas reservoirs. The effectiveness of altering capillary forces acting at the wetting front also depends on the injection pressure acting at the fracture-matrix boundary. The pressure at the interface between the fracture and the shale matrix is constrained between the reservoir pore pressure and formation pressure (rock fracture pressure, also known as breakdown pressure of the rock) and increases with depth to magnitudes that greatly exceed that of capillary pressures. Furthermore, the analyses presented here show that even maximum alteration of interfacial properties that result in strongly hydrophobic interactions between the fracturing fluid and reservoir rock is incapable of significantly reducing imbibition in deep reservoirs. Instead of using surfactants, this analysis points to decreases in wellbore shut-in pressures and shut-in times as practical options for reducing imbibition losses of water-based fluids.

03 NATURAL GAS↗
DOE OSTI2020

Mineral reaction kinetics constrain the length scale of rock matrix diffusion

Mass transport by aqueous fluids is a dynamic process in shallow crustal systems, redistributing nutrients as well as contaminants. Rock matrix diffusion into fractures (void space) within crystalline rock has been postulated to play an important role in the transient storage of solutes. The reacted volume of host rock involved, however, will be controlled by fluid-rock reactions. Here we present the results of a study which focusses on defining the length scale over which rock matrix diffusion operates within crystalline rock over timescales that are relevant to safety assessment of radioactive and other long-lived wastes. Through detailed chemical and structural analysis of natural specimens sampled at depth from an active system (Toki Granite, Japan), we show that, contrary to commonly proposed models, the length scale of rock matrix diffusion may be extremely small, on the order of centimetres, even over timescales of millions of years. This implies that in many cases the importance of rock matrix diffusion will be minimal. Additional analyses of a contrasting crystalline rock system (Carnmenellis Granite, UK) corroborate these results.

58 GEOSCIENCES↗
DOE OSTI2025

Invariant Forms of Dissolution Fingers

Dissolution of fractured and porous media introduces a positive feedback between fluid flow and reactant transport, leading to the emergence of pronounced, fingerlike channels. We investigate the formation of these structures using a microfluidic Hele-Shaw cell with a soluble bottom. Our experiments show that the shape of dissolution fingers is invariant and reveals itself over time as the fingers extend into the system. Here, by combining reactive-transport theory and conformal mapping techniques we derive these invariant forms. We relate these results to natural dissolution fingers in karst landscapes, and illustrate how to determine the groundwater flow rate responsible for their formation based on the finger shape.

Flows in porous media↗