Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Data Flow Reversal”

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 307 records · Page 17

Latitudinal Transport of Angular Momentum by Cellular Flows Observed with MDI

We have analyzed Doppler velocity images from the MDI instrument on SOHO to determine the latitudinal transport of angular momentum by the cellular photospheric flows. Doppler velocity images from 60-days in May to July of 1996 were processed to remove the p-mode oscillations, the convective blue shift, the axisymmetric flows, and any instrumental artifacts. The remaining cellular flows were examined for evidence of latitudinal angular momentum transport. Small cells show no evidence of any such transport. Cells the size of supergranules (30,000 km in diameter) show strong evidence for a poleward transport of angular momentum. This would be expected if supergranules are influenced by the Coriolis force, and if the cells are elongated in an east-west direction. We find good evidence for just such an east-west elongation of the supergranules. This elongation may be the result of differential rotation shearing the cellular structures. Data simulations of this effect support the conclusion that elongated supergranules transport angular momentum from the equator toward the poles, Cells somewhat larger than supergranules do not show evidence for this poleward transport. Further analysis of the data is planned to determine if the direction of angular momentum transport reverses for even larger cellular structures. The Sun's rapidly rotating equator must be maintained by such transport somewhere within the convection zone.

Hathaway, David H.↗

Estimates of the wavenumber wavelet power spectrum of magnetic fluctuations during magnetic reconnection

Fluctuation analyses of experimental observations generally lack high temporal resolution and are in frequency-space f, contrary to theoretical efforts in wavenumber-space k. This is due to the inherent limits of the Fourier transform, though it is prominent due to the ease of diagnostic implementation. Advances in wavelet-based analysis have provided relief due to its temporal resolution, but in its common use, is still hard to compare to theoretical models. By using the two-point correlation technique in conjunction with large data sets, a wavelet power spectrum in wavenumber-space can be created. Dubbed the wavenumber wavelet power spectrum, this spectrum relates wavenumber to power in time. Further, this analysis technique more closely connects characterizations of experimentally observed fluctuations with other system parameters and theoretical predictions. In this article, we develop the wavenumber wavelet power spectrum using magnetic fluctuations caused by tearing instability driven magnetic reconnection in reproducible, high temperature laboratory plasmas. These dynamic magnetic fluctuations generated in reversed field pinch plasmas are broadband, ranging from the low frequency, 10's of kHz, up to the ion gyroradii frequencies, 100's of kHz. The dominant fluctuations have poloidal and toroidal mode numbers (m,n)=(1,6−10) and can grow to 2%–3% of the mean magnetic field. During these reconnection events, ions, and electrons are energized, magnetic fluctuation amplitudes increase, plasma flow is halted, and the toroidal magnetic flux increases, all on a semi-periodic basis. The newly developed spectrum provides better temporal resolution of spectrum characteristics to correlate with these particle energization phenomena.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Aerosol Meteorology of the Maritime Continent for the 2012 7SEAS Southwest Monsoon Intensive Study - Part 1: Regional-scale Phenomena

The largest 7 Southeast Asian Studies (7SEAS) operation period within the Maritime Continent (MC) occurred in the August-September 2012 biomass burning season. Included was an enhanced deployment of Aerosol Robotic Network (AERONET) sun photometers, multiple lidars, and field measurements to observe transported smoke and pollution as it left the MC and entered the southwest monsoon trough. Here we describe the nature of the overall 2012 southwest monsoon (SWM) and biomass burning season to give context to the 2012 deployment. The MC in 2012 was in a slightly warm El Nino Southern Oscillation (ENSO) phase and with spatially typical burning activity. However, overall fire counts for 2012 were 10 lower than the Reid et al. (2012) baseline, with regions of significant departures from this norm, ranging from southern Sumatra (+30) to southern Kalimantan (42). Fire activity and monsoonal flows for the dominant burning regions were modulated by a series of intraseasonal oscillation events (e.g., Madden-Julian Oscillation, or MJO, and boreal summer intraseasonal oscillation, or BSISO). As is typical, fire activity systematically progressed eastward over time, starting with central Sumatran fire activity in June related to a moderately strong MJO event which brought drier air from the Indian Ocean aloft and enhanced monsoonal flow. Further burning in Sumatra and Kalimantan Borneo occurred in a series of significant events from early August to a peak in the first week of October, ending when the monsoon started to migrate back to its wintertime northeastern flow conditions in mid-October. Significant monsoonal enhancements and flow reversals collinear with tropical cyclone (TC) activity and easterly waves were also observed. Islands of the eastern MC, including Sulawesi, Java, and Timor, showed less sensitivity to monsoonal variation, with slowly increasing fire activity that also peaked in early October but lingered into November. Interestingly, even though fire counts were middling, resultant AERONET 500nm aerosol optical thickness (AOT) from fire activity was high, with maximums of 3.6 and 5.6 in the Sumatra and Kalimantan source regions at the end of the burning season and an average of approximately 1. AOTs could also be high at receptor sites, with a mean and maximum of 0.57 and 1.24 in Singapore and 0.61 and 0.8 in Kuching Sarawak. Ultimately, outside of the extreme 2015 El Nino event, average AERONET AOT values were higher than any other time since sites were established. Thus, while satellite fire data, models, and AERONET all qualitatively agree on the nature of smoke production and transport, the MC's complex environment resulted in clear differences in quantitative interpretation of these datasets.

