Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Jacobian”

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

Automatic guidance and control of a transport aircraft during a helical landing approach

A linear optimal regulator theory was applied to a nonlinear simulation of a transport aircraft performing a helical landing approach. A closed-form expression for the quasi-steady nominal flight path is presented along with the method for determining the corresponding constant nominal control inputs. The Jacobian matrices and the weighting matrices in the cost functional were time varying. A method of solving for the optimal feedback gains is reviewed. The control system was tested on several alternative landing approaches using both 3 deg and 6 deg flight path angles. On each landing approach, the aircraft was subjected to large random initial-state errors and to randomly directed crosswinds. The system was also tested for sensitivity to changes in the parameters of the aircraft and of the atmosphere. Results indicate that performance of the optimal controller on all the 3 deg approaches is very good. The control system proved to be reasonably insensitive to parametric uncertainties. Performance is not as good on the 6 deg approaches. A modification to the 6 deg flight path was proposed for the purpose of improving performance.

Crawford, D. J.↗

On differential transformations between Cartesian and curvilinear (geodetic) coordinates

Differential transformations are developed between Cartesian and curvilinear orthogonal coordinates. Only matrix algebra is used for the presentation of the basic concepts. After defining the reference systems used the rotation (R), metric (H), and Jacobian (J) matrices of the transformations between cartesian and curvilinear coordinate systems are introduced. A value of R as a function of H and J is presented. Likewise an analytical expression for J(-1) as a function of H(-2) and R is obtained. Emphasis is placed on showing that differential equations are equivalent to conventional similarity transformations. Scaling methods are discussed along with ellipsoidal coordinates. Differential transformations between elipsoidal and geodetic coordinates are established.

Soler, T.↗

Effects of radiation pressure on the equipotential surfaces in X-ray binaries

Equipotential surfaces incorporating the effect of radiation pressure were computed for the X-ray binaries Cen X-3, Cyg X-1 = HDE 226868, Vela XR-1 = 3U 0900-40 = HD 77581, and 3U 1700-37 = HD 153919. The topology of the equipotential surfaces is significantly affected by radiation pressure. In particular, the so-called critical Roche (Jacobian) lobes, the traditional figure 8's, do not exist. The effects of these results on modeling X-ray binaries are discussed.

Kondo, Y.↗

An improved method for the prediction of completely three-dimensional aerodynamic load distributions of configurations with leading edge vortex separation

The application of a higher-order subsonic potential flow panel method to the solution of three-dimensional flow about wing and wing-body combinations with leading-edge vortex separation is presented. The governing equations are the linear flow differential equation and nonlinear boundary conditions which require that the flow be parallel to the wing and body surfaces and that the free vortex sheet, springing from the leading and trailing edges, be aligned with the local flow and support no pressure jump. The vortex core is modeled as a simple line vortex which receives vorticity from the free sheet through a connecting sheet. The Kutta condition is imposed on all appropriate edges of the wing. This set of nonlinear equations is solved by an iterative procedure. The Goethert rule accounts for compressibility. The method has been programmed for the CDC 6600. Delta wings, gothic wings, arrow wings, cambered wings, and wing with body have been analyzed. Initial studies involving variations of panel density, vortex sheet sizing, Jacobian update, and initial geometry demonstrate that the present method generally exhibits good convergence characteristics.

Rubbert, p. E.↗

Mass flow in the O7f binary UW Canis Majoris. II

The far-ultraviolet spectrum of UW CMa was obtained from the telescope spectrometer on the Copernicus satellite at about phase 0.25 in order to supplement earlier observations of the same kind made near phase 0.75. Analysis of the data shows that the high-temperature gas giving rise to the P Cygni features in the far-ultraviolet spectrum is located primarily around the entire binary system and does not share in the motion of either component. The effects of radiation pressure on the Roche (Jacobian) equipotential surfaces and in generating a stellar wind are discussed.

Mccluskey, G. E., Jr.↗

RMS massless arm dynamics capability in the SVDS

