Search NASA⌕ Search

SEARCH · Search NASA

Results for “capillary”

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 163 records · Page 9

A generalized analytical energy balance model for evaluating agglomeration from a binary collision of wet particles

Agglomeration of wet particles, i.e., particles coated with a thin liquid layer, is a common phenomenon in many processes like fluidized bed combustion of low rank fuels. The availability of an agglomeration model that can evaluate the outcome of a binary collision between wet particles differing in solid particle properties, liquid layer thicknesses, and initial collision (impact) speeds is essential for obtaining a comprehensive understanding on the existing processes experiencing wet particle agglomeration or for a successful development of new processes with high chances of wet particle agglomeration. This study presents a generalized agglomeration model on the basis of energy conservation before and after collision when colliding wet particles may differ in solid particle properties, liquid layer thicknesses, and impact speeds. The model was established based on the approximate values of energy losses that may happen during the collision. It incorporates body forces, solid-solid contacting, liquid capillary, and viscous contributions, as well as the liquid bridge volume effect. Predictions of the new model for collision outcomes of identical wet particles were like those from an analytical energy balance model developed recently by the group for identical wet particles. We also validated the new model by experimental data from literature. The results of a collision direction analysis indicated that the direction often has a minimal effect on the collision outcome in many practical scenarios. The results of Monte Carlo uncertainty analyses with the new model revealed that proper estimations of impact speed, under capillary limiting conditions, and thickness of coating layers and asperity heights, under viscous limiting conditions, are critical for the realistic prediction of collision outcomes at impact speeds close to critical impact speed, i.e., the minimum particle speed required for the particles to rebound.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Informing field-scale CO 2 storage simulations with sandbox experiments: The effect of small-scale heterogeneities

Small-scale heterogeneities can significantly affect the fate of the CO 2 plume and trapping during CO 2 migration. We conducted geologic carbon storage field-scale simulations to investigate the impact of small-scale heterogeneities on plume dynamics and trapping performance. Small-scale heterogeneities have been shown to increase the amount of trapped CO 2 during buoyancy-driven flow. The trapped CO 2 saturation is validated by previous sandbox experimental work during buoyancy-driven flow in realistic heterogeneous domains and is implemented through the critical CO 2 saturation parameter (i.e., the first non-zero value in the drainage CO 2 relative permeability curve). Depending on the type and degree of heterogeneity, various critical CO 2 saturation values are exhibited. Furthermore, we investigated the effect of small-scale heterogeneities when multiple capillary pressure models are employed. This study demonstrates that an increase in critical CO 2 saturation reduces the CO 2 plume size and lateral extent, accompanying an increase in residual trapping and a decrease in solubility trapping. Lastly, we show that independent of the capillary pressure model used, an increase in critical saturation leads to similar CO 2 plume dynamics distribution and trapping performance. These results emphasize the importance of quantifying the effect of small-scale heterogeneity as they affect the large-scale behavior of the CO 2 plume.

58 GEOSCIENCES↗

Numerical analysis of coalescence-induced bubble departure for enhanced boiling heat transfer

Boiling heat transfer plays a crucial role in a wide range of applications, such as power generation, refrigeration, electronics cooling, and pharmaceutics. Among the various factors that influence boiling heat transfer, the dynamics of vapor bubble nucleation, growth, and departure from the heated surface stand out as particularly important. An emerging phenomenon that can promote the departure of bubbles smaller than the Fritz diameter is coalescence-induced departure. If the dynamics of this process are fully understood, then surfaces can be engineered to promote faster bubble departure and substantially increase the performance of boiling heat transfer. Further, this work expands on published results by presenting a detailed numerical analysis of bubble coalescence and departure for a range of initial bubble diameters and size ratios between coalescing bubbles. Analysis of the results is focused on explaining how the release of surface energy and bubble surface dynamics lead to bubble departure, as well as fundamentally distinguishing capillary–inertial jumping and buoyant–inertial departure mechanisms across different bubble sizes and size ratios. The results show that both the initial sizes of the coalescing bubbles and the ratio between their sizes can determine whether the merged bubble will leave the surface through capillary–inertial jumping or buoyant departure. Below a certain bubble size, the release of surface energy by the merger is not sufficient to propel the merged bubble from the surface.

