Engineering Papers⌕ Search

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.

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

At least 1,261 records · Page 70

On-line optimization of aircraft altitude and flight path angle dynamics

The optimization of aircraft altitude and flight path angle dynamics is addressed in a form suitable for on-line computation and control. The approach here is a direct extension of the work reported by Calise (1979), where singular perturbation methods were used to optimize position, energy and heading dynamics; it thus represents an optimal control solution that models all of the primary trajectory related dynamics. It is pointed out that the resulting algorithm can be regarded as a nonlinear feedback control law. The minimum time intercept of a fixed terminal point is used in setting the framework in which the analytical results are developed. The main theoretical result is that the dynamics, while not completely separable, can be approximated by singular perturbation methods when the control model includes relative position dynamics.

Calise, A. J.↗

Numerical techniques in radiative heat transfer for general, scattering, plane-parallel media

The study of radiative heat transfer with scattering usually leads to the solution of singular Fredholm integral equations. The present paper presents an accurate and efficient numerical method to solve certain integral equations that govern radiative equilibrium problems in plane-parallel geometry for both grey and nongrey, anisotropically scattering media. In particular, the nongrey problem is represented by a spectral integral of a system of nonlinear integral equations in space, which has not been solved previously. The numerical technique is constructed to handle this unique nongrey governing equation as well as the difficulties caused by singular kernels. Example problems are solved and the method's accuracy and computational speed are analyzed.

Sharma, A.↗

Wall interference correction improvements for the ONERA main wind tunnels

This paper describes improved methods of calculating wall interference corrections for the ONERA large windtunnels. The mathematical description of the model and its sting support have become more sophisticated. An increasing number of singularities is used until an agreement between theoretical and experimental signatures of the model and sting on the walls of the closed test section is obtained. The singularity decentering effects are calculated when the model reaches large angles of attack. The porosity factor cartography on the perforated walls deduced from the measured signatures now replaces the reference tests previously carried out in larger tunnels. The porosity factors obtained from the blockage terms (signatures at zero lift) and from the lift terms are in good agreement. In each case (model + sting + test section), wall corrections are now determined, before the tests, as a function of the fundamental parameters M, CS, CZ. During the windtunnel tests, the corrections are quickly computed from these functions.

Vaucheret, X.↗

An implicit, transonic, full-potential code for cascade flow on H-grid topology

A transonic, full-potential code is developed for computing the flow through two-dimensional cascades using an H-type grid topology that employs an implicit approximate-factorization scheme. The body-conforming H-grid is generated numerically by solving Poisson's equation. The flow-solution algorithm at the coordinate mapping singularity associated with this grid is investigated using two different types of finite-difference schemes. The grid-geometry effect on these schemes is also studied by noting free-stream capturing properties. It is found that by implementing a consistent spatial differencing scheme, the mapping singularities can be resolved numerically, and the grid-geometry-induced error minimized. The code is verified by computing model cascade flow problems.

Kwak, D.↗

Modified SUMT for structural synthesis

The concepts of the singular perturbation theory are employed to modify the SUMT (Sequential Unconstrained Minimization Technique) algorithm of Fiacco-McCormick (1968) in order to make it more robust and reliable. The algorithm, which uses a sequence of penalty parameters to convert a constrained problem into a sequence of unconstrained problems, suffers from the need to minimize an ill-conditioned penalty function. By using the modified SUMT algorithm on two different structural optimization problems, it is shown that the singular perturbation SUMT easily converges to accurate solutions.

Kamat, M. P.↗

Derivation of the fundamental equation of sound generated by moving aerodynamic surfaces

Ffowcs Williams and Hawkings (1969) based their derivation of the fundamental equation of the sound generated by arbitrarily moving aerodynamic surfaces on the study of mass and momentum balance of a control volume imbedding a mathematical surface(s) exactly corresponding to real surface(s). These investigators also sketched an alternative method, employing generalized functions, for its derivation. This latter method, which was later developed by Farassat (1975), is purely mathematical and formal. Goldstein (1976) used the free-space Green function to produce an implicit derivation of the Ffowcs Williams and Hawkings equation. In the study presented here, Lowson's (1965) concept of moving point singularities is generalized to moving surface singularities, and a new derivation is given of the fundamental equation. The derivation is based on topological considerations of the underlying space, the fluid medium, and the integral properties of the Dirac delta function.

Aggarwal, H. R.↗

Post instability in continuous systems. II - Post-instability models of continua

