Results for “channel estimation” — Page 7

DOE OSTI2024

Resilient Consensus State Observer for Nonlinear Systems and Against Attacks

Herein, we investigate the design of a resilient consensus observer that yields correct and continuous state estimates of a nonlinear control system whose state observations are subject to attacks. The observations are assumed to be obtained from multiple sensors and/or communication channels employed to provide redundancy. It is shown that, in an attack-free setting, individual globally convergent observers can be designed for each sensor/channel using the corresponding state-dependent Jacobian system. To counter attacks, a resilient consensus observer is then designed by merging the observers estimates to produce a smooth consensus estimate. Resilient convergence of the consensus estimate to the true state is ensured under the condition that the sparsity of the attacks being less than a half of the number of redundant channels. As an illustrative example, the proposed scheme is successfully applied to the chaotic circuit synchronization problem in which their synchronization has been attacked.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2026

Dynamics of the water-plasma interface in various discharge modes of atmospheric-pressure plasmas

This paper presents an experimental study of the dynamics of the water-plasma interface during the interaction high-voltage (HV) atmospheric pressure discharge with the water surface. Under the influence of the voltage applied to the pin type electrode, discharges are formed at the air layer between the sharp-tip HV electrode and deionized water. Three discharge regimes were identified from synchronized electrical waveforms and imaging: (i) a weak linear regime at 3–10.6 kV, (ii) a branching streamer regime at the maximum applied voltage of 12.6 kV, and (iii) a continuous arc-like regime occurring during the transition to lower voltage and higher current (e.g., U ≈ 3.2 kV, I ≈ 4.07 mA). High-speed shadowgraph images showed a symmetric interfacial cavity whose depth increased nonlinearly with voltage from h 0 ≈ 0.1 mm at 3 kV to h 0 ≈ 2.7 mm at 10.6 kV, reaching a maximum depth of h 0 ≈ 5.9 mm at 12.6 kV, while the cavity disappeared in the continuous-channel regime and was replaced by outward-propagating wave-like motion. The deformation of the water surface during the discharge-water interaction is governed by the balance between electric field forces, surface tension, gravitational forces, and electrohydrodynamic forces, whose relative contributions vary with the applied voltage. This effect was explained by quantitatively estimating the magnitudes of the acting forces and establishing the force balance. Optical emission spectroscopy in the continuous-channel regime indicated air plasma signatures and yielded a gas temperature of approximately 400 K, while prolonged operation (≈ 20 min) increased the water temperature to ~ 70 °C, reduced the water-layer thickness from 6 mm to 3 mm, and decreased pH from 7 to 4. These regime-resolved, quantitative results clarify how the dominant interfacial forcing shifts with discharge mode and provide a mechanistic basis for controlling plasma-water interactions in atmospheric-pressure applications.

atmospheric pressure plasma↗

DOE OSTI2020

Performance simulation of the soft gamma-ray concentrator

The soft gamma-ray concentrator is a telescope mission concept utilizing a suitable arrangement of bent multilayer structures of alternating low- and high-density materials. This lens is able to channel gamma-ray photons via total external reflection and concentrate the incident radiation to a point. The channeling technique offers the potential for concentrating gamma rays with focal lengths <10 m and energies >100 keV, beyond the reach of current grazing-incidence hard x-ray mirrors. For the performance estimation of such an instrument, we have developed a flexible set of computer modeling tools to compute the optical properties of multilayer structures, predict the channeling efficiency for a given multilayer configuration, and aid in the optimization of potential gamma-ray concentrator-based telescope designs. This modeling includes the multilayer optical properties calculated by the IMD software, the ray tracing using an IDL code, and the focal plane detector simulation by MEGAlib. We illustrate the potential of this approach by presenting simulated astronomical observations from a balloon-borne platform. The final result, including simulated effective area, instrument sensitivity, and polarization performance, shows that the gamma-ray concentrator will provide greatly increased sensitivity for next-generation soft gamma-ray missions with modest cost and complexity.

47 OTHER INSTRUMENTATION↗

DOE OSTI2021

DANTE as a primary temperature diagnostic for the NIF iron opacity campaign

