Search NASA⌕ Search

SEARCH · Search NASA

Results for “PRESSURE DISTRIBUTION - FLOAT BOTTOMS”

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 55 records · Page 3

Description and calibration of the Langley 6- by 19-inch transonic tunnel

A description and calibration is presented of the Langley 6- by 19-inch transonic tunnel which is a two-dimensional facility with top and bottom slotted walls used for testing two-dimensional airfoil sections. Basic tunnel-empty Mach number distributions and schlieren flow photographs as well as integrated normal-force coefficients, pitching-moment coefficients, surface-pressure distributions, and schlieren flow photographs of an NACA 0012 airfoil calibration model are presented. The Mach number capability of the facility is from 0.5 to about 1.1 with a corresponding Reynolds number range of 1.5 million to 3 million based on a 4.0-in. model chord. Comparisons of experimental results from the tests with previous data are also presented.

Ladson, C. L.↗

Jupiter's atmospheric composition and cloud structure deduced from absorption bands in reflected sunlight

The spectrum of sunlight reflected by Jupiter is analyzed by comparing observations of Woodman (1979) with multiple-scattering computations. The analysis yields information on the vertical cloud structure at several latitudes and on the abundance of CH4 and NH3 in the atmosphere of Jupiter. The abundances of CH4 and NH3 suggest that all ices and rocks are overabundant on Jupiter by a factor of 2 or more, providing an important constraint on models for the formation of Jupiter from the primitive solar nebula. The pressure level of the clouds, the gaseous NH3 abundance, the mean temperature profile, and the Clausius-Clapeyron relation suggest that these clouds are predominantly ammonia crystals with the cloud bottom at 600-700 mb. A diffuse distribution of aerosols exists between 150 and 500 mb, and the spectral variation of albedo reflects a changing bulk absorption coefficient of the material composing the aerosols and is diagnostic of the aerosol composition.

Sato, M.↗

A swept wing panel in a low speed flexible walled test section

The testing of two-dimensional airfoil sections in adaptive wall tunnels is relatively widespread and has become routine at all speeds up to transonic. In contrast, the experience with the three-dimensional testing of swept panels in adaptive wall test sections is very limited, except for some activity in the 1940's at NPL, London. The current interest in testing swept wing panels led to the work covered by this report, which describes the design of an adaptive-wall swept-wing test section for a low speed wind tunnel and gives test results for a wing panel swept at 40 deg. The test section has rigid flat sidewalls supporting the panel, and features flexible top and bottom wall with ribs swept at the same angle as the wing. When streamlined, the walls form waves swept at the same angle as the wing. The C sub L (-) curve for the swept wing, determined from its pressure distributions taken with the walls streamlined, compare well with reference data which was taken on the same model, unswept, in a test section deep enough to avoid wall interference.

Goodyer, M. J.↗

Geochemical Evidence for a Terrestrial Magma Ocean

The aftermath of phase separation and crystal-liquid fractionation in a magma ocean should leave a planet geochemically differentiated. Subsequent convective and other mixing processes may operate over time to obscure geochemical evidence of magma ocean differentiation. On the other hand, core formation is probably the most permanent, irreversible part of planetary differentiation. Hence the geochemical traces of core separation should be the most distinct remnants left behind in the mantle and crust, In the case of the Earth, core formation apparently coincided with a magma ocean that extended to a depth of approximately 1000 km. Evidence for this is found in high pressure element partitioning behavior of Ni and Co between liquid silicate and liquid iron alloy, and with the Ni-Co ratio and the abundance of Ni and Co in the Earth's upper mantle. A terrestrial magma ocean with a depth of 1000 km will solidify from the bottom up and first crystallize in the perovskite stability field. The largest effect of perovskite fractionation on major element distribution is to decrease the Si-Mg ratio in the silicate liquid and increase the Si-Mg ratio in the crystalline cumulate. Therefore, if a magma ocean with perovskite fractionation existed, then one could expect to observe an upper mantle with a lower than chondritic Si-Mg ratio. This is indeed observed in modern upper mantle peridotites. Although more experimental work is needed to fully understand the high-pressure behavior of trace element partitioning, it is likely that Hf is more compatible than Lu in perovskite-silicate liquid pairs. Thus, perovskite fractionation produces a molten mantle with a higher than chondritic Lu-Hf ratio. Arndt and Blichert-Toft measured Hf isotope compositions of Barberton komatiites that seem to require a source region with a long-lived, high Lu-Hf ratio. It is plausible that that these Barberton komatiites were generated within the majorite stability field by remelting a perovskite-depleted part of the upper mantle transition zone.

