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.

SourceResearch connectionFrom yearThrough year

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

At least 631 records · Page 35

All sources NASA NTRS DOE OSTI

NASA NTRS2013

The Force-Free Magnetosphere of a Rotating Black Hole

We revisit the Blandford-Znajek process and solve the fundamental equation that governs the structure of the steady-state force-free magnetosphere around a Kerr black hole. The solution depends on the distributions of the magnetic field angular velocity and the poloidal electric current. These are not arbitrary. They are determined self-consistently by requiring that magnetic field lines cross smoothly the two singular surfaces of the problem: the inner "light surface" located inside the ergosphere and the outer "light surface" which is the generalization of the pulsar light cylinder.We find the solution for the simplest possible magnetic field configuration, the split monopole, through a numerical iterative relaxation method analogous to the one that yields the structure of the steady-state axisymmetric force-free pulsar magnetosphere. We obtain the rate of electromagnetic extraction of energy and confirm the results of Blandford and Znajek and of previous time-dependent simulations. Furthermore, we discuss the physical applicability of magnetic field configurations that do not cross both "light surfaces."

Contopoulos, Ioannis↗

NASA NTRS1974

Fast direct numerical solution of the nonhomogeneous Cauchy-Riemann equations

A fast direct (noniterative) 'Cauchy-Riemann Solver' is developed for solving the finite-difference equations representing systems of first-order elliptic partial differential equations in the form of the nonhomogeneous Cauchy-Riemann equations. The method is second-order accurate and requires approximately the same computer time as a fast cyclic-reduction Poisson solver. The accuracy and efficiency of the direct solver are demonstrated in an application to solving an example problem in aerodynamics: subsonic inviscid flow over a biconvex airfoil. The analytical small-perturbation solution contains singularities, which are captured well by the computational technique. The algorithm is expected to be useful in nonlinear subsonic and transonic aerodynamics.

Lomax, H.↗

NASA NTRS1977

Numerical determination of the fundamental eigenvalue for the Laplace operator on a spherical domain

Methods for obtaining approximate solutions for the fundamental eigenvalue of the Laplace-Beltrami operator (i.e., the membrane eignevalue problem for the vibration equation) on the unit spherical surface are developed. Two types of spherical surface domains are considered: the interior of a spherical triangle, and the exterior of a great circle arc extending for less than pi radians (a spherical surface with a slit). In both cases, zero boundary conditions are imposed. In order to solve the resulting second-order elliptic partial differential equations in two independent variables, a finite difference approximation is employed. The fundamental eigenvalue is approximated by iteration utilizing the power method and point successive overrelaxation. Some numerical results are given and compared, in certain special cases, with analytical solutions to the eigenvalue problem. The significance of the numerical eigenvalue results is discussed in terms of the singularities in the solution of three-dimensional boundary-value problems near a polyhedral corner of the domain.

Walden, H.↗

NASA NTRS1978

A user's guide to the SUDAN computer program for determining the vibration modes of structural systems

The use of the SUDAN computer program for analyzing structural systems for their natural modes and frequencies of vibration is described. SUDAN is intended for structures which can be represented as an equivalent system of beam, spring, and rigid-body substructures. User-written constraint equations are used to analytically join the mass and stiffness matrices of the substructures to form the mass and stiffness matrices of the complete structure from which all the frequencies and modes of the system are determined. The SUDAN program can treat the case in which both the mass and stiffness matrices of the coupled system may be singular simultaneously. A general description of the FORTRAN IV program is given, the computer hardware and software specifications are indicated, and the input required by the program is described.

Kvaternik, R. G.↗

NASA NTRS1981

Numerical solution of transonic stream function equation

The stream function equation, in conservation form, looks similar to the full potential equation and existing methods (e.g. artificial compressibility) can be readily applied. Rotational flows can be calculated once the vorticity (due to shocks or nonuniformity) is evaluated. There are, however, two main difficulties: First, the density is not uniquely determined in terms of the flux (there are two solutions; the subsonic and the supersonic branch with a square root singularity at the sonic point). Methods to overcome this difficulty are studied and results are presented with some remarks on inviscid separation and closed stream lines. Second, the need of two stream functions for three dimensional calculations is briefly discussed.

Hafez, M.↗

NASA NTRS1976

Features of sound propagation through and stability of a finite shear layer

The plane wave propagation, the stability and the rectangular duct mode problems of a compressible inviscid linearly sheared parallel, but otherwise homogeneous flow, are shown to be governed by Whittaker's equation. The exact solutions for the perturbation quantities are essentially Whittaker M-functions. A number of known results are obtained as limiting cases of exact solutions. For the compressible finite thickness shear layer it is shown that no resonances and no critical angles exist for all Mach numbers, frequencies and shear layer velocity profile slopes except in the singular case of the vortex sheet.

Koutsoyannis, S. P.↗

