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 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
NASA NTRS1943

The Heat Transfer to a Plate in Flow at High Speed

The heat transfer in the laminar boundary layer of a heated plate in flow at high speed can be obtained by integration of the conventional differential equations of the boundary layer, so long as the material values can be regarded as constant. This premise is fairly well satisfied at speeds up to about twice the sonic speed and at not excessive temperature rise of the heated plate. The general solution of the equation includes Pohlhausen's specific cases of heat transfer to a plate at low speeds and of the plate thermometer. The solution shows that the heat transfer coefficient at high speed must be computed with the same equation as at low speed, when it is referred to the difference of the wall temperature of the heated plate in respect to its "natural temperature." Since this fact follows from the linear structure of the differential equation describing the temperature field, it is equally applicable to the heat transfer in the turbulent boundary layer.

Eckert, E.↗
NASA NTRS2006

Semi-Analytic Reconstruction of Flux in Finite Volume Formulations

Semi-analytic reconstruction uses the analytic solution to a second-order, steady, ordinary differential equation (ODE) to simultaneously evaluate the convective and diffusive flux at all interfaces of a finite volume formulation. The second-order ODE is itself a linearized approximation to the governing first- and second- order partial differential equation conservation laws. Thus, semi-analytic reconstruction defines a family of formulations for finite volume interface fluxes using analytic solutions to approximating equations. Limiters are not applied in a conventional sense; rather, diffusivity is adjusted in the vicinity of changes in sign of eigenvalues in order to achieve a sufficiently small cell Reynolds number in the analytic formulation across critical points. Several approaches for application of semi-analytic reconstruction for the solution of one-dimensional scalar equations are introduced. Results are compared with exact analytic solutions to Burger s Equation as well as a conventional, upwind discretization using Roe s method. One approach, the end-point wave speed (EPWS) approximation, is further developed for more complex applications. One-dimensional vector equations are tested on a quasi one-dimensional nozzle application. The EPWS algorithm has a more compact difference stencil than Roe s algorithm but reconstruction time is approximately a factor of four larger than for Roe. Though both are second-order accurate schemes, Roe s method approaches a grid converged solution with fewer grid points. Reconstruction of flux in the context of multi-dimensional, vector conservation laws including effects of thermochemical nonequilibrium in the Navier-Stokes equations is developed.

Gnoffo, Peter A.↗
NASA NTRS1975

Special Perturbations Using Back-Correction Methods of Numerical Integration

A new class of linear multistep methods for numerical integration of differential equations is reported that permits satellite computation solutions to be corrected at certain points in the past as the integration advances in time. Algorithms have been developed for the solution of both first- and second-order differential equations. The back correction method appears to be more efficient than classical methods when dominant and perturbing forces can be separated.

Feagin, T.↗
NASA NTRS1977

Transient wave propagation in composite plates due to impact

Stress-wave propagation in a composite plate due to an impact has been examined using anisotropic elasticity theory. The plate is modelled as a number of identical anisotropic layers whose elastic moduli are given. Mindlin's (1959) approximate theory of plates is then applied to each layer to obtain a set of difference-differential equations of motion by using the interlaminar stresses and displacements as explicit variables. Dispersion relationships for harmonic waves are found when traction-free boundary conditions are applied to both surfaces of the plate. The difference-differential equations are reduced to difference equations via integral transforms. With given impact boundary conditions these equations are solved for an arbitrary number of layers in the plate, and the transient propagation of a stress wave is calculated by means of a fast Fourier transform algorithm.

Sung Kim, B.↗
NASA NTRS1986

An approximate buckling analysis for rectangular orthotropic plates with centrally located cutouts

