Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “partial 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.

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

At least 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
NASA NTRS1987

Automating the parallel processing of fluid and structural dynamics calculations

The NASA Lewis Research Center is actively involved in the development of expert system technology to assist users in applying parallel processing to computational fluid and structural dynamic analysis. The goal of this effort is to eliminate the necessity for the physical scientist to become a computer scientist in order to effectively use the computer as a research tool. Programming and operating software utilities have previously been developed to solve systems of ordinary nonlinear differential equations on parallel scalar processors. Current efforts are aimed at extending these capabilities to systems of partial differential equations, that describe the complex behavior of fluids and structures within aerospace propulsion systems. This paper presents some important considerations in the redesign, in particular, the need for algorithms and software utilities that can automatically identify data flow patterns in the application program and partition and allocate calculations to the parallel processors. A library-oriented multiprocessing concept for integrating the hardware and software functions is described.

Arpasi, Dale J.↗
NASA NTRS1987

Self-consistent inclusion of space-charge in the traveling wave tube

It is shown how the complete field of the electron beam may be incorporated into the transmission line model theory of the traveling wave tube (TWT). The fact that the longitudinal component of the field due to the bunched beam is not used when formulating the beam-to-circuit coupling equation is not well-known. The fundamental partial differential equation for the traveling wave field is developed and compared with the older (now standard) one. The equation can be solved numerically using the same algorithms, but now the coefficients can be updated continuously as the calculation proceeds down the tube. The coefficients in the older equations are primarily derived from preliminary measurements and some trial and error. The newer coefficients can be found by a recursive method, since each has a well defined physical interpretation and can be calculated once a reasonable first trial solution is postulated. The results of the new expression were compared with those of the older forms, as well as to a field theory model to show the ease in which a reasonable fit to the field prediction is obtained. A complete summary of the existing transmission line modeling of the TWT is given to explain the somewhat vague ideas and techniques in the general area of drifting carrier-traveling circuit wave interactions. The basic assumptions and inconsistencies of the existing theory and areas of confusion in the general literature are examined and hopefully cleared up.

Freeman, Jon C.↗
NASA NTRS1987

Automating the parallel processing of fluid and structural dynamics calculations

The NASA Lewis Research Center is actively involved in the development of expert system technology to assist users in applying parallel processing to computational fluid and structural dynamic analysis. The goal of this effort is to eliminate the necessity for the physical scientist to become a computer scientist in order to effectively use the computer as a research tool. Programming and operating software utilities have previously been developed to solve systems of ordinary nonlinear differential equations on parallel scalar processors. Current efforts are aimed at extending these capabilties to systems of partial differential equations, that describe the complex behavior of fluids and structures within aerospace propulsion systems. This paper presents some important considerations in the redesign, in particular, the need for algorithms and software utilities that can automatically identify data flow patterns in the application program and partition and allocate calculations to the parallel processors. A library-oriented multiprocessing concept for integrating the hardware and software functions is described.

Arpasi, Dale J.↗
NASA NTRS1991

A recursive approach to the equations of motion for the maneuvering and control of flexible multi-body systems

Interest lies in a mathematical formulation capable of accommodating the problem of maneuvering a space structure consisting of a chain of articulated flexible substructures. Simultaneously, any perturbations from the 'rigid body' maneuvering and any elastic vibration must be suppressed. The equations of motion for flexible bodies undergoing rigid body motions and elastic vibrations can be obtained conveniently by means of Lagrange's equations in terms of quasi-coordinates. The advantage of this approach is that it yields equations in terms of body axes, which are the same axes that are used to express the control forces and torques. The equations of motion are nonlinear hybrid differential quations. The partial differential equations can be discretized (in space) by means of the finite element method or the classical Rayleigh-Ritz method. The result is a set of nonlinear ordinary differential equations of high order. The nonlinearity can be traced to the rigid body motions and the high order to the elastic vibration. Elastic motions tend to be small when compared with rigid body motions.

Kwak, Moon K.↗
NASA NTRS1991

PDEMOD: Software for control/structures optimization