The Opacity Platform on the National Ignition Facility (NIF) has been developed to measure iron opacities at varying densities and temperatures relevant to the solar interior and to verify recent experimental results obtained at the Sandia Z-machine, that diverge from theory. The first set of NIF experiments collected iron opacity data at ∼150 eV to 160 eV and an electron density of ∼7 × 1021 cm−3, with a goal to study temperatures up to ∼210 eV, with electron densities of up to ∼3 × 1022 cm−3. Among several techniques used to infer the temperature of the heated Fe sample, the absolutely calibrated DANTE-2 filtered diode array routinely provides measurements of the hohlraum conditions near the sample. However, the DANTE-2 temperatures are consistently low compared to pre-shot LASNEX simulations for a range of laser drive energies. We have re-evaluated the estimated uncertainty in the reported DANTE-2 temperatures and also the error generated by varying channel participation in the data analysis. An uncertainty of ±5% or better can be achieved with appropriate spectral coverage, channel participation, and metrology of the viewing slot.

Opachich, Y. P. (ORCID:0000000314323989)↗

DOE OSTI2020

Effects of matter density profiles on neutrino oscillations for T2HK and T2HKK

This paper explores the effects of changes in matter density profiles on neutrino oscillation probabilities, and whether these could potentially be seen by the future Hyper-Kamiokande long-baseline oscillation experiment (T2HK). The analysis is extended to include the possibility of having an additional detector in Korea (T2HKK). In both cases, we find that these effects will be immeasurable, as the magnitudes of the changes in the oscillation probabilities induced in all density profile scenarios considered here remain smaller than the estimated experimental sensitivity to the oscillation probabilities of each experiment, for both appearance and disappearance channels. Therefore, we conclude that using a constant density profile is sufficient for both the T2HK and T2HKK experiments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2023

Improving Signal-to-Noise Ratio (SNR) for Readout Signals Using Adaptive Filters on Reconfigurable Controls Hardware

This study investigates the optimization of Signal-to-Noise Ratio (SNR) in superconducting quantum computing readout signals through adaptive filtering. Quantum computing technology has the potential to revolutionize various fields by delivering exponential speedup in solving certain computational problems. However, the technology's practical implementation is hindered by the difficulty of extracting clean, reliable signals during the readout phase, with various sources of noise presenting a significant barrier to clean signals. This noise, often present in readout profiles due to imperfect isolation, degrades the system's overall SNR, thus impeding the ability to extract the quantum state accurately. The research leverages the power of adaptive filtering to improve the SNR of quantum computing readout signals. Specifically, an adaptive filter is implemented in a PYNQ overlay on an FPGA, and eventually will be connected to a quantum computing system. The system models the noise with a Least Mean Squares (LMS) adaptive filter, and then subtracts the estimated noise from the received signal to improve the SNR. A Direct Memory Access (DMA) channel is used to handle the signal processing, delivering efficient, high-speed data transfer between the PYNQ system and the hardware. The study explores the benefits of this adaptive filtering technique, potentially providing a significant contribution to practical and fast quantum computing.

Johnson, Hans↗

DOE OSTI2023

Improving Signal-to-Noise Ratio (SNR) for Readout Signals Using Adaptive Filters on Reconfigurable Controls Hardware

This study investigates the optimization of Signal-to-Noise Ratio (SNR) in superconducting quantum computing readout signals through adaptive filtering. Quantum computing technology has the potential to revolutionize various fields by delivering exponential speedup in solving certain computational problems. However, the technology's practical implementation is hindered by the difficulty of extracting clean, reliable signals during the readout phase, with various sources of noise presenting a significant barrier to clean signals. This noise, often present in readout profiles due to imperfect isolation, degrades the system's overall SNR, thus impeding the ability to extract the quantum state accurately. The research leverages the power of adaptive filtering to improve the SNR of quantum computing readout signals. Specifically, an adaptive filter is implemented in a PYNQ overlay on an FPGA, and eventually will be connected to a quantum computing system. The system models the n oise with a Least Mean Squares (LMS) adaptive filter, and then subtracts the estimated noise from the received signal to improve the SNR. A Direct Memory Access (DMA) channel is used to handle the signal processing, delivering efficient, high-speed data transfer between the PYNQ system and the hardware. The study explores the benefits of this adaptive filtering technique, potentially providing a significant contribution to practical and fast quantum computing.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2021

Analysis of Approximations in Modeling of BWR Bundle Void Distributions

