SEARCH · Search NASA

Results for “Divided differences”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,207 records · Page 67

All sources NASA NTRS DOE OSTI

NASA NTRS1993

A finite-element-based perturbation model for the rotordynamic analysis of shrouded pump impellers: Part 2: User's guide

This report describes the computational steps involved in executing a finite-element-based perturbation model for computing the rotor dynamic coefficients of a shrouded pump impeller or a simple seal. These arise from the fluid/rotor interaction in the clearance gap. In addition to the sample cases, the computational procedure also applies to a separate category of problems referred to as the 'seal-like' category. The problem, in this case, concerns a shrouded impeller, with the exception that the secondary, or leakage, passage is totally isolated from the primary-flow passage. The difference between this and the pump problem is that the former is analytically of the simple 'seal-like' configuration, with two (inlet and exit) flow-permeable stations, while the latter constitutes a double-entry / double-discharge flow problem. In all cases, the problem is that of a rotor clearance gap. The problem here is that of a rotor excitation in the form of a cylindrical whirl around the housing centerline for a smooth annular seal. In its centered operation mode, the rotor is assumed to give rise to an axisymmetric flow field in the clearance gap. As a result, problems involving longitudinal or helical grooves, in the rotor or housing surfaces, go beyond the code capabilities. Discarding, for the moment, the pre- and post-processing phases, the bulk of the computational procedure consists of two main steps. The first is aimed at producing the axisymmetric 'zeroth-order' flow solution in the given flow domain. Detailed description of this problem, including the flow-governing equations, turbulence closure, boundary conditions, and the finite-element formulation, was covered by Baskharone and Hensel. The second main step is where the perturbation model is implemented, with the input being the centered-rotor 'zeroth-order' flow solution and a prescribed whirl frequency ratio (whirl frequency divided by the impeller speed). The computational domain, in the latter case, is treated as three dimensional, with the number of computational planes in the circumferential direction being specified a priori. The reader is reminded that the deformations in the finite elements are all infinitesimally small because the rotor eccentricity itself is a virtual displacement. This explains why we have generically termed the perturbation model the 'virtually' deformable finite-element category. The primary outcome of implementing the perturbation model is the tangential and radial components, F(sub theta)(sup *) and F(sub r)(sup *) of the fluid-exerted force on the rotor surface due to the whirling motion. Repetitive execution of the perturbation model subprogram over a sufficient range of whirl frequency ratios, and subsequent interpolation of these fluid forces, using the least-square method, finally enable the user to compute the impeller rotor dynamic coefficients of the fluid/rotor interaction. These are the direct and cross-coupled stiffness, damping, and inertia effects of the fluid/rotor interaction.

Baskharone, Erian A.↗

NASA NTRS2017

Biosignatures in the Context of Low Energy Flux

Many of the features that are thought of as biosignatures - including the mediation of chemical and physical processes with speed, specificity, and selectivity - result directly or indirectly from life's unique capability to mediate and direct energy flux. As such, it is important to consider the impact that differences in energy flux may have on the quantity and quality of evidence for life. Earth differs from every other body in our solar system in the magnitude of biologically-usable energy flux into a liquid water environment. On a global basis, the capture of light energy into photosynthesis and the flux of chemical energy represented in the products of that photosynthesis (organic material + O2) are about six and four orders of magnitude larger, respectively, than the flux of energy represented in geochemical sources. Our conception of what an inhabited world "looks like" and our intuition about how to search for life are based in this high-energy context. Energy fluxes on worlds beyond Earth may be better approximated by the million-fold smaller flux provided to Earth's biosphere by geochemical sources. As a result, the nature, abundance, and quality of evidence for life that could be expected on an inhabited extraterrestrial world within our solar system may differ profoundly from that found on Earth. Understanding this potential difference in quantitative terms provides important context for the formulation of life detection strategies. The influence of energy flux on biosignatures can be evaluated through reference to the two basic purposes into which life partitions energy flux: (1) Life expends energy to sustain existing biomass in a metabolic steady state (metabolically functional but non-growing). The formal representation of this relationship in the traditional microbiology literature equates biomass directly with energy flux. The direct implication is that worlds having lower energy flux will have correspondingly lower potential to support biomass. Life detection strategies that directly target extant organisms should therefore be prepared to encounter average biomass densities that may be many orders of magnitude smaller than those found in most of Earth's surface environments (2) Life expends energy to synthesize new biomass. An end-member case in which new biomass is created at the energy-limited rate and the corresponding cells are immediately destroyed (so that the energy partitioned to cell maintenance is minimized) establishes an upper bound on the rate at which biological material can enter a bulk global pool. For a specified bulk concentration [i] of any particular biological compound, i, or for biologically produced matter overall, this synthesis rate, R (sub i), defines a characteristic time scale tau (sub i) equals [i] divided by R (sub i). tau (sub i) can be thought of as (a) the minimum time required for biosynthesis to yield a specific bulk concentration (e.g., a detection threshold) of i, and (b) the average residence time of i within a bulk pool when [i] is held in steady state through a balance between biosynthesis and attrition by physical, chemical, or biological consumption. tau (sub i) becomes an important quantity in considering the potential utility of enantiomeric excess (as a product of homochiral biosynthesis) as a biosignature. Spontaneous racemization of amino acids acts to "erase" the signature of homochiral synthesis over time scales that may range from hundreds to hundreds of thousands of years, depending on temperature. For environments in which low energy flux translates to low rates of biosynthesis, including the synthesis of homochiral amino acids, amino acid residence times in pools having detectable concentrations may compare to or significantly exceed the time scale for racemization. This and similar consequences of long residence times should be considered in the formulation of life detection strategies based on detection of biologically-produced species. Fluxes of biologically-useful energy on potentially habitable worlds within our solar system are, at present, not well constrained. Improving such constraint has the potential to inform priorities in the formulation and targeting of search-for-life strategies, based on the implications of energy flux for the abundance and quality of biosignatures overall, and in specific categories.