42 ENGINEERING↗

Derivation of physical equations for high-speed laser welding using large language models

It is challenging to formulate complex physical phenomena that occur in a manufacturing process, particularly when the available data are limited, rendering conventional data-driven approaches ineffective. This study aims to predict humping onset in high-speed laser welding by introducing a novel framework, namely text-to-equations generative pre-trained transformer (T2EGPT). This method leverages the capabilities of large language models (LLMs), in combination with sparse experimental data and enriched literature data, to derive an interpretable and generalizable equation for predicting humping initiation. By capturing key correlations among physical parameters, T2EGPT generates a compact and dimensionless expression that accurately predicts hump formation. The equation reveals that humping arises from the interplay between inertia-driven backward melt flow and capillary-driven surface stabilization, where inertial forces drive molten metal backward and capillary forces resist surface deformation. Furthermore, compared to traditional data-driven models, T2EGPT demonstrates enhanced predictive accuracy and cross-material transferability. More broadly, this study highlights the potential of LLMs to integrate textual information with data-driven discovery, enabling the extraction of physical laws in data-scarce scientific domains.

36 MATERIALS SCIENCE↗

Analysis of the laser welding keyhole using inline coherent imaging

Laser beam welding is a widely used fusion welding process in many industrial applications such as automotive, aerospace, energy, defense, and medical products. Industry has a fundamental need to model the laser welding process to minimize experimental testing and improve confidence in production welds. However, computational models attempting to predict weld formation are limited by an incomplete understanding of the beam-material interactions. As a consequence, these models do not accurately predict the mechanisms associated with laser weld formation. To improve the current state of prediction capabilities, it is vital to better detect/measure the physical aspects of the weld pool during high energy density welding. In this work, a novel real-time laser weld monitoring device using inline coherent imaging (ICI) was used to provide a fundamental understanding of laser weld formation via vaporization. The objective of this work was to investigate and quantify the relationship relating laser weld parameters and the vapor capillary (keyhole) through a state-of-the-art measurement technique. Bead-on-plate laser beam welds were produced with partial penetration on 304L stainless steel, 2205 duplex stainless steel, and Ti-6Al-4V. Keyhole monitoring was performed using a commercially available ICI system to collect keyhole penetration data in real time. These measurements were reconstructed to generate the vapor capillary shape at different welding parameters. Process parameters significantly influenced the keyhole shape and the keyhole root position relative to the process beam. Finally, the keyhole geometry showed distinct differences between the stainless steel alloys and Ti-6Al-4V.

2205↗

An experimental investigation on the characteristics of heat pipes with annular type composite wick structure

Here, an experimental investigation has been conducted on the hydraulic characteristics of annular type wick structures for heat pipes. An experimental facility which can measure porosity, permeability, and effective pore radius of the wick structures in a vacuum condition is established. Nine different types of multi-layered (6 layers in total) composite screen meshes are characterized. Based on the measurement results, a wick structure composed of one layer of 100×100 mesh, three layers of 400×400 mesh, and two layers of 60×60 mesh is determined to have the highest permeability to effective pore radius ratio $(K/r_{eff})$ and was selected as a targeted sample. The wick-to-wall gap effect on the hydraulic characteristics of annular type wick structure is also investigated by measuring the permeability and effective pore radius of the sample wick structure with varying gap widths. The result shows that the permeability increases as the gap increases from 0 mm to 1.2 mm. After a peak at 1.2 mm, the permeability decreases as the gap increases and converges to the value of the case measured without a wall structure. The effective pore radius becomes smaller as the gap increases, making a peak at a gap distance of 1.2 mm. This result implies that there is an optimal point in gap distance which is determined to be 1.2 mm for the selected composite mesh structure. A capillary limitation correlation with a multiplying factor which explains the enhanced performance of the wick due to the gap is suggested. An annular wick type heat pipe is constructed and tested. The capillary limitation of the heat pipe showed good agreement with the suggested correlation.

42 ENGINEERING↗

Preparation of illite coated geomaterial microfluidic surfaces: Effect of salinity and heat treatment

