Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “strain rate”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 631 records · Page 35

Symposium MT02: Statistical Mechanics-Based Computational Tools for the Study of Phase Transformation in Complex Materials (Final Report)

Symposium MT02 brought together a diverse and interdisciplinary community of scientists specializing in Statistical Mechanics-based computational modeling to investigate phase transformations in materials exhibiting complex disordered structures. As the demand for materials with extreme performance metrics grows—from aerospace components to next-generation optical fibers—the ability to predict microstructural evolution under non-equilibrium conditions has become paramount. The primary goal of this symposium was to identify, evaluate, and discuss advanced computational tools capable of designing precise manufacturing conditions to tailor material properties efficiently. By fostering a dialogue between computational theorists and experimentalists, the symposium sought to establish new protocols for predicting how processing history—such as cooling rates or strain paths—dictates the final microstructure.

36 MATERIALS SCIENCE↗

The Role of Helium on Ejecta Production in Copper

The effect of helium (He) concentration on ejecta production in OFHC-Copper was investigated using Richtmyer–Meshkov Instability (RMI) experiments. The experiments involved complex samples with periodic surface perturbations machined onto the surface. Each of the four target was implanted with a unique helium concentration that varied from 0 to 4000 appm. The perturbation’s wavelengths were λ ≈ 65 μ m, and their amplitudes h 0 were varied to determine the wavenumber ( 2 π / λ ) amplitude product k h 0 at which ejecta production beganfor Cu with and without He. The velocity and mass of the ejecta produced was quantified using Photon Doppler Velocimetry (PDV) and Lithium-Niobate (LN) pins, respectively. Our results show that there was an increase of 30% in the velocity at which the ejecta cloud was traveling in Copper with 4000 appm as compared to its unimplanted counterpart. Our work also shows that there was a finer cloud of ejecta particles that was not detected by the PDV probes but was detected by the early arrival of a “signal” at the LN pins. While the LN pins were not able to successfully quantify the mass produced due to it being in the solid state, they did provide information on timing. Our results show that ejecta was produced for a longer time in the 4000 appm copper.

36 MATERIALS SCIENCE↗

Considerations for Thermal Annealing of Zr-Alloy Cladding During Dry Storage

In January 2025, the U.S. Department of Energy (DOE) sponsored a technical workshop with fuel vendors that provided a forum for increased collaboration on backend fuel cycle research. The workshop reviewed observations of thermal annealing of irradiated zirconium alloys in simulated dry storage conditions and discussed implications to the industry’s interest in reducing wet storage to optimize operational flexibility and the expected trend of increasing discharge burnup (and decay heats), which could lead to an increase in dry storage temperatures. One approach to accommodate increased decay heats in dry storage is to increase the regulatory temperature limit during dry storage, which is now 400°C. However, this could lead to thermal annealing of the cladding, which will tend to increase creep rates, creep strains, and ductility while decreasing yield strength. These considerations are compounded by bonding between the pellet and cladding at high burnup, which results in the pellet carrying more load during fuel rod deformations, especially in bending, which is the dominant deformation phenomenon in canister drop analyses. Because the U.S. Nuclear Regulatory Commission (NRC) recommends yield strength as a primary failure criterion for accident analyses in storage and transportation, particularly for canister drop scenarios, the purpose of the workshop was to develop options for alternative failure criterion that could replace yield strength in canister drop analyses. Topics discussed were (1) the regulatory framework, (2) the effect of thermal annealing on microhardness and tensile properties of cladding materials, (3) the effect of thermal annealing on fuel rod failure during bending, bending fatigue, and pinch loading, and (4) an alternative failure criterion for canister drop analyses. This paper provides a summary of the key discussions and outcomes of the workshop

Cantonwine, Paul [ORNL] (ORCID:0009000522247033)↗

Infusible thermoplastic resin based sandwich structures for wind blade applications and the influence of scrim on facesheet to core interface debonding