ENERGY FLUX↗

DOE OSTI2024

The Design and Evaluation of Zero Trust Architecture for Electric Vehicle Charging Infrastructure: EVs @ Scale Series on EV Charging Station Cybersecurity

Implementing a zero trust architecture can significantly bolster the security of electric vehicle (EV) charging infrastructure. EV charging infrastructure includes numerous networked interfaces, each of which can present potential vulnerabilities. When these vulnerabilities are exploited, they can compromise the entire system, leading to severe operational and security risks. Zero trust is a security model that operates on the principle of "never trust, always verify," which helps manage the attack surface and limit the scope of any potential compromises. Fundamentally, this model ensures that no entity, whether inside or outside the network, is trusted by default. The design principles of zero trust include continuous verification, strict deny-by-default access controls, and micro-segmentation. Continuous verification ensures that every request is thoroughly checked, regardless of its origin. Strict access controls enforce the principle of least privilege, allowing users and devices only the minimum necessary access to perform their functions. Micro-segmentation involves dividing the network into smaller, isolated segments to prevent lateral movement in case of a breach. In the context of EV charging infrastructure, zero trust can be implemented through various strategies. For example, multi-factor authentication (MFA) can be required for engineers to access the management interfaces and control systems of charging stations. Real-time monitoring and analysis of network traffic can help detect and respond to anomalies. Systems that do not need to communicate with each other can be micro-segmented to enhance security. All communications should adhere to predefined policies to be permitted. Additionally, encrypting communications can protect sensitive information exchanged between chargers and management systems. This paper presents a zero trust architecture specifically designed for EV charging infrastructure. Implementing zero trust not only mitigates risks but also builds a resilient infrastructure capable of withstanding and quickly recovering from cyber threats. The architecture addresses six defined security objectives. A comprehensive test plan is developed to assess the architecture against these objectives, and the results of the evaluation are reported. This approach is essential for maintaining the reliability and integrity of EV charging services in an increasingly interconnected and vulnerable digital landscape. This is the first in a planned series of papers exploring the implementation of zero trust in EV charging infrastructure. Each paper will delve into different aspects and applications of zero trust, highlighting how various work processes and requirements can lead to distinct architectural designs. These architectures will be tailored to address specific security challenges and operational needs within the EV charging ecosystem, ensuring a robust and adaptable security framework.

33 ADVANCED PROPULSION SYSTEMS↗

NASA NTRS1994

Emission-line studies of young stars. 4: The optical forbidden lines

Optical forbidden line strengths and profiles are discussed for a sample of 30 T Tauri stars and 12 Herbig Ae-Be stars. Transitions of (C I), (N II), (O I), (O II), (S II), (Ca II), (Cr II), (Fe II), and (Ni II) are detected. Profile variability occurred in DG Tau and probably other sources. The ensemble profiles can be divided into four generic components that may represent distinct emitting regions; (1) narrow rest-velocity lines, (2) 'low'-velocity lines (peaking at less than or approximately +/- 50 km s(exp -1)), (3) 'high'-velocity (usually greater than or approximately +/- 100 km s(exp -1)) blueshifted peaks or wings, and (4) high-velocity redshifted peaks. Among T Tauri stars, the rest-velocity lines appear most often in sources with weak and narrow permitted lines, such as the Ca II triplet. The low- and high-velocity blueshifted components usually appear together in sources with strong and broad Ca II triplet lines. If the velocity-shifted lines form in jets, the smallest (full) opening angles required by the profiles are less than or approximately 20 deg for the narrow, blueshifted (Ca II) lines of DG Tau and HL Tau. Other lines in DG Tau are much broader, implying larger opening angles or greater velocity dispersions. The variability in DG Tau also implies significant changes in the collimation or velocity coherence on timescales of a few years. RW Aur and AS 353A have blue- and redshifted line peaks that could form in oppositely directed jets. The strong (S II) lambda 6716 and lambda 6731 lines in RW Aur are exclusively redshifted and require opening angles less than or approximately 60 deg. Measurements of different profiles in the same spectrum show that the physical conditions change with the line-of-sight velocities. The most persistent trends are for more (N II) and (O II) and less (O I) lambda 5577 flux at high velocities. Constraints on the physical conditions are derived by modeling the emission lines via multilevel ions in 'coronal ionization equilibrium.' A single temperature and density cannot fully describe the line spectra in any velocity interval. Temperatures in the (O I) region are 9000 less than or approximately T(sub e) less than 14,000 K, and the ionization fraction (of H) is less than 35%. The densities derived from (O I) include n(sub e) less than or approximately 5 x 10(exp 5) to approximately 10(exp 7) cm(exp -3), but n(sub e) greater than or approximately 10(exp 6) cm(exp -3) obtains only at low velocities. In the (S II) regions the densities are lower, 10(exp 3) less than or approximately n(sub e) less than or approximately 7 x 10(exp 4) cm(exp -3), and the temperatures are probably higher, T(sub e) greater than or approximately 13,000 K. At high velocities (only) there is additional hot gas that produces (N II) and (O II), possibly most of the (S II), and little (O I). This region is characterized by T(sub e) greater than or approximately 15,000 K, n(sub e) less than or approximately 10(exp 5) cm(exp -3), and an ionization fraction greater than or approximately 50%. When combined with the spatially segregated emitting regions observed by others by spectral imaging, these results suggest decreasing n(sub e) and increasing T(sub e) away from the star in at least the high velocity gas.