Unique characteristics of clay (e.g., high surface area, deficiency in positive charges) make any clay abundant shales a good source of adsorbed gas. However, the presence of clay in porous media can also affect shale rock's producibility. The clay-water interaction can cause significant swelling and fines migration in the formation and impair the overall hydrocarbon recovery from shale. Geomaterial micromodels, developed by functionalizing traditional glass or PDMS surfaces with geomaterials (e.g., calcite, Quartz, clay), can represent the physicochemical properties of the natural porous media and have been used during the last few years to understand and evaluate the solid-fluid physicochemical interactions. Here, this study focuses on developing a clay-coated geomaterial surface and investigates the effect of base fluid's salinity on clay adsorption to the glass surface. Glass capillary tubes and straight channel borosilicate glass micromodels are coated with Illite clay to represent the pore-scale clay chemistry of Caney Shale, a Mississippian unconventional play in Southern Oklahoma, USA. 10 wt% of Illite clay slurries made with brines of four different salinities are used to coat the glass-capillary tubes to understand the effect of salinity in clay adsorption on the glass surface and overall coating quality. To achieve a stable coating and evaluate the impact of heat treatment on coating's stability, straight channel glass flow cells coated with Illite clay are heat-treated at low (25 °C) and high temperature (125 °C), and the results are compared with an untreated coated surface. Flooding tests were carried out with brines of different salinities on the coated surface to evaluate the stability of the coating. The experimental results indicated a strong relationship between the brine's salinity and the adsorption of clay particles on the glass surface. An increase in brine concentration resulted in improved adsorption of clay particles on the glass surface. Experiments involving heat treating the glass surface following the coating demonstrated significant improvement in the stability of the coating.

02 PETROLEUM↗

Phenomena Identification and Ranking Table (PIRT) for heat pipes

Heat pipes are advanced passive thermal management devices that utilize phase change and capillary action to achieve efficient heat transfer. However, due to the complexity of the phenomena coupled in heat pipes, including capillary, phase change, turbulence, and compressibility effects, there are high uncertainties in the predictability of their operational regimes and performance. This PIRT exercise, conducted as a collaborative effort involving the Department of Energy (DOE) Microreactor Program (MRP), the Nuclear Regulatory Commission (NRC), and university partners systematically identifies, reviews, and prioritizes critical phenomena affecting the operation of heat pipes based on their importance and knowledge levels. Additional analyses and discussion are provided for phenomena with high importance and low knowledge, such as wick de-wetting, critical heat flux, contact angles, and pressure dynamics. The discussions included the recognizing challenges and proposing future research directions for both modeling and simulation and experimental efforts. Additionally, the report addresses phenomena with medium importance and low knowledge that could impact heat pipe operation during non-normal or transient operation, including frozen startup, laminar to turbulent transition, geyser boiling, wick priming, underfilling conditions, surface roughness of the wick, NCGs trapped in the wick, and the timescales of startup and shutdown. In conclusion, this comprehensive evaluation serves as a valuable resource for guiding future research and development efforts, supporting the successful integration of heat pipes into critical applications such as nuclear reactors, and contributing to the advancement of heat pipe technologies in safety-critical industries.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗

Rheological behavior and fiber spinning of polyacrylonitrile (PAN)/Carbon nanotube (CNT) dispersions at high CNT loading

Rheological studies play an important role in polymer processing including fiber spinning. In the current work, rheological behavior has been studied for polyacrylonitrile (PAN)/carbon nanotube (CNT) dispersions in dimethyl formamide (DMF), where CNT loading is as high as 15 wt% of the total solids (polymer + CNT). The presence of CNTs increased the elastic-like and shear thinning behavior of the dispersions. The terminal slope of the log G' versus log G" curve, which was used to characterize the homogeneity of solution or dispersion, is dependent on the CNT concentration but independent of temperature between 0 and 60 °C. PAN/CNT-10 (90 wt% PAN +10 wt% CNT) and PAN/CNT-15 were gel spun at different conditions, with small diameter capillary at high shear rate and large diameter capillary at relatively low shear rate. Additionally, for PAN/CNT-10, fiber spinning with lower molecular weight polymer showed increased continuous jetting time when compared to spinning with higher molecular weight polymer. In addition, the rheological behavior of PAN/cellulose nanocrystal (CNC) dispersion with 20–60 wt% CNC was also compared to that of the PAN/CNT dispersions. PAN/CNC and PAN/CNT share similarities in both rheological behavior and fiber spinning at high filler loading.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Stream lamination and rapid mixing in a microfluidic jet for X-ray spectroscopy studies

