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,297 records · Page 72

All sources NASA NTRS DOE OSTI

NASA NTRS1987

A new formulation of electromagnetic wave scattering using an on-surface radiation boundary condition approach

A new formulation of electromagnetic wave scattering by convex, two-dimensional conducting bodies is reported. This formulation, called the on-surface radiation condition (OSRC) approach, is based upon an expansion of the radiation condition applied directly on the surface of a scatterer. It is now shown that application of a suitable radiation condition directly on the surface of a convex conducting scatterer can lead to substantial simplification of the frequency-domain integral equation for the scattered field, which is reduced to just a line integral. For the transverse magnetic case, the integrand is known explicitly. For the transverse electric case, the integrand can be easily constructed by solving an ordinary differential equation around the scatterer surface contour. Examples are provided which show that OSRC yields computed near and far fields which approach the exact results for canonical shapes such as the circular cylinder, square cylinder, and strip. Electrical sizes for the examples are ka = 5 and ka = 10. The new OSRC formulation of scattering may present a useful alternative to present integral equation and uniform high-frequency approaches for convex cylinders larger than ka = 1. Structures with edges or corners can also be analyzed, although more work is needed to incorporate the physics of singular currents at these discontinuities. Convex dielectric structures can also be treated using OSRC.

Kriegsmann, Gregory A.↗

NASA NTRS1987

Asymptotic expansions in nonlinear rotordynamics

This paper is an examination of special nonlinearities of the Jeffcott equations in rotordynamics. The immediate application of this analysis is directed toward understanding the excessive vibrations recorded in the LOX pump of the SSME during hot-firing ground testing. Deadband, side force, and rubbing are three possible sources of inducing nonlinearity in the Jeffcott equations. The present analysis initially reduces these problems to the same mathematical description. A special frequency, named the nonlinear natural frequency, is defined and used to develop the solutions of the nonlinear Jeffcott equations as singular asymptotic expansions. This nonlinear natural frequency, which is the ratio of the cross-stiffness and the damping, plays a major role in determining response frequencies.

Day, William B.↗

NASA NTRS1987

Experience in distributed parameter modeling of the Spacecraft Control Laboratory Experiment (SCOLE) structure

The Spacecraft Control Laboratory Experiment (SCOLE) configuration is used to compare exact and approximate solutions of the partial differential equations which define its structural dynamics. The need for a proof model for evaluating competing control laws demands that solutions be generated which not only exhibit accurate modal characteristics, but precise static deflections as well. Because precise pointing is required, the motion of the end bodies of the Shuttle-attached antenna must be known with great accuracy. Modal models are attractive because of their stable solutions but require hundreds of modes to obtain a static deflection accuracy of only one percent. Although proportional damping in bending agrees well with experimental results using the SCOLE experimental apparatus, modes which involve both torsion and bending differ significantly from proportional damping. A lumped mass model is used to generate exact static deflections, but only approximate modal characteristics. Asymptotic solutions to the distributed parameter system approximate very accurately the modal characteristics at high mode numbers. Ways are examined for refining the approximate solutions by applying a first-order variation and by employing singular perturbation techniques which are usually limited to ordinary differential equations. The most accurate solutions of the distributed parameter model of SCOLE are obtained by combining exact and asymptotic solutions.

Taylor, L. W.↗

NASA NTRS1987

Evaluation of constraint stabilization procedures for multibody dynamical systems

Comparative numerical studies of four constraint treatment techniques for the simulation of general multibody dynamic systems are presented, and results are presented for the example of a classical crank mechanism and for a simplified version of the seven-link manipulator deployment problem. The staggered stabilization technique (Park, 1986) is found to yield improved accuracy and robustness over Baumgarte's (1972) technique, the singular decomposition technique (Walton and Steeves, 1969), and the penalty technique (Lotstedt, 1979). Furthermore, the staggered stabilization technique offers software modularity, and the only data each solution module needs to exchange with the other is a set of vectors plus a common module to generate the gradient matrix of the constraints, B.

Park, K. C.↗

NASA NTRS1987

Quaternions as astrometric plate constants

A new method for solving problems in relative astrometry is proposed. In it, the relationship between the measured quantities and the components of the position vector of a star is modeled using quaternions, in effect replacing the plate constants of a standard four-plate-constant solution with the four components of a quaternion. The method allows a direct solution for the position vectors of the stars, and hence for the equatorial coordinates. Distortions, magnitude, and color effects are readily incorporated into the formalism, and the method is directly applicable to overlapping-plate problems. The advantages of the method include the simplicity of the resulting equations, their freedom from singularities, and the fact that trigonometric functions and tangential point transformations are not needed to model the plate material. A global solution over the entire sky is possible.

