Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “data harmonization”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

Numerical experiments in geomagnetic modeling

Numerical tests were made, using least squares fitting of a spherical harmonic model, to a selection of Magsat data to determine the practical limits of this technique with modern computers. The resulting (M102189) model, whose coefficients were adjusted up to n = 50, was compared with M07AV6, a previous model which used least squares (on vector data) for coefficients up to n = 29, and Gauss-Legendre quadrature (on Z residuals) to adjust the coefficients up to n = 63. For the new least squares adjustment to n = 50 a condition number of 115 was obtained for the solution matrix, with a resulting precision of 11 significant figures. The M102189 model shows a lower and more Gaussian residual distribution than did M07AV6, though the Gaussian envelope fits to the residual distributions, even for the scalar field, gives "standard deviations' never lower than 6 nT, a factor of three higher than the estimated Magsat observational errors. Ionospheric currents are noted to have a significant effect on the coefficients of the internal potential functions.

Cain, Joseph C.↗
NASA NTRS2003

Measurement and Simulation of Spontaneous Raman Scattering Spectra in High-Pressure, Fuel-Rich H2-Air Flames

Rotational vibrational spontaneous Raman spectra (SRS) of H2, N2, and H2O have been measured in H2-air flames at pressures up to 30 atm as a first stem towards establishing a comprehensive Raman spectral database for temperatures and species in high-pressure combustion. A newly developed high-pressure burner facility provides steady, reproducible flames with a high degree of flow precision. We have obtained an initial set of measurements that indicate the spectra are of sufficient quality in terms of spectral resolution, wavelength coverage, and signal-to-noise ratio for use in future reference standards. The fully resolved Stokes and anti-Stokes shifted SRS spectra were collected in the visible wavelength range (400-700 nm) using pulse-stretched 532 nm excitation and a non-intensified CCD spectrograph with a high-speed shutter. Reasonable temperatures were determined via the intensity distribution of rotational H2 lines at stoichiometry and fuel-rich conditions. Theoretical Raman spectra of H2 were computed using a semi-classical harmonic-oscillator model with recent pressure broadening data and were compared with experimental results. The data and simulation indicated that high-J rotational lines of H2 might interfere with the N2 vibrational Q-branch lines, and this could lead to errors in N2-Raman thermometry based on the line-fitting method. From a comparison of N2 Q-branch spectra in lean H2 low-pressure (1.2 atm) and high-pressure (30 atm) flames, we found no significant line-narrowing or -broadening effects at the current spectrometer resolution of 0.04 nm.

Kojima, Jun↗
NASA NTRS1998

Multipath calibration in GPS pseudorange measurements

Novel techniques are disclosed for eliminating multipath errors, including mean bias errors, in pseudorange measurements made by conventional global positioning system receivers. By correlating the multipath signals of different satellites at their cross-over points in the sky, multipath mean bias errors are effectively eliminated. By then taking advantage of the geometrical dependence of multipath, a linear combination of spherical harmonics are fit to the satellite multipath data to create a hemispherical model of the multipath. This calibration model can then be used to compensate for multipath in subsequent measurements and thereby obtain GPS positioning to centimeter accuracy.

Kee, Changdon↗
DOE OSTI2023

Excitation of toroidally localized harmonics of global Alfvén eigenmodes

