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 181 records · Page 10

Multiple spacecraft observations of interplanetary shocks Four spacecraft determination of shock normals

ISEE 1, 2, 3, IMP 8, and Prognoz 7 observations of interplanetary shocks in 1978 and 1979 provide five instances where a single shock is observed by four spacecraft. These observations are used to determine best-fit normals for these five shocks. In addition to providing well-documented shocks for future investigations these data allow the evaluation of the accuracy of several shock normal determination techniques. When the angle between upstream and downstream magnetic field is greater than 20 deg, magnetic coplanarity can be an accurate single spacecraft method. However, no technique based solely on the magnetic measurements at one or multiple sites was universally accurate. Thus, the use of overdetermined shock normal solutions, utilizing plasma measurements, separation vectors, and time delays together with magnetic constraints, is recommended whenever possible.

Russell, C. T.↗

S-band Smallsat InSAR Constellation for Surface Deformation Science

Surface deformation studies using repeat-pass interferometric SAR have evolved into a powerful tool for geophysicists studying earthquake fault zones, volcanoes, ice sheet motion, and subterranean aquifers. Longer wavelengths (S-Band and L-Band) are preferred because they do not decorrelate as quickly as shorter wavelengths. Rapid revisit (1-3 days) is preferred because it allows the study of these phenomena at the timescales at which they commonly occur. Global access on such timescales is also required. Vector surface deformation measurements, taken from more than one direction, are a desired feature. This paper describes the conceptual architecture of a longer wave length, Smallsat SAR constellation of up to 12 satellites for rapid revisit surface deformation studies. The key to making such a constellation affordable is to lower launch costs, spacecraft costs, and instrument (SAR) costs. The first two objectives can be achieved using an ESPA-ring class, or Smallsat spacecraft. The third objective requires a SAR instrument sized to fit the mass and volume constraints imposed by such a spacecraft. Current state-of-the-art in miniaturization of electronics means that the radar transmit, receive and data handling functions can easily be implemented in a compact, low mass solution. The most significant challenge in designing a SAR to fit the Smallsat paradigm is in the dimensions of the antenna. The antenna sizing problem is addressed by adopting a smaller antenna than allowed by conventional SAR design rules. The baseline antenna design is simple, requiring no electronic beam-steering or beam-forming capability. Both reflectarray and microstrip patch antenna solutions are considered. The antenna structure is dual-purpose, to limit the overall system mass, with solar panels on the backplane providing power for the radar and spacecraft. The proposed solution easily accommodates radar squint angles of +/-30 degrees for repeat-pass interferometry measurements from multiple direct

Freeman, Anthony↗

Perturbations, TTVs & the (Un)reliability of Ephemerides of Kepler Planets

Fractional uncertainties quoted for orbital periods of the vast majority of planet candidates listed in recent Kepler catalogs are $<10^{-5}$, with values $\sim 10^{-6}$ (corresponding to 2 minutes per 4 years) being typical. These small uncertainties suggest ephemeris predictions for most Kepler planets are robust for decades to come. Tabulated periods list formal uncertainties in the mean times between midpoints of successive transits during the time interval in which transits were observed. However, transit timing variations (TTVs) produce errors in estimates of some planets’ orbital periods that need to be accounted for in certain dynamical investigations and ephemeris predictions. Periodic TTVs with timescales short compared to the interval of Kepler observations largely average out and do not produce significant errors in estimates of orbital periods. TTVs with timescales comparable to the interval of Kepler observations have been fit for dozens of Kepler planet candidates to estimate long-term average orbital periods (Holczer et al. 2016, ApJS 225, 9), and more detailed dynamical models have been used to estimate long-term average periods of a small number of well-studied planets. Most Kepler planets that show large TTVs are near mean-motion orbital resonances with other planets. The largest effect for planets moderately close to two-body resonances are due to rotation of the forced eccentricity vector by resonant perturbations; the timescale for this precession is typically short, so the variations tend to average out. Libration of planets locked in resonances typically occurs on timescales longer than the Kepler baseline, but few Kepler planets appear to be resonantly-locked. The more general but smaller (during the era of Kepler observations) effect is caused by secular precession of the planets' free eccentricities usually has a period much longer than the baseline of the Kepler observations, so it is not accounted for in estimates of mean periods or uncertainties. This precession causes a discrepancy between Kepler era mean orbital period and long-term mean orbital period can exist even for planets having TTVs that are too small to be observable during the epoch of Kepler observations.