The equations of motion for the remote manipulator system, assuming that the masses and inertias of the arm can be neglected, are developed for implementation into the space vehicle dynamics simulation (SVDS) program for the Orbiter payload system. The arm flexibility is incorporated into the equations by the computation of flexibility terms for use in the joint servo model. The approach developed in this report is based on using the Jacobian transformation matrix to transform force and velocity terms between the configuration space and the task space to simplify the form of the equations.

Flanders, H. A.↗

Mass flow in close binary systems

The manner of mass flow in close binary systems is examined with a special view to the role of the so-called critical Roche (or Jacobian) lobe, taking into consideration relevant physical conditions such as radiation pressure that may affect the restricted three-body problem treatment. The mass does not necessarily flow from component one to component two through the L1 point to form a gaseous ring surrounding the latter. These considerations are applied to X-ray binaries with early-type optical components, such as Cyg X-1 (HDE 226868) and 3U 1700 - 37 (HD 153919). In the two bright close binary systems Beta Lyr and UW CMa, which are believed to be undergoing dynamic mass transfer, recent Copernicus observations show that the gas giving rise to the prominent ultraviolet emission lines surrounds the entire binary system rather than merely component two. Implications of these observations are also discussed.

Kondo, Y.↗

A high order accurate finite element algorithm for high Reynolds number flow prediction

A Galerkin-weighted residuals formulation is employed to establish an implicit finite element solution algorithm for generally nonlinear initial-boundary value problems. Solution accuracy, and convergence rate with discretization refinement, are quantized in several error norms, by a systematic study of numerical solutions to several nonlinear parabolic and a hyperbolic partial differential equation characteristic of the equations governing fluid flows. Solutions are generated using selective linear, quadratic and cubic basis functions. Richardson extrapolation is employed to generate a higher-order accurate solution to facilitate isolation of truncation error in all norms. Extension of the mathematical theory underlying accuracy and convergence concepts for linear elliptic equations is predicted for equations characteristic of laminar and turbulent fluid flows at nonmodest Reynolds number. The nondiagonal initial-value matrix structure introduced by the finite element theory is determined intrinsic to improved solution accuracy and convergence. A factored Jacobian iteration algorithm is derived and evaluated to yield a consequential reduction in both computer storage and execution CPU requirements while retaining solution accuracy.

Baker, A. J.↗

Description of a computer program and numerical techniques for developing linear perturbation models from nonlinear systems simulations

A numerical technique was developed which generates linear perturbation models from nonlinear aircraft vehicle simulations. The technique is very general and can be applied to simulations of any system that is described by nonlinear differential equations. The computer program used to generate these models is discussed, with emphasis placed on generation of the Jacobian matrices, calculation of the coefficients needed for solving the perturbation model, and generation of the solution of the linear differential equations. An example application of the technique to a nonlinear model of the NASA terminal configured vehicle is included.

Dieudonne, J. E.↗

Coefficient matrices for implicit finite difference solution of the inviscid fluid conservation law equations

Although the Navier-Stokes equations describe most flows of interest in aerodynamics, the inviscid conservation law equations may be used for small regions with viscous forces. Thus, Euler equations and several time-accurate finite difference procedures, explicit and implicit, are discussed. Although implicit techniques require more computational work, they permit larger time steps to be taken without instability. It is noted that the Jacobian matrices for Euler equations in conservation-law form have certain eigenvalue-eigenvector properties which may be used to construct conservative-form coefficient matrices. This reduces the computation time of several implicit and semiimplicit schemes. Extensions of the basic approach to other areas are suggested.

Steger, J. L.↗

A high order accurate numerical solution algorithm for turbulent boundary layer flow

A fourth-order accurate numerical solution algorithm is derived using finite element interpolation theory for the non-linear parabolic equations governing turbulent boundary layer flow including a two-equation turbulence closure model. The results of carefully controlled numerical experiments firmly quantize for the first time performance differences between finite element and finite difference solution methodology for this type of equation. The developed algorithm takes advantage of the apparent semi-analytical formulational procedure, in establishment of a single, retarded-evaluation Jacobian matrix iterative solution algorithm. Numerical results document performance of solution economy features in terms of computer requirements and solution accuracy. The developed algorithm should find wide application in aerodynamics analysis.