Agee, Carl B.↗

Results of wind tunnel tests of an ASRM configured 0.03 scale Space Shuttle integrated vehicle model (47-OTS) in the AEDC 16-foot transonic wind tunnel, volume 2

An experimental Aerodynamic and Aero-Acoustic loads data base was obtained at transonic Mach numbers for the Space Shuttle Launch Vehicle configured with the ASRM Solid Rocket Boosters as an increment to the current flight configuration (RSRB). These data were obtained during transonic wind tunnel tests (IA 613A) conducted in the Arnold Engineering Development Center 16-Foot transonic propulsion wind tunnel from March 27, 1991 through April 12, 1991. This test is the first of a series of two tests covering the Mach range from 0.6 to 3.5. Steady state surface static and fluctuating pressure distributions over the Orbiter, External Tank and Solid Rocket Boosters of the Shuttle Integrated Vehicle were measured. Total Orbiter forces, Wing forces and Elevon hinge moments were directly measured as well from force balances. Two configurations of Solid Rocket Boosters were tested, the Redesigned Solid Rocket Booster (RSRB) and the Advanced Solid Rocket Motor (ASRM). The effects of the position (i.e., top, bottom, top and bottom) of the Integrated Electronics Assembly (IEA) box, mounted on the SRB attach ring, were obtained on the ASRM configured model. These data were obtained with and without Solid Plume Simulators which, when used, matched as close as possible the flight derived pressures on the Orbiter and External Tank base. Data were obtained at Mach numbers ranging from 0.6 to 1.55 at a Unit Reynolds Number of 2.5 million per foot through model angles of attack from -8 to +4 degrees at sideslip angles of 0, +4 and -4 degrees.

Marroquin, J.↗

Results of wind tunnel tests of an ASRM configured 0.03 scale Space Shuttle integrated vehicle model (47-OTS) in the AEDC 16-foot Transonic wind tunnel (IA613A), volume 1

An experimental Aerodynamic and Aero-Acoustic loads data base was obtained at transonic Mach numbers for the Space Shuttle Launch Vehicle configured with the ASRM Solid Rocket Boosters as an increment to the current flight configuration (RSRB). These data were obtained during transonic wind tunnel tests (IA 613A) conducted in the Arnold Engineering Development Center 16-Foot transonic propulsion wind tunnel from March 27, 1991 through April 12, 1991. This test is the first of a series of two tests covering the Mach range from 0.6 to 3.5. Steady state surface static and fluctuating pressure distributions over the Orbiter, External Tank and Solid Rocket Boosters of the Shuttle Integrated Vehicle were measured. Total Orbiter forces, Wing forces and Elevon hinge moments were directly measured as well from force balances. Two configurations of Solid Rocket Boosters were tested, the Redesigned Solid Rocket Booster (RSRB) and the Advanced Solid Rocket Motor (ASRM). The effects of the position (i.e. top, bottom, top and bottom) of the Integrated Electronics Assembly (IEA) box, mounted on the SRB attach ring, were obtained on the ASRM configured model. These data were obtained with and without Solid Plume Simulators which, when used, matched as close as possible the flight derived pressures on the Orbiter and External Tank base. Data were obtained at Mach numbers ranging from 0.6 to 1.55 at a Unit Reynolds Number of 2.5 million per foot through model angles of attack from -8 to +4 degrees at sideslip angles of 0, +4 and -4 degrees.

Marroquin, J.↗

Comparison of Water-Load Distributions Obtained during Seaplane Landings with Bureau of Aeronautics Specifications. TED No. NACA 2413