NASA NTRS1980

Forebody and base region real-gas flow in severe planetary entry by a factored implicit numerical method. I - Computational fluid dynamics

A new code for the simulation of full (forebody and base region) flowfields about bluff bodies in the hypersonic regime of severe planetary entry is described. The present 'maximally conservative, maximally differenced' formulation of the unsteady compressible Navier-Stokes equations for 2-D axisymmetric 3-D flow is contrasted for stability with previous formulations of Viviand, Kutler, et al, and Thomas and Lombard. Discrete metric relations peculiar to the axisymmetric finite volume formulation are presented along with a general discussion of their relations to and consequences of failure to close computational cells. A computational mesh of curvilinear coordinate topology singular in the flow regime is presented that permits aligned capturing of the major physical features of the complex flowfield.

Lombard, C. K.↗

NASA NTRS1981

Optimal maneuvers at supersonic speeds in a vertical plane

This paper presents several optimal maneuvers of supersonic aircraft in a vertical plane. The general dimensionless equations of motion are derived and the computation of the optimal trajectories are carried out using the aerodynamics and engine characteristics of a light-weight fighter called the supercruiser. Because of the normalizing of the control variables, namely, the load factor and the thrust-to-weight ratio, the results applied to any supersonic aircraft. The optimality of the singular thrust control and the optimal junction of different subarcs are discussed. The proposed method of computing the optimal trajectory is very efficient and makes explicit the selection of the optimal control. The technique should be useful for performance assessment of supersonic aircraft with potential for implementation of onboard flight control system.

Vinh, N. X.↗

NASA NTRS1980

Implicit Finite-Difference Simulations of Three-Dimensional Compressible Flow

An implicit finite-difference procedure for unsteady three-dimensional flow capable of handling arbitrary geometry through the use of general coordinate transformations is described. Viscous effects are optionally incorporated with a "thin-layer" approximation of the Navier-Stokes equations. An implicit approximate factorization technique is employed so that the small grid sizes required for spatial accuracy and viscous resolution do not impose stringent stability limitations. Results obtained from the program include transonic inviscid or viscous solutions about simple body configurations. Comparisons with existing theories and experiments are made. Numerical accuracy and the effect of three-dimensional coordinate singularities are also discussed.

Pulliam, Thomas H.↗

NASA NTRS1971

Exact solutions to Brans-Dicke cosmologies in flat Friedmann universes.

The Brans-Dicke cosmological equations for flat Friedmann-type expanding universes are solved parametrically for time, density, expansion parameter, and scalar field. These results reduce to a previously obtained exact solution to the radiation cosmology. Although the scalar field may be undetectable at the present epoch, it is felt that, if it exists, it must play an important role as one approaches the initial singularity of the cosmology.

Morganstern, R. E.↗

NASA NTRS1992

A non-iterative method for computing the infimum in H(infinity)-optimization

This paper presents a simple and non-iterative procedure for the computation of the exact value of the infimum in the singular H(infinity)-optimization problem, and is an extension of our earlier work. The problem formulation is general and does not place any restriction on the direct feedthrough terms between the control input and the controlled output variables, and between the disturbance input and the measurement output variables. Our method is applicable to a class of singular H(infinity)-optimization problems for which the transfer functions from the control input to the controlled output and from the disturbance input to the measurement output have no invariant zeros on the j-omega axis and also satisfy certain geometric conditions. The computation of the infimum in our method involves solving two well-defined Riccati and two Liapunov equations.

Chen, Ben M.↗

NASA NTRS1978

On implicit finite-difference simulations of three-dimensional flow

An implicit finite-difference procedure for unsteady three-dimensional flow capable of handling arbitrary geometry through the use of general coordinate transformations is described. Viscous effects are optionally incorporated with a 'thin layer' approximation of the Navier-Stokes equations. An implicit approximate factorization technique is employed so that the small grid sizes required for spatial accuracy and viscous resolution do not impose stringent stability limitations. Results obtained from the program include transonic inviscid and laminar-turbulent solutions about simple body configurations. Comparisons with existing theories and experiments are made. Numerical accuracy and the effect of three-dimensional coordinate singularities are also discussed.

Pulliam, T. H.↗

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 NTRS1993

Elastic interactions of a fatigue crack with a micro-defect by the mixed boundary integral equation method

In this paper, the mixed boundary integral equation method is developed to study the elastic interactions of a fatigue crack and a micro-defect such as a void, a rigid inclusion or a transformation inclusion. The method of pseudo-tractions is employed to study the effect of a transformation inclusion. An enriched element which incorporates the mixed-mode stress intensity factors is applied to characterize the singularity at a moving crack tip. In order to evaluate the accuracy of the numerical procedure, the analysis of a crack emanating from a circular hole in a finite plate is performed and the results are compared with the available numerical solution. The effects of various micro-defects on the crack path and fatigue life are investigated. The results agree with the experimental observations.

