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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 163 records · Page 9

All sources NASA NTRS DOE OSTI

NASA NTRS1985

Formulation and application of Russell's method

It is shown that the numerical technique of Russell's momentum approach can be derived by using Hamilton's principle and Vance's numerical scheme. It results in a set of first order differnce equations for solving the angular velocities. The numerical examples show that the method is reliable. The algorithm is modified next to perform the analysis of N-body systems with closed loop topology. To increase the formulation flexibility, the equations of motion are represented by using Cartesian coordinates and Lagrange multipliers. The algorithm consists of two parts, Vance's scheme and an unconstrained minimization. The Vance's scheme is used to find the angular velocities, and the unconstrained minimization is applied to provide the correct angular displacements. The proposed scheme is further extended to find the design sensitivity of an N-body system with closed loop configuration, and to carry out the design optimization as well. The numerical example of a small-scaled mechanical system is presented to verify the proposed formulation.

Hou, J. W.↗

NASA NTRS1986

An embedding method for the steady Euler equations

Certain difficulties arise in connection with the numerical solution of a direct finite difference representation of the steady Euler equations. Johnson (1979, 1981, 1982) has, therefore, proposed a surrogate-equation technique, in which the first-order steady Euler equations are embedded in a certain second-order system of equations. The present paper is concerned with the theoretical justification for such an embedding approach. For the numerical solution of the two-dimensional steady Euler equations, it is shown that, under a continuity restriction, it is possible to solve a second-order embedded system together with appropriate additional boundary conditions. The result indicates that a more direct and potentially more efficient approach to the steady solutions exists than the alternative of solving the unsteady equations.

Chang, S.-H.↗

NASA NTRS1986

The alpha dynamo parameter and measurability of helicities in magnetohydrodynamic turbulence

Alpha, an important parameter in dynamo theory, is shown to be proportional to either the kinetic, current, magnetic, or velocity helicities of the fluctuating magnetic field and fluctuating velocity field. The particular helicity to which alpha is proportional depends on the assumptions used in deriving the first-order smoothed equations that describe the alpha effect. In two cases, viz., when alpha is proportional to either the magnetic helicity or velocity helicity, alpha can be determined experimentally from two-point measurements of the fluctuating fields in incompressible, homogeneous turbulence with arbitrary rotational symmetry. For the other two possibilities, alpha can be determined if the turbulence is isotropic.

Matthaeus, W. H.↗

NASA NTRS1986

A spline-based parameter estimation technique for static models of elastic structures

The problem of identifying the spatially varying coefficient of elasticity using an observed solution to the forward problem is considered. Under appropriate conditions this problem can be treated as a first order hyperbolic equation in the unknown coefficient. Some continuous dependence results are developed for this problem and a spline-based technique is proposed for approximating the unknown coefficient, based on these results. The convergence of the numerical scheme is established and error estimates obtained.

Dutt, P.↗

NASA NTRS1987

Optimal cure cycle design of a resin-fiber composite laminate

A unified computed aided design method was studied for the cure cycle design that incorporates an optimal design technique with the analytical model of a composite cure process. The preliminary results of using this proposed method for optimal cure cycle design are reported and discussed. The cure process of interest is the compression molding of a polyester which is described by a diffusion reaction system. The finite element method is employed to convert the initial boundary value problem into a set of first order differential equations which are solved simultaneously by the DE program. The equations for thermal design sensitivities are derived by using the direct differentiation method and are solved by the DE program. A recursive quadratic programming algorithm with an active set strategy called a linearization method is used to optimally design the cure cycle, subjected to the given design performance requirements. The difficulty of casting the cure cycle design process into a proper mathematical form is recognized. Various optimal design problems are formulated to address theses aspects. The optimal solutions of these formulations are compared and discussed.

Hou, Jean W.↗

NASA NTRS1985

UNAERO: A package of FORTRAN subroutines for approximating unsteady aerodynamics in the time domain

This report serves as an instruction and maintenance manual for a collection of CDC CYBER FORTRAN IV subroutines for approximating the unsteady aerodynamic forces in the time domain. The result is a set of constant-coefficient first-order differential equations that approximate the dynamics of the vehicle. Provisions are included for adjusting the number of modes used for calculating the approximations so that an accurate approximation is generated. The number of data points at different values of reduced frequency can also be varied to adjust the accuracy of the approximation over the reduced-frequency range. The denominator coefficients of the approximation may be calculated by means of a gradient method or a least-squares approximation technique. Both the approximation methods use weights on the residual error. A new set of system equations, at a different dynamic pressure, can be generated without the approximations being recalculated.

