Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “vector 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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS1985

Evaluation of data obtained from atmospheric laser Doppler velocimeter

The Doppler lidar velocimeter produces a variety of kinds of data. Besides the line of sight velocity components, there are the reflected amplitudes and the Doppler widths. Computer graphics software for displaying these data was produced. Different methods of presentation are needed for the various items. A picture was produced as pictures are often the best way to gain understanding. The individual lidar scans form a crosshatch pattern. Thus one-dimensional measurements fit together to form a two-dimensional whole. A pair of velocity measurements at a point combine to form a wind measurement with direction and magnitude. This gives a forest of wind vectors. The goal is to find a recognizable pattern to these trees. Often it is necessary to show only part of the information. That is, show only the wind direction not its magnitude or reduce the wind to streamlines of air flow. In other cases data are reduced to contour plots. Just enough contour lines are included to show the picture described.

Keck, T. S.↗
NASA NTRS1982

Zonal harmonic model of Saturn's magnetic field from Voyager 1 and 2 observations

An analysis of the magnetic field of Saturn is presented which takes into account both the Voyager 1 and 2 vector magnetic field observations. The analysis is based on the traditional spherical harmonic expansion of a scale potential to derive the magnetic field within 8 Saturn radii. A third-order zonal harmonic model fitted to Voyager 1 and 2 observations is found to be capable of predicting the magnetic field characteristics at one encounter based on those observed at another, unlike models including dipole and quadrupole terms only. The third-order model is noted to lead to significantly enhanced polar surface field intensities with respect to dipole models, and probably represents the axisymmetric part of a complex dynamo field.

Connerney, J. E. P.↗
NASA NTRS1982

Sensitivity of constrained linear inversions to the selection of the Lagrange multiplier

The influence of the choice of the Lagrange multiplier on constrained linear inversions is explored, with reference made to applications in inferring the columnar aerosol size distributions from spectral aerosol optical depth measurements. A range of the Lagrange multiplier is examined to find all positive solutions for the solution vector, which represents modifying factors to the assumed form of the size distribution. An iterative method is devised to constrain the calculations to consideration of only positive quantities and a requirement that the regression fit to data be consistent with measurement errors. The determination of the variances and covariances is formulated and applied to existing data sets for optical depth. Variances in the solution are found to be large for particle radii when the information content of the data is small.

King, M. D.↗
NASA NTRS1988

Calculations of rate constants for the three-body recombination of H2 in the presence of H2

A new global potential energy hypersurface for H2 + H2 is constructed and quasiclassical trajectory calculations performed using the resonance complex theory and energy transfer mechanism to estimate the rate of three body recombination over the temperature range 100 to 5000 K. The new potential is a faithful representation of ab initio electron structure calculations, is unchanged under the operation of exchanging H atoms, and reproduces the accurate H3 potential as one H atom is pulled away. Included in the fitting procedure are geometries expected to be important when one H2 is near or above the dissociation limit. The dynamics calculations explicitly include the motion of all four atoms and are performed efficiently using a vectorized variable-stepsize integrator. The predicted rate constants are approximately a factor of two smaller than experimental estimates over a broad temperature range.

Schwenke, David W.↗
NASA NTRS1988

Calculations of rate constants for the three-body recombination of H2 in the presence of H2

A new global potential energy hypersurface for H2 + H2 is constructed and quasiclassical trajectory calculations performed using the resonance complex theory and energy transfer mechanism to estimate the rate of three body recombination over the temperature range 100 to 5000 K. The new potential is a faithful representation of ab initio electron structure calculations, is unchanged under the operation of exchanging H atoms, and reproduces the accurate H3 potential as one H atom is pulled away. Included in the fitting procedure are geometries expected to be important when one H2 is near or above the dissociation limit. The dynamics calculations explicitly include the motion of all four atoms and are performed efficiently using a vectorized variable-stepsize integrator. The predicted rate constants are approximately a factor of two smaller than experimental estimates over a broad temperature range.

Schwenke, David W.↗
NASA NTRS2006

A New Look at Some Solar Wind Turbulence Puzzles

