Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “diffusion length”

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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Experimental study of plate materials for evaporative air coolers

Following paper presents an experimental study on eight potential plate materials for evaporative air coolers. There are three groups of materials: currently used plate materials (kraft papers), new natural origin plate materials (cotton fibers), and new synthetic plate materials (synthetic fibers). Proposed materials were examined in terms of different properties such as upward and downward distribution height, water capacity, durability and diffusivity. The results were compared with data presented by other authors. It was established that, synthetic materials are characterized by superior water distribution qualities in comparison to the materials of natural origin. Basalt paper-like synthetic fiber allowed to rise water vertically to 40 cm during 120 min tests. The highest water rise obtained in other studies was 18 cm, hence basalt paper-like fiber is visibly superior in comparison with the currently known materials. In terms of downward water distribution, two new synthetic surfactant non-wovens show better performance than the existing materials by quickly and evenly distributing water over the entire length of the sample. The durability of synthetic fibers and non-wovens is higher comparing with the materials of natural origin. The results show that synthetic materials are better suited to the role of plate materials in evaporative air coolers.

Experiment↗
DOE OSTI2021

Light transport with weak angular dependence in fog

Random scattering and absorption of light by tiny particles in aerosols, like fog, reduce situational awareness and cause unacceptable down-time for critical systems or operations. Computationally efficient light transport models are desired for computational imaging to improve remote sensing capabilities in degraded optical environments. To this end, we have developed a model based on a weak angular dependence approximation to the Boltzmann or radiative transfer equation that appears to be applicable in both the moderate and highly scattering regimes, thereby covering the applicability domain of both the small angle and diffusion approximations. An analytic solution was derived and validated using experimental data acquired at the Sandia National Laboratory Fog Chamber facility. The evolution of the fog particle density and size distribution were measured and used to determine macroscopic absorption and scattering properties using Mie theory. A three-band (0.532, 1.55, and 9.68 µ m) transmissometer with lock-in amplifiers enabled changes in fog density of over an order of magnitude to be measured due to the increased transmission at higher wavelengths, covering both the moderate and highly scattering regimes. The meteorological optical range parameter is shown to be about 0.6 times the transport mean free path length, suggesting an improved physical interpretation of this parameter.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Two-fluid and kinetic transport physics of Kelvin–Helmholtz instabilities in nonuniform low-beta plasmas

Hall-magnetohydrodynamic (Hall-MHD) theory, two-fluid simulations, and kinetic simulations are used in this work to investigate the cross-field transport properties of Kelvin–Helmholtz instabilities in nonuniform low-beta collisionless plasmas. Hall-MHD analysis shows how the linear properties of the instability are modified by density gradients and magnetization. High-order accurate two-fluid and kinetic simulations, with complete dynamics of finite-mass electrons and ions, are applied to a suite of parameter cases to systematically assess the effects of diamagnetic drift, magnetization, charge separation, and finite Larmor motion. Initialization of exact two-species kinetic equilibria facilitates the study of isolated physical effects and enables detailed cross-comparisons between two-fluid and kinetic simulations, including for cases where ion gyroradii are comparable to gradient scale lengths. For nonuniform plasmas with significant space charge, the results of two-fluid and kinetic simulations are found to disagree with Hall-MHD predictions. Kelvin–Helmholtz instability growth rates, per unit shear, are shown to be smaller when ion diamagnetic drift and $E \times B$ drift are parallel and larger when the two drifts are antiparallel. The effect is attributed to polarization drift in the shear layer, which leads to redistribution of charge, alters the electric field that drives plasma advection, and consequently modifies growth rates. Instability-induced mass transport for different parameters is characterized in terms of the flux across the shear layer and a simplified diffusion model. Distribution functions from kinetic simulations are shown to deviate substantially from Maxwellian reconstructions, indicating the importance of kinetic physics during the nonlinear phase of the instability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Multi-modal imaging of a single mouse brain over five orders of magnitude of resolution