Abstract The spherical tokamak NSTX and the upgraded, higher toroidal field, NSTX-U, routinely create plasmas with large un-thermalized populations of super-Alfvénic fast ions. This population of non-thermal energetic ions excites a broad spectrum of Alfvénic waves, from the lower frequency (≈30 kHz–≈200 kHz) toroidal Alfvén eigenmodes (Cheng and Chance 1986 Phys. Fluids 29 3695), to the higher frequency (400 kHz–3 MHz) global Alfvén eigenmodes (GAEs) and compressional Alfvén eigenmodes (Goedbloed 1975 Phys. Fluids 18 1258). In this paper we present evidence that the GAE non-linearly excite modes, presumably GAE, at frequencies consistent with non-linear or 3-wave coupling. The observation of the excitation of 2 nd harmonic GAE through the intrinsic non-linearity of Global Alfvén modes demonstrates that the non-linear terms can act as an exciter-antenna inside the plasma, broadcasting at harmonics of the mode frequency and with concomitantly shorter wavelengths. As with experiments using external antenna to excite otherwise weakly stable Alfvénic modes, this data can provide information about the stability of modes at harmonics of the GAE. It may also provide information on the nature of the non-linearities in the wave dispersion equation or a direct measure of the mode amplitude (Smith et al 2006 Phys. Plasmas 13 042504). We also report that the short wavelength 2 nd harmonic GAE can be strongly toroidally localized. The observation that the shorter wavelength ( n ≈ 20) 2 nd harmonic modes can be toroidally localized potentially has implications for the impact of fast-ion driven instabilities on fast-ion confinement in ITER and future fusion reactors (Gorelenkov et al 2014 Nucl. Fusion 54 125001).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

First measurement of symmetric cumulants of hexagonal flow harmonics in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV

Correlations between event-by-event fluctuations of anisotropic flow harmonics are measured in Pb-Pb collisions at a center-of-mass energy per nucleon pair of 5.02 TeV, as recorded by the ALICE detector at the LHC. This study presents correlations up to the hexagonal flow harmonic 𝑣 6 , which was measured for the first time. The magnitudes of these higher-order correlations are found to vary as a function of collision centrality and harmonic order. These measurements are compared to viscous hydrodynamic model calculations with EKRT initial conditions and to the iEBE-VISHNU model with T R ⁢ENTo initial conditions. The observed discrepancies between the data and the model calculations vary depending on the harmonic combinations. Due to the sensitivity of model parameters estimated with Bayesian analyses to these higher-order observables, the results presented in this work provide new and independent constraints on the initial conditions and transport properties in theoretical models used to describe the system created in heavy-ion collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1993

Investigation of a geodesy coexperiment to the Gravity Probe B relativity gyroscope program

Geodesy is the science of measuring the gravitational field of and positions on the Earth. Estimation of the gravitational field via gravitation gradiometry, the measurement of variations in the direction and magnitude of gravitation with respect to position, is this dissertation's focus. Gravity Probe B (GP-B) is a Stanford satellite experiment in gravitational physics. GP-B will measure the precession the rotating Earth causes on the space time around it by observing the precessions of four gyroscopes in a circular, polar, drag-free orbit at 650 km altitude. The gyroscopes are nearly perfect niobium-coated spheres of quartz, operating at 1.8 K to permit observations with extremely low thermal noise. The permissible gyroscope drift rate is miniscule, so the torques on the gyros must be tiny. A drag-free control system, by canceling accelerations caused by nongravitational forces, minimizes the support forces and hence torques. The GP-B system offers two main possibilities for geodesy. One is as a drag-free satellite to be used in trajectory-based estimates of the Earth's gravity field. We described calculations involving that approach in our previous reports, including comparison of laser only, GPS only, and combined tracking and a preliminary estimate of the possibility of estimating relativistic effects on the orbit. The second possibility is gradiometry. This technique has received a more cursory examination in previous reports, so we concentrate on it here. We explore the feasibility of using the residual suspension forces centering the GP-B gyros as gradiometer signals for geodesy. The objective of this work is a statistical prediction of the formal uncertainty in an estimate of the Earth's gravitation field using data from GP-B. We perform an instrument analysis and apply two mathematical techniques to predict uncertainty. One is an analytical approach using a flat-Earth approximation to predict geopotential information quality as a function of spatial wavelength. The second estimates the covariance matrix arising in a least-squares estimate of a spherical harmonic representation of the geopotential using GP-B gradiometer data. The results show that the GP-B data set can be used to create a consistent estimate of the geopotential up to spherical harmonic degree and order 60. The formal uncertainty of all coefficients between degrees 5 and 50 is reduced by factors of up to 30 over current satellite-only estimates and up to 7 over estimates which include surface data. The primary conclusion resulting from this study is that the gravitation gradiometer geodesy coexperiment to GP-B is both feasible and attractive.