Because of the possibility of adverse interaction between the control system and the structural dynamics of large, flexible spacecraft, great care must be taken to ensure stability and system performance. Because of the high cost of insertion of mass into low earth orbit, it is prudent to optimize the roles of structure and control systems simultaneously. Because of the difficulty and the computational burden in modeling and analyzing the control structure system dynamics, the total problem is often split and treated iteratively. It would aid design if the control structure system dynamics could be represented in a single system of equations. With the use of the software PDEMOD (Partial Differential Equation Model), it is now possible to optimize structure and control systems simultaneously. The distributed parameter modeling approach enables embedding the control system dynamics into the same equations for the structural dynamics model. By doing this, the current difficulties involved in model order reduction are avoided. The NASA Mini-MAST truss is used an an example for studying integrated control structure design.

Taylor, Lawrence W., Jr.↗
NASA NTRS2017

Stabilization of Hypersonic Boundary Layers by Linear and Nonlinear Optimal Perturbations

The effect of stationary, finite-amplitude, linear and nonlinear optimal perturbations on the modal disturbance growth in a Mach 6 axisymmetric flow over a 7 deg. half-angle cone with 0:126 mm nose radius and 0:305 m length is investigated. The freestream parameters (M = 6, Re(exp 1) = 18 x 10(exp. 6) /m) are selected to match the flow conditions of a previous experiment in the VKI H3 hypersonic tunnel. Plane-marching parabolized stability equations are used in conjunction with a partial-differential equation based planar eigenvalue analysis to characterize the boundary layer instability in the presence of azimuthally periodic streaks. The streaks are observed to stabilize nominally planar Mack mode instabilities, although oblique Mack mode and first-mode disturbances are destabilized. Experimentally measured transition onset in the absence of any streaks correlates with an amplification factor of N = 6 for the planar Mack modes. For high enough streak amplitudes, the transition threshold of N = 6 is not reached by the Mack mode instabilities within the length of the cone; however, subharmonic first-mode instabilities, which are destabilized by the presence of the streaks, do reach N = 6 near the end of the cone. The highest stabilization is observed at streak amplitudes of approximately 20 percent of the freestream velocity. Because the use of initial disturbance profiles based on linear optimal growth theory may yield suboptimal control in the context of nonlinear streaks, the computational predictions are extended to nonlinear optimal growth theory. Results show that by using nonlinearly optimal perturbation leads to slightly enhanced stabilization of plane Mack mode disturbances as well as reduced destabilization of subharmonic first-mode disturbances.

Incompressible boundary layer↗
NASA NTRS1986

An efficient method for solving the steady Euler equations

An efficient numerical procedure for solving a set of nonlinear partial differential equations, the steady Euler equations, using Newton's linearization procedure is presented. A theorem indicating quadratic convergence for the case of differential equations is demonstrated. A condition for the domain of quadratic convergence Omega(2) is obtained which indicates that whether an approximation lies in Omega(2) depends on the rate of change and the smoothness of the flow vectors, and hence is problem-dependent. The choice of spatial differencing, of particular importance for the present method, is discussed. The treatment of boundary conditions is addressed, and the system of equations resulting from the foregoing analysis is summarized and solution strategies are discussed. The convergence of calculated solutions is demonstrated by comparing them with exact solutions to one and two-dimensional problems.

Liou, M.-S.↗
NASA NTRS1986

An efficient method for solving the steady Euler equations

An efficient numerical procedure for solving a set of nonlinear partial differential equations is given, specifically for the steady Euler equations. Solutions of the equations were obtained by Newton's linearization procedure, commonly used to solve the roots of nonlinear algebraic equations. In application of the same procedure for solving a set of differential equations we give a theorem showing that a quadratic convergence rate can be achieved. While the domain of quadratic convergence depends on the problems studied and is unknown a priori, we show that firstand second-order derivatives of flux vectors determine whether the condition for quadratic convergence is satisfied. The first derivatives enter as an implicit operator for yielding new iterates and the second derivatives indicates smoothness of the flows considered. Consequently flows involving shocks are expected to require larger number of iterations. First-order upwind discretization in conjunction with the Steger-Warming flux-vector splitting is employed on the implicit operator and a diagonal dominant matrix results. However the explicit operator is represented by first- and seond-order upwind differencings, using both Steger-Warming's and van Leer's splittings. We discuss treatment of boundary conditions and solution procedures for solving the resulting block matrix system. With a set of test problems for one- and two-dimensional flows, we show detailed study as to the efficiency, accuracy, and convergence of the present method.

Liou, M. S.↗
NASA NTRS2000

Polar Plate Theory for Orthogonal Anisotropy