Mammalian neurons operate at length scales spanning six orders of magnitude; they project millimeters to centimeters across brain regions, are composed of micrometer-scale-diameter myelinated axons, and ultimately form nanometer scale synapses. Capturing these anatomical features across that breadth of scale has required imaging samples with multiple independent imaging modalities. Translating between the different modalities, however, requires imaging the same brain with each. Here, we imaged the same postmortem mouse brain over five orders of spatial resolution using MRI, whole brain micrometer-scale synchrotron x-ray tomography ($\mu$CT), and large volume automated serial electron microscopy. Using this pipeline, we can track individual myelinated axons previously relegated to axon bundles in diffusion tensor MRI or arbitrarily trace neurons and their processes brain-wide and identify individual synapses on them. This pipeline provides both an unprecedented look across a single brain's multi-scaled organization as well as a vehicle for studying the brain's multi-scale pathologies.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Reshaping of Truncated Pd Nanocubes: Energetic and Kinetic Analysis Integrating Transmission Electron Microscopy with Atomistic-Level and Coarse-Grained Modeling

Stability against reshaping of metallic fcc nanocrystals synthesized with tailored far-from-equilibrium shapes is key to maintaining optimal properties for applications such as catalysis. Yet Arrhenius analysis of experimental reshaping kinetics, and appropriate theory and simulation, is lacking. As such, we use TEM to monitor the reshaping of Pd nanocubes of ~25 nm side length between 410 °C (over ~4.5 h) and 440 °C (over ~0.25 h), extracting a high effective energy barrier of E eff ≈ 4.6 eV. We also provide an analytic determination of the energy variation along the optimal pathway for reshaping that involves transfer of atoms across the nanocube surface from edges or corners to form new layers on side {100} facets. The effective barrier from this analysis is shown to increase strongly with the degree of truncation of edges and corners in the synthesized nanocube. Theory matches experiment for the appropriate degree of truncation. In addition, we perform simulations of a stochastic atomistic-level model incorporating a realistic description of diffusive hopping for undercoordinated surface atoms, thereby providing a visualization of the initial reshaping process.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Temperature dependence of resistivity of carbon micro/nanostructures: Microscale spatial distribution with mixed metallic and semiconductive behaviors

The temperature coefficient of resistivity (θ T ) of carbon-based materials is a critical property that directly determines their electrical response upon thermal impulses. It could have metal- (positive) or semiconductor-like (negative) behavior, depending on the combined temperature dependence of electron density and electron scattering. Its distribution in space is very difficult to measure and is rarely studied. Here, for the first time, we report that carbon-based micro/nanoscale structures have a strong non-uniform spatial distribution of θ T . This distribution is probed by measuring the transient electro-thermal response of the material under extremely localized step laser heating and scanning, which magnifies the local θ T effect in the measured transient voltage evolution. For carbon microfibers (CMFs), after electrical current annealing, θ T varies from negative to positive from the sample end to the center with a magnitude change of >130% over <1 mm. This θ T sign change is confirmed by directly testing smaller segments from different regions of an annealed CMF. For micro-thick carbon nanotube bundles, θ T is found to have a relative change of >125% within a length of ~2 mm, uncovering strong metallic to semiconductive behavior change in space. In conclusion, our θ T scanning technique can be readily extended to nm-thick samples with μm scanning resolution to explore the distribution of θ T and provide a deep insight into the local electron conduction.

36 MATERIALS SCIENCE↗
DOE OSTI2024

How does negative triangularity mitigate ITG turbulence and transport?

Improved confinement in negative triangularity (NT) experiments is attributed to reduced fluxes driven by micro-turbulence. The physical mechanism of why thermal confinement improves in NT relative to PT is unknown. This study employs gyrokinetic flux tube simulations using the GENE code with local Miller equilibrium to elucidate the physical mechanisms behind the beneficial effects of NT flux surface shapes. The focus is on collisionless ion temperature gradient (ITG) driven turbulence with adiabatic electrons. The kinetic profiles are held fixed across a scan of triangularity values, thus enabling comparisons on a level playing field. The reduced linear growth rates for NT is shown to be due to a reduced eigenmode averaged magnetic drift frequency and a wider, stronger negative local magnetic shear region about the outboard mid-plane. The nonlinear heat flux is lower for NT than that for PT, due to reduced radial correlation length and increased correlation time (τ c ) of fluctuations. These, in turn, are due to a comparatively higher level of self-generated zero-frequency E × B zonal shearing rate ω E in NT as compared to PT. Though the linear zonal potential residual is lower for NT, the nonlinearly generated E × B zonal shearing rate is higher for NT than for PT. This outcome is linked to the distinctive features of the radial wavenumber spectra of the zonal potential and the zonal shearing rate. The dimensionless parameter ω E τ c is suggested as a figure of merit. This is higher for NT than for PT. Thus, the reduced heat diffusivity for NT is linked to increased ω E τ c . Self-generated temperature corrugations (i.e. zonal temperature gradients) are much weaker than the background mean temperature gradient. Nevertheless, temperature corrugations are more pronounced in NT than in PT.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Resolving the Solvation Structure and Transport Properties of Aqueous Zinc Electrolytes from Salt-in-Water to Water-in-Salt Using Neural Network Potential

