SEARCH · Search NASA

Results for “Ordinary Differential 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 181 records · Page 10

All sources NASA NTRS DOE OSTI

NASA NTRS1983

Finite element or Galerkin type semidiscrete schemes

A finite element of Galerkin type semidiscrete method is proposed for numerical solution of a linear hyperbolic partial differential equation. The question of stability is reduced to the stability of a system of ordinary differential equations for which Dahlquist theory applied. Results of separating the part of numerical solution which causes the spurious oscillation near shock-like response of semidiscrete scheme to a step function initial condition are presented. In general all methods produce such oscillatory overshoots on either side of shocks. This overshoot pathology, which displays a behavior similar to Gibb's phenomena of Fourier series, is explained on the basis of dispersion of separated Fourier components which relies on linearized theory to be satisfactory. Expository results represented.

Durgun, K.↗

NASA NTRS1987

Green's function solution to heat transfer of a transparent gas through a tube

A heat transfer analysis of a transparent gas flowing through a circular tube of finite thickness is presented. This study includes the effects of wall conduction, internal radiative exchange, and convective heat transfer. The natural mathematical formulation produces a nonlinear, integrodifferential equation governing the wall temperature and an ordinary differential equation describing the gas temperature. This investigation proposes to convert the original system of equations into an equivalent system of integral equations. The Green's function method permits the conversion of an integrodifferential equation into a pure integral equation. The proposed integral formulation and subsequent computational procedure are shown to be stable and accurate.

Frankel, J. I.↗

NASA NTRS1994

Hydrodynamics of rotating stars and close binary interactions: Compressible ellipsoid models

We develop a new formalism to study the dynamics of fluid polytropes in three dimensions. The stars are modeled as compressible ellipsoids, and the hydrodynamic equations are reduced to a set of ordinary differential equations for the evolution of the principal axes and other global quantities. Both viscous dissipation and the gravitational radiation reaction are incorporated. We establish the validity of our approximations and demonstrate the simplicity and power of the method by rederiving a number of known results concerning the stability and dynamical oscillations of rapidly rotating polytropes. In particular, we present a generalization to compressible fluids of Chandrasekhar's classical results for the secular and dynamical instabilities of incompressible Maclaurin spheroids. We also present several applications of our method to astrophysical problems of great current interest, such as the tidal disruption of a star by a massive black hole, the coalescence of compact binaries driven by the emission of gravitational waves, and the development of instabilities in close binary systems.

Lai, Dong↗

NASA NTRS2003

Computation of Large-Scale Structure Jet Noise Sources With Weak Nonlinear Effects Using Linear Euler

An approximate technique is presented for the prediction of the large-scale turbulent structure sound source in a supersonic jet. A linearized Euler equations code is used to solve for the flow disturbances within and near a jet with a given mean flow. Assuming a normal mode composition for the wave-like disturbances, the linear radial profiles are used in an integration of the Navier-Stokes equations. This results in a set of ordinary differential equations representing the weakly nonlinear self-interactions of the modes along with their interaction with the mean flow. Solutions are then used to correct the amplitude of the disturbances that represent the source of large-scale turbulent structure sound in the jet.

Dahl, Milo D.↗

NASA NTRSDate not supplied

An ODE-Fitting Approach to Estimate Critical Aircraft Performance Parameters for Trajectory Prediction

Ground-based decision support tools (DST) in air traffic management (ATM) typically perform trajectory prediction based on aircraft performance model (APM) parameters, but some or all of these parameters might not be readily available. In particular, the three critical parameters required for trajectory prediction are the thrust setting, drag coefficients, and takeoff weight of the aircraft. Unfortunately, these parameters are coupled and appear together in physics-based kinetic models. Past approaches utilize data from a specific phase of flight (climb, level flight or descent), where one or more of the parameters are assumed to be known and estimate the remaining unknown parameters. This approach introduces bias/errors and also does not extend to scenarios where all of the above parameters are not known with sufficient accuracy. This paper is the first of its kind to propose a generalized framework for simultaneous estimation of all three critical APM parameters (thrust, drag, and mass). The proposed approach utilizes data from both the climb and descent phases and fits the ordinary differential equation for altitude in each phase using historical trajectory data available from radar tracks or ADS-B. The approach yields a set of optimized APM parameters that are best suited to fit each historical flight record. The methodology is applied on on sample flights from three different aircraft types, and the results demonstrate low fit error and consequently will yield a high level of prediction accuracy.

