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 37 records · Page 2

The effect of line damping, magneto-optics and parasitic light on the derivation of sunspot vector magnetic fields

The least square fitting of Stokes observations of sunspots using a Milne-Eddington-Unno model appears to lead, in many circumstances, to various inconsistencies such as anomalously large doppler widths and, hence, small magnetic fields which are significantly below those inferred solely from the Zeeman splitting in the intensity profile. It is found that the introduction of additional physics into the model such as the inclusion of damping wings and magneto-optic birefrigence significantly improves the fit to Stokes parameters. Model fits excluding the intensity profile, i.e., of both magnitude as well as spectral shape of the polarization parameters alone, suggest that parasitic light in the intensity profile may also be a source of inconsistencies. The consequences of the physical changes on the vector properties of the field derived from the Fe I lambda 6173 line for the 17 November 1975 spot as well as on the thermodynamic state are discussed. A Doppler width delta lambda (D) - 25mA is bound to be consistent with a low spot temperature and microturbulence, and a damping constant of a = 0.2.

Skumanich, A.↗

The modeling and prediction of multiple jet VTOL aircraft flow fields in ground effect

An engineering methodology based on an empirical data base and analytical fluid dynamic models was developed for the prediction of propulsive lift system induced aerodynamic effects for multiple lift jet VTOL aircraft operating in the hover mode in and out of ground effect. The effects of aircraft geometry, aircraft orientation (pitch, roll) as well as height above ground are considered. Lift jet vector and splay directions fit the airframe, lift jet exit flow conditions, and both axisymmetric and rectangular nozzle exit geometry are also accommodated. The induced suckdown flows are computed from the potential flowfield induced by the turbulent entrainment of both the free jets and wall jets in ground effect and from the free jets alone out of ground effect. The methodology emphasized geometric considerations, computation of stagnation lines and fountain upwash inclination, fountain upwash formation and development, and fountain impingement on the airframe.

Kotansky, D. R.↗

Oblique convergence and deformation along the Kuril and Japan trenches

The hypothesis that present-day deformation within the southern Kuril forearc is driven by oblique subduction of the Pacific plate is tested using 397 horizontal slip directions derived from shallow-thrust earthquakes from the Kuril and Japan trenches for the period 1963-1991. A simple two-plate model fits the 397 slip vectors significantly worse than a model that permits strike-slip motion of the southern Kuril forearc relative to the overlying plate. Weighted, mean slip directions along the southern Kuril trench are systematically rotated toward the direction orthogonal to the trench, which implies that the net convergence is partitioned into less oblique subduction and trench-parallel displacement of the southern Kuril forearc. The angular discrepancy between the observed slip direction and the direction predicted by the NUVEL-1 Pacific-North America Euler vector implies that the southern Kuril forearc translates 6-11 mm/yr to the southwest relative to the overlying North American plate. These results are consistent with geologically, geodetically, and seismologically observed convergence at the leading edge of the forearc sliver in southern Kokkaido and with previously inferred extension at the trailing edge of the sliver, which is located at the Bussol Strait at 46 deg N.

Demets, Charles↗

Autonomous frequency domain identification: Theory and experiment

The analysis, design, and on-orbit tuning of robust controllers require more information about the plant than simply a nominal estimate of the plant transfer function. Information is also required concerning the uncertainty in the nominal estimate, or more generally, the identification of a model set within which the true plant is known to lie. The identification methodology that was developed and experimentally demonstrated makes use of a simple but useful characterization of the model uncertainty based on the output error. This is a characterization of the additive uncertainty in the plant model, which has found considerable use in many robust control analysis and synthesis techniques. The identification process is initiated by a stochastic input u which is applied to the plant p giving rise to the output. Spectral estimation (h = P sub uy/P sub uu) is used as an estimate of p and the model order is estimated using the produce moment matrix (PMM) method. A parametric model unit direction vector p is then determined by curve fitting the spectral estimate to a rational transfer function. The additive uncertainty delta sub m = p - unit direction vector p is then estimated by the cross spectral estimate delta = P sub ue/P sub uu where e = y - unit direction vectory y is the output error, and unit direction vector y = unit direction vector pu is the computed output of the parametric model subjected to the actual input u. The experimental results demonstrate the curve fitting algorithm produces the reduced-order plant model which minimizes the additive uncertainty. The nominal transfer function estimate unit direction vector p and the estimate delta of the additive uncertainty delta sub m are subsequently available to be used for optimization of robust controller performance and stability.