Zn Cl 2 solutions are promising electrolytes for aqueous zinc-ion batteries. Here, we report a joint computational and experimental study of the structural and dynamic properties of aqueous Zn Cl 2 electrolytes with concentrations ranging from salt-in-water to water-in-salt (WIS). By developing a neural network potential (NNP) model, we perform molecular dynamics (MD) simulations with accuracy but at much larger lengths and longer timescales. The NNP predicted structures are validated by the structure factors measured by X-ray total scattering experiments. The MD trajectories provide a comprehensive and quantitative picture of the Zn 2 + solvation shell structures. Additionally, we find that the O − H covalent bonds in water are strengthened with increasing salt concentration, thus expanding the electrochemical stability window of aqueous electrolytes. In terms of dynamic properties, the calculated and experimentally measured conductivities are in good agreement. Through the analysis of the calculated cation transference number, we propose a three-stage charge carrier transport mechanism with increasing concentration: independent ion transport, strongly correlated ion transport, and small positive charge carrier diffusion through negatively charged polymeric clusters. Our study provides fundamental atomic scale insights into the structure and transport properties of the Zn Cl 2 electrolyte that can aid the optimization and development of WIS electrolytes. Published by the American Physical Society 2025

25 ENERGY STORAGE↗
DOE OSTI2022

Development of a creep model informed by lower-length scale simulations to simulate creep in doped UO 2

Using molecular dynamics, we predict information at the atomistic scale used to develop a mechanistic UO 2 creep model for use in higher length-scale fuel performance codes. The ultimate objective of the model is to better describe the grain size dependence and therefore impact of doping on creep rates in UO 2 . In a previous NEAMS milestone, we found that Nabarro-Herring (bulk diffusional) creep was too low to capture the experimentally observed creep rates in standard UO 2 . Moreover, in that milestone, other mechanisms were explored, such as Coble (grain boundary) creep and dislocation climb, with each mechanism exhibiting different grain size dependencies. Again, these were orders of magnitude too low to describe the experimental creep rates. In this work, we address the previous assumptions made for the Coble creep mechanism by investigating the diffusivity of various defects at grain boundaries in UO 2 and, critically, to determine if enhanced grain boundary diffusivity allows the model to better reproduce experimental results. The diffusivity as a function of temperature for different concentrations of uranium vacancies and interstitials for bulk UO 2 have also been examined using cluster dynamics. Furthermore, using a concentration dependent segregation model, the concentration of defects at the grain boundary were predicted. This atomistic data was then input into the various creep mechanisms and the creep rates compared to the empirical MATPRO correlation (used in BISON) and experiment.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Characterization of electromigration-induced short-range stress development in Al(0.25 at. % Cu) conductor line

Scanning x-ray microbeam topography and fluorescence experiments were conducted in situ to study the electromigration behavior of a 0.5 μm thick, 10 μm wide, and 200 μm long Al(0.25 at. % Cu) conductor line with 1.5 μm-thick SiO2 passivation on a single crystal Si substrate. The strain sensitivity of x-ray topography measurement allowed detailed examination of the electromigration-induced stress distribution and evolution in the conductor line in response to the depletion of Cu solute early in the electromigration process. Upon electromigration at 0.4 MA/cm2 and 303 °C, a short-range stress gradient was quickly induced by Al migration in the Cu-depleted cathode region to counteract further Al flow. The stress gradient was fully developed during the 5.3 h incubation time, extending over the critical Blech length of about 66 μm from the cathode end. Plastic deformation then occurred at the downstream end of the Cu-depleted region. The preferential electromigration of Cu did not cause detectable stress change outside the Cu-depleted region, except for the significant stress development from the Al2Cu precipitation at the anode end which appeared to initiate the fracture in the passivation. Preliminary finite difference modeling was undertaken to simulate the experimental observations, from which important parameters dictating electromigration in Al(Cu) line were extracted: an apparent effective valence of −5.6 and −1.9 for Cu and Al in Al(Cu), respectively, and a critical Cu concentration of 0.16 at. % above which Al grain boundary diffusion is effectively blocked.

