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 487 records · Page 27

Gamma-Ray Emission in Dissipative Pulsar Magnetospheres: from Theory to Fermi Observations

We compute the patterns of gamma-ray emission due to curvature radiation in dissipative pulsar magnetospheres. Our ultimate goal is to construct macrophysical models that are able to reproduce the observed gamma-ray light curve phenomenology recently published in the Second Fermi Pulsar Catalog. We apply specific forms of Ohm's law on the open field lines using a broad range for the macroscopic conductivity values that result in solutions ranging, from near-vacuum to near-force-free. Using these solutions, we generate model gamma-ray light curves by calculating realistic trajectories and Lorentz factors of radiating particles under the influence of both the accelerating electric fields and curvature radiation reaction. We further constrain our models using the observed dependence of the phase lags between the radio and gamma-ray emission on the gamma-ray peak separation. We perform a statistical comparison of our model radio-lag versus peak-separation diagram and the one obtained for the Fermi standard pulsars. We find that for models of uniform conductivity over the entire open magnetic field line region, agreement with observations favors higher values of this parameter. We find, however, significant improvement in fitting the data with models that employ a hybrid form of conductivity, specifically, infinite conductivity interior to the light cylinder and high but finite conductivity on the outside. In these models the gamma-ray emission is produced in regions near the equatorial current sheet but modulated by the local physical properties. These models have radio lags near the observed values and statistically best reproduce the observed light curve phenomenology. Additionally, they also produce GeV photon cut-off energies.

stars↗

Isolating Added Mass Load Components of CPAS Main Clusters

The current simulation for the Capsule Parachute Assembly System (CPAS) lacks fidelity in representing added mass for the 116 ft Do ringsail Main parachute. The availability of 3-D models of inflating Main canopies allowed for better estimation the enclosed air volume as a function of time. This was combined with trajectory state information to estimate the components making up measured axial loads. A proof-of-concept for an alternate simulation algorithm was developed based on enclosed volume as the primary independent variable rather than drag area growth. Databases of volume growth and parachute drag area vs. volume were developed for several flight tests. Other state information was read directly from test data, rather than numerically propagated. The resulting simulated peak loads were close in timing and magnitude to the measured loads data. However, results are very sensitive to data curve fitting and may not be suitable for Monte Carlo simulations. It was assumed that apparent mass was either negligible or a small fraction of enclosed mass, with little difference in results.

Ray, Eric S.↗

Optical Testing Using Portable Laser Coordinate Measuring Instruments

High precision, portable coordinate measuring instruments (CMI) such as laser radars (LR) and laser trackers (LT) have been used for optical system alignment and integration. The LRs ability to perform a non-contact scan of surfaces was previously utilized to characterize large spherical and aspheric mirrors. In this paper, we explore the use of a CMI as an accurate, fast, robust, and non-contact tool for prescription characterization of powered optical surfaces. Using Nikons MV-224350 LR and Leicas Absolute Tracker AT401402 instruments, proof of concept measurements were performed to characterize a variety of optical components by measuring the actual and apparent, or equivalently the direct and through (DT), coordinates of calibrated metrology targets. Custom macros in metrology software and other data reduction code were developed to compute surface-ray intercepts and surface slopes from the DT shots. The calculated data is fit to an aspheric surface formula to obtain the optimum prescription. The results were compared to the nominal parameters and were crosschecked using LR scans or other approaches. We discuss potential applications across the fields of optical component fabrication and system alignment and testing.

Khreishi, M.↗

Optical Testing Using Portable Laser Coordinate Measuring Instruments

High precision, portable coordinate measuring instruments (CMI) such as laser radars (LR) and laser trackers (LT) have been used for optical system alignment and integration. The LR's ability to perform a non-contact scan of surfaces was previously utilized to characterize large spherical and aspheric mirrors. In this paper, we explore the use of a CMI as an accurate, fast, robust, and non-contact tool for prescription characterization of powered optical surfaces. Using Nikon's MV-224/350 LR and Leica's Absolute Tracker AT401/402 instruments, proof of concept measurements were performed to characterize a variety of optical components by measuring the actual and apparent, or equivalently the "direct and through" (D&T), coordinates of calibrated metrology targets. Custom macros in metrology software and other data reduction code were developed to compute surface-ray intercepts and surface slopes from the D&T shots. The calculated data is fit to an aspheric surface formula to obtain the optimum prescription. The results were compared to the nominal parameters and were crosschecked using LR scans or other approaches. We discuss potential applications across the fields of optical component fabrication and system alignment and testing.

