Search NASA⌕ Search

SEARCH · Search NASA

Results for “height function”

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 1,513 records · Page 84
All sourcesNASA NTRSDOE OSTI
NASA NTRS2015

A Microwave Radiometric Method to Obtain the Average Path Profile of Atmospheric Temperature and Humidity Structure Parameters and Its Application to Optical Propagation System Assessment

The values of the key atmospheric propagation parameters Ct2, Cq2, and Ctq are highly dependent upon the vertical height within the atmosphere thus making it necessary to specify profiles of these values along the atmospheric propagation path. The remote sensing method suggested and described in this work makes use of a rapidly integrating microwave profiling radiometer to capture profiles of temperature and humidity through the atmosphere. The integration times of currently available profiling radiometers are such that they are approaching the temporal intervals over which one can possibly make meaningful assessments of these key atmospheric parameters. Since these parameters are fundamental to all propagation conditions, they can be used to obtain Cn2 profiles for any frequency, including those for an optical propagation path. In this case the important performance parameters of the prevailing isoplanatic angle and Greenwood frequency can be obtained. The integration times are such that Kolmogorov turbulence theory and the Taylor frozen-flow hypothesis must be transcended. Appropriate modifications to these classical approaches are derived from first principles and an expression for the structure functions are obtained. The theory is then applied to an experimental scenario and shows very good results.

Atmospheric turbulence↗
DOE OSTI2024

Proton discrimination in CLYC for fast neutron spectroscopy

The Cs 2 LiYCl 6 :Ce (CLYC) elpasolite scintillator is known for its response to fast and thermal neutrons along with good γ-ray energy resolution. While the 35 Cl(n,p) reaction has been identified as a potential means for CLYC-based fast neutron spectroscopy in the absence of time-of-flight (TOF), previous efforts to functionalize CLYC as a fast neutron spectrometer have been thwarted by the inability to isolate proton interactions from 6 Li(n,α) and 35 Cl(n,α) signals. This work introduces a new approach to particle discrimination in CLYC for fission spectrum neutrons using a multi-gate charge integration algorithm that provides excellent separation between protons and heavier charged particles. Neutron TOF data were collected using a 252 Cf source, an array of EJ-309 organic liquid scintillators, and a 6 Li-enriched CLYC scintillator outfitted with fast electronics. Modal waveforms were constructed corresponding to the different reaction channels, revealing significant differences in the pulse characteristics of protons and heavier charged particles at ultrafast, fast, and intermediate time scales. These findings informed the design of a pulse shape discrimination algorithm, which was validated using the TOF data. This study also proposes an iterative subtraction method to mitigate contributions from confounding reaction channels in proton and heavier charged particle pulse height spectra, opening the door for CLYC-based fast neutron and γ-ray spectroscopy while preserving sensitivity to thermal neutron capture signals.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS2003

Residual Strength Analyses of Monolithic Structures

