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 program, first written in REDUCE, was implemented in MACSYMA with several additional techniques which are explained. By utilizing many methods, the program can obtain closed-form and series solutions to a large class of linear and nonlinear problems. One of the techniques developed, reduction to a differential equation, has not previously appeared in the literature in the general form described.
The matrix solving algorithm of any finite element algorithm is extremely important since solution of the matrix equations requires a large amount of elapse time due to null calculations and excessive input/output operations. An alternate method of solving the matrix equations is presented. A symbolic processing step followed by numeric solution yields the solution very rapidly and is especially useful for nonlinear problems.
Linear acoustic calculations are compared with previously reported data for a small-scale hovering rotor operated at high tip Mach numbers. A detailed calculated description of the distributions of blade surface pressure and shear stress due to skin friction is presented. The noise due to skin friction and loading, in the rotor disk plane, is small compared to thickness noise. The basic conclusions of Boxwell et al about the importance of nonlinear effects are upheld. Some approximations involved in the current theories for the inclusion of nonlinear effects are discussed. Using a model nonlinear problem, it is shown that to use the acoustic analogy, good knowledge of the flowfield is required.
The hydrodynamic picture of the formation of primordial black holes (PBH) at the early stages of expansion of the Universe is considered. It is assumed that close to singularity, expansion occurs in a quasi-isotropic way. Using an EVM, a spherically symmetrical nonlinear problem of the evolution of primary strong deviation from the Fridman solution was solved. What these deviations must be, so that the formation of PBH occurred was clarified. Attention was devoted to the role of pressure gradients. It is pointed out that at the moment of formation of PBH, only a small part of matter enters into it, primarily the component of perturbation. It is also pointed out that at this moment, the mass of PBH essentially is smaller than the mass considered within the cosmic horizon. The possibility of changing the mass of the PBH as a result of accretion is analyzed.
An explicit time-advancing scheme for the spectral solution of parabolic equations is presented. Several two-dimensional examples are considered, including convection-diffusion and nonlinear problems, under various boundary conditions. Numerical evidence demonstrates the efficiency and accuracy of the spectral approach.
The development of a multiple-mass, nonlinear, finite-element analytical model - the Large AMplitude Propellant Slosh model (LAMPS) in two- and three-dimensional versions, is described. The model is used in predicting the forces on the Space Shuttle external tank caused by large amplitude motion of propellant fluid. Comparisons between measured propellant reorientation forces and those predicted by the LAMPS models are presented, and it is concluded that the model provides a cost effective, experimentally verified approach to a complicated nonlinear problem.
The theory, variational formulation and finite element analysis of a class of nonlinear problems arising in optimization, thin membrane deflection and compressible flows are developed. Error estimates are summarized, and representative numerical results are presented. Particular theoretical questions pertaining to the existence and uniqueness of solutions are addressed, and their implications for such applications of finite element analysis as the study of compressible flows indicated.
One of the objectives of the present investigation is to identify classes of equivalent mixed models and reduced/selective integration displacement models for curved shell structures. A second objective is concerned with the identification of the spurious modes exhibited by various mixed models and their equivalent reduced integration displacement models. The merits of using mixed models are also evaluated. The mixed elements developed in the current investigation differ from those presented by Noor and Hartley (1977) by the fact that the stress resultants are discontinuous at interelement boundaries and, therefore, are eliminated on the element level. The high accuracy and effectiveness of the elements developed is demonstrated by means of numerical examples, taking into account three geometrically nonlinear problems of shallow shells.
A slowly convergent stationary iterative process can be accelerated by explicitly annihilating (i.e., eliminating) the dominant eigenvector component of the error. The dominant eigenvalue or complex pair of eigenvalues can be estimated from the solution during the iteration. The corresponding eigenvector or complex pair of eigenvectors can then be annihilated by applying an explicit Richardson process over the basic iterative method. This can be done entirely in real arithmetic by analytically combining the complex conjugate annihilation steps. The technique is applied to an implicit algorithm for the calculation of two dimensional steady transonic flow over a circular cylinder using the equations of compressible inviscid gas dynamics. This demonstrates the use of explicit annihilation on a nonlinear problem.
A method is presented for performing efficient and stable finite element calculations of heat conduction with quadrilaterals using one-point quadrature. The stability in space is obtained by using a stabilization matrix which is orthogonal to all linear fields and its magnitude is determined by a stabilization parameter. It is shown that the accuracy is almost independent of the value of the stabilization parameter over a wide range of values; in fact, the values 3, 2, and 1 for the normalized stabilization parameter lead to the 5-point, 9-point finite difference, and fully integrated finite element operators, respectively, for rectangular meshes and have identical rates of convergence in the L2 norm. Eigenvalues of the element matrices, which are needed for stability limits, are also given. Numerical applications are used to show that the method yields accurate solutions with large increases in efficiency, particularly in nonlinear problems.
Difficulties in solution methodology to be used to deal with the potentially higher nonlinear rotor equations when dynamic coupling is included. A solution methodology is selected to solve the nonlinear differential equations. The selected method was verified to give good results even at large nonlinearity levels. The transfer matrix methodology is extended to the solution of nonlinear problems.
The mechanics of materials and material characterization are considered, taking into account micromechanics, the behavior of steel structures at elevated temperatures, and an anisotropic plasticity model for inelastic multiaxial cyclic deformation. Other topics explored are related to advances and trends in finite element technology, classical analytical techniques and their computer implementation, interactive computing and computational strategies for nonlinear problems, advances and trends in numerical analysis, database management systems and CAD/CAM, space structures and vehicle crashworthiness, beams, plates and fibrous composite structures, design-oriented analysis, artificial intelligence and optimization, contact problems, random waves, and lifetime prediction. Earthquake-resistant structures and other advanced structural applications are also discussed, giving attention to cumulative damage in steel structures subjected to earthquake ground motions, and a mixed domain analysis of nuclear containment structures using impulse functions.
Massively parallel computers, as considered in this investigation, are not yet available. However, a large-scale parallel computer cannot usefully be designed before the hypothetical algorithms which will employ it are studied. Most of the studies of parallel partial differential equations (PDE) solvers are based on solution techniques much slower (on sequential machines) than multigrid methods. Multigrid methods are highly parallelizable. Each of their processes can simultaneously be performed at all grid points. The present investigation is concerned with a preliminary exploration of the potential of multigrid, or, more generally, Multi-Level Adaptive Techniques (MLAT) on computers with many processors. Basic processes are considered, taking into account coarse-grid approximation, relaxation, coarse-grid corrections, full multigrid algorithms, nonlinear problems and eigenvalue problems, fine-to-coarse correction, and chains of problems. Details of parallel multigrid processing are also examined.
The globally convergent quasi-Newton minimization algorithm and the homotopy algorithms are discussed in detail and their effectiveness in solving certain classes of highly nonlinear problems of structural analysis is demonstrated. The application of the double dogleg strategy controls the directions and step-lengths of the quasi-Newtonian algorithm and overcomes the problem of nonpositive definite Hessians being produced during the iteration process. The algorithms are applied to a centrally loaded clamped beam, the snap-through of a shallow arch, and a shallow reticulated dome.
Conference paper discusses computational requirements for finiteelement analysis using quasi-linear approach to nonlinear problems. Paper evaluates computational efficiency of different computer architecturtural types in terms of relative cost and computing time.
A method is presented for performing efficient and stable finite element calculations of heat conduction with quadrilaterals using one-point quadrature. The stability in space is obtained by using a stabilization matrix which is orthogonal to all linear fields and its magnitude is determined by a stabilization parameter. It is shown that the accuracy is almost independent of the value of the stabilization parameter over a wide range of values; in fact, the values 3, 2 and 1 for the normalized stabilization parameter lead to the 5-point finite difference, 9-point finite difference and fully integrated finite element operators, respectively, for rectangular meshes; numerical experiments reported here show that the three have identical rates of convergence in the L2 norm. Eigenvalues of the element matrices, which are needed for stability limits, are also given. Numerical applications are used to show that the method yields accurate solutions with large increases in efficiency, particularly in nonlinear problems.
A computationally efficient method has been developed for highly nonlinear problems in which radiative transfer is an important aspect of the heating and cooling of the medium. An approximate probabilistic radiative transfer equation is derived for one-dimensional plane-parallel atmospheres of finite or semi-infinite extent, for both spectral lines and bound-free continua. Boundary conditions, accuracy, escape probabilities, and practical aspects of complete linearization, are discussed. The method is accurate to a few tens of percent of a wide variety of realistic problems in which frequency redistribution of scattered photons dominates the transfer and escape of radiation.
A slowly convergent stationary iterative process can be accelerated by explicitly annihilating (i.e., eliminating) the dominant eigenvector component of the error. The dominant eigenvalue or complex pair of eigenvalues can be estimated from the solution during the iteration. The corresponding eigenvector or complex pair of eigenvectors can then be annihilated by applying an explicit Richardson process over the basic iterative method. This can be done entirely in real arithmetic by analytically combining the complex conjugate annihilation steps. The technique is applied to an implicit algorithm for the calculation of two dimensional steady transonic flow over a circular cylinder using the equations of compressible inviscid gas dynamics. This demonstrates the use of explicit annihilation on a nonlinear problem.