Hamann, Fred↗

DOE OSTI2021

A consistent, Bayesian, approach to the cross section probability distribution in the unresolved resonance region

The cross sections of neutron-induced reactions can be divided into three energy ranges: the resolved resonance region (RRR), the unresolved resonance region (URR), and the fast region. In general, the cross sections in the URR show significant fluctuations that cannot be predicted and cannot be experimentally resolved, thus, it is commonly assumed that the cross section at a specific energy is given by a probability distribution function (PDF) over a range of values that can span several orders of magnitude. The current methodology used to describe such behavior is to construct the PDF by stochastically generating resonance ladders and numerically measuring the PDF. The resonance ladders are sampled using known resonance statistical properties and average resonance widths and spacings extrapolated from the RRR. Although this is a standard and widely used technique, it is computationally very expensive, therefore, an alternative, analytical, approach would be preferable due to the considerable speed up of the computational time in real life applications. Moreover, the current methodology does not take into account existing experimental data, such for total and capture cross sections, that are available for many nuclei. Finally, this approach was developed to be used in reactor-scale applications and it is not suited for use in single-event applications. In this work we will rethink the entire approach to the PDF construction using a Bayesian mindset. This will allow us to provide a different definition of the PDF that allows a much faster calculation of the higher-temperature PDFs and a proper combination of theoretical and experimental PDFs following the probability theory. We will also show that our definition is well suited for single-event applications and we will make an explicit connection between our method and the standard approach. We do this by showing that the central limit theorem applies and our method leads to the same PDF obtained with the standard methodology, for a large number of events per history.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2017

Facilitating the Use of Epidemiological Models for Infectious Disease Surveillance