Laminated fiber-reinforced (or filamentary) composites are used today for their high strength-to-weight and stiffness-to-weight ratios. However, because of the anisotropic behavior of composites, determining the response on a macroscopic scale is challenging. This is particularly evident in the evaluation of the governing differential equations of a circular disk with the fibers of the lamina oriented with rectilinear orthogonality. This includes any situation involving a composite plate of circular geometry in which out-of-plane displacements due to load are desired, such as fastener pull through loading of a composite plate. Current analysis techniques use numerical methods with rectilinear coordinate systems to solve problems with circular geometry. These analyses over predict plate stiffness by 20% and underpredict failure by 70%. Consequently, there is a need to transform classical composite plate theory to a polar coordinate system. In order to better analyze structures with circular geometries the classical composite plate equations are transformed into the plate equations for a rectilinearly anisotropic composite in polar coordinates. A composite plate is typically a laminate of fibers in rectilinear directions. Subsequent to the lay-tip the necessary geometry is cut out of a rectangular plate. In a similar manner, the derivation of the plate equation starts with the fundamental definitions of strain, displacement and curvature and incorporates the material property angular dependence into the equilibrium equations for a differential polar element. In the transformed state, the stiffness coefficients are no longer constant, adding to the complexity of the governing differential equations. This paper discusses the new derivation and evaluation of the plate equations for a circular composite disk with orthogonal rectilinear anisotropy. The resultant new three partial differential equations, which describe the circular anisotropic plate, can be used to evaluate out-of-plane displacements for given load conditions in design of composite Structures. Without the formation of these unique plate equations, design of such structures is not a precise engineering accomplishment due to the lack of precise design tools. Such structures are over designed to compensate and costly tests need to be performed after production to validate safety. Though impossible to present the mathematical solution within the limits of this paper, the application of the newly derived plate equation in its orthotropic form for a balanced symmetric laminate compared favorably with test results and surpassed the results of standard numerical methods. Within the limitations of the problem foundation the predictive model provides a mechanism for establishing out-of-plane deflection levels for a circular composite plate. Future work in this area should include incorporating the circular composite plate equations in Finite Element Models. The new anisotropic polar plate equations call be utilized as a design tool so that the design more accurately meets the requirements, thus taking full advantage of the weight savings of composites. These accomplishments have not previously been presented to the engineering community nor ever proposed to the IAF. Should the proposed paper be accepted for presentation, attendance is assured.

Bailey, Michelle D.↗
NASA NTRS1972

A finite difference method for predicting supersonic turbulent boundary layer flows with tangential slot injection

An implicit finite difference method has been applied to tangential slot injection into supersonic turbulent boundary layer flows. In addition, the effects induced by the interaction between the boundary layer displacement thickness and the external pressure field are considered. In the present method, three different eddy viscosity models have been used to specify the turbulent momentum exchange. One model depends on the species concentration profile and the species conservation equation has been included in the system of governing partial differential equations. Results are compared with experimental data at stream Mach numbers of 2.4 and 6.0 and with results of another finite difference method. Good agreement was generally obtained for the reduction of wall skin friction with slot injection and with experimental Mach number and pitot pressure profiles. Calculations with the effects of pressure interaction included showed these effects to be smaller than effects of changing eddy viscosity models.

Miner, E. W.↗
NASA NTRS1973

Supersonic turbulent boundary-layer flows with tangential slot injection.

An implicit finite-difference method has been developed for the solution of the compressible boundary-layer equations. This method is applied to tangential slot injection into supersonic turbulent boundary-layer flows. In addition, the effects induced by the interaction between the boundary-layer displacement thickness and the external pressure field are considered. Three different eddy viscosity models have been used to specify the turbulent momentum exchange. One model depends on the species concentration profile, and the species conservation equation has been included in the system of governing partial differential equations. For air injected into air, the freestream and injected gases are treated as separate species which have common fluid properties. Calculations were made and results were compared with experimental data at freestream Mach numbers of 2.4 and 6.0 and with results of another finite-difference method. Good agreement was obtained for the reduction of wall skin friction with slot injection.

Miner, E. W.↗
NASA NTRS1976

Implicit finite-difference procedures for the computation of vortex wakes