Khreishi, Manal↗

The JPL Neptune Radiation Model (NMOD)

The objective of this study is the development of a comprehensive radiation model of the Neptunian environment for JPL mission planning. The ultimate goal is to provide a description of the high-energy electron and proton environments and the magnetic field at Neptune that can be used for engineering design. The JPL Neptune Radiation Model (NMOD) models the high-energy electrons and protons between 0.025 MeV and 5 MeV based on the California Institute of Technology's Cosmic Ray Subsystem and the Applied Physics Laboratory's Low Energy Charged Particle Detector on Voyager 2. As in previous JPL radiation models, the form of the Neptunian model is based on magnetic field coordinates and requires a conversion from spacecraft coordinates to Neptunian-centered magnetic "B-L" coordinates. Two types of magnetic field models have been developed for Neptune: 1) simple "offset, tilted dipoles" (OTD), and 2) a complex, multi-pole expansion model ("O8"). A review of the existing data on Neptune and a search of the NASA Planetary Data System (PDS) were completed to obtain the most current descriptions of the Neptunian high-energy particle environment. These data were fit in terms of the O8 B-L coordinates to develop the electron and proton flux models. The flux predictions of the new model were used to estimate the total ionizing dose (TID) rate along the Neptunian equator, meridional flux contours for the electrons and protons, and for flux and dose comparisons with the other radiation belts in the Solar System.

Garrett, Henry↗

The JPL Neptune Radiation Model (NMOD)

The objective of this study is the development of a comprehensive radiation model of the Neptunian environment for JPL mission planning. The ultimate goal is to provide a description of the high-energy electron and proton environments and the magnetic field at Neptune that can be used for engineering design. The JPL Neptune Radiation Model (NMOD) models the high-energy electrons and protons between 0.025 MeV and 5 MeV based on the California Institute of Technology's Cosmic Ray Subsystem and the Applied Physics Laboratory's Low Energy Charged Particle Detector on Voyager 2. As in previous JPL radiation models, the form of the Neptunian model is based on magnetic field coordinates and requires a conversion from spacecraft coordinates to Neptunian-centered magnetic "B-L" coordinates. Two types of magnetic field models have been developed for Neptune: 1) simple "offset, tilted dipoles" (OTD), and 2) a complex, multi-pole expansion model ("O8"). A review of the existing data on Neptune and a search of the NASA Planetary Data System (PDS) were completed to obtain the most current descriptions of the Neptunian high-energy particle environment. These data were fit in terms of the O8 B-L coordinates to develop the electron and proton flux models. The flux predictions of the new model were used to estimate the total ionizing dose (TID) rate along the Neptunian equator, meridional flux contours for the electrons and protons, and for flux and dose comparisons with the other radiation belts in the Solar System.

magnetic field↗

Limits On Undetected Planets in the Six Transiting Planets Kepler-11 System

The Kepler-11 has five inner planets ranging from approx. 2 - 1 times as massive Earth in a tightly-packed configuration, with orbital periods between 10 and 47 days. A sixth planet, Kepler-11 g, with a period of118 days, is also observed. The spacing between planets Kepler-11 f and Kepler-11 g is wide enough to allow room for a planet to orbit stably between them. We compare six and seven planet fits to measured transit timing variations (TTVs) of the six known planets. We find that in most cases an additional planet between Kepler-11 f and Kepler-11 g degrades rather than enhances the fit to the TTV data, and where the fit is improved, the improvement provides no significant evidence of a planet between Kepler-11 f and Kepler-11 g. This implies that any planet in this region must be low in mass. We also provide constraints on undiscovered planets orbiting exterior to Kepler-11 g. representations will be described.

Lissauer, Jack↗

A Fundamental Plane for Gamma-Ray Pulsars

