Search NASA⌕ Search

SEARCH · Search NASA

Results for “FIRST ORDER EQUATION”

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 91 records · Page 5

Minimal parameter solution of the orthogonal matrix differential equation

As demonstrated in this work, all orthogonal matrices solve a first order differential equation. The straightforward solution of this equation requires n sup 2 integrations to obtain the element of the nth order matrix. There are, however, only n(n-1)/2 independent parameters which determine an orthogonal matrix. The questions of choosing them, finding their differential equation and expressing the orthogonal matrix in terms of these parameters are considered. Several possibilities which are based on attitude determination in three dimensions are examined. It is shown that not all 3-D methods have useful extensions to higher dimensions. It is also shown why the rate of change of the matrix elements, which are the elements of the angular rate vector in 3-D, are the elements of a tensor of the second rank (dyadic) in spaces other than three dimensional. It is proven that the 3-D Gibbs vector (or Cayley Parameters) are extendable to other dimensions. An algorithm is developed employing the resulting parameters, which are termed Extended Rodrigues Parameters, and numerical results are presented of the application of the algorithm to a fourth order matrix.

Bar-Itzhack, Itzhack Y.↗

Minimal parameter solution of the orthogonal matrix differential equation

As demonstrated in this work, all orthogonal matrices solve a first order differential equation. The straightforward solution of this equation requires n sup 2 integrations to obtain the element of the nth order matrix. There are, however, only n(n-1)/2 independent parameters which determine an orthogonal matrix. The questions of choosing them, finding their differential equation and expressing the orthogonal matrix in terms of these parameters are considered. Several possibilities which are based on attitude determination in three dimensions are examined. It is shown that not all 3-D methods have useful extensions to higher dimensions. It is also shown why the rate of change of the matrix elements, which are the elements of the angular rate vector in 3-D, are the elements of a tensor of the second rank (dyadic) in spaces other than three dimensional. It is proven that the 3-D Gibbs vector (or Cayley Parameters) are extendable to other dimensions. An algorithm is developed employing the resulting parameters, which are termed Extended Rodrigues Parameters, and numerical results are presented of the application of the algorithm to a fourth order matrix.

Baritzhack, Itzhack Y.↗

Minimal parameter solution of the orthogonal matrix differential equation

As demonstrated in this work, all orthogonal matrices solve a first order differential equation. The straightforward solution of this equation requires n sup 2 integrations to obtain the element of the nth order matrix. There are, however, only n(n-1)/2 independent parameters which determine an orthogonal matrix. The questions of choosing them, finding their differential equation and expressing the orthogonal matrix in terms of these parameters are considered. Several possibilities which are based on attitude determination in three dimensions are examined. It is shown that not all 3-D methods have useful extensions to higher dimensions. It is also shown why the rate of change of the matrix elements, which are the elements of the angular rate vector in 3-D, are the elements of a tensor of the second rank (dyadic) in spaces other than three dimensional. It is proven that the 3-D Gibbs vector (or Cayley Parameters) are extendable to other dimensions. An algorithm is developed emplying the resulting parameters, which are termed Extended Rodrigues Parameters, and numerical results are presented of the application of the algorithm to a fourth order matrix.

Bar-Itzhack, Itzhack Y.↗

Fast methods incorporating direct elliptic solvers for nonlinear applications in fluid dynamics

Semidirect methods are discussed, their present role, as well as some developments for their application in computational fluid dynamics. A semidirect method is a computational scheme that uses a fast, direct, elliptic solver as the driving algorithm for the iterative solution of finite difference equations. Specific subtopics include: (1) direct Cauchy Riemann solvers for first order elliptic equations; (2) application of the semidirect method to the mixed elliptic hyperbolic problem of steady, inviscid transonic flow; and (3) the treatment of interior conditions, such as those on an airfoil or wing, in semidirect methods.

Martin, E. D.↗

Three-dimensional boundary layer calculation by a characteristic method

A numerical method for solving the three-dimensional boundary layer equations for bodies of arbitrary shape is presented. In laminar flows, the application domain extends from incompressible to hypersonic flows with the assumption of chemical equilibrium. For turbulent boundary layers, the application domain is limited by the validity of the mixing length model used. In order to respect the hyperbolic nature of the equations reduced to first order partial derivative terms, the momentum equations are discretized along the local streamlines using of the osculator tangent plane at each node of the body fitted coordinate system. With this original approach, it is possible to overcome the use of the generalized coordinates, and therefore, it is not necessary to impose an extra hypothesis about the regularity of the mesh in which the boundary conditions are given. By doing so, it is possible to limit, and sometimes to suppress, the pre-treatment of the data coming from an inviscid calculation. Although the proposed scheme is only semi-implicit, the method remains numerically very efficient.

Houdeville, R.↗