Kepler↗

Thrust-induced effects on low-speed aerodynamics of fighter aircraft

Results of NASA Langley has conducted wind-tunnel investigations of several fighter configurations conducted to determine the effects of both thrust vectoring and spanwise blowing are reviewed. A recent joint NASA/Grumman Aerospace Corporation/U.S. Air Force Wright Aeronautical Laboratory wind-tunnel investigation was conducted to examine the effects of spanwise blowing on the trailing-edge flap system. This application contrasts with the more familiar method of spanwise blowing near the wing leading edge. Another joint program among NASA/McDonnell Aircraft Company/U.S. Air Force Wright Aeronautical Laboratory investigated the effects of reverse thrust on the low-speed aerodynamics of an F-15 configuration. The F-15 model was fitted with a rotating van thrust reverser concept which could simulate both in-flight reversing for approach and landing or full reversing for ground roll reduction. The significant results of these two joint programs are reported.

Banks, D. W.↗

A conservative approach for flow field calculations on multiple grids

In the computation of flow fields about complex configurations, it is very difficult to construct body-fitted coordinate systems. An alternative approach is to use several grids at once, each of which is generated independently. This procedure is called the 'multiple grids' or 'zonal grids' approach and its applications are investigated in this study. The method follows the conservative approach and provides conservation of fluxes at grid interfaces. The Euler equations are solved numerically on such grids for various configurations. The numerical scheme used is the finite-volume technique with a three-stage Runge-Kutta time integration. The code is vectorized and programmed to run on the CDC VPS-32 computer. Some steady state solutions of the Euler equations are presented and discussed.

Kathong, M.↗

A conservative approach for flow field calculations on multiple grids

In the computation of flow fields about complex configurations, it is very difficult to construct body-fitted coordinate systems. An alternative approach is to use several grids at once, each of which is generated independently. This procedure is called the multiple grids or zonal grids approach and its applications are investigated in this study. The method follows the conservative approach and provides conservation of fluxes at grid interfaces. The Euler equations are solved numerically on such grids for various configurations. The numerical scheme used is the finite-volume technique with a three-state Runge-Kutta time integration. The code is vectorized and programmed to run on the CDC VPS-32 computer. Some steady state solutions of the Euler equations are presented and discussed.

Kathong, Monchai↗

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↗

A Post-Pathfinder Evaluation of Areocentric Solar Coordinates with Improved Timing Recipes for Mars Seasonal/Diurnal Climate Studies

