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

Finite element solution of turbulent and reacting compressible three-dimensional flows

Results of detailed numerical evaluation of turbulent and reacting three-dimensional compressible flows are presented. These solutions were obtained using a finite element algorithm to solve the system of explicitly non-linear partial differential equations governing three-dimensional, confined unidirectional turbulent flow of a compressible, reacting fluid. Assessments are made of solution accuracy and convergence, including detailed comparisons to experimental data.

Baker, A. J.↗
NASA NTRS1977

Nonlinear propagation of a wave packet in a two-dimensional acoustically lined duct

The method of multiple scales is used to analyze the nonlinear effects of the gas motion and the acoustic lining material on the propagation and attenuation of a wave packet in a two-dimensional duct of uniform cross section. The partial differential equations describing the space and time variation of amplitudes and phases are obtained and used to show that both the monochromatic waves and the pure amplitude modulated waves are stable. The spatial attenuation of the pure amplitude modulated waves is found to be lower than that of monochromatic waves, while the temporal attenuation of the former waves has a minimum value near the resonant frequency. The nonlinearity shifts the wavenumber and frequency to higher values without changing phase speed for the pure phase modulated waves.

Tsai, M.-S.↗
NASA NTRS1978

Body-fitted coordinates systems transformations

Two computer programs generate two-dimensional body-fitted coordinate systems and coordinate transformation. Programs are useful in fields requiring accurate numerical representation of boundary conditions and accurate numerical solutions of partial differential equations.

Mastin, C. W.↗
NASA NTRS1978

Two-dimensional computer code for design assessment of SRM aft closure environment

A general interpolants method for constructing numerical analogs of the partial differential equations of continuum mechanics and which combines the best features of both finite element and finite difference methods was used in a two dimensional axisymmetric analysis of the space shuttle solid motor aft closure region flowfield. The numerical technique used is discussed as well as the results of the analysis. A steady state solution for the submerged nozzle nose region indicates the development of an area gas flow recirculation near the propellant burning surface boundary and in the region of the underside of the nozzle nose. This recirculation region remains spatially fixed for successive iterations (greater than 800) of the solution.

Freeman, J. A.↗
NASA NTRS1978

Calculation of the flow field in supersonic mixed-compression inlets at angle of attack using the three-dimensional method of characteristics with discrete shock wave fitting

The influence of molecular transport is included in the computation by treating viscous and thermal diffusion terms in the governing partial differential equations as correction terms in the method of characteristics scheme. The development of a production type computer program is reported which is capable of calculating the flow field in a variety of axisymmetric mixed-compression aircraft inlets. The results agreed well with those produced by the two-dimensional method characteristics when axisymmetric flow fields are computed. For three-dimensional flow fields, the results agree well with experimental data except in regions of high viscous interaction and boundary layer removal.

Vadyak, J.↗
NASA NTRS1978

A high order accurate finite element algorithm for high Reynolds number flow prediction

A Galerkin-weighted residuals formulation is employed to establish an implicit finite element solution algorithm for generally nonlinear initial-boundary value problems. Solution accuracy, and convergence rate with discretization refinement, are quantized in several error norms, by a systematic study of numerical solutions to several nonlinear parabolic and a hyperbolic partial differential equation characteristic of the equations governing fluid flows. Solutions are generated using selective linear, quadratic and cubic basis functions. Richardson extrapolation is employed to generate a higher-order accurate solution to facilitate isolation of truncation error in all norms. Extension of the mathematical theory underlying accuracy and convergence concepts for linear elliptic equations is predicted for equations characteristic of laminar and turbulent fluid flows at nonmodest Reynolds number. The nondiagonal initial-value matrix structure introduced by the finite element theory is determined intrinsic to improved solution accuracy and convergence. A factored Jacobian iteration algorithm is derived and evaluated to yield a consequential reduction in both computer storage and execution CPU requirements while retaining solution accuracy.

Baker, A. J.↗
NASA NTRS1978

Prediction of hydrodynamics and chemistry of confined turbulent methane-air frames in a two concentric tube combustor

A formulation of the governing partial differential equations for fluid flow and reacting chemical species in a two-concentric-tube 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 considered consists of two coaxial ducts. Concentric streams of gaseous fuel and air enter the inlet duct at one end; the flow then reverses and flows out through the outer duct. Two sample cases with specified inlet and boundary conditions are considered and the results are discussed.