1. To develop a comprehensive model characterization framework to describe epidemiological models in an operational context. 2. To apply the framework to characterize “operational” models for specific infectious diseases and provide a web-based directory, the biosurveillance analytics resource directory (BARD) to the global infectious disease surveillance community. Epidemiological modeling for infectious disease is useful for disease management and routine implementation needs to be facilitated through better description of models in an operational context. A standardized model characterization process that allows selection or making manual comparisons of available models and their results is currently lacking. Los Alamos National Laboratory(LANL) has developed a comprehensive framework that can be used to characterize an infectious disease model in an operational context. We offer this framework and an associated database to stakeholders of the infectious disease modeling field as a tool for standardizing model description and facilitating the use of epidemiological models. Such a framework could help the understanding of diverse models by various stakeholders with different preconceptions, backgrounds, expertise ,and needs, and can foster greater use of epidemiological models as tools in infectious disease surveillance. We define, “operational” as the application of an epidemiological model to a real-world event for decision support and can be used by experts and non-experts alike. The term “model” covers three major types, risk mapping, disease dynamics and anomaly detection. To develop a framework for characterizing epidemiological models we collected information via a three-step process: a literature search of model characteristics, a review of current operational infectious disease epidemiological models, and subject matter expert (SME)panel consultation. We limited selection of operational models to five infectious diseases: influenza, malaria, dengue, cholera and foot-and-mouth disease (FMD). These diseases capture a variety of transmission modes, represent high or potentially high epidemic or endemic burden, and are well represented in the literature. We also developed working criteria for what attributes can be used to comprehensively describe an operational model including a model’s documentation, accessibility, and sustainability. To apply the model characterization framework, we built the BARD, which is publicly available (http://brd.bsvgateway.org).A document was also developed to describe the usability requirements for the BARD; potential users (and non-users) and use cases are formally described to explain the scope of use. 1. Framework for model characterization. The framework is divided into six major components (Figure 1):Model Purpose, Model Objective, Model Scope, Biosurveillance(BSV) goals, Conceptual Model and Model Utility; each of which has several sub-categories for characterizing each aspect of a model.2. Application to model characterization Models for five infectious diseases—cholera, malaria, influenza, FMD and dengue were characterized using the framework and are included in the BARD database. Our framework characterized disparate models in a streamlined fashion. Model information could be binned into the same categories, allowing easy manual comparison and understanding of the models.3. Development of the BARD. Our model characterization framework was implemented into an actionable tool which provides specific information about a model that has been systematically categorized. It allows manual category-to-category comparison of multiple models for a single disease and while the tool does not rank models it provides model information in a format that allows a user to make a ranking or an assessment of the utility of the model. With the model characterization framework we hope to encourage model developers to start describing the many features of their models using a common format. We illustrate the application of the framework through the development of the BARD which is a scientific and non-biased tool for selecting an appropriate epidemiological model for infectious disease surveillance. Epidemiological models are not necessarily being developed with decision makers in mind. This gap between model developers and decision makers needs to be narrowed before modeling becomes routinely implemented in decision making. The characterization framework and the tool developed (BARD) area first step towards addressing this gap.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1999

Experimental Investigation of Pool Boiling Heat Transfer Enhancement in Microgravity in the Presence of Electric Fields

In boiling high heat fluxes are possible driven by relatively small temperature differences, which make its use increasingly attractive in aerospace applications. The objective of the research is to develop ways to overcome specific problems associated with boiling in the low gravity environment by substituting the buoyancy force with the electric force to enhance bubble removal from the heated surface. Previous studies indicate that in terrestrial applications nucleate boiling heat transfer can be increased by a factor of 50, as compared to values obtained for the same system without electric fields. The goal of our research is to experimentally explore the mechanisms responsible for EHD heat transfer enhancement in boiling in low gravity conditions, by visualizing the temperature distributions in the vicinity of the heated surface and around the bubble during boiling using real-time holographic interferometry (HI) combined with high-speed cinematography. In the first phase of the project the influence of the electric field on a single bubble is investigated. Pool boiling is simulated by injecting a single bubble through a nozzle into the subcooled liquid or into the thermal boundary layer developed along the flat heater surface. Since the exact location of bubble formation is known, the optical equipment can be aligned and focused accurately, which is an essential requirement for precision measurements of bubble shape, size and deformation, as well as the visualization of temperature fields by HI. The size of the bubble and the frequency of bubble departure can be controlled by suitable selection of nozzle diameter and mass flow rate of vapor. In this approach effects due to the presence of the electric field can be separated from effects caused by the temperature gradients in the thermal boundary layer. The influence of the thermal boundary layer can be investigated after activating the heater at a later stage of the research. For the visualization experiments a test cell was developed. All four vertical walls of the test cell are transparent, and they allow transillumination with laser light for visualization experiments by HI. The bottom electrode is a copper cylinder, which is electrically grounded. The copper block is heated with a resistive heater and it is equipped with 6 thermocouples that provide reference temperatures for the measurements with HI. The top electrode is a mesh electrode. Bubbles are injected with a syringe into the test cell through the bottom electrode. The working fluids presently used in the interferometric visualization experiments, water and PF 5052, satisfy requirements regarding thermophysical, optical and electrical properties. A 30kV power supply equipped with a voltmeter allows to apply the electric field to the electrodes during the experiments. The magnitude of the applied voltage can be adjusted either manually or through the LabVIEW data acquisition and control system connected to a PC. Temperatures of the heated block are recorded using type-T thermocouples, whose output is read by a data acquisition system. Images of the bubbles are recorded with 35mm photographic and 16mm high-speed cameras, scanned and analyzed using various software packages. Visualized temperature fields HI allows the visualization of temperature fields in the vicinity of bubbles during boiling in the form of fringes. Typical visualized temperature distributions around the air bubbles injected into the thermal boundary layer in PF5052 are shown. The temperature of the heated surface is 35 C. The temperature difference for a pair of fringes is approximately 0.05 C. The heat flux applied to the bottom surface is moderate, and the fringe patterns are regular. In the image a bubble penetrating the thermal boundary layer is visible. Because of the axial symmetry of the problem, simplified reconstruction techniques can be applied to recover the temperature field. The thermal plume developing above the heated surface for more intensive heating is shown. The temperature distribution in the liquid is clearly 3D, and tomographic techniques have to be applied to recover the temperature distribution in such a physical situation. A sequence of interferometric images showing the temperature distribution around the rising bubble, recorded with a high-speed camera is shown. Again, the temperature distribution is 3D, and a more complex approach to the evaluation, the tomographic reconstruction has to be taken. Measurement of the temperature distribution from the fringe pattern temperature distributions that yield important information regarding heat transfer are determined. Two algorithms that allow the quantitative evaluation of interferometric fringe patterns and the reconstruction of temperature fields during boiling have been developed at the Heat Transfer Laboratory of the Johns Hopkins University. In the first algorithm the bubble is assumed to be axially symmetrical, which significantly reduces the computational effort for quantifying the temperature distribution around the bubble. For this purpose the thermal boundary layer around the bubble is divided into equidistant concentric shells, and the refractive index is assumed to be constant in each of the shells. Since large temperature gradients are expected in the vicinity of the bubble during boiling, the deflection of the light beam cannot be neglected in boiling experiments. Since the exit angle of the light beam is known, this allows to account for the deflections and phase shifts outside the boundary layer (in the bulk fluid and in the windows of the test cell). Three dimensional temperature distributions in the vicinity of the bubble are reconstructed using tomographic techniques. In tomography, the measurement volume is sliced into 2D planes. In the present study these planes are parallel to the heated surface. The objective is to determine the values of the field parameter of interest in form of the field function in these 2D planes. The field parameter is the change of the refractive index of the liquid in the measurement volume caused by temperature changes. By superimposing data for many 2D planes recorded at the same time instant, the 3D temperature distribution in the measurement volume is recovered.

Herman, Cila↗

NASA NTRS2011

An Overview of NASA's Integrated Design and Engineering Analysis (IDEA) Environment

Historically, the design of subsonic and supersonic aircraft has been divided into separate technical disciplines (such as propulsion, aerodynamics and structures), each of which performs design and analysis in relative isolation from others. This is possible, in most cases, either because the amount of interdisciplinary coupling is minimal, or because the interactions can be treated as linear. The design of hypersonic airbreathing vehicles, like NASA's X-43, is quite the opposite. Such systems are dominated by strong non-linear interactions between disciplines. The design of these systems demands that a multi-disciplinary approach be taken. Furthermore, increased analytical fidelity at the conceptual design phase is highly desirable, as many of the non-linearities are not captured by lower fidelity tools. Only when these systems are designed from a true multi-disciplinary perspective, can the real performance benefits be achieved and complete vehicle systems be fielded. Toward this end, the Vehicle Analysis Branch at NASA Langley Research Center has been developing the Integrated Design and Engineering Analysis (IDEA) Environment. IDEA is a collaborative environment for parametrically modeling conceptual and preliminary designs for launch vehicle and high speed atmospheric flight configurations using the Adaptive Modeling Language (AML) as the underlying framework. The environment integrates geometry, packaging, propulsion, trajectory, aerodynamics, aerothermodynamics, engine and airframe subsystem design, thermal and structural analysis, and vehicle closure into a generative, parametric, unified computational model where data is shared seamlessly between the different disciplines. Plans are also in place to incorporate life cycle analysis tools into the environment which will estimate vehicle operability, reliability and cost. IDEA is currently being funded by NASA?s Hypersonics Project, a part of the Fundamental Aeronautics Program within the Aeronautics Research Mission Directorate. The environment is currently focused around a two-stage-to-orbit configuration with a turbine-based combined cycle (TBCC) first stage and a reusable rocket second stage. IDEA will be rolled out in generations, with each successive generation providing a significant increase in capability, either through increased analytic fidelity, expansion of vehicle classes considered, or by the inclusion of advanced modeling techniques. This paper provides the motivation behind the current effort, an overview of the development of the IDEA environment (including the contents and capabilities to be included in Generation 1 and Generation 2), and a description of the current status and detail of future plans.

Robinson, Jeffrey S.↗

DOE OSTI2017

Fundamentals of Condensed Matter Physics Marvin L. Cohen and Steven G. Louie

This graduate level textbook on Condensed Matter Physics is written lucidly by two leading luminaries in this field. The volume draws its material from the graduate course in condensed matter physics that has been offered by the authors for several decades at the University of California, Berkeley. Cohen and Louie have done an admirable job of guiding the reader gradually from elementary concepts to advanced topics. The book is divided into four main parts that have four chapters each. Chapter 1 presents models of solids in terms of interacting atoms, which is appropriate for the ground state, and excitations to describe collective effects. Chapter 2 deals with the properties of electrons in crystalline materials. The authors introduce the Born-Oppenheimer approximation and then proceed to the periodic potential approximation. Chapter 3 discusses energy bands in materials and covers concepts from the free electron model to the tight binding model and periodic boundary conditions. Chapter 4 starts with fixed atomic cores and introduces lattice vibrations, phonons, and the concept of density of states. By the end of this part, the student should have a basic understanding of electrons and phonons in materials. Part II presents electron dynamics and the response of materials to external probes. Chapter 5 covers the effective Hamiltonian approximation and the motion of the electron under a perturbation, such as an external field. The discussion moves to many-electron interactions and the exchange-correlation energy in Chapter 6, the widely-used Density Functional Theory (DFT) in chapter 7, and the dielectric response function in Chapter 8. The next two parts of the book cover advanced topics. Part III begins with a discussion of the response of materials to photons in Chapter 9. Chapter 10 goes into the details of electron-phonon interactions in different materials and introduces the polaron. Chapter 11 presents electron dynamics in a magnetic field and Chapter 12 discusses electrical and thermal transport in materials. Part IV takes the reader further into many body effects, superconductivity, and nanoscale materials. The authors introduce Feynman diagrams and many-body perturbation theory in Chapter 13, theories of superconductivity in Chapter 14, magnetism in Chapter 15, and low dimensional systems in Chapter 16. The first two parts are required reading for the beginner planning to perform DFT calculations. The advanced student interested in conducting research in condensed matter physics will benefit from continuing on to the last two parts. There is a set of problems at the end of each part. The narrative is aided by equations and detailed figures. References at the end of the book direct the reader to relevant books and review articles for each chapter. The inside covers include a periodic table and a useful list of fundamental physical constants. The authors present the underlying mathematics elegantly, which makes the textbook quite readable for those with a good mathematical background. Students lacking a firm footing in math will find the terrain rough after Chapter 1. This field has seen many good undergraduate textbooks including those by Kittel and by Ashcroft and Mermin. This volume fills the need for a rigorous graduate level textbook, and is a required addition to the bookshelf of every condensed matter physicist. Cohen and Louie have brought refreshing clarity to a challenging subject and made it eminently accessible to the motivated student.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2019

Cs-137 Body-burdens in Wild Pigs on the Department of Energy's Savannah River Site: Determination of an Appropriate Framework for Spatiotemporal Analysis of Hunter Dose Tracking System Data

The Savannah River Site (SRS) is a U.S. Department of Energy facility with an area of 310 square miles. Since its opening in 1951, SRS has produced and processed nuclear materials for national defense and industrial purposes. These activities have resulted in the release of radiocesium (Cs-137) which has led to contamination of associated streams, floodplains, man-made reservoirs, soil and vegetation. Contamination is spread to wild pigs as they move about SRS, consume vegetation, drink water and wallow in soil. SRS provides annual deer hunts open to the public in which wild pigs may be taken as well. The Cs-137 concentration of each animal is determined before it is released to the hunter for consumption and a cumulative dose is retained for each individual hunter. These tasks are accomplished using the Hunter Dose Tracking System (HDTS). The data collected on SRS by the HDTS for White-Tailed deer has been extensively analyzed. A repeatable method for calculating the background radiocesium body burden for the SRS deer population was established by Gaines and Novak (2009,2016). This model was framed using Cs-137 contamination obtained via gamma over-flight data, the decay rate of radiocesium and the historical behavior observed in the SRS deer population. These parameters were used in a spatiotemporal analysis to estimate background Cs-137 body burdens (from global fallout, etc.), as well as the effect SRS activities may have. The purpose of this project is to modify the framework to represent the behavior of the wild pig population on SRS such that HDTS data collected may be analyzed similarly. The 2018 data were collected via the HDTS in the field to screen wild pigs dispatched by Forest Service - Savannah River. These data will provide insight for an analysis of seasonal variance in the body burdens of the SRS wild pig population. The HDTS quickly provided Cs-137 concentrations in the field, as well as automatically updated the electronic database. Space: Historically, the average home range of wild pigs on the Savannah River Site has been reported at approximately six square miles. This range was confirmed accurate as of July 2018, with the observed minimum and maximum at one and ten square miles respectively. This range was geoprocessed using ArcMap to project buffers based on 1991 gamma over-flight data. The buffers were constructed to originate at all Cs-137 contamination point coordinates and protrude outward in all directions for 1.5, 2.45 and 6 square miles. While the 2.45 mile buffer is theoretically representative of potential hog exposure based on the average home range, it does not consider the possibility of home ranges not uniform in shape. Therefore, the 1.5 and 6 mile buffers were also constructed as visual aids for determining the appropriate hunt compartments for background level considerations. The SRS is divided into numbered hunt compartments of varying size and shape. These compartment numbers are used during sponsored hunts to record the location in which a wild pig was dispatched. Time: The Gaines report determined deer analysis should be carried out at the half-decade level for the best representation of true body burdens. This is due to the impact natural factors, such as drought or flood, can have on the bioaccumulation of Cs-137. To determine the appropriate time scale for hog analysis, all available data was totaled and averaged at 1, 2, 3 and 7 year increments. The annual averages and running average are provided. The appropriate compartments for the determination of the SRS hog population background level should be: 2, 3, 4, 9, 10, 14 and 15. Hogs dispatched within other compartments are most likely exposed to Cs-137 contamination and should not be included in the background level. The current dataset is small, but was still able to demonstrate the relationship between where hogs on SRS and the Cs-137 body burdens they incur. Geoprocessing the data using the average home range of wild pigs highlighted the importance of collecting additional samples from areas surrounding SRS if an accurate background is to be determined. When the average concentration of Cs-137 in pigs from background compartments is compared against that of pigs from contaminated compartments, a significant difference can be seen. A Table provides a summary of the analyzed data from the 15 years processed thus far.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2021

Geometry versus growth: Internal consistency of the flat ΛCDM model with KiDS-1000

Here, we carry out a multi-probe self-consistency test of the flat Lambda Cold Dark Matter (ΛCDM) model with the aim of exploring potential causes of the reported tensions between high- and low-redshift cosmological observations. We divide the model into two theory regimes determined by the smooth background (geometry) and the evolution of matter density fluctuations (growth), each governed by an independent set of ΛCDM cosmological parameters. This extended model is constrained by a combination of weak gravitational lensing measurements from the Kilo-Degree Survey, galaxy clustering signatures extracted from Sloan Digital Sky Survey campaigns and the Six-Degree Field Galaxy Survey, and the angular baryon acoustic scale and the primordial scalar fluctuation power spectrum measured in Planck cosmic microwave background (CMB) data. For both the weak lensing data set individually and the combined probes, we find strong consistency between the geometry and growth parameters, as well as with the posterior of standard ΛCDM analysis. In the non-split analysis, for which one single set of parameters was used, tension in the amplitude of matter density fluctuations as measured by the parameter S 8 persists at around 3σ, with a 1.5% constraint of S 8 = 0.776$_{–0.008}^{+0.016}$ for the combined probes. We also observe a less significant preference (at least 2σ) for higher values of the Hubble constant, H 0 = 70.5$_{–1.5}^{+0.7}$ km s –1 Mpc –1 , as well as for lower values of the total matter density parameter Ω m = 0.289$_{–0.005}^{+0.007}$ compared to the full Planck analysis. Including the subset of the CMB information in the probe combination enhances these differences rather than alleviate them, which we link to the discrepancy between low and high multipoles in Planck data. Our geometry versus growth analysis does not yet yield clear signs regarding whether the origin of the discrepancies lies in ΛCDM structure growth or expansion history but holds promise as an insightful test for forthcoming, more powerful data.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2010

Packing Optimization of an Intentionally Stratified Sorbent Bed Containing Dissimilar Media Types

The Fire Cartridge is a packed bed air filter with two different and separate layers of media designed to provide respiratory protection from combustion products after a fire event on the International Space Station (ISS). The first layer of media is a carbon monoxide catalyst made from gold nanoparticles dispersed on iron oxide. The second layer of media is universal carbon, commonly used in commercial respirator filters. Each layer must be optimally packed to effectively remove contaminants from the air. Optimal packing is achieved by vibratory agitations. However, if post-packing movement of the media within the cartridge occurs, mixing of the bed layers, air voids, and channeling could cause preferential air flow and allow contaminants to pass. Several iterations of prototype fire cartridges were developed to reduce post-packing movement of the media within each layer (settling), and to prevent mixing of the two media types. Both types of movement of the media contribute to decreased fire cartridge performance. Each iteration of the fire cartridge design was tested to demonstrate mechanical loads required to cause detrimental movement within the bed, and resulting level of functionality of the media beds after movement was detected. In order to optimally pack each layer, vertical, horizontal, and orbital agitations were tested and a final packed bulk density was calculated for each method. Packed bulk density must be calculated for each lot of catalyst to accommodate variations in particle size, shape, and density. In addition, a physical divider sheet between each type of media was added within the fire cartridge design to further inhibit intermixing of the bed layers.

Kidd, Jessica↗

NASA NTRS1993

Large cosmic ray transient decreases observed in the heliosphere at the Earth, Voyagers 1 and 2 and Pioneer 10 from 1987 to 1991

Using the greater than 70-MeV cosmic ray rate data from the Interplanetary Monitoring Platform (IMP), Voyager, and Pioneer spacecraft, we have examined the role of large transient decreases in the overall cosmic ray decrease from 1987 to 1990. At least 12 separate transient decreases with magnitudes greater than 5% were observed at Voyager 2, which moved from a heliocentric radius of approximately 23 to 35 AU during the time period. Many, but not all, of these decreases were observed at IMP 8 near the Earth and at the Pioneer 10 between 42 and 54 AU near the ecliptic plane. However, only five of these decreases were observed at Voyager 1 at north heliographic latitude approximately 30 deg and radial distance between 32 and 46 AU. The period of decreasing cosmic ray intensity from 1987 to 1990 can be divided into two phases. During phase 1, between mid-1987 and early 1989, the average heliospheric current sheet tilt rapidly increased from approximately 8 deg to 65 deg as the solar activity was also increasing rapidly and the intensity at the Earth decreased by about 35%. Four large transient decreases were observed at the Earth during this time; however, transient decreases beyond the Earth were observed only near the ecliptic plane. In phase 2, from early 1989 to mid-1990 when the intensity at the Earth decreased by an additional approximately 60%, five large transient decreases were observed at all latitudes. During this time the current sheet tilt was always greater than 60 deg. Using a procedure to separate these large transient decreases from any possible slower long-term variations, we argue that these decreases appear to be responsible for much of the overall intensity decrease of 20 - 35% observed at all spacecraft near the equatorial plane in phase 1 as well as most of the larger intensity decrease of about 45 - 60% observed at all spacecraft during phase 2. The recovery of the cosmic ray intensity that started at the Earth in early 1990 moved outward in the heliosphere reaching the Voyager and Pioneer spacecraft about 6 months later. This is in contrast to what happened as the intensity began to recover was observed to begin almost simultaneously at all radii less than 25 AU. This difference suggests that during the current cycle, the mechanism responsible for the recovery must have propagated outward from the Sun. This could have been related to the reversal of the solar polar magnetic field which occurred in early 1990. The effects of large transients propagating outward in the heliosphere and becoming ineffective when they reach the heliospheric boundary could not have played a major role in the onset of the recovery in 1990, since in this case the recovery should have been observed first in the outer heliosphere.

Webber, W. R.↗

NASA NTRS2017

A Case for Nebula Scale Mixing Between Non-Carbonaceous and Carbonaceous Chondrite Reservoirs: Testing the Grand Tack Model with Chromium Isotopic Composition of Almahata Sitta Stone 91A

There is an increasing number of Cr-O-Ti isotope studies that show solar system materials are divided into two main populations, one carbonaceous chondrite (CC)-like and the other is non-carbonaceous (NC)-like, with minimal mixing attributed to a gap opened in the protoplanetary disk due to Jupiter's formation. The Grand Tack model suggests there should be large-scale mixing between S- and C-type asteroids, an idea supported by our recent work on chondrule (Delta)17O-ε54Cr isotope systematics. The Almahata Sitta (AhS) meteorite provides a unique opportunity to test the Grand Tack model. The meteorite fell to Earth in October 2008 and has been linked to the asteroid 2008 TC3 which was discovered just prior to the fall of the AhS stones. The AhS meteorite is composed of up to 700 individual pieces with approx.140 of those pieces having some geochemical and/or petrologic studies. Almahata Sitta is an anomalous polymict ureilite with other meteorite components, including enstatite, ordinary, and carbonaceous chondrites with an approximate abundance of 70% ureilites and 30% chondrites. This observation has lead to the suggestion that TC3 2008 was a loosely aggregated rubble pile-like asteroid with the non-ureilite sample clasts within the rubble-pile. Due to the loosely-aggregated nature of AhS, the object disintegrated during atmospheric entry resulting in the weakly held clasts falling predominantly as individual stones in the AhS collection area. However, recent work has identified one sample of AhS, sample 91A, which may represent two different lithologies coexisting within a single stone. The predominate lithology type in 91A appears to be that of a C2 chondrite based on mineralogy but also contains olivine, pyroxene, and albite that have ureilite-like compositions. Previous Cr isotope investigations into AhS stones are sparse and what data is available show nearly uniform isotopic composition similar to that of typical ureilites with negative ε54Cr values.

Sanborn, M. E.↗

NASA NTRS2005

Experiments and Models for Polymeric Microsphere Foams

The current project was performed under the direction of Dr. Byron Pipes as its lead investigator from January 2001 to August 2004. With the permission of the NASA, the project was transferred to Dr. Thein Kyu as the principle investigator for the period of September 2004 - June 2005. There were two major thrust areas in the original proposal; (1) experimental characterization and kinematics of foam structure formation and (2) determination of the mechanical, physical, and thermal properties, although these thrust areas were further sub- divided into 7 tasks. The present project has been directed primarily to elucidate kinematics of micro-foam formation (tasks 1 and 3) and to characterize micro-foam structures, since the control of the micro-structure of these foams is of paramount importance in determining their physical, mechanical and thermal properties. The first thrust area was accomplished in a timely manner; however, the second thrust area of foam properties (tasks 2,4-7) has yet to be completed because the area of kinematics of foam structure formation turned out to be extremely complex and thus consumed more time than what have been anticipated. As will be reported in what follows, the present studies have greatly enhances the in-depth understanding of mechanisms and kinematics of the micro-foam formation from solid powders. However, in order to implement all objectives of the second thrust areas regarding investigations of mechanical, physical, and thermal properties and establishment of the correlation of structure - properties of the foams, the project needs additional time and resources. The technical highlights of the accomplishment are summarized as follows. The present study represents a first approach to understanding the complexities that act together in the powder foaming process to achieve the successful inflation of polyimide microstructures. This type of study is novel as no prior work had dissected the fundamentals that govern the inflation process in this type of systems. The systematic approach to each of the different phenomena (i.e. morphological, diffusive, kinetic and dynamic) brings into context each of them in a way that allows separate understanding and analysis. Of the different phenomena studied, probably the one that gives a higher level of control over the inflation process has been shown to be the morphological aspects of the precursor particles. It is a major contribution of the present work to isolate and identify this phenomenon and highlight the features that with careful control during the synthesis of the precursor material can lead to a highly optimized and specialized final product (neat foam or microstructure). Some of these accomplishments have been presented in various national meetings and some of which are either published in refereed journals or still in various stages of publications. One of the presentations was selected for "Best of ANTEC 2004" Online Presentation Series of the Society of Plastics Engineers (SPE) (September 2004)

Pipes, R. Byrona↗

NASA NTRS1973

On Lidar Sounding of the Atmosphere to Estimate Static and Dynamic Characteristics of Aerosol Inhomogeneities

1. A possible application of intensity fluctuations of a pulse light signal reflected by atmospheric aerosols is analyzed by the correlation method to evaluate static (medium sizes, shape) and dynamic (speed and direction of movement, lifetime) characteristics of aerosol inhomogeneities. The aerosol inhomogeneities are assumed to be expanded, pressed, disintegrated and originated constantly in accordance with random laws, the set of inhomogeneities as a whole traveling together with air masses and having predominant movement in wind direction. It is shown that the characteristics of aerosol inhomogeneities considered can be expressed by the coefficients of the correlation function expansion of the reflected signal fluctuation intensity in Tailor series. 2/ Correlation systems for evaluating static and dynamic characteristics of driving objects can be divided into two types according to the kind and quantity of used information: the systems with coordinates of the information removal "points" to be fixed in space, and the systems with a parallel simultaneous information removal at discrete moments of time. The systems for determination of wind direction considered in are the examples of the first type system. However, the operating information removal for two points is insufficient to estimate completely static and dynamic characteristics of inhomogenities, their quantity ought to be increased up to three of them for two-dimensional problem and up to four of them for three-dimensional problem as it is usually done in the ionospheric studies. The second type systems are used for the investigation of a medium shape and speed of the clouds according to photographs made from satellites. These systems are also used for solution of navigation problems. The use of optical quantum generators with a scanning beam is seen to increase greatly the working information removal in comparison with the first type systems. Nevertheless, scanning rate is not sufficient sometimes in order· to consider a general picture of aerosol inhomogeneities to be stationary. In this connection the use of the systems of second type treatment becomes a matter of essential difficulty. 3. Aerosol inhomogeneities simulation has been carried out on the basis of the digital computer experiments with the aim of estimating static and dynamic characteristics of inhomogeneities by an optical beam in the atmosphere at different scanning procedures. The dependence of determination accuracy of these characteristics on the type of chosen laws of aerosol particle distribution in the atmosphere, the parameters of inhomogeneities geometry, their speed and the law of scanning have been obtained.

Zuev, V. E.↗

NASA NTRS2017

Effects of Acute Exposures to Carbon Dioxide upon Cognitive Functions

Carbon dioxide (CO2) originates from human metabolism and typically remains about 10-fold higher in concentration on the International Space Station (ISS) than at the earth's surface. There have been recurring complaints by crew members of episodes of "mental viscosity" adversely affecting their performance, and there is evidence from the ISS that associates CO2 levels with reports of headaches by crewmembers. Consequently, flight rules have been employed to control CO2 below 3 mm Hg, which is well below the existing Spacecraft Maximum Allowable Concentration (SMAC) of 10 mm Hg for 24-hour exposures, and 5.3 mm Hg for exposures of 7 to 180 days. Headaches, while sometime debilitating themselves, are also symptoms that can provide evidence that physiological defense mechanisms have been breached, and there is evidence that CO2 has effects at levels below the threshold for headaches. This concern appears to be substantiated in reports that CO2 at concentrations below 2 mm Hg substantially reduced some cognitive functions that are associated with the ability to make complex decisions in conditions that are characterized by volatility, uncertainty, complexity, ambiguity, and delayed feedback. These are conditions that could be encountered by crews in off-nominal situations or during the first missions beyond low earth orbit. Therefore, we set out to determine if decision-making under volatile, uncertain, confusing and ambiguous circumstances, where feedback is delayed or absent, is correlated with low levels of CO2 during acute exposures (several hours) in crew-like subjects and to determine if additional cognitive domains are sensitive to concentrations of CO2 at, or below, current ISS levels by using a test battery that is currently available onboard ISS. We enrolled 22 volunteers (8 females, 14 males) between the ages of 30-55 (38.8 +/- 7.0) years whose training and professional experience reflect that of the astronaut corps. Subjects were divided among 4 study groups. Exposures occurred in the 20-foot environmental chamber in Bldg 7 at JSC. Each group was exposed to each of four concentrations of CO2 (600, 1200, 2500, and 5000 ppm), and the exposure order was balanced across the groups which were randomly assigned to 4 possible complete exposure orders. Study participants and investigators were blinded to the exposure order. Each exposure lasted 4 hours and occurred on the same day each week for 4 consecutive weeks. Testing included the Strategic Management Simulation (SMS) methodology and the Cognition battery of psychometric measures that are being utilized aboard the ISS. Subjects participated in a familiarization session one week prior to the start of exposure sessions. Each exposure session followed the schedule noted in Figure 1 below. On the morning of each exposure, subjects completed the Cognition battery in a conference room at indoor ambient CO2 levels. Subjects then entered the exposure chamber and were acclimated to the environment for 15 minutes before completing another round of the Cognition battery. After a 20 minute rest, subject completed the approx.80 minute SMS. Subjects were allowed another 20 minute rest period before completing the final Cognition test battery in the exposure chamber. Subjects then returned to the conference room where they completed one final round of Cognition. Data from both testing methodologies are currently being analyzed and will be presented. Study investigators have not yet been unblinded to the coded exposure concentrations, but preliminary results suggest differences.

Ryder, V. E.↗

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↗