SEARCH · Engineering Papers
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.
The Reaction kinetics ablation program /REKAP/, developed to simulate ablation of various materials, provides mathematical formulations for computer programs which can simulate certain industrial processes. The programs are based on the use of nonsymmetrical difference equations that are employed to solve complex partial differential equation systems.
Wave-wave interactions are examined, particularly for significance in relation to Langmuir turbulence associated with type III solar radio emissions. The Zakharov equations are considered as two coupled nonlinear partial differential equations for the low-frequency density and the Langmuir envelope, and the electron and ion temperatures are regarded as equal. The Zhakharov equations are solved in one dimension, and display no collapse of a driven broadband pump, indicating that scattering instabilities are important in the saturation of a beam-plasma instability. A two-dimensional numerical simulation is developed to interpret the role of parametric instabilities, and it is concluded that the distortion of a wave packet is due to scattering instability, and not a change in the physics of direct collapse due to an artificial numerical damping. Further discussions are presented for comparisons of models of Langmuir turbulence in type III bursts.
We present an abstract approximation framework for estimation of parameters in nonlinear nonautonomous distributed systems. Specific examples involving linear spline approximations for retarded delay equations and cubic spline approximations for parabolic partial differential equations are discussed and shown to be included as special cases of our general framework.
Discussed is the implementation of a single algorithm on three parallel-vector computers. The algorithm is a relaxation scheme for the solution of the Cauchy-Riemann equations; a set of coupled first order partial differential equations. The computers were chosen so as to encompass a variety of architectures. They are: the MPP, and SIMD machine with 16K bit serial processors; FLEX/32, an MIMD machine with 20 processors; and CRAY/2, an MIMD machine with four vector processors. The machine architectures are briefly described. The implementation of the algorithm is discussed in relation to these architectures and measures of the performance on each machine are given. Simple performance models are used to describe the performance. These models highlight the bottlenecks and limiting factors for this algorithm on these architectures. Conclusions are presented.
Three new Runge-Kutta methods are presented for numerical integration of systems of linear inhomogeneous ordinary differential equations (ODES) with constant coefficients. Such ODEs arise in the numerical solution of the partial differential equations governing linear wave phenomena. The restriction to linear ODEs with constant coefficients reduces the number of conditions which the coefficients of the Runge-Kutta method must satisfy. This freedom is used to develop methods which are more efficient than conventional Runge-Kutta methods. A fourth-order method is presented which uses only two memory locations per dependent variable, while the classical fourth-order Runge-Kutta method uses three. This method is an excellent choice for simulations of linear wave phenomena if memory is a primary concern. In addition, fifth- and sixth-order methods are presented which require five and six stages, respectively, one fewer than their conventional counterparts, and are therefore more efficient. These methods are an excellent option for use with high-order spatial discretizations.
The Magnetic Field Solver computer program calculates the magnetic field generated by a group of collinear, cylindrical axisymmetric electromagnet coils. Given the current flowing in, and the number of turns, axial position, and axial and radial dimensions of each coil, the program calculates matrix coefficients for a finite-difference system of equations that approximates a two-dimensional partial differential equation for the magnetic potential contributed by the coil. The program iteratively solves these finite-difference equations by use of the modified incomplete Cholesky preconditioned-conjugate-gradient method. The total magnetic potential as a function of axial (z) and radial (r) position is then calculated as a sum of the magnetic potentials of the individual coils, using a high-accuracy interpolation scheme. Then the r and z components of the magnetic field as functions of r and z are calculated from the total magnetic potential by use of a high-accuracy finite-difference scheme. Notably, for the finite-difference calculations, the program generates nonuniform two-dimensional computational meshes from nonuniform one-dimensional meshes. Each mesh is generated in such a way as to minimize the numerical error for a benchmark one-dimensional magnetostatic problem.
Finite difference analysis of hyperbolic partial differential equations
As is well known, the full scattering amplitude can be expressed as an integral involving the complete scattering wave function. We have shown that the integral can be simplified and used in a practical way. Initial application to electron-hydrogen scattering without exchange was highly successful. The Schrodinger equation (SE) can be reduced to a 2d partial differential equation (pde), and was solved using the finite element method. We have now included exchange by solving the resultant SE, in the static exchange approximation. The resultant equation can be reduced to a pair of coupled pde's, to which the finite element method can still be applied. The resultant scattering amplitudes, both singlet and triplet, as a function of angle can be calculated for various energies. The results are in excellent agreement with converged partial wave results.
The differential equations describing the flap-lag-torsional motion of a flexible rotor blade including third-order nonlinearities were derived for hover and forward flight. Making use of the two boundary conditions, those equations were reduced to a set of three integro partial differential equations written in terms of the flexural deflections and the torsional variable.
The differential equations of motion, and boundary conditions, describing the flap-lead/lag-torsional motion of a flexible rotor blade with a precone angle and a variable pitch angle, which incorporates a pretwist, are derived via Hamilton's principle. The meaning of inextensionality is discussed. The equations are reduced to a set of three integro partial differential equations by elimination of the extension variable. The generalized aerodynamic forces are modelled using Greenberg's extension of Theodorsen's strip theory. The equations of motion are systematically expanded into polynomial nonlinearities with the objective of retaining all terms up to third degree. The blade is modeled as a long, slender, of isotropic Hookean materials. Offsets from the blade's elastic axis through its shear center and the axes for the mass, area and aerodynamic centers, radial nonuniformaties of the blade's stiffnesses and cross section properties are considered and the effect of warp of the cross section is included in the formulation.
As is well known, the full scattering amplitude can be expressed as an integral involving the complete scattering wave function. We have shown that the integral can be simplified and used in a practical way. Initial application to electron-hydrogen scattering without exchange was highly successful. The Schrodinger equation (SE), which can be reduced to a 2d partial differential equation (pde), was solved using the finite element method. We have now included exchange by solving the resultant SE, in the static exchange approximation, which is reducible to a pair of coupled pde's. The resultant scattering amplitudes, both singlet and triplet, calculated as a function of energy are in excellent agreement with converged partial wave results.
Research was undertaken to determine the mass required to protect a surface during atmospheric reentry by injecting a cold gas tangentially along the surface. A consideration of this problem revealed that the boundary layer equations could be employed to approximate the most important flow configurations which occur in film cooling. However, even with this simplification, the equations to be solved are complex, nonlinear, coupled partial differential equations. Due to the relative weakness of classical methods of boundary layer analysis, finite difference methods were chosen as the best approach to the formidable problem under consideration. Instead of blindly adopting a proposed technique and trying to force it to the problem, a more systematic study was deemed to be not only desirable but necessary. A relatively basic study of the application of finite difference methods was undertaken with the objective of developing an algorithm capable of obtaining accurate solutions to film cooling problems.
A fast direct (noniterative) 'Cauchy-Riemann Solver' is developed for solving the finite-difference equations representing systems of first-order elliptic partial differential equations in the form of the nonhomogeneous Cauchy-Riemann equations. The method is second-order accurate and requires approximately the same computer time as a fast cyclic-reduction Poisson solver. The accuracy and efficiency of the direct solver are demonstrated in an application to solving an example problem in aerodynamics: subsonic inviscid flow over a biconvex airfoil. The analytical small-perturbation solution contains singularities, which are captured well by the computational technique. The algorithm is expected to be useful in nonlinear subsonic and transonic aerodynamics.
A second order viscous layer solution procedure has been developed that does in a consistent way include curvature effects and the corresponding normal pressure gradients. In the present system, the normal momentum equation is retained. The parabolic system of nonlinear partial differential equations is converted by linear finite differencing procedures to a system of linear algebraic equations and solved in primitive coordinates. The solutions have been shown to give smooth stable distributions for all the variables, most particularly the normal velocity which plays an important role in the interaction procedure. An algorithm for matching the viscous layer solution with a rotational characteristics outer solution has been developed.
It is noted that the nonlinear partial differential equation for the perturbation velocity potential and boundary conditions describing steady inviscid compressible transonic flow past a thin two-dimensional airfoil can be transformed into a singular integrodifferential equation and that differentiation of the latter yields an integral equation. Two forms of this integral equation currently exist: one for the singularity that is enclosed in an infinitely long strip of vanishing thickness and the other for the singularity that is enclosed in a vanishing circle. In the present article, a more general integral equation is derived by enclosing the singularity in a vanishing rectangular cavity of arbitrary aspect ratio. The two existing forms of this equation are deduced as special cases distinguished by the respective values for the aspect ratio (infinity for the first form and unity for the second).
A rigorous mathematical model was obtained for the boundary-layer free-edge stress singularity in angleplied and crossplied fiber composite laminates. The solution was obtained using a method consisting of complex-variable stress function potentials and eigenfunction expansions. The required order of the boundary-layer stress singularity is determined by solving the transcendental characteristic equation obtained from the homogeneous solution of the partial differential equations. Numerical results obtained show that the boundary-layer stress singularity depends only upon material elastic constants and fiber orientation of the adjacent plies. For angleplied and crossplied laminates the order of the singularity is weak in general.
The application of conforming mixed finite element methods to obtain approximate solutions of linearized Navier-Stokes equations is examined. Attention is given to the convergence rates of various finite element approximations of the pressure and the velocity field. The optimality of the convergence rates are addressed in terms of comparisons of the approximation convergence to a smooth solution in relation to the best approximation available for the finite element space used. Consideration is also devoted to techniques for efficient use of a Gaussian elimination algorithm to obtain a solution to a system of linear algebraic equations derived by finite element discretizations of linear partial differential equations.
Basic concepts associated with the numerical solution of elliptic partial differential equations are introduced, and procedures used to solve the full potential equation for transonic flow fields are discussed. Governing equations, classical relaxation schemes and concepts regarding transonic, full potential equation algorithms are covered. The equation transformation and grid generation procedures; full potential spatial differencing schemes; full potential iteration schemes, emphasizing convergence acceleration; and three dimensional applications are presented.