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 research is reported concerning the development of a three-dimensional nonlinear nozzle admittance relation to be used as a boundary condition in the nonlinear combustion instability theories for liquid propellant rocket engines. The derivation of the nozzle wave equation and the application of the Galerkin method are discussed along with the nozzle response.
It is demonstrated that the problem of elasto-plastic finite deformation is governed by a quasi-linear model irrespective of deformation magnitude. This feature follows from the adoption of a rate viewpoint toward finite deformation analysis in an Eulerian reference frame. Objectivity of the formulation is preserved by introduction of a frame-invariant stress rate. Equations for piece-wise linear incremental finite element analysis are developed by application of the Galerkin method to the instantaneously linear governing differential equations of the quasi-linear model. Numerical solution capability has been established for problems of plane strain and plane stress. The accuracy of the numerical analysis is demonstrated by consideration of a number of problems of homogeneous finite deformation admitting comparative analytic solution. It is shown that accurate, objective numerical solutions can be obtained for problems involving dimensional changes of an order of magnitude and rotations of a full radian.
An approximate analytical technique has been developed for the solution of nonlinear three-dimensional, transverse and axial combustion-instability problems that are frequently observed in liquid-propellant rocket motors. This theory is an extension and generalization of previous analyses, which could analyze either transverse or axial instabilities in liquid combustors with quasi-steady nozzles, to the practical situations of three-dimensional instabilities in combustors with conventional DeLaval nozzles. Unlike the quasi-steady nozzle, the presence of a conventional nozzle imposes restrictions upon the behavior of both the amplitudes and phases of the oscillations at the nozzle entrance plane. The applicability of the solution technique, which is based on the Galerkin method, is demonstrated by analyzing the nonlinear stability of a cylindrical liquid-rocket combustor with uniform injection of propellants at one end and a conventional nozzle at the other end. Crocco's pressure-sensitive time-lag model is used to describe the unsteady combustion process.
The quasi-static, pressure characteristics of a gas-lubricated thrust bearing with shrouded, Rayleigh-step pads are determined for a time-varying film thickness. The axial response of the thrust bearing to an axial forcing function or an axial rotor disturbance is investigated by treating the gas film as a spring having nonlinear restoring and damping forces. These forces are related to the film thickness by a power relation. The nonlinear equation of motion in the axial mode is solved by the Ritz-Galerkin method as well as the direct, numerical integration. Results of the nonlinear response by both methods are compared with the response based on the linearized equation. Further, the gas-film instability of an infinitely wide Rayleigh step thrust pad is determined by solving the transient Reynolds equation coupled with the equation of the motion of the pad. Results show that the Rayleigh-step geometry is very stable for bearing number A up to 50. The stability threshold is shown to exist only for ultrahigh values of Lambda equal to or greater than 100, where the stability can be achieved by making the mass heavier than the critical mass.
A solution procedure was developed using linear small deflection theory for the flutter of simply supported laminated plates. For such plates, the bending and extensional governing equations are coupled and have cross-stiffness terms which do not appear in classical plate theory. An extended Galerkin method is used to obtain approximate solutions to the governing equations, and the aerodynamic pressure loading used in the analysis is that given by linear piston theory with flow at arbitrary cross-flow angle. A limited parametric study was conducted for typical laminated composite plates. The calculations show that both the bending-extensional coupling and the cross-stiffness terms have a large destabilizing effect on flutter. Since classical plate theory does not consider bending-extensional coupling and cross stiffness terms, it usually gives inaccurate and nonconservative flutter boundaries for laminated plates.
A three dimensional, nonlinear nozzle admittance relation is developed by solving the wave equation describing finite amplitude oscillatory flow inside the subsonic portion of a choked, slowly convergent axisymmetric nozzle. This nonlinear nozzle admittance relation is then used as a boundary condition in the analysis of nonlinear combustion instability in a cylindrical liquid rocket combustor. In both nozzle and chamber analyses solutions are obtained using the Galerkin method with a series expansion consisting of the first tangential, second tangential, and first radial modes. Using Crocco's time lag model to describe the distributed unsteady combustion process, combustion instability calculations are presented for different values of the following parameters: (1) time lag, (2) interaction index, (3) steady-state Mach number at the nozzle entrance, and (4) chamber length-to-diameter ratio. In each case, limit cycle pressure amplitudes and waveforms are shown for both linear and nonlinear nozzle admittance conditions. These results show that when the amplitudes of the second tangential and first radial modes are considerably smaller than the amplitude of the first tangential mode the inclusion of nozzle nonlinearities has no significant effect on the limiting amplitude and pressure waveforms.
An anisotropic buckling and flutter analysis is developed with allowance for both bending-extensional coupling and bending-twisting coupling within the framework of linear small deflection theory for simply supported general laminated plates. The extended Galerkin method is used to obtain approximate solutions to the coupled governing equations. The effects of various anisotropic stiffness parameters on the static and dynamic stability of laminated plates are evaluated, with particular emphasis on assessing the range of applicability of classical orthotropic plate theory. It is shown that bending-extensional coupling and bending-twisting stiffness terms have a destabilizing effect on buckling and flutter, the effect being more pronounced for a small number of layers. For symmetric plates, the number of layers required for orthotropic plate theory to be applicable is generally less for the buckling problem than for flutter. For square plates, aligning the fibers with the direction of airflow over the plate surface results in the highest flutter dynamic pressure.
The Galerkin method is used to investigate the large-amplitude free vibrations of simply supported and clamped tapered beams of rectangular cross-section which are frequently encountered in practical structures. Two types of linear tapers are considered: breadth and depth tapers. Two solutions are obtained for each type of taper, using trigonometric and polynomial displacement distributions. Frequency-amplitude relationships are obtained for all these cases for the fundamental flexural mode. The results indicate good agreement between trigonometric and polynomial solutions for all cases. The nonlinearity is noted to be always of the hardening type and, as expected, it is severe for beams with depth taper as compared with beams with breadth taper.
The stability of an electrically conducting fluid subjected to two dimensional disturbance was investigated. A physical system consisting of two parallel infinite vertical plates which are thermally insulated was studied. An external magnetic field of constant strength was applied to normal plates. The fluid was heated from below so that a steady temperature gradient was maintained in the fluid. The governing equations were derived by perturbation technique, and solutions were obtained by a modified Galerkin method. It was found that the presence of the magnetic field increases the stability of the physical system and instability can occur in the form of neutral or oscillatory instability.
An integral-eigenvalue problem is formulated for a crossed rectangular waveguide and solved numerically by applying the Ritz-Galerkin method. Theoretical formulas for determining cutoff frequencies and modal-field expressions are obtained for the specific case of a symmetrical rectangular waveguide, cutoff frequencies are calculated numerically, and the results are verified by comparison with available experimental data. The modal fields are expressed in terms of Fourier series for both TE and TM modes. It is found that the bandwidth can be increased to a maximum of 38% when the waveguide dimensions are properly selected and that the numerical results are in agreement with those computed by the method of partial regions. Some practical applications of the modal-field equations are briefly noted.
The generalized Galerkin method (or the method of integral relations) is applied to the type of problem described by quasilinear parabolic equations. As an example the problem of nonlinear transient slab diffusion with a general reservoir boundary condition is worked out. The integral relations are given for an arbitrary number of strips, and solutions using up to seven strips have been obtained in order to investigate the convergence of the method.
A multimodal finite-element analysis for two-dimensional and axisymmetric ducts with flow has been adopted to predict acoustic transmission in nonuniform ducts carrying a nonuniform flow (i.e., turbofan aircraft engine conditions). The analysis is based on a Galerkin method with a finite-element discretization, and multimodal matching at the duct ends. The element mesh and the number of modes required at each end are found to have a strong influence on the analysis. The numerical results compare well with results obtained from other techniques, including the method of weighted residuals.
The problems of combustion instability in an annular combustion chamber are investigated. A modified Galerkin method was used to produce a set of modal amplitude equations from the general nonlinear partial differential acoustic wave equation. From these modal amplitude equations, the two variable perturbation method was used to develop a set of approximate equations of a given order of magnitude. These equations were modeled to show the effects of velocity sensitive combustion instabilities by evaluating the effects of certain parameters in the given set of equations. By evaluating these effects, parameters which cause instabilities to occur in the combustion chamber can be ascertained. It is assumed that in the annular combustion chamber, the liquid propellants are injected uniformly across the injector face, the combustion processes are distributed throughout the combustion chamber, and that no time delay occurs in the combustion processes.
High frequency combustion instability problems in a liquid fuel annular combustion chamber are examined. A modified Galerkin method was used to produce a set of modal amplitude equations from the general nonlinear partial differential acoustic wave equation in order to analyze the problem of instability. From these modal amplitude equations, the two variable perturbation method was used to develop a set of approximate equations of a given order of magnitude. These equations were modeled to show the effects of velocity sensitive combustion instabilities by evaluating the effects of certain parameters in the given set of equations.
A study of the response of a clamped panel to a supersonic turbulent boundary layer, based on a recently developed theoretical model, is presented. It is demonstrated that the model incorporates the effect of coupling between the panel motion and the flow of the surrounding fluid. Further, a Ritz-Galerkin method is used to obtain approximate solutions for the statistics of the panel response to the turbulence. The numerical results are compared with previous experimental data and the theoretical model is assessed.
The aeroelastic stability and response problem of the coupled flap-lag-torsional dynamics of a hingeless rotor blade in forward flight is treated in a comprehensive manner. The spatial dependence of the partial differential, nonlinear, equations of motion is discretized using a multimodal Galerkin method. The aeroelastic problem is coupled with the trim state of the helicopter obtained from improved, representative, trim procedures. The nonlinear time dependent equilibrium position, or response, about which the equations are linearized is obtained by solving a sequence of linear periodic response problems, using quasi-linearization. Numerous results illustrating blade behavior in forward flight are presented.
Definite triple integrals encountered in applying the Galerkin method to the problem of heat and mass transfer across rectangular enclosures are discussed. Rather than evaluating them numerically, the technique described by Reid and Harris (1958) was extended to obtain the exact solution of the integrals. In the process, four linear simultaneous equations with triple integrals as unknowns were obtained. These equations were then solved exactly to obtain the closed form solution. Since closed form representations of this type have been shown to be useful in solving nonlinear hydrodynamic problems by series expansion, the integrals are presented here in general form.
Recent developments with finite-difference techniques are emphasized. The quotation marks reflect the fact that any finite discretization procedure can be included in this category. Many so-called finite element collocation and galerkin methods can be reproduced by appropriate forms of the differential equations and discretization formulas. Many of the difficulties encountered in early Navier-Stokes calculations were inherent not only in the choice of the different equations (accuracy), but also in the method of solution or choice of algorithm (convergence and stability, in the manner in which the dependent variables or discretized equations are related (coupling), in the manner that boundary conditions are applied, in the manner that the coordinate mesh is specified (grid generation), and finally, in recognizing that for many high Reynolds number flows not all contributions to the Navier-Stokes equations are necessarily of equal importance (parabolization, preferred direction, pressure interaction, asymptotic and mathematical character). It is these elements that are reviewed. Several Navier-Stokes and parabolized Navier-Stokes formulations are also presented.