Bureau of Aeronautics Design Specifications SS-IC-2 for water loads in sheltered water are compared with experimental water loads obtained during a full--scale landing investigation. This investigation was conducted with a JRS-1 flying boat which has a 20 degrees dead-rise V-bottom with a partial chine flare. The range of landing conditions included airspeeds between 88 and 126 feet per second, sinking speeds between 1.6 and 9.1 feet per second, flight angles less than 6 degrees, and trims between 2 degrees and 12 degrees. Landings were moderate and were made in calm water. Measurements were obtained of maximum over-all loads, maximum pitching moments, and pressure distributions. Maximum experimental loads include over-all load factors of 2g, moments of 128,000 pound-feet, and maximum local pressures greater than 40 pounds per square inch. Experimental over-all loads are approximately one-half the design values, while local pressures are of the same order as or larger than pressures calculated from specifications for plating, stringer, floor, and frame design. The value of this comparison is limited, to some extent, by the moderate conditions of the test and by the necessary simplifying assumptions used in comparing the specifications with the experimental loads.

Smiley, Robert F.↗

In-Flight Boundary-Layer Transition on a Large Flat Plate at Supersonic Speeds

A flight experiment was conducted to investigate the pressure distribution, local flow conditions, and boundary-layer transition characteristics on a large flat plate in flight at supersonic speeds up to Mach 2.0. The primary objective of the test was to characterize the local flow field in preparation for future tests of a high Reynolds number natural laminar flow test article. The tests used a F-15B testbed aircraft with a bottom centerline mounted test fixture. A second objective was to determine the boundary-layer transition characteristics on the flat plate and the effectiveness of using a simplified surface coating for future laminar flow flight tests employing infrared thermography. Boundary-layer transition was captured using an onboard infrared imaging system. The infrared imagery was captured in both analog and digital formats. Surface pressures were measured with electronically scanned pressure modules connected to 60 surface-mounted pressure orifices. The local flow field was measured with five 5-hole conical probes mounted near the leading edge of the test fixture. Flow field measurements revealed the local flow characteristics including downwash, sidewash, and local Mach number. Results also indicated that the simplified surface coating did not provide sufficient insulation from the metallic structure, which likely had a substantial effect on boundary-layer transition compared with that of an adiabatic surface. Cold wall conditions were predominant during the acceleration to maximum Mach number, and warm wall conditions were evident during the subsequent deceleration. The infrared imaging system was able to capture shock wave impingement on the surface of the flat plate in addition to indicating laminar-to-turbulent boundary-layer transition.

Banks, Daniel W.↗

Estimation of Hydrodynamic Impact Loads and Pressure Distributions on Bodies Approximating Elliptical Cylinders with Special Reference to Water Landings of Helicopters

An approximate method for computing water loads and pressure distributions on lightly loaded elliptical cylinders during oblique water impacts is presented. The method is of special interest for the case of emergency water landings of helicopters. This method makes use of theory developed and checked for landing impacts of seaplanes having bottom cross sections of V and scalloped contours. An illustrative example is given to show typical results obtained from the use of the proposed method of computation. The accuracy of the approximate method was evaluated through comparison with limited experimental data for two-dimensional drops of a rigid circular cylinder at a trim of 0 degrees and a flight -path angle of 90 degrees. The applicability of the proposed formulas to the design of rigid hulls is indicated by the rough agreement obtained between the computed and experimental results. A detailed computational procedure is included as an appendix.

Schnitzer, Emanuel↗

Ocean foam generation and modeling

A laboratory investigation was conducted to determine the physical and microwave properties of ocean foam. Special foam generators were designed and fabricated, using porous glass sheets, known as glass frits, as the principal element. The glass frit was sealed into a water-tight vertical box, a few centimeters from the bottom. Compressed air, applied to the lower chamber, created ocean foam from sea water lying on the frit. Foam heights of 30 cm were readily achieved, with relatively low air pressures. Special photographic techniques and analytical procedures were employed to determine foam bubble size distributions. In addition, the percentage water content of ocean foam was determined with the aid of a particulate sampling procedure. A glass frit foam generator, with pore diameters in the range 70 - 100 micrometers, produced foam with bubble distributions very similar to those found on the surface of natural ocean foam patches.

Porter, R. A.↗

Numerical Simulation of No Vent Chill and Fill of a Large Liquid Hydrogen Tank