Everitt, C. W. F.↗
NASA NTRS2000

Test Cases for a Rectangular Supercritical Wing Undergoing Pitching Oscillations

Steady and unsteady measured pressures for a Rectangular Supercritical Wing (RSW) undergoing pitching oscillations have been presented. From the several hundred compiled data points, 27 static and 36 pitching oscillation cases have been proposed for computational Test Cases to illustrate the trends with Mach number, reduced frequency, and angle of attack. The wing was designed to be a simple configuration for Computational Fluid Dynamics (CFD) comparisons. The wing had an unswept rectangular planform plus a tip of revolution, a panel aspect ratio of 2.0, a twelve per cent thick supercritical airfoil section, and no twist. The model was tested over a wide range of Mach numbers, from 0.27 to 0.90, corresponding to low subsonic flows up to strong transonic flows. The higher Mach numbers are well beyond the design Mach number such as might be required for flutter verification beyond cruise conditions. The pitching oscillations covered a broad range of reduced frequencies. Some early calculations for this wing are given for lifting pressure as calculated from a linear lifting surface program and from a transonic small perturbation program. The unsteady results were given primarily for a mild transonic condition at M = 0.70. For these cases the agreement with the data was only fair, possibly resulting from the omission of viscous effects. Supercritical airfoil sections are known to be sensitive to viscous effects (for example, one case cited). Calculations using a higher level code with the full potential equations have been presented for one of the same cases, and with the Euler equations. The agreement around the leading edge was improved, but overall the agreement was not completely satisfactory. Typically for low-aspect-ratio rectangular wings, transonic shock waves on the wing tend to sweep forward from root to tip such that there are strong three-dimensional effects. It might also be noted that for most of the test, the model was tested with free transition, but a few points were taken with an added transition strip for comparison. Some unpublished results of a rigid wing of the same airfoil and planform that was tested on the pitch and plunge apparatus mount system (PAPA) showed effects of the lower surface transition Strip on flutter at the lower subsonic Mach numbers. Significant effects of a transition strip were also obtained on a wing with a thicker supercritical section on the PAPA mount system. Both of these flutter tests on the PAPA resulted in very low reduced frequencies that may be a factor in this influence of the transition strip. However, these results indicate that correlation studies for RSW may require some attention to the estimation of transition location to accurately treat viscous effects. In this report several Test Cases are selected to illustrate trends for a variety of different conditions with emphasis on transonic flow effects. An overview of the model and tests is given and the standard formulary for these data is listed. Sample data points are presented in both tabular and graphical form. A complete tabulation and plotting of all the Test Cases is given. Only the static pressures and the real and imaginary parts of the first harmonic of the unsteady pressures are available. All the data for the test are available in electronic file form. The Test Cases are also available as separate electronic files.

Bennett, Robert M.↗
NASA NTRS1985

Gravity Field of Venus and Comparison with Earth

The acceleration (gravity) anomaly estimates by spacecraft tracking, determined from Doppler residuals, are components of the gravity field directed along the spacecraft Earth line of sight (LOS). These data constitute a set of vector components of a planet's gravity field, the specific component depending upon where the Earth happened to be at the time of each measurement, and they are at varying altitudes above the planet surface. From this data set the gravity field vector components were derived using the method of harmonic splines which imposes a smoothness criterion to select a gravity model compatible with the LOS data. Given the piecewise model it is now possible to upward and downward continue the field quantities desired with a few parameters unlike some other methods which must return to the full dataset for each desired calculation.

Bowin, C.↗
DOE OSTI2023

Optimal Frequency-domain Analysis for Spacecraft Time Series: Introducing the Missing-data Multitaper Power Spectrum Estimator

