Equilibrium computations for multicomponent plasmas.
Perturbation method analysis of thermodynamic equilibrium properties of multicomponent plasma
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.
Perturbation method analysis of thermodynamic equilibrium properties of multicomponent plasma
An approach is developed to the application of the perturbation method for the solution of problems with essential external nonlinearities, based on identification in the boundary condition of a small nonlinear complex which is considered a perturbing function. The solutions obtained in the first approximation with error of 1 to 2% in calculating the unsteady temperature fields are then used to determine the temperature stresses and deformations in solid bodies of classical form.
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.
Physical phenomena involving rapid and sudden transitions, such as snap buckling of elastic shells, explosions, and earthquakes, are characterized mathematically as a small disturbance causing a large-amplitude response. Because of this, standard asymptotic and perturbation methods are ill-suited to these problems. In the present paper, a new method of analyzing jump phenomena is proposed. The principal feature of the method is the representation of the response in terms of rational functions. For illustration, the method is applied to the snap buckling of an elastic arch and to a simple combustion problem.
Two methods for analyzing the time-scale properties of aircraft trajectory optimization problems are presented. Time-scale properties must be identified before solutions can be obtained by using singular perturbation methods. Both methods only require a knowledge of the state equations, the aircraft characteristics, and the bounds on the state and control variables. Although these methods give only rough estimates of time-scale separation, they do not require that an 'exact' optimal trajectory be known, as do the more rigorous methods, and they are an improvement on the ad hoc methods currently in use. The two methods are applied to an example problem for a high performance aircraft.
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.
The three most successful models for describing scattering from random rough surfaces are the Kirchhoff approximation (KA), the small-perturbation method (SPM), and the two-scale-roughness (or composite roughness) surface-scattering (TSR) models. In this paper it is shown how these three models can be derived rigorously from one perturbation expansion based on the extinction theorem for scalar waves scattering from perfectly rigid surface. It is also shown how corrections to the KA proportional to the surface curvature and higher-order derivatives may be obtained. Using these results, the scattering cross section is derived for various surface models.
The nonlinear evolution of a two-dimensional second mode unstable wave in a Mach 4.5 wall-bounded flow is computed by solving the full time-dependent compressible Navier-Stokes equations. A highly accurate solution is obtained using spectral collocation methods. It is shown that departure from linearity first occurs in the critical layer due to the cubic nonlinearities in the momentum equation. This is a direct result of the large density perturbations in this regime. Time evolution studies of the growth rate as a function of normal distance from the plate suggests that the mode is evolving toward a nonlinear saturated state, and that this problem is possibly amenable to standard weakly nonlinear perturbation methods.
A one-dimensional model for bulk motion induced by a transient volumetric heat source in a confined gas at zero gravity is considered. Rational approximation methods are used to derive a quantitative theory for the gas response to a spatiallly distributed, time-dependent internal power deposition. The resulting low Mach number compressible flow equations are solved by using perturbation methods. Solutions are given for a conduction-free core and thin conductive boundary layers adjacent to the end walls. It is found that may any spatially nonuniform power deposition will cause fluid motion. Net mass transport in the closed container will occur for certain spatially distributed heating. The model mimics the thermal effects of an exothermic gas phase reaction in vapor transport experiments conducted in space. The solutions demonstrate that thermally induced mass transport can be as large as diffusive mass transport in a typical experiment.
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.
The paper describes a perturbation method for turbomachinery calculations, particularly where it is necessary to carry out a number of calculations for closely-related flows such as are needed in a parametric study. The method is applied for solving a model problem involving blade-to-blade surfaces of a transonic compressor. Basically, the method makes use of a previously calculated base solution to determine first-order changes in the flow field due to variations in one or more of a variety of geometrical or flow field parameters. The fundamental assumption associated with the perturbation solution is that the magnitude of the deviations from the base solution lies within the range of a linear perturbation analysis. Comparisons are made for results obtained, by varying the thickness ratio of an unstaggered nonlifting cascade composed of biconvex profiles in a flow with an oncoming freestream Mach number of 0.60.
An accurate method was developed for predicting effects of streamline curvature and coordinate system rotation on turbulent boundary layers. A new two-equation model of turbulence was developed which serves as the basis of the study. In developing the new model, physical reasoning is combined with singular perturbation methods to develop a rational, physically-based set of equations which are, on the one hand, as accurate as mixing-length theory for equilibrium boundary layers and, on the other hand, suitable for computing effects of curvature and rotation. The equations are solved numerically for several boundary layer flows over plane and curved surfaces. For incompressible boundary layers, results of the computations are generally within 10% of corresponding experimental data. Somewhat larger discrepancies are noted for compressible applications.
A numerical solution technique is developed for computing the flow field around an isolated helicopter rotor in hover. The flow is governed by the compressible Euler equations which are integrated using a finite volume approach. The Euler equations are coupled to a free wake model of the rotary wing vortical wake. This wake model is incorporated into the finite volume solver using a prescribed flow, or perturbation, technique which eliminates the numerical diffusion of vorticity due to the artificial viscosity of the scheme. The work is divided into three major parts: (1) comparisons of Euler solutions to experimental data for the flow around isolated wings show good agreement with the surface pressures, but poor agreement with the vortical wake structure; (2) the perturbation method is developed and used to compute the interaction of a streamwise vortex with a semispan wing. The rapid diffusion of the vortex when only the basic Euler solver is used is illustrated, and excellent agreement with experimental section lift coefficients is demonstrated when using the perturbation approach; and (3) the free wake solution technique is described and the coupling of the wake to the Euler solver for an isolated rotor is presented. Comparisons with experimental blade load data for several cases show good agreement, with discrepancies largely attributable to the neglect of viscous effects. The computed wake geometries agree less well with experiment, the primary difference being that too rapid a wake contraction is predicted for all the cases.
A procedure is provided for generalizing the method of harmonic balance to obtain higher-order approximations to the periodic solutions of differential equations for two systems, one a conservative system and the other nonconservative. The procedure is currently being applied to investigate the possible solution behaviors of singular, nonlinear oscillators, where the usual perturbation methods do not work.
Linear equations are solved by a method of superposition of solutions of a sequence of initial value problems. For nonlinear equations and/or boundary conditions, the solution is iterative and in each iteration a problem like the linear case is solved. A simple Taylor series expansion is used for the linearization of both nonlinear equations and nonlinear boundary conditions. The perturbation method of solution is used in preference to quasilinearization because of programming ease, and smaller storage requirements; and experiments indicate that the desired convergence properties exist although no proof or convergence is given.
A procedure cast in a form appropriate to the finite element method is presented for geometrically nonlinear prebuckling and postbuckling structural analysis, including the identification of snap-through type of buckling. The principal features of this procedure are the use of direct iteration for solution of the nonlinear algebraic equations in the prebuckling range, an interpolation scheme for determination of the initial bifurcation point, a perturbation method in definition of the load-displacement behavior through the postbuckling regime, and extrapolation in determination of the limit point for snap-through buckling. Three numerical examples are presented in illustration of the procedure and in comparison with alternative approaches.
The latest trends and theoretical developments involved with the modeling and control of Large Flexible Space Structures (LFSS) are described. The paper addresses first the basic problems, characteristics, and difficulties inherent in modeling and control of LFSS. Major sources of difficulties and errors are the stiffness and damping operators of the dynamic model. Extensions of Linear Quadratic Gaussian (LQG) theory as applied to LFSS are presented, including frequency-shaped cost functionals and perturbation methods. The minimum data/maximum entropy approach which uses a stochastic design model to overcome difficulties found in the LQG-based methods is described. Latest trends in system theory including balanced realization and singular-value analysis are used to determine reduced order controllers and models. Ad hoc methods such as component cost analysis and modal cost analysis are discussed in context with the closed-loop reduction problem of controller order versus performance. The minimum data/maximum entropy approach also addresses controller order versus performance. Those areas of control science and large scale systems that appear to have an important role in understanding and solving LFSS modeling and control are also identified.
The present evaluation of the method of asymptotic expansions (MAE) indicates that the various terms of the common solution of MAE can be generated as polynomials in stretched variables, without actually solving them from the outer solution, as is currently the practice. It is also noted that the common solution of the MAE and the intermediate solution of the singular-perturbation method are the same; these methods therefore yield identical results for a certain class of problems. Two illustrative problems are treated.