Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “differential equations”

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 955 records · Page 53
All sourcesNASA NTRSDOE OSTI
NASA NTRS2004

Variable Step Integration Coupled with the Method of Characteristics Solution for Water-Hammer Analysis, A Case Study

One-dimensional water-hammer modeling involves the solution of two coupled non-linear hyperbolic partial differential equations (PDEs). These equations result from applying the principles of conservation of mass and momentum to flow through a pipe, and usually the assumption that the speed at which pressure waves propagate through the pipe is constant. In order to solve these equations for the interested quantities (i.e. pressures and flow rates), they must first be converted to a system of ordinary differential equations (ODEs) by either approximating the spatial derivative terms with numerical techniques or using the Method of Characteristics (MOC). The MOC approach is ideal in that no numerical approximation errors are introduced in converting the original system of PDEs into an equivalent system of ODEs. Unfortunately this resulting system of ODEs is bound by a time step constraint so that when integrating the equations the solution can only be obtained at fixed time intervals. If the fluid system to be modeled also contains dynamic components (i.e. components that are best modeled by a system of ODEs), it may be necessary to take extremely small time steps during certain points of the model simulation in order to achieve stability and/or accuracy in the solution. Coupled together, the fixed time step constraint invoked by the MOC, and the occasional need for extremely small time steps in order to obtain stability and/or accuracy, can greatly increase simulation run times. As one solution to this problem, a method for combining variable step integration (VSI) algorithms with the MOC was developed for modeling water-hammer in systems with highly dynamic components. A case study is presented in which reverse flow through a dual-flapper check valve introduces a water-hammer event. The predicted pressure responses upstream of the check-valve are compared with test data.

Turpin, Jason B.↗
NASA NTRS1994

The Dynamics of Some Iterative Implicit Schemes

The global asymptotic nonlinear behavior of some standard iterative procedures in solving nonlinear systems of algebraic equations arising from four implicit linear multistep methods (LMMs) in discretizing 2 x 2 systems of first-order autonomous nonlinear ordinary differential equations is analyzed using the theory of dynamical systems. With the aid of parallel Connection Machines (CM-2 and CM-5), the associated bifurcation diagrams as a function of the time step, and the complex behavior of the associated 'numerical basins of attraction' of these iterative implicit schemes are revealed and compared. Studies showed that all of the four implicit LMMs exhibit a drastic distortion and segmentation but less shrinkage of the basin of attraction of the true solution than standard explicit methods. The numerical basins of attraction of a noniterative implicit procedure mimic more closely the basins of attraction of the differential equations than the iterative implicit procedures for the four implicit LMMs.

Yee, H. C.↗
NASA NTRS1981

An analytic solution to the classical two-body problem with drag

The two body problem was analyzed with a specific drag model. The model treats drag as a force proportional to the vector velocity and inversely proportional to the distance to the center of attraction. The solution is expressed in terms of known functions and is of a simple and compact form. The time of flight is expressed as a quadrature in the true anomaly. The results are: (1) development of a vector differential equation which allows analysis of an infinite number of gravitational and drag models; and (2) obtaining the solution of a linear differential equation using the inverse method of laplace transforms.

Mittleman, D.↗
NASA NTRS1984

Singular asymptotic expansions in nonlinear rotordynamics

During hot firing ground testing of the Space Shuttle's Main Engine, vibrations of the liquid oxygen pump occur at frequencies which cannot be explained by the linear Jeffcott model of the rotor. The model becomes nonlinear after accounting for deadband, side forces, and rubbing. Two phenomena present in the numerical solutions of the differential equations are unexpected periodic orbits of the rotor and tracking of the nonlinear frequency. A multiple scale asymptotic expansion of the differential equations is used to give an analytic explanation of these characteristics.

Day, W. B.↗
NASA NTRS1985

Singular asymptotic expansions in nonlinear rotordynamics

During hot firing ground testing of the Space shuttle's Main Engine, vibrations of the liquid oxygen pump occur at frequencies which cannot be explained by the linear Jeffcott model of the rotor. The model becomes nonlinear after accounting for deadband, side forces, and rubbing. Two phenomena present in the numerical solutions of the differential equations are unexpected periodic orbits of the rotor and tracking of the nonlinear frequency. A multiple scale asymptotic expansion of the differential equations is used to give an analytic explanation of these characteristics.

Day, W. B.↗
NASA NTRS2000

LSENS, The NASA Lewis Kinetics and Sensitivity Analysis Code