While the Lomb–Scargle periodogram is foundational to astronomy, it has a significant shortcoming: the variance in the estimated power spectrum does not decrease as more data are acquired. Statisticians have a 60 yr history of developing variance-suppressing power spectrum estimators, but most are not used in astronomy because they are formulated for time series with uniform observing cadence and without seasonal or daily gaps. Here we demonstrate how to apply the missing-data multitaper power spectrum estimator to spacecraft data with uniform time intervals between observations but missing data during thruster fires or momentum dumps. The F-test for harmonic components may be applied to multitaper power spectrum estimates to identify statistically significant oscillations that would not rise above a white noise–based false alarm probability. Multitapering improves the dynamic range of the power spectrum estimate and suppresses spectral window artifacts. We show that the multitaper–F-test combination applied to Kepler observations of KIC 6102338 detects differential rotation without requiring iterative sinusoid fitting and subtraction. Significant signals reside at harmonics of both fundamental rotation frequencies and suggest an antisolar rotation profile. Next we use the missing-data multitaper power spectrum estimator to identify the oscillation modes responsible for the complex "scallop-shell" shape of the K2 light curve of EPIC 203354381. We argue that multitaper power spectrum estimators should be used for all time series with regular observing cadence.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1985

Venus topography - A harmonic analysis

A model of Venusian global topography has been obtained by fitting an eighteenth-degree harmonic series to Pioneer Venus orbiter radar altimeter data. The mean radius is (6051.45 + or - 0.04) km. The corresponding mean density is (5244.8 + or 0.5) kg/cu m. The center of figure is displaced from the center of mass by (0.339 + or - 0.088) km towards (6.6 + or 10.1) deg N, (148. 8 + or - 7.7) deg. The figure of Venus is distinctly triaxial, but the orientation and magnitudes of the principal topographic axes correlate rather poorly with the gravitational principal axes. However, the higher-degree harmonics of topography and gravity are significantly correlated. The topographic variance spectrum of Venus is very similar in form to those of the moon, Mars, and especially earth. It is suggested that this spectral similarity simply reflects a statistical balance between constructional and degradational geomorphic proceses. Venus and earth are particularly similar (and differ from the moon and Mars) in that the larger bodies both exhibit a significant low degree deficit (relative to the extrapolated trend of the higher harmonics).

Bills, B. G.↗
NASA NTRS2004

Vibration-Based Data Used to Detect Cracks in Rotating Disks

Rotor health monitoring and online damage detection are increasingly gaining the interest of aircraft engine manufacturers. This is primarily due to the fact that there is a necessity for improved safety during operation as well as a need for lower maintenance costs. Applied techniques for the damage detection and health monitoring of rotors are essential for engine safety, reliability, and life prediction. Recently, the United States set the ambitious goal of reducing the fatal accident rate for commercial aviation by 80 percent within 10 years. In turn, NASA, in collaboration with the Federal Aviation Administration, other Federal agencies, universities, and the airline and aircraft industries, responded by developing the Aviation Safety Program. This program provides research and technology products needed to help the aerospace industry achieve their aviation safety goal. The Nondestructive Evaluation (NDE) Group of the Optical Instrumentation Technology Branch at the NASA Glenn Research Center is currently developing propulsion-system-specific technologies to detect damage prior to catastrophe under the propulsion health management task. Currently, the NDE group is assessing the feasibility of utilizing real-time vibration data to detect cracks in turbine disks. The data are obtained from radial blade-tip clearance and shaft-clearance measurements made using capacitive or eddy-current probes. The concept is based on the fact that disk cracks distort the strain field within the component. This, in turn, causes a small deformation in the disk's geometry as well as a possible change in the system's center of mass. The geometric change and the center of mass shift can be indirectly characterized by monitoring the amplitude and phase of the first harmonic (i.e., the 1 component) of the vibration data. Spin pit experiments and full-scale engine tests have been conducted while monitoring for crack growth with this detection methodology. Even so, published data are extremely limited, and the basic foundation of the methodology has not been fully studied. The NDE group is working on developing this foundation on the basis of theoretical modeling as well as experimental data by using the newly constructed subscale spin system shown in the preceding photograph. This, in turn, involved designing an optimal sub-scale disk that was meant to represent a full-scale turbine disk; conducting finite element analyses of undamaged and damaged disks to define the disk's deformation and the resulting shift in center of mass; and creating a rotordynamic model of the complete disk and shaft assembly to confirm operation beyond the first critical concerning the subscale experimental setup. The finite element analysis data, defining the center of mass shift due to disk damage, are shown. As an example, the change in the center of mass for a disk spinning at 8000 rpm with a 0.963-in. notch was 1.3 x 10(exp -4) in. The actual vibration response of an undamaged disk as well as the theoretical response of a cracked disk is shown. Experiments with cracked disks are continuing, and new approaches for analyzing the captured vibration data are being developed to better detect damage in a rotor. In addition, the subscale spin system is being used to test the durability and sensitivity of new NDE sensors that focus on detecting localized damage. This is designed to supplement the global response of the crack-detection methodology described here.

