Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Perturbative methods”

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.

At least 91 records · Page 5

Scattering of Pruppacher-Pitter raindrops at 30 GHz

Optimum design of modern ground-satellite communication systems requires the knowledge of rain-induced differential attenuation, differential phase shift, and cross polarization factors. Different available analytical techniques for raindrop scattering problems were assessed. These include: (1) geometrical theory of diffraction; (2) method of moment; (3) perturbation method; (4) point matching methods; (5) extended boundary condition method; and (6) global-local finite element method. The advantages and disadvantages of each are listed. The extended boundary condition method, which was determined to yield the most scattering results, is summarized. The scattered fields for Pruppacher-Pitter raindrops with sizes ranging from 0.5 mm to 3.5 mm at 20 C and at 30 GHz for several incidence angles are tabulated.

Yeh, C.↗

The Effect of Stochastic Perturbation of Fuel Distribution on the Criticality of a One Speed Reactor and the Development of Multi-Material Multinomial Line Statistics

The effect of random fuel redistribution on the eigenvalue of a one-speed reactor is investigated. An ensemble of such reactors that are identical to a homogeneous reference critical reactor except for the fissile isotope density distribution is constructed such that it meets a set of well-posed redistribution requirements. The average eigenvalue, , is evaluated when the total fissile loading per ensemble element, or realization, is conserved. The perturbation is proven to increase the reactor criticality on average when it is uniformly distributed. The various causes of the change in reactivity, and their relative effects are identified and ranked. From this, a path towards identifying the causes. and relative effects of reactivity fluctuations for the energy dependent problem is pointed to. The perturbation method of using multinomial distributions for representing the perturbed reactor is developed. This method has some advantages that can be of use in other stochastic problems. Finally, some of the features of this perturbation problem are related to other techniques that have been used for addressing similar problems.

Jahshan, S. N.↗

Application of the Method of Coordinate Perturbation to Unsteady Duct Flow

The method of coordinate perturbation is applied to the unsteady flow of a compressible fluid in ducts of variable cross section. Solutions, in the form of perturbation series, are obtained for unsteady flows in ducts for which the logarithmic derivative of area variation with respect to the space coordinate is a function of the 'smallness' parameter of the perturbation series. This technique is applied to the problem of the interaction of a disturbance and a shock wave in a diffuser flow. It is found that, for a special choice of the function describing the disturbance, the path of the shock wave can be expressed in closed form to first order. The method is then applied to the determination of the flow field behind a shock wave moving on a prescribed path in the x,t-plane. Perturbation series solutions for quite general paths are developed. The perturbation series solutions are compared with the more exact solutions obtained by the application of the method of characteristics. The approximate solutions are shown to be in reasonably accurate agreement with the solutions obtained by the method of characteristics.

Himmel, Seymour C.↗

Control optimization, stabilization and computer algorithms for aircraft applications

Research related to reliable aircraft design is summarized. Topics discussed include systems reliability optimization, failure detection algorithms, analysis of nonlinear filters, design of compensators incorporating time delays, digital compensator design, estimation for systems with echoes, low-order compensator design, descent-phase controller for 4-D navigation, infinite dimensional mathematical programming problems and optimal control problems with constraints, robust compensator design, numerical methods for the Lyapunov equations, and perturbation methods in linear filtering and control.

Source record↗

Improvement of modal scaling factors using mass additive technique

A general investigation into the improvement of modal scaling factors of an experimental modal model using additive technique is discussed. Data base required by the proposed method consists of an experimental modal model (a set of complex eigenvalues and eigenvectors) of the original structure and a corresponding set of complex eigenvalues of the mass-added structure. Three analytical methods,i.e., first order and second order perturbation methods, and local eigenvalue modification technique, are proposed to predict the improved modal scaling factors. Difficulties encountered in scaling closely spaced modes are discussed. Methods to compute the necessary rotational modal vectors at the mass additive points are also proposed to increase the accuracy of the analytical prediction.

Zhang, Qiang↗

Efficient optimization of integrated aerodynamic-structural design

The introduction of composite materials is having a profound effect on the design process. Because these materials permit the designer to tailor material properties to improve structural, aerodynamic and acoustic performance, they require a more integrated multidisciplinary design process. Because of the complexity of the design process numerical optimization methods are required. The present paper is focused on a major difficulty associated with the multidisciplinary design optimization process - its enormous computational cost. We consider two approaches for reducing this computational burden: (1) development of efficient methods for cross-sensitivity calculation using perturbation methods; and (2) the use of approximate numerical optimization procedures. Our efforts are concentrated upon combined aerodynamic-structural optimization. Results are presented for the integrated design of a sailplane wing. The impact of our computational procedures on the computational costs of integrated costs of integrated designs are discussed.

