Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “singular equations”

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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
NASA NTRS1979

An evaluation of four single element airfoil analytic methods

A comparison of four computer codes for the analysis of two-dimensional single element airfoil sections is presented for three classes of section geometries. Two of the computer codes utilize vortex singularities methods to obtain the potential flow solution. The other two codes solve the full inviscid potential flow equation using finite differencing techniques, allowing results to be obtained for transonic flow about an airfoil including weak shocks. Each program incorporates boundary layer routines for computing the boundary layer displacement thickness and boundary layer effects on aerodynamic coefficients. Computational results are given for a symmetrical section represented by an NACA 0012 profile, a conventional section illustrated by an NACA 65A413 profile, and a supercritical type section for general aviation applications typified by a NASA LS(1)-0413 section. The four codes are compared and contrasted in the areas of method of approach, range of applicability, agreement among each other and with experiment, individual advantages and disadvantages, computer run times and memory requirements, and operational idiosyncrasies.

Freuler, R. J.↗
NASA NTRS1979

Reduction of computer usage costs in predicting unsteady aerodynamic loadings caused by control surface motions: Computer program description

A digital computer program was developed to calculate unsteady loadings caused by motions of lifting surfaces with leading edge and trailing edge controls based on the subsonic kernel function approach. The pressure singularities at hinge line and side edges were extracted analytically as a preliminary step to solving the integral equation of collocation. The program calculates generalized aerodynamic forces for user supplied deflection modes. Optional intermediate output includes pressure at an array of points, and sectional generalized forces. From one to six controls on the half span can be accomodated.

Petrarca, J. R.↗
NASA NTRS1983

On obtaining the forward phase functions of Saturn ring features from radio occultation observations

It is noted that the near-forward scattering functions of particles in Saturn ring features are related to 3.6 cm radio occultation power spectra by a Fredholm integral equation of the first kind. The equation reduces to an algebraic system of equation whose solution by usual inversion techniques (that is, least mean squares) is ruled out by the near singularity of the forward transformation matrix. A combination of constrained linear inversion and a filtering algorithm based on eigenvector decomposition of the matrix reduces the instabilities; this yields derived phase functions valid over the range of zero to about 12 mrad. It is noted that these functions represent the collective forward diffraction lobe of particles greater than about 1 m in radius. Since multiple scattering of the signal is a significant effect, the measured phase functions must be adjusted to obtain the singly scattered component. This single-scattering correction is examined for two physical models, namely the monolayer and the classical discrete random slab; in addition, the fraction of opacity in submeter particles for each model for particular ring features is estimated.

Zebker, H. A.↗
NASA NTRS1986

Existence of needle crystals in local models of solidification

The way in which surface tension acts as a singular perturbation to destroy the continuous family of needle-crystal solutions of the steady-state growth equations is analyzed in detail for two local models of solidification. All calculations are performed in the limit of small surface tension or, equivalently, small velocity. The basic mathematical ideas are introduced in connection with a quasilinear, isotropic version of the geometrical model of Brower et al., in which case the continuous family of solutions dissappears completely. The formalism is then applied to a simplified boundary-layer model with an anisotropic kinetic attachment coefficient. In the latter case, the solvability condition for the existence of needle crystals can be satisfied whenever the coefficient of anisotropy is arbitrarily small but nonzero.

Langer, J. S.↗
NASA NTRS1989

Semi-Lagrangian integration of a gridpoint shallow water model on the sphere

A stable, semi-Lagrangian, semi-implicit, two-time-level, gridpoint integration scheme for the shallow water equations on the sphere is presented. A rotated spherical coordinate system is used to integrate the equations of motion at each gridpoint poleward of a certain latitude, thus overcoming problems associated with the polar singularity. The results of medium term integrations of large scale test patterns using a long time step are presented.

Mcdonald, A.↗
NASA NTRS2004

Position-Finding Instrument Built Around a Magnetometer

A coarse-positioning instrument is built around a three-axis magnetometer. The magnetometer is of a type that is made of inexpensive hardware and is suitable for use aboard spacecraft orbiting no more than 1,000 km above the surface of the Earth. A data processor programmed with suitable software and equipped with a central processing unit, random-access memory, programmable read-only memory, and interface circuitry for communication with external equipment are added to the basic magnetometer to convert it into a coarse-positioning instrument. Although the instrument was conceived for use aboard spacecraft, it could be useful for navigation on Earth under some circumstances. A major feature of the proposed instrument is an ability to generate a coarse estimate of its position in real time (that is, without start-up delay). Algorithms needed to solve the position equations have been developed. These include algorithms to work around gaps in measurement data that arise from a singularity near the minimum in the magnetic field of the Earth. Some work has been done to develop a prototype of this instrument incorporating a standard three-axis flux-gate magnetometer and a Pentium P-5 (or equivalent) processor with a clock frequency of 120 MHz. Alternatively, the processor could be of the 486 class. A computer model of the instrument has been completed and tested.