Implicit finite-difference procedures for the primitive form of the incompressible Navier-Stokes and the compressible Euler equations are used to compute vortex wake flows. The partial differential equations in strong conservation-law form are transformed to cluster grid points in regions with large changes in vorticity. In addition to clustering, fourth-order accurate, spatial difference operators are used to help resolve the flow-field gradients. The use of implicit time-differencing permits large time steps to be taken since temporal variations are typically small. Computational efficiency is achieved by approximate factorization. Both two-dimensional and preliminary three-dimensional calculations are described and qualitatively compared with existing experimental data.

Steger, J. L.↗
NASA NTRS1980

An analysis of turbulent diffusion flame in axisymmetric jet

The kinetic theory of turbulent flow was employed to study the mixing limited combustion of hydrogen in axisymmetric jets. The integro-differential equations in two spatial and three velocity coordinates describing the combustion were reduced to a set of hyperbolic partial differential equations in the two spatial coordinates by a binodal approximation. The MacCormick's finite difference method was then employed for solution. The flame length was longer than that predicted by the flame-sheet analysis, and was found to be in general agreement with a recent experimental result. Increase of the turbulence energy and scale resulted in an enhancement of the combustion rate and, hence, in a shorter flame length. Details of the numerical method as well as of the physical findings are discussed.

Chung, P. M.↗
NASA NTRS1980

The multispecies modeling of the premixed, laminar steady-state ozone flame

Species dependent kinetic, transport and thermodynamic coefficients were employed in a one dimensional model of the premixed, laminar, steady state ozone flame. Convenient expressions for these coefficients are reported. They are based on independent measurements, no arbitrary parameters are used. The governing equations are solved using a relaxation technique and the partial differential equation package, PDECOL. Species and temperature profiles and the burning velocities are found over the range of initial ozone mole fraction of 0.25 to 1.00. The computed burning velocities are no more than 30% greater than the measurements of Streng and Grosses. Comparison with the computed results of Warnatz shows agreement within + or - 12%, even though quite different expressions for some of the kinetic coefficients were used. These differences are most obvious in the atomic oxygen and temperature profiles at an initial ozone mole fraction of unity.

Heimerl, J. M.↗
NASA NTRS1993

The telegraph equation in charged particle transport

We present a new derivation of the telegraph equation which modifies its coefficients. First, an infinite order partial differential equation is obtained for the velocity space solid angle-averaged phase-space distribution of particles which underwent at least a few collisions. It is shown that, in the lowest order asymptotic expansion, this equation simplifies to the well-known diffusion equation. The second-order asymptotic expansion for isotropic small-angle scattering results in a modified telegraph equation with a signal propagation speed of v(5/11) exp 1/2 instead of the usual v/3 exp 1/2. Our derivation of a modified telegraph equation follows from an expansion of the Boltzmann equation in the relevant smallness parameters and not from a truncation of an eigenfunction expansion. This equation is consistent with causality. It is shown that, under steady state conditions in a convecting plasma, the telegraph equation may be regarded as a diffusion equation with a modified transport coefficient, which describes a combination of diffusion and cosmic-ray inertia.

Gombosi, T. I.↗
NASA NTRS2010

Analysis of Instabilities in Non-Axisymmetric Hypersonic Boundary Layers Over Cones

Hypersonic flows over circular cones constitute one of the most important generic configurations for fundamental aerodynamic and aerothermodynamic studies. In this paper, numerical computations are carried out for Mach 6 flows over a 7-degree half-angle cone with two different flow incidence angles and a compression cone with a large concave curvature. Instability wave and transition-related flow physics are investigated using a series of advanced stability methods ranging from conventional linear stability theory (LST) and a higher-fidelity linear and nonlinear parabolized stability equations (PSE), to the 2D eigenvalue analysis based on partial differential equations. Computed N factor distribution pertinent to various instability mechanisms over the cone surface provides initial assessments of possible transition fronts and a guide to corresponding disturbance characteristics such as frequency and azimuthal wave numbers. It is also shown that strong secondary instability that eventually leads to transition to turbulence can be simulated very efficiently using a combination of advanced stability methods described above.

Li, Fei↗
NASA NTRS1986

Iterative methods for elliptic finite element equations on general meshes

Iterative methods for arbitrary mesh discretizations of elliptic partial differential equations are surveyed. The methods discussed are preconditioned conjugate gradients, algebraic multigrid, deflated conjugate gradients, an element-by-element techniques, and domain decomposition. Computational results are included.

Nicolaides, R. A.↗