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 109 records · Page 6

Three-dimensional inviscid analysis of the scramjet inlet flow field

A computer code has been developed to analyze the inviscid flow field in a supersonic combustion ramjet (scramjet) inlet. The code uses the three-dimensional Euler equations in full conservation form to describe the inlet flow. An algebraic numerical coordinate transformation is used to generate a set of boundary-fitted curvilinear coordinates. The governing equations are solved by a time-asymptotic, unsplit, two-step, finite-difference method. This method is highly efficient on the vector processing computers for which the current code is written. Detailed results are presented for two scramjet inlet configurations over a range of Mach numbers. The calculated results are compared with the available experimental and theoretical results.

Kumar, A.↗

Numerical experiments in geomagnetic modeling

Numerical tests were made, using least squares fitting of a spherical harmonic model, to a selection of Magsat data to determine the practical limits of this technique with modern computers. The resulting (M102189) model, whose coefficients were adjusted up to n = 50, was compared with M07AV6, a previous model which used least squares (on vector data) for coefficients up to n = 29, and Gauss-Legendre quadrature (on Z residuals) to adjust the coefficients up to n = 63. For the new least squares adjustment to n = 50 a condition number of 115 was obtained for the solution matrix, with a resulting precision of 11 significant figures. The M102189 model shows a lower and more Gaussian residual distribution than did M07AV6, though the Gaussian envelope fits to the residual distributions, even for the scalar field, gives "standard deviations' never lower than 6 nT, a factor of three higher than the estimated Magsat observational errors. Ionospheric currents are noted to have a significant effect on the coefficients of the internal potential functions.

Cain, Joseph C.↗

The Z3 model of Saturns magnetic field and the Pioneer 11 vector helium magnetometer observations

Magnetic field observations obtained by the Pioneer 11 vector helium magnetometer are compared with the Z(sub 3) model magnetic field. These Pioneer 11 observations, obtained at close-in radial distances, constitute an important and independent test of the Z(sub 3) zonal harmonic model, which was derived from Voyager 1 and Voyager 2 fluxgate magnetometer observations. Differences between the Pioneer 11 magnetometer and the Z(sub 3) model field are found to be small (approximately 1%) and quantitatively consistent with the expected instrumental accuracy. A detailed examination of these differences in spacecraft payload coordinates shows that they are uniquely associated with the instrument frame of reference and operation. A much improved fit to the Pioneer 11 observations is obtained by rotation of the instrument coordinate system about the spacecraft spin axis by 1.4 degree. With this adjustment, possibly associated with an instrumental phase lag or roll attitude error, the Pioneer 11 vector helium magnetometer observations are fully consistent with the Voyager Z(sub 3) model.

Connerney, J. E. P.↗

Neural-Network Simulator

F77NNS (FORTRAN 77 Neural Network Simulator) computer program simulates popular back-error-propagation neural network. Designed to take advantage of vectorization when used on computers having this capability, also used on any computer equipped with ANSI-77 FORTRAN Compiler. Problems involving matching of patterns or mathematical modeling of systems fit class of problems F77NNS designed to solve. Program has restart capability so neural network solved in stages suitable to user's resources and desires. Enables user to customize patterns of connections between layers of network. Size of neural network F77NNS applied to limited only by amount of random-access memory available to user.

Mitchell, Paul H.↗

The Z3 model of Saturn's magnetic field and the Pioneer 11 vector helium magnetometer observations

Magnetic field observations obtained by the Pioneer 11 vector helium magnetometer are compared with the Z(sub 3) model magnetic field. These Pioneer 11 observations, obtained at close-in radial distances, constitute an important and independent test of the Z(sub 3) zonal harmonic model, which was derived from Voyager 1 and Voyager 2 fluxgate magnetometer observations. Differences between the Pioneer 11 magnetometer and the Z(sub 3) model field are found to be small (approximately 1 percent) and quantitatively consistent with the expected instrumental accuracy. A detailed examination of these differences in spacecraft payload coordinates shows that they are uniquely associated with the instrument frame of reference and operation. A much improved fit to the Pioneer 11 observations is obtained by rotation of the instrument coordinate system about the spacecraft spin axis by 1.4 degree. With this adjustment, possibly associated with an instrumental phase lag or roll attitude error, the Pioneer 11 vector helium magnetometer observations are fully consistent with the Voyager Z(sub 3) model.