Lua, Yuan J.↗

NASA NTRS1983

Mechanics of composite materials: Recent advances; Proceedings of the Symposium, Virginia Polytechnic Institute and State University, Blacksburg, VA, August 16-19, 1982

The present conference on the mechanics of composites discusses microstructure's influence on particulate and short fiber composites' thermoelastic and transport properties, the elastoplastic deformation of composites, constitutive equations for viscoplastic composites, the plasticity and fatigue of metal matrix composites, laminate damping mechanisms, the micromechanical modeling of Kevlar/epoxy composites' time-dependent failure, the variational characterization of waves in composites, and computational methods for eigenvalue problems in composite design. Also discussed are the elastic response of laminates, elastic coupling nonlinear effects in unsymmetrical laminates, elasticity solutions for laminate problems having stress singularities, the mechanics of bimodular composite structures, the optimization of laminated plates and shells, NDE for laminates, the role of matrix cracking in the continuum constitutive behavior of a damaged composite ply, and the energy release rates of various microcracks in short fiber composites.

Hashin, Z.↗

NASA NTRS2003

Many-Body Effects in Quantum-Well Intersubband Transitions

Intersubband polarization couples to collective excitations of the interacting electron gas confined in a semiconductor quantum well (Qw) structure. Such excitations include correlated pair excitations (repellons) and intersubband plasmons (ISPs). The oscillator strength of intersubband transitions (ISBTs) strongly varies with QW parameters and electron density because of this coupling. We have developed a set of kinetic equations, termed the intersubband semiconductor Bloch equations (ISBEs), from density matrix theory with the Hartree-Fock approximation, that enables a consistent description of these many-body effects. Using the ISBEs for a two-conduction-subband model, various many-body effects in intersubband transitions are studied in this work. We find interesting spectral changes of intersubband absorption coefficient due to interplay of the Fermi-edge singularity, subband renormalization, intersubband plasmon oscillation, and nonparabolicity of bandstructure. Our results uncover a new perspective for ISBTs and indicate the necessity of proper many-body theoretical treatment in order for modeling and prediction of ISBT line shape.

Li, Jian-Zhong↗

NASA NTRS2023

Relationship Between the Sub-Micron Fraction (SMF) and Fine-Mode Fraction (FMF) in the Context of AERONET Retrievals

The sub-micron (SM) aerosol optical depth (AOD) is an optical separation based on the fraction of particles below a specified cutoff radius of the particle size distribution (PSD) at a given particle radius. It is fundamentally different from spectrally separated FM (fine-mode) AOD. We present a simple (AOD-normalized) SM fraction versus FM fraction (SMF vs. FMF) linear equation that explains the well-recognized empirical result of SMF generally being greater than the FMF. The AERONET inversion (AERinv) products (combined inputs of spectral AOD and sky radiance) and the spectral deconvolution algorithm (SDA) products (input of AOD spectra) enable, respectively, an empirical SMF vs. FMF comparison at similar (columnar) remote sensing scales across a variety of aerosol types. SMF (AERinv-derived) vs. FMF (SDA-derived) behavior is primarily dependent on the relative truncated portion (ε c ) of the coarse-mode (CM) AOD associated with the cutoff portion of the CM PSD and, to a second order, the cutoff FM PSD and FM AOD (ε f ). The SMF vs. FMF equation largely explains the SMF vs. FMF behavior of the AERinv vs. SDA products as a function of PSD cutoff radius (“inflection point”) across an ensemble of AERONET sites and aerosol types (urban-industrial, biomass burning, dust, maritime and a mixed class of Arctic aerosols). The overarching dynamic was that the linear SMF vs. FMF relation pivots clockwise about the approximate (SMF, FMF) singularity of (1, 1) in a “linearly inverse” fashion (slope and intercept of approximately 1−ε c and ε c ) with increasing cutoff radius. SMF vs. FMF slopes and intercepts derived from AERinv and SDA retrievals confirmed the general domination of ε c over ε f in controlling that dynamic. A more general conclusion is the apparent confirmation that the optical impact of truncating modal (whole) PSD features can be detected by an SMF vs. FMF analysis.

Norman T. O'Neill↗

NASA NTRS1977

NASA Ames three-dimensional potential flow analyses system (POTFAN) boundary condition code (BCDN), version 1

This segment of the POTFAN system is used to generate right hand sides (boundary conditions) of the system of equations associated with the flow field under consideration. These specified flow boundary conditions are encountered in the oblique derivative boundary value problem (boundary value problem of the third kind) and contain the Neumann boundary condition as a special case. Arbitrary angle of attack and/or sideslip and/or rotation rates may be specified, as well as an arbitrary, nonuniform external flow field and the influence of prescribed singularity distributions.

Davis, J. E.↗