In boiling water reactors, complex heterogeneous bundle designs, control blades adjacent to the corner of bundles, and the presence of boiling can lead to complex internal void distributions. A few approximations exist to model these void distributions. They could be modeled using a 1D axial solver in which each axial node is assumed to be at an average void, or each pin cell could be modeled with its own void concentration. In the latter case, the void could be discretized in pin-centered or coolant-centered channels. The goal of this project was to quantify the effect of using the different approximations for modeling internal void distributions on neutronics calculations. Using 3D void distributions calculated with CTF, Monte Carlo Neutral Particle (MCNP) transport code models were created for GE-9 and GE-14 lattices. For each model, the internal void distribution from CTF at a given axial node was selected, and a lattice calculation was carried out with MCNP. Comparisons between models using a lattice-averaged void, or using a void distribution in coolant-centered channels, showed large differences in reactivity which in some cases were well above 1,000 pcm, and it also showed differences in normalized fission rates greater than 20%. It was also found that using a lattice average void can lead to a significant difference in the estimation of the worth of a control blade. The differences found when comparing results from models using pin-centered and coolant-centered channels were up to 200 pcm in reactivity and up to 1.4% in the normalized fission rates. In addition to these two sets of comparisons, MCNP models were set up so that each subchannel had a saturated liquid component around the fuel pins and a saturated vapor component in the center to approximate annular flow. In comparison to the models using coolant-centered subchannels, up to 1–3% differences in normalized fission rates could be found.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2023

Optimizing cloud motion estimation on the edge with phase correlation and optical flow

Abstract. Phase correlation (PC) is a well-known method for estimating cloud motion vectors (CMVs) from infrared and visible spectrum images. Commonly, phase shift is computed in the small blocks of the images using the fast Fourier transform. In this study, we investigate the performance and the stability of the blockwise PC method by changing the block size, the frame interval, and combinations of red, green, and blue (RGB) channels from the total sky imager (TSI) at the United States Atmospheric Radiation Measurement user facility's Southern Great Plains site. We find that shorter frame intervals, followed by larger block sizes, are responsible for stable estimates of the CMV, as suggested by the higher autocorrelations. The choice of RGB channels has a limited effect on the quality of CMVs, and the red and the grayscale images are marginally more reliable than the other combinations during rapidly evolving low-level clouds. The stability of CMVs was tested at different image resolutions with an implementation of the optimized algorithm on the Sage cyberinfrastructure test bed. We find that doubling the frame rate outperforms quadrupling the image resolution in achieving CMV stability. The correlations of CMVs with the wind data are significant in the range of 0.38–0.59 with a 95 % confidence interval, despite the uncertainties and limitations of both datasets. A comparison of the PC method with constructed data and the optical flow method suggests that the post-processing of the vector field has a significant effect on the quality of the CMV. The raindrop-contaminated images can be identified by the rotation of the TSI mirror in the motion field. The results of this study are critical to optimizing algorithms for edge-computing sensor systems.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Nucleon clustering at kinetic freezeout of heavy-ion collisions via path-integral Monte Carlo

Clustering of the four-nucleon system at kinetic freezeout conditions is studied using path-integral Monte Carlo techniques. This method seeks to improve upon previous calculations which relied on approximate semiclassical methods or few-body quantum mechanics. Estimates are given for the decay probabilities of the 4N system into various light nuclei decay channels and the strength of spatial correlations is characterized. Additionally, a simple model is presented to describe the impact of this clustering on nucleon multiplicity distributions. Additionally, the effects of a possible modification of the inter-nucleon interaction due to the close critical line (and hypothetical QCD critical point) on the clustering are also studied.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2023

Air‐Sea CO 2 Fluxes Localized by Topography in a Southern Ocean Channel

Abstract Air‐sea exchange of carbon dioxide (CO 2 ) in the Southern Ocean plays an important role in the global carbon budget. Previous studies have suggested that flow around topographic features of the Southern Ocean enhances the upward supply of carbon from the deep to the surface, influencing air‐sea CO 2 exchange. Here, we investigate the role of seafloor topography on the transport of carbon and associated air‐sea CO 2 flux in an idealized channel model. We find elevated CO 2 outgassing upstream of a seafloor ridge, driven by anomalous advection of dissolved inorganic carbon. Argo‐like Lagrangian particles in our channel model sample heterogeneously in the vicinity of the seafloor ridge, which could impact float‐based estimates of CO 2 flux.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2020

Stastically determined experimental uncertainties in photon Doppler velocimetry (PDV) measurements

