Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “loss function regularization”

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.

53 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Physics-guided dual implicit neural representations for source separation

Significant challenges exist in efficient data analysis of most advanced experimental and observational techniques because the collected signals often include unwanted contributions, such as background and signal distortions, that can obscure the physically relevant information of interest. To address this, we have developed a self-supervised machine-learning approach for source separation using a dual implicit neural representation framework that jointly trains two neural networks: one for approximating distortions of the physical signal of interest and the other for learning the effective background contribution. Our method learns directly from the raw data by minimizing a reconstruction-based loss function without requiring labeled data or pre-defined dictionaries. We demonstrate the effectiveness of our framework by considering a challenging case study involving large-scale simulated, as well as experimental, momentum-energy-dependent inelastic neutron scattering data in a four-dimensional parameter space, characterized by heterogeneous background contributions and unknown distortions to the target signal. The method is found to successfully separate physically meaningful signals from a complex or structured background even when the signal characteristics vary across all four dimensions of the parameter space. An analytical approach that informs the choice of the regularization parameter is presented. Our method offers a versatile framework for addressing source separation problems across diverse domains, ranging from superimposed signals in astronomical measurements to structural features in biomedical image reconstructions.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Enhancing Interpretability in Generative Modeling: Statistically Disentangled Latent Spaces Guided by Generative Factors in Scientific Datasets

This study addresses the challenge of statistically extracting generative factors from complex, high-dimensional datasets in unsupervised or semi-supervised settings. We investigate encoder-decoder-based generative models for nonlinear dimensionality reduction, focusing on disentangling low-dimensional latent variables corresponding to independent physical factors. Introducing Aux-VAE, a novel architecture within the classical Variational Autoencoder framework, we achieve disentanglement with minimal modifications to the standard VAE loss function by leveraging prior statistical knowledge through auxiliary variables. These variables guide the shaping of the latent space by aligning latent factors with learned auxiliary variables. We validate the efficacy of Aux-VAE through comparative assessments on multiple datasets, including astronomical simulations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Inverter Reliability Estimation for Advanced Inverter Functionality

In the near future, grid operators are expected to regularly use advanced distributed energy resource (DER) functions, defined in IEEE 1547-2018, to perform a range of grid-support operations. Many of these functions adjust the active and reactive power of the device through commanded or autonomous modes, which will produce new stresses on the grid-interfacing power electronics components, such as DC/AC inverters. In previous work, multiple DER devices were instrumented to evaluate additional component stress under multiple reactive power setpoints. We utilize quasi-static time-series simulations to determine voltage-reactive power mode (volt-var) mission profile of inverters in an active power system. Mission profiles and loss estimates are then combined to estimate the reduction of the useful life of inverters from different reactive power profiles. It was found that the average lifetime reduction was approximately 0.15% for an inverter between standard unity power factor operation and the IEEE 1547 default volt-var curve based on thermal damage due to switching in the power transistors. For an inverter with an expected 20-year lifetime, the 1547 volt-var curve would reduce the expected life of the device by 12 days. This framework for determining an inverter's useful life from experimental and modeling data can be applied to any failure mechanism and advanced inverter operation.

component degradation↗
NASA NTRS1980

Long nonlinear waves in stratified shear flows

The propagation of finite-amplitude internal waves in a shear flow is considered for wavelengths that are long compared to the shear-layer thickness. Both singular and regular modes are investigated, and the equation governing the amplitude evolution is derived. The theory is generalized to allow for a radiation condition when the region outside the stratified shear layer is unbounded and weakly stratified. In this case, the evolution equation contains a damping term describing energy loss by radiation which can be used to estimate the persistence of solitary waves or nonlinear wave packets in realistic environments. A continuous three-layer model is studied in detail and closed-form expressions are obtained for the phase speed and the coefficients of the nonlinear and dispersive terms in the amplitude equation as a function of Richardson number.

Maslowe, S. A.↗
DOE OSTI2022

Automating Component-Level Stress Measurements for Inverter Reliability Estimation