Ketchum, Eleanor↗
NASA NTRS1981

General relaxation schemes in multigrid algorithms for higher order singularity methods

Relaxation schemes based on approximate and incomplete factorization technique (AF) are described. The AF schemes allow construction of a fast multigrid method for solving integral equations of the second and first kind. The smoothing factors for integral equations of the first kind, and comparison with similar results from the second kind of equations are a novel item. Application of the MD algorithm shows convergence to the level of truncation error of a second order accurate panel method.

Oskam, B.↗
NASA NTRS1975

Prediction of unsteady aerodynamic loadings caused by leading edge and trailing edge control surface motions in subsonic compressible flow: Computer program description

A digital computer program has been developed to calculate unsteady loadings caused by motions of lifting surfaces with leading edge or trailing edge controls based on the subsonic kernel function approach. The pressure singularities at hinge line and side edges have been extracted analytically as a preliminary step to solving the integral equation by collocation. The program calculates generalized aerodynamic forces for user supplied deflection modes. Optional intermediate output includes pressure at an array of points, and sectional generalized forces. From one to six controls on the half span can be accommodated.

Redman, M. C.↗
NASA NTRS1987

A new approach to the effect of sound on vortex dynamics

Analytical results are presented on the effect of acoustic radiation on three-dimensional vortex motions in a homogeneous, slightly compressible, inviscid fluid. The flow is considered as linear and irrotational everywhere except inside a very thin cylindrical core region around the vortex filament. In the outside region, a velocity potential is introduced that must be multivalued, and it is shown how to compute this scalar potential if the motion of the vortex filament is prescribed. To find the motion of this singularity in an external potential flow, a variational principle involving a volume integral that must exclude the singular region is considered. A functional of the external potential and vortex filament position is obtained whose extrema give equations to determine the sought-after evolution. Thus, a generalization of the Biot-Savart law to flows with constant sound speed at low Mach number is obtained.

Lund, Fernando↗
NASA NTRS1992

Aerodynamics of loaded cascades in subsonic flows subject to unsteady three-dimensional vortical disturbances

A highly efficient numerical method is developed for three-dimensional, periodic, vortical flows around a cascade of loaded airfoils. The method linearizes Euler's equations about the mean flow of the cascade and thus fully accounts for the effects of distortion of the vortical disturbances as they propogate and interact with the cascade mean flow. The numerical scheme is based on splitting the unsteady velocity into vortical and potential parts. The latter is governed by a non-constant coefficient inhomogeneous convective wave equation. A new and computationally suitable out-flow conditions are derived and avoid the difficulties associated with the singular velocity downstream. Solutions were obtained in the frequency domain by using a body-fitted coordinate system. Results are presented to demonstrate the effects of the out-flow boundary conditions, cascade spacing, mean blade loading and gust upstream conditions on the aerodynamics response and unsteady pressure field of a cascade.

Fang, J.↗
NASA NTRS1995

Ball bearing vibrations amplitude modeling and test comparisons

Bearings generate disturbances that, when combined with structural gains of a momentum wheel, contribute to induced vibration in the wheel. The frequencies generated by a ball bearing are defined by the bearing's geometry and defects. The amplitudes at these frequencies are dependent upon the actual geometry variations from perfection; therefore, a geometrically perfect bearing will produce no amplitudes at the kinematic frequencies that the design generates. Because perfect geometry can only be approached, emitted vibrations do occur. The most significant vibration is at the spin frequency and can be balanced out in the build process. Other frequencies' amplitudes, however, cannot be balanced out. Momentum wheels are usually the single largest source of vibrations in a spacecraft and can contribute to pointing inaccuracies if emitted vibrations ring the structure or are in the high-gain bandwidth of a sensitive pointing control loop. It is therefore important to be able to provide an a priori knowledge of possible amplitudes that are singular in source or are a result of interacting defects that do not reveal themselves in normal frequency prediction equations. This paper will describe the computer model that provides for the incorporation of bearing geometry errors and then develops an estimation of actual amplitudes and frequencies. Test results were correlated with the model. A momentum wheel was producing an unacceptable 74 Hz amplitude. The model was used to simulate geometry errors and proved successful in identifying a cause that was verified when the parts were inspected.

Hightower, Richard A., III↗
NASA NTRS1994

The Nonlinear Evolution of the Inviscid Secondary Instability of Streamwise Vortex Structures

The weakly nonlinear evolution of an inviscid marginally unstable wave growing on a boundary layer supporting a streamwise vortex structure is investigated. The nonlinear growth of the wave is found to be controlled by the diffusion layer located at the edge of the critical layer associated with the wave. The evolution equation is found to depend on the upstream history of the wave and the solution of the equation suggests that the wave either restructures the mean state so as to make it stable or develops a singularity at a finite distance downstream of the point of neutral stability.