A general chemical kinetics and sensitivity analysis code for complex, homogeneous, gas-phase reactions is described. The main features of the code, LSENS (the NASA Lewis kinetics and sensitivity analysis code), are its flexibility, efficiency and convenience in treating many different chemical reaction models. The models include: static system; steady, one-dimensional, inviscid flow; incident-shock initiated reaction in a shock tube; and a perfectly stirred reactor. In addition, equilibrium computations can be performed for several assigned states. An implicit numerical integration method (LSODE, the Livermore Solver for Ordinary Differential Equations), which works efficiently for the extremes of very fast and very slow reactions, is used to solve the "stiff" ordinary differential equation systems that arise in chemical kinetics. For static reactions, the code uses the decoupled direct method to calculate sensitivity coefficients of the dependent variables and their temporal derivatives with respect to the initial values of dependent variables and/or the rate coefficient parameters. Solution methods for the equilibrium and post-shock conditions and for perfectly stirred reactor problems are either adapted from or based on the procedures built into the NASA code CEA (Chemical Equilibrium and Applications).

Radhakrishnan, K.↗
NASA NTRS1982

Hybrid state vector methods for structural dynamic and aeroelastic boundary value problems

A computational technique is developed that is suitable for performing preliminary design aeroelastic and structural dynamic analyses of large aspect ratio lifting surfaces. The method proves to be quite general and can be adapted to solving various two point boundary value problems. The solution method, which is applicable to both fixed and rotating wing configurations, is based upon a formulation of the structural equilibrium equations in terms of a hybrid state vector containing generalized force and displacement variables. A mixed variational formulation is presented that conveniently yields a useful form for these state vector differential equations. Solutions to these equations are obtained by employing an integrating matrix method. The application of an integrating matrix provides a discretization of the differential equations that only requires solutions of standard linear matrix systems. It is demonstrated that matrix partitioning can be used to reduce the order of the required solutions. Results are presented for several example problems in structural dynamics and aeroelasticity to verify the technique and to demonstrate its use. These problems examine various types of loading and boundary conditions and include aeroelastic analyses of lifting surfaces constructed from anisotropic composite materials.

Lehman, L. L.↗
DOE OSTI2025

ORMATEX

The Oak Ridge Matrix Exponential (ORMATEX) software library contains methods to compute the matrix exponential and the action of the matrix exponential on a vector. Additionally, this package contains the related methods for the phi-functions which commonly appear in a wide class of exponential time integration methods. Krylov methods are provided to evaluate the matrix exponential-vector and phi-vector products for cases where the matrix is large and sparse. Utilizing these methods, ORMATEX implements performant exponential integrators for large systems of coupled ordinary differential equations (ODEs). The exponential time integration routines in ORMATEX are particularly suitable to large, stiff systems of equations. These routines may be utilized as a competitive alternative to classical implicit and explicit time integration schemes for certain classes of differential equations where the problem stiffness can be predominately explained by the linear terms.