An approximate analysis for predicting buckling of rectangular orthotropic composite plates with centrally located cutouts is presented. In this analysis, prebuckling and buckling problems are converted from a two-dimensional to a one-dimensional system of linear differential equations with variable coefficients. The conversion is accomplished by expressing the displacements as series with each element containing a trigonometric function of one coordinate and a coefficient that is an arbitrary function of the other coordinate. Ordinary differential equations are then obtained from a variational principle. Analytical results obtained from the approximate analysis are compared with finite element analyses for isotropic plates and for specially orthotropic plates with central circular cutouts of various sizes. Experimental results for the specially orthotropic plates are also presented. In nearly all cases, the approximate analysis predicts the buckling mode shapes correctly and predicts the buckling loads to within a few percent of the finite element and experimental results.

Nemeth, M. P.↗
NASA NTRS1994

Gradient drift eigenmodes in the equatorial electrojet

The problem of kilometer-scale irregularities in the daytime equatorial electrojet is revisited by means of an eigenmode analysis of the gradient drift instability. Realistic physical parameters are used, including the modeled altitude variations of ion and electron collision frequencies and mobilities. The full fourth-order system of two coupled differential equations (each of second order) for the denisty and electrostatic potential perturbations is solved numerically by a relaxation technique. Under some approximations, the fourth-order system can be shown to reduce to a second-order differential equation for the perturbed potential or density. The latter is solved using a shooting technique and provides initial guesses for numerical solutions to the full problem. It is shown that the linear growth rate peaks for kilometer-scale waves, contrary to the findings of recent initial-value studies. This occurs because the equilibrium velocity shear is much more effective as a damping mechanism for short-wavelength modes than it is for the longer, kilometer-scale modes. These results provide a natural qualitative explanation for the observed dominance of kilometer-scale structures in the daytime electrojet spectrum.

Wang, X.-H.↗
NASA NTRS1978

Description of a computer program and numerical techniques for developing linear perturbation models from nonlinear systems simulations

A numerical technique was developed which generates linear perturbation models from nonlinear aircraft vehicle simulations. The technique is very general and can be applied to simulations of any system that is described by nonlinear differential equations. The computer program used to generate these models is discussed, with emphasis placed on generation of the Jacobian matrices, calculation of the coefficients needed for solving the perturbation model, and generation of the solution of the linear differential equations. An example application of the technique to a nonlinear model of the NASA terminal configured vehicle is included.

Dieudonne, J. E.↗
NASA NTRS1980

Costas loop tracking performance in the presence of antenna switching transients

The present paper examines the impact of phase, amplitude and Doppler transients on the tracking performance of a Costas loop, both in the presence and absence of noise; the loop's lock maintenance capability is also examined. Specific considerations involve analyzing second order loop transient behavior in the absence of noise and the presence of the three transients noted above; here, the operational differential equations are formulated and numerically solved. In the context of noise presence, a first order loop is treated in conjunction with amplitude and phase transients; toward this end the time dependent version of the Fokker-Planck partial differential equation is developed and numerically solved. Selected performance curves describe second order loop behavior in the presence of transients, whose values are of varying severity. First order loop behavior in the presence of noise is characterized by the illustration of suitably defined time dependent 'truncated probability density functions' and associated loss of lock probabilities.

Weinberg, A.↗
NASA NTRS1990

Influence of convection on microstructure

In eutectic growth, as the solid phases grow they reject atoms to the liquid. This results in a variation of melt composition along the solid/liquid interface. In the past, mass transfer in eutectic solidification, in the absence of convection, was considered to be governed only by the diffusion induced by compositional gradients. However, mass transfer can also be generated by a temperature gradient. This is called thermotransport, thermomigration, thermal diffusion or the Soret effect. A theoretical model of the influence of the Soret effect on the growth of eutectic alloys is presented. A differential equation describing the compositional field near the interface during unidirectional solidification of a binary eutectic alloy was formulated by including the contributions of both compositional and thermal gradients in the liquid. A steady-state solution of the differential equation was obtained by applying appropriate boundary conditions and accounting for heat flow in the melt. Following that, the average interfacial composition was converted to a variation of undercooling at the interface, and consequently to microstructural parameters. The results obtained show that thermotransport can, under certain circumstances, be a parameter of paramount importance.