Infusible thermoplastic Elium® family of resins from Arkema have garnered much attention in recent years as a possible replacement for thermoset resins in laminate and sandwich composite manufacturing for wind blade applications due to its ease of recyclability and the ability to utilize existing manufacturing processes without imposing complicated variations. However, physical and mechanical properties of the proposed Elium® based thermoplastic composites must be comparable to existing epoxy (thermoset) based composites using manufacturing processes relevant for large wind turbine blades. A 13-meter-long demonstration blade was manufactured for that purpose and sandwich samples were obtained from that project for a detailed study. This paper details three-point flexural properties of unidirectional E-glass fiber reinforced acrylic and epoxy based sandwich panels with identical balsa wood core materials. In addition, to evaluate the relative merit considering debond failure mode, the interfacial critical strain energy release rate, predominantly in mode-1, was compared via single cantilever beam testing. In sandwich composites constructed with balsa wood core material, resin uptake by the balsa core is traditionally impeded via the insertion of a scrim material at the facesheet to core interface. Results revealed that inclusion of scrim mesh layer at the facesheet to core interface reduced flexural properties and strain energy release rates in panels infused with acrylic resin but did not significantly reduce these properties in epoxy infused facesheets.

Engineering↗

A macro-micro approach for identifying crystal plasticity parameters for necking and failure in nickel-based alloy haynes 282

Here, this work develops a two-scales macro-micro approach to address the challenge in calibrating crystal plasticity microstructural models when samples undergo necking prior to fracture. The crystal plasticity models are crucial for predicting the materials’ plastic deformation and failure at the microstructure level, identifying the materials’ intrinsic properties as well as investigating the microstructure-properties relationships. However, after necking occurs, the experimentally measured stress-strain curves fail to reflect the materials ‘true’ stress-strain behavior and cannot be directly fitted into crystal plasticity models. The proposed macro-micro approach employs a top-down strategy to address this challenge, which has been studied with experimental tests on precipitation-strengthened Ni-based superalloy Haynes® 282®. In this approach, a macro rate-dependent anisotropic plasticity model with Voce-type hardening and Rice-Tracey damage law is first utilized to model the deformation and failure of the tensile bar, and calibrated by matching the stress-strain curves, necking strain, and reduction of area. Especially, to match the testing results under different applied strain rates, the rate-sensitivity parameter m and saturation stress in the elasticity model are modified to incorporate dependence on the local strain rate. Then, the ‘true’ stress-strain behaviors are extracted from the necking zone of the macro-model, which are used to calibrate a micro-model with explicit microstructures and governed by an extended crystal plasticity law. The consistency between the micro-model and macro-model are enforced during calibration. The calibration outcomes from the crystal plasticity model elucidate the materials intrinsic properties for slip, hardening, and failure, which is vital for further investigations into the microstructure-properties relationship and for accurate prediction of the material behavior under various test and service conditions.

36 MATERIALS SCIENCE↗

Error Analysis on Numerical Integration Algorithms in a Hypoelasticity Framework

This report determines local truncation errors for common stress integration algorithms used in explicit finite element codes with hypoelastic material models. The hypoelastic integration algorithms in question utilize an operator splitting procedure in a rotation neutralized configuration, where the stress response is determined from de- coupling the total deformation into rotational and strain dependent components. This document analyzes the error in evolving the stress given a one-step time increment Δt and compares the errors associated with both the rotational and strain components of the operator splitting method. A slight modification to a traditional algorithm is proposed and studied, where the rate of deformation is appropriately rotated from the midstep configuration at t n+1/2 to the end step configuration at t n+1 before the constitutive evaluation. The proposed modification either completely eliminates the error associated with the rotation rate or is of the same order of magnitude as the original algorithm for the three test cases considered in this report. These cases consist of an unaxial stretch with a constant true strain rate with a rigid body rotation, an uniaxial stretch with a constant engineering strain rate with a rigid body rotation, and a simple shear deformation. All three cases are compared to a closed form solution, and in almost every test case the alternative algorithm yields the most accurate one-step local truncation error.

97 MATHEMATICS AND COMPUTING↗

Mechanochemical Regulation of Oxidative Addition to a Palladium(0) Bisphosphine Complex

Here, we report the effect of force applied to the biaryl backbone of a bisphosphine ligand on the rate of oxidative addition of bromobenzene to a ligand-coordinated palladium center. Local compressive and tensile forces on the order of 100 pN were generated using a stiff stilbene force probe. A compressive force increases the rate of oxidative addition, whereas a tensile force decreases the rate, relative to that of the parent complex of strain-free ligand. Rates vary by a factor of ~6 across ~340 pN of force applied to the complexes. The crystal structures and DFT calculations support that force-induced perturbation of the geometry of the reactant is negligible. The force-rate relationship observed is mainly attributed to the coupling of force to nuclear motion comprising the reaction coordinate. These observations inform the development of catalysts whose activity can be tuned by an external force that is adjusted within a catalytic cycle.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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

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