Jefferys, William H.↗

NASA NTRS1986

Approximate methods for black hole collisions

The head-on collision of two nonrotating black holes, initially falling toward each other with speed of approach equal to escape velocity, is studied. The initial-value problem for the two-black-hole configuration is solved by giving the extrinsic curvature on a Euclidean space slice, with the equal masses of the two holes and their separation as free parameters. An approximating space-time is then constructed by stacking up a whole sequence of initial-value slices with decreasing separation between the black holes. The space-time singularities at the centers of the two black holes collide at some instant of time, and after that time the space-time so constructed is just the Schwarzschild solution with mass equal to the total mass of both black holes. The approximating space-time describes the coalescence of the event horizons, but contains no gravitational waves.

Eardley, D. M.↗

NASA NTRS1987

The synoptic setting and possible energy sources for mesoscale wave disturbances

Published data on 13 cases of mesoscale wave disturbances and their environment were examined to isolate common features for these cases and to determine possible energy sources for the waves. These events are characterized by either a singular wave of depression or wave packets with periods of 1-4 h, horizontal wavelengths of 50-500 km, and surface-pressure perturbation amplitudes of 0.2-7.0 mb. These wave events are shown to be associated with a distinct synoptic pattern (including the existence of a strong inversion in the lower troposphere and the propagation of a jet streak toward a ridge axis in the upper troposphere) while displaying little correlation with the presence of convective storm cells. The observed development of the waves is consistent with the hypothesis that the energy source needed to initiate and sustain the wave disturbances may be related to a geostrophic adjustment process associated with upper-tropospheric jet streaks.

Uccellini, Louis W.↗

NASA NTRS1987

A new algorithm for the Navier-Stokes equations applied to transonic flows over wings

A new algorithm has been developed for the Navier-Stokes equations. For the convective terms, this algorithm employs flux vector splitting together with a locally rotated coordinate system that is aligned with the streamwise direction. The algorithm has been implemented into a Navier-Stokes code, the NASA Ames Research Center's TNS code, and several cases of steady flow have been calculated. The results show improvements over the original method in the code. First, in a case of separated flow, there is better agreement with the experimental results and also there is an absence of numerical oscillations that occur with the original method. Second, in a case with a strong supersonic to subsonic shock wave, comparisons show an improvement in the resolution of the shock wave and the reexpansion singularity.

Goorjian, Peter M.↗

NASA NTRS1987

Interaction of a decaying vortex ring with a rotational background flow bounded by a solid wall

The interaction between a vortex ring of finite strength and an axisymmetric rotational background flow bounded by a solid wall is studied by a singular perturbation method. The analysis is carried out by combining a composite solution of a vortex ring and an unsteady Navier-Stokes solution for the background rotational flow. Using the method of averaging, numerical scheme is developed to obtain a Navier-Stokes solution in which the grid and time-step sizes depend solely on the length and velocity scales of the background flow. Numerical results are presented to illustrate the separation of the boundary layer on a solid wall and its interaction with the vortex ring.

Ishii, K.↗

NASA NTRS1987

Spatial correlations in multifractals

Spatial correlations within multifractals or fractal measures are considered. It is shown how scaling dimensions can be used to predict scaling laws for spatial correlation functions. It is concluded that the overlaid singularity picture does not provide an adequate description of the spatial properties of multifractals, and that blob concepts seem to provide a much more suitable description of multifractal correlations.

Cates, M. E.↗

NASA NTRS1987

Gravitational lensing by an ensemble of isothermal galaxies

Calculation of 28,000 models of gravitational lensing of a distant quasar by an ensemble of randomly placed galaxies, each having a singular isothermal mass distribuiton, is reported. The average surface mass density was 0.2 of the critical value in all models. It is found that the surface mass density averaged over the area of the smallest circle that encompasses the multiple images is 0.82, only slightly smaller than expected from a simple analytical model of Turner et al. (1984). The probability of getting multiple images is also as large as expected analytically. Gravitational lensing is dominated by the matter in the beam; i.e., by the beam convergence. The cases where the multiple imaging is due to asymmetry in mass distribution (i.e., due to shear) are very rare. Therefore, the observed gravitational-lens candidates for which no lensing object has been detected between the images cannot be a result of asymmetric mass distribution outside the images, at least in a model with randomly distributed galaxies. A surprisingly large number of large separations between the multiple images is found: up to 25 percent of multiple images have their angular separation 2 to 4 times larger than expected in a simple analytical model.

Katz, Neal↗

NASA NTRS1987

The thermal stability of coronal loops by nonlinear diffusion asymptotics