Yam, Yeung↗

Ion Diffusive Transport Across the Separatrix Between the Low-Latitude Mantle and the Plasma Sheet By Kinetic Alfvén Waves: MMS Observation

To understand the entry of the cool low‐latitude mantle ions into the tail plasma sheet near the flanks under persistent interplanetary magnetic field B y , we evaluate the role of the cross‐field diffusive transport by kinetic Alfvén waves (KAWs) by investigating two events observed by multiscale (MMS) spacecraft. Around the separatrix between the open and closed field‐line regions, a two‐component mixing of hot plasma sheet ions of a few keV with cool mantle ions of a few hundred eV was observed, indicating transport across the separatrix. The waves observed between 0.01 and 10 Hz around the separatrix had characteristics consistent with those of KAWs. The consistency allowed us to estimate the wave vectors as a function of frequency by fitting KAW dispersion to the observations. Using the observed wave powers, plasma moments, and the estimated wave vectors, we computed the cross‐field diffusion rates associated with KAWs. The diffusion rates were found to be comparable to or larger than the Bohm diffusion rates during the intervals when the two‐component mixing was observed, indicating that the KAW diffusive transport can play a role in the entry of low‐latitude mantle ions into the plasma sheet.

Kinetic Alfvén waves↗

Polar decomposition for attitude determination from vector observations

This work treats the problem of weighted least squares fitting of a 3D Euclidean-coordinate transformation matrix to a set of unit vectors measured in the reference and transformed coordinates. A closed-form analytic solution to the problem is re-derived. The fact that the solution is the closest orthogonal matrix to some matrix defined on the measured vectors and their weights is clearly demonstrated. Several known algorithms for computing the analytic closed form solution are considered. An algorithm is discussed which is based on the polar decomposition of matrices into the closest unitary matrix to the decomposed matrix and a Hermitian matrix. A somewhat longer improved algorithm is suggested too. A comparison of several algorithms is carried out using simulated data as well as real data from the Upper Atmosphere Research Satellite. The comparison is based on accuracy and time consumption. It is concluded that the algorithms based on polar decomposition yield a simple although somewhat less accurate solution. The precision of the latter algorithms increase with the number of the measured vectors and with the accuracy of their measurement.

Bar-Itzhack, Itzhack Y.↗

MAGSAT and aeromagnetic data in the North American continent

Problems were encountered in deriving a proper reference field to be subtracted from the aeromagnetic data obtained from Project MAGNET. Field models tried thus far do not seem to eliminate properly the main field. The MAGSAT data in the North American continent for the period November 1 to December 22, 1979 are being compiled and compared with MAGNET data. Efforts are being made to eliminate the orbital bias errors. A computer program was developed and successfully tested which computes a topographic profile of the Curie depth isotherm which fits best to the observed vector or scalar field magnetic data.

Source record↗

Luneburg lens and optical matrix algebra research

Planar, as opposed to channelized, integrated optical circuits (IOCs) were stressed as the basis for computational devices. Both fully-parallel and systolic architectures are considered and the tradeoffs between the two device types are discussed. The Kalman filter approach is a most important computational method for many NASA problems. This approach to deriving a best-fit estimate for the state vector describing a large system leads to matrix sizes which are beyond the predicted capacities of planar IOCs. This problem is overcome by matrix partitioning, and several architectures for accomplishing this are described. The Luneburg lens work has involved development of lens design techniques, design of mask arrangements for producing lenses of desired shape, investigation of optical and chemical properties of arsenic trisulfide films, deposition of lenses both by thermal evaporation and by RF sputtering, optical testing of these lenses, modification of lens properties through ultraviolet irradiation, and comparison of measured lens properties with those expected from ray trace analyses.