Markatos, N. C.↗
NASA NTRS1978

Analysis of turbulent free jet hydrogen-air diffusion flames with finite chemical reaction rates

The nonequilibrium flow field resulting from the turbulent mixing and combustion of a supersonic axisymmetric hydrogen jet in a supersonic parallel coflowing air stream is analyzed. Effective turbulent transport properties are determined using the (K-epsilon) model. The finite-rate chemistry model considers eight reactions between six chemical species, H, O, H2O, OH, O2, and H2. The governing set of nonlinear partial differential equations is solved by an implicit finite-difference procedure. Radial distributions are obtained at two downstream locations of variables such as turbulent kinetic energy, turbulent dissipation rate, turbulent scale length, and viscosity. The results show that these variables attain peak values at the axis of symmetry. Computed distributions of velocity, temperature, and mass fraction are also given. A direct analytical approach to account for the effect of species concentration fluctuations on the mean production rate of species (the phenomenon of unmixedness) is also presented. However, the use of the method does not seem justified in view of the excessive computer time required to solve the resulting system of equations.

Sislian, J. P.↗
NASA NTRS1978

Computation of three-dimensional nozzle-exhaust flow fields with the GIM code

A methodology is introduced for constructing numerical analogs of the partial differential equations of continuum mechanics. A general formulation is provided which permits classical finite element and many of the finite difference methods to be derived directly. The approach, termed the General Interpolants Method (GIM), can combined the best features of finite element and finite difference methods. A quasi-variational procedure is used to formulate the element equations, to introduce boundary conditions into the method and to provide a natural assembly sequence. A derivation is given in terms of general interpolation functions from this procedure. Example computations for transonic and supersonic flows in two and three dimensions are given to illustrate the utility of GIM. A three-dimensional nozzle-exhaust flow field is solved including interaction with the freestream and a coupled treatment of the shear layer. Potential applications of the GIM code to a variety of computational fluid dynamics problems is then discussed in terms of existing capability or by extension of the methodology.

Spradley, L. W.↗
NASA NTRS1978

Boundary-value problem for a counterrotating electrical discharge in an axial magnetic field

An electrical discharge between two ring electrodes embedded in the mantle of a cylindrical chamber is considered, in which the plasma in the anode and cathode regions rotates in opposite directions under the influence of an external axial magnetic field. The associated boundary-value problem for the coupled partial differential equations describing the azimuthal velocity and radial current-density fields is solved in closed form. The velocity, current density, induced magnetic induction, and electric fields are presented for typical Hartmann numbers, magnetic Reynolds numbers, and geometry parameters. The discharge is shown to produce anodic and cathodic plasma sections rotating at speeds of the order 1,000,000 cm/sec for conventional magnetic field intensities. Possible application of the magnetoactive discharge as a plasma centrifuge for isotope separation is discussed.

Hong, S. H.↗
NASA NTRS1978

The relation between temporal and spatial stability in three-dimensional flows

An analysis is presented of the nonparallel spatial or temporal stability of three-dimensional incompressible, isothermal boundary-layer flows taking into account the transverse velocity component as well as the axial and crossflow variations of the mean flow. The method of multiple scales is used to derive partial differential equations that describe the axial and crossflow variations of the disturbance amplitude, phase and wavenumbers. This equation is used to derive the expressions that relate the temporal and spatial instabilities. These relations are functions of the complex group velocities. Moreover, this equation is used to derive the expression that relates the spatial amplification in any direction to a calculated amplification in any other direction. These relations are verified by numerical results obtained for two- and three-dimensional disturbances in two- and three-dimensional flows.

Nayfeh, A. H.↗
NASA NTRS1978

Fracture and crack growth in orthotropic laminates

An approximate solution is developed for the determination of the interlaminar normal and shear stresses in the vicinity of a crack in a three dimensional composite containing unidirectional linearly elastic fibers in an infinite linearly elastic matrix. In order to reduce the complexity of the formulation, certain assumptions are made as to the physically significant stresses to be retained. These simplifications reduce the partial differential equations of elasticity to differential-difference equations which are tractable using Fourier transform techniques. The potential for damaged or debonded zones to be generated by an embedded crack is discussed, and stress concentration factors for fibers near the crack are given. Detailed comparisons are made between the present solution, the analogous two dimensional problem, and corresponding shear-lag models.