In the near future, grid operators are expected to regularly use advanced distributed energy resource (DER) functions, defined in IEEE 1547-2018, to perform a range of grid-support operations. Many of these functions adjust the active and reactive power of the device through commanded or autonomous operating modes which induce new stresses on the power electronics components. In this work, an experimental and theoretical framework is introduced which couples laboratory-measured component stress with advanced inverter functionality and derives a reduction in useful lifetime based on an applicable reliability model. Multiple DER devices were instrumented to calculate the additional component stress under multiple reactive power setpoints to estimate associated DER lifetime reductions. A clear increase in switch loss was demonstrated as a function of irradiance level and power factor. This is replicated in the system-level efficiency measurements, although magnitudes were different—suggesting other loss mechanisms exist. Using an approximate Arrhenius thermal model for the switches, the experimental data indicate a lifetime reduction of 1.5% when operating the inverter at 0.85 PF—compared to unity PF—assuming the DER failure mechanism thermally driven within the H-bridge. If other failure mechanisms are discovered for a set of power electronics devices, this testing and calculation framework can easily be tailored to those failure mechanisms.

Flicker, Jack (ORCID:0000000212232154)↗
NASA NTRS1998

High sustained +Gz acceleration: physiological adaptation to high-G tolerance

Since the early 1940s, a significant volume of research has been conducted in an effort to describe the impact of acute exposures to high-G acceleration on cardiovascular mechanisms responsible to maintaining cerebral perfusion and conscious in high performance aircraft pilots during aerial combat maneuvers. The value of understanding hemodynamic characteristics that underlie G-induced loss of consciousness has been instrumental in the evolution of optimal technology development (e.g., G-suits, positive pressure breathing, COMBAT EDGE, etc.) and pilot training (e.g., anti-G straining maneuvers). Although the emphasis of research has been placed on the development of protection against acute high +Gz acceleration effects, recent observations suggest that adaptation of cardiovascular mechanism associated with blood pressure regulation may contribute to a protective 'G-training' effect. Regular training at high G enhances G tolerance in humans, rats, guinea pigs, and dogs while prolonged layoff from exposure in high G profiles (G-layoff) can result in reduced G endurance. It seems probable that adaptations in physiological functions following chronically-repeated high G exposure (G training) or G-layoff could have significant impacts on performance during sustained high-G acceleration since protective technology such as G-suits and anit-G straining maneuvers are applied consistently during these periods of training. The purpose of this paper is to present a review of new data from three experiments that support the notion that repeated exposure on a regular basis to high sustained +Gz acceleration induces significant physiological adaptations which are associated with improved blood pressure regulation and subsequent protection of cerebral perfusion during orthostatic challenges.

Review, Tutorial↗
DOE OSTI2022

Microbial community shifts correspond with suppression of decomposition 25 years after liming of acidic forest soils

Microbial community structure and function regularly covary with soil pH, yet effects of these interactions on soil carbon are rarely tested experimentally within natural ecosystems. We investigated the enduring (25 year) impacts of liming on microbial community structure and decomposition at an acidic northern hardwood forest, where experimental liming increased pH one unit and surprisingly doubled the organic carbon stocks of the forest floor. Here we show that this increase in carbon storage corresponded with restructuring of the bacterial and fungal communities that drive decomposition. In the Oe horizon, liming reduced the activities of five extracellular enzymes that mediate decomposition, while the Oa horizon showed an especially large (64%) reduction in the activity of a sixth, peroxidase, which is an oxidative enzyme central to lignocellulose degradation. Decreased enzyme activities corresponded with loss of microbial taxa important for lignocellulose decay, including large reductions in the dominant ectomycorrhizal genera Russula and Cenococcum, saprotrophic and wood decaying fungi, and Actinobacteria (Thermomonosporaceae). These results demonstrate the importance of pH as a dominant regulator of microbial community structure and illustrate how changes to this structure can produce large, otherwise unexpected increases in carbon storage in forest soils.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2003

Universal Responses of Cyclic-Oxidation Models Studied

Oxidation is an important degradation process for materials operating in the high-temperature air or oxygen environments typical of jet turbine or rocket engines. Reaction of the combustion gases with the component material forms surface layer scales during these oxidative exposures. Typically, the instantaneous rate of reaction is inversely proportional to the existing scale thickness, giving rise to parabolic kinetics. However, more realistic applications entail periodic startup and shutdown. Some scale spallation may occur upon cooling, resulting in loss of the protective diffusion barrier provided by a fully intact scale. Upon reheating, the component will experience accelerated oxidation due to this spallation. Cyclic-oxidation testing has, therefore, been a mainstay of characterization and performance ranking for high-temperature materials. Models simulate this process by calculating how a scale spalls upon cooling and regrows upon heating (refs. 1 to 3). Recently released NASA software (COSP for Windows) allows researchers to specify a uniform layer or discrete segments of spallation (ref. 4). Families of model curves exhibit consistent regularity and trends with input parameters, and characteristic features have been empirically described in terms of these parameters. Although much insight has been gained from experimental and model curves, no equation has been derived that can describe this behavior explicitly as functions of the key oxidation parameters.