Gurecky, William [Oak Ridge National Laboratory (O↗
NASA NTRS1977

A user's guide for V174, a program using a finite difference method to analyze transonic flow over oscillating wings

The design and usage of a pilot program using a finite difference method for calculating the pressure distributions over harmonically oscillating wings in transonic flow are discussed. The procedure used is based on separating the velocity potential into steady and unsteady parts and linearizing the resulting unsteady differential equation for small disturbances. The steady velocity potential which must be obtained from some other program, is required for input. The unsteady differential equation is linear, complex in form with spatially varying coefficients. Because sinusoidal motion is assumed, time is not a variable. The numerical solution is obtained through a finite difference formulation and a line relaxation solution method.

Butler, T. D.↗
NASA NTRS1965

Necessary conditions for a multistage bolza- mayer problem involving control variables and having inequality and finite equation constraints

A multiplier rule and analogues of the Weierstrass and Clebsch conditions are developed for a multistage Bolza-Meyer calculus of variations problems. The number of stages is fixed, but partition points defining state boundaries are variable. Discontinuities are allowed in variables finite equations and inequalities, as well as differential equations, all of which involve control variables. An appendix summarizes some of the results obtained by C. H. Denbow, as modified by R. W. hunt, for a generalized Bolza problem. The appendix is independent of the rest of the paper.

Differential equation↗
NASA NTRS1988

Coupled bending-torsion steady-state response of pretwisted, nonuniform rotating beams using a transfer-matrix method

Using the Newtonian method, the equations of motion are developed for the coupled bending-torsion steady-state response of beams rotating at constant angular velocity in a fixed plane. The resulting equations are valid to first order strain-displacement relationships for a long beam with all other nonlinear terms retained. In addition, the equations are valid for beams with the mass centroidal axis offset (eccentric) from the elastic axis, nonuniform mass and section properties, and variable twist. The solution of these coupled, nonlinear, nonhomogeneous, differential equations is obtained by modifying a Hunter linear second-order transfer-matrix solution procedure to solve the nonlinear differential equations and programming the solution for a desk-top personal computer. The modified transfer-matrix method was verified by comparing the solution for a rotating beam with a geometric, nonlinear, finite-element computer code solution; and for a simple rotating beam problem, the modified method demonstrated a significant advantage over the finite-element solution in accuracy, ease of solution, and actual computer processing time required to effect a solution.

Gray, Carl E., Jr.↗
NASA NTRS1972

Longitudinal dynamic stability of a shuttle vehicle.

Analytical study of the longitudinal dynamic stability of a nonrolling, lifting vehicle gliding at hypersonic speeds. The analysis applies to shuttle vehicles designed for operating up to the rim of a planetary atmosphere. A general nondimensional time transformation is introduced to derive a unified second-order linear differential equation for the angle of attack, valid for all types of reentry of a general type of vehicle. The stability of motion is discussed for two fundamental regimes of flight that are based on widely different assumptions. For near ballistic entry along a straight line trajectory, the equation reduces to a confluent hypergeometric equation, the solution of which can be expressed in terms of Whittaker's function. Using a theorem in the theory of stability of differential equations, criteria for damped oscillations are derived. It is shown that the aerodynamic criteria for stability are the same as for the case of ballistic entry. In addition, for each vehicle configuration, and specified planetary atmosphere, there exists an altitude range where the angle of attack frequency is nearly equal to the orbital frequency causing instability in pitch. This resonance instability is due to the ellipticity of the orbit. Criteria for eccentricity instability are derived.

Vinh, N. X.↗
NASA NTRS1977

Prediction of hydrodynamics and chemistry of confined turbulent methane-air flames with attention to formation of oxides of nitrogen

A formulation of the governing partial differential equations for fluid flow and reacting chemical species in a tubular combustor is presented. A numerical procedure for the solution of the governing differential equations is described, and models for chemical equilibrium and chemical kinetics calculations are presented. The chemical equilibrium model is used to characterize the hydrocarbon reactions. The chemical kinetics model is used to predict the concentrations of the oxides of nitrogen. The combustor consists of a cylindrical duct of varying cross sections with concentric streams of gaseous fuel and air entering the duct at one end. Four sample cases with specified inlet and boundary conditions are considered, and the results are discussed

Elghobashi, S.↗
NASA NTRS1978

Long range orbital error estimation for applications satellites

A method of optimum orbital averaging was employed to study the long range accuracy potential of polar orbiting applications satellites. This approach involved the determination of the boundary conditions of one set of differential equations of motion by adjusting the initial conditions in a least square sense with the use of data generated by another set of differential equations of motion.

Bonavito, N. L.↗
NASA NTRS1980

Acoustic transmission matrix of a variable area duct or nozzle carrying a compressible subsonic flow

The differential equations governing the propagation of sound in a variable area duct or nozzle carrying a one dimensional subsonic compressible fluid flow are derived and put in state variable form using acoustic pressure and particle velocity as the state variables. The duct or nozzle is divided into a number of regions. The region size is selected so that in each region the Mach number can be assumed constant and the area variation can be approximated by an exponential area variation. Consequently, the state variable equation in each region has constant coefficients. The transmission matrix for each region is obtained by solving the constant coefficient acoustic state variable differential equation. The transmission matrix for the duct or nozzle is the product of the individual transmission matrices of each region. Solutions are presented for several geometries with and without mean flow.

Miles, J. H.↗
NASA NTRS1980

Acoustic transmission matrix of a variable area duct or nozzle carrying a compressible subsonic flow

The differential equations governing the propagation of sound in a variable area duct or nozzle carrying a one-dimensional subsonic compressible fluid flow are derived and put in state variable form using acoustic pressure and particle velocity as the state variables. The duct or nozzle is divided into a number of regions. The region size is selected so that in each region the Mach number can be assumed constant and the area variation can be approximated by an exponential area variation. Consequently, the state variable equation in each region has constant coefficients. The transmission matrix for each region is obtained by solving the constant coefficient acoustic state variable differential equation. The transmission matrix for the duct or nozzle is the product of the individual transmission matrices of each region. Solutions are presented for several geometries with and without mean flow.

Miles, J. H.↗
NASA NTRS1999

Quasi-Static Viscoelastic Finite Element Model of an Aircraft Tire

An elastic large displacement thick-shell mixed finite element is modified to allow for the calculation of viscoelastic stresses. Internal strain variables are introduced at the element's stress nodes and are employed to construct a viscous material model. First order ordinary differential equations relate the internal strain variables to the corresponding elastic strains at the stress nodes. The viscous stresses are computed from the internal strain variables using viscous moduli which are a fraction of the elastic moduli. The energy dissipated by the action of the viscous stresses is included in the mixed variational functional. The nonlinear quasi-static viscous equilibrium equations are then obtained. Previously developed Taylor expansions of the nonlinear elastic equilibrium equations are modified to include the viscous terms. A predictor-corrector time marching solution algorithm is employed to solve the algebraic-differential equations. The viscous shell element is employed to computationally simulate a stair-step loading and unloading of an aircraft tire in contact with a frictionless surface.

Johnson, Arthur R.↗