Physics↗
DOE OSTI2021

Fast multiscale contrast independent preconditioners for linear elastic topology optimization problems

The goal of this work is to present a fast and viable approach for the numerical solution of the high-contrast state problems arising in topology optimization. The optimization process is iterative, and the gradients are obtained by an adjoint analysis, which requires the numerical solution of large high-contrast linear elastic problems with features spanning several length scales. The size of the discretized problems forces the utilization of iterative linear solvers with solution time dependent on the quality of the preconditioner. The lack of clear separation between the scales, as well as the high-contrast, imposes severe challenges on the standard preconditioning techniques. Thus, here we propose new methods for the high-contrast elasticity equation with performance independent of the high-contrast and the multi-scale structure of the elasticity problem. The solvers are based on two-levels domain decomposition techniques with a carefully constructed coarse level to deal with the high-contrast and multi-scale nature of the problem. The construction utilizes spectral equivalence between scalar diffusion and each displacement block of the elasticity problems and, in contrast to previous solutions proposed in the literature, is able to select the appropriate dimension of the coarse space automatically. The new methods inherit the advantages of domain decomposition techniques, such as easy parallelization and scalability. Finally, the presented numerical experiments demonstrate the excellent performance of the proposed methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Dry electrodes with a printed cellulose–graphene ink for low-profile strain sensors in electromyography

Dihydrolevoglucosenone, commonly known as Cyrene, is a renewable and fully biodegradable cellulose-waste derived, environmentally friendly solvent, presenting a non-toxic alternative to N-methyl-2-pyrrolidone (NMP). Currently, solution-based processing of graphene and other similar van der Waals solids favor toxic solvents such as NMP, limiting their use for biosensing. However, with the use of Cyrene, bio-compatible printable devices are possible, and studies have already demonstrated its use in temperature and other biosensing methods through screen-printing. Screen-printing unfortunately often requires masks that constrain the minimum acquirable feature size to be above hundreds of centimeters and wastes material, adding to process complexity and cost. Conversely, inkjet-printing is an attractive alternative for the maskless patterning of hierarchically assembled structures, with micron length scales attainable. Graphene's high conductivity positions it ideally for long-wear sensors such as dry electrodes or respiration monitors. Here, we demonstrate the potential of Cyrene-based graphene inks through few-layer inkjet printing on flexible substrates for the first time, to produce non-toxic conductors toward a strain-mediated mechanism for biosensing, used to detect bodily motion for wearable electronics. The challenges overcome in this study include engineering ink chemistry and printing parameters such that Cyrene's relatively high viscosity compared to typical inkjet solvents, still allows for droplet ejection in a conventional material printer, yielding well-resolved clean line-edges in contrast to other solvents that exhibit diffuse line-edges possibly from stray droplets and ink-splashing. Temperature-dependent transport measurements on the inkjet-printed Cyrene-based graphene films showed the conductivity to be largely temperature-invariant but at lower temperatures below 100 K, conductivity decreased, likely as a result of increased inter-membrane separation arising from thermal contraction. Additionally, temperature-dependent Raman spectroscopy showed the red-shift in the G-band, 2D-band and D-band peaks, as temperature increased. As a result, by validating flexion motion detection of the proximal interphalangeal joint demonstrated in this study, our work is the first of its kind to successfully additively manufacture inkjet-printed Cyrene-based graphene strain sensors on flexible substrates for bio-sensing and wearables.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Intermediate scattering functions of a rigid body monoclonal antibody protein in solution studied by dissipative particle dynamic simulation