Gyekenyesi, Andrew L.↗
NASA NTRS2000

Test Cases for Flutter of the Benchmark Models Rectangular Wings on the Pitch and Plunge Apparatus

The supercritical airfoil was chosen as a relatively modem airfoil for comparison. The BOO12 model was tested first. Three different types of flutter instability boundaries were encountered, a classical flutter boundary, a transonic stall flutter boundary at angle of attack, and a plunge instability near M = 0.9 and for zero angle of attack. This test was made in air and was Transonic Dynamics Tunnel (TDT) Test 468. The BSCW model (for Benchmark SuperCritical Wing) was tested next as TDT Test 470. It was tested using both with air and a heavy gas, R-12, as a test medium. The effect of a transition strip on flutter was evaluated in air. The B64AOlO model was subsequently tested as TDT Test 493. Some further analysis of the experimental data for the BOO12 wing is presented. Transonic calculations using the parameters for the BOO12 wing in a two-dimensional typical section flutter analysis are given. These data are supplemented with data from the Benchmark Active Controls Technology model (BACT) given and in the next chapter of this document. The BACT model was of the same planform and airfoil as the BOO12 model, but with spoilers and a trailing edge control. It was tested in the heavy gas R-12, and was instrumented mostly at the 60 per cent span. The flutter data obtained on PAPA and the static aerodynamic test cases from BACT serve as additional data for the BOO12 model. All three types of flutter are included in the BACT Test Cases. In this report several test cases are selected to illustrate trends for a variety of different conditions with emphasis on transonic flutter. Cases are selected for classical and stall flutter for the BSCW model, for classical and plunge for the B64AOlO model, and for classical flutter for the BOO12 model. Test Cases are also presented for BSCW for static angles of attack. Only the mean pressures and the real and imaginary parts of the first harmonic of the pressures are included in the data for the test cases, but digitized time histories have been archived. The data for the test cases are available as separate electronic files. An overview of the model and tests is given, the standard formulary for these data is listed, and some sample results are presented.

Robert M Bennett↗
NASA NTRS2003

Characterization of Viscoelastic Properties of Polymeric Materials Through Nanoindentation

Nanoindentation is used to determine the dynamic viscoelastic properties of six polymer materials. It is shown that varying the harmonic frequency of the nanoindentation does not have any significant effect on the measured storage and loss moduli of the polymers. Agreement is found between these results and data from DMA testing of the same materials. Varying the harmonic amplitude of the nanoindentation does not have a significant effect on the measured properties of the high performance resins, however, the storage modulus of the polyethylene decreases as the harmonic amplitude increases. Measured storage and loss moduli are also shown to depend on the density of the polyethylene.

Odegard, G. M.↗
NASA NTRS1994

Isostatic compensation of equatorial highlands on Venus

