SEARCH · Engineering Papers

Results for “Singularity”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,117 records · Page 62

All sources NASA NTRS DOE OSTI

NASA NTRS1983

Unsteady laminar boundary-layer calculations on oscillating configurations including backflow. Part 1: Flat plate, oscillating in its own plane

A finite difference method has been developed to calculate the unsteady boundary layer over an oscillating flat plate. Low- and high frequency approximations were used for comparison with numerical results. Special emphasis was placed on the behavior of the flow and on the numerical calculation procedure as soon as reversed flow has occurred over part of the oscillation cycle. The numerical method displayed neither problems nor singular behavior at the beginning of or within the reversed flow region. Calculations, however, came to a limit where the back-flow region reached the plate's leading edge in the case of high oscillation amplitudes. It is assumed that this limit is caused by the special behavior of the flow at the plate's leading edge where the boundary layer equations are not valid.

Geissler, W.↗

NASA NTRS1983

A class of nonideal solutions. 1: Definition and properties

A class of nonideal solutions is defined by constructing a function to represent the composition dependence of thermodynamic properties for members of the class, and some properties of these solutions are studied. The constructed function has several useful features: (1) its parameters occur linearly; (2) it contains a logarithmic singularity in the dilute solution region and contains ideal solutions and regular solutions as special cases; and (3) it is applicable to N-ary systems and reduces to M-ary systems (M or = N) in a form-invariant manner.

Zeleznik, F. J.↗

NASA NTRS1983

Numerical methods for stiff systems of two-point boundary value problems

Numerical procedures are developed for constructing asymptotic solutions of certain nonlinear singularly perturbed vector two-point boundary value problems having boundary layers at one or both endpoints. The asymptotic approximations are generated numerically and can either be used as is or to furnish a general purpose two-point boundary value code with an initial approximation and the nonuniform computational mesh needed for such problems. The procedures are applied to a model problem that has multiple solutions and to problems describing the deformation of thin nonlinear elastic beam that is resting on an elastic foundation.

Flaherty, J. E.↗

NASA NTRS1983

An investigation of the accuracy of finite difference methods in the solution of linear elasticity problems

The accuracy of the finite difference method in the solution of linear elasticity problems that involve either a stress discontinuity or a stress singularity is considered. Solutions to three elasticity problems are discussed in detail: a semi-infinite plane subjected to a uniform load over a portion of its boundary; a bimetallic plate under uniform tensile stress; and a long, midplane symmetric, fiber reinforced laminate subjected to uniform axial strain. Finite difference solutions to the three problems are compared with finite element solutions to corresponding problems. For the first problem a comparison with the exact solution is also made. The finite difference formulations for the three problems are based on second order finite difference formulas that provide for variable spacings in two perpendicular directions. Forward and backward difference formulas are used near boundaries where their use eliminates the need for fictitious grid points.

Bauld, N. R., Jr.↗

NASA NTRS1983

Analysis of structural perturbations in systems via cost decomposition methods

It has long been common practice to analyze linear dynamic systems by decomposing the total response in terms of individual contributions which are easier to analyze. Examples of this philosophy include the expansion of transfer functions using: (1) the superposition principle, (2) residue theory and partial fraction expansions, (3) Markov parameters, Hankel matrices, and (4) regular and singular perturbations. This paper summarizes a new and different kind of expansion designed to decompose the norm of the response vector rather than the response vector itself. This is referred to as "cost-decomposition' of the system. The notable advantages of this type of decomposition are: (a) easy application to multi-input, multi-output systems, (b) natural compatibility with Linear Quadratic Gaussian Theory, (c) applicability to the analysis of more general types of structural perturbations involving inputs, outputs, states, parameters. Property (c) makes the method suitable for problems in model reduction, measurement/actuator selections, and sensitivity analysis.

Skelton, R. E.↗

NASA NTRS1982

Pole placement and order reduction in two-time-scale control systems through Riccati iteration

A transformation of variables taken from singular perturbations may be applied to two-time-scale linear systems in state space form to reduce the system to block-diagonal form with slow and fast modes decoupled. The transformation is easily computed by applying the new Riccati iteration. The iteration yields a solution to the nonsymmetric algebraic Riccati equation obtained by partitioning the original system matrix A. The numerical procedure is initiated with the trivial iterate L(0) = 0, and is globally convergent to the desired unique time scale decoupling solution. After transformation, the decoupled system may be used in controller design to achieve exact closed loop pole placement in the slow subsystem without altering the poles of the fast subsystem. The decoupled form may also be used to reduce system order by wetting a small parameter to zero. Provided the fast subsystem is stable, the order reduction can be expected to yield a good approximation to the original system. These methods are demonstrated using the 16th order linear model of a turbofan engine.