Smialek, James L.↗
NASA NTRS2015

Comparison of Spares Logistics Analysis Techniques for Long Duration Human Spaceflight

As the durations and distances involved in human exploration missions increase, the logistics associated with the repair and maintenance becomes more challenging. Whereas the operation of the International Space Station (ISS) depends upon regular resupply from the Earth, this paradigm may not be feasible for future missions. Longer mission durations result in higher probabilities of component failures as well as higher uncertainty regarding which components may fail, and longer distances from Earth increase the cost of resupply as well as the speed at which the crew can abort to Earth in the event of an emergency. As such, mission development efforts must take into account the logistics requirements associated with maintenance and spares. Accurate prediction of the spare parts demand for a given mission plan and how that demand changes as a result of changes to the system architecture enables full consideration of the lifecycle cost associated with different options. In this paper, we utilize a range of analysis techniques - Monte Carlo, semi-Markov, binomial, and heuristic - to examine the relationship between the mass of spares and probability of loss of function related to the Carbon Dioxide Removal System (CRS) for a notional, simplified mission profile. The Exploration Maintainability Analysis Tool (EMAT), developed at NASA Langley Research Center, is utilized for the Monte Carlo analysis. We discuss the implications of these results and the features and drawbacks of each method. In particular, we identify the limitations of heuristic methods for logistics analysis, and the additional insights provided by more in-depth techniques. We discuss the potential impact of system complexity on each technique, as well as their respective abilities to examine dynamic events. This work is the first step in an effort that will quantitatively examine how well these techniques handle increasingly more complex systems by gradually expanding the system boundary.

Owens, Andrew↗
DOE OSTI2024

Robust wind farm layout optimization

Wake interactions in wind farms cause losses in annual energy production (AEP) on the order of 10%. Wind farm designers optimize the layout of the farm to mitigate wake losses, especially in the dominant site-specific wind directions. As wind turbines and wind farms grow in scale, optimization becomes more complex. Offshore wind farms regularly comprise more than 100 wind turbines and are characterized by complex boundaries due to shipping lanes, neighboring wind farms, and other constraints. Layout optimization methods are broadly split between gradient-based and gradient-free approaches. Gradient-based approaches can converge quickly and perform well for smaller, academic problems but are often sensitive to initial conditions and tuning parameters and require expert knowledge to use. On the other hand, gradient-free approaches can be more robust to problem complexities. We present a robust layout optimization approach based on a random search algorithm. The algorithm is intended for those who are not optimization experts and has few tuning parameters that need specification to achieve satisfactory results. Unlike off-the-shelf methods, which use generally available, non-domain-specific optimization routines that accept as inputs an optimization function and constraint definitions, this approach takes advantage of the relative computational costs of the different evaluations by evaluating cheaper computations first (boundary and minimum distance constraints) and running expensive AEP evaluations only if all other checks pass. Moreover, an outer genetic algorithm allows multiple solutions to evolve in parallel, enabling rapid solution development on high-performance computers. We discuss the relative ease of selecting necessary tuning parameters and demonstrate the efficacy of the genetic random search on a complex layout problem consisting of placing 70 turbines in a nonconvex and unconnected boundary region.

17 WIND ENERGY↗
NASA NTRSDate not supplied

Developing A Continuous Ozone Record Through the SAGE and Aura Missions With NASA Reanalysis Products