Connerney, J. E. P.↗

Multi-spacecraft observations of magnetopause surface waves - ISEE 1 and 2 determinations of amplitude, wavelength and period

The multispacecraft ISEE mission made it possible to study propagation of surface waves on the magnetopause. One case studied in detail was near local noon at a latitude of 24.4 deg. The inferred instantaneous magnetopause normal vectors oscillated about the model magnetopause normal vector in a manner expected for a wave propagating along the magnetopause from the subsolar region. By assuming a model monochromatic wave and a constant magnetopause thickness, the best-fit that yields the three wave parameters( amplitude, wavelength, and period) is determined.

Song, PU↗

In Situ Surface Characterization

Operation of in situ space assets, such as rovers and landers, requires operators to acquire a thorough understanding of the environment surrounding the spacecraft. The following programs help with that understanding by providing higher-level information characterizing the surface, which is not immediately obvious by just looking at the XYZ terrain data. This software suite covers three primary programs: marsuvw, marsrough, and marsslope, and two secondary programs, which together use XYZ data derived from in situ stereo imagery to characterize the surface by determining surface normal, surface roughness, and various aspects of local slope, respectively. These programs all use the Planetary Image Geometry (PIG) library to read mission-specific data files. The programs themselves are completely multimission; all mission dependencies are handled by PIG. The input data consists of images containing XYZ locations as derived by, e.g., marsxyz. The marsuvw program determines surface normals from XYZ data by gathering XYZ points from an area around each pixel and fitting a plane to those points. Outliers are rejected, and various consistency checks are applied. The result shows the orientation of the local surface at each point as a unit vector. The program can be run in two modes: standard, which is typically used for in situ arm work, and slope, which is typically used for rover mobility. The difference is primarily due to optimizations necessary for the larger patch sizes in the slope case. The marsrough program determines surface roughness in a small area around each pixel, which is defined as the maximum peak-to-peak deviation from the plane perpendicular to the surface normal at that pixel. The marsslope program takes a surface normal file as input and derives one of several slope-like outputs from it. The outputs include slope, slope rover direction (a measure of slope radially away from the rover), slope heading, slope magnitude, northerly tilt, and solar energy (compares the slope with the Sun s location at local noon). The marsuvwproj program projects a surface normal onto an arbitrary plane in space, resulting in a normalized 3D vector, which is constrained to lie in the plane. The marsuvwrot program rotates the vectors in a surface normal file, generating a new surface normal file. It also can change coordinate systems for an existing surface normal file. While the algorithms behind this suite are not particularly unique, what makes the programs useful is their integration into the larger in situ image processing system via the PIG library. They work directly with space in situ data, understanding the appropriate image metadata fields and updating them properly. The secondary programs (marsuvwproj, marsuvwrot) were originally developed to deal with anomalous situations on Opportunity and Spirit, respectively, but may have more general applicability.

Deen, Robert G.↗

The Extreme Spin of the Black Hole in Cygnus X-1

The compact primary in the X-ray binary Cygnus X-1 was the first black hole to be established via dynamical observatIOns. We have recently determined accurate values for its mass and distance, and for the orbital inclination angle of the binary. Building on these.results, which are based on our favored (asynchronous) dynamical model, we have measured the radius of the inner edge of the black hole's accretion disk by fitting its thermal continuum.spectrum to a fully relativistic model of a thin accretion disk. Assuming that the spin axis of the black hole is aligned with the orbital angular momentum vector, we have determined that Cygnus X-I contains a near-extreme Kerr black hole with a spin parameter a* > 0.95 (3(sigma)). For a less probable (synchronous) dynamIcal model, we find a* > 0.92 (3(sigma)). In our analysis, we include the uncertainties in black hole mass orbital inclination angle and distance, and we also include the uncertainty in the calibration of the absolute flux via the Crab. These four sources of uncertainty totally dominate the error budget. The uncertainties introduced by the thin-disk model we employ are particularly small in this case given the extreme spin of the black hole and the disk's low luminosity.