Haftka, R. T.↗

Dynamical Evolution of Meteoroid Streams, Developments Over the Last 30 Years

As soon as reliable methods for observationally determining the heliocentric orbits of meteoroids and hence the mean orbit of a meteoroid stream in the 1950s and 60s, astronomers strived to investigate the evolution of the orbit under the effects of gravitational perturbations from the planets. At first, the limitations in the capabilities of computers, both in terms of speed and memory, placed severe restrictions on what was possible to do. As a consequence, secular perturbation methods, where the perturbations are averaged over one orbit became the norm. The most popular of these is the Halphen- Goryachev method which was used extensively until the early 1980s. The main disadvantage of these methods lies in the fact that close encounter can be missed, however they remain useful for performing very long-term integrations. Direct integration methods determine the effects of the perturbing forces at many points on an orbit. This give a better picture of the orbital evolution of an individual meteoroid, but many meteoroids have to be integrated in order to obtain a realistic picture of the evolution of a meteoroid stream. The notion of generating a family of hypothetical meteoroids to represent a stream and directly integrate the motion of each was probably first used by Williams Murray & Hughes (1979), to investigate the Quadrantids. Because of computing limitations, only 10 test meteoroids were used. Only two years later, Hughes et. al. (1981) had increased the number of particles 20-fold to 200 while after a further year, Fox Williams and Hughes used 500 000 test meteoroids to model the Geminid stream. With such a number of meteoroids it was possible for the first time to produce a realistic cross-section of the stream on the ecliptic. From that point on there has been a continued increase in the number of meteoroids, the length of time over which integration is carried out and the frequency with which results can be plotted so that it is now possible to produce moving images of the stream. As a consequence, over recent years, emphasis has moved to considering stream formation and the role fragmentation plays in this.

Williams, I. P.↗

Perturbation solutions for transonic flow on the blade-to-blade surface of compressor blade rows

A preliminary investigation was conducted to establish the theoretical basis of perturbation techniques with the objective of minimizing computational requirements associated with parametric studies of transonic flows in turbomachines. The theoretical analysis involved the development of perturbation methods for determining first order changes in the flow solution due to variations of one or more geometrical or flow parameters. The formulation is primarily directed toward transonic flows on the blade to blade surface of a single blade row compressor. Two different perturbation approaches were identified and studied. Applications and results of these methods for various perturbations are presented for selected two dimensional transonic cascade flows to illustrate the advantages and disadvantages of each technique. Additionally, it was found that, for flows with shock waves, proper account of shock displacement was crucial.

Stahara, S. S.↗

Comparison of remote sensing measurements with a two-scale polarimetric emission and scattering model for sea surfaces

Recently, it has been observed that the brightness temperatures of sea surfaces correlate with the azimuth angle of the ocean wind vector, indicating that ocean wind direction can not only be retrieved from the microwave backscatter, but also from the brightness temperature measurements. In this paper, comparison of the theoretical result calculated from a two-scale emission and scattering model with the Seasat and SSM/I model functions and aircraft measurements is presented and potential applications of microwave polarimetry are discussed. In our two-scale model, the modified reflectivities of large scale surfaces are calculated by extending the small perturbation method to the second order for small scale perturbation with anisotropic directional spectrum. It was found that the modified reflectivities derived from the second-order scattered field agree excellently well with the results obtained from a Monte Carlo simulation technique which numerically calculates the polarimetric reflectivities of one-dimensional random rough surfaces with a power-law spectrum. Without the second-order correction, the modified reflectivities of the rough surfaces are significantly over-estimated and sign errors are observed in the third and fourth Stokes parameters for thermal emissions. The surface spectrum parameters and two-scale cutoff are selected so that the calculated scattering and emission signatures agree with the reported model functions for Seasat and SSM/I. Subsequently, the polarimetric signatures of sea surfaces are illustrated to indicate the possibility of reducing the number of azimuthal looks required for spaceborne sensors in the remote sensing of ocean wind by using the polarimetric information. Furthermore, it is found that contrary to the dependence of backscattering coefficients on incidence angles, polarimetric brightness temperatures display a stronger wind direction dependence in the near nadir-looking direction than away from nadir. Finally, we discuss the potential of a single-look polarimetric radiometer for wind retrieval.

Yueh, S. H.↗

Determination of the Rate Coefficients of the SO2 plus O plus M yields SO3 plus M Reaction