We show that the γ-ray pulsar observables, i.e., their total γ-ray luminosity, L(sub γ), spectral cutoff energy, ϵ (sub cut), stellar surface magnetic field, B(sub ⋆), and spin-down power ε, obey a relation of the form L(sub γ) = f (ϵ (sub cut), B(sub ⋆), ε ), which represents a 3D plane in their 4D logspace. Fitting the data of 88 pulsars of the second Fermi pulsar catalog, we show this relation to be L(sub γ) ∝ ϵ (sub cut)(sup 1.18±0.24 B(sub ⋆)(sup 0.17± 0.05) ε(sup 0.41±0.08), a pulsar fundamental plane (FP). We show that the observed FP is remarkably close to the theoretical relation L(sub γ) ∝ ϵ (sub cut)(sup 4/3) B(sub ⋆)(sup 1/6) ε(sup 5/12) obtained assuming that the pulsar γ-ray emission is due to curvature radiation by particles accelerated at the pulsar equatorial current sheet just outside the light cylinder. Interestingly, the FP seems incompatible with emission by synchrotron radiation. The corresponding scatter about the FP is ∼0.35 dex and can only partly be explained by the observational errors while the rest is probably due to the variation of the inclination and observer angles. We predict also that ϵ (sub cut) ∝ ε (sup 7/16) toward low ε for both young and millisecond pulsars implying that the observed death line of γ-ray pulsars is due to ϵ (sub cut) dropping below the Fermi band. Our results provide a comprehensive interpretation of the observations of γ-ray pulsars, setting requirements for successful theoretical modeling.

Kalapotharakos, Constantinos↗

Numerical Investigation and Optimization of a Flushwall Injector for Scramjet Applications at Hypervelocity Flow Conditions

An investigation utilizing Reynolds-averaged simulations (RAS) was performed in order to find optimal designs for an interdigitated flushwall injector suitable for scramjet applications at hypervelocity conditions. The flight Mach number, duct height, spanwise width, and injection angle were the design variables selected to maximize two objective functions: the thrust potential and combustion efficiency. A Latin hypercube sampling design-of-experiments method was used to select design points for RAS. A methodology was developed that automated building geometries and generating grids for each design. The ensuing RAS analysis generated the performance database from which the two objective functions of interest were computed using a one-dimensional performance utility. The data were fitted using four surrogate models: an artificial neural network (ANN) model, a cubic polynomial, a quadratic polynomial, and a Kriging model. Variance-based decomposition showed that both objective functions were primarily driven by changes in the duct height. Multiobjective design optimization was performed for all four surrogate models via a genetic algorithm method. Optimal solutions were obtained at the upper and lower bounds of the flight Mach number range. The Kriging model obtained an optimal solution set that predicted high values for both objective functions. Additionally, three challenge points were selected to assess the designs on the Pareto fronts. Further sampling among the designs of the Pareto fronts are required in order to lower the errors and perform more accurate surrogate-based optimization. sed optimization.

Shenoy, Rajiv R.↗

Numerical Investigation and Optimization of a Flushwall Injector for Scramjet Applications at Hypervelocity Flow Conditions

An investigation utilizing Reynolds-averaged simulations (RAS) was performed in order to demonstrate the use of design and analysis of computer experiments (DACE) methods in Sandia’s DAKOTA software package for surrogate modeling and optimization. These methods were applied to a flow- path fueled with an interdigitated flushwall injector suitable for scramjet applications at hyper- velocity conditions and ascending along a constant dynamic pressure flight trajectory. The flight Mach number, duct height, spanwise width, and injection angle were the design variables selected to maximize two objective functions: the thrust potential and combustion efficiency. Because the RAS of this case are computationally expensive, surrogate models are used for optimization. To build a surrogate model a RAS database is created. The sequence of the design variables comprising the database were generated using a Latin hypercube sampling (LHS) method. A methodology was also developed to automatically build geometries and generate structured grids for each design point. The ensuing RAS analysis generated the simulation database from which the two objective functions were computed using a one-dimensionalization (1D) of the three-dimensional simulation data. The data were fitted using four surrogate models: an artificial neural network (ANN), a cubic polynomial, a quadratic polynomial, and a Kriging model. Variance-based decomposition showed that both objective functions were primarily driven by changes in the duct height. Multiobjective design optimization was performed for all four surrogate models via a genetic algorithm method. Optimal solutions were obtained at the upper and lower bounds of the flight Mach number range. The Kriging model predicted an optimal solution set that exhibited high values for both objective functions. Additionally, three challenge points were selected to assess the designs on the Pareto fronts. Further sampling among the designs of the Pareto fronts may be required to lower the surrogate model errors and perform more accurate surrogate-model-based optimization.

Shenoy, Rajiv R.↗

Automatic Design of Missions to Small Bodies

The new JPL Small-Body Automatic Mission-Design System comprises two main elements: a database of pre-computed mission options to all known asteroids and comets, and an interactive web interface that can be used to design transfers to each small body. The system is kept current with the JPL Small Body orbit catalog, is publicly available, and can be accessed from the JPL Solar System Dynamics Group website. The missions computed by the automatic system are impulsive. However, a low-thrust v estimate is also provided. The database of pre-computed missions can be filtered to find potential targets with certain orbital and physical properties, and that meet specific mission-design constraints. In addition to describing the system in detail, this paper presents empirical analytic expressions to approximate the impulsive v requirements of missions to comets and to each family of asteroids, obtained by fitting statistical data. To show how the interactive interface works, we consider mission options to asteroid 99942 Apophis. We also include optimal mission opportunities to other selected small bodies.

Farnocchia, Davide↗

Juno Orbit Determination Experience During First Year At Jupiter

The Juno spacecraft successfully inserted into a polar orbit around Jupiter on 5- July-2016. Since the Jupiter Orbit Insertion (JOI) maneuver, Juno has completed six orbits around Jupiter. The mission plan at the time of JOI was for Juno to perform two 53.5-day capture orbits before executing a Period Reduction Maneuver (PRM) to place the spacecraft into its intended 14-day science orbit. This maneuver was canceled due to a concern with the propulsion system. As a result, the Juno spacecraft will remain in its longer orbit period for rest of its mission. This paper discusses the Navigation Team’s experience: the orbit determination strategy and how it changed due to the cancellation of the PRM, challenges fitting the data during perijove, and how we reconstructed the trajectory during Juno’s first year in orbit.

Thompson, Paul↗

Recommended Direct Simulation Monte Carlo Collision Model Parameters for Reacting Methane Flows

While hydrocarbon combustion and gas-surface interactions have been simulated for a variety of applications in continuum solvers, such processes have been rarely investigated in Direct Simulation Monte Carlo (DSMC). Thus Variable Hard/Soft Sphere (VHS/VSS)collision parameters are not available for most molecules of interest in the decomposition of methane, a significant pyrolysis gas product and potential atmospheric species for Titan re-entries. As these properties are essential for accurate modeling of gas transport and thermochemical effects, a method must be devised to compute them for use in DSMC. In this study, the methodology outlined in Stephani et al is utilized to compute the VHS/VSS parameters by fitting the Ω(1,1)and Ω(2,2) collision integrals [1]. These collision integrals are typically computed from intermolecular potentials and are used in CFD simulations to calculate transport properties such as viscosity, thermal conductivity, and diffusion. This approach ensures consistency between micro-scale collisions and macro-scale transport properties. To produce the necessary collision integrals needed for the fitting routine, we utilize the method described by Laricchiuta et al, where collision integrals are computed based on the constituent particles’ polarizability and effective electron number [2].Preliminary computations of pure-species Ω(1,1)and Ω(2,2)collision integrals produced from the coupled approach of Stephani et al and Laricchiuta et al for methane show good agreement with those computed via available ab-initio Potential Energy Surface Computations, as seen in Figure 1. Above 300K, no more than 12% deviation is observed between the two results over the temperature range of the fitted PES data [3].Current results of this study include full collision-specificVSS transport properties for a complex-species mixture involving major byproducts of pure methane decomposition. Targeted final results will include a comprehensive list of VSS parameters relevant toDSMC simulations of Titan atmospheric entry, based on the reduced kinetic model of Savajano et al [4].

M R Gosma↗

Model for Dielectric Constant of Seawater based on L-band Measurements with Conductivity by Definition

This article reports an improvement in the model for the dielectric constant of seawater used to fit laboratory measurements of the dielectric constant at the L-band. The new model (dielectric constant as a function of salinity, temperature, and frequency) is based on the response of a polar molecule proposed by Debye and fits the same measurement as reported in earlier work but uses a functional form for conductivity, σ(S,T) , that is given by the definition of salinity. The new version of this model fits the data well and has the advantages that the relaxation time constant is allowed to be a function of temperature and salinity and is well behaved when extrapolated to high salinities.

L-band↗

Quantifying urban, industrial, and background changes in NO2 during the COVID-19 lockdown period based on TROPOMI satellite observations

The COVID-19 lockdown had a large impact on anthropogenic emissions of air pollutants and particularly on nitrogen dioxide (NO 2 ). While the overall NO 2 decline over some large cities is well-established, understanding the details remains a challenge since multiple source categories contribute. In this study, a new method of isolation of three components (background NO 2 , NO 2 from urban sources, and NO 2 from industrial point sources) is applied to estimate the impact of the COVID-19 lockdown on each of them. The approach is based on fitting satellite data by a statistical model with empirical plume dispersion functions driven by a meteorological reanalysis. Population density and surface elevation data as well as coordinates of industrial sources were used in the analysis. The tropospheric NO 2 vertical column density (VCD) values measured by the Tropospheric Monitoring Instrument (TROPOMI) on board the Sentinel-5 Precursor over 261 urban areas for the period from 16 March to 15 June 2020 were compared with the average VCD values for the same period in 2018 and 2019. While the background NO 2 component remained almost unchanged, the urban NO 2 component declined by −18 % to −28 % over most regions. India, South America, and a part of Europe (particularly, Italy, France, and Spain) demonstrated a −40 % to −50 % urban emission decline. In contrast, the decline over urban areas in China, where the lockdown was over during the analysed period, was, on average, only -4.4± 8 %. Emissions from large industrial sources in the analysed urban areas varied greatly from region to region from -4.8 ± 6 % for China to -40 ± 10 % for India. Estimated changes in urban emissions are correlated with changes in Google mobility data (the correlation coefficient is 0.62) confirming that changes in traffic were one of the key elements in the decline in urban NO2 emissions. No correlation was found between changes in background NO2 and Google mobility data. On the global scale, the background and urban components were remarkably stable in 2018, 2019, and 2021, with averages of all analysed areas all being within ±2.5 % and suggesting that there were no substantial drifts or shifts in TROPOMI data. The 2020 data are clearly an outlier: in 2020, the mean background component for all analysed areas (without China) was -6.0 % ± 1.2 % and the mean urban component was -26.7 ± 2.6 % or 20σ below the baseline level from the other years.

COVID-19↗

Accounting for Dose Uncertainty in Dose-Response Curve Estimation Using Hierarchical Bayes Models

As part of their development of the technology for low-boom supersonic flight, the National Aeronautics and Space Administration (NASA) is planning to conduct a set of supersonic aircraft annoyance surveys in select communities in the United States, to measure public perception of the reduced sonic boom. The relationship between a noise exposure level of a supersonic flight event and the probability of an individual being highly annoyed by the event is quantified by a dose-response curve, which is usually based on logistic regression modeling. Unavoidable amounts of estimation error are expected in the upcoming noise measurements, which can bias the results of the logistic regression if ignored. Hence, the development of an estimation approach that accounts for such error when estimating the dose-response curve is imperative. In this paper, an evaluation study is conducted to assess the impact of measurement error on dose-response curve estimation. For this, hierarchical Bayes models using different specifications are fit to data collected by NASA in 2018, as part of a risk-reduction study of supersonic flights affecting Galveston, Texas. It is observed that the estimated dose-response curve appears sensitive to uncertainty in dose measurements, being subject to attenuation bias.

community response↗

Uncertainty Estimates for Fitting Zernike Polynomials to Discrete Data

Zernike polynomials are a widely used metric in modern optical analysis. They conveniently represent surfaces as a series of weighted terms corresponding to various optical aberrations. Ideally, each term is independent of others in the series, but Zernike polynomials lose this property when working with sets of discrete data. This gives rise to uncertainty in each polynomial’s actual contribution and affects metrology and simulation estimates of their relative weights. Several factors influencing these estimates are the number and arrangement of sample locations, the method for calculating the weights, and the total number of Zernike terms used in the calculation. Discussed is the uncertainty associated with linear regression using random sampling. Other topics reviewed are complex Zernike polynomials and vector spaces of functions.

Zernike Polynomials↗

Calibration of the COBE FIRAS instrument

The Far-Infrared Absolute Spectrophotometer (FIRAS) instrument on the Cosmic Background Explorer (COBE) satellite was designed to accurately measure the spectrum of the cosmic microwave background radiation (CMBR) in the frequency range 1-95/cm with an angular resolution of 7 deg. We describe the calibration of this instrument, including the method of obtaining calibration data, reduction of data, the instrument model, fitting the model to the calibration data, and application of the resulting model solution to sky observations. The instrument model fits well for calibration data that resemble sky condition. The method of propagating detector noise through the calibration process to yield a covariance matrix of the calibrated sky data is described. The final uncertainties are variable both in frequency and position, but for a typical calibrated sky 2.6 deg square pixel and 0.7/cm spectral element the random detector noise limit is of order of a few times 10(exp -7) ergs/sq cm/s/sr cm for 2-20/cm, and the difference between the sky and the best-fit cosmic blackbody can be measured with a gain uncertainty of less than 3%.

Fixsen, D. J.↗