Soliman, M. O.↗

Matrix perturbation for analytical model improvement

A matrix perturbation method is proposed to calculate the Jacobian matrix and to compute the new eigendata for the parameter estimation procedure. The advantages of the method are the applicability to large complex structures without knowing the analytical expressions for the mass and stiffness matrices, and a cost effective approach for the re-computation of the eigendata. This method also allows the use of other measurements such as modal forces, kinetic energy distribution, and strain energy distributions in the estimation procedure. A realistic sample problem is presented to demonstrate the effectiveness of the proposed method.

Chen, J. C.↗

Computer coordination of limb motion for a three-legged walking robot

Coordination of the limb motion of a vehicle which could perform assembly and maintenance operations on large structures in space is described. Manipulator kinematics and walking robots are described. The basic control scheme of the robot is described. The control of the individual arms are described. Arm velocities are generally described in Cartesian coordinates. Cartesian velocities are converted to joint velocities using the Jacobian matrix. The calculation of a trajectory for an arm given a sequence of points through which it is to pass is described. The free gait algorithm which controls the lifting and placing of legs for the robot is described. The generation of commanded velocities for the robot, and the implementation of those velocities by the algorithm are discussed. Suggestions for further work in the area of robot legged locomotion are presented.

Klein, C. A.↗

Computer-aided analysis of nonlinear problems in transport phenomena

The paper describes algorithms for equilibrium and steady-state problems with coefficients in the expansions derived by the Galerkin weighted residual method and calculated from the resulting sets of nonlinear algebraic equations by the Newton-Raphson method. Initial approximations are obtained from nearby solutions by continuation techniques as parameters are varied. The Newton-Raphson technique is preferred because the Jacobian of the solution is useful for continuation, for analyzing the stability of solutions, for detecting bifurcation of solution families, and for computing asymptotic estimates of the effects on any solution of small changes in parameters, boundary conditions, and boundary shape.

Brown, R. A.↗

Two-boundary grid generation for the solution of the three dimensional compressible Navier-Stokes equations

A grid generation technique called the two boundary technique is developed and applied for the solution of the three dimensional Navier-Stokes equations. The Navier-Stokes equations are transformed from a cartesian coordinate system to a computational coordinate system, and the grid generation technique provides the Jacobian matrix describing the transformation. The two boundary technique is based on algebraically defining two distinct boundaries of a flow domain and the distribution of the grid is achieved by applying functions to the uniform computational grid which redistribute the computational independent variables and consequently concentrate or disperse the grid points in the physical domain. The Navier-Stokes equations are solved using a MacCormack time-split technique. Grids and supersonic laminar flow solutions are obtained for a family of three dimensional corners and two spike-nosed bodies.

Smith, R. E.↗

Visual tracking in stereo

A method is described for visual object tracking by a computer vision system using TV cameras and special low-level image processing hardware. The tracker maintains an internal model of the location, orientation, and velocity of the object in three-dimensional space. This model is used to predict where features of the object will lie on the two-dimensional images produced by stereo TV cameras. The differences in the locations of features in the two-dimensional images as predicted by the internal model and as actually seen create an error signal in the two-dimensional representation. This is multiplied by a generalized inverse Jacobian matrix to deduce the error in the internal model. The procedure repeats to update the internal model of the object's location, orientation and velocity continuously.

Saund, E.↗

Grid spacing control with variation diminishing splines

Methods used to specify and control two and three dimensional grids on which numerical solutions of partial differential equations may be obtained were studied. While initially focusing on grid generation, the research will evolve into a consideration of the interaction of grid generation with the solution of a partial differential equation. The multisurface method of grid generation was used to continuously patch a grid onto an existing grid. In the resulting grid the elements of the Jacobian matrix must be continuous across the boundary between the original grid and the patched grid. Programs were written which accept as input the coordinates of the original grid and the desired new boundary, and then use the three-surface or four-surface version of the multisurface method to extend the original grid out to a new boundary.

Smith, P. W.↗