Reid, Jeffrey S.↗

Network-Aware and Welfare-Maximizing Dynamic Pricing for Energy Sharing: Preprint

The proliferation of behind-the-meter (BTM) distributed energy resources (DER) within the electrical distribution network presents significant supply and demand flexibilities, but also introduces operational challenges such as voltage spikes and reverse power flows. In response, this paper proposes a network-aware dynamic pricing framework tailored for energy-sharing coalitions that aggregate small, but ubiquitous, BTM DER downstream of a distribution system operator's (DSO) revenue meter that adopts a generic net energy metering (NEM) tariff. By formulating a Stackelberg game between the energy-sharing market leader and its prosumers, we show that the dynamic pricing policy induces the prosumers toward a network-safe operation and decentrally maximizes the energysharing social welfare. The dynamic pricing mechanism involves a combination of a locational ex-ante dynamic price and an ex-post allocation, both of which are functions of the energy sharing's BTM DER. The ex-post allocation is proportionate to the price differential between the DSO NEM price and the energy sharing locational price. Simulation results using real DER data and the IEEE 13-bus test systems illustrate the dynamic nature of network-aware pricing at each bus, and its impact on voltage.

energy communities↗

Network-Aware and Welfare-Maximizing Dynamic Pricing for Energy Sharing

The proliferation of behind-the-meter (BTM) distributed energy resources (DER) within the electrical distribution network presents significant supply and demand flexibilities, but also introduces operational challenges such as voltage spikes and reverse power flows. In response, this paper proposes a network-aware dynamic pricing framework tailored for energy-sharing coalitions that aggregate small, but ubiquitous, BTM DER downstream of a distribution system operator's (DSO) revenue meter that adopts a generic net energy metering (NEM) tariff. By formulating a Stackelberg game between the energy-sharing market leader and its prosumers, we show that the dynamic pricing policy induces the prosumers toward a network-safe operation and decentrally maximizes the energysharing social welfare. The dynamic pricing mechanism involves a combination of a locational ex-ante dynamic price and an ex-post allocation, both of which are functions of the energy sharing's BTM DER. The ex-post allocation is proportionate to the price differential between the DSO NEM price and the energy-sharing locational price. Simulation results using real DER data and the IEEE 13-bus test systems illustrate the dynamic nature of network-aware pricing at each bus, and its impact on voltage.

aggregates↗

Instantaneous velocity field imaging instrument for supersonic reacting flows