Dunn, H. J.↗

NASA NTRS1988

Development of an unsteady wake theory appropriate for aeroelastic analyses of rotors in hover and forward flight

The purpose of this research is the development of an unsteady aerodynamic model for rotors such that it can be used in conventional aeroelastic analysis (e.g., eigenvalue determination and control system design). For this to happen, the model must be in a state-space formulation such that the states of the flow can be defined, calculated and identified as part of the analysis. The fluid mechanics of the problem is given by a closed-form inversion of an acceleration potential. The result is a set of first-order differential equations in time for the unknown flow coefficients. These equations are hierarchical in the sense that they may be truncated at any number of radial or azimuthal terms.

Peters, David A.↗

NASA NTRS1988

Consistent boundary conditions for reduced Navier-Stokes (RNS) scheme applied to three-dimensional internal viscous flows

Three-dimensional internal viscous flow problems are presently addressed by a consistent and efficient set of boundary conditions for the multisweep space-marching, pressure-elliptic Reduced Navier-Stokes scheme. Since continuity is by this means directly satisfied at all points in the flow domain, the first-order momentum equations are directly solvable for pressure without the requirement for a Poisson pressure-correction equation. Incompressible flow solutions are obtained for straight and curved ducts of square cross section, in order to validate the procedure. Usefulness is demonstrated for internal flows with strong interactions, as would be found in turbomachine geometries.

Reddy, D. R.↗

NASA NTRS1989

A spline-based parameter estimation technique for static models of elastic structures

Dutt, P.↗

NASA NTRS1990

Cell averaging Chebyshev methods for hyperbolic problems

A cell averaging method for the Chebyshev approximations of first order hyperbolic equations in conservation form is described. Formulas are presented for transforming between pointwise data at the collocation points and cell averaged quantities, and vice-versa. This step, trivial for the finite difference and Fourier methods, is nontrivial for the global polynomials used in spectral methods. The cell averaging methods presented are proven stable for linear scalar hyperbolic equations and present numerical simulations of shock-density wave interaction using the new cell averaging Chebyshev methods.

Wei, Cai↗

NASA NTRS1990

Theory of biaxial graded-index optical fiber

The problem of wave propagation in a biaxial graded-index fiber with circular symmetry is considered. The problem is formulated in terms of four first-order differential equations for the tangential components of the electric and magnetic fields. A general solution method for solving systems of differential equations is presented. This solution method is then used to solve the system of equations for a particular example of a biaxial graded-index fiber. Numerical results for the propagation constant in the fiber are also given.

Kawalko, Stephen F.↗

NASA NTRS1989

Interpolation solutions for the problem of synthesis of dual-shaped offset reflector antennas

Synthesis of dual-shaped offset reflector antennas to control the exit aperture distribution of amplitude and phase has received considerable attention in recent years. For a given feed illumination and desired aperture field distribution, an exact formulation of the problem of simultaneously synthesizing the shapes of the sub and main reflectors was presented recently by Galindo-Israel et al. (1987) in terms of a set of nonlinear first-order differential equations. In this paper, a numerical approach to solving these equations is discussed which circumvents some of the difficulties encountered by Galindo-Israel et al., particularly for small values of theta.

Jervase, Joseph A.↗

NASA NTRS1990

Analysis of the single-blow transient testing technique for perforated plate heat exchangers

An analysis of the single-blow transient testing technique for heat exchanger surfaces is made for perforated plate and similar discontinuous surfaces. The model assumes that there is no temperature variation across each plate, and allows for axial conduction in spacers: the resulting axial fluid temperature profile is discontinuous, rather than the usual continuous profile. Numerical solutions are obtained to the resulting set of coupled first-order differential equations for a step change in inlet fluid temperature. Results are presented in tabular form, which allow the heat transfer coefficients to be calculated from test data using the maximum slope technique.

Rodriguez, Jose I.↗

NASA NTRS1991

A method for the analysis of biaxial graded-index optical fibers

The problem of wave propagation is a biaxial graded-index fiber with circular symmetry is considered. The problem is formulated in terms of four first-order differential equations for the tangential components of the electric and magnetic fields. A general solution method for solving systems of differential equations is presented. This solution method is then used to solve the system of equations for a particular example of a biaxial graded-index fiber. Numerical results for the propagation constant in the fiber are given.