During the last quarter of the 20th century, the Stratospheric Aerosol and Gas Experiment (SAGE) missions were crucial in monitoring the loss and subsequent recovery of the stratospheric ozone layer. Due to the employed solar occultation and self-calibration method, the SAGE monitors have produced stable data throughout the lifetime of each instrument. However, over ten years passed between the end of the SAGE II and SAGE III/M3M missions in 2005 and the launch of SAGE III/ISS instrument in 2017, leaving a gap in the data that must be bridged in order to assess trends in the ozone record. The Modern-Era Retrospective analysis for Research and Applications, version 2 (MERRA-2) reanalysis product, with output available starting in 1980, is an attractive candidate for trend analysis due to the statistically optimized combination of multiple observing systems and the regular temporal and spatial coverage. However, changes in the assimilated observation systems can introduce discontinuities within the MERRA-2 ozone record, such as in 2004 when the MERRA-2 system shifted from assimilating ozone retrievals collected by SBUV instruments to those collected by instruments onboard the Aura satellite. In this study, we explore using the SAGE II record as a transfer function to develop a stable reanalysis data product, suitable for trend analysis, from the start of the SAGE II record in 1984 through the present. We follow the procedure outlined by Wargan et al. (2018) to address discontinuities in the MERRA-2 ozone dataset at the 2004 transition and during the Aura record. SAGE II ozone profiles are used to correct discontinuities in upper stratospheric ozone associated with changes in the MERRA-2 meteorological observing system in 1998 and 1995. We will then assess the relative performance of the data from different SAGE sensors using the resulting bias-corrected MERRA-2 ozone fields.

Pamela Wales↗
DOE OSTI2022

A fourth-order phase-field fracture model: Formulation and numerical solution using a continuous/discontinuous Galerkin method

Modeling crack initiation and propagation in brittle materials is of great importance to be able to predict sudden loss of load-carrying capacity and prevent catastrophic failure under severe dynamic loading conditions. Second-order phase-field fracture models have gained wide adoption given their ability to capture the formation of complex fracture patterns, e.g. via crack merging and branching, and their suitability for implementation within the context of the conventional finite element method. Higher-order phase-field models have also been proposed to increase the regularity of the exact solution and thus increase the spatial convergence rate of its numerical approximation. However, they require special numerical techniques to enforce the necessary continuity of the phase field solution. In this paper, we derive a fourth-order phase-field model of fracture in two independent ways; namely, from Hamilton’s principle and from a higher-order micromechanics-based approach. The latter approach is novel, and provides a physical interpretation of the higher-order terms in the model. In addition, we propose a continuous/discontinuous Galerkin (C/DG) method for use in computing the approximate phase-field solution. This method employs Lagrange polynomial shape functions to guarantee -continuity of the solution at inter-element boundaries, and enforces the required regularity with the aid of additional variational and interior penalty terms in the weak form. Finally, the phase-field equation is coupled with the momentum balance equation to model dynamic fracture problems in hyper-elastic materials. Two benchmark problems are presented to compare the numerical behavior of the C/DG method with mixed finite element methods.

42 ENGINEERING↗
NASA NTRS2000

Multiresolution With Super-Compact Wavelets

The solution data computed from large scale simulations are sometimes too big for main memory, for local disks, and possibly even for a remote storage disk, creating tremendous processing time as well as technical difficulties in analyzing the data. The excessive storage demands a corresponding huge penalty in I/O time, rendering time and transmission time between different computer systems. In this paper, a multiresolution scheme is proposed to compress field simulation or experimental data without much loss of important information in the representation. Originally, the wavelet based multiresolution scheme was introduced in image processing, for the purposes of data compression and feature extraction. Unlike photographic image data which has rather simple settings, computational field simulation data needs more careful treatment in applying the multiresolution technique. While the image data sits on a regular spaced grid, the simulation data usually resides on a structured curvilinear grid or unstructured grid. In addition to the irregularity in grid spacing, the other difficulty is that the solutions consist of vectors instead of scalar values. The data characteristics demand more restrictive conditions. In general, the photographic images have very little inherent smoothness with discontinuities almost everywhere. On the other hand, the numerical solutions have smoothness almost everywhere and discontinuities in local areas (shock, vortices, and shear layers). The wavelet bases should be amenable to the solution of the problem at hand and applicable to constraints such as numerical accuracy and boundary conditions. In choosing a suitable wavelet basis for simulation data among a variety of wavelet families, the supercompact wavelets designed by Beam and Warming provide one of the most effective multiresolution schemes. Supercompact multi-wavelets retain the compactness of Haar wavelets, are piecewise polynomial and orthogonal, and can have arbitrary order of approximation. The advantages of the multiresolution algorithm are that no special treatment is required at the boundaries of the interval, and that the application to functions which are only piecewise continuous (internal boundaries) can be efficiently implemented. In this presentation, Beam's supercompact wavelets are generalized to higher dimensions using multidimensional scaling and wavelet functions rather than alternating the directions as in the 1D version. As a demonstration of actual 3D data compression, supercompact wavelet transforms are applied to a 3D data set for wing tip vortex flow solutions (2.5 million grid points). It is shown that high data compression ratio can be achieved (around 50:1 ratio) in both vector and scalar data set.

