Search NASA⌕ Search

SEARCH · Search NASA

Results for “INITIAL VALUE PROBLEM”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
NASA NTRS1976

On the Navier-Stokes equations with constant total temperature

For various applications in fluid dynamics, one can assume that the total temperature is constant. Therefore, the energy equations can be replaced by an algebraic relation. The resulting set of equations in the inviscid case is analyzed in this paper. It is shown that the system is strictly hyperbolic and well posed for the initial-value problem. Boundary conditions are described such that the linearized system is well posed. The hopscotch method is investigated and numerical results are presented.

Gottlieb, D.↗
NASA NTRS1976

Thermal instability in post-flare plasmas

The cooling of post-flare plasmas is discussed and the formation of loop prominences is explained as due to a thermal instability. A one-dimensional model was developed for active loop prominences. Only the motion and heat fluxes parallel to the existing magnetic fields are considered. The relevant size scales and time scales are such that single-fluid MHD equations are valid. The effects of gravity, the geometry of the field and conduction losses to the chromosphere are included. A computer code was constructed to solve the model equations. Basically, the system is treated as an initial value problem (with certain boundary conditions at the chromosphere-corona transition region), and a two-step time differencing scheme is used.

Antiochos, S. K.↗
NASA NTRS1973

Program for solution of ordinary differential equations

A program for the solution of linear and nonlinear first order ordinary differential equations is described and user instructions are included. The program contains a new integration algorithm for the solution of initial value problems which is particularly efficient for the solution of differential equations with a wide range of eigenvalues. The program in its present form handles up to ten state variables, but expansion to handle up to fifty state variables is being investigated.

Sloate, H.↗
NASA NTRS1979

Numerical spatial marching techniques in duct acoustics

Direct calculation of the internal structure of a ducted noise source from farfield pressure measurements is regarded as an initial value problem, where the pressure and pressure gradient (farfield impedance) are assumed to be known along a line in the farfield. If pressure and impedance are known at the boundary of the farfield, the pressure can be uniquely determined in the vicinity of the inlet and inside the inlet ducting. A marching procedure is developed which, with this information obtained from measurements, enables a description of a ducted noise source. The technique uses a finite difference representation of the homogeneous Helmholtz equation.

Baumeister, K. J.↗
NASA NTRS1979

Analysis and computation of viscous subsonic primary and secondary flows

A new approximate flow analysis, designed to enable numerical solution as an initial value problem, is developed for a wide class of viscous subsonic flows at high Reynolds number and in straight or smoothly curved three-dimensional flow geometries. The analysis is coordinate-independent and corrects an a priori known inviscid primary flow for viscous and thermal effects, secondary flows, total pressure distortion, internal flow blockage and pressure drop. Computed results include laminar solutions for three-dimensional boundary layer flow, fully viscous flow in circular arc ducts, and also flow in a curved duct shaped like a turbine blade passage.

Briley, W. R.↗
NASA NTRS1980

Mountain winds (revisited)

The prediction of extremely high wind speeds, at ground level on the downstream side of a mountain range, is possible by solving the initial value problem for a two-layered nonlinear shallow water model of the atmosphere. Three different numerical methods are described to find the solutions which may involve shocks: (1) the vonNeumann-Richtmyer artificial viscosity method, (2) a filtering scheme, and (3) a hybrid method.

Isaacson, E.↗
NASA NTRS1980

The dynamics of oceanic fronts. I - The Gulf Stream

The establishment and maintenance of the mean hydrographic properties of large-scale density fronts in the upper ocean is considered. The dynamics is studied by posing an initial value problem starting with a near-surface discharge of buoyant water with a prescribed density deficit into an ambient stationary fluid of uniform density; full time dependent diffusion and Navier-Stokes equations are then used with constant eddy diffusion and viscosity coefficients, together with a constant Coriolis parameter. Scaling analysis reveals three independent scales of the problem including the radius of deformation of the inertial length, buoyancy length, and diffusive length scales. The governing equations are then suitably scaled and the resulting normalized equations are shown to depend on the Ekman number alone for problems of oceanic interest. It is concluded that the mean Gulf Stream dynamics can be interpreted in terms of a solution of the Navier-Stokes and diffusion equations, with the cross-stream circulation responsible for the maintenance of the front; this mechanism is suggested for the maintenance of the Gulf Stream dynamics.

Kao, T. W.↗
NASA NTRS1970

The dynamics of oceanic fronts. Part 1: The Gulf Stream