Microfluidic mixers offer new possibilities for the study of fast reaction kinetics down to the microsecond time scale, and methods such as soft X-ray absorption spectroscopy are powerful analysis techniques. These systems impose challenging constraints on mixing time scales, sample volume, detection region size and component materials. The current work presents a novel micromixer and jet device which aims to address these limitations. The system uses a so-called ‘theta’ mixer consisting of two sintered and fused glass capillaries. Sample and carrier fluids are injected separately into the inlets of the adjacent capillaries. At the downstream end, the two streams exit two micron-scale adjoining nozzles and form a single free-standing jet. The flow-rate difference between the two streams results in the rapid acceleration and lamination of the sample stream. This creates a small transverse dimension and induces diffusive mixing of the sample and carrier stream solutions within a time scale of 0.9 microseconds. The reaction occurs at or very near a free surface so that reactants and products are more directly accessible to interrogation using soft X-ray. We use a simple diffusion model and quantitative measurements of fluorescence quenching (of fluorescein with potassium iodide) to characterize the mixing dynamics across flow-rate ratios.

Huyke, Diego A. (ORCID:0000000283356613)↗

Rapid Droplet Sampling Interface for Low-Volume, High-Throughput Mass Spectrometry Analysis

Here, we present a rapid droplet sampling interface (RDSI) electrospray ionization mass spectrometry (ESI-MS) system as a high-throughput, low-volume, noncontact, and minimal-carryover approach for characterization of liquids. Liquid characterization was achieved by combining droplet ejection with an open-face microflow capillary with a 2.5 μL/min continuous flow of carrier solvent. Through this implementation, single 0.3 nL droplets containing the analyte effectively mix with 4–8 nL of carrier solvent and create a combined electrospray plume. The carrier solvent continuously cleaned the system, eliminating carryover. A sampling rate of 5 Hz was achieved for droplets containing 1 μM propranolol or 5 μM leu-enkephalin with each droplet fully baseline-resolved (138 ± 32 ms baseline-to-baseline). Using a SCIEX API4000 mass spectrometer, a lower limit of quantification (LLOQ) of propranolol was 15 nM, corresponding to 1.16 fg of propranolol in the droplet, and was linear across 3 orders of magnitude. Further, quantitation could be achieved by adding an isotopically labeled internal standard, as done in conventional ESI. Signal transients were faster than the acquisition speed of the mass spectrometer, resulting in artificially high reproducibility of 15–30% RSD droplet-to-droplet. Analyte–solvent mixing ratios could be controlled by adjusting droplet positioning along the open-face capillary with an optimal position about 0.4 mm from the tip end. The range of analyte coverage was exemplified by measures of peptides and drugs in methanol, water, and buffer solutions. In a comparison to the Open Port Sampling Interface (OPSI) implemented on the same system, the RDSI had 78× greater sensitivity, 6× greater throughput and used significantly less carrier solvent.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Optimization of Delignified Wood Rods for Interfacial Solar Evaporation

Interfacial solar evaporation using three-dimensional evaporator materials has shown promise to achieve high evaporation rates exceeding the photothermal limit for water treatment and resource recovery processes. However, challenges remain in optimizing the material geometry to balance the vertical capillary uptake of water and the rate of evaporation to achieve both stable and high evaporation rates. Delignified wood can serve as a highly porous and hydrophilic substrate material to achieve high evaporation rates. In this work, we designed a suspended biochar-coated delignified wood rod evaporator and investigated the influence of its geometric parameters on its evaporation performance. The height of the evaporator exposed to air for evaporation was observed to be an important parameter in governing the evaporation rate, as shorter evaporators limited the available surface area for evaporation, but taller evaporators could not achieve sufficient rates of capillary uptake to maintain stable evaporation rates over 24 h under 1-sun illumination (1 kW m –2 ). In conclusion, the exposed height of the evaporator was optimized to achieve an average hourly evaporation rate of 3.05 ± 0.23 kg m –2 h –1 , which could be maintained in saline solutions containing up to 5 wt % NaCl.