Lee, Dohyung↗
DOE OSTI2021

Shake loss intensities in x-ray photoelectron spectroscopy: Theory, experiment, and atomic composition accuracy for MgO and related compounds

The relative intensities of XPS core levels, scaled by their photoionization cross sections, are regularly used to determine sample atomic composition. Cross sections, however, give the intensity to all possible final states for the core ionizations, not just to the main peak. This includes all intrinsic satellite structure (shake states and, for open shell systems, the different ionic multiplets). In practice, for solids, this is usually experimentally impossible to determine accurately because such a satellite structure sits on the inelastically scattered electron background and cannot be easily separated. Therefore, usually, only the intensity of the main peak is used. This limits the ultimate possible accuracy of XPS composition determination. The purpose of the present paper is to examine the contributions that a theoretical analysis of losses of intensity can make to improve quantitation. For an MgO single crystal, we show that the correct stoichiometry of 1:1 can be recovered using the theoretical analysis of the experimental MgO peak ratio intensities. For materials with a sufficient bandgap for the XPS main peaks to be separated from the scattered background, the intensity of main peaks can often be accurately determined. Thus, if one uses theory to calculate that fraction of the total intensity lost from a main peak into all its satellite structure, the intensity of just main peaks could then be used to more accurately determine relative atom % composition. This work tests this approach using a single crystal MgO (50% Mg, 50% O) standard. Ab initio electronic structure theory of representative MgO clusters is used to determine Hartree–Fock wave functions for the ground state and final ionized states corresponding to the main Mg 2p and O1s XPS peaks of the oxide. The sudden approximation, SA, is used to determine the fractional losses from these main peaks to shake satellites, which is found to be greater for O1s than Mg2p. This results in predicted “apparent composition” for stoichiometric MgO of 55.2% Mg, 44.8% O instead of the true 50% Mg, 50% O. Equivalent theory for CaO results in a predicted apparent Ca value of 53.4%. Experimentally, using Mg2s or 2p intensity ratio to O1s, we find values between 52.2% and 56.0% Mg using two crystals and four different instrument electron pass energies. The average value of the measurements is 54.5% Mg when corrected for the presence of an adventitious carbon overlayer and slight surface hydroxide. Though this agreement with theory may be somewhat fortuitous, given the potential experimental errors, which are fully discussed, it is similar to that in our earlier study on LiF. We also present preliminary experimental data on Mg(OH) 2 and MgSO 4 , which show a similar trend of apparently higher than 50% Mg, but we have no theory values. We are not yet able to experimentally test for validation of the difference between apparent composition for MgO (55.2% Mg) and CaO (53.4% Ca), owing to significant carbonate formation at the surface of the single crystal CaO. Here, an important conclusion is that the theoretical determination of shake losses, obtained with ab initio wavefunctions and the SA, is likely to be a useful way to calibrate the accuracy and reliability of compositions obtained from XPS intensities and merits further study.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Relation of Observable Stellar Parameters to Mass-loss Rate of AGB Stars in the LMC

Using the data set of Riebel et al. for 6889 pulsating AGB stars in the LMC, we have derived formulae for mass-loss rate as a function of luminosity and pulsation period, or luminosity and mass, in three ways, for each of five subsets of data: fundamental-mode oxygen-rich stars, first-overtone-mode oxygen-rich stars, fundamental-mode carbon stars, first-overtone-mode carbon stars, and extreme carbon stars. Using the distribution of the stars in period versus luminosity and mass versus luminosity, we are able to derive a power-law fit to the dependence of mass-loss rate on those quantities. This results in formulae that reproduce observed mass-loss rates and are in general agreement with the expectation from mass-loss models that the mass-loss rate is highly sensitive to luminosity, mass, and pulsation period. In the process of carrying out this analysis we have found radius–mass–luminosity relations and examined pulsation–mass–radius relations using published evolutionary and pulsation models. These allow us to derive mass and radius from the observed quantities luminosity and pulsation period. We also derived new mass-loss rate versus color relations.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2011

Molecular Mechanisms of Circadian Regulation During Spaceflight