Wood, V. E.↗

Evaluation of X-Ray Reflectors by Optical Diffraction Patterns

Performance of X-ray reflectors affects that of X-ray mirrors. Modern X-ray mirrors have thousands of reflectors to gain large effective area. Evaluation of the reflectors is an important process in production of the mirrors. A diffraction pattern dominates reflector image when the parallel optical beam illuminates the reflector along its optical axis because the reflectors are used at grazing incident angles of around 1 deg and their effective width are 1–10 mm. A diffraction pattern from the entire reflector surface can be acquired at once with the aid of a lens. The diffraction pattern holds information of the surface profiles of the reflectors. To quantitatively evaluate the reflectors with the diffraction pattern, we created a diffraction pattern model by the wave optics with the ideal surface profile and fitted it to data. As a result, a correlation between fitting residual and the normal vector distribution of the surface profile was found. With our method, the reflectors can be evaluated and sorted out more efficiently.

Hayashi, Takayuki↗

DFW microburst model based on AA-539 data

Analysis of the August 2, 1985 crash for an L-1011 jumbo jet (DL-191) on approach to the Dallas-Ft. Worth International Airport (DFW) in a thunderstorm indicates that the severe windshear microburst that caused the crash was composed not only of a strong downflow and outflow but also included several large-scale vortex rings entrained in the flowfield. This paper presents a detailed two-dimensional model of the DFW microburst based on data from the MD-80 (AA-539) that followed behind DL-191 and flew through the microburst about two minutes after the crash of DL-191. The model was developed using wind-vector and flight-path data reconstructed by NASA Ames Research Center and a combination of interactive graphics and least-squares error best fit between the modeled and measured wind vectors along the AA-539 flight path. The model indicates that the flowfield contains some significant elements and vortices not previously reported. The alternating direction of rotation of the vortices in the model suggests a microburst structure based on a von Karman vortex street rather than on a Kelvin-Helmholtz instability. The model also indicates that the reconstructed wind-vector data contain a time lag of at least one second in the horizontal winds.

Grantham, Walter J.↗

A DFW microburst model based on DL-191 data

Analysis of the August 2, 1985 crash of an L-1011 jumbo jet (DL-191), on approach to the Dallas-Ft. Worth International Airport (DFW) in a thunderstorm indicates that the severe windshear microburst which caused the crash was composed not only of a strong downflow and outflow but also included several large-scale vortex rings entrained in the flowfield. A detailed two-dimensional model of the DFW microburst, based on data from DL-191, is presented. The model uses wind vector and flight path data reconstructed by NASA-Ames Research Center from the L-1011's digital flight data records and ground-based radar measurements. The model was developed using a combination of interactive graphics and a least-square-error best fit between the modeled and measured wind vectors along the DL-191 flight path. The model confirms a microburst structure based on a von Karman vortex street rather than on a Kelvin-Helmholtz instability and indicates that the L-1011 flew near or directly through several strong vortices. In addition, the results also confirm that the reconstructed wind vector data contain a time lag in the horizontal winds. For the DL-191 this lag is approximately three seconds.

Grantham, Walter J.↗

REQUEST: A Recursive QUEST Algorithm for Sequential Attitude Determination

In order to find the attitude of a spacecraft with respect to a reference coordinate system, vector measurements are taken. The vectors are pairs of measurements of the same generalized vector, taken in the spacecraft body coordinates, as well as in the reference coordinate system. We are interested in finding the best estimate of the transformation between these coordinate system.s The algorithm called QUEST yields that estimate where attitude is expressed by a quarternion. Quest is an efficient algorithm which provides a least squares fit of the quaternion of rotation to the vector measurements. Quest however, is a single time point (single frame) batch algorithm, thus measurements that were taken at previous time points are discarded. The algorithm presented in this work provides a recursive routine which considers all past measurements. The algorithm is based on on the fact that the, so called, K matrix, one of whose eigenvectors is the sought quaternion, is linerly related to the measured pairs, and on the ability to propagate K. The extraction of the appropriate eigenvector is done according to the classical QUEST algorithm. This stage, however, can be eliminated, and the computation simplified, if a standard eigenvalue-eigenvector solver algorithm is used. The development of the recursive algorithm is presented and illustrated via a numerical example.