Goree, J. G.↗
NASA NTRS1978

Programming effort analysis of the ELLPACK language

ELLPACK is a problem statement language and system for elliptic partial differential equations which is implemented by a FORTRAN preprocessor. ELLPACK's principal purpose is as a tool for the performance evaluation of software. However, it is used here as an example with which to study the programming effort required for problem solving. It is obvious that problem statement languages can reduce programming effort tremendously; the goal is to quantify this somewhat. This is done by analyzing the lengths and effort (as measured by Halstead's software science technique) of various approaches to solving these problems.

Rice, J. R.↗
NASA NTRS1978

Galactic civilizations: Population dynamics and interstellar diffusion

The interstellar diffusion of galactic civilizations is reexamined by potential theory; both numerical and analytical solutions are derived for the nonlinear partial differential equations which specify a range of relevant models, drawn from blast wave physics, soil science, and, especially, population biology. An essential feature of these models is that, for all civilizations, population growth must be limited by the carrying capacity of the environment. Dispersal is fundamentally a diffusion process; a density-dependent diffusivity describes interstellar emigration. Two models are considered: the first describing zero population growth (ZPG), and the second which also includes local growth and saturation of a planetary population, and for which an asymptotic traveling wave solution is found.

Newman, W. I.↗
NASA NTRS1979

Recent developments in finite element analysis for transonic airfoils

The prediction of aerodynamic forces in the transonic regime generally requires a flow field calculation to solve the governing non-linear mixed elliptic-hyperbolic partial differential equations. Finite difference techniques were developed to the point that design and analysis application are routine, and continual improvements are being made by various research groups. The principal limitation in extending finite difference methods to complex three-dimensional geometries is the construction of a suitable mesh system. Finite element techniques are attractive since their application to other problems have permitted irregular mesh elements to be employed. The purpose of this paper is to review the recent developments in the application of finite element methods to transonic flow problems and to report some recent results.

Hafez, M. M.↗
NASA NTRS1979

Development of a Model and Computer Code to Describe Solar Grade Silicon Production Processes

Mathematical models and computer codes based on these models, which allow prediction of the product distribution in chemical reactors for converting gaseous silicon compounds to condensed-phase silicon were developed. The following tasks were accomplished: (1) formulation of a model for silicon vapor separation/collection from the developing turbulent flow stream within reactors of the Westinghouse (2) modification of an available general parabolic code to achieve solutions to the governing partial differential equations (boundary layer type) which describe migration of the vapor to the reactor walls, (3) a parametric study using the boundary layer code to optimize the performance characteristics of the Westinghouse reactor, (4) calculations relating to the collection efficiency of the new AeroChem reactor, and (5) final testing of the modified LAPP code for use as a method of predicting Si(1) droplet sizes in these reactors.

Srivastava, R.↗
NASA NTRS1979

An artificial viscosity method for the design of supercritical airfoils

A numerical technique is presented for the design of two-dimensional supercritical wing sections with low wave drag. The method is a design mode of the analysis code H which gives excellent agreement with experimental results and is widely used in the aircraft industry. Topics covered include the partial differential equations of transonic flow, the computational procedure and results; the design procedure; a convergence theorem; and description of the code.

Mcfadden, G. B.↗
NASA NTRS1979

A simplified computer program for the prediction of the linear stability behavior of liquid propellant combustors

A program for predicting the linear stability of liquid propellant rocket engines is presented. The underlying model assumptions and analytical steps necessary for understanding the program and its input and output are also given. The rocket engine is modeled as a right circular cylinder with an injector with a concentrated combustion zone, a nozzle, finite mean flow, and an acoustic admittance, or the sensitive time lag theory. The resulting partial differential equations are combined into two governing integral equations by the use of the Green's function method. These equations are solved using a successive approximation technique for the small amplitude (linear) case. The computational method used as well as the various user options available are discussed. Finally, a flow diagram, sample input and output for a typical application and a complete program listing for program MODULE are presented.

Mitchell, C. E.↗