In the past decade, there was increased research interest in studying internal motions of flexible proteins in solution using Neutron Spin Echo (NSE) as NSE can simultaneously probe the dynamics at the length and time scales comparable to protein domain motions. However, the collective intermediate scattering function (ISF) measured by NSE has the contributions from translational, rotational, and internal motions, which are rather complicated to be separated. Widely used NSE theories to interpret experimental data usually assume that the translational and rotational motions of a rigid particle are decoupled and independent to each other. To evaluate the accuracy of this approximation for monoclonal antibody (mAb) proteins in solution, dissipative particle dynamic computer simulation is used here to simulate a rigid-body mAb for up to about 200 ns. The total ISF together with the ISFs due to only the translational and rotational motions as well as their corresponding effective diffusion coefficients is calculated. The aforementioned approximation introduces appreciable errors to the calculated effective diffusion coefficients and the ISFs. For the effective diffusion coefficient, the error introduced by this approximation can be as large as about 10% even though the overall agreement is considered reasonable. Thus, we need to be cautious when interpreting the data with a small signal change. In addition, the accuracy of the calculated ISFs due to the finite computer simulation time is also discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Structure and Electronic Properties of InSb Nanowires Grown in Flexible Polycarbonate Membranes

A dense array of vertically aligned indium antimonide (InSb) nanowires with high aspect ratio (diameter 150 nm, length 20 μ m) were grown in the pores of a track-etched polycarbonate membrane via a one-step electrochemical method. There are several reports on InSb nanowire growth in the pores of a mechanically rigid, nano-channel alumina template (NCA), where nanowire growth occurs in the pores of the NCA. This work on InSb nanowire growth in pores of track-etched polycarbonate (PC) membrane sheds light on the various factors that affect nucleation and nanowire growth. The average length and diameter of the as-grown nanowires was about 10 μ m and 150 nm, respectively. Two possible mechanisms accounting for two different morphologies of the as-grown nanowires are proposed. The polycrystallinity observed in some of the nanowires is explained using the 3D ‘nucleation-coalescence’ mechanism. On the other hand, single crystal nanowires with a high density of twin defects and stacking faults grow epitaxially by a two-dimensional (2D) nucleation/growth mechanism. To assess the electrical quality of the nanowires, two- and four-terminal devices were fabricated using a single InSb nanowire contacted by two Ni electrodes. It was found that, at low bias, the ohmic current is controlled by charge diffusion from the bulk contacts. On the other hand, at high bias, the effects of space charge limited current (SCLC) are evident in the current–voltage behavior, characteristic of transport through structures with reduced electrostatic screening. A cross-over from ohmic to SCLC occurs at about 0.14 V, yielding a free carrier concentration of the order of 10 14 cm -3 .

36 MATERIALS SCIENCE↗
DOE OSTI2023

Computational Study on the Effect of Thermal Boundary Conditions and Axial Aspect Ratio on Catalytic Oxidative Coupling of Methane

Packed-bed reactor models with a large reactor-diameter-to-particle-diameter ratio (N > 10) are commonly used in catalysis research. Their use is underpinned by strong simplifying assumptions including unidirectional plug flow, negligible radial gradients, small axial diffusion, and isothermal operation. However, the complex exothermic–endothermic nature of reactions like oxidative coupling of methane (OCM) coupled with heat transfer from the catalyst surface and convective heat losses from the reactor wall provides an interesting case study to evaluate packed-bed assumptions and understand axial and radial temperature variations that may impact conversion and selectivity. This computational study investigates the effect of thermal boundary conditions on catalyst performance in laboratory OCM reactors (N = ~13) through computational fluid dynamics (CFD) simulation. Further, the 1-D model typically used for modeling plug flow reactors assumes negligible radial temperature and species concentration gradients. This assumption is evaluated for OCM reactors over a range of geometric axial aspect ratios, defined as the reactor length divided by the diameter, comparing the 1-D with 3-D CFD results. 1-D model predictions of catalyst performance start to deviate from the 3-D model predictions as the axial aspect ratio decreases, specifically as diameter increases at constant length. In the work, the 1-D model assumption is only found to be valid for axial aspect ratios higher than 1.6, a result considered to be relevant for all exothermic catalytic reactors.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Nanometer-Scale Acoustic Wave Packets Generated by Stochastic Core-Level Photoionization Events