Mixed finite-difference scheme for analysis of simply supported thick plates.

A mixed finite-difference scheme is presented for the stress and free vibration analysis of simply supported nonhomogeneous and layered orthotropic thick plates. The analytical formulation is based on the linear, three-dimensional theory of orthotropic elasticity and a Fourier approach is used to reduce the governing equations to six first-order ordinary differential equations in the thickness coordinate. The governing equations possess a symmetric coefficient matrix and are free of derivatives of the elastic characteristics of the plate. In the finite difference discretization two interlacing grids are used for the different fundamental unknowns in such a way as to reduce both the local discretization error and the bandwidth of the resulting finite-difference field equations. Numerical studies are presented for the effects of reducing the interior and boundary discretization errors and of mesh refinement on the accuracy and convergence of solutions. It is shown that the proposed scheme, in addition to a number of other advantages, leads to highly accurate results, even when a small number of finite difference intervals is used.

Noor, A. K.↗

Analytical description of charged particle transport along arbitrary guiding field configurations

Aspects of charged particle transport theory have been studied by Earl (1973, 1974, 1976) in a number of investigations. The investigations were conducted under the assumption that particles propagate along a DC guiding field with negligible perpendicular transport and that pitch-angle scattering by random fields can be described by the Fokker-Planck formalism. The present study provides a generalized and simplified treatment of focused transport. The first step is taken with respect to a perturbation analysis by which the Boltzmann equation is solved in terms of successively improved approximations to its exact solution. The first approximation is the sum of isotropic pseudodiffusive and anisotropic supercoherent components whose coefficients must satisfy two coupled partial differential equations. The normalization constant K and the pseudodiffusive velocity, which appear in the first-order transport equations, are given by simple integrals over the pitch angle.

Earl, J. A.↗

Algorithms For Integrating Nonlinear Differential Equations

Improved algorithms developed for use in numerical integration of systems of nonhomogenous, nonlinear, first-order, ordinary differential equations. In comparison with integration algorithms, these algorithms offer greater stability and accuracy. Several asymptotically correct, thereby enabling retention of stability and accuracy when large increments of independent variable used. Accuracies attainable demonstrated by applying them to systems of nonlinear, first-order, differential equations that arise in study of viscoplastic behavior, spread of acquired immune-deficiency syndrome (AIDS) virus and predator/prey populations.

Freed, A. D.↗

The spectrum of steady state turbulent convection.

Based on Heisenberg's statistical theory of turbulence, a model for steady state turbulent convection is herein proposed, and on the basis of this model, equations for the energy spectrum for steady state turbulent convection are derived. The spectrum is obtained from the solution of a nonlinear integral equation. After the integral equation is brought into a universally valid nondimensional form, it is transformed into a nonlinear first order differential equation to be solved numerically, with the Rayleigh number appearing as the only parameter. The energy spectrum has a substantial deviation from the Kolmogoroff law, as a result of the buoyancy force acting on the rising and falling eddies. The presented theory may be applicable to convection in planetary and stellar atmospheres wherein the radiative heat transport is small.

Winterberg, F.↗

Nonlinear behavior of acoustic waves in combustion chambers

The nonlinear growth and limiting amplitude of acoustic waves in a combustion chamber are considered. A formal framework is provided within which practical problems can be treated with a minimum of effort and expense. The general conservation equations were expanded in two small parameters, one characterizing the mean flow field and one measuring the amplitude of oscillations, and then combined to yield a nonlinear inhomogeneous wave equation. The unsteady pressure and velocity fields were expressed as syntheses of the normal modes of the chamber, but with unknown time-varying amplitudes. This procedure yielded a representation of a general unsteady field as a system of coupled nonlinear oscillators. The system of nonlinear equations was treated by the method of averaging to produce a set of coupled nonlinear first order differential equations for the amplitudes and phases of the modes. The analysis is applicable to any combustion chamber. The most interesting applications are probably to solid rockets, liquid rockets, or thrust augmentors on jet engines.

Culick, F. E. C.↗

Nonlinear behavior of acoustic waves in combustion chambers. I, II

The general problem of the nonlinear growth and limiting amplitude of acoustic waves in a combustion chamber is treated in three parts: (1) the general conservation equations are expanded in two small parameters, and then combined to yield a nonlinear inhomogeneous wave equation, (2) the unsteady pressure and velocity fields are expressed as a synthesis of the normal modes of the chamber, but with unknown time-varying amplitudes, and (3) the system of nonlinear equations is treated by the method of averaging to produce a set of coupled nonlinear first order differential equations for the amplitudes and phases of the modes. This approximate analysis is applied to the investigation of the unstable motions in a solid propellant rocket engine and in a T burner.

Culick, F. E. C.↗

Numerical solutions of the complete Navier-Stokes equations

