Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fluid migration”

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

Regional-scale analysis of karst underground flow deduced from tracing experiments: examples from carbonate aquifers in Malaga province, southern Spain

Tracer concentration data from field experiments conducted in several carbonate aquifers (Malaga province, southern Spain) were analyzed following a dual approach based on the graphical evaluation method (GEM) and solute transport modeling to decipher flow mechanisms in karst systems at regional scale. The results show that conduit system geometry and flow conditions are the principal factors influencing tracer migration through the examined karst flow routes. Solute transport is mainly controlled by longitudinal advection and dispersion throughout the conduit length, but also by flow partitioning between mobile and immobile fluid phases, while the matrix diffusion process appears to be less relevant. The simulation of tracer breakthrough curves (BTCs) suggests that diffuse and concentrated flow through the unsaturated zone can have equivalent transport properties under extreme recharge, with high flow velocities and efficient mixing due to the high hydraulic gradients generated. Tracer mobilization within the saturated zone under low flow conditions mainly depends on the hydrodynamics (rather than on the karst conduit development), which promote a lower longitudinal advection and retardation in the tracer migration, resulting in a marked tailing effect of BTCs. The analytical advection-dispersion equation better approximates the effective flow velocity and longitudinal dispersion estimations provided by the GEM, while the non-equilibrium transport model achieves a better adjustment of most asymmetric and long-tailed BTCs. The assessment of karst underground flow properties from tracing tests at regional scale can aid design of groundwater management and protection strategies, particularly in large hydrogeological systems (i.e. transboundary carbonate aquifers) and/or in poorly investigated ones.

58 GEOSCIENCES↗

Nonwoven Fabric Uses and Prospects in Human Space Flight

The US space shuttle fleet has been flying for over 20 years. Although the shuttle operates in a unique exterior environment, the interior is intentionally made to be as close to the "normal" human environment as possible. The filtration needs of the shuttle are not substantially different from those of a large mobile home or camper, supporting the needs of a family of seven for up to two weeks. Therefore, most of the materials that are used to filter the air, water, and other fluids on the Shuttle are similar or identical to those employed in other sectors of the transportation industry. The only significantly different feature of the space environment is the unique "three-phase" nature of the air (with suspended liquids and solids ranging in size from aerosol droplets to binoculars). Such suspended debris contributes to the air filtration and waste management problem. Careful flow management and cleanliness practices help to mitigate the effect of debris, and liquid spills are rare, seldom making it to the filters. (It has been common on all spacecraft to look first for lost items on the air intake filters, since all objects ultimately migrate there in the flow. Liquids tend to seep rather than "spill", and so tend to aggregate in a ball near the source.) In addition to the basic fluids of the interior environment (water and water wastes, air, and its constituent supply gasses) the shuttle also has unfiltered fluid systems for Freon, hydrogen, helium, ammonia, hydraulic fluid, and propellants. Only the propellant system, owing to its uncommon chemistry, represents a fluid system that is not typical of household or medical applications. Careful external filtration prior to flight assures the cleanliness in these closed systems.

Bacon, Jack↗

Induced migration of non-magnetic particles and fabrication of metallic-based graded materials by applying a strong magnetic field

Strong magnetic fields e.g., 10 T are now frequently used during materials preparation. In the present paper, a method of metal-ceramics graded materials preparation is proposed by applying a strong magnetic field with a high gradient. Experimental and theoretical analyses are performed to investigate the effect of the strong magnetic field on the migration and interaction behavior of diamagnetic oxide particles in a liquid metal. The migration of micrometer or sub-micrometer sized particles is clearly enhanced by the magnetic field gradient. After being treated by the strong magnetic field, a self-assembly structure of the particles is achieved. Various factors such as the magnetic dipole-dipole interaction and chain-chain interaction, are governing the particles assembly. The present results provide direct evidence that nonmetallic particles of micrometer size or less can be manipulated in a conductive melt by applying a strong magnetic field and as a consequence graded materials can be prepared.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Experimental Results for Large Particle Separation from Non-Newtonian Slurries Using Tapered Bump Arrays