well test, sinusoidal, periodic pumping, strain↗

Universal behavior in fragmenting brittle, isotropic solids across material properties

A bonded particle model is used to explore how variations in the material properties of brittle, isotropic solids affect critical behavior in fragmentation. To control material properties, a model is proposed which includes breakable two- and three-body particle interactions to calibrate elastic moduli and mode I and mode II fracture toughnesses. In the quasistatic limit, fragmentation leads to a power-law distribution of grain sizes which is truncated at a maximum grain mass that grows as a nontrivial power of system size. In the high-rate limit, truncation occurs at a mass that decreases as a power of increasing rate. A scaling description is used to characterize this behavior by collapsing the mean-square grain mass across rates and system sizes. Consistent scaling persists across all material properties studied, although there are differences in the evolution of grain size distributions with strain as the initial number of grains at fracture and their subsequent rate of production depend on Poisson's ratio. Importantly, this evolving granular structure is found to induce a unique rheology where the ratio of the shear stress to pressure, an internal friction coefficient, decays approximately as the logarithm of increasing strain rate. The stress ratio also decreases at all rates with increasing strain as fragmentation progresses and depends on elastic properties of the solid.

36 MATERIALS SCIENCE↗

Hot deformation behavior of an industrially cast large grained low density austenitic steel

Hot compression testing was performed on specimens of an Fe–30Mn–9Al–1Si-0.9C-0.5Mo wt.% low density steel which were extracted from an industrial casting with large grain size (320 ± 50 μm). Grain boundaries act as nucleation sites for dynamic recrystallization in the present material and the large grain size results in fewer nucleation sites, which delays dynamic recrystallization to larger strains and/or higher temperatures and causes higher rates of work hardening relative to similar materials but with smaller grain sizes. Here, a relatively high activation energy for hot deformation of 470 ± 90 kJ mol –1 was determined from a Zener-Hollomon analysis of the flow stresses at different temperatures and strain rates, attributed in part to the large grain size.

36 MATERIALS SCIENCE↗

Cyclic strain enhances the early stage mineral nucleation and the modulus of demineralized bone matrix

The adaptive response of bones to mechanical loading is essential for musculoskeletal development. Despite the importance of collagen in bone mineralization, little is known about how cyclic strain influences physicochemical responses of collagen, especially at the early stage of mineralization when the levels of strain are higher than those in mature bones. The findings in this study show that, without any cell-mediated activity, cyclic strain increases nucleation rates of calcium phosphate (CaP) nanocrystals in highly-organized collagen matrices. The cyclic strain enhances the transport of mineralization fluids with nucleation precursors into the matrix, thus forming more CaP nanocrystals and increasing the elastic modulus of the collagen matrix. The results also suggest that the multiscale spatial distribution of nanocrystals in the fibrous collagen network determines tissue-level mechanical properties more critically than the total mineral content. Furthermore, by linking nano- and micro-scale observations with tissue-level mechanical properties, we provide new insights into designing better biomaterials.

36 MATERIALS SCIENCE↗

Inviscid damping of an elliptical vortex subject to an external strain flow

Inviscid spatial Landau damping is studied experimentally for the case of oscillatory motion of a two-dimensional vortex about its elliptical equilibrium in the presence of an applied strain flow. Here, the experiments are performed using electron plasmas in a Penning–Malmberg trap. They exploit the isomorphism between the two-dimensional Euler equations for an ideal fluid and the drift-Poisson equations for the plasma, where plasma density is the analog of vorticity. Perturbed elliptical vortex states are created using E x B strain flows, which are generated by applying voltages to electrodes surrounding the plasma. Measurements of spatial Landau damping (also called critical-layer damping) are in agreement with previous studies in the absence of an applied strain, where the damping is due to a resonance between the local fluid motion and the vortex oscillations. Interestingly, the damping rate does not change significantly over a wide range of applied strain rates. This can be accurately predicted from the initial vorticity profile, even though the resonant frequency is reduced substantially due to the applied strain. For higher amplitude perturbations, nonlinear trapping oscillations also exhibit behavior similar to the strain-free case. In principle, higher-order effects of the applied strain, such as separatrix crossing of peripheral vorticity and interactions with harmonics of the fundamental resonance, are expected to change the damping rate. However, this occurs only for conditions that are not realized in the experiments described here. Vortex-in-cell simulations are used to investigate the possible roles of these effects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Viscometric flow of dense granular materials under controlled pressure and shear stress