trajectory prediction↗

NASA NTRS1976

Nonlinear Equations of Motion for Cantilever Rotor Blades in Hover with Pitch Link Flexibility, Twist, Precone, Droop, Sweep, Torque Offset, and Blade Root Offset

Nonlinear equations of motion for a cantilever rotor blade are derived for the hovering flight condition. The blade is assumed to have twist, precone, droop, sweep, torque offset and blade root offset, and the elastic axis and the axes of center of mass, tension, and aerodynamic center coincident at the quarter chord. The blade is cantilevered in bending, but has a torsional root spring to simulate pitch link flexibility. Aerodynamic forces acting on the blade are derived from strip theory based on quasi-steady two-dimensional airfoil theory. The equations are hybrid, consisting of one integro-differential equation for root torsion and three integro-partial differential equations for flatwise and chordwise bending and elastic torsion. The equations are specialized for a uniform blade and reduced to nonlinear ordinary differential equations by Galerkin's method. They are linearized for small perturbation motions about the equilibrium operating condition. Modal analysis leads to formulation of a standard eigenvalue problem where the elements of the stability matrix depend on the solution of the equilibrium equations. Two different forms of the root torsion equation are derived that yield virtually identical numerical results. This provides a reasonable check for the accuracy of the equations.

Hodges, D. H.↗

NASA NTRS1967

Heat transfer in a plane incompressible laminar jet

Partial differential equation for heat transfer in plane incompressible laminar jet is reduced by similarity transformation to ordinary differential equation

HEAT TRANSFER↗

NASA NTRS1983

Modeling and control of distributed structures

An inherent incongruity in the modeling of controls for structures is discussed. Structures are basically distributed-parameter systems, described by partial differential equations, and control theory is concerned almost exclusively with discrete (in space) systems, described by ordinary differential equations. The standard approach to solving this dilemma is to discretize the system in space, which precludes the use of distributed controls. A different approach, known as the independent modal-space control method, is designed to eliminate the incongruity by bringing about a closer correspondence between modeling and control theory. The independent modal-space control method can treat distributed structures as well as discretized models and it permits design of both distributed and discrete-point controls.

Meirovitch, L.↗

NASA NTRS1986

A three-dimensional axisymmetric calculation procedure for turbulent flows in a radial vaneless diffuser

An analytical model is proposed to calculate the three-dimensional axisymmetric turbulent flowfield in a radial vneless diffuser. The model assumes that the radial and tangential boundary layer profiles are approximated by power law profiles. Then, using the integrated radial and tangential momentum and continuity equations for the boundary layer and corresponding inviscid equations for the core flow, there results six ordinary differential equation in six unknowns which are easily solved using a Runge-Kutta technique. A model is also proposed for fully developed flow. The results using this technique were compared with the results from a three-dimensional viscous, axisymmetric duct code and with experimental data and good quantitative agreement was obtained.

Schumann, L. F.↗

NASA NTRS1987

Time-dependent approximation schemes for some problems of parameter estimation in distributed systems

A parameter estimation method that can be used to estimate functional parameters in delay differential equations and moving boundary problems is discussed. In either problem, the original model equation (which is infinite-dimensional) is approximated by a system of ordinary differential equations that can be solved numerically in an efficient way. The approximation scheme is based on time-dependent spline elements. For the delay equation with time-varying delay, convergence results are presented that indicate the estimates obtained using the approximating system. Numerical test examples converge in some sense to a best-fit parameter for the original system are included by means of which time-varying and state-dependent delays and a time-varying diffusion coefficient in a one-phase, one-dimensional Stefan problem are estimated.