Lateral displacement arrays are useful for separating particles such as blood cells and sand from carrier fluids. These arrays consist of staggered posts, which allow smaller particles to follow streamlines and larger particles to flow around the posts and migrate to one side. This migration increases the particle concentration in one direction and depletes the particle concentration in the other direction allowing particle separation to occur. Experiments were conducted to separate large particles in non-Newtonian yield stress slurries using tapered bump arrays. The non-Newtonian slurry used was composed of a bentonite kaolin blend with the inclusion of larger diameter particles. These experiments were conducted to evaluate the performance of particle separation using a tapered array of staggered posts as configured in a bump array for non-Newtonian yield stress slurries, an application that has not been explored experimentally. A theory for large particle separation was developed for slurries that exhibit Bingham or Cross rheology. The results of these experiments are described. This information could be applied in industrial settings such as separation of particles from nuclear waste slurries including those to be processed at the Hanford site, where removing large particles from waste streams is important to processing.

slurry, solids removal, deterministic lateral disp↗

Evidence for organic synthesis in high temperature aqueous media - facts and prognosis

Hydrothermal systems are common along the active tectonic areas of the earth. Potential sites being studied for organic matter alteration and possible organic synthesis are spreading ridges, off-axis systems, back-arc activity, hot spots, volcanism, and subduction. Organic matter alteration, primarily reductive and generally from immature organic detritus, occurs in these high temperature and rapid fluid flow hydrothermal regimes. Hot circulating water (temperature range - warm to greater than 400 C) is responsible for these molecular alterations, expulsion and migration. Compounds that are obviously synthesized are minor components because they are generally masked by the pyrolysis products formed from contemporary natural organic precursors. The reactivity of organic compounds in hot water (200-350 C) has been studied in autoclaves, and supercritical water as a medium for chemistry has also been evaluated. This high temperature aqueous organic chemistry and the strong reducing conditions of the natural systems suggest this as an important route to produce organic compounds on the primitive earth. Thus a better understanding of the potential syntheses of organic compounds in hydrothermal systems will require investigations of the chemistry of condensation, autocatalysis, catalysis and hydrolysis reactions in aqueous mineral buffered systems over a range of temperatures from warm to greater than 400 C.

Simoneit, Bernd R. T.↗

Simulating the Effect of Fracture Connectivity on Repository Performance with GDSA Framework - 18589