Coating materials↗

Topology-Dependent Polymer Stretching and Scission in Solution at Extreme Shear Rates

There has been significant interest in the engineering of polymer topology to control rheology and mechanical stability in dilute solutions for applications involving extreme shear rate flows. However, methods to experimentally probe properties at relevant shear rates (≈ 10 4 –10 6 s –1 ) have remained ex situ, obscuring access to measures of polymer deformation and rheology that would otherwise provide mechanistic insight into the topology-dependent properties that control their behavior in extreme shear flows. In this study, we used novel in situ small angle neutron scattering measurements in a capillary rheometer (capillary rheo-SANS) to simultaneously measure solution viscosities and polymer deformations in high shear on a series of chemically homologous topology-defined polymers including linear, randomly branched, and star-shaped molecules. We demonstrate that differences in the onset of chain stretching and shear thinning of these polymers in dilute solution are controlled primarily by differences in their molecular relaxation time. These differences correlate with differences in chain scission inferred from ex situ measurements at more extreme shear rates. Together, the results demonstrate a direct coupling between chain deformation and scission, and suggest that the dominant effect of branching as a means to impart resilience against mechanical degradation is through differences in relaxation dynamics due to branching. We anticipate that these results will provide key insights to engineer topology-controlled polymers for rheological modification, mechanical stability, and controlled mechano-chemistry.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Design and Evaluation of a Tethered, Open Port Sampling Interface for Liquid Extraction-Mass Spectrometry Chemical Analysis

Presented in this work is a tethered, liquid-extraction-sampling interface designed for the mass spectrometric surface sampling/analysis of 3D objects. The tethered, open port sampling interface (TOPSI) incorporates a vacuum line between the sampling probe and ionization source, which enables the ability for an extended, tethered sample transfer line. Herein, several designs of the hand-held TOPSI are presented and evaluated on the basis of the analytical metrics of analyte transport time, peak width, and analyte sensitivity. The best analytical metrics were obtained with capillary flow resistances arranged in a particular order and the vacuum region set at 6.2 kPa. This TOPSI design incorporated a transfer capillary 1 m in length, while retaining a fast analyte transport time (12 s), short signal peak width (5 s baseline-to-baseline), and high analyte signal at 90% of that obtained with a regular open port sampling interface (OPSI). The hand-held TOPSI was demonstrated for the characterization of extracted small molecules and metabolites from the surface of mint and rosemary leaves.

47 OTHER INSTRUMENTATION↗

Link Between CO 2 –Induced Wettability and Pore Architecture Alteration

Changes in pore (throat) size, surface roughness and mineralogy induced by supercritical-CO 2 -water-rock reactions impact petrophysical properties such as porosity, permeability, and especially wettability. Herein, we show that these changes directly impact relative permeability and capillary pressure curves, a fact rarely studied in the literature. In this work, we show that CO 2 contact angle changes emerge after Madison Limestone samples were soaked for 400 hours in CO 2 -enriched brine. Coreflooding results show that the water production rate and cumulative water production increased after the rock was exposed to carbonic acid. Moreover, the mercury capillary pressure decreased in meso- and macro-pores, indicating the increase of size in these pores due to reactions. This compounded wettability and pore network alteration can directly affect CO 2 injectivity, migration and storage capacity. Furthermore, this fundamental insight into CO 2 geological storage processes should aid practitioners to reduce uncertainties in forecasting CO 2 distribution via injection simulation.

58 GEOSCIENCES↗

Three-Dimensional Free Gas Flow Focuses Basin-Wide Microbial Methane to Concentrated Methane Hydrate Reservoirs in Geological System