Rate coefficients of the title reaction R(sub 31) (SO2 +O+M yields SO3 +M) and R(sub 56) (SO2 + HO2 yields SO3 +OH), important in the conversion of S(IV) to S(VI),were obtained at T =970-1150 K and rho (sub ave) = 16.2 micro mol/cubic cm behind reflected shock waves by a perturbation method. Shock-heated H2/ O2/Ar mixtures were perturbed by adding small amounts of SO2 (1%, 2%, and 3%) and the OH temporal profiles were then measured using laser absorption spectroscopy. Reaction rate coefficients were elucidated by matching the characteristic reaction times acquired from the individual experimental absorption profiles via simultaneous optimization of k(sub 31) and k(sub 56) values in the reaction modeling (for satisfactory matches to the observed characteristic times, it was necessary to take into account R(sub 56)). In the experimental conditions of this study, R(sub 31) is in the low-pressure limit. The rate coefficient expressions fitted using the combined data of this study and the previous experimental results are k(sub 31,0)/[Ar] = 2.9 10(exp 35) T(exp ?6.0) exp(?4780 K/T ) + 6.1 10(exp 24) T(exp ?3.0) exp(?1980 K/T ) cm(sup 6) mol(exp ?2)/ s at T = 300-2500 K; k(sub 56) = 1.36 10(exp 11) exp(?3420 K/T ) cm(exp 3)/mol/s at T = 970-1150 K. Computer simulations of typical aircraft engine environments, using the reaction mechanism with the above k(sub 31,0) and k(sub 56) expressions, gave the maximum S(IV) to S(VI) conversion yield of ca. 3.5% and 2.5% for the constant density and constant pressure flow condition, respectively. Moreover, maximum conversions occur at rather higher temperatures (?1200 K) than that where the maximum k(sub 31,0) value is located (approximately 800 K). This is because the conversion yield is dependent upon not only the k(sup 31,0) and k(sup 56) values (production flux) but also the availability of H, O, and HO2 in the system (consumption flux).

Hwang, S. M.↗

The B/K/ method - Application to methylene

The B(K) method is compared with conventional Configuration Interaction (CI) and other perturbation methods for the 3B1 and 1A1 states of methylene. Use of Rayleigh-Schroedinger perturbation theory for construction of an effective Hamiltonian is seen to be more accurate for estimating the full CI energy limit.

Davidson, E. R.↗

Review of probabilistic analysis of dynamic response of systems with random parameters

The various methods that have been studied in the past to allow probabilistic analysis of dynamic response for systems with random parameters are reviewed. Dynamic response may have been obtained deterministically if the variations about the nominal values were small; however, for space structures which require precise pointing, the variations about the nominal values of the structural details and of the environmental conditions are too large to be considered as negligible. These uncertainties are accounted for in terms of probability distributions about their nominal values. The quantities of concern for describing the response of the structure includes displacements, velocities, and the distributions of natural frequencies. The exact statistical characterization of the response would yield joint probability distributions for the response variables. Since the random quantities will appear as coefficients, determining the exact distributions will be difficult at best. Thus, certain approximations will have to be made. A number of techniques that are available are discussed, even in the nonlinear case. The methods that are described were: (1) Liouville's equation; (2) perturbation methods; (3) mean square approximate systems; and (4) nonlinear systems with approximation by linear systems.

Kozin, F.↗

A hybrid Pade-Galerkin technique for differential equations

A three-step hybrid analysis technique, which successively uses the regular perturbation expansion method, the Pade expansion method, and then a Galerkin approximation, is presented and applied to some model boundary value problems. In the first step of the method, the regular perturbation method is used to construct an approximation to the solution in the form of a finite power series in a small parameter epsilon associated with the problem. In the second step of the method, the series approximation obtained in step one is used to construct a Pade approximation in the form of a rational function in the parameter epsilon. In the third step, the various powers of epsilon which appear in the Pade approximation are replaced by new (unknown) parameters (delta(sub j)). These new parameters are determined by requiring that the residual formed by substituting the new approximation into the governing differential equation is orthogonal to each of the perturbation coordinate functions used in step one. The technique is applied to model problems involving ordinary or partial differential equations. In general, the technique appears to provide good approximations to the solution even when the perturbation and Pade approximations fail to do so. The method is discussed and topics for future investigations are indicated.

Geer, James F.↗

Methods for the calculation of axial wave numbers in lined ducts with mean flow

A survey is made of the methods available for the calculation of axial wave numbers in lined ducts. Rectangular and circular ducts with both uniform and non-uniform flow are considered as are ducts with peripherally varying liners. A historical perspective is provided by a discussion of the classical methods for computing attenuation when no mean flow is present. When flow is present these techniques become either impractical or impossible. A number of direct eigenvalue determination schemes which have been used when flow is present are discussed. Methods described are extensions of the classical no-flow technique, perturbation methods based on the no-flow technique, direct integration methods for solution of the eigenvalue equation, an integration-iteration method based on the governing differential equation for acoustic transmission, Galerkin methods, finite difference methods, and finite element methods.

Eversman, W.↗