Shock-ramp experiments place high demands on the precision, accuracy, and temporal resolution of transformed velocimetry signals. To minimize uncertainty after total propagation of a shockramp experiment, ideally multiple PDV traces are used. However, each trace has some uncertainty that is a function of the PDV system and analysis of the velocimetry data. In order to determine a reasonable experimental estimate of minimum uncertainty, we calculate observed uncertainties in as-fielded PDV systems using multiple channels to sample the velocity on each shot. We also calculate the uncertainty’s sensitivity to choice of analysis methods.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

An Algorithm to Retrieve Total Precipitable Water Vapor in the Atmosphere from FengYun 3D Medium Resolution Spectral Imager 2 (FY-3D MERSI-2) Data

The atmosphere has substantial effects on optical remote sensing imagery of the Earth’s surface from space. These effects come through the functioning of atmospheric particles on the radiometric transfer from the Earth’s surface through the atmosphere to the sensor in space. Precipitable water vapor (PWV), CO 2 , ozone, and aerosol in the atmosphere are very important among the particles through their functioning. This study presented an algorithm to retrieve total PWV from the Chinese second-generation polar-orbiting meteorological satellite FengYun 3D Medium Resolution Spectral Imager 2 (FY-3D MERSI-2) data, which have three near-infrared (NIR) water vapor absorbing channels, i.e., channel 16, 17, and 18. The algorithm was improved from the radiance ratio technique initially developed for Moderate-Resolution Imaging Spectroradiometer (MODIS) data. MODTRAN 5 was used to simulate the process of radiant transfer from the ground surfaces to the sensor at various atmospheric conditions for estimation of the coefficients of ratio technique, which was achieved through statistical regression analysis between the simulated radiance and transmittance values for FY-3D MERSI-2 NIR channels. The algorithm was then constructed as a linear combination of the three-water vapor absorbing channels of FY-3D MERSI-2. Measurements from two ground-based reference datasets were used to validate the algorithm: the sun photometer measurements of Aerosol Robotic Network (AERONET) and the microwave radiometer measurements of Energy’s Atmospheric Radiation Measurement Program (ARMP). The validation results showed that the algorithm performs very well when compared with the ground-based reference datasets. The estimated PWV values come with root mean square error (RMSE) of 0.28 g/cm 2 for the ARMP and 0.26 g/cm 2 for the AERONET datasets, with bias of 0.072 g/cm 2 and 0.096 g/cm 2 for the two reference datasets, respectively. The accuracy of the proposed algorithm revealed a better consistency with ground-based reference datasets. Thus, the proposed algorithm could be used as an alternative to retrieve PWV from FY-3D MERSI-2 data for various remote sensing applications such as agricultural monitoring, climate change, hydrologic cycle, and so on at various regional and global scales.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Supercritical CO2 Recuperators, Presented at The ASME Turbo Expo 2017, June 29, 2017, by Dr. John Kelly, President, Altex Technologies Corporation

Closed Brayton super-critical CO2 power cycles are well suited to waste heat bottoming cycles, due to their increased efficiency and compactness, relative to Rankine steam bottoming cycles. Since waste heat applications would be retrofits, the power system compactness is important. To achieve high efficiency, these power cycles require high pressure recuperative type heat exchangers, of substantial heat duty. Current Printed Circuit Heat Exchangers (PCHE), originally developed for high pressure gas and oil applications, can be used as recuperators, but costs are higher than desired. Altex is developing a purpose-built high pressure and effectiveness recuperator to provide the reliability, compactness, performance and pressure capability of current recuperators, but at a reduced cost. The High Effectiveness Low Cost (HELC) recuperative heat exchanger design yields volume and weight metrics of .0024 m3/UA and 10.2 kg/UA, which are 4% and 77.3% below recuperator target metrics, respectively. A 50 kW test article was designed and fabricated. Performance tests on water and oil showed that the HELC design model could predict heat transfer, to within 10% of the measured value. This model was then used to project HELC performance, when operating on supercritical CO2. Besides performance tests, the test article was hydrostatically tested, for integrity at up to 4,000 psi pressure. At these conditions, some distortion of channels was encountered. To mitigate distortion, the inserts were redesigned and these results were used to project the cost of 500 kW and 5,083 kW HELC units. The cost metric for the 5,083 kW unit was determined to be $1,349/UA, which is 10% lower than the recuperator desired cost metric desired target of $1,500/UA. In addition, the 5,083 kW unit HELC cost of $61.09/kW is 33.6% lower than the $92/kW estimated cost for a PCHE. Hydrostatic pressure tests, at up to 4,000psi, showed that the unit did not leak. However, channels were distorted at this pressure differential. Design updates to minimize stress concentrations and distortion were prepared and analyzed, to show that distortion could be controlled, but to date tests have not been run to prove the design. Project results show the potential of the HELC approach, but more work is required to confirm this potential, at the larger scales of interest.