Bar-Itzhack, Itzhack Y.↗

Adjoint Methods for Adjusting Three-Dimensional Atmosphere and Surface Properties to Fit Multi-Angle/Multi-Pixel Polarimetric Measurements

This paper derives an efficient procedure for using the three-dimensional (3D) vector radiative transfer equation (VRTE) to adjust atmosphere and surface properties and improve their fit with multi-angle/multi-pixel radiometric and polarimetric measurements of scattered sunlight. The proposed adjoint method uses the 3D VRTE to compute the measurement misfit function and the adjoint 3D VRTE to compute its gradient with respect to all unknown parameters. In the remote sensing problems of interest, the scalar-valued misfit function quantifies agreement with data as a function of atmosphere and surface properties, and its gradient guides the search through this parameter space. Remote sensing of the atmosphere and surface in a three-dimensional region may require thousands of unknown parameters and millions of data points. Many approaches would require calls to the 3D VRTE solver in proportion to the number of unknown parameters or measurements. To avoid this issue of scale, we focus on computing the gradient of the misfit function as an alternative to the Jacobian of the measurement operator. The resulting adjoint method provides a way to adjust 3D atmosphere and surface properties with only two calls to the 3D VRTE solver for each spectral channel, regardless of the number of retrieval parameters, measurement view angles or pixels. This gives a procedure for adjusting atmosphere and surface parameters that will scale to the large problems of 3D remote sensing. For certain types of multi-angle/multi-pixel polarimetric measurements, this encourages the development of a new class of three-dimensional retrieval algorithms with more flexible parametrizations of spatial heterogeneity, less reliance on data screening procedures, and improved coverage in terms of the resolved physical processes in the Earth's atmosphere.

linearization↗

Mars physical parameters as determined from Mariner 9 observations of the natural satellites

Derivation of values for the Mars gravity field, as well as the mass and spin axis direction, by processing combined radio and optical data taken with the Mariner 9 spacecraft. The optical data consists of 62 TV photographs of Phobos and Deimos taken during the Mars orbiting phase of Mariner 9. The radio data consists of apoapsis state vectors obtained from 195 one-revolution fits of the Doppler data. A first-order analytic theory is used to model both the satellites and the spacecraft's motion. However, the shallow resonance condition of Mariner 9's orbit necessitated development of a partial second-order theory to model the long-period perturbations which derive from Mars even-order tesseral harmonics. Use of an analytic theory allows rapid and inexpensive processing of both data types simultaneously. The results are in good agreement with those previously published from Mariner 9 Doppler and landmark data.

Born, G. H.↗

Comparing the jerk with other global models of the geomagnetic field from 1960 to 1978

About 3300 satellite values of geomagnetic intensity and about 700 observatory values of annual mean magnetic vector components from 1960 to 1978 were fitted by three global models of the geomagnetic field B. Each model includes a spatially constant external field whose time dependence is a constant plus another constant times the Dst index, and each model accepts a time-independent station correction at each observatory. The time dependence of the internal Gauss coefficients is either cubic, quintic, or biquadratic (two independent quadratics, one before and one after January 1, 1970); and g1(0) also has an induced term proportional to the Dst index. The rms residual of the data fit is the same for the cubic and biquadratic models and insignificantly smaller for the quintic model. The quintic and biquadratic models have 1164 adjustable parameters, and the cubic has 1038. At a high level of significance the parameters of the best fitting biquadratic rule out a physical model for the magnetic impulse of 1969 in which the level surfaces of electrical conductivity in the lower mantle are approximately spherical, and the radial magnetic field at the core-mantle boundary goes from one quadratic time dependence to another in a year or less.

