SEARCH · Search NASA
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.
Explicit finite difference schemes for the computation of weak solutions of nonlinear scalar conservation laws is presented and analyzed. These schemes are uniformly second-order accurate and nonoscillatory in the sense that the number of extrema of the discrete solution is not increasing in time.
Report gives information on use of NASTRAN computer program in finite-element analysis of rotating flexible blades like in compressors and on turboprop engines. Predicts steady-state components of deflections and stresses under centrifugal forces, generates data for plots of natural frequency versus rotational speed, and provides vibration-mode data for calculations of flutter. Describes use of NASTRAN solution sequence 64 for geometrical nonlinear analysis and solution sequence 63 for determination of frequencies and vibrational-mode shapes. Includes sample problem with NASTRAN input data. Emphasizes key factors in analysis of rotating blades, such as setting angle and centrifugal softening effects. Combined analyses of solution sequences 64 and 63 reduces computer time and number of output listings, in comparison with separate analyses. In central-processing-unit time cut in half.
A method is presented for the approximate solution of the nonlinear equations transonic flow theory. Solutions are found for two-dimensional flows at a Mach number of 1 and for purely subsonic and purely supersonic flows. Results are obtained in closed analytic form for a large and significant class of nonlifting airfoils. At a Mach number of 1 general expressions are given for the pressure distribution on an airfoil of specified geometry and for the shape of an airfoil having a prescribed pressure distribution. Extensive comparisons are made with available data, particularly for a Mach number of 1, and with existing solutions.
Approximate solutions for three nonlinear orthotropic plate problems are presented: (1) a thick plate attached to a pad having nonlinear material properties which, in turn, is attached to a substructure which is then deformed; (2) a long plate loaded in inplane longitudinal compression beyond its buckling load; and (3) a long plate loaded in inplane shear beyond its buckling load. For all three problems, the two dimensional plate equations are reduced to one dimensional equations in the y-direction by using a one dimensional trigonometric approximation in the x-direction. Each problem uses different trigonometric terms. Solutions are obtained using an existing algorithm for simultaneous, first order, nonlinear, ordinary differential equations subject to two point boundary conditions. Ordinary differential equations are derived to determine the variable coefficients of the trigonometric terms.
The finite-element method is used to establish a numerical solution algorithm for the Navier-Stokes equations for two-dimensional flows of a viscous compressible fluid. Numerical experiments confirm the advection property for the finite-element equivalent of the nonlinear convection term for both unidirectional and recirculating flowfields. For linear functionals, the algorithm demonstrates good accuracy using coarse discretizations and h squared convergence with discretization refinement.
High order accurate centered flux approximations used in the computation of numerical solutions to nonlinear partial differential equations produce large oscillations in regions of sharp transitions. Here, we present a new class of filtering methods denoted by Essentially Nonoscillatory Least Squares (ENOLS), which constructs an upgraded filtered solution that is close to the physically correct weak solution of the original evolution equation. Our method relies on the evaluation of a least squares polynomial approximation to oscillatory data using a set of points which is determined via the ENO network. Numerical results are given in one and two space dimensions for both scalar and systems of hyperbolic conservation laws. Computational running time, efficiency, and robustness of method are illustrated in various examples such as Riemann initial data for both Burgers' and Euler's equations of gas dynamics. In all standard cases, the filtered solution appears to converge numerically to the correct solution of the original problem. Some interesting results based on nonstandard central difference schemes, which exactly preserve entropy, and have been recently shown generally not to be weakly convergent to a solution of the conservation law, are also obtained using our filters.
The paper analyzes the boundary layers formed in a rotating fluid by an oscillating flow over an infinite half plate, with particular attention paid to the effects of unsteadiness, the critical latitude effect and the structure of the solution to the boundary layer equations at resonance. The Navier-Stokes boundary layer equations are obtained through an asymptotic expansion with the incorporation of the Rossby and Ekman numbers and are analyzed as the sum of a nonlinear steady solution and a linearized unsteady solution. The solution is predominantly composed of two inertial wave vector components, one circularly polarized to the left and the other circularly polarized to the right. The problem considered here has relevance in oceanography and meteorology, with special reference to the unsteady atmospheric boundary layer.
A solution to second order in thickness is derived for harmonically oscillating two-dimensional airfoils in supersonic flow. For slow oscillations of an arbitrary profile, the result is found as a series including the third power of frequency. For arbitrary frequencies, the method of solution for any specific profile is indicated, and the explicit solution derived for a single wedge. Nonlinear thickness effects are found generally to reduce the torsional damping, and so enlarge the range of Mach numbers within which torsional instability is possible.
A proof is given of the existence and uniqueness of the solution to the automatic control problem with a nonlinear state equation of the form y' = f(t,y,u) and nonlinear operator controls u = U(y) acting onto the state function y which satisfies the initial condition y(t) = x(t) for t or = 0.
Approach to solution of nonlinear ordinary differential equations involves transfer functions based on recently-introduced Laplace-Borel and Fourier-Borel transforms. Main theorem gives transform of response of nonlinear system as Cauchy product of transfer function and transform of input function of system, together with memory effects. Used to determine responses of electrical circuits containing variable inductances or resistances. Also possibility of doing all noncommutative algebra on computers in such symbolic programming languages as Macsyma, Reduce, PL1, or Lisp. Process of solution organized and possibly simplified by algebraic manipulations reducing integrals in solutions to known or tabulated forms.
This paper presents details of a flutter and stability analysis of aerospace structures such as hypersonic vehicles. Both structural and aerodynamic domains are discretized by the common finite element technique. A vibration analysis is first performed by the STARS code employing a block Lanczos solution scheme. This is followed by the generation of a linear aerodynamic grid for subsequent linear flutter analysis within subsonic and supersonic regimes of the flight envelope; the doublet lattice and constant pressure techniques are employed to generate the unsteady aerodynamic forces. Flutter analysis is then performed for several representative flight points. The nonlinear flutter solution is effected by first implementing a CFD solution of the entire vehicle. Thus, a 3-D unstructured grid for the entire flow domain is generated by a moving front technique. A finite element Euler solution is then implemented employing a quasi-implicit as well as an explicit solution scheme. A novel multidisciplinary analysis is next effected that employs modal and aerodynamic data to yield aerodynamic damping characteristics. Such analyses are performed for a number of flight points to yield a large set of pertinent data that define flight flutter characteristics of the vehicle. This paper outlines the finite-element-based integrated analysis procedures in detail, which is followed by the results of numerical analyses of flight flutter simulation.
The successful high angle of attack (HAOA) operation of fighter aircraft will necessarily require the introduction of a new onboard control methodology that address the nonlinearity of the system when flown at the stall/poststall limits of the craft's flight envelope. As a precursor to this task, a researcher endeavored to familarize himself with the dynamics of one specific aircraft, the F-18, when it is flown at HAOA. This was accomplished by conducting a number of real time flight sorties using the NASA-Langley Research Center's F-18 simulator, which was operated with a pilot in the loop. In addition to developing a first hand familarity with the aircraft's dynamic characteristic at HAOA, work was also performed to identify the input/output operational footprint of the F-18's control surfaces. This investigator proposes to employ the nonlinear models of the plant identified this summer in a subsequent research effort that will make it possible to fly the F-18 effectively at poststall angles of attack. The controller design used there will rely on a new technique proposed by this investigator that provides for the automatic generation of online optimal control solutions for nonlinear dynamic systems.
A suitable grid system for complex three dimensional configurations such as wing/body/nacelle shape for the solution of nonlinear transonic flow problems was constructed. Two approaches were explored based on Thompson's body fitted coordinate concept. The most general approach was to divide the computational domain into multiple rectangular blocks where the configuration itself was also represented by a set of blocks whose structure follows the natural lines of the configuration. The block structured grid system was adaptable to complex configurations and gives good grid quality near physical corners. However, it introduced algorithm issues for the flow solution concerning the treatment of nonanalytic grid block boundaries and nonstandard grid cells. These issues were explored in relation to the grid generation. A more limited approach treats a wing/body configuration with only a single rectangular block in computational space. The issues involving nonstandard cells were avoided, but other limitations on grid resolution appear. Both a linear and a nonlinear system of grid generation equations were developed including methods of grid control.
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.
The development of an analytical technique used in the solution of nonlinear velocity-sensitive combustion instability problems is presented. The Galerkin method was used and proved successful. The pressure wave forms exhibit a strong second harmonic distortion and a variety of behaviors are possible depending on the nature of the combustion process and the parametric values involved. A one dimensional model provides insight into the problem by allowing a comparison of Galerkin solutions with more exact finite difference computations.
A two step hybrid analysis technique to predict the nonlinear steady state temperature distribution in structures and solids is presented. The technique is based on the regular perturbation expansion and the classical Bubnov-Galerkin approximation. The functions are obtained by using the regular perturbation method. These functions are selected as coordinate functions and the classical Bubnov-Galerkin technique is used to compute their amplitudes. The potential of the proposed hybrid technique for the solution of nonlinear thermal problems is discussed. The effectiveness of this technique is demonstrated by the effects of conduction, convection, and radiation modes of heat transfer. It is indicated that the hybrid technique overcomes the two major drawbacks of the classical techniques: (1) the requirement of using a small parameter in the regular perturbation method; and (2) the arbitrariness in the choice of the coordinate functions in the Bubnov-Galerkin technique. The proposed technique extends the range of applicability of the regular perturbation method and enhances the effectiveness of the Bubnov-Galerkin technique.
A method for improving the efficiency of nonlinear structural analysis by the use of global displacement functions is presented. The computer programs include options to define the global functions as input or let the program automatically select and update these functions. The program was applied to a number of structures: (1) 'pear-shaped cylinder' in compression, (2) bending of a long cylinder, (3) spherical shell subjected to point force, (4) panel with initial imperfections, (5) cylinder with cutouts. The sample cases indicate the usefulness of the procedure in the solution of nonlinear structural shell problems by the finite element method. It is concluded that the use of global functions for extrapolation will lead to savings in computer time.
Some results of studies of convergence and accuracy of finite element approximations of certain nonlinear problems encountered in finite elasticity are presented. A general technique for obtaining error bounds is also described together with an existence theorem. Numerical results obtained by solving a representative problem are also included.