COMPARISON OF EXACT AND APPROXIMATE SOLUTIONS FOR NONEQUILIBRIUM NOZZLE FLOWS
Comparison of exact and approximate solutions for nonequilibrium nozzle flows
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.
Comparison of exact and approximate solutions for nonequilibrium nozzle flows
Internal and distributed damping as possible methods for reducing the vibration response of turbomachine blades and theoretical methods for analyzing damped vibration were studied. It is demonstrated how the Ritz-Galerkin methods may be used to straightforwardly to analyze forced vibrations with damping. This is done directly without requiring the free vibration eigenfunctions. The Galerkin method is an effective technique for these types of problems. The Ritz method has the further advantage of not needing to satisfy the force type boundary conditions, which is particularly important for plates and shells. But proper functionals representing the forcing and damping terms must be developed, and this is done. Two types of damping--viscous and material (hysteretic) are considered. Both distributed and concentrated exciting forces are treated. Numerical results are obtained for cantilevered beams and rectangular plates. Studies showing the rates of convergence of the solutions are made. In the case of the cantilever beam, approximate solutions from the present methods are compared with the exact solutions.
The Ritz-Galerkin methods were used to obtain approximate solutions for free undamped, vibration problems. It is demonstrated that these same methods may be used straightforwardly to analyze forced vibrations with damping without requiring the free vibration eigenfunctions. It was shown that the Galerkin method is an effective technique for these types of problems. The Ritz method has the advantage that it does not need to satisfy the force-type boundary conditions, which is particularly important for plates and shells. Proper functionals representing the forcing and damping terms were developed. Two types of damping--viscous and material (hysteretic) are discussed. Distributed and concentrated exciting forces are treated. Numerical results are obtained for cantilevered beams and rectangular plates. The rates of convergence of the solutions are shown. Approximate solutions from the present methods are compared with the exact solutions for the cantilever beam.
The two-dimensional steady-state shape of a solidified region, such as a frost layer, was determined analytically for formation on a plate that is convectively cooled. The nonuniform shape of the layer is produced by exposure to a spatially nonuniform distribution of radiant energy. For high convective cooling the cooled wall approaches a uniform temperature, and an exact solution is obtained for the free boundary shape. For a lesser amount of convective cooling, the variation in temperature along the cooled boundary is treated by a boundary perturbation method. Some illustrative examples are given that show the effects of nonuniform heating and the magnitude of convective heat transfer at the cooled wall. Only one boundary condition is approximated by the perturbation solution; all of the other boundary conditions are satisfied exactly. The calculated results given here were found to satisfy the approximate boundary condition within a very small error.
This paper presents for the first time an exact solution to the fuel optimal propulsive reboost problem for flexible spacecraft. The spacecraft undergoes rigid-body motion and flexible-body motion and reboost is achieved propulsively through the use of reaction control jets. The exact fuel optimal solution to the associated minimization problem is found numerically based on an adaptive grid bisection search. The reboost of the floating harmonic oscillator reveals properties of the fuel optimal solution. Nondimensional plots of minimum fuel versus maneuver time expose the nature of the solution classes depending on the maneuver time. Some very interesting properties are observed, among them, the shifting of impulses to the smaller oscillator mass, and comparisons are made with near fuel optimal solutions obtained by other investigators.
The time domain code ASSPIN gives acousticians a powerful technique of advanced propeller noise prediction. Except for nonlinear effects, the code uses exact solutions of the Ffowcs Williams-Hawkings equation with exact blade geometry and kinematics. By including nonaxial inflow, periodic loading noise, and adaptive time steps to accelerate computer execution, the development of this code becomes complete.
A recently reported solution for stationary stability of a thermosolutal system with Soret diffusion is re-derived and examined using a symbolic computational package. Symbolic computational languages are well suited for such an analysis and facilitate a pragmatic approach that is adaptable to similar problems. Linearization of the equations, normal mode analysis, and extraction of the final solution are performed in a Mathematica notebook format. An exact solution is obtained for stationary stability in the limit of zero gravity. A closed form expression is also obtained for the location of asymptotes in relevant parameter, (Sm(sub c), Mac(sub c)), space. The stationary stability behavior is conveniently examined within the symbolic language environment. An abbreviated version of the Mathematica notebook is given in the Appendix.
This project is concerned with the possibility of using computers for the simulation of structural systems with the same degree of reliability as full scale physical experiments. Reliable numerical simulation will make it possible to reduce the costs of engineering and improve the quality of engineering decisions based on computed information. An error of idealization is an error between the actual physical quantities on which engineering decisions are based (e.g., maximum principal stress, first natural frequency, etc.) and the same data corresponding to the exact solution of the mathematical model. An error of discretization is an error between the quantities of interest corresponding to the exact and approximate solutions of a mathematical model. A high degree of reliability can be achieved in numerical simulation only if both the errors of idealization and errors of discretization can be shown to be small.
The theory of Ritz is applied to the equation that Hamilton called the 'Law of Varying Action.' Direct analytical solutions are obtained for the transient motion of beams, both conservative and nonconservative. The results achieved are compared to exact solutions obtained by the use of rigorously exact free-vibration modes in the differential equations of Lagrange and to an approximate solution obtained through the application of Gurtin's principles for linear elastodynamics. A brief discussion of Hamilton's law and Hamilton's principle is followed by examples of results for both free-free and cantilever beams with various loadings.
The theory of Ritz is applied to the equation that Hamilton called the 'Law of Varying Action'. Direct analytical solutions are obtained for the transient motion of beams, both conservative and nonconservative. The results obtained are compared to exact solutions obtained by the use of rigorously exact free-vibration modes in the differential equations of Lagrange and to an approximate solution obtained through the application of Gurtin's principles for linear elastodynamics. A brief discussion of Hamilton's law and Hamilton's principle is followed by examples of results for both free-free and cantilever beams with various loadings.
Derivation of an exact solution for the problem of line transfer by resonant scattering (assumed coherent in the fluid frame) in a differentially expanding atmosphere. The probabilistic method developed is a generalization of the method of addition of layers plus a decomposition of the radiation field into partial radiation fields due to those photons that scatter exactly n times before escaping (n = 0, 1, 2, 3,
The calculation of the outer inviscid flow about a rectangular wing moving at supersonic speeds is reported. The inviscid equations of motion governing the flow generated by the wing form a set of hyperbolic differential equations. The flow field about the rectangular wing is separated into three regions consisting of the forebody, the afterbody, and the wing wake. Solutions for the forebody are obtained using conical flow techniques while the afterbody and the wing wake regions are treated as initial value problems. The numerical solutions are compared in the two dimensional regions with known exact solutions.
Matched asymptotic expansions (MAE) are used to obtain a first order approximation to the solution of a singularly perturbed second order system. A special case is considered in which the uniform asymptotic solution obtained by MAE is shown to converge to the exact solution. Ways in which the method can be used to sole higher-order linear systems, including those which are not singularly perturbed, are also discussed.
The random element method is a generalization of the random vortex method that was developed for the numerical modeling of momentum transport processes as expressed in terms of the Navier-Stokes equations. The method is based on the concept that random walk, as exemplified by Brownian motion, is the stochastic manifestation of diffusional processes. The algorithm based on this method is grid-free and does not require the diffusion equation to be discritized over a mesh, it is thus devoid of numerical diffusion associated with finite difference methods. Moreover, the algorithm is self-adaptive in space and explicit in time, resulting in an improved numerical resolution of gradients as well as a simple and efficient computational procedure. The method is applied here to an assortment of problems of diffusion of momentum and energy in one-dimension as well as heat conduction in two-dimensions in order to assess its validity and accuracy. The numerical solutions obtained are found to be in good agreement with exact solution except for a statistical error introduced by using a finite number of elements, the error can be reduced by increasing the number of elements or by using ensemble averaging over a number of solutions.
We describe a simple and highly efficient and accurate radiative transfer technique for computing bidirectional reflectance of a macroscopically flat scattering layer composed of nonabsorbing or weakly absorbing, arbitrarily shaped, randomly oriented and randomly distributed particles. The layer is assumed to be homogeneous and optically semi-infinite, and the bidirectional reflection function (BRF) is found by a simple iterative solution of the Ambartsumian's nonlinear integral equation. As an exact Solution of the radiative transfer equation, the reflection function thus obtained fully obeys the fundamental physical laws of energy conservation and reciprocity. Since this technique bypasses the computation of the internal radiation field, it is by far the fastest numerical approach available and can be used as an ideal input for Monte Carlo procedures calculating BRFs of scattering layers with macroscopically rough surfaces. Although the effects of packing density and coherent backscattering are currently neglected, they can also be incorporated. The FORTRAN implementation of the technique is available on the World Wide Web at http://ww,,v.giss.nasa.gov/-crmim/brf.html and can be applied to a wide range of remote sensing, engineering, and biophysical problems. We also examine the potential effect of ice crystal shape on the bidirectional reflectance of flat snow surfaces and the applicability of the Henyey-Greenstein phase function and the 6-Eddington approximation in calculations for soil surfaces.
As more and more nonlinear estimation techniques become available, our interest is in finding out what performance improvement, if any, they can provide for practical nonlinear problems that have been traditionally solved using linear methods. In this paper we examine the problem of estimating spacecraft position using conical scan (conscan) for NASA's Deep Space Network antennas. We show that for additive disturbances on antenna power measurement, the problem can be transformed into a linear one, and we present a general solution to this problem, with the least square solution reported in literature as a special case. We also show that for additive disturbances on antenna position, the problem is a truly nonlinear one, and we present two approximate solutions based on linearization and Unscented Transformation respectively, and one 'exact' solution based on Markov Chain Monte Carlo (MCMC) method. Simulations show that, with the amount of data collected in practice, linear methods perform almost the same as MCMC methods. It is only when we artificially reduce the amount of collected data and increase the level of noise that nonlinear methods show significantly better accuracy than that achieved by linear methods, at the expense of more computation.
This paper presents some exact solutions of problems in random function theory for the purpose of testing the validity of an approximate method known variously in the many different fields of its application as first order smoothing theory, first order cumulant discard, quasilinear theory, or the adiabatic approximation. The hydromagnetic dynamo equations are used here, as particularly appropriate for such an investigation. The calculations show that in one case the exact and approximate solutions agree. In the other case the approximate solution is wrong. Hence, in the absence of a general criterion for validity, a result based on first order smoothing theory is a conjecture rather than a fact. This impacts strongly on much of the recent work on hydromagnetic dynamos.
A relatively simple equation is presented for estimating the induced drag ideal efficiency factor e for arbitrary cross sectional wing forms. This equation is based on eight basic but varied wing configurations which have exact solutions. The e function which relates the basic wings is developed statistically and is a continuous function of configuration geometry. The basic wing configurations include boxwings shaped as a rectangle, ellipse, and diamond; the V-wing; end-plate wing; 90 degree cruciform; circle dumbbell; and biplane. Example applications of the e equations are made to many wing forms such as wings with struts which form partial span rectangle dumbbell wings; bowtie, cruciform, winglet, and fan wings; and multiwings. Derivations are presented in the appendices of exact closed form solutions found of e for the V-wing and 90 degree cruciform wing and for an asymptotic solution for multiwings.