Some aspects of solar wind turbulence have defied explanation. While it seems likely that the evolution of Alfvenicity and power spectra are largely explained by the shearing of an initial population of solar-generated Alfvenic fluctuations, the evolution of the anisotropies of the turbulence does not fit into the model so far. A two-component model, consisting of slab waves and quasi-two-dimensional fluctuations, offers some ideas, but does not account for the turning of both wave-vector-space power anisotropies and minimum variance directions in the fluctuating vectors as the Parker spiral turns. We will show observations that indicate that the minimum variance evolution is likely not due to traditional turbulence mechanisms, and offer arguments that the idea of two-component turbulence is at best a local approximation that is of little help in explaining the evolution of the fluctuations. Finally, time-permitting, we will discuss some observations that suggest that the low Alfvenicity of many regions of the solar wind in the inner heliosphere is not due to turbulent evolution, but rather to the existence of convected structures, including mini-clouds and other twisted flux tubes, that were formed with low Alfvenicity. There is still a role for turbulence in the above picture, but it is highly modified from the traditional views.

Roberts, Aaron↗
NASA NTRS2002

Ground Operations of the ISS GNC Babb-Mueller Atmospheric Density Model

The ISS GNC system was updated recently with a new software release that provides onboard state determination capability. Prior to this release, only the Russian segment maintained and propagated the onboard state, which was periodically updated through Russian ground tracking. The new software gives the US segment the capability for maintaining the onboard state, and includes new GPS and state vector propagation capabilities. Part of this software package is an atmospheric density model based on the Babb-Mueller algorithm. Babb-Mueller efficiently mimics a full analytical density model, such as the Jacchia model. While lacchia is very robust and is used in the Mission Control Center, it is too computationally intensive for use onboard. Thus, Babb-Mueller was chosen as an alternative. The onboard model depends on a set of calibration coefficients that produce a curve fit to the lacchia model. The ISS GNC system only maintains one set of coefficients onboard, so a new set must be uplinked by controllers when the atmospheric conditions change. The onboard density model provides a real-time density value, which is used to calculate the drag experienced by the ISS. This drag value is then incorporated into the onboard propagation of the state vector. The propagation of the state vector, and therefore operation of the BabbMueller algorithm, will be most critical when GPS updates and secondary state vector sources fail. When GPS is active, the onboard state vector will be updated every ten seconds, so the propagation error is irrelevant. When GPS is inactive, the state vector must be updated at least every 24 hours, based on current protocol. Therefore, the Babb-Mueller coefficients must be accurate enough to fulfill the state vector accuracy requirements for at least one day. A ground operations concept was needed in order to manage both the on board Babb-Mueller density model and the onboard state quality. The Babb-Mueller coefficients can be determined operationally in two ways. The first method is to calibrate the coefficients in real-time, where a set of custom coefficients is generated for the real-time atmospheric conditions. The second approach is to generate pre-canned sets of coefficients that encompass the expected atmospheric conditions over the lifetime of the vehicle. These predetermined sets are known as occurrences. Even though a particular occurrence will not match the true atmospheric conditions, the error will be constrained by limiting the breadth of each occurrence. Both methods were investigated and the advantages and disadvantages of each were considered. The choice between these implementations was a trade-off between the additional accuracy of the real-time calibration and the simpler development for the approach using occurrences. The operations concept for the frequency of updates was also explored, and depends on the deviation in solar flux that still achieves the necessary accuracy of the coefficients. This was determined based on historical solar flux trends. This analysis resulted in an accurate and reliable implementation of the Babb-Mueller coefficients and how flight controllers use them during realtime operations.

Brogan, Jonathan↗
NASA NTRS1978

The radial dependences of the interplanetary magnetic field between 1 and 5 AU Pioneer 10

Pioneer 10 vector helium magnetometer data acquired in 1972-1973 during Bartels solar rotations 1896-1918 are used to investigate the radial dependences of the distant interplanetary magnetic field (IMF) between 1 and 5 AU. Least-square fits were determined for the radial dependences of the averages of the magnitudes of IMF components and total field and plane projections, and radial fits were prepared for the standard deviations of these variables over the solar rotation, one day, and three-hour intervals. The variation of the weighted averages of the radial component of the field with respect to the heliocentric distance, the variation of the tangential component of the field, and the characteristics of a subset corresponding to a relatively low average solar wind velocity are reported.