Wilcox, William R.↗
NASA NTRS1974

Contribution to the theory of tidal oscillations of an elastic earth. External tidal potential

The differential equations of the tidal oscillations of the earth were established under the assumption that the interior of the earth is laterally inhomogeneous. The theory was developed using vectorial and dyadic symbolism to shorten the exposition and to reduce the differential equations to a symmetric form convenient for programming and for numerical integration. The formation of tidal buldges on the surfaces of discontinuity and the changes in the internal density produce small periodic variations in the exterior geopotential which are reflected in the motion of artificial satellites. The analoques of Love elastic parameters in the expansion of exterior tidal potential reflect the asymmetric and inhomogeneous structure of the interior of the earth.

Musen, P.↗
NASA NTRS1981

Transient dynamic analysis of high-speed lightly loaded cylindrical roller bearings. 2: Computer program and results

The governing differential equations of motion for a high speed cylindrical roller bearing are programmed for numerical solution and plotted output. This computer program has the capability of performing a two dimensional or three dimensional simulation. Two numerical solutions of the governing differential equations were obtained to simulate the motion of a roller bearing, one for the two dimensional equations of motion and one for the three dimensional equations of motion. Computer generated plots were obtained and present such data as roller/cage interaction forces, roller/race traction forces, roller/race relative slip velocities and cage angular speed over a nondimensional time equivalent to 1.2 revolutions of the inner race. Roller axial displacement, roller skew angle, and skew moment are also plotted for the three dimensional solution. The trajectory of the cage center is plotted for both the two dimensional and three dimensional solutions.

Conry, T. F.↗
NASA NTRS1986

Theoretical Foundation for Weld Modeling

Differential equations describe physics of tungsten/inert-gas and plasma-arc welding in aluminum. Report collects and describes necessary theoretical foundation upon which numerical welding model is constructed for tungsten/inert gas or plasma-arc welding in aluminum without keyhole. Governing partial differential equations for flow of heat, metal, and current given, together with boundary conditions relevant to welding process. Numerical estimates for relative importance of various phenomena and required properties of 2219 aluminum included

Traugott, S.↗
NASA NTRS1992

Examples Of Synthesis Of Dual-Shaped Reflectors

Report presents examples to demonstrate validity and utility of method of synthesis of offset dual-shape reflectors. Method of synthesis described by the authors in previous journal article. Current report reviews derivation of partial differential equations and iterative method of numerical solution. Discusses significance of starting point of numerical integration on each reflector surface; this point could be at center, on outer rim, or at interior point. Emphasizes that one of notable attributes of partial differential equations is speed with which they can be solved.

Galindo, Victor↗
NASA NTRS1995

CCD/CID Processors Would Offer Greater Precision

Charge-coupled-device/charge-injection-device (CCD/CID) data processors of proposed type offer advantages of massively parallel computational architecture and high computational speed typical of older CCD/CID data processors, but with increased precision. Useful in performing matrix vector multiplications in variety of applications, including solving partial differential equations, processing signal and image data, control computations, and neural-network simulations. Greater precision of proposed devices help to ensure accuracy in CCD/CID implementations of pseudospectral neural networks - particular class of artificial neural networks especially suited to solving nonlinear differential equations.

Barhen, Jacob↗
NASA NTRS2014

Badhwar-O'Neill 2011 Galactic Cosmic Ray Model Update and Future Improvements