We demonstrate that the absorption of femtosecond hard x-ray pulses excites quasispherical, high-amplitude, and high-wave-vector coherent acoustic phonon wave packets using an all hard-x-ray pump-probe scattering experiment. The time- and momentum-resolved diffuse scattering signal is consistent with an ensemble of 3D strain wave packets induced by the rapid electron cascade dynamics following photoionization at uncorrelated excitation centers. We quantify key parameters of this process, including the localization size of the stress field and the photon energy conversion efficiency into elastic energy. The parameters are determined by the photoelectron and Auger electron cascade dynamics, as well as the electron-phonon interaction. In particular, we obtain the localization size of the observed strain wave packet to be 1.5 and 2.5 nm for bulk SrTiO 3 and KTaO 3 single crystals, respectively. The results provide crucial information on the mechanism of x-ray energy deposition into matter and shed light on the shortest collective length scales accessible to coherent acoustic phonon generation using x-ray excitation. Published by the American Physical Society 2024

36 MATERIALS SCIENCE↗
DOE OSTI2019

Big Wheel Farm: Farmland Scour Reduction

Flood irrigation benefits from low infrastructure costs and maintenance but the scour near the weirs can cause channeling of the flow preventing the water from evenly dispersing across the field. Using flow obstructions in front of the weir could reduce be a low cost solution to reduce the scour. The mitigation strategy was to virtually simulate the effects of various geometric changes to the morphology (e.g. holes and bumps) in front of the weir as a means to diffuse the high intensity flow coming from the gate. After running a parametric study for the dimensions of the shapes that included a Gaussian, semi-circle, and rectangle; a Gaussian-hole in front of the gates showed the most promise to reduce farm field shear-stresses with the added benefit of being easy to construct and implement in practice. Further the simulations showed that the closer the Gaussian-hole could be placed to the gate the sooner the high shear stress could be reduced. To realize the most benefit from this mitigation strategy, it was determined that the maximum depth of the Gaussian-hole should be 0.5 m. The width of the hole in the flow direction and the length of the Gaussian-hole normal to the flow should be 0.5 m and 3 m respectively as measured by the full width at half maximum.

42 ENGINEERING↗
DOE OSTI2020

Using Strain in the Vadose Zone During Sinusoidal-Rate Well Tests to Characterize an Aquifer

Well tests with sinusoidal variation in pumping rate are appealing because of the ability to filter competing signals and they also can be set up with no net removal of water, which is an advantage at contaminated sites. Current methods of interpreting sinusoidal-rate well tests analyze the phase lag and amplitude of pressure signals in the aquifer. This approach is restricted by the availability of monitoring wells. An alternative approach could be to measure the stain in the vadose zone instead of pressure in the aquifer. The goal of this research is to evaluate the feasibility of using signals from in-situ sensors in the vadose zone to characterize sinusoidal-rate well tests. Well tests with a variety of periods and amplitudes were performed using a well completed in an unconfined saprolite aquifer with a water table at 9m depth at our test site near Clemson, South Carolina. Vertical strain was measured every meter between 1 to 7 m using a distributed optical fiber sensing technique called coherence length gated microwave photonics interferometry (CMPI). A novel packaging system was designed to couple the optical fiber to the saprolite. The results show that strain measured in the vadose zone during sinusoidal-rate pumping varies periodically with amplitudes from 0.1 to 3 microstrain, depending on the amplitude of the sinusoidal pumping rate. Hydraulic heads at similar radial distances were also periodic with changing amplitudes ranging up to ±0.2 meters. The strain in the vadose zone responds to pressure in the underlying aquifer—the strain is tensile when the pressure is falling, and compressive when the pressure is rising. The phase lag between the pumping rate and the strain was 3 min with a pumping period of 15 minutes, and this gives a hydraulic diffusivity of 0.06 m2/s. This is the same diffusivity obtained by analyzing a constant rate pumping test, and it is slightly larger than the diffusivity (0.044 m2/s) obtained by analyzing the phase lag of the pressure in a nearby monitoring well. Strain measurements in the vadose zone could be made with much less cost and effort than installing a deep monitoring well. These results suggest that strain measurements in the vadose zone could be used to supplement monitoring wells, thereby increasing the resolution, and reducing costs of analyzing well tests.

well test, sinusoidal, periodic pumping, strain↗