A nonlinear reaction-diffusion equation and some additional constraints are derived which describe the time-dependent behavior of the temperature structure of the plasma in coronal loops. The equation is analyzed using nonlinear diffusion asymptotics, in particular singular perturbation techniques, and the results are interpreted in the context of the physical problem of the thermal stability and temporal behavior of the plasma. The results are consistent with the possibility of cyclic thermal behavior of the plasma, as suggested by Kuin and Martens (1982).

Pakkert, J. W.↗

NASA NTRS1987

Orbital changes during hypersonic aerocruise

A novel mathematical approach that allows the analysis of orbital changes occurring during an aerocruise maneuver to be conducted in two distinct stages is presented. In the first stage, the aerodynamic turn is determined using a nondimensional form of the equations of motion that is free of singularities, and the way in which speed, altitude, angle of attack, and thrust direction should be chosen to maximize the aerodynamic turn for a given propellant expenditure is demonstrated. In the second analysis stage, the aerodynamic turn is translated into changes in the orbital elements with respect to the equatorial plane; analytic solutions for the initial arguments of latitude that maximize the change in inclination and in the longitude of the ascending node are given. As the initial inclination decreases toward zero, the optimal location moves from the apex toward the node.

Mease, Kenneth D.↗

NASA NTRS1987

Optimal heading change with minimum energy loss for a hypersonic gliding vehicle

A three state model is presented for analyzing the problem of optimal changes in heading with minimum energy loss for a hypersonic gliding vehicle. A further model order reduction to a single state model is examined using singular perturbation theory. The optimal solution for the reduced problem defines an optimal altitude profile dependent on the current energy of the vehicle, and the corresponding optimal lift and bank angle. A separate boundary layer analysis, based on an expansion of the necessary conditions about the reduced solution, is used to account for altitude and flight path angle dynamics and to derive a guidance law in feedback form. The guidance law is evaluated for a hypothetical vehicle.

Calise, Anthony J.↗

NASA NTRS1987

Robust eigensystem assignment for flexible structures

An improved method is developed for eigenvalues and eigenvectors placement of a closed-loop control system using either state or output feedback. The method basically consists of three steps. First, the singular value of QR decomposition is used to generate an orthonormal basis that spans admissible eigenvector space corresponding to each assigned eigenvalue. Secondly, given a unitary matrix, the eigenvector set which best approximates the given matrix in the least-square sense and still satisfy eigenvalue cosntraints is determined. Thirdly, a unitary matrix is sought to minimize the error between the unitary matrix and the assignable eigenvector matrix. For use as the desired eigenvector set, two matrices, namely, the open-loop eigenvector matrix and its closest unitary matrix are proposed. The latter matrix generally encourages both minimum conditioning and control gains. In addition, the algorithm is formulated in real arithmetic for efficient implementation. To illustrate the basic concepts, numerical examples are included.

Juang, Jer-Nan↗

NASA NTRS1986

Multigrid-sinc methods

A Galerkin method using Whittaker cardinal or 'sinc' functions as basis functions is described for the solution of boundary-value problems. When the solution is analytic in the interior of the domain, the error of approximation using 2N + 1 points is O(e exp /-gamma sq rt N/) even if derivatives of the solution are singular at the boundaries. A multigrid method with overall complexity O(N log N) is used to solve the discrete equations. This paper contains a description of the multigrid-sinc algorithm along with some preliminary numerical results for two-point boundary-value problems.

Schaffer, Steve↗

NASA NTRS1987

Three dimensional numerical scattering from ocean-like surfaces

A new method of calculating the electric field scattered from two-dimensional conducting rough surfaces is presented. This method uses the operator conjugate gradients technique and the extended boundary condition equation. This avoids the build-up of round-off errors encountered in the inversion of large matrices. It also avoids the singularities of the integral kernels usually encountered in the method of moments. Techniques for improving the rate of convergence by applying physical constraints are discussed. This method is presently being used to calculate the polarization signatures from simulated ocean-like surfaces.

Rodriguez, Ernesto↗

NASA NTRS1984

Cellular interface morphologies in directional solidification. II - The effect of grain boundaries

A singular perturbation analysis valid for small grain-boundary slopes is used with the one-sided model for solidification to show that grain boundaries introduce imperfections into the symmetry of the developing cellular interfaces which rupture the junction between the family of planar shapes and the bifurcating cellular families. Undulating interfaces are shown to develop first near grain boundaries, and to evolve with decreasing temperature gradient either by a smooth transition from the almost planar family or by a sudden jump to moderate-amplitude cellular forms, depending on the growth rate.

Ungar, Lyle H.↗