This work examines the possible migration of radionuclides from a deep geologic nuclear waste repository sited in fractured crystalline host rock. The key safety concern is the potential effect on waste isolation of the inter-connectivity of the fracture network, which is primarily established by the temporal evolution of the temperature and stress fields at the time of original rock intrusion and cooling. Two end members are considered in these simulations, one with a high degree of connectivity to the biosphere, such that advective transport through the fracture and fault network controls radionuclide migration, and the other with a low degree of connectivity, such that slow diffusive transport through the crystalline rock matrix is the controlling process for radionuclide migration to the biosphere. Both end members have several areally extensive, high-transmissivity deformation zones, whose distance to the repository is an important factor for waste isolation capability. Uncertainties in fracture properties (transmissivity, orientation, radius) and fracture distribution also give rise to uncertainty in waste isolation capability (i.e., the overall connectivity of the repository with the surface biosphere). These types of natural system uncertainties are generally present in any geologic site-characterization program and thus represent an important factor in assessing repository performance in hard rock environments. The waste considered in these simulations are spent fuel rod assemblies from the U.S. commercial nuclear reactor fleet. Thermal output of the waste must be considered in the simulations and can influence the rate and timing of waste package failure, waste form degradation, and fluid flux due to thermal expansion around the repository horizon. To achieve a high-fidelity representation of radionuclide transport in fractures and rock matrix, combined with thermal energy transport and fluid flow in fractures and matrix, the mathematical model is solved numerically in a parallel high-performance computing (HPC) environment on a finite volume unstructured grid consisting of approximately 5 million cells, using Geologic Disposal Safety Assessment (GDSA) Framework (https://pa.sandia.gov), an open-source performance assessment tool for deep underground disposal of nuclear waste. GDSA Framework uses PFLOTRAN to solve the balance equations on a three-dimensional grid with heterogeneous properties, using multiple processors in a parallel configuration based on domain decomposition. The fracture networks in these simulations are originally generated as discrete fracture networks (DFNs), which are sets of two-dimensional planes distributed in a three-dimensional domain. The method used in GDSA Framework maps the stochastically generated DFN to an equivalent continuous porous medium (ECPM) domain that allows for the simulation of coupled heat flow, fluid flow, and radionuclide transport, including heat conduction through the matrix of the fractured rock, which is a process not easily modeled in a DFN representation. Computational efficiency is also greatly enhanced using the ECPM method, allowing for a realistic representation and analysis of uncertainties in a multi-realization performance assessment of a deep geologic repository. The effect of fracture connectivity on the waste isolation safety function, as brought to light by these GDSA Framework simulations, points to the importance of including a realistic representation of uncertainties in fracture properties and distribution (effectively, uncertainty in spatial heterogeneity) in repository safety assessment simulations. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Lubricated immersed boundary method in two dimensions

Many biological examples of fluid-structure interaction, including the transit of red blood cells through the narrow slits in the spleen and the intracellular trafficking of vesicles into dendritic spines, involve the near-contact of elastic structures separated by thin layers of fluid. Motivated by such problems, we introduce an immersed boundary method that uses elements of lubrication theory to resolve thin fluid layers between immersed boundaries. We demonstrate 2nd-order accurate convergence for simple two-dimensional flows with known exact solutions to showcase the increased accuracy of this method compared to the standard immersed boundary method. Motivated by the phenomenon of wall-induced migration, we apply the lubricated immersed boundary method to simulate an elastic vesicle near a wall in shear flow. We also simulate the dynamics of a vesicle traveling through a narrow channel and observe the ability of the lubricated method to capture the vesicle motion on relatively coarse fluid grids.

59 BASIC BIOLOGICAL SCIENCES↗

Characterization of induced seismicity patterns derived from internal structure in event clusters

Microseismicity induced by CO 2 injection at Decatur, Illinois, occurs in distinct clusters and shows no obvious correlation with the proceeding pressure front. Here, we analyze some of these clusters in more depth by using a waveform cross-correlation approach. With this approach we can associate about 1400 events from two clusters, with moment magnitudes between 1.1 and -1.7, with specific formations of much smaller vertical dimensions (tens of meters) than the depth resolution of traveltime-based event locations. The differentiation of reservoir and basement events, and the definition of subclusters by waveform correlation, rather than by location, helps to better analyze the spatiotemporal evolution of the events within a cluster. In the Decatur case, this is characterized by event migration from the reservoir into the adjacent basement. The spatial variation of Brune stress drop and Gutenberg b value exhibits signs of a fluid-driven triggering mechanism at the cluster level, revealing a punctual hydraulic connection between reservoir and basement, most likely associated with basement faults cutting into the reservoir. The observed clustering of microseismicity can thus be explained by the lateral heterogeneity of permeability and crustal strength and is overall consistent with a pressure-induced triggering mechanism. Hence, proper long-term risk mitigation for large-scale fluid injection close to the basement requires prior mapping of small subseismic basement-connected faults.

58 GEOSCIENCES↗

Nanoconfined Grain Boundaries Increase the Conductivity of Polycrystalline Molecular Crystals

Soft-solid molecular crystals consist of crystalline grains and fluid grain boundaries (GBs) that enhance the grain binding and transport of Li + ions between the grains. The total ionic conductivity consists of ion migration in both the grains and GBs. To unravel these contributions in adiponitrile (Adpn):LiPF 6 molecular crystals, the GB volume fraction was varied by changing the size of the crystals and the Adpn:LiPF 6 molar ratio. Molecular dynamics (MD) simulations indicate that ion motion was subdiffusive in the grains and “well-diffusive” in the GBs, with GBs characterized as disordered nanoconfined regions of higher charge carrier concentration (∼1 M) than in saturated Adpn:LiPF 6 solutions (0.04 M), and Li + ions predominantly solvated by cyano groups with few contact ion pairs. The diffusivity in the GBs is at least an order of magnitude higher than that in the crystalline grains. The emergent picture is the grains as a reservoir of ions that migrate to faster-conducting GBs.

Crystals↗

The role of fault-zone architectural elements and basal altered zones on pore pressure propagation and induced seismicity

Here, we used hydrogeologic models to assess how fault-zone properties promote or inhibit the downward propagation of fluid overpressures from a basal reservoir injection well (150 m from fault zone, Q = 5000 m 3 /day) into the underlying crystalline basement rocks. We varied the permeability of the fault-zone architectural components and a crystalline basement weathered layer as part of a numerical sensitivity study. Realistic conduit-barrier style fault zones effectively transmit elevated pore pressures associated with four years of continuous injection to depths of ~2.5 km within the crystalline basement while compartmentalizing fluid flow within the injection reservoir. The presence of a laterally continuous, relatively low-permeability altered/weathered basement horizon (k altered layer = 0.1 × k basement ) can limit the penetration depth of the pressure front to ~500 m. On the other hand, the presence of a discontinuous altered/weathered horizon that partially confines the injection reservoir without blocking the fault fluid conduit promotes downward propagation of pressures. Permeability enhancement via hydromechanical failure was found to increase the depth of early-time pressure front migration by a factor of 1.3 to 1.85. Dynamic permeability models may help explain seismicity at depths of >10 km such as is observed within the Permian Basin, NM.

58 GEOSCIENCES↗

Application of Generic Disposal System Models

This report describes specific GDSA activities in fiscal year 2015 (FY2015) toward the development of the enhanced disposal system modeling and analysis capability for geologic disposal of nuclear waste. The GDSA framework employs the PFLOTRAN thermal-hydrologic-chemical multi-physics code (Hammond et al., 2011) and the Dakota uncertainty sampling and propagation code (Adams et al., 2013). Each code is designed for massively-parallel processing in a high-performance computing (HPC) environment. Multi-physics representations in PFLOTRAN are used to simulate various coupled processes including heat flow, fluid flow, waste dissolution, radionuclide release, radionuclide decay and ingrowth, precipitation and dissolution of secondary phases, and radionuclide transport through the engineered barriers and natural geologic barriers to a well location in an overlying or underlying aquifer. Dakota is used to generate sets of representative realizations and to analyze parameter sensitivity.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Understanding Multicomponent Multiphase Hydrocarbon Fluid Interaction and Transport in Shale Matrix: From Nanoscale to Wellbore Production

One critical issue related to shale oil/gas production is a rapid decline in wellbore production and a low recovery rate. Therefore, maximizing wellbore production and extending the production life cycle are crucial for the sustainability of shale oil/gas production. Shale oil/gas production starts with creating a fracture network by injecting a pressurized fluid in a wellbore. The induced fractures are then held open by proppant particles. During production, oil and gas release from the shale matrix, migrate to nearby fractures, and ultimately reach a production wellbore. Given the relatively high permeability of the induced fractures, oil/gas release and transport in low-permeability shale matrix are likely a limiting step for long-term wellbore production. This project is aimed to (1) fundamentally understand the disposition and release of complex hydrocarbon mixtures in nanopore networks of shale matrix and their transport from low-permeability matrix to hydrofracking-induced fractures under various reservoir conditions ranging from black oil to dry gas and (2) use machine learning to upscale and integrate the nanometer-scale understanding into reservoir-scale model simulations. The work will help develop new stimulation strategies to enable efficient resource recovery from fewer and less environmentally impactful wells.

03 NATURAL GAS↗

rHealth One Demonstration Aboard ISS: A Microfluidic Bioanalyzer Based on Sheath/Hydrodynamic Focusing Flow Cytometry

The Exploration Medical Capability (ExMC) element aims to provide astronauts with the means for their own health monitoring, diagnosis, and treatment during exploration missions. As space flight ventures further from earth, the need for autonomous medical care increases under greater constraints on size, mass, and resources. One pillar for diagnosis that crew would be separated from is laboratory analysis. Even now on ISS, blood samples must be collected and returned to earth for testing. In response, ExMC is assessing how assays for hematology, bone health, radiation exposure etc. could be addressed through commercial-off-the-shelf (COTS) and Small Business Innovation Research (SBIR) funded bioanalyzers that are miniaturizing lab technology. Since missions will reach distances where there are no timely replacements, validation on the International Space Station (ISS) is a necessary part of that assessment. In partnership with NASA Johnson Space Center (JSC) Immunology Lab and the Research Operations and Integration (ROI) element, ExMC conducted a technology demonstration on ISS of the rHEALTH ONE, a flow cytometry based bioanalyzer, to assess future devices based on this design. In flow cytometry, there are predominantly three ways to focus the cells (or particles) into a single file stream: hydrodynamic focusing, microcapillary, and acoustic focusing. The rHEALTH ONE utilizes the commercial standard, sheath-based hydrodynamic focusing. rHEALTH itself represents both a company and a suite of medical tools NASA has funded through SBIR grants towards the development of a promising diagnosis instrument for exploration missions. rHEALTH ONE is the interim version of the technology, functional as a benchtop analyzer and test bed for the next generation of rHEALTH in development. Several modifications were made to the rHEALTH ONE analyzer for operation in microgravity. For function, fluid management was key. A sheath-based analyzer uses sheath fluid to flow the sample, cleaning fluid to prevent biological contamination, and a reservoir to collect the liquid waste. The rHEALTH ONE analyzer uses bottles dependent on gravity to separate the air and liquid pathways and keep the liquids contained. It uses only 1 psig of air pressure to directly push the liquid through the device, requiring little-to-no resistance at inlet and outlet. A microtubing assembly with self-sealing luer connectors was designed – featuring 0.014 mm thick durable medical balloons to hold water inside the supply bottles – to create safe containment and easy access for the crew while maintaining the analyzer’s function. For safety, copper tape was added to the interior of the plastic housing to reduce electromagnetic interference, fluid and electrical connections were secured against vibration and leaks, and gaps were further sealed to ensure containment of the optical block and lasers. Water as the sheath and cleaning fluid and TOX 0 samples were used to reduce the biohazard risk to the crew. In May 2022, European Space Agency astronaut Samantha Cristoforetti demonstrated the rHEALTH ONE aboard ISS for its sample loading, flow cytometry, and data collection capability in microgravity. The JSC Immunology Lab provided flow cytometry expertise, designed the sample test protocol, and manufactured and benchmarked the flight samples on a ‘gold-standard’ flow cytometer. The analyzer was primed with water, purged of air, and four calibration solutions of polystyrene microparticles were tested to characterize its performance. Tightly controlled procedures and excellent execution prevented air bubble interference during air/water separation, fluid transfer, sample mixing and loading. Data showed a slight increase in signal noise on 3 of 5 channels and an anomaly of fluorophores migrating in one sample (confirmed by JSC post-flight). Flight results correctly detected the change in sample, matched the analyzer’s ground performance, and were consistent with the gold-standard’s ground results. Although areas were noted for improvement, these outcomes signified complete mission success.

R. S. Miller↗

Analysis of Geologic CO 2 Migration Pathways in Farnsworth Field, NW Anadarko Basin

This study reports on analyses of natural, geologic CO 2 migration paths in Farnsworth Oil Field, northern Texas, where CO 2 was injected into the Pennsylvanian Morrow B reservoir as part of enhanced oil recovery and carbon sequestration efforts. We interpret 2D and 3D seismic reflection datasets of the study site, which is located on the western flank of the Anadarko basin, and compare our seismic interpretations with results from a tracer study. Petroleum system models are developed to understand the petroleum system and petroleum- and CO 2 -migration pathways. We find no evidence of seismically resolvable faults in Farnsworth Field, but interpret a karst structure, erosional structures, and incised valleys. These interpretations are compared with results of a Morrow B well-to-well tracer study that suggests that inter-well flow is up-dip or lateral. Southeastward fluid flow is inhibited by dip direction, thinning, and draping of the Morrow B reservoir over a deeper, eroded formation. Petroleum system models predict a deep basin-ward increase in temperature and maturation of the source rocks. In the northwestern Anadarko Basin, petroleum migration was generally up-dip with local exceptions; the Morrow B sandstone was likely charged by formations both below and overlying the reservoir rock. Based on this analysis, we conclude that CO 2 escape in Farnsworth Field via geologic pathways such as tectonic faults is unlikely. Abandoned or aged wellbores remain a risk for CO 2 escape from the reservoir formation and deserve further monitoring and research.

58 GEOSCIENCES↗

Remobilizing the Interfaces of Thermocapillary Driven Bubbles Retarded by the Adsorption of a Surfactant Impurity on the Bubble Surface

Thermocapillary migration is a method for moving bubbles in space in the absence of buoyancy. A temperature gradient is applied to the continuous phase in which a bubble is situated, and the applied gradient impressed on the bubble surface causes one pole of the drop to be cooler than the opposite pole. As the surface tension is a decreasing function of temperature, the cooler pole pulls at the warmer pole, creating a flow which propels the bubble in the direction of the warmer fluid. A major impediment to the practical use of thermocapillarity to direct the movement of bubbles in space is the fact that surfactant impurities which are unavoidably present in the continuous phase can significantly reduce the migration velocity. A surfactant impurity adsorbed onto the bubble interface is swept to the trailing end of the bubble. When bulk concentrations are low (which is the case with an impurity), diffusion of surfactant to the front end is slow relative to convection, and surfactant collects at the back end of the bubble. Collection at the back lowers the surface tension relative to the front end setting up a reverse tension gradient. For buoyancy driven bubble motions in the absence of a thermocapillarity, the tension gradient opposes the surface flow, and reduces the surface and terminal velocities (the interface becomes more solid-like). When thermocapillary forces are present, the reverse tension gradient set up by the surfactant accumulation reduces the temperature tension gradient, and decreases to near zero the thermocapillary velocity. The objective of our research is to develop a method for enhancing the thermocapillary migration of bubbles which have been retarded by the adsorption onto the bubble surface of a surfactant impurity, Our remobilization theory proposes to use surfactant molecules which kinetically rapidly exchange between the bulk and the surface and are at high bulk concentrations. Because the remobilizing surfactant is present at much higher concentrations than the impurity, it adsorbs to the bubble much faster than the impurity when the bubble is formed, and thereby prevents the impurity from adsorbing onto the surface. In addition the rapid kinetic exchange and high bulk concentration maintain a saturated surface with a uniform surface concentrations. This prevents retarding surface tension gradients and keeps the velocity high. In our first report last year, we detailed experimental results which verified the theory of remobilization in ground based experiments in which the steady velocity of rising bubbles was measured in a continuous phase consisting of a glycerol/water mixture containing a polyethylene glycol surfactant C12E6 (CH3(CH2)11(OCH2CH2)6OH). In our report this year, we detail our efforts to describe theoretically the remobilization observed. We construct a model in which a bubble rises steadily by buoyancy in a continuous (Newtonian) viscous fluid containing surfactant with a uniform far field bulk concentration. We account for the effects of inertia as well as viscosity in the flow in the continuous phase caused by the bubble motion (order one Reynolds number), and we assume that the bubble shape remains spherical (viscous and inertial forces are smaller than capillary forces, i e. small Weber and capillary numbers). The surfactant distribution is calculated by solving the mass transfer equations including convection and diffusion in the bulk, and finite kinetic exchange the bulk and the surface. Convective effects dominate diffusive mass transfer in the bulk of the liquid (high Peclet numbers) except in a thin boundary layer near the surface. A finite volume method is used to numerically solve the hydrodynamic and mass transfer equations on a staggered grid which accounts specifically for the thin boundary layer. We present the results of the nondimensional drag as a function of the bulk concentration of surfactant for different rates of kinetic exchange, from which we develop criteria for the concentration necessary to develop a prescribed degree of remobilization. The criteria compare favorably with the experimental results.

Palaparthi, Ravi↗

Evaluation of Electrical Resistivity Tomography to Monitor the Transport of Past Releases Beneath Tank Farms

Underground storage tanks at the Hanford Site, in southeastern Washington State, hold radioactive waste generated from four decades of plutonium production. The 149 single-shell tanks and the 28 double-shell tanks have all exceeded their initial design life of approximately 25 years. At least 67 tanks are assumed to have leaked in the past, resulting in radioactive releases into the vadose zone. Gamma ray logging within dry monitoring wells is currently the primary method for tracking the migration of leaked tank waste through the vadose zone. While this approach provides an accurate assessment of radioactive contamination, that information is only provided near (within ~1m) the borehole, leaving most of the vadose zone unmonitored, particularly the important region directly beneath the tank. This report describes a numerical study that investigates the feasibility and performance of time-lapse 3D electrical resistivity tomography (ERT) for long-term monitoring of a hypothetical tank waste location and migration through the vadose zone. ERT is a method of remotely imaging the bulk electrical conductivity (EC) of the subsurface, which is significantly impacted by the presence of conductive solid and liquid tank waste. The release of liquid tank wastes increases subsurface fluid conductivity and saturation over time, creating a target to use time-lapse ERT for long-term monitoring. Although the presence of metallic infrastructure can cause ERT interference, recent advancements in ERT data processing enable the deleterious effects of buried metallic infrastructure (e.g. pipes, wellbore casings, tanks) to be removed to better determine the liquid tank waste migration over time. Three hypothetical realistic scenarios were simulated in the ERT evaluation. The first two scenarios assume the same leak amount and rate (i.e., between 1/1/1951 and 12/31/1951 at the rate of 347 m 3 per year) but different leaky tanks. Scenario 1 assumes leaks under tank B-102, which is located on the edge of the B-tank farm and surrounded by a few metallic infrastructure including cased pipes/wells/tanks. Scenario 2 assumes leaks under tank B-108, which is located near the center of the B-tank farm and surrounded by larger amount of metallic infrastructure than B-102. Scenario 3 assumes the same metallic infrastructure as B-102, with a more recent contaminant leak that was simulated to have occurred between 1/1/2018 and 12/31/2023 at a rate of 1.89 m 3 per year. The leak time in Scenarios 1 and 2 corresponds to a historical overfill event in 1951 and Scenario 3 corresponds to a recent found tank leak in 2019. In each scenario, a “true” bulk EC model vs. time reflecting contaminant migration was generated. ERT data was simulated from these “true” bulk EC models and a time-lapse ERT inversion produced “imaged” bulk EC vs. time. Three electrode configurations in two, four and eight boreholes surrounding the leak tank were used in the ERT simulations in each scenario. These borehole configurations were considered logistically feasible and cost-effective for monitoring. The hypothetical ERT boreholes are assumed to have non-metallic casing. By comparing the “imaged” bulk EC with the “true” bulk EC, it was demonstrated that the three configurations of wells used (two, four, and eight wells) were able to successfully monitor the migration of tank leaks through the vadose zone, with bulk EC resolution increasing with the number of down borehole ERT arrays for all scenarios. Therefore, the use of eight boreholes to perform ERT monitoring beneath the tanks provided the best spatiotemporal information.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Remobilizing the Interface of Thermocapillary Driven Bubbles Retarded By the Adsorption of a Surfactant Impurity on the Bubble Surface

Thermocapillary migration is a method for moving bubbles in space in the absence of buoyancy. A temperature gradient is the continuous phase in which a bubble is situated, and the applied gradient impressed on the bubble surface causes one pole of the drop to be cooler than the opposite pole. As the surface tension is a decreasing function of temperature, the cooler pole pulls at the warmer pole, creating a flow that propels the bubble in the direction of the warmer fluid. A major impediment to the practical use of thermocapillary to direct the movement of bubbles in space is the fact that surfactant impurities, which are unavoidably present in the continuous phase, can significantly reduce the migration velocity. A surfactant impurity adsorbed onto the bubble interface is swept to the trailing end of the bubble. When bulk concentrations are low (which is the case with an impurity), diffusion of surfactant to the front end is slow relative to convection, and surfactant collects at the back end of the bubble. Collection at the back lowers the surface tension relative to the front end setting up a reverse tension gradient. (This can also be the case if kinetic desorption of surfactant at the back end of the bubble is much slower than convection.) For buoyancy driven bubble motions in the absence of a thermocapillarity, the tension gradient opposes the surface flow, and reduces the surface and terminal velocities (the interface becomes more solid-like and bubbles translate as solid particles). When thermocapillary forces are present, the reverse tension gradient set up by the surfactant accumulation reduces the temperature-induced tension gradient, and can decrease to near zero the bubble's thermocapillary velocity. The objective of our research is to develop a method for enhancing the thermocapillary migration of bubbles which have be retarded by the adsorption onto the bubble surface of a surfactant impurity. Our remobilization theory proposes to use surfactant molecules which kinetically rapidly exchange between the bulk and the surface and are at high bulk concentrations. Because the remobilizing surfactant is present at much higher concentrations than the impurity, it adsorbs to the bubble surface much faster than the impurity when the bubble is formed, and thereby prevents the impurity from adsorbing onto the surface. In addition, the rapid kinetic exchange and high bulk concentration maintain a saturated surface with uniform surface concentrations. This prevents retarding surface tension gradients and keeps the thermocapillary velocity high. In our reports over the first 2 years, we presented numerical simulations of the bubble motion and surfactant transport which verified theoretically the concept of remobilization, and the development of an apparatus to track and measure the velocity of rising bubbles in a glycerol/water surfactant solution. This year, we detail experimental observations of remobilization. Two polyethylene oxide surfactants were studied, C12E6 (CH3(CH2)11(OCH2)6OH) and C10E8 (CH3(CH2)4(OCH2CH2)8OH). Measurements of the kinetic exchange for these surfactants show that the one with the longer hydrophobe chain C12E6 has a lower rate of kinetic exchange. In addition, this surfactant is much less soluble in the glycerol/water mixture because of the shorter ethoxylate chain. As a result, we found that C12E6 had only a very limited ability to remobilize rising bubbles because of the limited kinetic exchange and reduced solubility. However, C10E8, with its higher solubility and more rapid exchange was found to dramatically remobilize rising bubbles. We also compared our theoretical calculations to the experimental measurements of velocity for both the non-remobilizing and remobilizing surfactants and found excellent agreement. We further observed that for C10E8 at high concentrations, which exceeded the critical micelle concentrations, additional remobilization was measured. In this case the rapid exchange of monomer between micelle and surfactant provides an additional mechanism for maintaining a uniform surface concentrations.

Palaparthi, Ravi↗

Fundamental Studies of Fluid Mechanics: Stability in Porous Media

We summarize our research results in three main areas: coating flows; electrohydrodynamics of drops; and wetting and spreading of drops. Experimental, computational and analytical methods are used to address a variety of issues. Coating flow studies include the effect of roughness, surfactants, and adsorbed particles on the dynamics of dip-coating. Electrohydrodynamic studies include drop deformation in uniform electric fields, shape distortion due to charge convection in sedimenting drops, and driving chaotic advection by either an electric field inclined to the direction of drop motion or time-periodic changes in the direction of the electric field. Heat and mass transport from chaotically mixed droplets exhibit unexpected and remarkable increases in the rates of transport. Finally, we develop an analytical solution to the problem of a static droplet, and use numerical techniques to predict its migration due to surface tension gradients.

02 PETROLEUM↗