The establishment and maintenance of the mean hydrographic properties of large scale density fronts in the upper ocean is considered. The dynamics is studied by posing an initial value problem starting with a near surface discharge of buoyant water with a prescribed density deficit into an ambient stationary fluid of uniform density. The full time dependent diffusion and Navier-Stokes equations for a constant Coriolis parameter are used in this study. Scaling analysis reveals three independent length scales of the problem, namely a radius of deformation or inertial length scale, Lo, a buoyance length scale, ho, and a diffusive length scale, hv. Two basic dimensionless parameters are then formed from these length scales, the thermal (or more precisely, the densimetric) Rossby number, Ro = Lo/ho and the Ekman number, E = hv/ho. The governing equations are then suitably scaled and the resulting normalized equations are shown to depend on E alone for problems of oceanic interest. Under this scaling, the solutions are similar for all Ro. It is also shown that 1/Ro is a measure of the frontal slope. The governing equations are solved numerically and the scaling analysis is confirmed. The solution indicates that an equilibrium state is established. The front can then be rendered stationary by a barotropic current from a larger scale along-front pressure gradient. In that quasisteady state, and for small values of E, the main thermocline and the inclined isopycnics forming the front have evolved, together with the along-front jet. Conservation of potential vorticity is also obtained in the light water pool. The surface jet exhibits anticyclonic shear in the light water pool and cyclonic shear across the front.

Kao, T. W.↗
NASA NTRS1981

Prediction of laminar and turbulent primary and secondary flows in strongly curved ducts

The analysis is based on a primary secondary velocity decomposition in a given coordinate system, and leads to approximate governing equations which correct an a priori inviscid solution for viscous effects, secondary flows, total pressure distortion, heat transfer, and internal flow blockage and losses. Solution of the correction equations is accomplished as an initial value problem in space using an implicit forward marching technique. The overall solution procedure requires significantly less computational effort than Navier-Stokes algorithms. The solution procedure is effective even with the extreme local mesh resolution which is necessary to solve near wall sublayer regions in turbulent flow calculations. Computed solutions for both laminar and turbulent flow compared very favorably with available analytical and experimental results. The overall method appears very promising as an economical procedure for making detailed predictions of viscous primary and secondary flows in highly curved passages.

Kreskovsky, J. P.↗
NASA NTRS1980

Transmitted sound field due to an impulsive line acoustic source bounded by a plate followed by a vortex sheet

The propagation of sound due to a line acoustic source in the moving stream across a semiinfinite vortex sheet which trails from a rigid plate is examined in a linear theory for the subsonic case. A solution for the transmitted sound field is obtained with the aid of multiple integral transforms and the Wiener-Hopf technique for both the steady state (time harmonic) and initial value (impulsive source) situations. The contour of inverse transform and hence the decomposition of the functions are determined through causality and radiation conditions. The solution obtained satisfies causality and the full Kutta conditions. The transmitted sound field is composed of two waves in both the stady state and initial value problems. One is the wave scattered from the edge of the plate which is associated with the bow wave and the instability wave. These waves exist in the downstream sectors. The other is the wave transmitted through the vortex sheet which is also associated with the instability wave. Regional divisions of the transmitted sound field are identified.

Miura, T.↗
NASA NTRS1982

Design of supercritical cascades with high solidity

The method of complex characteristics of Garabedian and Korn was successfully used to design shockless cascades with solidities of up to one. A code was developed using this method and a new hodograph transformation of the flow onto an ellipse. This code allows the design of cascades with solidities of up to two and larger turning angles. The equations of potential flow are solved in a complex hodograph like domain by setting a characteristic initial value problem and integrating along suitable paths. The topology that the new mapping introduces permits a simpler construction of these paths of integration.

Sanz, J. M.↗
NASA NTRS1982

Two-and-a-half-dimensional magnetohydrodynamic turbulence

The homogeneous turbulence for which fluctuating magnetic fields and velocity fields are independent of one spatial coordinate but still possess all three components are studied to form a generalized two-and-half dimensional geometry. The integral of the z-vector potential and the magnetic helicity are shown to be ideal invariants, and the basic dynamical variables and the equations of uniform-density incompressible magnetohydrodynamics are defined. The possibility of simultaneous inverse cascades is considered, and Kolmogoroff dimensional analysis is employed to infer omnidirectional inverse cascade spectra. Implications of a selective decay hypothesis, where a tendency exists in the initial value problem for the ideal invariants directly cascaded to higher wavenumbers to be selectively dissipated, and the possibility of a dynamo action in the considered geometry, are examined.

Montgomery, D.↗
NASA NTRS1981

Analysis and computation of three-dimensional flow in strongly curved ducts

