Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Data fitting”

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,315 records · Page 73

A Synthesis of VIIRS Solar and Lunar Calibrations

The NASA VIIRS Ocean Science Team (VOST) has developed two independent calibrations of the SNPP VIIRS moderate resolution reflective solar bands using solar diffuser and lunar observations through June 2013. Fits to the solar calibration time series show mean residuals per band of 0.078-0.10%. There are apparent residual lunar libration correlations in the lunar calibration time series that are not accounted for by the ROLO photometric model of the Moon. Fits to the lunar time series that account for residual librations show mean residuals per band of 0.071-0.17%. Comparison of the solar and lunar time series shows that the relative differences in the two calibrations are 0.12-0.31%. Relative uncertainties in the VIIRS solar and lunar calibration time series are comparable to those achieved for SeaWiFS, Aqua MODIS, and Terra MODIS. Intercomparison of the VIIRS lunar time series with those from SeaWiFS, Aqua MODIS, and Terra MODIS shows that the scatter in the VIIRS lunar observations is consistent with that observed for the heritage instruments. Based on these analyses, the VOST has derived a calibration lookup table for VIIRS ocean color data based on fits to the solar calibration time series.

Calibrations↗

Specular reflectance of optical-black coatings in the far infrared

Far-infrared specular reflectance spectra of seven optically black coatings near normal incidence are presented. Seven photometric spectra were obtained using eleven bandpass transmission filters in the wavelength range between 12 and 500 microns, and three interferometric spectra were obtained for corroboration. Data on the construction, thickness, and rms surface roughness of the coatings are also presented. The chemical composition of three coatings can be distinguished from that of the others by a strong absorption feature between 20 and 40 microns, which can be largely attributed to amorphous silicate material. At 100 microns, the most and least reflective coatings differ by nearly 3 orders of magnitude. Inverse relationships observed between the spectra and the roughness and thickness of the coatings led to development of a reflecting-layer model for the measured reflectance. The model successfully describes the spectra at wavelengths outside the silicate absorption, and optical constants are deduced from a nonlinear least squares fit to the data. Parametric errors are estimated by chi-square analysis, and sensitivity tests are performed to determine which parameters control reflectance in different spectral regions.

Smith, S. M.↗

Properties of dust in the Martian atmosphere and its effect on temperature structure

Data from the Viking lander imaging experiment and the IRIS experiment on board the Mariner 9 Orbiter are used to develop a one-dimensional radiative-convective model to describe the Mars temperature profile. Photography was employed to determine the optical depth, mean size, shape, visible absorption coefficient, and vertical distribution of aerosols in the Mars atmosphere. The calculation procedures are outlined, particularly for the particle size distribution function, which indicated a cross-sectional average particle radius of 2.5 microns. Dust-free conditions in the atmosphere near the equator and at lower altitudes permit full surface absorption of sunlight and subsequent adiabatic lapse rates. Difficulties were encountered in accurately modeling temperature profiles to fit Viking lander data, and resultant variations are attributed to effects of both suspended dust particles and large scale circulation.

Pollack, J. B.↗

Improvements in Modeling 90 degree Bleed Holes for Supersonic Inlets

The modeling of porous bleed regions as boundary conditions in computational fluid dynamics (CFD) simulations of supersonic inlet flows has been improved through a scaling of sonic flow coefficient data for 90deg bleed holes. The scaling removed the Mach number as a factor in computing the sonic flow coefficient and allowed the data to be fitted with a quadratic equation, with the only factor being the ratio of the plenum static pressure to the surface static pressure. The implementation of the bleed model into the Wind-US CFD flow solver was simplified by no longer requiring the evaluation of the flow properties at the boundary-layer edge. The quadratic equation can be extrapolated to allow the modeling of small amounts of blowing, which can exist when recirculation of the bleed flow occurs within the bleed region. The improved accuracy of the bleed model was demonstrated through CFD simulations of bleed regions on a flat plate in supersonic flow with and without an impinging oblique shock. The bleed model demonstrated good agreement with experimental data and three-dimensional CFD simulations of bleed holes.

Slater, John W.↗

Thermal remote sensing of surface soil water content with partial vegetation cover for incorporation into climate models