Anderson, L. R.↗

NASA NTRS1983

A new consistent spatial differencing scheme for the transonic full-potential equation

A new spatial differencing scheme for the transonic full-potential equation in conservative form has been developed. This scheme guarantees zero truncation error on any curvilinear mesh for freestream flows in either two- or three-space dimensions. Solutions obtained with this new differencing scheme, away from freestream regions, exhibit greatly improved accuracy, especially for nonsmooth or singular meshes.

Flores, J.↗

NASA NTRS1983

Flow in a differentially rotated cylindrical drop at low Reynolds number

A liquid drop held captive between parallel disks that are differentially rotated is a model for the swirling flows induced by crystal rotation in the floating-zone process for growing semiconductor materials. An asymptotic analysis for a cylindrical drop is presented that elucidates the structure of the axisymmetric cellular motions caused by disk rotation at low Reynolds number. Variations of meniscus shape induced by these flows are described in the limit of small capillary number. Most cellular flow fields break the bifurcation point that corresponds to the Plateau-Rayleigh limit for the length of a static drop into two disjoint shape families and lower the maximum stable drop length. This effect is studied by a singular bifurcation analysis.

Harriott, G. M.↗

NASA NTRS1983

An analysis of the two-dimensional acoustic field in a nonuniform duct carrying compressible flow

An analytical/numerical linear acoustic solution in a nonuniform two-dimensional duct carrying a compressible mean flow is developed. A quasi-one dimensional mean flow model is employed together with a consistent expression for the cross-flow velocity. The acoustic solution is obtained using the Wave Envelope Method. Numerical results are compared with those of an existing Wave Envelope solution which includes a more general mean flow model. In addition, the singular behavior of the solutions derived in the current work for nearly sonic flow conditions is examined explicitly.

Uenishi, K.↗

NASA NTRS1982

Altitude transitions in energy climbs

The aircraft energy-climb trajectory for configurations with a sharp transonic drag rise is well known to possess two branches in the altitude/Mach-number plane. Transition in altitude between the two branches occurs instantaneously, a 'corner' in the minimum-time solution obtained with the energy-state model. If the initial and final values of altitude do not lie on the energy-climb trajectory, then additional jumps (crude approximations to dives and zooms) are required at the initial and terminal points. With a singular-perturbation approach, a 'boundary-layer' correction is obtained for each altitude jump, the transonic jump being a so-called 'internal' boundary layer, different in character from the initial and terminal layers. The determination of this internal boundary layer is examined and some computational results for an example presented.

Weston, A. R.↗

NASA NTRS1983

Measurements of the scattering of sound by a line vortex

This paper presents measurements of the phase and magnitude of the scattered field arising from the incidence of a monochromatic plane sound field as a steady vortex. The amplitude of the scattered field was found to vary linearly with the vortex strength, and with the incident wave amplitude and frequency as predicted by solutions based on the Born approximation. The scattered field was observed to be nonsingular in the incidence direction, and this was similar to predictions by the Parabolic Equation Method (PEM) rather than the Born approximation, which predicts singular behavior in the incidence direction.

Horne, W. C.↗

NASA NTRS1983

Incidence angle effects on convected gust airfoil noise

An analysis is developed which predicts the influence of airfoil mean loading on noise generation due to convected gusts. The theory is based on a linearization of the exact inviscid equations about a nonuniform compressible mean flow and the solution is developed using singular perturbation techniques. The case of a flat plate airfoil, at incidence angle alpha, interacting with three-dimensional disturbances is analyzed. It is found that in the vicinity of the airfoil leading and trailing edges, local regions are present which scale on the disturbance wavelength, with the noise generation concentrated in these regions. Away from the airfoil edges, the mean flow variation is found to be slow compared to the disturbance wavelength and no significant noise generation occurs. The mean flow variation near the leading edge generates additional noise by distorting the convected gust. The cumulative effect of the airfoil mean loading in the trailing edge region produces a 0(1) phase shift between the disturbances on the upper and lower surfaces of the airfoil. A corresponding 0(1) decrease, compared to the alpha = 0 case, is found in the noise generated at the trailing edge.

Kerschen, E. J.↗

NASA NTRS1983

Compact finite difference schemes for the Euler and Navier-Stokes equations