14 SOLAR ENERGY↗

DOE OSTI2020

High Effectiveness, Compact, High Pressure and Low-Cost Recuperator, Presented at The Fifth International Symposium – Super-Critical CO2 Power Cycles, San Antonio, Texas, March 28-31, 2016, by Dr. John Kelly, President, Altex Technologies Corporation

Closed Brayton super-critical CO2 power cycles are well suited to waste heat bottoming cycles, due to their increased efficiency and compactness, relative to Rankine steam bottoming cycles. Since waste heat applications would be retrofits, the power system compactness is important. To achieve high efficiency, these power cycles require high pressure recuperative type heat exchangers, of substantial heat duty. Current Printed Circuit Heat Exchangers (PCHE), originally developed for high pressure gas and oil applications, can be used as recuperators, but costs are higher than desired. Altex is developing a purpose-built high pressure and effectiveness recuperator to provide the reliability, compactness, performance and pressure capability of current recuperators, but at a reduced cost. The High Effectiveness Low Cost (HELC) recuperative heat exchanger design yields volume and weight metrics of .0024 m3/UA and 10.2 kg/UA, which are 4% and 77.3% below recuperator target metrics, respectively. A 50 kW test article was designed and fabricated. Performance tests on water and oil showed that the HELC design model could predict heat transfer, to within 10% of the measured value. This model was then used to project HELC performance, when operating on supercritical CO2. Besides performance tests, the test article was hydrostatically tested, for integrity at up to 4,000 psi pressure. At these conditions, some distortion of channels was encountered. To mitigate distortion, the inserts were redesigned and these results were used to project the cost of 500 kW and 5,083 kW HELC units. The cost metric for the 5,083 kW unit was determined to be $1,349/UA, which is 10% lower than the recuperator desired cost metric desired target of $1,500/UA. In addition, the 5,083 kW unit HELC cost of $61.09/kW is 33.6% lower than the $92/kW estimated cost for a PCHE. Hydrostatic pressure tests, at up to 4,000psi, showed that the unit did not leak. However, channels were distorted at this pressure differential. Design updates to minimize stress concentrations and distortion were prepared and analyzed, to show that distortion could be controlled, but to date tests have not been run to prove the design. Project results show the potential of the HELC approach, but more work is required to confirm this potential, at the larger scales of interest.

14 SOLAR ENERGY↗

DOE OSTI2020

High Effectiveness, Compact, High Pressure and Low Cost Recuperator for Super-Critical CO2 Power Cycles, Paper Presented at The Fifth International Symposium – Super-Critical CO2 Power Cycles, San Antonio, Texas, March 28-31, 2016, by Dr. John Kelly, President, Altex Technologies Corporation

Closed Brayton super-critical CO2 power cycles are well suited to waste heat bottoming cycles, due to their increased efficiency and compactness, relative to Rankine steam bottoming cycles. Since waste heat applications would be retrofits, the power system compactness is important. To achieve high efficiency, these power cycles require high pressure recuperative type heat exchangers, of substantial heat duty. Current Printed Circuit Heat Exchangers (PCHE), originally developed for high pressure gas and oil applications, can be used as recuperators, but costs are higher than desired. Altex is developing a purpose-built high pressure and effectiveness recuperator to provide the reliability, compactness, performance and pressure capability of current recuperators, but at a reduced cost. The High Effectiveness Low Cost (HELC) recuperative heat exchanger design yields volume and weight metrics of .0024 m3/UA and 10.2 kg/UA, which are 4% and 77.3% below recuperator target metrics, respectively. A 50 kW test article was designed and fabricated. Performance tests on water and oil showed that the HELC design model could predict heat transfer, to within 10% of the measured value. This model was then used to project HELC performance, when operating on supercritical CO2. Besides performance tests, the test article was hydrostatically tested, for integrity at up to 4,000 psi pressure. At these conditions, some distortion of channels was encountered. To mitigate distortion, the inserts were redesigned and these results were used to project the cost of 500 kW and 5,083 kW HELC units. The cost metric for the 5,083 kW unit was determined to be $1,349/UA, which is 10% lower than the recuperator desired cost metric desired target of $1,500/UA. In addition, the 5,083 kW unit HELC cost of $61.09/kW is 33.6% lower than the $92/kW estimated cost for a PCHE. Hydrostatic pressure tests, at up to 4,000psi, showed that the unit did not leak. However, channels were distorted at this pressure differential. Design updates to minimize stress concentrations and distortion were prepared and analyzed, to show that distortion could be controlled, but to date tests have not been run to prove the design. Project results show the potential of the HELC approach, but more work is required to confirm this potential, at the larger scales of interest.