This study examines the flow of dense granular materials under external shear stress and pressure using discrete element method simulations. In this method, the material is allowed to strain along all periodic directions and adapt its solid volume fraction in response to an imbalance between the internal state of stress and the external applied stress. By systematically varying the external shear stress and pressure, the steady rheological response is simulated for: (1) rate-independent quasi-static flow; and (2) rate-dependent inertial flow. The simulated flow is viscometric with non-negligible first and second normal stress differences. While both normal stress differences are negative in inertial flows, the first normal stress difference switches from negative to slightly positive, and second normal stress difference tends to zero in quasi-static flows. The first normal stress difference emerges from a lack of coaxiality between a second-rank contact fabric tensor and strain rate tensor in the flow plane, while the second normal stress difference is linked to an excess of contacts in the shear plane compared with the vorticity direction. A general rheological model of second order (in terms of strain rate tensor) is proposed to describe the two types of flow, and the model is calibrated for various values of interparticle friction from simulations on nearly monodisperse spheres. The model incorporates normal stress differences in both regimes of flow and provides a complete viscometric description of steady dense granular flows.

42 ENGINEERING↗

Comparison of wild-type KT2440 and genome-reduced EM42 Pseudomonas putida strains for muconate production from aromatic compounds and glucose

Pseudomonas putida KT2440 is a robust, aromatic catabolic bacterium that has been widely engineered to convert bio-based and waste-based feedstocks to target products. Towards industrial domestication of P. putida KT2440, rational genome reduction has been previously conducted, resulting in P. putida strain EM42, which exhibited characteristics that could be advantageous for production strains. Here, we compared P. putida KT2440-and EM42-derived strains for cis,cis-muconic acid production from an aromatic compound, p-coumarate, and in separate strains, from glucose. To our surprise, the EM42-derived strains did not outperform the KT2440-derived strains in muconate production from either substrate. In bioreactor cultivations, KT2440-and EM42-derived strains produced muconate from p-coumarate at titers of 45 g/L and 37 g/L, respectively, and from glucose at 20 g/L and 13 g/L, respectively. To provide additional insights about the differences in the parent strains, we analyzed growth profiles of KT2440 and EM42 on aromatic compounds as the sole carbon and energy sources. In general, the EM42 strain exhibited reduced growth rates but shorter growth lags than KT2440. We also observed that EM42-derived strains resulted in higher growth rates on glucose compared to KT2440-derived strains, but only at the lowest glucose concentrations tested. Transcriptomics revealed that genome reduction in EM42 had global effects on transcript levels and showed that the EM42-derived strains that produce muconate from glucose exhibit reduced modulation of gene expression in response to changes in glucose concentrations. Overall, our results highlight that additional studies are warranted to understand the effects of genome reduction on microbial metabolism and physiology, especially when intended for use in production strains.

09 BIOMASS FUELS↗

Understanding the transient large amplitude oscillatory shear behavior of yield stress fluids

A full understanding of the sequence of processes exhibited by yield stress fluids under large amplitude oscillatory shearing is developed using multiple experimental and analytical approaches. A novel component rate Lissajous curve, where the rates at which strain is acquired unrecoverably and recoverably are plotted against each other, is introduced and its utility is demonstrated by application to the analytical responses of four simple viscoelastic models. Using the component rate space, yielding and unyielding are identified by changes in the way strain is acquired, from recoverably to unrecoverably and back again. The behaviors are investigated by comparing the experimental results with predictions from the elastic Bingham model that is constructed using the Oldroyd–Prager formalism and the recently proposed continuous model by Kamani, Donley, and Rogers in which yielding is enhanced by rapid acquisition of elastic strain. The physical interpretation gained from the transient large amplitude oscillatory shear (LAOS) data is compared to the results from the analytical sequence of physical processes framework and a novel time-resolved Pipkin space. The component rate figures, therefore, provide an independent test of the interpretations of the sequence of physical processes analysis that can also be applied to other LAOS analysis frameworks. Each of these methods, the component rates, the sequence of physical processes analysis, and the time-resolved Pipkin diagrams, unambigiously identifies the same material physics, showing that yield stress fluids go through a sequence of physical processes that includes elastic deformation, gradual yielding, plastic flow, and gradual unyielding.

Kamani, Krutarth M. (ORCID:0000000338975420)↗

A cavitation and dynamic void growth model for a general class of strain-softening amorphous materials