Rosenberg, R. L.↗
NASA NTRS1988

Three-dimensional multigrid algorithms for the flux-split Euler equations

The Full Approximation Scheme (FAS) multigrid method is applied to several implicit flux-split algorithms for solving the three-dimensional Euler equations in a body fitted coordinate system. Each of the splitting algorithms uses a variation of approximate factorization and is implemented in a finite volume formulation. The algorithms are all vectorizable with little or no scalar computation required. The flux vectors are split into upwind components using both the splittings of Steger-Warming and Van Leer. The stability and smoothing rate of each of the schemes are examined using a Fourier analysis of the complete system of equations. Results are presented for three-dimensional subsonic, transonic, and supersonic flows which demonstrate substantially improved convergence rates with the multigrid algorithm. The influence of using both a V-cycle and a W-cycle on the convergence is examined.

Anderson, W. Kyle↗
NASA NTRS1995

Calculation of three-dimensional (3-D) internal flow by means of the velocity-vorticity formulation on a staggered grid

A method has been developed to accurately compute the viscous flow in three-dimensional (3-D) enclosures. This method is the 3-D extension of a two-dimensional (2-D) method developed for the calculation of flow over airfoils. The 2-D method has been tested extensively and has been shown to accurately reproduce experimental results. As in the 2-D method, the 3-D method provides for the non-iterative solution of the incompressible Navier-Stokes equations by means of a fully coupled implicit technique. The solution is calculated on a body fitted computational mesh incorporating a staggered grid methodology. In the staggered grid method, the three components of vorticity are defined at the centers of the computational cell sides, while the velocity components are defined as normal vectors at the centers of the computational cell faces. The staggered grid orientation provides for the accurate definition of the vorticity components at the vorticity locations, the divergence of vorticity at the mesh cell nodes and the conservation of mass at the mesh cell centers. The solution is obtained by utilizing a fractional step solution technique in the three coordinate directions. The boundary conditions for the vorticity and velocity are calculated implicitly as part of the solution. The method provides for the non-iterative solution of the flow field and satisfies the conservation of mass and divergence of vorticity to machine zero at each time step. To test the method, the calculation of simple driven cavity flows have been computed. The driven cavity flow is defined as the flow in an enclosure driven by a moving upper plate at the top of the enclosure. To demonstrate the ability of the method to predict the flow in arbitrary cavities, results will he shown for both cubic and curved cavities.

Stremel, Paul M.↗
NASA NTRS2013

Improving Photometry and Stellar Signal Preservation with Pixel-Level Systematic Error Correction

The Kepler Mission has demonstrated that excellent stellar photometric performance can be achieved using apertures constructed from optimally selected CCD pixels. The clever methods used to correct for systematic errors, while very successful, still have some limitations in their ability to extract long-term trends in stellar flux. They also leave poorly correlated bias sources, such as drifting moiré pattern, uncorrected. We will illustrate several approaches where applying systematic error correction algorithms to the pixel time series, rather than the co-added raw flux time series, provide significant advantages. Examples include, spatially localized determination of time varying moiré pattern biases, greater sensitivity to radiation-induced pixel sensitivity drops (SPSDs), improved precision of co-trending basis vectors (CBV), and a means of distinguishing the stellar variability from co-trending terms even when they are correlated. For the last item, the approach enables physical interpretation of appropriately scaled coefficients derived in the fit of pixel time series to the CBV as linear combinations of various spatial derivatives of the pixel response function (PRF). We demonstrate that the residuals of a fit of soderived pixel coefficients to various PRF-related components can be deterministically interpreted in terms of physically meaningful quantities, such as the component of the stellar flux time series which is correlated with the CBV, as well as, relative pixel gain, proper motion and parallax. The approach also enables us to parameterize and assess the limiting factors in the uncertainties in these quantities.

Kolodzijczak, Jeffrey J.↗
NASA NTRS1997

Steady-State Pursuit Is Driven by Object Motion Rather Than the Vector Average of Local Motions