The technical tasks conducted to develop and demonstrate a new gas velocity measurement technique for high enthalpy reacting flows is described. The technique is based on Doppler-shifted Planar Laser-induced Fluorescence (PLIF) imaging of the OH radical. The imaging approach permits, in principle, single-shot measurements of the 2-D distribution of a single velocity component in the measurement plane, and is thus a technique of choice for applications in high enthalpy transient flow facilities. In contrast to previous work in this area, the present program demonstrated an approach which modified the diagnostic technique to function under the constraints of practical flow conditions of engineering interest, rather than vice-versa. In order to accomplish the experimental demonstrations, the state-of-the-art in PLIF diagnostic techniques was advanced in several ways. Each of these tasks is described in detail and is intended to serve as a reference in supporting the transition of this new capability to the fielded PLIF instruments now installed at several national test facilities. Among the new results of general interest in LlF-based flow diagnostics, a detailed set of the first measurements of the collisional broadening and shifting behavior of OH (1,0) band transitions in H7-air combustion environments is included. Such measurements are critical in the design of a successful strategy for PLIF velocity imaging; they also relate to accurate concentration and temperature measurements, particularly in compressible flow regimes. Furthermore, the results shed new light on the fundamental relationship between broadening and energy transfer collisions in OH A(sup 2)Sigma(+)v(sup ') = 1. The first single-pulse, spectrally-resolved measurements of the output of common pulsed dye lasers were also produced during the course of this effort. As with the OH broadening measurements, these data are a significant aspect of a successful velocity imaging strategy, and also have potential implications for many other LIF measurement techniques. Our results indicated the need to modify the commercially available laser cavity in order to accommodate the constraints imposed by typical SCRAMJET combustion characteristics as well as to increase the instrument's velocity dynamic range to span an intra-image range in excess of 2 km/s. The various technical efforts were brought together in a series of experiments demonstrating the applicability of the technique in a high pressure, high temperature H2-air combustion system. The resultant images were compared with 2-D flow simulations in order to determine the accuracy of the instrument. Mean velocity imaging in flows with an axis of symmetry was demonstrated with an accuracy of +/- 50 m/s out of an intra-image dynamic range of 1600 m/s, including reversed flow. A more complex configuration amenable to single-shot imaging in flows without an axis of symmetry was also demonstrated. Limitations imposed by available equipment resulted in an accuracy of about +/- 200 m/s out of 1750 m/s in these demonstrations. Minor modifications to the present configuration were suggested to improve this performance. Each technical task is described in detail, along with significance of the results for the overall imaging velocimeter configuration. This report should allow the user community to integrate this new measurement capability in their existing instrumentation platforms.

Allen, M. G.↗

Observations of the reappearance of polar coronal holes and the reversal of the polar magnetic field

Observational data relating to the evolution of the polar magnetic field around sunspot maximum is examined. Particular emphasis is given to coronal hole observations performed during the last two solar maxima. Long-term averages of the latitudinal dependence of the photospheric magnetic field and the evolutionary pattern of the polar crown filaments are used to trace the poleward motion of the reversal of the large-scale field and are compared to the redevelopment of polar holes. Within the context of phenomenological models of the solar cycle, it is concluded that: (1) the process of polarity reversal and redevelopment of polar holes is discontinuous, with surges of flux of old-cycle polarity interrupting the poleward migration of new-cycle flux; (2) contrary to the Babcock (1961) hypothesis, the polar crown disappears months after the magnetic pole reversal; and (3) the observations support suggestions of a poleward meridional flow around solar maximum that cannot be accounted for by Leighton-type (1964) diffusion.

Webb, D. F.↗

Rayleigh Scattering Density Measurements, Cluster Theory, and Nucleation Calculations at Mach 10

In an exploratory investigation, quantitative unclustered laser Rayleigh scattering measurements of density were performed in the air in the NASA Langley Research Center's 31 in. Mach 10 wind tunnel. A review of 20 previous years of data in supersonic and Mach 6 hypersonic flows is presented where clustered signals typically overwhelmed molecular signals. A review of nucleation theory and accompanying nucleation calculations are also provided to interpret the current observed lack of clustering. Data were acquired at a fixed stagnation temperature near 990Kat five stagnation pressures spanning 2.41 to 10.0 MPa (350 to 1454 psi) using a pulsed argon fluoride excimer laser and double-intensified charge-coupled device camera. Data averaged over 371 images and 210 pixels along a 36.7mmline measured freestream densities that agree with computed isentropic-expansion densities to less than 2% and less than 6% at the highest and lowest densities, respectively. Cluster-free Mach 10 results are compared with previous clustered Mach 6 and condensation-free Mach 14 results. Evidence is presented indicating vibrationally excited oxygen and nitrogen molecules are absorbed as the clusters form, release their excess energy, and inhibit or possibly reverse the clustering process. Implications for delaying clustering and condensation onset in hypersonic and hypervelocity facilities are discussed.

Balla, R. Jeffrey↗

DC Electric Fields and Associated Plasma Drifts Observed with the C/NOFS Satellite

Initial DC electric field observations and associated plasma drifts are presented from the Vector Electric Field Investigation (VEFI) on the Air Force Communication/Navigation Outage Forecasting System (C/NOFS) satellite. We present statistical averages of the vector fields for the first year of operations that include both the zonal and radial components of the resulting E x B plasma flows at low latitudes. Magnetic field data from the VEFI science magnetometer are used to compute the plasma flows. The DC electric field detector reveals zonal and radial electric fields that undergo strong diurnal variations, typically displaying eastward and outward-directed fields during the day and westward and downward-directed fields at night. There is considerable variation in the large scale DC electric field data, in both the daytime and nighttime cases, with enhanced structures typically observed at night. In general, the measured zonal DC electric field amplitudes include excursions that extend within the 0.4 - 2 m V/m range, corresponding to E x B drifts of the order of 30-150 m/s. The average vertical or radial electric fields may exceed the zonal fields in amplitude by a factor of 1.5 to 2. Although the data compare well, in a general sense, with previous satellite observations and statistical patterns of vertical ion drifts, the E x B drifts we report from C/NOFS rarely show a pronounced pre-reversal enhancement after sunset. We attribute this to a combination of extreme solar minimum conditions and the fact that the C/NOFS orbit of 401 by 867 km carries the probes essentially above the lower altitude regions where the wind-driven dynamo might be expected to create enhanced upwards drifts in the early evening. Evidence for wavenumber 4 tidal effects and other longitudinal signatures have been detected and will be presented. We also discuss off-equatorial electric fields and their relation to the ambient plasma density.

Pfaff, R.↗

Microtopographic evolution of lava flows at Cima volcanic field, Mojave Desert, California

Microtopographic profiles were measured and power spectra calculated for dated lava flow surfaces at Cima volcanic field in the eastern Mojave Desert of California in order to quantify changes in centimeter- to meter-scale roughness as a function of age. For lava flows younger than about 0.8 m.y., roughness over all spatial scales decreases with age, with meter-scale roughness decreasing slightly more than centimeter scales. Flows older than about 0.8 m.y. show a reversal of this trend, becoming as rough as young flows at these scales. Modeling indicates that eolian deposition can explain most of the change observed in the offset, or roughness amplitude, of power spectra of flow surface profiles up to 0.8 m.y. Other processes, such as rubbing and stone pavement development, appear to have a minor effect in this age range. Changes in power spectra of surfaces older than about 0.8 m.y. are consistent with roughening due to fluvial dissection. These results agree qualitatively with a process-response model that attributes systematic changes in flow surface morphology to cyclic changes in the rates of eolian, soil formation, and fluvial processes. Identification of active surficial processes and estimation of the extent of their effects, or stage of surficial evolution, through measurement of surface roughness will help put the correlation of surficial units on a quantitative basis. This may form the basis for the use of radar remote sensing data to help in regional correlations of surficial units.

Farr, Tom G.↗

Hybrid Reduced Order Model for N2-N2 Interactions for Application to Dissociation and Energy Transfer Processes

Recent work in the aerothermodynamics community has focused on the development of reduced-order models for thermo-chemical non-equilibrium which avoid the restrictive assumptions of multi-temperature models and the prohibitive cost associated with State-to-State (StS) models. In the present work, this is accomplished by lumping energy states together and assuming groups of states are roughly in equilibrium. As a result, the non-equilibrium behavior of a gas can be captured at a reduced computational cost from a full StS simulation. In this work, we present a hybrid grouping model for studying energy transfer and dissociation in a mixture of nitrogen molecules due to N2-N2 reactions. This is accomplished by making use of a grouping strategy informed by data from the N2-N StS kinetic data. However, due to the massive computational cost associated with constructing StS data for the N2-N2 system, the kinetic data for the hybrid grouping model are calculated using the quasi-classical trajectory (QCT) method by sampling states for trajectory within the groups. This general framework is called the Maximum-Entropy Quasi-Classical Trajectory (ME-QCT) method. The primary challenge associated with this method is that rates for reverse grouped reactions cannot be obtained through detailed balance at a group level, due to the variation of group internal temperatures. To construct the full model for N2-N2 grouped kinetics using the ME-QCT method, detailed balance is invoked at the microscopic level, allowing for the calculation of the full kinetic data from QCT. Results will be presented using the full ME-QCT model for the N2-N2 system in an isothermal and isochoric reactor simulation. In addition, simple CFD test cases for a one-dimensional standing shock and for a quasi-one-dimensional nozzle will be used for demonstration of the ME-QCT method. This method allows for the calculation of non-equilibrium behavior for the N2-N2 system without the prohibitive cost of a full StS simulation. Moreover, it enables the construction of a unified model for the dissociating and recombining non-equilibrium flows.

nonequilibrium modeling↗

Reaction Temperature Manipulation as a Process Intensification Approach for CO 2 Absorption

Reactor temperature manipulation to increase product yields of chemical reactions is a known technique used in many industrial processes. In the case of exothermic chemical reactions, the well-known Le Chatelier’s principle predicts that a decrease in temperature will displace the chemical reaction toward the formation of products by increasing the value of the equilibrium constant. The reverse is true for endothermic reactions. Reactor temperature manipulation in an industrial system, however, affects the values of many variables, including physical properties, transport parameters, reaction kinetic parameters, etc. In the case of reactive absorption, some variables change with increasing temperatures due to solute absorption, while others change in such a way that the solute absorption rate decreases. For example, temperature drop increases product formation for exothermic reactions but reduces the value of transport parameters, leading to decreasing interfacial concentrations and absorption rates. Therefore, temperature manipulation strategies must be designed carefully to achieve the process goals. In this work, we theoretically study the use of temperature as a tool to increase CO 2 absorption by solvents in a semi-batch reactor. A computer code has been developed and validated using reported experimental data. Calculated results demonstrate an increase in absorbed CO 2 of more than 28% with respect to the highest temperature used. Despite high agitation and high gas flow rate, the system is mass transfer controlled at short times, becoming kinetically controlled as time increases. An operating strategy to decrease cooling energy costs is also proposed. This study reveals that reactor temperature manipulation can be an effective process to improve CO 2 absorption by solvents in two-phase semi-batch reactors.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Coherent fluctuations in a partially magnetized small magnetron discharge operating on molecular nitrogen

We report experimental data on rotating plasma disturbances in a quasi-steady direct-current magnetron discharge operating with molecular nitrogen at pressures of 27, 40, 67, and 133 Pa. Experiments with a segmented anode measuring the local current flow reveal up to 12 coherent structures (“spokes”) propagating in the +E×B and −E×B azimuthal directions. The current–voltage characteristics are studied and found to follow a negative resistance profile at the lower discharge currents. The existence of coherent structures and the number of spokes observed depend on the discharge current and operating pressure. Lower pressures are more favorable to the retrograde −E×B spoke rotation, which is observed in the 27, 40, and 67 Pa discharges. The results are interpreted with a gradient-driven drift wave model, and using the dispersion relation, we show that the propagation direction depends on the orientation of the local electric field within the plasma. Using an electron fluid model, we show how experimental conditions can lead to field reversals that generate −E×B spoke rotation in the lower pressure discharges.

Physics↗

Reduction of rotor-turbulence interaction noise in static fan noise testing

Three effective inlet cleanup methods are described. A flared reverse cone inlet is used to eliminate wakes from the fan casings and/or probe supports. Boundary layer suction is employed ahead of the fan rotor and on the outer flare of the cone to reduce the boundary layer turbulence and remove any residual wakes. To reduce the midstream turbulent intensity and length scales, a turbulence control structure is used which is constructed with both a layer of honeycomb and a fine mesh screen. The effects of these cleanup methods are quantified by measuring the far-field noise in an anechoic chamber, using a high-speed 20-in.-diam fan of the current high passband type. The changes in the turbulent field impinging on this rotor are quantified by mapping the streamwise and transverse turbulent properties (spectra, intensity and length scale) with crossed hot film probes. It is convincingly shown that it is possible to clean up the inlet flow of a static fan noise test facility to a point where the static acoustic data simulate flight data.

Kantola, R. A.↗

Seasonal Cycles of Meridional Overturning and Heat Transport of the Indian Ocean

A general circulation model of the Indian Ocean is fitted to monthly averaged climatological temperatures, salinities, and surface fluxes using the adjoint method. Interannual variability is minimized by penalizing the temporal drift from one seasonal cycle to another during a two-year integration. The resultant meridional overturning and heat transport display large seasonal variations, with maximum amplitudes of 18 and 22 (x 10(exp 6) cubic m/s) for the overturning and 1.8 and 1.4 (x 10(exp 15) W) for heat transport near 10 S and 10 N, respectively. A dynamical decomposition of the overturning and heat transport shows that the time-varying Ekman How plus its barotropic compensation can explain a large part of the seasonal variations in overturning and heat transport. The maximum variations at 10 deg N and 10 deg S are associated with monsoon reversal over the northern Indian Ocean and changes of the easterlies over the southern Indian Ocean. An external mode with variable topography has a moderate contribution where the Somali Current and the corresponding gyre reverse direction seasonally. Contribution front vertical shear (thermal wind and ageostrophic shear) is dominant near the southern boundary and large near the Somali Current latitudes. The dominant balance in the zonally integrated heat budget is between heat storage change and heat transport convergence except south of 15 S. Optimization with seasonal forcings improves estimates of sea surface temperatures, but the annual average overturning and heat transport are very similar to previous results with annual mean forcings. The annual average heat transport consists of roughly equal contributions from time-mean and time-varying fields of meridional velocities and temperatures in the northern Indian Ocean. indicating a significant rectification to the heat transport due to the time-varying fields. The time-mean and time-varying contributions are primarily due to the overturning and horizontal gyre, respectively. Inclusion of TOPEX data enhances the seasonal cycles of the estimated overturning and heat transport in the central Indian Ocean significantly and improves the estimated equatorial zonal flows but leads 10 unrealistic estimates of the velocity structure near the Indonesian Throughflow region, most likely owing to the deficiencies in the lateral boundary conditions.

Lee, Tong↗

Generative models on phase space

Deep generative models such as diffusion and flow matching are powerful machine learning tools capable of learning and sampling from high-dimensional distributions. They are particularly useful when the training data appears to be concentrated on a submanifold of the data embedding space. For high-energy physics data, consisting of collections of relativistic energy-momentum 4-vectors, this submanifold can enforce extremely strong physically-motivated priors, such as energy and momentum conservation. If these constraints are learned only approximately, rather than exactly, this can inhibit the interpretability and reliability of such generative models. To remedy this deficiency, we introduce generative models which are, by construction, confined at every step of their sampling trajectory to the manifold of massless N-particle Lorentz-invariant phase space in the center-of-momentum frame. In the case of diffusion models, the "pure noise" forward process endpoint corresponds to the uniform distribution on phase space, which provides a clear starting point from which to identify how correlations among the particles emerge during the reverse (de-noising) process. We demonstrate that our models are able to learn both few-particle and many-particle distributions with various singularity structures, paving the way for future interpretability studies using generative models trained on simulated jet data.

Bogorad, Zachary [Fermilab]↗

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↗

Seven hole probe measurement of leading edge vortex flows

This paper discusses the use of a seven-hole probe on measurements of leading edge vortices of highly sweep delta wing planforms. Intrusive probe data taken with the pressure probe were compared with nonintrusive measurements made with laser Doppler anemometry system. In addition to probe size, the natural position of breakdown and the sweep angle of the wing are also factors in determining sensitivity of the flow to probe interference. At low angles of attach vortex breakdown does not occur in the vicinity of the model and the seven hole probe was found to yield reasonably accurate measurements. When the angle of attack of the model was increased so that vortex breakdown was near the trailing edge, introducing the probe over the wing would cause the breakdown position to move ahead of the probe. However, when breakdown naturally occurred ahead of the mid-chord of the wing the vortices were found to be less sensitive to a probe placed behind the breakdown point. Vortex breakdown on a lower swept wing is found to be more sensitive to interference. Near the breakdown region, seven hole probe measurement is less accurate due to a combination of probe interference and flow reversal.

Payne, F. M.↗