Spherical harmonic models for Venus' global topography and gravity incorporating Magellan data are used to test isostatic compensation models in five 30 deg x 30 deg regions representative of the main classes of equatorial highlands. The power spectral density for the harmonic models obeys a power-law scaling with spectral slope Beta approximately 2 (Brown noise) for the topography and Beta approximately 3 (Kaula's law) for the geoid, similar to what is observed for Earth. The Venus topography spectrum has lower amplitudes than Earth's which reflects the dominant lowland topography on Venus. Observed degree geoid to topography ratios (GTRs) on Venus are significantly smaller than degree GTRs for uncompensated topography, indicative of substantial compensation. Assuming a global Airy compensation, most of the topography is compensated at depths greater than 100 km, suggesting a thick lithosphere on Venus. For each region considered we obtain a regional degree of compensation C from a linear regression of Bouguer anomaly versus Bouguer gravity data. Geoid anomaly (N) versus topography variation (h) data for each sample were compared, in the least-squares sense, to theoretical correlations for Pratt, Airy, and thermal thinning isostasy models yielding regional GTR, zero-elevation crustal thickness (H), and zero elevation thermal lithosphere thickness (y(sub L(sub 0)), respectively. We find the regional compensation to be substantial (C approximately 52-80%), and the h, N data correlations in the chosen areas can be explained by isostasy models applicable on the Earth and involving variations in crustal thickness (Airy) and/or lithospheric (thermal thinning) thickness. However, a thick crust and lithosphere (y(sub L(sub 0)) approximately 300 km) must be assumed for Venus.

Kucinskas, Algis B.↗
NASA NTRS1990

State-of-the-art of high-speed propeller noise prediction - A multidisciplinary approach and comparison with measured data

The results of NASA's Propeller Test Assessment program involving extensive flight tests of a large-scale advanced propeller are presented. This has provided the opportunity to evaluate the current capability of advanced propeller noise prediction utilizing principally the exterior acoustic measurements for the prediction of exterior noise. The principal object of this study was to evaluate the state-of-the-art of noise prediction for advanced propellers utilizing the best available codes of the disciplines involved. The effects of blade deformation on the aerodynamics and noise of advanced propellers were also studied. It is concluded that blade deformation can appreciably influence propeller noise and aerodynamics, and that, in general, centrifugal and blade forces must both be included in the calculation of blade forces. It is noted that the present capability for free-field noise prediction of the first three harmonics for advanced propellers is fairly good. Detailed data and diagrams of the test results are presented.

Dunn, Mark H.↗
DOE OSTI2024

Data for Grogan et al. "Bringing Hydrologic Realism to Water Markets"