Disruption of the regular environmental circadian cues in addition to stringent and demanding operational schedules are two main factors that undoubtedly impact sleep patterns and vigilant performance in the astronaut crews during spaceflight. Most research is focused on the behavioral aspects of the risk of circadian desynchronization, characterized by fatigue and health and performance decrement. A common countermeasure for circadian re-entrainment utilizes blue-green light to entrain the circadian clock and mitigate this risk. However, an effective countermeasure targeting the photoreceptor system requires that the basic circadian molecular machinery remains intact during spaceflight. The molecular clock consists of sets of proteins that perform different functions within the clock machinery: circadian oscillators (genes whose expression levels cycle during the day, keep the pass of cellular time and regulate downstream effector genes), the effector or output genes (those which impact the physiology of the tissue or organism), and the input genes (responsible for sensing the environmental cues that allow circadian entrainment). The main environmental cue is light. As opposed to the known photoreceptors (rods and cones), the non-visual light stimulus is received by a subset of the population of retinal ganglion cells called intrinsically photosensitive retinal ganglion cells (ipRGC) that express melanopsin (opsin 4 -Opn4-) as the photoreceptor. We hypothesize that spaceflight may affect ipRGC and melanopsin expression, which may be a contributing cause of circadian disruption during spaceflight. To answer this question, eyes from albino Balb/cJ mice aboard STS-133 were collected for histological analysis and gene expression profiling of the retina at 1 and 7 days after landing. Both vivarium and AEM (animal enclosure module) mice were used as ground controls. Opn4 expression was analyzed by real time RT/qPCR and retinal sections were stained for Opn4 immunofluorescence. Opn4 was decreased (abrogated in one case) in retinas that concurrently showed higher evidence of oxidative stress. We propose that oxidative stress can lead to a decrease in melanopsin expression, likely via ipRGC loss or impairment, and thus, it can be a contributing factor to circadian disruption during spaceflight. Countermeasures contemplating the use of light should therefore be complemented with melanopsin expression maintenance and/or reduction in oxidative stress.

Zanello, Susana↗
DOE OSTI2023

Evaluation of Zero-Net-Rate Pumping Tests