Murphy, K. A.↗

NASA NTRS1990

Mathematical modeling of the flow field and particle motion in a rotating bioreactor at unit gravity and microgravity

The biotechnology group at NASA Johnson Space Center is developing systems for culturing mammalian cells that stimulate some aspect of microgravity and provide a low shear environment for microgravity-based studies on suspension and anchorage dependent cells. The design of these vessels for culturing cells is based on the need to suspend cells and aggregates of cells and microcarrier beads continually in the culturing medium. The design must also provide sufficient circulation for adequate mass transfer of nutrients to the cells and minimize the total force on the cells. Forces, resulting from sources such as hydrodynamic fluid shear and collisions of cells and walls of the vessels, may damage delicate cells and degrade the formation of three dimensional structures. This study examines one particular design in both unit gravity and microgravity based on two concentric cylinders rotating in the same direction at different speeds to create a Couette flow between them. A numerical simulation for the flow field and the trajectories of particles in the vessel. The flow field for the circulation of the culturing medium is modeled by the Navier-Stokes equations. The forces on a particle are assumed to be drag from the fluid's circulation, buoyancy from the gravitational force and centrifugal force from the rotation of the vessel. The problem requires first solving the system of partial differential equations for the fluid flow by a finite difference method and then solving the system of ordinary differential equations for the trajectories by Gear's stiff method. Results of the study indicate that the trajectories in unit gravity and microgravity are very similar except for small spatial deviations on the fast time scale in unit gravity. The total force per unit cross sectional area on a particle in microgravity, however, is significantly smaller than the corresponding value in unit gravity, which is also smaller than anticipated. Hence, this study indicates that this design for a bioreactor with optimal rates of rotation can provide a good environment for culturing cells in microgravity with adequate circulation and minimal force on the cells.

Boyd, Ernest J.↗

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

Burnett solutions along the stagnation line of a cooled cylinder in low-density hypersonic flows

The Burnett equations in the cylindrical coordinate system have been applied to the hypersonic flow over a cylinder in low density environments. Through a local similarity transformation, the Burnett equations are degenerated into nonlinear third-order ordinary differential equations which are only valid along the stagnation line. These equations are solved by the multi-staged Runge-Kutta integration technique. Steady state solutions are asymptotically reached by the time-marching procedure. A newly-developed computer code is used to solve both the Navier-Stokes and the Burnett equations within the shock layer. The computed heat transfer coefficients are consistent with wind-tunnel and STS flight data. Comparisons with the DSMC solutions and the Boltzmann solutions show that the Burnett equations are more appropriate as the governing equations than the Navier-Stokes equations in the transitional flow regime.

Liaw, G. S.↗

NASA NTRS1995

Hydrodynamics of coalescing binary neutron stars: Ellipsoidal treatment

We employ an approximate treatment of dissipative hydrodynamics in three dimensions to study the coalescence of binary neutron stars driven by the emission of gravitational waves. The stars are modeled as compressible ellipsoids obeying a polytropic equation of state; all internal fluid velocities are assumed to be linear functions of the coordinates. The hydrodynamics equations then reduce to a set of coupled ordinary differential equations for the evolution of the principal axes of the ellipsoids, the internal velocity parameters, and the binary orbital parameters. Gravitational radiation reaction and viscous dissipation are both incorporated. We set up exact initial binary equilibrium configurations and follow the transition from the quasi-static, secular decay of the orbit at large separation to the rapid dynamical evolution of the configurations just prior to contact. A hydrodynamical instability resulting from tidal interactions significantly accelerates the coalescence at small separation, leading to appreciable radial infall velocity and tidal lag angles near contact. This behavior is reflected in the gravitational waveforms and may be observable by gravitational wave detectors under construction. In cases where the neutron stars have spins which are not aligned with the orbital angular momentum, the spin-induced quadrupole moment can lead to precession of the orbital plane and therefore modulation of the gravitational wave amplitude even at large orbital radius. However, the amplitude of the modulation is small for typical neutron star binaries with spins much smaller than the orbital angular momentum.