14 SOLAR ENERGY↗

DOE OSTI2020

Active- and transfer-learning applied to microscale-macroscale coupling to simulate viscoelastic flows

Active- and transfer-learning are applied to microscale dynamics of polymer flows for the multiscale discovery of effective constitutive approximations required in viscoelastic flow simulation. The result is macroscopic rheology directly connected to a microstructural model. Micro and macroscale simulations are adaptively coupled by means of Gaussian process regression (GPR) to run the expensive microscale computations only as necessary. This multiscale method is demonstrated with flows of a polymer solution as a model system. At the microscale level dissipative particle dynamics (DPD) is employed to model the fluid as a suspension of bead-spring micro-structures subjected to steady shear flow. The results yield the non-Newtonian viscosity and the first normal stress difference at strain rates as training data used in a GPR model. DPD parameters are calibrated with respect to experimental data for a real polymer solution. Compliance with these data requires adjustment of the DPD model's cutoff radius, which then becomes a function of the second invariant of the strain rate tensor. The FENE-P model is chosen for the macroscale description using the spectral element method (SEM) to simulate channel flow and flow past a circular cylinder. The DPD results at the lowest possible shear strain rate yield an estimate of the zero-shear rate viscosity, which allows the initiation of the macroscale flow by SEM as a Newtonian fluid. The resulting strain-rate field is surveyed to determine additional shear strain rate sampling points for the DPD system. This new information allows an initial fitting of parameters of the constitutive equation followed by new SEM simulations at the macroscale. Additionally, guided by active-learning GPR to select new sampling points, this process continues until convergence is achieved. The effectiveness of this new simulation paradigm for viscoelastic flows is tested with different macroscale operating conditions. The effective closure learned in the channel simulation is then transferred directly to the flow past a circular cylinder at low Reynolds number, where the results show that only two additional DPD simulations are required to achieve a satisfactory constitutive model. With an increase of the Reynolds number, the active-learning scheme automatically detects the inaccuracy of the learned constitutive model, and initiates additional DPD simulations for the extra data needed to once again close the microscale-macroscale coupled system. This new paradigm of active- and transfer-learning for multiscale modeling is readily applicable to other microscale-macroscale coupled simulations of complex fluids and other materials. Furthermore, the coupling between microscale and macroscale solvers can be seamlessly implemented with our open source multiscale universal interface (MUI) library.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

Fast Advective Water Flow through Nanochannels in Clay Interlayers: Implications for Moisture Transport in Soils and Unconventional Oil/Gas Production

Water flow in nanometer or sub-nanometer hydrophilic channels bears special importance in diverse fields of science and engineering. However, the nature of such water flow remains elusive. In this work, we report our molecular-modeling results on water flow in a sub-nanometer clay interlayer between two montmorillonite layers. We show that a fast advective flow can be induced by evaporation at one end of the interlayer channel, that is, a large suction pressure created by evaporation (~818 MPa) is able to drive the fast water flow through the channel (~0.88 m/s for a 46 Å-long channel). Scaled up for the pressure gradient to a 2 μm particle, the velocity of water is estimated to be about 95 μm/s, indicating that water can quickly flow through a μm-sized clay particle within seconds. The prediction seems to be confirmed by our thermogravimetric analysis of bentonite hydration and dehydration processes, which indicates that water transport at the early stage of the dehydration is a fast advective process, followed by a slow diffusion process. The possible occurrence of a fast advective water flow in clay interlayers prompts us to reassess water transport in a broad set of natural and engineered systems such as clay swelling/shrinking, moisture transport in soils, water uptake by plants, water imbibition/release in unconventional hydrocarbon reservoirs, and cap rock integrity of supercritical CO 2 storage.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