This study outlines a method for the estimation of regional patterns of surface moisture availability (M(sub 0)) and fractional vegetation (Fr) in the presence of spatially variable vegetation cover. The method requires relating variations in satellite-derived (NOAA, Advanced Very High Resolution Radiometer (AVHRR)) surface radiant temperature to a vegetation index (computed from satellite visible and near-infrared data) while coupling this association to an inverse modeling scheme. More than merely furnishing surface soil moisture values, the method constitues a new conceptual and practical approach for combining thermal infrared and vegetation index measurements for incorporating the derived values of M(sub 0) into hydrologic and atmospheric prediction models. Application of the technique is demonstrated for a region in and around the city of Newcastle upon Tyne situated in the northeast of England. A regional estimate of M(sub 0) is derived and is probabbly good for fractional vegetation cover up to 80% before errors in the estimated soil water content become unacceptably large. Moreover, a normalization scheme is suggested from which a nomogram, `universal triangle,' is constructed and is seen to fit the observed data well. The universal triangle also simplifies the inclusion of remotely derived M(sub 0) in hydrology and meteorological models and is perhaps a practicable step toward integrating derived data from satellite measurements in weather forecasting.

Gillies, Robert R.↗

Observation of double J / ψ meson production in p Pb collisions at s NN = 8.16 TeV

The first observation of the concurrent production of two J / ψ mesons in proton-nucleus collisions is presented. The analysis is based on a proton-lead ( p Pb ) data sample recorded at a nucleon-nucleon center-of-mass energy of 8.16 TeV by the CMS experiment at the CERN LHC and corresponding to an integrated luminosity of 174.6 nb − 1 . The two J / ψ mesons are reconstructed in their μ + μ − decay channels with transverse momenta p T > 6.5 GeV and rapidity | y | < 2.4 . Events where one of the J / ψ mesons is reconstructed in the dielectron channel are also considered in the search. The p Pb → J / ψ J / ψ + X process is observed with a significance of 5.3 standard deviations. The measured inclusive fiducial cross section, using the four-muon channel alone, is σ ( p Pb → J / ψ J / ψ + X ) = 22.0 ± 8.9 ( stat ) ± 1.5 ( syst ) nb . A fit of the data to the expected rapidity separation for pairs of J / ψ mesons produced in single (SPS) and double (DPS) parton scatterings yields σ SPS p Pb → J / ψ J / ψ + X = 16.5 ± 10.8 ( stat ) ± 0.1 ( syst ) nb and σ DPS p Pb → J / ψ J / ψ + X = 5.4 ± 6.2 ( stat ) ± 0.4 ( syst ) nb , respectively. This latter result can be transformed into a lower bound on the effective DPS cross section, closely related to the squared average interparton transverse separation in the collision, of σ eff > 1.0 mb at 95% confidence level. © 2024 CERN, for the CMS Collaboration 2024 CERN

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A determination of the composition of the Saturnian stratosphere using the IUE

Ultraviolet spectra of Saturn from the IUE satellite was reduced to produce a geometric albedo of the planet from 1500 to 3000 A. By matching computer models to the albedo a chemical composition consistent with the data was determined. This model includes C2H2 and C2H6 with mixing ratios and distributions of 9 + or - 3 x 10 to the -8th in the top 20 mbar of the atmosphere with none below for C2H2 and 6 + or 1 x 10 to the -6th also in the top 20 mbar with none below for C2H6. The C2H2 and C2H6 distributions and the C2H6 mixing ratio are taken directly from the Voyager IRIS model (R. Courtin et al., 1981). The Voyager IRIS model also includes PH3, which is not consistent with the UV albedo from 1800 to 2400 A. This model requires a previously unidentified absorber to explain the albedo near 1600 A. After considering several candidates, it is found that the best fit to the data is obtained with H2O, having a column density of 6 + or - 1 x 10 to the -3 cm-am.

Winkelstein, P.↗

Data Analysis & Statistical Methods for Command File Errors

This paper explains current work on modeling for managing the risk of command file errors. It is focused on analyzing actual data from a JPL spaceflight mission to build models for evaluating and predicting error rates as a function of several key variables. We constructed a rich dataset by considering the number of errors, the number of files radiated, including the number commands and blocks in each file, as well as subjective estimates of workload and operational novelty. We have assessed these data using different curve fitting and distribution fitting techniques, such as multiple regression analysis, and maximum likelihood estimation to see how much of the variability in the error rates can be explained with these. We have also used goodness of fit testing strategies and principal component analysis to further assess our data. Finally, we constructed a model of expected error rates based on the what these statistics bore out as critical drivers to the error rate. This model allows project management to evaluate the error rate against a theoretically expected rate as well as anticipate future error rates.

Correlation Analysis↗

Parameters Inference and Model Reduction for the Single-Particle Model of Li Ion Cells

The Single-Particle Model (SPM) of Li ion cell is a computationally efficient model for simulating Li ion cell for weak to moderate currents. The model depends n a number of parameters describing the geometry and material properties of a cell components. In order to apply the model to simulating a cell, the best-fit parametric values have to be inferred from a constant discharge data. We report our efforts to determine the best-fit set for 18650 LP batteries. We found that rather than being best-fit by a particular point in the parametric space the data is fit equally well by an ensemble of points clustering about an effective multidimensional manifold in the parametric space. This property of the SPM is known to be shared by a multitude of the so-called "sloppy models" of complex systems, characterized by a few stiff directions in the parametric space, in which the predicted behavior varies significantly, and a number of sloppy directions in which the behavior doesn't change appreciably. Only the stiff parameters combinations are identifiable. Geometrical features of the BFM give insights to possible reduction of the SPM to a model having fewer sloppy parameters. We have constructed a hierarchy of such models. The fully reduced model depends on only stiff effective parameters which are identifiable and can be used for characterization of the battery's state of health.

Khasin, Michael↗

Projecting Power Converter Specific Power Through 2050 for Aerospace Applications

In order to analyze the potential fuel burn benefit from the electrification of aircraft powertrains, it is important to quantify the amount of weight that will be added to the aircraft for each additional component of the electric powertrain. This paper provides a projection of the specific power and efficiency of power converters, (AC-DC, DC-AC, or DC-DC), through the year 2050. Data was first collected on state of the art power converters in multiple application areas, creating a power converter database. Relevant specific powers were added to a set of historical data from 1976-2020, and then three different logistic curves were fit through the historical data to represent S-curve shaped growth through the year 2050. The three curves were differentiated by conservative, nominal, and aggressive assumptions for the year in which the logistic curve begins to bend down towards slower growth. With a 30% knockdown factor accounting for the additional weight required for a high altitude converter, projections range from the aggressive specific power projection of 52.9 kW/kg in 2050 to a much more conservative specific power of 12 kW/kg in which growth is limited due to certifiability concerns. Little historical data was found on converter efficiencies to project efficiency based on historical trends. Projections are based on expert opinion on yearly decreases in converter losses. 2050 projections range from 0.987 to 0.997.

fuel burn↗

Public Reference Data for Megawatt-Scale Hydrogen Electrolysis - NLR Historical Wind

The U.S. Department of Energy and the National Laboratory of the Rockies (NLR) demonstrate hydrogen electrolysis from variable sources, hydrogen compression and storage, and hydrogen fuel cell power production using megawatt-scale equipment at NLR’s Flatirons Campus as part of the Advanced Research on Integrated Energy Systems (ARIES) initiative. This dataset represents part of that effort and is intended for academic, national laboratory, industrial, and other stakeholders to plan, design, and validate models of megawatt-scale hydrogen technologies and diverse energy infrastructure nationwide. These data provide a baseline for how existing hydrogen electrolysis technologies perform when coupled with various energy technologies. Future datasets will demonstrate how existing hydrogen fuel cell technologies can provide controllable, dispatchable, and variable power output for artificial intelligence (AI) data centers and other variable loads. This dataset entry describes hydrogen production by conducting a statistical analysis of historical wind data over a five-year period (2020-2025) from a single 1.5MW turbine manufactured by General Electric (GE) located at NLR’s Flatirons Campus, to generate an experimental test profile that was deployed on a 1.25-MW proton exchange membrane type MC250 electrolyzer system manufactured by Nel Hydrogen . [1] While the electrolyzer balance-of-plant supports up to 2.5 MW of electrolysis, NLR only has a single 1.25-MW electrolysis stack. The historical wind data provided several metrics, however, the analysis particularly focused on the measured power output by the wind turbine. The power output time series of data for each day was categorized by total energy generation and standard deviation, and the day that represented the highest combination of these two metrics was chosen – December 25th, 2022. This process was then repeated for a moving four-hour window within this day to identify the most statistically variable period. Finally, this four-hour period was scaled by 65% to match the 1.25 MW electrolyzer. The electrolysis system controls hydrogen production by varying DC current applied to the stack, from a maximum of 3000 A to a minimum safe operation of 300 A, or 10%. Because the current – voltage characteristic changes as the stack ages and efficiency degrades, the actual minimum safe operating power changes over time. The historical wind profiles were translated from power (kilowatts) to current (amperes) using a curve fit with calibration data and sent to the electrolyzer power supply at 1 Hz frequency. For more details on the statistical analysis process, see the presentation labeled “ Public Reference Data for Megawatt-Scale Hydrogen Electrolysis” provided with each data entry. These datasets report relevant hydrogen balance-of-plant and system data, all captured at 1 Hz, including hydrogen mass production measured with an Emerson Coriolis flow meter. Each .zip file represents a single wind turbine electrolysis experiment and is formatted as follows: {technology}_{scaling factor}-{electrolyzer ramp rate in amperes/second} For instance, “wind-GE1.5MW_0.65-400.zip” represents the hour-long experiment using historical data from the wind-GE1.5MW turbine, scaled to 65%, with the electrolyzer power supply set to a maximum ramp rate (gain and slew) of 400 A/s. Each .zip folder contains the following files: A .csv file containing raw data An .xlsx file explaining all the fields in the raw data. A .png plot showing the time series of hydrogen production, electrolysis power consumption, and wind power input. A PDF file detailing the historical wind data statistical analysis used to generate the wind profile. An experiment labeled “characterization_200.zip” demonstrates the MC250 electrolyzer steady-state response with 30-minute load steps for a total duration of 5 hours. Finally, a .csv file is provided with all simulated wind experiments combined into one dataset labeled "combined_historical_wind_experiments.csv". NLR also built an AI/machine-learning predictive model based on these datasets. The model ingests the electrolyzer current command in amperes, as well as various pressures and temperatures across the system, and predicts hydrogen output in kilograms per hour. The complete model can be found at https://huggingface.co/NatLabRockies/ptmelt-hydrogen-electrolysis [1] nelhydrogen.com/product/mc-series-electrolyser .

08 HYDROGEN↗

Application of Support Vector Regression to Derive Crater Depth/Diameter From Satellite Images

Through the study of impact crater shapes, one can draw important conclusions about the nature and evolution of planetary surfaces [e.g., 1-4].In particular, studying the depth (d) to diameter (D)ratio (d/D) of a population of impact craters, in combination with crater count statistics, can yield valuable insights regarding rates of erosion and burial[5]. Motivated by the great abundance of available planetary surface image data, the goal of this project is to develop an efficient way to estimate d/D from satellite images of impact craters for which stereo information is not available [6]. We set out to develop and train a machine learning algorithm to extract d/D from a dataset of synthetic impact crater images for which model d/D is known. The applications of machine learning to planetary science are numerous and diverse [7], including automatic planetary surface mapping [8] and the detection of impact craters [9]. Our algorithm makes use of Support Vector Regression (SVR), which is a type of Support Vector Machine (SVM) [10, 11].SVMs are a branch of supervised machine learning valued for their straightforward implementation and versatility in solving both classification and regression problems. In regression analysis, an SVR algorithm produces a hyperplane function to fit the training data points, as well as an ε-tube that surrounds the hyperplane. Tunable hyperparameters include the width of the ε-tube (ε) and the amount an algorithm is penalized for points which fall outside the ε-tube.

L R Chin↗

Response Surface Modeling Tolerance and Inference Error Risk Specifications: Proposed Industry Standards

This paper reviews the derivation of an equation for scaling response surface modeling experiments. The equation represents the smallest number of data points required to fit a linear regression polynomial so as to achieve certain specified model adequacy criteria. Specific criteria are proposed which simplify an otherwise rather complex equation, generating a practical rule of thumb for the minimum volume of data required to adequately fit a polynomial with a specified number of terms in the model. This equation and the simplified rule of thumb it produces can be applied to minimize the cost of wind tunnel testing.

DeLoach, Richard↗

Retrievals of Profiles of Fine And Coarse Aerosols Using Lidar And Radiometric Space Measurements

In couple of years we expect the launch of the CALIPSO lidar spaceborne mission designed to observe aerosols and clouds. CALIPSO will collect profiles of the lidar attenuated backscattering coefficients in two spectral wavelengths (0.53 and 1.06 microns). Observations are provided along the track of the satellite around the globe from pole to pole. The attenuated backscattering coefficients are sensitive to the vertical distribution of aerosol particles, their shape and size. However the information is insufficient to be mapped into unique aerosol physical properties and vertical distribution. Infinite number of physical solutions can reconstruct the same two wavelength backscattered profile measured from space. CALIPSO will fly in formation with the Aqua satellite and the MODIS spectro-radiometer on board. Spectral radiances measured by MODIS in six channels between 0.55 and 2.13 microns simultaneously with the CALIPSO observations can constrain the solutions and resolve this ambiguity, albeit under some assumptions. In this paper we describe the inversion method and apply it to aircraft lidar and MODIS data collected over a dust storm off the coast of West Africa during the SHADE experiment. It is shown that the product of the single scattering albedo, omega, and the phase function, P, for backscattering can be retrieved from the synergism between measurements avoiding a priori hypotheses required for inverting lidar measurements alone. The resultant value of (omega)P(180 deg.) = 0.016/sr are significantly different from what is expected using Mie theory, but are in good agreement with recent results obtained from lidar observations of dust episodes. The inversion is robust in the presence of noise of 10% and 20% in the lidar signal in the 0.53 and 1.06 pm channels respectively. Calibration errors of the lidar of 5 to 10% can cause an error in optical thickness of 20 to 40% respectively in the tested cases. The lidar calibration errors cause degradation in the ability to fit the MODIS data. Therefore the MODIS measurements can be used to identify the calibration problem and correct for it. The CALIPSO-MODIS measurements of the profiles of fine and coarse aerosols, together with CALIPSO measurements of clouds vertical distribution, is expected to be critically important in understanding aerosol transport across continents and political boundaries, and to study aerosol-cloud interaction and its effect on precipitation and global forcing of climate.

Kaufman, Yoram↗

Misclassification in Workers’ Telecommuting Frequency Choices Using a Generalized Extreme Value Model

Telecommuting frequency is a response variable collected in travel surveys and is, therefore, prone to errors leading to mismeasurements or misclassification. Misclassification of explanatory variables is a common risk when using statistical modeling techniques. We define “misclassification” as a response reported or recorded in the wrong category; for example, a variable is recorded as a 1 when it should be 0. Here, in this context, this study aims to develop a statistical model to analyze telecommuting data which accounts for potential misclassification errors by building on existing literature in econometrics. The empirical analysis was undertaken using the 2017 National Household Travel Survey (NHTS) and the general extreme value (GEV) models available in the literature. Specifically, the frequency of telecommuting days was analyzed using the negative binomial (NB) model recast as the multinomial logit (MNL) model. By nature—and consistent with other studies—NHTS data are prone to errors that can be classified as intentional or unintentional misinformation provided by the person being interviewed. Ignoring these errors while modeling telecommuting frequencies using standard discrete count models can result in biased parameter estimates. The misclassification parameter was calculated for both over-reporting and under-reporting scenarios. The misclassification errors can be as high as 14% over-reported and 10% under-reported, particularly for the neighboring values. Statistical fit comparison between the models shows that models that ignore misclassification have worse data fit and biased parameter estimates with significant policy implications.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Searching for axionlike particles from core-collapse supernovae with Fermi LAT’s low-energy technique

Light axionlike particles (ALPs) are expected to be abundantly produced in core-collapse supernovae (CCSNe), resulting in a ~ 10 -second long burst of ALPs. These particles subsequently undergo conversion into gamma rays in external magnetic fields to produce a long gamma-ray burst (GRB) with a characteristic spectrum peaking in the 30–100-MeV energy range. At the same time, CCSNe are invoked as progenitors of ordinary long GRBs, rendering it relevant to conduct a comprehensive search for ALP spectral signatures using the observations of long GRBs with the Fermi Large Area Telescope (LAT). We perform a data-driven sensitivity analysis to determine CCSN distances for which a detection of an ALP signal is possible with the LAT’s low-energy technique which, in contrast to the standard LAT analysis, allows for a a larger effective area for energies down to 30 MeV. Assuming an ALP mass m a ≲ 10 - 10 eV and ALP-photon coupling g a γ = 5.3 × 10 - 12 GeV - 1 , values considered and deduced in ALP searches from SN1987A, we find that the distance limit ranges from ~ 0.5 to ~ 10 Mpc , depending on the sky location and the CCSN progenitor mass. Furthermore, we select a candidate sample of 24 GRBs and carry out a model comparison analysis in which we consider different GRB spectral models with and without an ALP signal component. We find that the inclusion of an ALP contribution does not result in any statistically significant improvement of the fits to the data. Finally, we discuss the statistical method used in our analysis and the underlying physical assumptions, the feasibility of setting upper limits on the ALP-photon coupling, and give an outlook on future telescopes in the context of ALP searches.

79 ASTRONOMY AND ASTROPHYSICS↗

Direct evaluation of the phase diagrams of dense multicomponent plasmas by integration of the Clapeyron equations

Accurate phase diagrams of multicomponent plasmas are required for the modeling of dense stellar plasmas, such as those found in the cores of white dwarf stars and the crusts of neutron stars. Those phase diagrams have been computed using a variety of standard techniques, which suffer from physical and computational limitations. Here we present an efficient and accurate method that overcomes the drawbacks of previously used approaches. In particular, finite-size effects are avoided as each phase is calculated separately; the plasma electrons and volume changes are explicitly taken into account; and arbitrary analytic fits to simulation data as well as particle insertions are avoided. Furthermore, no simulations at “uninteresting” state conditions, i.e., away from the phase coexistence curves, are required, which improves the efficiency of the technique. The method consists of an adaptation of the so-called Gibbs-Duhem integration approach to electron-ion plasmas, where the coexistence curve is determined by direct numerical integration of its underlying Clapeyron equation. The thermodynamics properties of the coexisting phases are evaluated separately using Monte Carlo simulations in the isobaric semigrand canonical ensemble (NPT Δ μ ). We describe this Monte Carlo-based Clapeyron integration method, including its basic physical and numerical principles, our extension to electron-ion plasmas, and our numerical implementation. We illustrate its applicability and benefits with the calculation of the melting curve of dense carbon-oxygen plasmas under conditions relevant for the cores of white dwarf stars and provide analytic fits to implement this new melting curve in white dwarf models. While this work focuses on the liquid-solid phase boundary of dense two-component plasmas, a wider range of physical systems and phase boundaries are within the scope of the Clapeyron integration method, which had until now only been applied to simple model systems of neutral particles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Analyzing Electric Field Morphology Through Data-Model Comparisons of the GEM IM/S Assessment Challenge Events

The storm-time inner magnetospheric electric field morphology and dynamics are assessed by comparing numerical modeling results of the plasmasphere and ring current with many in situ and remote sensing data sets. Two magnetic storms are analyzed, April 22,2001 and October 21-23,2001, which are the events selected for the Geospace Environment Modeling (GEM) Inner Magnetosphere/Storms (IM/S) Assessment Challenge (IMSAC). The IMSAC seeks to quantify the accuracy of inner magnetospheric models as well as synthesize our understanding of this region. For each storm, the ring current-atmosphere interaction model (RAM) and the dynamic global core plasma model (DGCPM) were run together with various settings for the large-scale convection electric field and the nightside ionospheric conductance. DGCPM plasmaspheric parameters were compared with IMAGE-EUV plasmapause extractions and LANL-MPA plume locations and velocities. RAM parameters were compared with Dst*, LANL-MPA fluxes and moments, IMAGE-MENA images, and IMAGE-HENA images. Both qualitative and quantitative comparisons were made to determine the electric field morphology that allows the model results to best fit the plasma data at various times during these events. The simulations with self-consistent electric fields were, in general, better than those with prescribed field choices. This indicates that the time-dependent modulation of the inner magnetospheric electric fields by the nightside ionosphere is quite significant for accurate determination of these fields (and their effects). It was determined that a shielded Volland-Stern field description driven by the 3-hour Kp index yields accurate results much of the time, but can be quite inconsistent. The modified Mcllwain field description clearly lagged in overall accuracy compared to the other fields, but matched some data sets (like Dst*) quite well. The rankings between the simulations varied depending on the storm and the individual data sets, indicating that each field description did well for some place, time, and energy range during the events, as well as doing less well in other places, times, and energies. Several unresolved issues regarding the storm-time inner magnetospheric electric field are discussed.

Liemohn, Michael W.↗