Gou, Lijun↗

An approximate factorization solution of the Navier-Stokes equations for transonic flow using body-fitted coordinates with application to NACA 64A010 airfoils

The implementation of the approximate factorization algorithm and its ability to efficiently and accurately describe transonic flow about an NACA 64A010 airfoil section is examined. The approximate factorization algorithm is developed from the nondimensional, conservative, vectorized Navier-Stokes equations expressed in curvilinear coordinates. Equations of state and transport coefficient relations appropriate to atmospheric air are appended to close the system of partial differential equations. An algebraic turbulence model is also incorporated into the equation set. This algorithm was verified by investigating the flow about an NACA 64A010 airfoil at 0, 2, and 3.5 deg angle of attack for free-stream conditions of 2,000,000 Reynolds number and 0.8 Mach number. Overall results were in good qualitative agreement with wind tunnel data sets. However, while nondimensional times of six were attained, numerical difficulties prevented any case from reaching a true steady state.

Copper, G. K.↗

Improved image decompression for reduced transform coding artifacts

The perceived quality of images reconstructed from low bit rate compression is severely degraded by the appearance of transform coding artifacts. This paper proposes a method for producing higher quality reconstructed images based on a stochastic model for the image data. Quantization (scalar or vector) partitions the transform coefficient space and maps all points in a partition cell to a representative reconstruction point, usually taken as the centroid of the cell. The proposed image estimation technique selects the reconstruction point within the quantization partition cell which results in a reconstructed image which best fits a non-Gaussian Markov random field (MRF) image model. This approach results in a convex constrained optimization problem which can be solved iteratively. At each iteration, the gradient projection method is used to update the estimate based on the image model. In the transform domain, the resulting coefficient reconstruction points are projected to the particular quantization partition cells defined by the compressed image. Experimental results will be shown for images compressed using scalar quantization of block DCT and using vector quantization of subband wavelet transform. The proposed image decompression provides a reconstructed image with reduced visibility of transform coding artifacts and superior perceived quality.

Orourke, Thomas P.↗

Recent update of the RPLUS2D/3D codes

The development of the RPLUS2D/3D codes is summarized. These codes utilize LU algorithms to solve chemical non-equilibrium flows in a body-fitted coordinate system. The motivation behind the development of these codes is the need to numerically predict chemical non-equilibrium flows for the National AeroSpace Plane Program. Recent improvements include vectorization method, blocking algorithms for geometric flexibility, out-of-core storage for large-size problems, and an LU-SW/UP combination for CPU-time efficiency and solution quality.

Tsai, Y.-L. Peter↗

The Extreme Spin of the Black Hole in Cygnus X-1

The compact primary in the X-ray binary Cygnus X-1 was the first black hole to be established via dynamical observations. We have recently determined accurate values for its mass and distance, and for the orbital inclination angle of the binary. Building on these results, which are based on our favored (asynchronous) dynamical model, we have measured the radius of the inner edge of the black hole s accretion disk by fitting its thermal continuum spectrum to a fully relativistic model of a thin accretion disk. Assuming that the spin axis of the black hole is aligned with the orbital angular momentum vector, we have determined that Cygnus X-1 contains a near-extreme Kerr black hole with a spin parameter a* > 0.95 (3(sigma)). For a less probable (synchronous) dynamical model, we find a. > 0.92 (3 ). In our analysis, we include the uncertainties in black hole mass, orbital inclination angle, and distance, and we also include the uncertainty in the calibration of the absolute flux via the Crab. These four sources of uncertainty totally dominate the error budget. The uncertainties introduced by the thin-disk model we employ are particularly small in this case given the extreme spin of the black hole and the disk s low luminosity.

Gou, Lijun↗