Numerical solutions are presented for three-dimensional laminar and turbulent flow in curved ducts of rectangular cross section and significant curvature. The analysis is based on a primary-secondary velocity decomposition in a given coordinate system, and leads to approximate governing equations which correct an a priori inviscid solution for viscous effects, secondary flows, total pressure distortion, heat transfer, and internal flow blockage and losses. Solution of the correction equations is accomplished as an initial-value problem in space using an implicit forward-marching technique. The overall solution procedure requires significantly less computational effort than Navier-Stokes algorithms. The present solution procedure is effective even with the extreme local mesh resolution which is necessary to resolve near-wall sublayer regions in turbulent flow calculations. Computed solutions for both laminar and turbulent flow compare very favorably with available analytical and experimental results.

Kreskovsky, J. P.↗
NASA NTRS1982

Design of supercritical cascades with high solidity

The method of complex characteristics of Garabedian and Korn has been successfully used to design shockless cascades with solidities of up to one. A new code has been developed using this method and a new hodograph transformation of the flow onto an ellipse. This new code allows the design of cascades with solidities of up to two and larger turning angles. The equations of potential flow are solved in a complex hodograph-like domain by setting a characteristic initial value problem and integrating along suitable paths. The topology that the new mapping introduces permits a simpler construction of these paths of integration.

Sanz, J. M.↗
NASA NTRS1983

Stress and Fracture Analyses Under Elastic-plastic and Creep Conditions: Some Basic Developments and Computational Approaches

A new hybrid-stress finite element algorith, suitable for analyses of large quasi-static deformations of inelastic solids, is presented. Principal variables in the formulation are the nominal stress-rate and spin. A such, a consistent reformulation of the constitutive equation is necessary, and is discussed. The finite element equations give rise to an initial value problem. Time integration has been accomplished by Euler and Runge-Kutta schemes and the superior accuracy of the higher order schemes is noted. In the course of integration of stress in time, it has been demonstrated that classical schemes such as Euler's and Runge-Kutta may lead to strong frame-dependence. As a remedy, modified integration schemes are proposed and the potential of the new schemes for suppressing frame dependence of numerically integrated stress is demonstrated. The topic of the development of valid creep fracture criteria is also addressed.

Reed, K. W.↗
NASA NTRS1983

Prediction of a propellant tank pressure history using state space methods

An analysis of the time response of a propellant supply system operating in the blowdown mode is presented. The supply system is part of a pump-fed propulsion system intended for use on interplanetary spacecraft. As such, the supply system must provide the pump with propellant at sufficient pressure to avoid pump cavitation. The system, consisting of the tank, the liquid propellant, the pressurant gas and propellant vapor mixture, and a film layer separating the liquid and vapor phases, is analyzed using the principles of mass and energy conservation. The resulting set of ordinary, coupled, nonlinear differential equations for the thermodynamic state variables is integrated as an initial value problem. The resulting histories of total pressure, propellant vapor pressure, propellant liquid temperature, film layer temperature, propellant vapor/pressurant gas temperature, propellant vapor mass, and propellant liquid mass enable the calculation of the net positive suction head available at the pump which determines the viability of the pump-fed system concept when operated in the blowdown mode.

Estey, P. N.↗
NASA NTRS1983

Computer simulation of auroral kilometric radiation

An investigation is carried out of the linear amplification and nonlinear saturation of electromagnetic waves associated with two components of electrons consisting of a cold Maxwellian background and a weakly relativistic population of electrons possessing a loss-cone distribution. The goal is to understand the mechanism responsible for the production of auroral kilometric radiation (AKR). A relativistic 1-2/2 dimensional electromagnetic particle simulation code is used in studying an initial value problem modeling the AKR source region. The simulation makes it possible to follow the growth of all modes of radiation, including the ordinary mode (O-mode), slow extraordinary mode (Z-mode), whistler mode and fast extraordinary mode (X-mode), past the point of saturation.

Wagner, J. S.↗
NASA NTRS1984

Boolean difference equations. I - Formulation and dynamic behavior

In many biological and physical systems, feedback mechanisms depend on a set of thresholds associated with the state variables. Each feedback has a characteristic time scale. It is suggested that delay-difference equations for Boolean-valued variables are an appropriate mathematical framework for such situations: the feedback thresholds result in the discrete, on-off character of the variables, and the interaction time scales of the feedbacks are expressed as delays. The initial-value problem for Boolean delay equations (B-Delta-Es) is formulated, and shown to have unique solutions for all times. Examples of periodic and aperiodic solutions are given. Aperiodic solutions have increasing complexity which depends on time t roughly as t to the l-1 power, l being the number of delays. Stability of solutions is defined, and some examples of stability analysis are given; additional stability questions are raised. The present formulation of (B-Delta-Es) is compared with related work and generalizations are suggested.

Dee, D.↗