The Badhwar-O'Neill Galactic Cosmic Ray (GCR) Model based on actual GCR measurements is used by deep space mission planners for the certification of microelectronic systems and the analysis of radiation health risks to astronauts in space missions. The BO GCR Model provides GCR flux in deep space (outside the earth's magnetosphere) for any given time from 1645 to present. The energy spectrum from 50 MeV/n - 20 GeV/n is provided for ions from hydrogen to uranium. This work describes the most recent version of the BO GCR model (BO'11). BO'11 determined the GCR flux at a given time applying an emperical time delay function to past sunspot activity. We describe the GCR measurement data used in the BO'11 update - modern data from BESS, PAMELA, CAPRICE, and ACE emphasized more than the older balloon data used for the previous BO model (BO'10). We look at the GCR flux for the last 24 solar minima and show how much greater the flux was for the cycle 24 minimum in 2010. The BO'11 Model uses the traditional, steady-state Fokker-Planck differential equation to account for particle transport in the heliosphere due to diffusion, convection, and adiabatic deceleration. It assumes a radially symmetrical diffusion coefficient derived from magnetic disturbances caused by sunspots carried outward by a constant solar wind. A more complex differential equation is now being tested to account for particle transport in the heliosphere in the next generation BO model. This new model is time-dependent (no longer a steady state model). In the new model, the dynamics and anti-symmetrical features of the actual heliosphere are accounted for so emperical time delay functions will no longer be required. The new model will be capable of simulating the more subtle features of modulation - such as the Sun's polarity and modulation dependence on gradient and curvature drift. This improvement is expected to significantly improve the fidelity of the BO GCR model. Preliminary results of its performance will be presented.

O'Neill, Pat M.↗
NASA NTRS2014

Badhwar-O'Neill 2011 Galactic Cosmic Ray Model Update and Future Improvements

The Badhwar-O'Neill Galactic Cosmic Ray (GCR) Model based on actual GR measurements is used by deep space mission planners for the certification of micro-electronic systems and the analysis of radiation health risks to astronauts in space missions. The BO GCR Model provides GCR flux in deep space (outside the earth's magnetosphere) for any given time from 1645 to present. The energy spectrum from 50 MeV/n-20 GeV/n is provided for ions from hydrogen to uranium. This work describes the most recent version of the BO GCR model (BO'11). BO'11 determines the GCR flux at a given time applying an empirical time delay function to past sunspot activity. We describe the GCR measurement data used in the BO'11 update - modern data from BESS, PAMELA, CAPRICE, and ACE emphasized for than the older balloon data used for the previous BO model (BO'10). We look at the GCR flux for the last 24 solar minima and show how much greater the flux was for the cycle 24 minimum in 2010. The BO'11 Model uses the traditional, steady-state Fokker-Planck differential equation to account for particle transport in the heliosphere due to diffusion, convection, and adiabatic deceleration. It assumes a radially symmetrical diffusion coefficient derived from magnetic disturbances caused by sunspots carried onward by a constant solar wind. A more complex differential equation is now being tested to account for particle transport in the heliosphere in the next generation BO model. This new model is time-dependent (no longer a steady state model). In the new model, the dynamics and anti-symmetrical features of the actual heliosphere are accounted for so empirical time delay functions will no longer be required. The new model will be capable of simulating the more subtle features of modulation - such as the Sun's polarity and modulation dependence on the gradient and curvature drift. This improvement is expected to significantly improve the fidelity of the BO GCR model. Preliminary results of its performance will be presented.

O'Neill, Pat M.↗
NASA NTRS1994

Nonlinear Asymptotic Integration Algorithms for One-dimensional Autonomous Dissipative First-order Odes

Nonlinear asymptotic integrators are applied to one-dimensional, nonlinear, autonomous, dissipative, ordinary differential equations. These integrators, including a one-step explicit, a one-step implicit, and a one- and two-step midpoint algorithm, are designed to follow the asymptotic behavior of a system approaching a steady state. The methods require that the differential equation be written in a particular asymptotic form. This is always possible for a one-dimensional equation with a globally asymptotic steady state. In this case, conditions are obtained to guarantee that the implicit algorithms are well defined. Further conditions are determined for the implicit methods to be contractive. These methods are all first order accurate, while under certain conditions the midpoint algorithms may also become second order accurate. The stability of each method is investigated and an estimate of the local error is provided.

Haslach, Henry W., Jr.↗