On orbit transfer of cryogenic liquids is necessary for many future space missions. The techniques required to transfer cryogens in micro gravity are quite different from ground system applications. In normal gravity, cryogenic tanks are usually filled from the bottom at nearly atmospheric pressure. The vapor generated due to heat transfer from warm tank walls, is allowed to vent from the top of the tank while the tank is being filled. If the normal gravity technique is used on orbit, the uncertainty of liquid and vapor distributions in unsettled condition may result in dumping of large amounts of liquid propellant overboard. The Vented Chill and No Vent Fill (VCNVF) process is a methodology [1] used to reduce propellant loss by holding the vent valve open during initial wall chilldown, preventing over-pressurization of the tank and stalling of the incoming flow, but closing the tank vent valve after chilling the tank wall sufficiently. An accurate simulation of no vent tank chill and fill is necessary to estimate the amount of propellant necessary to chill the tank, determine the maximum sustainable flowrate to fill the tank without stalling, and establish a timeline for the filling process. Numerical modeling of filling a cryogenic tank is complex because the flow is two-phase and the process involves both boiling and condensation with heat transfer between solid and fluid. In recent years, progress has been made to develop a numerical modeling technique [2,3,4] that has been verified by comparing with the test data from a small flight tank using liquid nitrogen [5]. The purpose of the present paper is to apply this technique in a larger tank with liquid hydrogen and compare with the test data of Chato [1]. The K-site Test Facility [1], shown in Figure 1, consists of a Test Tank, spray system, instrumentation, and the vacuum chamber. The test tank selected was ellipsoidal with an 87 inch major diameter and a 1.2 to 1 major to minor axis ratio. The two ends are joined by a short 1.5 inch cylindrical section. The tank is made of 2219 aluminum chemically milled to a nominal thickness of 0.087 inches. The tank weighs 329.25 pounds, and the tank volume is 175 ft 3 . The GFSSP model of the Test Tank is shown in Figure 2. Node 13 is a boundary node that represents the supply tank which is supplying liquid hydrogen at a given pressure and temperature. The total measured flowrate is evenly distributed through Branches 131 through 139. Nodes 1 through 9 represent the test tank and they are connected with metal solid nodes 14 through 22 through fluid-to-solid conductors that allow convective heat transfer between the solid and fluid nodes. The predicted tank pressure, filling rate, wall temperature and fluid temperature will be compared with test data.

Alok Majumdar↗

The abundances of methane and ortho/para hydrogen on Uranus and Neptune: Implications of New Laboratory 4-0 H2 quadrupole line parameters