In this work, we derive a homogenized framework for cavitation and dynamic void growth in a general class of strain-softening materials with particular emphasis on amorphous materials. Analytic solutions are derived by considering an assemblage of spherical shells subject to dynamic hydrostatic tensile pressure. The isochoric response of the matrix material is captured by a free volume theory based amorphous viscoplasticity model. The framework accounts for the coupled rate-dependent, strain-softening behavior typically exhibited in amorphous materials. Validation of the homogenized theory is carried out against numerous direct numerical simulations with excellent agreement. Furthermore, we find that a reduction in cavitation strength and significant acceleration of dynamic void growth rates is induced as a consequence of strain-softening in the amorphous matrix material. Lastly, the model is utilized to interpret some non-intuitive experimental observations of amorphous materials subject to shock compression and subsequent spall failure.

36 MATERIALS SCIENCE↗

Hybrid combustion modeling approach for turbulent jet ignition in natural-gas pre-chamber spark-ignition engines at high EGR

Here, this study presented a hybrid modeling approach for simulating turbulent jet ignition and combustion processes in a natural-gas pre-chamber spark-ignition engine operating under exhaust gas recirculation (EGR) diluted conditions. In-depth analyses of experimental data and simulation results from previous work [Chinnathambi et al., ICEF2021-67836; Kim et al., Fuel 409: 137815, 2026] revealed two key findings: (i) the magnitude of pressure difference between the pre-chamber and main chamber ($∆P_{PC-MC}$) was positively correlated with the combustion duration from the moment of $∆P_{PC-MC}=0$ to the point of 5% mass fraction burned, with larger $∆P_{PC-MC}$ associated with longer duration; and (ii) the turbulent combustion regime in the main chamber transitioned from the broken reaction zone to the corrugated flamelet regime, with the Karlovitz number exceeding 100 immediately after turbulent hot jets were ejected from nozzles, coinciding with observed local extinction events. To accurately simulate the entire combustion process, a hybrid approach was developed under Reynolds-Averaged Navier Stokes framework, combining the G-equation model for pre-chamber combustion with the multi-zone well-stirred reactor approach and a turbulence-chemistry interaction (TCI) submodel for main chamber combustion. The TCI submodel accounted for the attenuation of reaction rates due to turbulent strain and modeled local extinction by suppressing reaction rates under certain flow and flame conditions. When applied to three EGR rate conditions toward the dilution limit, the hybrid modeling approach accurately reproduced experimental data in terms of cylinder pressure, apparent heat release rate, and the observed positive correlation, including the delayed onset of main chamber combustion—a feature not captured by existing combustion models.

computational fluid dynamics simulation↗

A Novel Technique for Distributed Measuring of Strain in the Vadose Zone