We have previously shown that humans can pursue the motion of objects whose trajectories can be recovered only by spatio-temporal integration of local motion signals. We now explore the integration rule used to derive the target-motion signal driving pursuit. We measured the pursuit response of 4 observers (2 naive) to the motion of a line-figure diamond viewed through two vertical bar apertures (0.2 cd/square m). The comers were always occluded so that only four line segments (93 cd/square m) were visible behind the occluding foreground (38 cd/square m). The diamond was flattened (40 & 140 degree vertex angles) such that vector averaging of the local normal motions and vertical integration (e.g. IOC) yield very I or different predictions, analogous to using a Type II plaid. The diamond moved along Lissajous-figure trajectories (Ax = Ay = 2 degrees; TFx = 0.8 Hz; TFy = 0.4 Hz). We presented only 1.25 cycles and used 6 different randomly interleaved initial relative phases to minimize the role of predictive strategies. Observers were instructed to track the diamond and reported that its motion was always coherent (unlike type II plaids). Saccade-free portions of the horizontal and vertical eye-position traces sampled at 240 Hz were fit by separate sinusoids. Pursuit gain with respect to the diamond averaged 0.7 across subjects and directions. The ratio of the mean vertical to horizontal amplitude of the pursuit response was 1.7 +/- 0.7 averaged across subjects (1SD). This is close to the prediction of 1.0 from vertical motion-integration rules, but far from 7.7 predicted by vector averaging and infinity predicted by segment- or terminator-tracking strategies. Because there is no retinal motion which directly corresponds to the diamond's motion, steady-state pursuit of our "virtual" diamond is not closed-loop in the traditional sense. Thus, accurate pursuit is unlikely to result simply from local retinal negative feedback. We conclude that the signal driving steady-state pursuit is not the vector average of local motion signals, but rather a more vertical estimate of object motion, derived in extrastriate cortical areas beyond V1, perhaps NIT or MST.

Stone, Leland S.↗
NASA NTRS1976

Ocean color spectra measured off the Oregon coast - Characteristic vectors

The ocean color spectrum is defined as the ratio of the spectrum of light upwelled from the sea to the spectrum of light incident on the sea surface. Ocean color spectra, observed from an airplane flown over waters off Oregon, are analyzed. The original spectra are resolved into fifty-five wavelength bands, each 5 nm wide. The shapes of these spectra are parameterized by, and shown to be accurately recoverable from, their first four principal components. These components are the scalar projections of each spectrum on the first four characteristic vectors of the sample covariance matrix. Regression equations are found with which phytoplankton pigment concentration and water transparency may be estimated as linear functions of the principal components. Pigment concentration estimates thus obtained are imprecise. The poor fit is at least partly due to the inappropriateness of the linear regression model and the neglect of other optically important substances typically present in sea water.

Mueller, J. L.↗
NASA NTRS1996

Spatial orientation of semicircular canals and afferent sensitivity vectors in pigeons

Rotational head motion in vertebrates is detected by the semicircular canal system, whose innervating primary afferent fibers carry information about movement in specific head planes. The semicircular canals have been qualitatively examined over a number of years, and the canal planes have been quantitatively characterized in several animal species. The present study first determined the geometric relationship between individual semicircular canals and between the canals and the stereotactic head planes in pigeons. Stereotactic measurements of multiple points along the circumference of the bony canals were taken, and the measured points fitted with a three-dimensional planar surface. Direction normals to the plane's surface were calculated and used to define angles between semicircular canal pairs. Because of the unusual shape of the anterior semicircular canals in pigeons, two planes, a major and a minor, were fitted to the canal's course. Calculated angle values for all canals indicated that the horizontal and posterior semicircular canals are nearly orthogonal, but the anterior canals have substantial deviations from orthogonality with other canal planes. Next, the responses of the afferent fibers that innervate each of the semicircular canals to 0.5 Hz sinusoidal rotation about an earth-vertical axis were obtained. The head orientation relative to the rotation axis was systematically varied so that directions of maximum sensitivity for each canal afferent could be determined. These sensitivity vectors were then compared with the canal plane direction normals. The afferents that innervated specific semicircular canals formed homogeneous clusters of sensitivity vectors in different head planes. The horizontal and posterior afferents had average sensitivity vectors that were largely co-incident with the innervated canal plane direction normals. Anterior canal afferents, however, appeared to synthesize contributions from the major and minor plane components of the bony canal structure to produce a resultant sensitivity vector that was positioned between the canal planes. Calculated angles between the average canal afferent sensitivity vectors revealed that direction orthogonality is preserved at the afferent signal level, even though deviations from canal plane orthogonality exist.