The tropospheric methane molar fraction (f(sub CH4, t) and the ortho/para hydrogen ratio are derived for Uranus and Neptune based on new determinations of spectroscopic parameters for key hydrogen features as reported by D. W./ Ferguson et al. (1993). For each planet, the relatively weak laboratory linestrengths (approximately 30 and 15% less than the theoretical 4-0 S(0) and S(1) linestrengths, respectively) results, when compared to analyses adopting theroetical values, in a approximately 30% decrease in the tropospheric methane ratio and a comparable increase in the pressure level of the optically thick cloudtop marking the bottom of the visible atmosphere (P(sub c/d)). The increase in the ratio of S(1)/S(0) linestrengths from 4.4 (theoretical) to approximately 5.9 (measured) results in a decrease in the range of viable ortho/para ratios; an equilibrium hydrogen distribution is now the best fit for both planets. The methane mixing ratios reported here are in agreement with the value of 0.023 derived by the Voyager Radio Occultation Experiment (G. F. Lindal, 1992) for Neptune, but slightly lower than the Voyager Uranus measurement of 0.023 reported by G. F. LIndel et al. (1987). The relative carbon-to-hydrogen abundances for Uranus and Neptune support planetary formation mechanisms involving the dissolution of carbon-bearing planetesimals in the atmospheres of both planets during their early stages of formation (e.g., J. B. Pollack et al., 1986).

Baines, Kevin H.↗

Buoyancy-Driven Heat Transfer During Application of a Thermal Gradient for the Study of Vapor Deposition at Low Pressure Using and Ideal Gas

A mathematical model has been developed to determine heat transfer during vapor deposition of source materials under a variety of orientations relative to gravitational accelerations. The model demonstrates that convection can occur at total pressures as low as 10-2 mm Hg. Through numerical computation, using physical material parameters of air, a series of time steps demonstrates the development of flow and temperature profiles during the course of vapor deposition. These computations show that in unit gravity vapor deposition occurs by transport through a fairly complicated circulating flow pattern when applying heat to the bottom of the vessel with parallel orientation with respect to the gravity vector. The model material parameters for air predict the effect of kinematic viscosity to be of the same order as thermal diffusivity, which is the case for Prandtl number approx. 1 fluids. Qualitative agreement between experiment and the model indicates that 6-(2-methyl-4-nitroanilino)-2,4-hexadiyn-l-ol (DAMNA) at these pressures indeed approximates an ideal gas at the experiment temperatures, and may validate the use of air physical constants. It is apparent that complicated nonuniform temperature distribution in the vapor could dramatically affect the homogeneity, orientation, and quality of deposited films. The experimental test i's a qualitative comparison of film thickness using ultraviolet-visible spectroscopy on films generated in appropriately oriented vapor deposition cells. In the case where heating of the reaction vessel occurs from the top, deposition of vapor does not normally occur by convection due to a stable stratified medium. When vapor deposition occurs in vessels heated at the bottom, but oriented relative to the gravity vector between these two extremes, horizontal thermal gradients induce a complex flow pattern. In the plane parallel to the tilt axis, the flow pattern is symmetrical and opposite in direction from that where the vessel is positioned vertically. The ground-based experiments are sufficient preliminary tests of theory and should be of significant interest regarding vapor deposited films in microgravity.

Frazier, D. O.↗

Electric current in a unipolar sunspot with an untwisted field

The return flux (RF) sunspot model is applied to a round, unipolar sunspot observed by H. Kawakami (1983). Solving the magnetohydrostatic problem using the gas pressure deficit between the umbral and quiet-sun atmospheres as a source function, a distribution of electric current density in an untwisted, unipolar sunspot as a function of height and radial distance from the sunspot center is observed. Maximum electric current density is about 32 mA/sq m at the bottom of the sunspot.

Osherovich, V. A.↗

Distributed Capacitive Sensor for Sample Mass Measurement

Previous robotic sample return missions lacked in situ sample verification/ quantity measurement instruments. Therefore, the outcome of the mission remained unclear until spacecraft return. In situ sample verification systems such as this Distributed Capacitive (DisC) sensor would enable an unmanned spacecraft system to re-attempt the sample acquisition procedures until the capture of desired sample quantity is positively confirmed, thereby maximizing the prospect for scientific reward. The DisC device contains a 10-cm-diameter pressure-sensitive elastic membrane placed at the bottom of a sample canister. The membrane deforms under the weight of accumulating planetary sample. The membrane is positioned in close proximity to an opposing rigid substrate with a narrow gap. The deformation of the membrane makes the gap narrower, resulting in increased capacitance between the two parallel plates (elastic membrane and rigid substrate). C-V conversion circuits on a nearby PCB (printed circuit board) provide capacitance readout via LVDS (low-voltage differential signaling) interface. The capacitance method was chosen over other potential approaches such as the piezoelectric method because of its inherent temperature stability advantage. A reference capacitor and temperature sensor are embedded in the system to compensate for temperature effects. The pressure-sensitive membranes are aluminum 6061, stainless steel (SUS) 403, and metal-coated polyimide plates. The thicknesses of these membranes range from 250 to 500 m. The rigid substrate is made with a 1- to 2-mm-thick wafer of one of the following materials depending on the application requirements glass, silicon, polyimide, PCB substrate. The glass substrate is fabricated by a microelectromechanical systems (MEMS) fabrication approach. Several concentric electrode patterns are printed on the substrate. The initial gap between the two plates, 100 m, is defined by a silicon spacer ring that is anodically bonded to the glass substrate. The fabricated proof-of-concept devices have successfully demonstrated tens to hundreds of picofarads of capacitance change when a simulated sample (100 g to 500 g) is placed on the membrane.

Toda, Risaku↗

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↗

A parametric study of effect of forebody shape on flow angularity at Mach 8

Flow angularity and static pressure measurements have been made on the lower surface of nine forebody models that simulate the bottom forward surface of a hypersonic aircraft. Measurements were made in an area of the forebody that represents the location of an inlet of a scramjet engine. A parametric variation of the forebody surface investigated the effect of: (1) spanwise curvature; (2) longitudinal curvature; and (3) planform shape on both flow angularity and static pressure distribution. Results of each of the three parametric variations of geometry were compared to those for the same flat delta forebody. Spanwise curvature results showed that a concave shape and the flat delta had the lowest flow angularity and lowest rate of increase in flow angularity with angle of attack. Longitudinal curvature results showed a convex surface to give the better flow at the higher angles of attack. The better of the two planform shapes tested was a convex elliptical shape. Limited flow field calculations were made at angles of attack using a three dimensional, method-of-characteristics program. In general, at all angles of attack there was agreement between data and theory.

Johnson, C. B.↗

LLNL Kimberlina 1.2 NUFT Simulations June 2018 (v2)

This dataset contains the output 6,000, 3-dimensional reactive multi-phase flow and transport aquifer simulations of brine and CO2 leakage into a protective aquiver in California’s San Joaquin Valley and input data files detailing the geologic mesh, aquifer physical properties and CO2 and brine injection rates. This data set was generated as an ongoing effort with the US DOE National Risk Assessment Partnership (NRAP) to evaluate the effectiveness of monitoring techniques to detect brine and CO2 leakage from legacy wells into underground sources of drinking water overlaying a CO2 storage reservoir. Each simulation contains a unique set of input parameters, generated stochastically. The outputs consist of these upper three geologic layers (from top): the Etchegoin, Macoma-Chanac, Santa Margarita-McLure formations. These simulations span the several distances (1, 3 and 6 km or wells W31-0.2, W31-0.5 and W31-1.0, respectively) from the CO2 injector, initiated from bottom hole pressure and saturation to calculate wellbore leakage from the storage reservoir, with low and high regional groundwater gradients and wellbore leakage into 5 leaky nodes. The dataset includes 1,000 unique simulations for each distance, which each contain a unique aquifer heterogeneity, aquifer and caprock permeability, and two model generations are included with a high permeability (prod07) and hybrid permeability (prod09). The range of permeability distributions is listed in Table 1. Each model generation consists of 3,000 simulations. Included in the dataset are the leakage rates determined from 2D wellbore models which utilize the pressure and CO2 saturation from LBL's reservoir simulations, NUFT mesh files with distributed lithology, NUFT rocktab files which describe the material properties for the geologic layers and the NUFT input files and post-processed output 'ntab' files. Each ntab file contains spatial (rows) and temporal (columns) model output tables for each model cell, the locations (x,y,z) and dimensions for each cells (dx, dy, dz). Table 1. Permeability distribution ranges for prod07 and prod09 model generations Geologic Layer: Permeability Range (log10 m^2) prod07 prod09 Etchegoin -12.92 to -10.92 -13.70 to -11.44 Macoma-Chanac -12.72 to -10.72 -13.50 to -11.24 Santa Margarita-McLure -12.70 to -10.70 -13.48 to -11.22 The input files used to generate the model include which are included in the dataset are: Time series of CO2 leakage input into the model (ex: Q_brn.W31-0.2.sim1000.layers123.tab) Time series of CO2 leakage input into the model (ex: Q_CO2.W31-0.2.sim1000.layers123.tab) Physical properties of the aquifer materials detailing the aquifer porosity, solid density, partitioning coefficients, permeabilities and van-Genuchten parameters detailed in a NUFT rocktab file: (ex: sim1000.usnt.rocktab) Numerical mesh and geologic data assigned to each model cell detailed in a NUFT genmsh format (ex: sim1000.mesh_k16.prod07.trans.genmsh) The primary output parameters are: pH (use absolute value) Change in TDS (mg/kg) Change in Pressure (Pa) Change CO2 gas saturation (fraction range 0.0-1.0) for example, the directory /p/lscratchh/mansoor1/nrap/kimberlina/prod09/mainfiles/sim1000/W31- 0.2 contains: sim1000.W31-0.2.trans.pH.red.ntab sim1000.W31-0.2.no_bg.trans.TDS.red.ntab sim1000.W31-0.2.usnt.P.deltabg.red.ntab sim1000.W31-0.2.usnt.CO2_sat.deltabg.red.ntab Each row in the NTAB files consist of model output per numerical grid cell. Each output file contains 33 columns (variables), including the information of numerical records, geologic location and sizes and the simulated parameter values over time. The first 13 variables are about numerical records and relative geologic information for a simulation grid: 1. index: simulation index 2. i: the ith grid of x-axis 3. j: the ith grid of y-axis 4. k: the ith grid of z-axis 5. element_ref: element reference 6. nuft_ind: nuft index 7. x: grid location in the x axis direction 8. y: grid location in the y axis direction 9. z: grid location in the z axis direction 10. dx: grid length in the x axis direction 11. dy: grid length in the y axis direction 12. dz: grid length in the z axis direction 13. volume: volume of the simulation grid The remainder (14, 15, 16...) variables are the simulated parameter values over time, take Pressure as an example, are: 14. 0.0y: initial pressure per cell. 15. 10.0y: simulated pressure at the end of the 10th year. 16. 20.0y: simulated pressure at the end of the 20th year. ... (repeated for every 10 years until 200 years)... The model extends 10,000 m, 5,000 m and 1,411 m in the x,y and z dimensions, respectively. The mesh consists of 164,832 cells with mesh dimensions of 101 x 51 x 32 (nx, ny, nz), with cell dimensions ranging from 100 m laterally (along x and y-axis) and model layers are as designated in the z-axis: Layer 1: atmosphere (1e-30 m thick) Layer 2: upper caprock (10 m thick) Layers 3-13: Etchegoin (536.23 m thck) Layers 14-27: Macoma-Chanac (679.04 m thick) Layers 28-32: Santa Margarita-McLure (185.94 m thick) The wellbore is placed along node i=51, j=26, and extends vertically along 5 nodes from the top to the bottom of the model. Special instructions when extracting files: Each Gzip archive (ex: prod07.sim1000-sim00099.tar.gz) contains 100 simulations. Gzip archives should be transferred into base directories (ie. In Linux: mkdir prod07; mv prod07.*.tar.gz prod07/.) before extracting, or files will be overwritten. Each sub-simulation tree should have the following file structure pattern (using the linux 'tree' command): |-- prod07 | |-- sim0001 | |-- W31-0.2 | | |-- Q_brn.W31-0.2.sim0001.layers123.tab | | |-- Q_co2.W31-0.2.sim0001.layers123.tab | | |-- sim0001.W31-0.2.no_bg.trans.TDS.red.ntab | | |-- sim0001.W31-0.2.trans.pH.red.ntab | | |-- sim0001.W31-0.2.usnt.CO2_sat.deltabg.red.ntab | | |-- sim0001.W31-0.2.usnt.P.deltabg.red.ntab | |-- W31-0.5 | | |-- Q_brn.W31-0.5.sim0001.layers123.tab | | |-- Q_co2.W31-0.5.sim0001.layers123.tab | | |-- sim0001.W31-0.5.no_bg.trans.TDS.red.ntab | | |-- sim0001.W31-0.5.trans.pH.red.ntab | | |-- sim0001.W31-0.5.usnt.CO2_sat.deltabg.red.ntab | | |-- sim0001.W31-0.5.usnt.P.deltabg.red.ntab | |-- W31-1.0 | | |-- Q_brn.W31-1.0.sim0001.layers123.tab | | |-- Q_co2.W31-1.0.sim0001.layers123.tab | | |-- sim0001.W31-1.0.no_bg.trans.TDS.red.ntab | | |-- sim0001.W31-1.0.trans.pH.red.ntab | | |-- sim0001.W31-1.0.usnt.CO2_sat.deltabg.red.ntab | | |-- sim0001.W31-1.0.usnt.P.deltabg.red.ntab | |-- sim0001.mesh_k16.prod07.trans.genmsh Disclaimer This document was prepared as an account of work sponsored by an agency of the United States government. Neither the United States government nor Lawrence Livermore National Security, LLC, nor any of their employees makes any warranty, expressed or implied, or assumes any legal liability or responsibility for the accuracy, completeness, or usefulness of any information, apparatus, product, or process disclosed, or represents that its use would not infringe privately owned rights. Reference herein to any specific commercial product, process, or service by trade name, trademark, manufacturer, or otherwise does not necessarily constitute or imply its endorsement, recommendation, or favoring by the United States government or Lawrence Livermore National Security, LLC. The views and opinions of authors expressed herein do not necessarily state or reflect those of the United States government or Lawrence Livermore National Security, LLC, and shall not be used for advertising or product endorsement purposes. Lawrence Livermore National Laboratory is operated by Lawrence Livermore National Security, LLC, for the U.S. Department of Energy, National Nuclear Security Administration under Contract DE-AC52-07NA27344. This report was reviewed and released as LLNL-MI-753464.

aquifer↗