Characterizing water content and pressure changes in the vadose zone is important to understanding a variety of geologic processes, ranging from permeability, to evapotranspiration, and aquifer recharge. Changes in water content, or pressure, cause strain in the solid porous medium of the vadose zone, and it is possible that measuring those strains could be used a characterization tool. The Coherence-length-gated Microwave Photonics Interferometry (CMPI) technique measures strain at high resolution along many intervals defined by pairs of reflectors distributed along an optical fiber. This technique has recently been developed at Clemson University, and it has the spatial, and temporal resolution to characterize the strains that are expected to occur with hydrologic changes in the vadose zone. However, the technique has never been used to measure strain in porous media, so its capabilities remain uncertain. The objective of this thesis is to evaluate the ability of using CMPI to measure strain changes in the vadose zone. The research approach consists of conducting laboratory tests using a column filled with sand that was subjected to changes in water content and pressure. An optical fiber with CMPI reflectors was deployed in a high surface area ribbon and used to measure strain along the axis of the column. The column was made from 8-inch, Schedule 40 PVC pipe (75 cm tall, 20 cm inner diameter) and filled with medium-grained sand (K = 2.5x10-6 m/s, porosity = 0.28, Coefficient of Uniformity = 2, van Genuchten (n = 4.34 and  = 0.00039 m-1), Young’s Modulus = 18 – 47 MPa). Five pressure and four temperature sensors attached to probes inserted into the wall. The optical fiber includes five reflectors spaced 10 cm apart on the inside and outside of the column along the vertical axis. Strain measurements are made between pairs of reflectors. The reflectors are created in 250-micron-diameter acrylate-coated single mode Corning SMF28e+ optical fiber using a femtosecond laser. The optical fiber was laminated between two pieces of polyester film creating a large surface area to transfer strain from the porous media to the optical fiber. The experiments were conducted by filling the column from the bottom or infiltrating water from the top. The water was allowed to equilibrate to room temperature prior to each test in order to limit thermoelastic strain. Five injection tests with a rate of 250 ml/min and three infiltration tests at varying rates were conducted, and the results show patterns of strain and pressure are generally similar. Hydrologic conditions define three zones based on the pressure magnitude and distribution. 1.) Ambient Zone where the pressure heads are quasi-static and the pressure gradient is roughly unity (head gradient of zero). This is the uppermost zone and is characterized by negative pressures. 2.) Transition Zone where the pressures increase from ambient to zero, are changing relatively rapidly and the pressure gradient is relatively steep (pressure head gradients of 2). The Transition zone is 10 to 15 cm thick. 3.) Positive Pressure Zone where the pressure is positive, the rate of change is slower than in the transition zone and the gradient is flatter (pressure head gradient 1.1 to 1.2). This is the lowest zone in the column. Injection of water causes the pressure to increase and the three zones to move upward. The Transition zone moves at a rate of approximately 0.0002 m/s +/- .00005, according to analyses of pressure profiles. This is the same as the average velocity of the water calculated as volumetric flux/effective porosity =5.4x10-5 m/s/0.28. Strain signals in the range of 10s  were observed with a noise level of generally less than 0.1  (signal to noise ratio of greater than 100) during injection and drainage. The spatial and temporal distributions of strain caused by injection depend on the location of the moving pressure zones. The locations of the different pressure zones were determined from pressure profiles at different times and these data were transferred to strain time series. This showed that strain in the Ambient pressure zone is either unchanged or slightly compressive, whereas strain in Transition zone is tensile and roughly proportional to the pressure change, and strain ranges from tensile to compressive in the Positive Pressure zone. In general, the strain is variable at the top of the Positive Pressure zone, but it appears to be consistently tensile lower in the zone. The spatial distribution of strain causes three distinct stages in the strain times series measured between pairs of CMPI reflectors. The strain is unchanged or slightly compressive during Stage 1, it increases (positive strain is tensile) approximately linearly at a rate of 0.02 to 0.03 μ/s during the Stage 2, and then flattens out in Stage 3. The strain increases again during Stage 4 of the time series. The pressure also changes in stages. It is unchanged during Stage 1 and then increases relatively rapidly at a rate of between 4 and 10 Pa/s during Stage 2 and continues to increase during Stages 3 and 4, but at a rate slightly slower than during Stage 2 (from 2.5 to 3.5 Pa/s). The rates of pressure and strain change are consistent with basic analyses. The rate of pressurization is similar to the calculated velocity of flow in the saturated zone, and this is similar to the measured velocity of the strain increase. The ratio of the pressurization rate and strain rate during Stage 2 (2.5 (Pa/s) /0.03 (μ/s)) is 80 MPa, which is approximately equal to the upper range expected considering the uniaxial compression and the Young’s Modulus of the sand. An interesting effect occurs when the upper surface of the sand becomes saturated. Significant compression (several 10s of ) occurs as the pressure and saturation increase at the upper surface. Compression occurs throughout the column, but the effect is greatest at the top of the column. The rate of compression is fastest slightly before ponding occurs, but it slows markedly and nearly stops when water starts to accumulate at the surface (ponding). This effect reverses when the surface of the sand is drained, resulting in tension throughout the column. The compression caused by this effect can be as large or larger than the tensile strain that accumulated during filling. This effect was unexpected because increasing pressure is normally associated with tensile strain. Nevertheless, this effect was observed consistently in all tests where the pressure changes at the upper surface, including tests where water was injected from below or infiltrated from above. This effect behaves as if the pore pressure at the upper surface of the sand exerts a normal force on the boundary (increasing pore pressure exerts a downward compression on the boundary). The results of these laboratory experiments indicate that CMPI can measure strain caused by fluid pressure changes in the vadose zone with a signal to noise ratio of 100 or more. Injection and drainage cause a strain signal that is complex, but repeatable. The magnitudes of the strain signal are consistent with magnitudes that are expected based on poroelastic calculations using independently measured properties of the sand. These results indicate that the CMPI technique with an optical fiber laminated in a polyester ribbon generates data that represent the strain distribution during hydrologic processes of imbibition and drainage in variably saturated sand. This suggests that distributed strain measurements using CMPI could be a viable approach for evaluating processes in the vadose zone, laying the groundwork for future field implementation.

47 OTHER INSTRUMENTATION↗