Backus, G. E.↗

Full potential methods for analysis/design of complex aerospace configurations

The steady form of the full potential equation, in conservative form, is employed to analyze and design a wide variety of complex aerodynamic shapes. The nonlinear method is based on the theory of characteristic signal propagation coupled with novel flux biasing concepts and body-fitted mapping procedures. The resulting codes are vectorized for the CRAY XMP and the VPS-32 supercomputers. Use of the full potential nonlinear theory is demonstrated for a single-point supersonic wing design and a multipoint design for transonic maneuver/supersonic cruise/maneuver conditions. Achievement of high aerodynamic efficiency through numerical design is verified by wind tunnel tests. Other studies reported include analyses of a canard/wing/nacelle fighter geometry.

Shankar, Vijaya↗

The NASA OBPG 2020 On-Orbit Calibration of SNPP VIIRS for Ocean Color Applications

The NASA Ocean Biology Processing Group (OBPG) has continued monitoring the SNPP VIIRS on-orbit calibration for bands M1-M11 over its mission to optimize the calibration for ocean color applications. The OBPG has recently implemented several changes to the calibration scheme: using solar-derived f-factors to detrend the lunar observations; using long-term exponentials of time as basis vectors (along with libration angles) for radiometric fits to any resulting lunar temporal drifts; deriving gain adjustments to the solar f-factors from these exponentials; and deriving gain adjustments due to modulated RSRs outside of the solar/lunar calibration using TOA reference spectra. These calibration changes minimize the impact of uncertainties in any one component of the calibration on the derived f-factors. The final f-factors incorporate VIIRS solar diffuser measurements, h-factor BRDF corrections, lunar-derived gains, and modulated RSR gains. The combined BRDF corrections, lunar gain adjustments, and mRSR gain adjustments define effective h-factors for each band. The improvements in the on-orbit calibration are validated by evaluation of globally-derived anomaly plots of remote sensing reflectance for the ocean color bands. The ultimate goal of the OBPG calibration effort is incorporation of a consistent SNPP VIIRS ocean color data set into the NASA multi-mission ocean color climate data record.

VIIRS↗

Adjoint Methods for Adjusting Three-Dimensional Atmosphere and Surface Properties to Fit Multi-Angle Multi-Pixel Polarimetric Measurements

This paper derives an efficient procedure for using the three-dimensional (3D) vector radiative transfer equation (VRTE) to adjust atmosphere and surface properties and improve their fit with multi-angle/multi-pixel radiometric and polarimetric measurements of scattered sunlight. The proposed adjoint method uses the 3D VRTE to compute the measurement misfit function and the adjoint 3D VRTE to compute its gradient with respect to all unknown parameters. In the remote sensing problems of interest, the scalar-valued misfit function quantifies agreement with data as a function of atmosphere and surface properties, and its gradient guides the search through this parameter space. Remote sensing of the atmosphere and surface in a three-dimensional region may require thousands of unknown parameters and millions of data points. Many approaches would require calls to the 3D VRTE solver in proportion to the number of unknown parameters or measurements. To avoid this issue of scale, we focus on computing the gradient of the misfit function as an alternative to the Jacobian of the measurement operator. The resulting adjoint method provides a way to adjust 3D atmosphere and surface properties with only two calls to the 3D VRTE solver for each spectral channel, regardless of the number of retrieval parameters, measurement view angles or pixels. This gives a procedure for adjusting atmosphere and surface parameters that will scale to the large problems of 3D remote sensing. For certain types of multi-angle/multi-pixel polarimetric measurements, this encourages the development of a new class of three-dimensional retrieval algorithms with more flexible parametrizations of spatial heterogeneity, less reliance on data screening procedures, and improved coverage in terms of the resolved physical processes in the Earth׳s atmosphere.

adjoint methods↗