The accurate determination of the Mars pole vector derived from Pathfinder and Viking Lander radio data, together with the VSOP87 representation of planetary orbits, have been applied to a new evaluation of the right ascension of the "fictitious mean sun" (FMS) at Mars. With DELTA t (sub J2000) the elapsed time in days from the J2000 epoch (J.D.2451545.0 (sup TT), alpha FMS = 270 degrees.3863 + 0.52403840(degrees/d) (raised dot) DELTA T (sub j2000) - 4 x 10 (exp -13) (degrees/d (sup 2)) (raised dot) DELTA t (sup 2) (sub J2000) represents a best least-squares quadratic fit of the FMS, including aberration, to each instance of the four equinox and solstice passages for each of 134 Mars orbits spanning the calendar years 1874-2127. The implied tropical orbit period for Mars, 686.9726 (sup d), closely agrees with the recent evaluations. Together with the Pathfinder radio determination of the Mars sidereal rotation, the derived FMS rate corresponds to a mean solar day (or "sol") of 1.027491251 (sup d). The new FMS determination would serve to define the Mean Solar Time at Mars to the nearest tenth-second, according to historical conventions originally established for terrestrial time keeping, once the Mars prime meridian defined by the crater Airy-O is navigated to the same accuracy. For convenient reference to current epochs, 2000 Jan 06 00:00 UTC (= MJD 51549.000 (sup UTC)) corresponds to a coincidence of (alpha (sub FMS)) and the rotation angle of the crater Airy-O measured with respect to the Mars equinox (i.e. "mean solar midnight" on the planet's prime meridian), to within the current uncertainty of several seconds in the locational definition of the planet's cartographic grid. As a further result of the analysis, the consistently derived Mars obliquity of date is epsilon = 25 degrees.192 + 3.45 x l0 (exp -7)(degrees/d)(raised dot) DELTA t (sub J2000). An improved analytic recipe for the calculation of the solar areocentric longitude (L (sub s)) of Mars to an accuracy of 0 degrees.01 is also provided, accounting for the primary perturbations of Earth, Jupiter, and Venus, which may in turn be applied to an efficient evaluation of Mars local true solar time (LTST) to within the uncertainty of the inertial position of the Mars prime meridian. For specific applications to the data archives for landed Mars spacecraft, simple conversion formulae are given for the determination of the Viking "Local Lander Time" and the Pathfinder "Local True Solar Time" in terms of the terrestrial calendar date and UTC.

Allison, Michael↗

Aeroelastic Model Structure Computation for Envelope Expansion

Structure detection is a procedure for selecting a subset of candidate terms, from a full model description, that best describes the observed output. This is a necessary procedure to compute an efficient system description which may afford greater insight into the functionality of the system or a simpler controller design. Structure computation as a tool for black-box modeling may be of critical importance in the development of robust, parsimonious models for the flight-test community. Moreover, this approach may lead to efficient strategies for rapid envelope expansion that may save significant development time and costs. In this study, a least absolute shrinkage and selection operator (LASSO) technique is investigated for computing efficient model descriptions of non-linear aeroelastic systems. The LASSO minimises the residual sum of squares with the addition of an l(Sub 1) penalty term on the parameter vector of the traditional l(sub 2) minimisation problem. Its use for structure detection is a natural extension of this constrained minimisation approach to pseudo-linear regression problems which produces some model parameters that are exactly zero and, therefore, yields a parsimonious system description. Applicability of this technique for model structure computation for the F/A-18 (McDonnell Douglas, now The Boeing Company, Chicago, Illinois) Active Aeroelastic Wing project using flight test data is shown for several flight conditions (Mach numbers) by identifying a parsimonious system description with a high percent fit for cross-validated data.

Kukreja, Sunil L.↗

Drift and observations in cosmic-ray modulation, 1

It is illustrated that a relative simple drift model can, in contrast with no drift models, simultaneously fit proton and electron spectra observed in 1965-66 and 1977, using a single set of modulation parameters except for a change in the IMF polarity. This result is interpreted together with the observation of Evenson and Meyer that electrons are recovering more rapidly than protons after 1980, in contrast with what Burger and Swanenburg observed in 1968-72, as a charge sign dependent effect due to the occurrence of drift in cosmic ray modulation. The same set of parameters produces a shift in the phase and amplitude of the diurnal anisotropy vector, consistent with observations in 1969-71 and 1980-81.

Potgieter, M. S.↗

The probabilistic neural network architecture for high speed classification of remotely sensed imagery

In this paper we discuss a neural network architecture (the Probabilistic Neural Net or the PNN) that, to the best of our knowledge, has not previously been applied to remotely sensed data. The PNN is a supervised non-parametric classification algorithm as opposed to the Gaussian maximum likelihood classifier (GMLC). The PNN works by fitting a Gaussian kernel to each training point. The width of the Gaussian is controlled by a tuning parameter called the window width. If very small widths are used, the method is equivalent to the nearest neighbor method. For large windows, the PNN behaves like the GMLC. The basic implementation of the PNN requires no training time at all. In this respect it is far better than the commonly used backpropagation neural network which can be shown to take O(N6) time for training where N is the dimensionality of the input vector. In addition the PNN can be implemented in a feed forward mode in hardware. The disadvantage of the PNN is that it requires all the training data to be stored. Some solutions to this problem are discussed in the paper. Finally, we discuss the accuracy of the PNN with respect to the GMLC and the backpropagation neural network (BPNN). The PNN is shown to be better than GMLC and not as good as the BPNN with regards to classification accuracy.

Chettri, Samir R.↗

Conservation equations and physical models for hypersonic air flows in thermal and chemical nonequilibrium

The conservation equations for simulating hypersonic flows in thermal and chemical nonequilibrium and details of the associated physical models are presented. These details include the curve fits used for defining thermodynamic properties of the 11 species air model, curve fits for collision cross sections, expressions for transport properties, the chemical kinetics models, and the vibrational and electronic energy relaxation models. The expressions are formulated in the context of either a two or three temperature model. Greater emphasis is placed on the two temperature model in which it is assumed that the translational and rotational energy models are in equilibrium at the translational temperature, T, and the vibrational, electronic, and electron translational energy modes are in equilibrium at the vibrational temperature, T sub v. The eigenvalues and eigenvectors associated with the Jacobian of the flux vector are also presented in order to accommodate the upwind based numerical solutions of the complete equation set.

Peter A. Gnoffo↗

GPS-Based Reduced Dynamic Orbit Determination Using Accelerometer Data

Currently two gravity field satellite missions, CHAMP and GRACE, are equipped with high sensitivity electrostatic accelerometers, measuring the non-conservative forces acting on the spacecraft in three orthogonal directions. During the gravity field recovery these measurements help to separate gravitational and non-gravitational contributions in the observed orbit perturbations. For precise orbit determination purposes all these missions have a dual-frequency GPS receiver on board. The reduced dynamic technique combines the dense and accurate GPS observations with physical models of the forces acting on the spacecraft, complemented by empirical accelerations, which are stochastic parameters adjusted in the orbit determination process. When the spacecraft carries an accelerometer, these measured accelerations can be used to replace the models of the non-conservative forces, such as air drag and solar radiation pressure. This approach is implemented in a batch least-squares estimator of the GPS High Precision Orbit Determination Software Tools (GHOST), developed at DLR/GSOC and DEOS. It is extensively tested with data of the CHAMP and GRACE satellites. As accelerometer observations typically can be affected by an unknown scale factor and bias in each measurement direction, they require calibration during processing. Therefore the estimated state vector is augmented with six parameters: a scale and bias factor for the three axes. In order to converge efficiently to a good solution, reasonable a priori values for the bias factor are necessary. These are calculated by combining the mean value of the accelerometer observations with the mean value of the non-conservative force models and empirical accelerations, estimated when using these models. When replacing the non-conservative force models with accelerometer observations and still estimating empirical accelerations, a good orbit precision is achieved. 100 days of GRACE B data processing results in a mean orbit fit of a few centimeters with respect to high-quality JPL reference orbits. This shows a slightly better consistency compared to the case when using force models. A purely dynamic orbit, without estimating empirical accelerations thus only adjusting six state parameters and the bias and scale factors, gives an orbit fit for the GRACE B test case below the decimeter level. The in orbit calibrated accelerometer observations can be used to validate the modelled accelerations and estimated empirical accelerations computed with the GHOST tools. In along track direction they show the best resemblance, with a mean correlation coefficient of 93% for the same period. In radial and normal direction the correlation is smaller. During days of high solar activity the benefit of using accelerometer observations is clearly visible. The observations during these days show fluctuations which the modelled and empirical accelerations can not follow.

VanHelleputte, Tom↗

Remanent magnetization and 3-dimensional density model of the Kentucky anomaly region

A three-dimensional model of the Kentucky body was developed to fit surface gravity and long wavelength aeromagnetic data. Magnetization and density parameters for the model are much like those of Mayhew et al (1982). The magnetic anomaly due to the model at satellite altitude is shown to be much too small by itself to account for the anomaly measured by Magsat. It is demonstrated that the source region for the satellite anomaly is considerably more extensive than the Kentucky body sensu stricto. The extended source region is modeled first using prismatic model sources and then using dipole array sources. Magnetization directions for the source region found by inversion of various combinations of scalar and vector data are found to be close to the main field direction, implying the lack of a strong remanent component. It is shown by simulation that in a case (such as this) where the geometry of the source is known, if a strong remanent component is present its direction is readily detectable, but by scalar data as readily as vector data.

Mayhew, M. A.↗

Computations of Unsteady Viscous Compressible Flows Using Adaptive Mesh Refinement in Curvilinear Body-fitted Grid Systems

A methodology for accurate and efficient simulation of unsteady, compressible flows is presented. The cornerstones of the methodology are a special discretization of the Navier-Stokes equations on structured body-fitted grid systems and an efficient solution-adaptive mesh refinement technique for structured grids. The discretization employs an explicit multidimensional upwind scheme for the inviscid fluxes and an implicit treatment of the viscous terms. The mesh refinement technique is based on the AMR algorithm of Berger and Colella. In this approach, cells on each level of refinement are organized into a small number of topologically rectangular blocks, each containing several thousand cells. The small number of blocks leads to small overhead in managing data, while their size and regular topology means that a high degree of optimization can be achieved on computers with vector processors.

Steinthorsson, E.↗

Photoionization of the hydrogen atom in strong magnetic fields

The photoionization of the hydrogen atom in magnetic fields B about 10 exp 11 - 10 exp 13 G typical of the surface layers of neutron stars is investigated analytically and numerically. We consider the photoionization from various tightly bound and hydrogen-like states of the atom for photons with arbitrary polarizations and wave-vector directions. It is shown that the length form of the interaction matrix elements is more appropriate in the adiabatic approximation than the velocity form, at least in the most important frequency range omega much less than omega(B), where omega(B) is the electron cyclotron frequency. Use of the length form yields nonzero cross sections for photon polarizations perpendicular to the magnetic field at omega less than omega(B); these cross sections are the ones that most strongly affect the properties of the radiation escaping from an optically thick medium, e.g., from the atmosphere of a neutron star. The results of the numerical calculations are fitted by simple analytical formulas.

Potekhin, Aleksandr IU.↗

Time-domain analysis of planar microstrip devices using a generalized Yee-algorithm based on unstructured grids

The generalized Yee-algorithm is presented for the temporal full-wave analysis of planar microstrip devices. This algorithm has the significant advantage over the traditional Yee-algorithm in that it is based on unstructured and irregular grids. The robustness of the generalized Yee-algorithm is that structures that contain curved conductors or complex three-dimensional geometries can be more accurately, and much more conveniently modeled using standard automatic grid generation techniques. This generalized Yee-algorithm is based on the the time-marching solution of the discrete form of Maxwell's equations in their integral form. To this end, the electric and magnetic fields are discretized over a dual, irregular, and unstructured grid. The primary grid is assumed to be composed of general fitted polyhedra distributed throughout the volume. The secondary grid (or dual grid) is built up of the closed polyhedra whose edges connect the centroid's of adjacent primary cells, penetrating shared faces. Faraday's law and Ampere's law are used to update the fields normal to the primary and secondary grid faces, respectively. Subsequently, a correction scheme is introduced to project the normal fields onto the grid edges. It is shown that this scheme is stable, maintains second-order accuracy, and preserves the divergenceless nature of the flux densities. Finally, for computational efficiency the algorithm is structured as a series of sparse matrix-vector multiplications. Based on this scheme, the generalized Yee-algorithm has been implemented on vector and parallel high performance computers in a highly efficient manner.

Gedney, Stephen D.↗

Modal vector estimation for closely spaced frequency modes

Techniques for obtaining improved modal vector estimates for systems with closely spaced frequency modes are discussed. In describing the dynamical behavior of a complex structure modal parameters are often analyzed: undamped natural frequency, mode shape, modal mass, modal stiffness and modal damping. From both an analytical standpoint and an experimental standpoint, identification of modal parameters is more difficult if the system has repeated frequencies or even closely spaced frequencies. The more complex the structure, the more likely it is to have closely spaced frequencies. This makes it difficult to determine valid mode shapes using single shaker test methods. By employing band selectable analysis (zoom) techniques and by employing Kennedy-Pancu circle fitting or some multiple degree of freedom (MDOF) curve fit procedure, the usefulness of the single shaker approach can be extended.

Craig, R. R., Jr.↗