Lai, Dong↗

NASA NTRS1995

Nonlinear evolution of magnetic flux ropes. 2: Finite beta plasma

In this second paper on the evolution of magnetic flux ropes we study the effects of gas pressure. We assume that the energy transport is described by a polytropic relationship and reduce the set of ideal MHD equations to a single, second-order, nonlinear, ordinary differential equation for the evolution function. For this conservative system we obtain a first integral of motion. To analyze the possible motions, we use a mechanical analogue -- a one-dimensional, nonlinear oscillator. We find that the effective potential for such an oscillator depends on two parameters: the polytropic index gamma and a dimensionless quantity kappa the latter being a function of the plasma beta, the strength of the azimuthal magnetic field relative to the axial field of the flux rope, and gamma. Through a study of this effective potential we classify all possible modes of evolution of the system. In the main body of the paper, we focus on magnetic flux ropes whose field and gas pressure increase steadily towards the symmetry axis. In this case, for gamma greater than 1 and all values of kappa, only oscillations are possible. For gamma less than 1, however, both oscillations and expansion are allowed. For gamma less than 1 and kappa below a critical value, the energy of the nonlinear oscillator determines whether the flux rope will oscillate or expand to infinity. For gamma less than 1 and kappa above critical, however, only expansion occurs. Thus by increasing kappa while keeping gamma fixed (less than 1), a phase transition occurs at kappa = kappa(sub critical) and the oscillatory mode disappears. We illustrate the above theoretical considerations by the example of a flux rope of constant field line twist evolving self-similarly. For this example, we present the full numerical MHD solution. In an appendix to the paper we catalogue all possible evolutions when (1) either the magnetic field or (2) the gas pressure decreases monotonically toward the axis. We find that in these cases critical conditions can occur for gamma greater than 1. While in most cases the flux rope collapses, there are notable exceptions when, for certain ranges of kappa and gamma, collapse may be averted.

Osherovich, V. A.↗

NASA NTRS2008

Simulink Model of the Ares I Upper Stage Main Propulsion System

A numerical model of the Ares I upper stage main propulsion system is formulated based on first principles. Equation's are written as non-linear ordinary differential equations. The GASP fortran code is used to compute thermophysical properties of the working fluids. Complicated algebraic constraints are numerically solved. The model is implemented in Simulink and provides a rudimentary simulation of the time history of important pressures and temperatures during re-pressurization, boost and upper stage firing. The model is validated against an existing reliable code, and typical results are shown.

Burchett, Bradley T.↗

NASA NTRS2009

Simulating Bone Loss in Microgravity Using Mathematical Formulations of Bone Remodeling

Most mathematical models of bone remodeling are used to simulate a specific bone disease, by disrupting the steady state or balance in the normal remodeling process, and to simulate a therapeutic strategy. In this work, the ability of a mathematical model of bone remodeling to simulate bone loss as a function of time under the conditions of microgravity is investigated. The model is formed by combining a previously developed set of biochemical, cellular dynamics, and mechanical stimulus equations in the literature with two newly proposed equations; one governing the rate of change of the area of cortical bone tissue in a cross section of a cylindrical section of bone and one governing the rate of change of calcium in the bone fluid. The mechanical stimulus comes from a simple model of stress due to a compressive force on a cylindrical section of bone which can be reduced to zero to mimic the effects of skeletal unloading in microgravity. The complete set of equations formed is a system of first order ordinary differential equations. The results of selected simulations are displayed and discussed. Limitations and deficiencies of the model are also discussed as well as suggestions for further research.

Pennline, James A.↗