Blackaby, Nicholas D.↗
NASA NTRS1974

Quantum descriptions of singularities leading to pair creation

A class of cosmological models is analyzed which provide a mathematically convenient (but idealized) description of a cosmological singularity that develops into a pair creation epoch and terminates in an adiabatic expansion with redshifting particle energies. This class of models was obtained by Gowdy (1971, 1974) as a set of exact solutions of the classical empty space Einstein equations describing inhomogeneous universes populated only by gravitational waves. It is shown that these models can be used to exhibit simplified models of quantized gravitational fields, and that a quantum description can be given arbitrarily near a cosmological singularity. Graviton pair creation occurs, and can be seen to convert anisotropic expansion rates into the energy of graviton pairs.

Misner, C. W.↗
NASA NTRS2014

Comptonization in Ultra-Strong Magnetic Fields: Numerical Solution to the Radiative Transfer Problem

We consider the radiative transfer problem in a plane-parallel slab of thermal electrons in the presence of an ultra-strong magnetic field (B approximately greater than B(sub c) approx. = 4.4 x 10(exp 13) G). Under these conditions, the magnetic field behaves like a birefringent medium for the propagating photons, and the electromagnetic radiation is split into two polarization modes, ordinary and extraordinary, that have different cross-sections. When the optical depth of the slab is large, the ordinary-mode photons are strongly Comptonized and the photon field is dominated by an isotropic component. Aims. The radiative transfer problem in strong magnetic fields presents many mathematical issues and analytical or numerical solutions can be obtained only under some given approximations. We investigate this problem both from the analytical and numerical point of view, provide a test of the previous analytical estimates, and extend these results with numerical techniques. Methods. We consider here the case of low temperature black-body photons propagating in a sub-relativistic temperature plasma, which allows us to deal with a semi-Fokker-Planck approximation of the radiative transfer equation. The problem can then be treated with the variable separation method, and we use a numerical technique to find solutions to the eigenvalue problem in the case of a singular kernel of the space operator. The singularity of the space kernel is the result of the strong angular dependence of the electron cross-section in the presence of a strong magnetic field. Results. We provide the numerical solution obtained for eigenvalues and eigenfunctions of the space operator, and the emerging Comptonization spectrum of the ordinary-mode photons for any eigenvalue of the space equation and for energies significantly lesser than the cyclotron energy, which is on the order of MeV for the intensity of the magnetic field here considered. Conclusions. We derived the specific intensity of the ordinary photons, under the approximation of large angle and large optical depth. These assumptions allow the equation to be treated using a diffusion-like approximation.

acceleration of particles↗
NASA NTRS1973

Maximally slicing a black hole.

Analytic and computer-derived solutions are presented of the problem of slicing the Schwarzschild geometry into asymptotically flat, asymptotically static, maximal spacelike hypersurfaces. The sequence of hypersurfaces advances forward in time in both halves (u greater than or equal to 0, u less than or equal to 0) of the Kruskal diagram, tending asymptotically to the hypersurface r = 3/2 M and avoiding the singularity at r = 0. Maximality is therefore a potentially useful condition to impose in obtaining computer solutions of Einstein's equations.

Estabrook, F.↗
NASA NTRS1975

Second-order small-disturbance solutions for hypersonic flow over power-law bodies

Similarity solutions were found which give the adiabatic flow of an ideal gas about two-dimensional and axisymmetric power-law bodies at infinite Mach number to second order in the body slenderness parameter. The flow variables were expressed as a sum of zero-order and perturbation similarity functions for which the axial variations in the flow equations separated out. The resulting similarity equations were integrated numerically. The solutions, which are universal functions, are presented in graphic and tabular form. To avoid a singularity in the calculations, the results are limited to body power-law exponents greater than about 0.85 for the two-dimensional case and 0.75 for the axisymmetric case. Because of the entropy layer induced by the nose bluntness (for power-law bodies other than cones and wedges), only the pressure function is valid at the body surface. The similarity results give excellent agreement with the exact solutions for inviscid flow over wedges and cones having half-angles up to about 20 deg. They give good agreement with experimental shock-wave shapes and surface-pressure distributions for 3/4-power axisymmetric bodies, considering that Mach number and boundary-layer displacement effects are not included in the theory.

Townsend, J. C.↗
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 NTRS1986

MHD bending waves in a current sheet

Transverse MHD bending waves are considered in an isothermal and compressible two-dimensional current sheet of finite thickness in which the magnetic field changes direction and strength. The general form of the wave equation is obtained. It is shown that rotation of the magnetic field across the current sheet prevents the existence of singular points so that continuous spectrum solutions and the concomitant wave decay disappear. Instead, normal modes exist and closed integral solution for arbitrary current sheet structure are found. The results are discussed in terms of small-scale waves on the heliospheric current sheet.

Musielak, Z. E.↗