NASA Discipline Neuroscience↗
NASA NTRS1979

CLASSY: An adaptive maximum likelihood clustering algorithm

The CLASSY clustering method alternates maximum likelihood iterative techniques for estimating the parameters of a mixture distribution with an adaptive procedure for splitting, combining, and eliminating the resultant components of the mixture. The adaptive procedure is based on maximizing the fit of a mixture of multivariate normal distributions to the observed data using its first through fourth central moments. It generates estimates of the number of multivariate normal components in the mixture as well as the proportion, mean vector, and covariance matrix for each component. The basic mathematical model for CLASSY and the actual operation of the algorithm as currently implemented are described. Results of applying CLASSY to real and simulated LANDSAT data are presented and compared with those generated by the iterative self-organizing clustering system algorithm on the same data sets.

Lennington, R. K.↗
NASA NTRS1982

Estimation of wave vector characteristics

A brief summary is presented of the techniques used to receive, transmit, and display frequency-time information. The mathematical basis for extracting wave vector information from the three electric and three magnetic wave fields is stated for the simple plane wave, single-source case. For the more realistic multiple-wave, multiple source case, basic correlation schemes and model-fitting techniques are described. Examples of results are given from various satellites for which two or more wave components could be treated. Finally, expectations for the upcoming OPEN mission are presented.

Shawhan, S. D.↗
NASA NTRS1995

Galileo flybys of Earth: The nature of the distant shock

We report on observations of the Earth's bow shock at unprecedentedly large downtail distances, some as remote as 360 R(sub E). Suprisingly, we find that even at these large distances, the bow shock signature remains clear. The cases we report are among the weakest shocks ever clearly identified. These shocks reveal patterns of field changes remarkably similar to those observed for stronger shocks. Indeed, several of the shocks could serve as textbook examples because they occur in unusually quiet and steady solar wind conditions. The quasi- perpendicular shocks are, in some cases preceded by whistler wave trains. Several of the shocks which have normal vectors in the transitional region between quasiperpendicular and quasiparallel are associated with large amplitude wave disturbances in the downstream plasma. Although low Mach number, these shocks are the source of large amplitude turbulence. We show a range of fits to the locations of the shock crossings, and argue that, for steady solar wind conditions, there is no sign of tail flapping. The displacement of the shock from the tail axis appears to be governed by the interplanetary magnetic field orientation as previously reported for Venus.

Kivelson, M. G.↗
NASA NTRS2019

Three CA-Related Covariance Issues and Their Solutions

This presentation addresses three issues that arise in the use of DoD-produced satellite state estimate covariances in the conjunction assessment process: the realism of the provided covariances, how to address correlated error between two satellites' covariance matrices, and how to proceed when a furnished covariance is non-positive-definite. To address the first, DoD has implemented a set of two "consider parameters" with which to alter/expand the covariance to account for atmospheric density forecast error and satellite frontal area uncertainty, the two largest sources of unmodeled position prediction error for LEO orbits; these values are governed by a satellite's orbital parameters and ballistic coefficient histories, as well as the current and predicted space weather situation. This approach substantially improves the realism of the covariance by accounting for expected prediction errors that are not part of the fit process that generates the covariance. The second issue of covariance correlation is a complicated one, but the main shared error source—global atmospheric density error that is common to both satellites can be characterized and, through sensitivity vectors, quantified at the conjunction's time of closest approach and removed from the joint covariance, which is used in the CA probability of collision (Pc) calculation. Finally, while according to the orbit determination theory a non-positive-definite (NPD) covariance is not possible, numerical truncation and covariance interpolation can conspire to produce NPD results, which render the covariance unusable for certain CA risk assessment calculations, such as Monte Carlo Pc determination. Three similar techniques are profiled and the simplest of them recommended as a reasonable remediation technique when NPD covariances are encountered by CA practitioners.

Dynamic Consider Parameter↗