This data set provides model output and post-processing files required to reproduce the results, tables, and figures in the paper "Bringing Hydrologic Realism to Water Markets" by Grogan et al. (in review). Other input data used in this study includes: Lisk, M., Grogan, D., Zuidema, S., Caccese, R., Peklak, D., Zheng, J., Fisher-Vanden, K., Lammers, R., Olmstead, S., & Fowler, L. (2023). Harmonized Database of Western U.S. Water Rights (HarDWR) (Version v1) [Data set]. MSD-LIVE Data Repository. https://doi.org/10.57931/2205619 Two models were used in this study: (1) The University of New Hampshire Water Balance Model WBM, and (2) a Water Market Model. Market model code and model output post-processing code that make use of these data can be found here Model output files are: 1. WBM output files: scenario[x]_wbm_output.zip Where [x] is one of 1, 2, 2a, 3, and 3a Each zipped directory contains 7 gridded NetCDF files, each reporting the 10-year annual average value of a given variable, in units of average mm/day: File Name: wbm_indUseGross_yc.nc; Description: Water withdrawals by industry (part of the urban sector) File Name: wbm_domUseGross_yc.nc; Description: Water withdrawals by the domestic sector (part of the urban sector) File Name: wbm_irrigationGross_yc.nc; Description: Water withdrawals for agriculture File Name: wbm_irrigationExtra_yc.nc; Description: Water withdrawals from unsustainable groundwater for agriculture File Name: wbm_indUseEvap_yc.nc; Description: Consumptive water use by industry File Name: wbm_domUseEvap_yc.nc; Description: Consumptive water use by the domestic sector File Name: wbm_irrigationNet_yc.nc; Description: Consumptive water use by agriculture The file full_cell_area.nc gives the area of each grid cell in km2, which is used for converting water depth to water volume. 2. Water market model output & post processing output Folder: marketTrdSummaries/ Description: Files in this folder are used as input to code 1_WelfareCalculation_actual_trades.R. They summarize historical water right trade transactions in each state. File Name: welfare_gain_by_state_sector.csv; Description: Welfare gains by state and sector, as shown in Figure 3F. Used in code Figure3.R and produced (as a .xlsx file) by code 2_DemandCurves_simulated_trades.R File Name: welfare_data_actual.rdata; Description: welfare gains by WMA from actual historical trades, as shown in Figure 3A. This data is the output of code 1_WelfareCalculation_actual_trades.R File Name: welfare_summary_simulated.xlsx; Description: Welfare gains by state as simulated by the market model in Scenario 1. Produced by code 2_DemandCurves_simulated_trades.R, and used in code 4_WelfareCalculation.R. File Name: welfare_summary_cutoffs.xlsx; Description: Welfare gains by state as simulated by the market model in Scenario 2. Produced by code 3_DemandCurves_simulated_trades_cutoffs.R, and used in code 4_WelfareCalculation.R. File Name: welfare_summary_cutoffs_SGMS.xlsx; Description: Welfare gains by state as simulated by the market model in Scenario 2a. Produced by code 3_DemandCurves_simulated_trades_cutoffs.R, and used in code 4_WelfareCalculation.R. File Name: welfare_data_actual.rdata; Description: Spatial data, actual historical welfare gains by WMA as shown in Figure 3A. Produced by code 4_WelfareCalculation.R and used by code Figure3.R. File Name: welfare_data_simulated.rdata; Description: Spatial data, simulated Scenario 1 welfare gains by WMA as shown in Figure 3B. Produced by code 4_WelfareCalculation.R and used by code Figure3.R. File Name: welfare_data_simulated_cutoffs.rdata; Description: Spatial data, simulated Scenario 2 welfare gains by WMA. Produced by code 4_WelfareCalculation.R and used by code Figure3.R. File Name: welfare_data_simulated_cutoffs_SGMA.rdata; Description: Spatial data, simulated Scenario 2a welfare gains by WMA. Produced by code 4_WelfareCalculation.R and used by code Figure3.R. Additional files are provided for efficient reproduction of tables and figures. These include: File Name: wma_thresold_dates_Scenario2(a).csv; Description: Wet vs. paper right threshold dates for each WMA. Shown in Figure 2A,B. Produced and used by code calculate_thresolds_Figure2.R File Name: WWRTradeBounds (directory); Description: Trade boundary shapefile required to reproduce Figure 3A-D. Used in code Figure3.R File Name: welfare_region_totals.csv; Description: Welfare gains for the entire study region, as shown in Figure 3E. Used in code Figure3.R File Name: Welfare_gain_by_state_sector.csv; Description: Welfare gains by state and sector, as shown in Figure 3F. Used in code Figure3.R and produced (as a .xlsx file) by code 2_DemandCurves_simulated_trades.R File Name: WECC_MERIT_5min_v3b_mask.nc; Description: Gridded file that identified which land grid cells are in the WBM model domain, used for processing in code Figure4.py File Name: Table_1.csv; Description: All data in Table 1, reproducible from WBM output files using code table_1.R

Economics↗
NASA NTRS1976

A possible experiment with two counter-orbiting drag-free satellites to obtain a new test of Einstein's general theory of relativity and improved measurements in geodesy

In 1918, Lense and Thirring calculated that a moon in orbit around a massive rotating planet would experience a nodal dragging effect due to general relativity. We describe an experiment to measure this effect by means of two counter-orbiting drag-free satellites in polar orbit about the earth. For a 2-1/2 year experiment, the measurement should approach an accuracy of 1%. An independent measurement of the geodetic precession of the orbit plane due to the motion about the sun may also be possible to about 10% accuracy. In addition to precision tracking data from existing ground stations, satellite-to-satellite Doppler data are taken at points of passing near the poles to yield an accurate measurement of the separation distance between the two satellites. New geophysical information on both earth harmonics and tidal effects is inherent in this polar ranging data.

Van Patten, R. A.↗