Accurately estimating the distribution of aquifer properties is key to understanding contaminant movement in the subsurface. The distribution of aquifer properties is typically addressed using slug or constant-rate well tests, and the pros and cons of these tests are well known. Slug tests are appealing because they avoid removing contaminated water, but their results are affected by well skin and the small volume of displaced water limits the volume of aquifer that can be evaluated. Constant-rate well tests have the disadvantage of requiring disposal of potentially contaminated water, but they can generate properties that are more representative than slug tests, and they can be used to estimate well efficiency and storativity, which are difficult to characterize using slug tests. Periodic pumping tests are appealing because they have many of the advantages and few of the disadvantages of slug and constant-rate well tests. Periodic pumping tests involve cycling the pumping rate with a regular period and measuring the resulting response in monitoring wells. Some of these tests involve imposing a periodic rate on a constant mean rate. However, other tests involve moving water out and back into the well at rates that are balanced so the net rate after each period is zero. Some zero-net-rate (ZNR) well tests use rates that follow a sinusoidal pattern, whereas others use a square wave pattern that switches between a constant rate of pumping and a constant rate of injection to achieve a zero-net rate. ZNR well tests appear to be a useful compromise between slug and constant rate tests, but methods for conducting and analyzing the results from these tests have received limited evaluation. Periodic pumping tests are typically evaluated through recorded pressure responses to a disturbance source. Recent work has demonstrated that the strain field in the vadose zone shows a response to disturbances in the underlying aquifer. Various instrument designs are available that can record vertical and horizontal strain to very precise resolutions (10-9 ε). Measuring the strain in the vadose zone can present lower costs than measuring pressure in a monitoring well, so using strain could improve the resolution of ZNR tests, as well as make these tests potentially cheaper. A ZNR periodic pumping system was constructed and used to perform tests in the Clemson, South Carolina area. The system was designed to pump and inject with a periodic square-wave, so pumping occurs at a constant rate for half the period and is followed by injection at the same rate for half the period resulting in no net gain or loss of water from the aquifer. The system is designed to generate flow rates from 1 to 3.5 gallons per minute (gpm) with a capacity of 875-gallon per half period. Pressure and strain monitoring points are located in the vicinity of the pumping well. A series of 10 ZNR periodic pumping tests were conducted at the field site, with periods ranging from 2 to 540 minutes. Pressure data were collected and recorded in four monitoring wells in various locations around the pumping well. In addition, vertical and horizontal strain and tilt data were recorded in various locations around the pumping well. Traditional aquifer tests, including slug tests and constant-rate pumping tests were conducted at the field site to provide a baseline of aquifer parameter estimates that are compared the results of ZNR tests. Pressure data measured at monitoring wells during ZNR tests were analyzed utilizing an analytical solution (Streltsova, 1988) that assumes a confined aquifer and estimates hydraulic diffusivity and transmissivity using the time lag and amplitude of the pressure signal. The storativity can be separated out from the hydraulic diffusivity using the transmissivity estimate. The time-lag and amplitude of the pressure signals were estimated using a Fourier transform. The time-lag of the primary period of the pressure increases as a linear function of distance from the pumping well, and it increases as the square root of the period of the ZNR test over a range of periods spanning two orders of magnitudes. The first few harmonics follow similar trends, but the higher frequency harmonics diverge from this trend. Strain and tilt data from various instruments in the vadose and saturated zones are periodic during ZNR tests. The time lags of the strain components increase roughly linearly with distance from the well. Many of the time lags of the strain data from the vadose zone are similar to the time lags observed in the pressure data from a similar distance and pumping period. Theoretical experiments were performed to understand how pressure and strain responds to ZNR periodic pumping. The time-lags of the simulated pressure responses match the field data trends. The results were used to validate the Streltsova (1988) solution for periodic pumping to estimate hydraulic diffusivity. Hydraulic diffusivity was estimated to be 1.8 x 10-2 m2s-1 < Dh < 9.0 x 10-2 m2s-1 using pressure data from all the wells during ZNR tests. Transmissivity was estimated to be 0.8 x 10-4 m2s-1 < T < 4.3 x 10-4 m2s-1. Assuming Storativity = T/Dh gives 1.1 x 10-3 < S < 19 x10-3. Transmissivity is used to estimate the hydraulic conductivity, K, by dividing by the assumed aquifer thickness. Two constant-rate pumping tests were conducted and analyzed using two conceptual models: The Hantush (1961) solution was used to analyze data assuming confined conditions and the Neuman (1974) solution was used to analyze data assuming unconfined conditions. The Hantush (1961) solution estimated transmissivity to be 1.6 x 10-4 m2s-1 < T < 1.8 x 10-4 m2s-1 when using data from all the wells. Storativity was estimated to be 4.4 x 10-3 < S < 8.2 x 10-3. Assuming Dh = T/S gives 0.22 x 10-2 m2s-1 < Dh < 3.6 x 10-2 m2s-1. The Neuman (1974) solution estimated transmissivity to be 0.80 x 10-4 m2s-1 < T < 0.84 x 10-4 m2s-1. Total storativity (S + Sy) was estimated to be 22 x 10-3 < S < 110 x 10-3. Storativity excluding the Sy term was estimated to be 2.2 x 10-3 < S < 2.5 x 10-3. Assuming Dh = T/S gives 0.073 x 10-2 m2s-1 < S < 0.36 x 10-2 m2s-1. Two slug tests were performed on the pumping well, PW-2. The Bouwer and Rice (1976) solution for slug tests was used to analyze the data. This solution works for both confined and unconfined conditions, so both conceptual models were covered through the analysis. Hydraulic conductivity was estimated to be 1.95 x 10-6 ms-1 < K < 2.49 x 10-6 ms-1. Statistical comparison show generally no significant difference in the estimates of Dh, K, or S made using ZNR and conventional pumping tests.. However, the total storage estimated by the Neuman (1974) solution for unconfined settings is larger than that estimate using the Streltsove analysis applied to the ZNR tests. This is likely because the Neuman solution considers delayed yield from storage that the water table, whereas the Streltsova (1988) solution assumes confined conditions. The analysis also indicates that there is a statistically significant difference between parameters estimated with conventional slug tests and those measured with either constant-rate pumping tests or ZNR tests. Hydraulic diffusivity is estimated from the time-lag and the distance from the well. The time-lags of the strain data from the vadose zone are similar to the lags of the pressure data from a similar distance. This indicates that hydraulic diffusivities estimated from the strain data measured in the vadose zone would be similar to the estimates from the pressure data measured in the aquifer. There appear to be sign reversals in some of the strain data that were corrected to estimate the time lag. The dissolved oxygen (CO) concentration was measured during several ZNR tests to evaluate the feasibility of using the procedure to increase contaminant degradation kinetics that are related to DO. The results indicate that DO can be increased during the cycling of ZNR tests in some cases.

Smith-Jones, Austin↗