Implicit compact finite difference schemes for the Euler equations are described which furnish equivalent treatment of the conservation and nonconservation forms; a simple modification yields an entropy-producing scheme. An extension of the scheme also treats the compressible Navier-Stokes equations; when the viscosity and heat conduction coefficients are negligible only the boundary data appropriate to the Euler equation influence the solution to any significant extent, a result consistent with singular perturbation theory.

Rose, M. E.↗

NASA NTRS1983

Transformations of the perturbed two-body problem to unperturbed harmonic oscillators

Singular, nonlinear, and Liapunov unstable equations are made regular and linear through transformations that change the perturbed planar problem of two bodies into unperturbed and undamped harmonic oscillators with constant coefficients, so that the stable solution may be immediately written in terms of the new variables. The use of arbitrary and special functions for the transformations allows the systematic discussion of previously introduced and novel anomalies. For the case of the unperturbed two-body problem, it is proved that if transformations are power functions of the radial variable, only the eccentric and the true anomalies (with the corresponding transformations of the radial variable) will result in harmonic oscillators. The present method significantly reduces computation requirements in autonomous space operations.

Szebehely, V.↗

NASA NTRS1983

A study of H-alpha velocities in NGC 1499, NGC 7000, and IC 1318B/C

Multiple slit echelle spectrograph observations of the H-alpha emission line are used to map the radial velocities of the California Nebula (NGC 1499), the North American Nebula complex (NGC 7000 and IC 5070), and IC 1318B/C. The California Nebula is singularly constant in velocity, considering its geometry. The North American Nebula complex reflects a very simple, classical dynamical picture. The expansion discovered earlier in IC 1318B/C is confirmed, detailed, and the model refined. The new data, along with that in earlier papers of this series, show that stellar wind acceleration and champagne flow mechanisms both play important roles in determining the evolution of H II regions.

Fountain, W. F.↗

NASA NTRS1983

Accuracy and resolution of Earth Radiation Budget Estimates

A numerical filter inversion technique that reduces wide-angle satellite measurements to top-of-the-atmosphere radiant exitances has been proposed for the Earth Radiation Budget Experiment (ERBE). The matrix formulation of this technique is presented, and the design of the numerical filter is discussed. The filter is smoothed with a singular value decomposition. The inversion process is simulated by generating synthetic measurements from a 24 degree spherical harmonic radiation field derived from Nimbus 6 ERB data. The numerical filter is applied to these measurements after they are corrupted with instrument error. The results are curves of expected error versus resolution area.

Green, R. N.↗

NASA NTRS1983

Asymptotic unbounded root loci - Formulas and computation

A new geometric way of computing the asymptotic behavior of unbounded root loci of a strictly proper linear time-invariant control system as loop gain goes to infinity is presented. Properties of certain restricted linear maps and nested restrictions of linear maps are developed, and formulas are obtained for the leading coefficient of the asymptotic values of the unbounded multivariable root loci are obtained in terms of eigenvalues of those maps. Published results and a certain simple null structure assumption are used to relate these asymptotic values to the structure at infinity of the Smith-McMillan form of the open loop transfer function. Explicit matrix formulas for the more abstract derived formulas are given and additional geometric insights are developed with orthogonal projections and singular value decomposition. Formulas for the pivots of the unbounded root loci are calculated and shown to have the same form as the coefficients of the unbounded asymptotic root loci.

Sastry, S. S.↗

NASA NTRS1983

Dynamic fields near a crack tip growing in an elastic-perfectly-plastic solid

A full asymptotic solution is presented for the fields in the neighborhood of the tip of a steadily advancing crack in an incompressible elastic-perfectly-plastic solid. There are four findings for mode I crack growth in the plane strain condition. The first is that the entire crack tip in steady crack growth is surrounded by a plastic region and that no elastic unloading is predicted by the complete dynamic asymptotic solution. The second is that, in contrast to the quasi-static solution, the dynamic solution yields strain fields with a logarithmic singularity everywhere near the crack tip. The third is that whereas the stress field varies throughout the entire crack tip neighborhood, it does not exhibit behavior that can be approximated by a constant field followed by an essentially centered-fan field and then by another constant field, especially for small crack growth speeds. The fourth finding is that there are two shock fronts emanating from the crack tip across which certain stress and strain components undergo jump discontinuities. After reviewing the mode III steady-state crack growth, it is concluded that ductile fracture criteria for nonstationary cracks must be based on solutions that include the inertia effects and that for this purpose quasi-static solutions may be inadequate.

Nemat-Nasser, S.↗