Roll tracking effects of G-vector tilt and various types of motion washout

In a dogfight scenario, the task was to follow the target's roll angle while suppressing gust disturbances. All subjects adopted the same behavioral strategies in following the target while suppressing the gusts, and the MFP-fitted math model response was generally within one data symbol width. The results include the following: (1) comparisons of full roll motion (both with and without the spurious gravity tilt cue) with the static case. These motion cues help suppress disturbances with little net effect on the visual performance. Tilt cues were clearly used by the pilots but gave only small improvement in tracking errors. (2) The optimum washout (in terms of performance close to real world, similar behavioral parameters, significant motion attenuation (60 percent), and acceptable motion fidelity) was the combined attenuation and first-order washout. (3) Various trends in parameters across the motion conditions were apparent, and are discussed with respect to a comprehensive model for predicting adaptation to various roll motion cues.

Jex, H. R.↗

Modification to shock fitting program

A modified form of the Lepping - Argentiero single spacecraft, shock normal determination procedure is presented. The modified method incorporates a simple predictor-corrector algorithm which allows a faster convergence rate and the use of average values of the parameters for the starting vector.

Acuna, M. H.↗

Numerical solution of the Navier-Stokes equations for laminar, transonic flows

An implicit finite difference solution of the Navier-Stokes equations yielded time histories of the transonic laminar flow development about a circular cylinder and NACA-0018 airfoil. Reynolds numbers ranged from those corresponding to purely laminar flow to those corresponding to significant turbulence in the boundary layer. Body thermal conditions of an adiabatic wall and a specified body temperature were considered. Versatility in treating arbitrary bodies was incorporated by using numerically generated, body-fitted coordinate transformations. Solution of the simultaneous difference equations for the dependent variables was obtained using an accelerated Gauss-Seidel iterative scheme. Computational results are presented in the form of velocity vector fields, Mach number contours, aerodynamic coefficients, heat transfer rates at the body surface, and body temperature distributions. Truncation analyses of first and second derivative difference approximations resulted in general criteria for numerically generated coordinates so that flow near a body is more accurately represented.

Turner, L., III↗

Implicit multigrid algorithms for the three-dimensional flux split Euler equations

The full approximation scheme multigrid method is applied to several implicit flux-split algorithms for solving the three-dimensional Euler equations in a body fitted coordinate system. Each uses a variation of approximate factorization and is implemented in a finite volume formulation. The algorithms are all vectorizable with little or no scalar computations required. The flux vectors are split into upwind components using both the splittings of Steger-Warming and Van Leer. Results comparing pressure distributions with experimental data using both splitting types are shown. 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. Using the multigrid method on both subsonic and transonic wing calculations, the final lift coefficient is obtained to within 0.1 percent of its final value in a few as 15 cycles for a mesh with over 210,000 points. A spectral radius of 0.89 is achieved for both subsonic and transonic flow over the ONERA M6 wing while a spectral radius of 0.83 is obtained for supersonic flow over an analytically defined forebody. Results compared with experiment for all cases show good agreement.

Anderson, W. K.↗

Interactive Spectral Analysis and Computation (ISAAC)

Isaac is a task in the NSO external package for IRAF. A descendant of a FORTRAN program written to analyze data from a Fourier transform spectrometer, the current implementation has been generalized sufficiently to make it useful for general spectral analysis and other one dimensional data analysis tasks. The user interface for Isaac is implemented as an interpreted mini-language containing a powerful, programmable vector calculator. Built-in commands provide much of the functionality needed to produce accurate line lists from input spectra. These built-in functions include automated spectral line finding, least squares fitting of Voigt profiles to spectral lines including equality constraints, various filters including an optimal filter construction tool, continuum fitting, and various I/O functions.

Lytle, D. M.↗

Modification to shock fitting program

A modified form of the Lepping-Argentiero single spacecraft, shock normal determination procedure is presented. The modified method incorporates a simple predictor-corrector algorithm which allows a faster convergence rate and the use of average values of the parameters for the starting vector.

Acuna, M. H.↗