New models simulating postinstability behavior of continuous systems are discussed. Multivalued velocity fields of the motion of a medium are analyzed, deriving concepts which are applied to a new model of turbulence. Postinstability models of solids are examined, emphasizing the multivaluedness of the parameters. The consequences are addressed of the simplifying assumption that in solids the rate of dissipation of energy of the pulsations is much higher than in fluids. The principal stress formulation of postinstability models is considered, including the geometry, kinematics, and dynamics of the principal stress coordinates. The fundamental properties of enlarged models of a solid are illustrated by the theory of wrinkling films. Results are presented for the statics of such films, the effects of shocks and shock waves, and the behavior of singular effects and singular thick wrinkles.

Zak, M.↗

Energy state revisited

Kaiser (1944) has introduced the concept of 'resultant height' in connection with aircraft minimum-time climbs. Its use as a state variable in trajectory work is attractive because it is a 'slower' variable than either altitude or velocity. Kelley (1972, 1973) has made an attempt to synthesize 'slow' state variables in connection with singular-perturbation procedures. In the present investigation, attempts are made to synthesize both 'fast' and 'slow' variables for the minimum-time-to-climb problem along lines explored by Kelley. Attention is given to climb equations, energy-modeling simplifications, 'slow'-variable choice, 'fast'-variable-choice considerations, a singular-perturbation analysis, the choice of a 'fast' variable, and the climb-dash problem.

Kelley, H. J.↗

A least squares finite element scheme for transonic flow around harmonically oscillating airfoils

The present investigation shows that a finite element scheme with a weighted least squares variational principle is applicable to the problem of transonic flow around a harmonically oscillating airfoil. For the flat plate case, numerical results compare favorably with the exact solution. The obtained numerical results for the transonic problem, for which an exact solution is not known, have the characteristics of known experimental results. It is demonstrated that the performance of the employed numerical method is independent of equation type (elliptic or hyperbolic) and frequency. The weighted least squares principle allows the appropriate modeling of singularities, which such a modeling of singularities is not possible with normal least squares.

Cox, C. L.↗

The mechanics of delamination in fiber-reinforced composite materials. Part 2: Delamination behavior and fracture mechanics parameters

Based on theories of laminate anisotropic elasticity and interlaminar fracture, the complete solution structure associated with a composite delamination is determined. Fracture mechanics parameters characterizing the interlaminar crack behavior are defined from asymptotic stress solutions for delaminations with different crack-tip deformation configurations. A numerical method employing singular finite elements is developed to study delaminations in fiber composites with any arbitrary combinations of lamination, material, geometric, and crack variables. The special finite elements include the exact delamination stress singularity in its formulation. The method is shown to be computationally accurate and efficient, and operationally simple. To illustrate the basic nature of composite delamination, solutions are shown for edge-delaminated (0/-0/-0/0) and (+ or - 0/+ or - 0/90/90 deg) graphite-epoxy systems under uniform axial extenstion. Three-dimensional crack-tip stress intensity factors, associated energy release rates, and delamination crack-closure are determined for each individual case. The basic mechanics and mechanisms of composite delamination are studied, and fundamental characteristics unique to recently proposed tests for interlaminar fracture toughness of fiber composite laminates are examined.

Wang, S. S.↗

Energy state revisited

Kaiser (1944) has introduced the concept of "resultant height' in connection with aircraft minimum-time climbs. Its use as a state variable in trajectory work is attractive because it is a "slower' variable than either altitude or velocity. Kelley (1972, 1973) has made an attempt to synthesize "slow' state variables in connection with singular-perturbation procedures. In the present investigation, attempts are made to synthesize both "fast' and "slow' variables for the minimum-time-climb problem along lines explored by Kelley. Attention is given to climb equations, energy-modeling simplifications, "slow'-variable choice, "fast'-variable-choice considerations, a singular-perturbation analysis, the choise of a "fast' variable, and the climb-dash problem.

Kelley, H. J.↗

Kinematics, influence functions and field quantities for disturbance propagation from moving disturbance sources

A unified method is presented for deriving the influence functions of moving singularities which determine the field quantities in aerodynamics and aeroacoustics. The moving singularities comprise volume and surface distributions having arbitrary orientations in space and to the trajectory. Hence one generally valid formula for the influence functions which reveal some universal relationships and remarkable properties in the disturbance fields. The derivations used are completely consistent with the physical processes in the propagation field, such that treatment renders new descriptions for some standard concepts. The treatment is uniformly valid for subsonic and supersonic Mach numbers.

Das, A.↗

Crack problems for bonded nonhomogeneous materials under antiplane shear loading

The singular nature of the crack tip stress field in a nonhomogeneous medium with a shear modulus with a discontinuous derivative was investigated. The simplest possible loading and geometry, the antiplane shear loading of two bonded half spaces in which the crack is perpendicular to the interface is considered. It is shown that the square root singularity of the crack tip stress field is unaffected by the discontinuity in the derivative of the shear modulus. The problem is solved for a finite crack and results for the stress intensity factors are presented.