Using ideas from the kinetic theory, the Navier-Stokes equations are modified in such a way that they can be cast as a set of first order hyperbolic equations. This is achieved by incorporating time dependent terms into the definition of the stress tensor and the heat flux vectors. The boundary conditions are then determined from the theory of characteristics. Because the resulting equations reduce to the traditional Navier-Stokes equations when the steady state is reached, the present approach provides a straightforward scheme for the determination of inflow and outflow boundary conditions. The method is validated by comparing its predictions with known exact solutions of the steady Navier-Stokes equations.

Hassan, H. A.↗

Boundary conditions for the Navier-Stokes equations

Using ideas from kinetic theory, the Navier-Stokes equations are modified in such a way that they can be cast as a set of first order hyperbolic equations. This is achieved by incorporating time dependent terms into the definition of the stress tensor and the heat flux vectors. The boundary conditions are then determined from the theory of characteristics. Because the resulting equations reduce to the traditional Navier-Stokes equations when the steady state is reached, the present approach provides a straightforward scheme for the determination of inflow and outflow boundary conditions. The method is validated by comparing its predictions with exact and numerical solutions of the steady Navier-Stokes equations.

Rudisill, Edgar N., Jr.↗

An Accuracy Evaluation of Unstructured Node-Centred Finite Volume Methods

Node-centred edge-based finite volume approximations are very common in computational fluid dynamics since they are assumed to run on structured, unstructured and even on mixed grids. We analyse the accuracy properties of both first and second derivative approximations and conclude that these schemes can not be used on arbitrary grids as is often assumed. For the Euler equations first-order accuracy can be obtained if care is taken when constructing the grid. For the Navier-Stokes equations, the grid restrictions are so severe that these finite volume schemes have little advantage over structured finite difference schemes. Our theoretical results are verified through extensive computations.

Svard, Magnus↗

Static and dynamic characteristics of parallel-grooved seals

Presented is an analytical method to determine static and dynamic characteristics of annular parallel-grooved seals. The governing equations were derived by using the turbulent lubrication theory based on the law of fluid friction. Linear zero- and first-order perturbation equations of the governing equations were developed, and these equations were analytically investigated to obtain the reaction force of the seals. An analysis is presented that calculates the leakage flow rate, the torque loss, and the rotordynamic coefficients for parallel-grooved seals. To demonstrate this analysis, we show the effect of changing number of stages, land and groove width, and inlet swirl on stability of the boiler feed water pump seals. Generally, as the number of stages increased or the grooves became wider, the leakage flow rate and rotor-dynamic coefficients decreased and the torque loss increased.

Iwatsubo, Takuzo↗

Wave equations and computational models for sonic boom propagation through a turbulent atmosphere

The improved simulation of sonic boom propagation through the real atmosphere requires greater understanding of how the transient acoustic pulses popularly termed sonic booms are affected by atmospheric turbulence. A nonlinear partial differential equation that can be used to simulate the effects of smaller-scale atmospheric turbulence on sonic boom waveforms is described. The equation is first order in the time derivative and involves an extension of geometrical acoustics to include diffraction phenomena. Various terms in the equation are explained in physical terms. Such terms include those representing convection at the wave speed, diffraction, molecular relaxation, classical dissipation, and nonlinear steepening. The atmospheric turbulence enters through an effective sound speed, which varies with all three spatial coordinates, and which is the sum of the local sound speed and the component of the turbulent flow velocity projected along a central ray that connects the aircraft trajectory with the listener.

Pierce, Allan D.↗

On the integrability cases of the equation of motion for a satellite in an axially symmetric gravitational field.

Investigation of two cases of integrability of a second-order differential equation describing the projection of an axisymmetric satellite orbit on to a plane perpendicular to the rotation axis. It is demonstrated that for these two cases the integration can be carried out either by quadratures or reduced to a first-order differential equation. Analytical and physical properties are expressed, and it is shown that the equation can be derived from the classical plane eikonal equation of geometric optics.

Ghaffari, A.↗

Rotordynamic coefficients for compressible flow in tapered annular seals

Derivation of the governing equations for compressible flow in a tapered annular seal is based on Hirs' turbulent bulk-flow model. Zeroth and first-order perturbation equations are developed by an expansion in the eccentricity ratio. These equations are numerically integrated to obtain the leakage, and the direct and cross-coupled stiffness and damping coefficients. Seal parameters similar to the Space Shuttle Main Engine High Pressure Oxidizer Turbopump are used to demonstrate output from the analysis procedure. The effects of preswirl and seal taper are shown for three different length-to-diameter ratios. Generally the results indicate that prerotating the fluid significantly increases the cross-coupled stiffness but has little effect on the other coefficients, and increasing the convergent taper increases the direct stiffness while decreasing the direct damping and cross-coupled stiffness.

Nelson, C. C.↗