Kawalko, Stephen F.↗

NASA NTRS1992

Multi-wave interaction theory for wind-generated surface gravity waves

Consistently employing the assumption of localness of wave-wave interactions in the wavenumber space, the Kolmogorov treatment of the energy cascade is applied to the case of wind-generated surface gravity waves. The effective number v of resonantly interacting wave harmonics is not limited to four but is found as a solution of a coupled system of equations expressing: (i) the dependence of the spectrum shape on the degree of the wave nonlinearity, and (ii) the continuity of the wave action flux through the spectrum given a continuous positive input from wind. The latter is specified in a Miles-type fashion, and a simple scaling relationship based on the concept of the turnover time is derived in place of the kinetic equation. The mathematical problem is reduced to an ordinary differential equation of first order. The exponent in the 'power law' for the spectral density of the wave potential energy and the effective number of resonantly interacting wave harmonics are found as functions of the wave frequency and of external factors of wind-wave interaction. The solution is close to the Zakharov-Filonenko spectrum at low frequencies and low wind input while approaching the Phillips spectrum at high frequencies and sufficiently high wind.

Glazman, Roman E.↗

NASA NTRS1993

On the superconvergence of Galerkin methods for hyperbolic IBVP

Finite element Galerkin methods for periodic first order hyperbolic equations exhibit superconvergence on uniform grids at the nodes, i.e., there is an error estimate 0(h(sup 2r)) instead of the expected approximation order 0(h(sup r)). It will be shown that no matter how the approximating subspace S(sup h) is chosen, the superconvergence property is lost if there are characteristics leaving the domain. The implications of this result when constructing compact implicit difference schemes is also discussed.

Gottlieb, David↗

NASA NTRS1997

Multi-Body Analysis of a Tiltrotor Configuration

The paper describes the aeroelastic analysis of a tiltrotor configuration. The 1/5 scale wind tunnel semispan model of the V-22 tiltrotor aircraft is considered. The analysis is performed by means of a multi-body code, based on an original formulation. The differential equilibrium problem is stated in terms of first order differential equations. The equilibrium equations of every rigid body are written, together with the definitions of the momenta. The bodies are connected by kinematic constraints, applied in form of Lagrangian multipliers. Deformable components are mainly modelled by means of beam elements, based on an original finite volume formulation. Multi-disciplinar problems can be solved by adding user-defined differential equations. In the presented analysis the equations related to the control of the swash-plate of the model are considered. Advantages of a multi-body aeroelastic code over existing comprehensive rotorcraft codes include the exact modelling of the kinematics of the hub, the detailed modelling of the flexibility of critical hub components, and the possibility to simulate steady flight conditions as well as wind-up and maneuvers. The simulations described in the paper include: 1) the analysis of the aeroelastic stability, with particular regard to the proprotor/pylon instability that is peculiar to tiltrotors, 2) the determination of the dynamic behavior of the system and of the loads due to typical maneuvers, with particular regard to the conversion from helicopter to airplane mode, and 3) the stress evaluation in critical components, such as the pitch links and the conversion downstop spring.

Ghiringhelli, G. L.↗

NASA NTRS2006

Analytical and Numerical Results for an Adhesively Bonded Joint Subjected to Pure Bending

A one-dimensional, semi-analytical methodology that was previously developed for evaluating adhesively bonded joints composed of anisotropic adherends and adhesives that exhibit inelastic material behavior is further verified in the present paper. A summary of the first-order differential equations and applied joint loading used to determine the adhesive response from the methodology are also presented. The method was previously verified against a variety of single-lap joint configurations from the literature that subjected the joints to cases of axial tension and pure bending. Using the same joint configuration and applied bending load presented in a study by Yang, the finite element analysis software ABAQUS was used to further verify the semi-analytical method. Linear static ABAQUS results are presented for two models, one with a coarse and one with a fine element meshing, that were used to verify convergence of the finite element analyses. Close agreement between the finite element results and the semi-analytical methodology were determined for both the shear and normal stress responses of the adhesive bondline. Thus, the semi-analytical methodology was successfully verified using the ABAQUS finite element software and a single-lap joint configuration subjected to pure bending.

Smeltzer, Stanley S., III↗