Finite-element fracture simulation methodology predicts the residual strength of damaged aircraft structures. The methodology uses the critical crack-tip-opening-angle (CTOA) fracture criterion to characterize the fracture behavior of the material. The CTOA fracture criterion assumes that stable crack growth occurs when the crack-tip angle reaches a constant critical value. The use of the CTOA criterion requires an elastic- plastic, finite-element analysis. The critical CTOA value is determined by simulating fracture behavior in laboratory specimens, such as a compact specimen, to obtain the angle that best fits the observed test behavior. The critical CTOA value appears to be independent of loading, crack length, and in-plane dimensions. However, it is a function of material thickness and local crack-front constraint. Modeling the local constraint requires either a three-dimensional analysis or a two-dimensional analysis with an approximation to account for the constraint effects. In recent times as the aircraft industry is leaning towards monolithic structures with the intention of reducing part count and manufacturing cost, there has been a consistent effort at NASA Langley to extend critical CTOA based numerical methodology in the analysis of integrally-stiffened panels.In this regard, a series of fracture tests were conducted on both flat and curved aluminum alloy integrally-stiffened panels. These flat panels were subjected to uniaxial tension and during the test, applied load-crack extension, out-of-plane displacements and local deformations around the crack tip region were measured. Compact and middle-crack tension specimens were tested to determine the critical angle (wc) using three-dimensional code (ZIP3D) and the plane-strain core height (hJ using two-dimensional code (STAGS). These values were then used in the STAGS analysis to predict the fracture behavior of the integrally-stiffened panels. The analyses modeled stable tearing, buckling, and crack branching at the integral stiffener using different values of critical CTOA for different material thicknesses and orientation. Comparisons were made between measured and predicted load-crack extension, out-of-plane displacements and local deformations around the crack tip region. Simultaneously, three-dimensional capabilities to model crack branching and to monitor stable crack growth of multiple cracks in a large thick integrally-stiffened flat panels were implemented in three-dimensional finite element code (ZIP3D) and tested by analyzing the integrally-stiffened panels tested at Alcoa. The residual strength of the panels predicted from STAGS and ZP3D code compared very well with experimental data. In recent times, STAGS software has been updated with new features and now one can have combinations of solid and shell elements in the residual strength analysis of integrally-stiffened panels.

Forth, Scott↗
DOE OSTI2019

Properties of HPK UFSD after neutron irradiation up to 6e15 n/cm 2

In this article we report results from a neutron irradiation campaign of Ultra-Fast Silicon Detectors (UFSD) with fluences of 1e14, 3e14, 6e14, 1e15, 3e15 and 6e15 neq/cm 2 . The UFSD used in this study are circular 50 μm thick Low-Gain Avalanche Detectors (LGAD), with a 1.0 mm diameter active area. Hamamatsu Photonics (HPK), Japan, produced the UFSD with pre-irradiation internal gain in the range 5–70 depending on the bias voltage. The sensors were tested pre-irradiation and post-irradiation with minimum ionizing particles (MIPs) from a 90 Sr β-source. The leakage current, internal gain and the timing resolution were measured as a function of bias voltage at -20 °C and -30 °C. The timing resolution of each device under test was extracted from the time difference with a second calibrated UFSD in coincidence, using the constant fraction discriminator (CFD) method for both. The dependence of the gain upon the irradiation fluence is consistent with the acceptor removal mechanism; at -20 °C the highest gain decreases from 70 before radiation to 2 after a fluence of 6e15 n/cm 2 . Consequently, the timing resolution was found to deteriorate from 20 ps to 50 ps. The results indicate that the most accurate time resolution is obtained varying with fluence the CFD value used to determine the time of arrival, from 0.1 for pre-irradiated sensors to 0.6 at the highest fluence. Key changes to the pulse shape induced by irradiation, i.e. (i) the contribution of charge multiplication not limited to the gain layer zone, (ii) the shortening of the rise time and (iii) the reduced pulse height, were compared with the WF2 simulation program and were found to be in agreement.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS1999

Ames Optimized TCA Configuration

Configuration design at Ames was carried out with the SYN87-SB (single block) Euler code using a 193 x 49 x 65 C-H grid. The Euler solver is coupled to the constrained (NPSOL) and the unconstrained (QNMDIF) optimization packages. Since the single block grid is able to model only wing-body configurations, the nacelle/diverter effects were included in the optimization process by SYN87's option to superimpose the nacelle/diverter interference pressures on the wing. These interference pressures were calculated using the AIRPLANE code. AIRPLANE is an Euler solver that uses a unstructured tetrahedral mesh and is capable of computations about arbitrary complete configurations. In addition, the buoyancy effects of the nacelle/diverters were also included in the design process by imposing the pressure field obtained during the design process onto the triangulated surfaces of the nacelle/diverter mesh generated by AIRPLANE. The interference pressures and nacelle buoyancy effects are added to the final forces after each flow field calculation. Full details of the (recently enhanced) ghost nacelle capability are given in a related talk. The pseudo nacelle corrections were greatly improved during this design cycle. During the Ref H and Cycle 1 design activities, the nacelles were only translated and pitched. In the cycle 2 design effort the nacelles can translate vertically, and pitch to accommodate the changes in the lower surface geometry. The diverter heights (between their leading and trailing edges) were modified during design as the shape of the lower wing changed, with the drag of the diverter changing accordingly. Both adjoint and finite difference gradients were used during optimization. The adjoint-based gradients were found to give good direction in the design space for configurations near the starting point, but as the design approached a minimum, the finite difference gradients were found to be more accurate. Use of finite difference gradients was limited by the CPU time limit available on the Cray machines. A typical optimization run using finite difference gradients can use only 30 to 40 design variables and one optimization iteration within the 8 hour queue limit for the chosen grid size and convergence level. The efficiency afforded by the adjoint method allowed for 50-120 design variables and 5-10 optimization iterations in the 8 hour queue. Geometric perturbations to the wing and fuselage were made using the Hicks/Henne (HH) shape functions. The HH functions were distributed uniformly along the chords of the wing defining sections and lofted linearly. During single-surface design, constraints on thickness and volume at selected wing stations were imposed. Both fuselage camber and cross-sectional area distributions were permitted to change during design. The major disadvantage to the use of these functions is the inherent surface waviness produced by repeated use of such functions. Many smoothing operations were required following optimization runs to produce a configuration with reasonable smoothness. Wagner functions were also used on the wing sections but were never used on the fuselage. The Wagner functions are a family of increasingly oscillatory functions that have also been used extensively in airfoil design. The leading and trailing edge regions of the wing were designed by use of polynomial and monomial functions respectively. Twist was attempted but was abandoned because of little performance improvement available from changing the baseline twist.

Cliff, Susan E.↗
DOE OSTI2021

Implementation of Perdew–Zunger self-interaction correction in real space using Fermi–Löwdin orbitals

Most widely used density functional approximations suffer from self-interaction error, which can be corrected using the Perdew–Zunger (PZ) self-interaction correction (SIC). We implement the recently proposed size-extensive formulation of PZ-SIC using Fermi–Löwdin Orbitals (FLOs) in real space, which is amenable to systematic convergence and large-scale parallelization. Here, we verify the new formulation within the generalized Slater scheme by computing atomization energies and ionization potentials of selected molecules and comparing to those obtained by existing FLOSIC implementations in Gaussian based codes. The results show good agreement between the two formulations, with new real-space results somewhat closer to experiment on average for the systems considered. We also obtain the ionization potentials and atomization energies by scaling down the Slater statistical average of SIC potentials. The results show that scaling down the average SIC potential improves both atomization energies and ionization potentials, bringing them closer to experiment. Finally, we verify the present formulation by calculating the barrier heights of chemical reactions in the BH6 dataset, where significant improvements are obtained relative to Gaussian based FLOSIC results.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2017

Aero-MINE (Motionless INtegrated Energy) for Distributed Scalable Wind Power.

The proposed Aero-MINE technology will extract energy from wind without any exterior moving parts. Aero-MINEs can be integrated into buildings or function stand-alone, and are scalable. This gives them advantages similar to solar panels, but with the added benefit of operation in cloudy or dark conditions. Furthermore, compared to solar panels, Aero-MINEs can be manufactured at lower cost and with less environmental impact. Power generation is isolated internally by the pneumatic transmission of air and the outlet air-jet nozzles amplify the effectiveness. Multiple units can be connected to one centrally located electric generator. Aero-MINEs are ideal for the built-environment, with numerous possible configurations ranging from architectural integration to modular bolt-on products. Traditional wind turbines suffer from many fundamental challenges. The fast-moving blades produce significant aero-acoustic noise, visual disturbances, light-induced flickering and impose wildlife mortality risks. The conversion of massive mechanical torque to electricity is a challenge for gears, generators and power conversion electronics. In addition, the installation, operation and maintenance of wind turbines is required at significant height. Furthermore, wind farms are often in remote locations far from dense regions of electricity customers. These technical and logistical challenges add significantly to the cost of the electricity produced by utility-scale wind farms. In contrast, distributed wind energy eliminates many of the logistical challenges. However, solutions such as micro-turbines produce relatively small amounts of energy due to the reduction in swept area and still suffer from the motion-related disadvantages of utility-scale turbines. Aero-MINEs combine the best features of distributed generation, while eliminating the disadvantages.

17 WIND ENERGY↗
NASA NTRS1997

Topography of the Moon from the Clementine Lidar

Range measurements from the lidar instrument carried aboard the Clementine spacecraft have been used to produce an accurate global topographic model of the Moon. This paper discusses the function of the lidar; the acquisition, processing, and filtering of observations to produce a global topographic model; and the determination of parameters that define the fundamental shape of the Moon. Our topographic model: a 72nd degree and order spherical harmonic expansion of lunar radii, is designated Goddard Lunar Topography Model 2 (GLTM 2). This topographic field has an absolute vertical accuracy of approximately 100 m and a spatial resolution of 2.5 deg. The field shows that the Moon can be described as a sphere with maximum positive and negative deviations of approx. 8 km, both occurring on the farside, in the areas of the Korolev and South Pole-Aitken (S.P.-Aitken) basins. The amplitude spectrum of the topography shows more power at longer wavelengths as compared to previous models, owing to more complete sampling of the surface, particularly the farside. A comparison of elevations derived from the Clementine lidar to control point elevations from the Apollo laser altimeters indicates that measured relative topographic heights generally agree to within approx. 200 in over the maria. While the major axis of the lunar gravity field is aligned in the Earth-Moon direction, the major axis of topography is displaced from this line by approximately 10 deg to the cast and intersects the farside 24 deg north of the equator. The magnitude of impact basin topography is greater than the lunar flattening (approx. 2 km) and equatorial ellipticity (approx. 800 m), which imposes a significant challenge to interpreting the lunar figure. The floors of mare basins are shown to lie close to an equipotential surface, while the floors of unflooded large basins, except for S.P.-Aitken, lie above this equipotential. The radii of basin floors are thus consistent with a hydrostatic mechanism for the absence of significant farside maria except for S.P.-Aitken, whose depth and lack of mare require significant internal compositional and/or thermal heterogeneity. A macroscale surface roughness map shows that roughness at length scales of 10(exp 1) - 10(exp 2) km correlates with elevation and surface age.

Smith, David E.↗
NASA NTRS2004

Residual Strength Characterization of a Curved Integrally-Stiffened Panel

Over the years, Finite-element fracture simulation methodology has been very well established at NASA Langley to predict the residual strength of damaged aircraft structures. This methodology has been experimentally verified at NASA Langley for structures ranging from laboratory coupons up to full-scale built-up structural components with single and multiple-site damage cracking. The methodology uses the critical crack-tip-opening-angle (CTOA) fracture criterion to characterize the fracture behavior of the material. The CTOA fracture criterion assumes that stable crack growth occurs when the crack-tip angle reaches a constant critical value. The use of the CTOA criterion requires an elastic-plastic, finite-element analysis. The critical CTOA value is determined by simulating fracture behavior in laboratory specimens, such as a compact specimen, to obtain the angle that best fits the observed test behavior. The critical CTOA value appears to be independent of loading, crack length, and in-plane dimensions. However, it is a function of material thickness and local crack-front constraint. Modeling the local constraint requires either a three-dimensional analysis or a two-dimensional analysis with an approximation to account for the constraint effects. In recent times as the aircraft industry is leaning towards monolithic structures with the intension of reducing part count and manufacturing cost, there has been a consistent effort at NASA Langley to extend critical CTOA based numerical methodology in the analysis of integrally-stiffened panels. In this regard, a series of fracture tests were conducted on curved aluminum-alloy integrally-stiffened panels. These curved panels were subjected to uniaxial tension and pressure loading. During the test, applied load-crack extension, out-of-plane displacements and local deformations around the crack tip region were measured. Compact and middle-crack tension specimens were tested to determine the critical angle (psi(sub c) using three-dimensional code (ZIP3D) and the plane-strain core height (h(sub c)) using two-dimensional code (STAGS).

Seshadri, B. R.↗
NASA NTRS1937

Full-scale Wind-tunnel and Flight Tests of a Fairchild 22 Airplane Equipped with External-airfoil Flaps

Wind-tunnel and flight tests have been made of a Fairchild 22 airplane equipped with a wing having external-airfoil flaps that also perform the function of ailerons. Lift, drag, and pitching-moment coefficients of the airplane with several flap settings, and the rolling- and yawing-moment coefficients with the flaps deflected as ailerons were measured in the full-scale tunnel with the horizontal tail surfaces and propeller removed. The effect of the flaps on the low speed and on the take-off and landing characteristics, the effectiveness of flaps when used as ailerons, and the forces required to operate them as ailerons were determined in flight. The wind-tunnel tests showed that the flaps increased the maximum lift coefficient of the airplane from 1.51 with the flap in the minimum drag position to 2.12 with the flap in the minimum drag position to 2.12 with the flap deflected 30 degrees. In the flight tests the minimum speed decreased from 46.8 miles per hour with the flaps up to 41.3 miles per hour with the flaps deflected. The required take-off run to attain a height of 50 feet was reduced from 820 to 750 feet and the landing run from a height of 50 feet was reduced from 930 to 480 feet. The flaps for this installation gave lateral control that was not entirely satisfactory. Their rolling action was good but the adverse yaw resulting from their use was greater than is considerable, and the stick forces required to operate them increased too rapidly with speed.

Reed, Warren D↗
DOE OSTI2021

Fermi-Löwdin-orbital self-interaction correction using the optimized-effective-potential method within the Krieger-Li-Iafrate approximation

Perdew-Zunger self-interaction correction (PZ-SIC) offers a route to remove self interaction errors on an orbital-by-orbital basis. A recent formulation of PZ-SIC by Pederson, Ruzsinszky and Perdew proposes restricting the unitary transformation to localized orbitals called Fermi-Löwdin orbitals. This formulation, called the FLOSIC method, simplifies PZ-SIC calculations and was implemented self-consistently using a Jacobi-like (FLOSIC-Jacobi) iteration scheme. In this work we implement the FLOSIC approach using the Krieger-Li-Iafrate (KLI) approximation to the optimized effective potential (OEP). We compare the results of present FLOSIC-KLI approach with FLOSIC-Jacobi scheme for atomic energies, atomization energies, ionization energies, barrier heights, polarizability of chains of hydrogen molecules etc. to validate the FLOSIC-KLI approach. The FLOSICKLI approach, which is within the realm of Kohn-Sham theory, predicts smaller energy gaps between frontier orbitals due to the lowering of eigenvalues of the lowest unoccupied orbitals. Results show that atomic energies, atomization energies, ionization energy as an absolute of highest occupied orbital eigenvalue, and polarizability of chains of hydrogen molecules between the two methods agree within 2%. Lastly the FLOSIC-KLI approach is used to determine the vertical ionization energies of water clusters.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS2012

How Did I Get Here?

What country was responsible for much of the early work in laser altimetry, who was the first scientist to look at a vegetation profile generated by a laser, who first related laser ranges to forest structure? The purpose of this paper is to look back 3+ decades ago when forestry lidar was in its formative years and to give credit to those who first thought to employ a laser to measure the Earth's surface. As with much scientific research, advances were often made independently and concurrently in a number of countries. Functioning lasers were first demonstrated in 1960 in the USA and in 1961 in the USSR, but research into the use of lasers as forest measurement tools did not begin for another IS years. Initially, with respect to Earth resources, lasers were employed to measure sea ice surface roughness, to make near-shore bathymetric measurements, to penetrate forests to make detailed topographic measurements, and to fluoresce oceanic phytoplankton for surface current studies. Some of these early studies noted that forest profiles were evident but in fact added noise to topographic retrievals. Trees became the signal rather than noise in 1976 when researchers in the Canadian Forest Service went to an IUFRO conference in Oslo and reported on efforts to better estimate timber volume in tropical forests. They employed large-scale airphotos (1:500 up to 1:4000) to measure the top of the tropical canopy, an avionic radar to measure the location of the ground beneath the canopy, and a barometric sensor to record aircraft/radar height above a datum (e.g., sea level). Their radar did not work well in dense vegetation and mentioned in passing that a laser altimeter had worked well in Canadian forests. But in a report 2 years later they reported that, on a study area in Costa Rica, the radar ground line was much improved and no mention is made of a lidar altimeter. In the USSR in 1977, Russian researchers felled a birch and a spruce, aimed a helium-neon, 0.63 micron laser with a spot size of approx 25mm, at the horizontal trees, produced a profilograph, compared it to tape measurements, and concluded that, with increased power, such a laser could be mounted on an aircraft to remotely measure forest canopies. In 1979, they mounted their He-Ne laser on an AN-2 biplane and acquired their fIrSt airborne profiles. Scientists with TRANARG, a mapping/surveying company in Caracas, Venezuela, reported on 1976 flights employing a helium-neon lidar to collect over 11,000 km of lidar profiles spaced 1.5 km apart in order to construct a topographic map to help site a new reservoir. They noted 35-40 m median canopy heights with emergents up to 55 m in their profiles. These studies and others that utilize these height and canopy density metrics for forest mensuration are reviewed in this glance backwards at the history of forestry lidar.

Nelson, Ross↗
NASA NTRS2023

Extragalactic Magnetism with SOFIA (SALSA Legacy Program). V. First Results on the Magnetic Field Orientation of Galaxies

We present the analysis of the magnetic field ( B -field) structure of galaxies measured with far-infrared (FIR) and radio (3 and 6 cm) polarimetric observations. We use the first data release of the Survey of extragALactic magnetiSm with SOFIA of 14 nearby (<20 Mpc) galaxies with resolved (5"-18" ; 90 pc–1 kpc) imaging polarimetric observations using SOFIA/HAWC+ from 53 to 214 μm . We compute the magnetic pitch-angle (𝜓 B ) profiles as a function of the galactocentric radius. We introduce a new magnetic alignment parameter ( ζ ) to estimate the disordered-to-ordered ratio of spiral B -fields. We find FIR and radio wavelengths to not generally trace the same B -field morphology in galaxies. The 𝜓 B profiles tend to be more ordered across all galactocentric radii in radio ( ζ 6cm = 0.93 + 0.03) than in FIR ( ζ 154μm = 0.84 + 0.14). For spiral galaxies, FIR B -fields are 2%–75% more turbulent than the radio B -fields. For starburst galaxies, we find that FIR polarization is a better tracer of the B -fields along the galactic outflows than radio polarization. Our results suggest that the B -fields associated with dense, dusty, turbulent star-forming regions (those traced at FIR) are less ordered than warmer, less dense regions (those traced at radio) of the interstellar medium. The FIR B -fields seem to be more sensitive to the activity of the star-forming regions and molecular clouds within a vertical height of a few hundred parsecs in the disk of spiral galaxies than the radio B -fields.

magnetic field (B-field)↗
NASA NTRS2023

Extragalactic Magnetism with SOFIA (SALSA Legacy Program) - V. First Results on the Magnetic Field Orientation of Galaxies

We present the analysis of the magnetic field (B-field) structure of galaxies measured with far-infrared (FIR) and radio (3 and 6 cm) polarimetric observations. We use the first data release of the Survey of extragALactic magnetiSm with SOFIA (SALSA) of 14 nearby (< 20 Mpc) galaxies with resolved (5′′ − 18′′; 90 pc–1 kpc) imaging polarimetric observations using SOFIA/HAWC+ from 53 to 214 μm. We compute the magnetic pitch angle (Ψ B ) profiles as a function of the galactocentric radius. We introduce a new magnetic alignment parameter (ζ) to estimate the disordered-to-ordered ratio of spiral B -fields. We find FIR and radio wavelengths to not generally trace the same B -field morphology in galaxies. The Ψ B profiles tend to be more ordered across all galactocentric radii in radio (ζ 6cm = 0.93 ± 0.03) than in FIR (ζ 154μm = 0.84 ± 0.14). For spiral galaxies, FIR B-fields are 2−75% more turbulent than the radio B -fields. For starburst galaxies, we find that FIR polarization is a better tracer of the B -fields along the galactic outflows than radio polarization. Our results suggest that the B -fields associated with dense, dusty, turbulent star-forming regions, those traced at FIR, are less ordered than warmer, less-dense regions, those traced at radio, of the interstellar medium. The FIR B -fields seem to be more sensitive to the activity of the star-forming regions and molecular clouds within a vertical height of few hundred pc in the disk of spiral galaxies than the radio B -fields.

SALSA↗
NASA NTRS1997

Implications of Large Scale Shifts in Tropospheric NO(x) Levels in the Remote Tropical Pacific

A major observation recorded during NASA's western Pacific Exploratory Mission (PEM-West B) was the large shift in tropical NO levels as a function of geographical location. High-altitude NO levels exceeding 100 pptv were observed during portions of tropical flights 5-8, while values almost never exceeded 20 pptv during tropical flights 9 and 10. The geographical regions encompassing these two flight groupings are here labeled "high" and "low" NOx regimes. A comparison of these two regimes, based on back trajectories and chemical tracers, suggests that air parcels in both were strongly influenced by deep convection. The low NO(x) regime appears to have been predominantly impacted by marine convection, whereas the high NO(x) regime shows evidence of having been more influenced by deep convection over a continental land mass. DMSP satellite observations point strongly toward lightning as the major source of NOx in the latter regime. Photochemical ozone formation in the high NOx regime exceeded that for low NO(x) by factors of 2 to 6, whereas O3 destruction in the low NO(x)regime exceeded that for high NOx by factors of up to 3. Taking the tropopause height to be 17 km, estimates of the net photochemical effect on the O3 column revealed that the high NO(x) regime led to a small net production. By contrast, the low NOx regime was shown to destroy O3 at the rate of 3.4 % per day. One proposed mechanism for off-setting this projected large deficit would involve the transport of 03 rich midlatitude air into the tropics. Alternatively, it is suggested that O3 within the tropics may be overall near self-sustaining with respect to photochemical activity. This scenario would require that some tropical regions, unsampled at the time of PEM-B, display significant net column O3 production, leading to an overall balanced budget for the "greater" tropical Pacific basin. Details concerning the chemical nature of such regimes are discussed.

Crawford, J. H.↗
DOE OSTI2022

Experiments and Modeling of Proppant Embedment and Fracture Conductivity for the Caney Shale, Oklahoma, USA

ABSTRACT: The ultimate aim of hydraulic fracturing is to have a long and conductive flow path that extends from the wellbore into the formation. The effective fracture length is part of a hydraulically propped fracture which contributes to production. The difficulty in achieving economical production targets from shale reservoirs is at the forefront in many exploration companies. Fracture conductivity loss is related to; proppant embedment under depletion, proppant crushing, damage as a result of fracturing fluid, fines migration and proppant-pack permeability-damage are some of the factors that contribute to production decline after hydraulic fracturing in shale reservoirs. The Caney Shale is a calcareous organic-rich mudrock. Various studies have investigated the effect that clay on shale well productivity, however, there is currently no literature on the Caney shale in relation to horizontal wells; all available literature exists in vertical wells as well as on formations of the Caney that are shallow in comparison to an emerging play which is twice the depth. In this paper we investigate stress-dependent fracture conductivity of proppant-filled fractures and proppant embedment in Caney shale through laboratory and modeling studies. API fracture conductivity tests were conducted using 2% KCl on five locations within the Caney shale that consisted of selecting three brittle(reservoir) zones and two ductile zones. Confining pressures range from 1,000 psi to 12,000 psi at 210°F. Conductivity, permeability as well as embedment were measured during the test. Our experimental results have confirmed that improved fracture conductivity is attributed to; proppant size, the increase in porosity of the proppant pack, closure pressure changes and the reduction in fracture conductivity are a function of many factors such as fracture closure stress. The findings from this study could help the stimulation design by providing new insights into the critical factors that are to be determined to facilitate the choice of proppants as well as fracturing fluids for long term production and recovery from shale reservoirs. 1 INTRODUCTION The development of low permeability formations, like shales, has been aided by hydraulic fracturing of horizontal wells (Radonjic et al., 2020). Hydraulic fracturing fluid is injected at a high pressure to induce tensile fractures that can link to and stimulate natural fractures (Katende et al., 2021a,b). Preserving adequate conductivity in hydraulic fractures over the life of the wells is required for economic production; nevertheless, conserving such conductivity can be difficult in some circumstances, particularly in soft, clay-rich formations (Wang et al., 2021). Proppant particles help to keep the fractures open when the pumping stops and the fracturing fluid returns to the wellbore, producing one or more propped hydraulic fractures of varying length, breadth, and height (Katende et al., 2021a). The proppant pack within the hydraulic fracture boosts well output by providing a greater permeability flowpath for hydrocarbons (Duenckel et al., 2016). Proppant in the fracture is under complicated stress conditions, and the interplay between the rock formation and the proppant pack has a significant impact on proppant-pack permeability (Karazincir et al., 2019). Proppant may be embedded (Katende et al., 2021a) in the rock or crushed into small pieces if the proppant size and strength characteristics are not specified appropriately, resulting in a loss in proppant-pack permeability and fracture aperture, and consequently a fall in well output.

Katende, A.↗
NASA NTRS2003

Regional Variability in Convection and Rain Retrievals from the TRMM Microwave Imager (TMI)

Precipitation Radar (PR) on board the TRMM satellite shows that the average height of 30 dBz in convective rain areas of the tropics varies significantly from one region to the other. When the convection is weak this height is shallow and when convection is strong this height extends deeper into the troposphere. The brightness temperature (Tb) measured by the microwave radiometer by itself does not reflect this nature of convection satisfactorily. Radiative transfer simulations of Tbs reveal that this could be due to the variations in the vertical distribution of optically active water and ice hydrometeors and their density, shape, and size. These variations are not coupled uniquely to the strength of the convective updrafts, and as a result the Tbs do not reflect properly the convective strength indicated by PR. Because of this deficiency in the Tbs the rain rate deduced from them differs from that of PR. For this reason, to improve the estimation of rain rate we have developed an empirical method. In this method a parameter based on the areal extent of the Tbs that exceed a certain magnitude is included along with the Tbs. Rain rate deduced with this approach is better correlated with that of PR when compared to the current Version 5 operational algorithm. The percentage of rain volume as a function of rain rate, for a given region of 5deg lat. X 5deg long. over a period of three months, deduced from this method, is also in better agreement with that of the PR.

Prabhakara, C.↗
DOE OSTI2026

Room temperature sodium diffusion in sodium iron phosphate investigated with neutron scattering

Sodium iron phosphate, Na x FePO 4 , is currently one of the most promising cathode materials for sodium ion batteries. Here we employ a systematic approach using three neutron scattering techniques to investigate temperature dependent properties with a focus on the charge transfer processes. Two distinct Na + diffusive motions are observed with quasielastic neutron scattering, including room temperature motion. These processes correspond to rapid double jumps through extrinsic vacancy sites of ∼6 Å along with the more fundamental nearest neighbor jumps of 3.1 Å. Experimental observations of room temperature ionic diffusion are highly desirable, as the diffusion pathways are naturally well suited for efficient transport. Inelastic neutron scattering was used to probe the vibrational density of states that assist in the ionic and electronic (polaronic) transport processes, which are further supported with density functional theory. Neutron diffraction examines the disordering of the sodium sublattice upon transition to the solid-solution phase, which is believed to enhance the overall kinetics during sodium insertion/extraction. Finally, 57 Fe Mössbauer spectroscopy measurements probes polaron hopping and the local structural evolution that modifies the electron density surrounding the redox active iron sites. This study sheds light on the temperature induced atomic scale dynamics that occur in olivine Na x FePO 4 .

Novak, Eric [Swarthmore College, Swarthmore, PA (U↗