Here, we present a systematic model that links the generation, migration, phase partitioning, and accumulation of methane into a closed loop as the sediment is deposited from the seafloor and buried through the base of hydrate stability zone (BHSZ). In our model, methane is generated by biodegradation of organic carbon in muds. Hydrate does not form and methane is not trapped until a coarse-grained layer is deposited, because the small pores prevent hydrate formation in muds. Instead, methane diffuses into sands/silts where methane solidifies into hydrate. As hydrate-bearing sands/silts pass through the BHSZ during sediment burial, methane hydrate dissociates, and releases free gas. The released and the newly generated free gas below the BHSZ concentrates into a vertical/dipping zone with low capillary entry pressure and high permeability and flows upward driven its buoyancy. When free gas reaches the hydrate stability zone (HSZ), capillary forces drive free gas to flow laterally, preferentially enter sands/silts, feed hydrate growth, and elevate hydrate saturation. With three-dimensional focused free gas flow, microbial methane that is generated from a much larger fetch area of the entire basin, both above and below the BHSZ, is concentrated into coarse-grained layers at structural closures for hydrate formation. Our model illustrates how geological evolution, microbial methane generation, and gas flow by buoyancy couple to generate concentrated hydrate deposits in geological system. These insights can be used to explore for high-concentration methane hydrate and are important for understanding the methane budget and carbon cycle under the seafloor.

58 GEOSCIENCES↗

Diffusion-to-Imbibition Transition in Water Sorption in Nanoporous Media: Theoretical Studies

The ability to predict multiphase fluid transport in nanoporous rocks such as shales is critical for many geoscience applications, for example unconventional hydrocarbon production, geologic carbon sequestration, and nuclear waste disposal. When the pore sizes approach nanoscales, the impact of the molecular interaction forces between fluids and solids becomes increasingly important. These forces can alter macroscopic fluid phase behavior and control transport. Recent experimental studies have shown that capillary condensation and subsequent imbibition of liquid water can occur in hydrophilic nanoporous media even if the vapor phase is at a critical relative humidity (rhcrit) well below vapor saturation. This study presents a theoretical investigation of the processes controlling adsorption, capillary condensation and imbibition in nanoporous media, using the square-gradient classical density functional theory. The proposed theoretical model explicitly includes the relevant interaction forces among fluids and solids in macroscopic porous media. Application of the model to a relative-humidity-controlled water adsorption experiment is presented to demonstrate the impact of water-pore wall attractive forces on multiphase water behavior in a hydrophilic silicon nanoporous medium. Here, the model represents well the measured time-dependent evolution of the water imbibition front inside the nanoporous medium and also explains the diffusion-like water transport regimes observed at rh < rh crit and the imbibition-like flow regimes observed at rh > rh crit . The study furthermore gives an insight on hysteresis phenomenon in adsorption and desorption isotherms.

58 GEOSCIENCES↗

Impacts of Pore Network-Scale Wettability Heterogeneity on Immiscible Fluid Displacement: A Micromodel Study

The mixed-wet nature of reservoir formations imposes a wide range of rock wettability from strong resident-fluid wetting to strong invading-fluid wetting. The characteristics of two-phase flow in porous media composed of mixed-wetting surfaces remain poorly understood. In this study, we investigated the displacement of resident ethylene glycol (EG) by hexane in two mixed-wet micromodels of identical 2.5-D geometry heterogeneity, with uniformly or heterogeneously distributed patches strongly wetting to hexane. These patches are mixed among pores with unaltered EG-wetting surfaces. Along with control tests in the originally EG-wet micromodel, here we show the classic fingering and transitions in flow regimes at $logCa$ (capillary number) from -7.2 to -3.9. Moreover, pore-scale distributions of wettability and their spatial correlation influence displacement efficiency. In the two mixed-wet micromodels, we found (a) an increase of steady-state hexane saturation at the end of experiments by up to 0.12 in the capillary fingering regime and a decrease of at most by 0.06 in the viscous fingering regime, compared to the EG-wet micromodel, and (b) dispersed and fragmented hexane distribution after displacement. Brine drainage during supercritical CO 2 (scCO 2 ) injections in these micromodels occurs with lower wettability contrasts, and under similar viscosity ratios and interfacial tensions resulted in higher displacement efficiency relative to displacement of EG by hexane. While mixed-wettability can enhance displacement efficiency compared to uniform wettability, the dynamics of immiscible fluids in strong mixed-wet reservoirs are expected to be less pronounced in contributing to the efficiency of geological CO 2 sequestration, oil recovery, and remediation of hydrocarbon-contaminated aquifers.

2.5-D micromodel↗