Erdogan, F.↗

The mechanics of delamination in fiber-reinforced composite materials. II - The delamination behavior and fracture mechanics parameters

Based on theories of laminate anisotropic elasticity and interlaminar fracture, the complete solution structure associated with a composite delamination is determined. Fracture mechanics parameters characterizing the interlaminar crack behavior are defined from asymptotic stress solutions for delaminations with different crack-tip deformation configurations. A numerical method employing singular finite elements is developed to study delaminations in fiber composites with any arbitrary combinations of lamination, material, geometric, and crack variables. The special finite elements include the exact delamination stress singularity in its formulation. The method is shown to be computationally accurate and efficient, and operationally simple. To illustrate the basic nature of composite delamination, solutions are shown for edge-delaminated (0/-0/-0/0) and (+ or - 0/+ or - 0/90/90 deg) graphite-epoxy systems under uniform axial extension. Three-dimensional crack-tip stress intensity factors, associated energy release rates, and delamination crack-closure are determined for each individual case. The basic mechanics and mechanisms of composite delamination are studied, and fundamental characteristics unique to recently proposed tests for interlaminar fracture toughness of fiber composite laminates are examined. Previously announced in STAR as N84-13222

Wang, S. S.↗

Thermal ripples in a resistive and radiative instability

The development of the resistive tearing instability in the case of sheared magnetic fields is considered, taking into account also the occurrence of a radiatively driven thermal instability. It is pointed out that thermal conduction has generally been neglected in theories similar to those discussed. The present investigation is concerned with a consideration of both parallel and perpendicular thermal conduction, in addition to finite resistivity and radiative loss. Attention is given to the equations and the model, the spatial singularity which arises with consideration of only the parallel heat conduction, the removal of this singularity and the formation of temperature oscillations (thermal ripples) by inclusion of the perpendicular heat-flux component, and details regarding the numerical procedure. A brief explanation is provided of the conditions required for the oscillations, and potential implications of the results with respect to the solar flare are discussed.

Steinolfson, R. S.↗

Progress in wind tunnel wall interference assessment/correction procedures at the NAE

Wall corrections as a function of wall porosity in the transonic wall interference problem was assessed. Effective porosities primarily for the two dimensional case were established as follows: (1) comparison of experimental data for two geometrically similar models of different chord/height ratio, an overall value of wall porosity could be deduced; (2) theoretical development which allows for unequal porosity for the floor and ceiling and wall boundary pressure measurements, porosities for floor and ceiling could be deduced; (3) a scheme was developed which allowed unequal porosity of floor and ceiling and streamwise varying porosity. The boundary layer development along the perforated floor and ceiling under the influence of the model pressure field, variations in boundary layer thickness underlining the difficulties in deducing meaningful values of wall porosity were determined. Wall boundary pressure measurement, in combination with singularity modelling of the airfoil, was sufficient to yield required information on the wall interference flow without having to establish some value for wall porosity. The singularity modelling of the airfoil initially covered only lift and volume but was extended to include drag and pitching moment, and second order volume term. It is shown by asymptotic transonic small disturbance analysis, that the derived corrections to angle of attack and free stream Mach number are correct to the first order.

Ohman, L. H.↗

Comparison of uniform perturbation solutions and numerical solutions for some potential flows past slender bodies

Approximate solutions for potential flow past an axisymmetric slender body and past a thin airfoil are calculated using a uniform perturbation method and then compared with either the exact analytical solution or the solution obtained using a purely numerical method. The perturbation method is based upon a representation of the disturbance flow as the superposition of singularities distributed entirely within the body, while the numerical (panel) method is based upon a distribution of singularities on the surface of the body. It is found that the perturbation method provides very good results for small values of the slenderness ratio and for small angles of attack. Moreover, for comparable accuracy, the perturbation method is simpler to implement, requires less computer memory, and generally uses less computation time than the panel method. In particular, the uniform perturbation method yields good resolution near the regions of the leading and trailing edges where other methods fail or require special attention.

Wong, T. C.↗

PROGRAM VSAERO: A computer program for calculating the non-linear aerodynamic characteristics of arbitrary configurations: User's manual

VSAERO is a computer program used to predict the nonlinear aerodynamic characteristics of arbitrary three-dimensional configurations in subsonic flow. Nonlinear effects of vortex separation and vortex surface interaction are treated in an iterative wake-shape calculation procedure, while the effects of viscosity are treated in an iterative loop coupling potential-flow and integral boundary-layer calculations. The program employs a surface singularity panel method using quadrilateral panels on which doublet and source singularities are distributed in a piecewise constant form. This user's manual provides a brief overview of the mathematical model, instructions for configuration modeling and a description of the input and output data. A listing of a sample case is included.

Maskew, B.↗