SEARCH · Search NASA

Results for “general relativity equations & solutions”

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 217 records · Page 12

All sources NASA NTRS DOE OSTI

NASA NTRS1991

Plasma expansion characteristics of ionized clouds in the ionosphere - Macroscopic formulation

A macroscopic plasma expansion model, based on a numerical solution of the type-dependent nonlinear coupled continuity and momentum equations for background O(+) ions and several released ion species, was used to examine plasma expansion characteristics of ionized clouds in the ionosphere. Information is obtained on Ba(+), Li(+), and Ba(+) - Li(+) clouds; different cloud sizes (Gaussian half widths of 0.1, 0.5, and 1 km); cloud/background ion density ratios covering two orders of magnitude; electron/ion temperature ratios of 1, 5, and 10; and several cloud-background relative velocities along B (0, 1, 2, 4, 6, and 8 km/sec). The macroscopic expansion features obtained were found to be in general agreement with those obtained from the small-scale numerical simulations.

Schunk, R. W.↗

NASA NTRS1991

Hybrid accretion disks in active galactic nuclei. I - Structure and spectra

A unified treatment is presented of the two distinct states of vertically thin AGN accretion disks: a cool (about 10 to the 6th K) optically thick solution, and a hot (about 10 to the 9th K) optically thin solution. A generalized formalism and a new radiative cooling equation valid in both regimes are introduced. A new luminosity limit is found at which the hot and cool alpha solutions merge into a single solution of intermediate optical depth. Analytic solutions for the disk structure are given, and output spectra are computed numerically. This is used to demonstrate the prospect of fitting AGN broadband spectra containing both the UV bump as well as the hard X-ray and gamma-ray tail, using a single accretion disk model. Such models are found to make definite predictions about the observed spectrum, such as the relation between the hard X-ray spectral index, the UV-to-X-ray luminosity ratio, and a feature of about 1 MeV.

Wandel, Amri↗

NASA NTRS1950

The solution of the laminar-boundary-layer equation for the flat plate for velocity and temperature fields for variable physical properties and for the diffusion field at high concentration

In connection with Pohlhausen's solution for the temperature field on the flat plate, a series of formulas were indicated by means of which the velocity and temperature field for variable physical characteristics can be computed by an integral equation and an iteration method based on it. With it, the following cases were solved: On the assumption that the viscosity simply varies with the temperature while the other fluid properties remain constant, the velocity and temperature field on the heated and cooled plate, respectively, was computed at the Prandtl numbers 12.5 and 100 (viscous fluids). A closer study of these two cases resulted in general relations: The calculations for a gas of Pr number 0.7 (air) were conducted on the assumption that all fluid properties vary with the temperature, and the velocities are low enough for the heat of friction to be discounted. The result was a thickening of the boundary layers, but no appreciable modification in shearing stress or heat-transfer coefficient.

BOUNDARY LAYER, LAMINAR - PLATES, FLAT↗

NASA NTRS1982

Jacchia-Lineberry upper atmosphere density model

A series of increasingly accurate models which are a careful blend of empirical and theoretical formulae were developed. The exospheric temperature is assumed to be a function of: (1) the average and daily variations in the solar flux, (2) the average and three hourly variations in the geomagnetic index, (3) the angle between the position vector and the axis of the unsymmetric atmospheric bulge, and (4) the angle between the position vector and the geomagnetic pole. The exospheric temperature is related to the density by the solution of the diffusion equilibrium equations for the different constituents of the atmosphere as a function of altitude. Other variations are modeled directly as changes in the density. They are: (1) changes due to the semiannual effect, and (2) changes due to the seasonal latitudinal effect. The causes for these variations are not exactly known but may be modeled sufficiently by empirical formulae. The Jacchia model is assumed to be valid over the altitude range of 90 to 2500 km. The residuals between the observed density from satellite drag observations and the computed densities show the mean relative error to be generally less than 10 percent with occasional peak errors near 50 percent.

Mueller, A. C.↗

NASA NTRS1986

Satellite attitude motion models for capture and retrieval investigations

The primary purpose of this research is to provide mathematical models which may be used in the investigation of various aspects of the remote capture and retrieval of uncontrolled satellites. Emphasis has been placed on analytical models; however, to verify analytical solutions, numerical integration must be used. Also, for satellites of certain types, numerical integration may be the only practical or perhaps the only possible method of solution. First, to provide a basis for analytical and numerical work, uncontrolled satellites were categorized using criteria based on: (1) orbital motions, (2) external angular momenta, (3) internal angular momenta, (4) physical characteristics, and (5) the stability of their equilibrium states. Several analytical solutions for the attitude motions of satellite models were compiled, checked, corrected in some minor respects and their short-term prediction capabilities were investigated. Single-rigid-body, dual-spin and multi-rotor configurations are treated. To verify the analytical models and to see how the true motion of a satellite which is acted upon by environmental torques differs from its corresponding torque-free motion, a numerical simulation code was developed. This code contains a relatively general satellite model and models for gravity-gradient and aerodynamic torques. The spacecraft physical model for the code and the equations of motion are given. The two environmental torque models are described.

Cochran, John E., Jr.↗

NASA NTRS1988

Numerical Simulation of Subsonic and Transonic Propeller Flow

The numerical simulation of 3-D transonic flow about a system of propeller blades is investigated. In particular, it is shown that the use of helical coordinates significantly simplifies the form of the governing equation when the propeller system is assumed to be surrounded by an irrotational flow field of an inviscid fluid. The unsteady small disturbance equation, valid for lightly loaded blades and expressed in helical coordinates, is derived from the general blade-fixed potential equation, given for an arbitrary coordinate system. The use of a coordinate system which inherently adapts to the mean flow results in a disturbance equation requiring relatively few terms to accurately model the physics of the flow. Furthermore, the helical coordinate system presented here is novel in that it is periodic in the circumferential direction while, simultaneously, maintaining orthogonal properties at the mean blade locations. The periodic characteristic allows a complete cascade of blades to be treated, and the orthogonality property affords straightforward treatment of blade boundary conditions. An ADI numerical scheme is used to compute the solution of the steady flow as an asymptotic limit of an unsteady flow. As an example of the method, solutions are presented for subsonic and transonic flow about a 5 percent thick bicircular arc blade of an 8-bladed cascade. Both high and low advance ratio cases are computed and include a lifting as well as nonlifting cases. The nonlifting solutions obtained are compared to solutions from a Euler code.

Snyder, Aaron↗

NASA NTRS1993

On the stability of nongyrotropic ion populations - A first (analytic and simulation) assessment

The wave and dispersion equations for perturbations propagating parallel to an ambient magnetic field in magnetoplasmas with nongyrotropic ion populations show, in general, the occurrence of coupling between the parallel (left- and right-hand circularly polarized electromagnetic and longitudinal electrostatic) eigenmodes of the associated gyrotropic medium. These interactions provide a means to driving linearly one mode with free-energy sources of other modes in homogeneous media. Different types of nongyrotropy bring about distinct classes of coupling. The stability of a hydrogen magnetoplasma with anisotropic, nongyrotropic protons that only couple the electromagnetic modes to each other is investigated analytically (via solution of the derived dispersion equation) and numerically (via simulation with a hybrid code). Nongyrotropy enhances growth and enlarges the unstable spectral range relative to the corresponding gyrotropic situation. The relevance of the properties of nongyrotropic populations to space plasma environments is also discussed.

Brinca, A. L.↗

NASA NTRS2021

Transition Analysis for the CRM-NLF Wind Tunnel Configuration

This paper presents the results of an ongoing study into the linear stability characteristics of the boundary layer flow over the common research model with natural laminar flow (CRMNLF) aircraft configuration. The flow conditions match selected test conditions from a recent wind tunnel experiment in the National Transonic Facility at the NASA Langley Research Center. Previous work involving parallel stability computations of a boundary layer flow based on the conical flow approximation has shown that the measured onset of laminar-turbulent transition during the experiments can be correlated with the linear amplification of Tollmien- Schlichting (TS) and stationary crossflow (CF) instabilities in the swept wing boundary layer. Here, we examine the effects of the simplifying approximations in both basic state computation and the stability analysis, with the goal of quantifying the resulting changes in the N-factor correlations. Specifically, the basic states are computed by using full Navier-Stokes equations and the stability analysis is performed by using a nonorthogonal coordinate system that allows a clear distinction between planar TS and CF instabilities. Furthermore, the effects of curvature and nonparallel mean flow have been included in the stability computations based on the parabolized stability equations (PSE). The fully turbulent Reynolds-Averaged-Navier-Stokes (RANS) mean flow solutions show good agreement with the measured wall pressure distribution. Viscous-inviscid interactive effects are observed to be important because the shock fronts along the suction surface are influenced by the imposed transition front. The stability results confirm the previous findings related to TS amplification within the inboard region of the wing and the dominance of stationary CF modes in the outboard region. However, given the close proximity of the measured transition front and the dual shock system within the outer part of the wing, the onset of transition may well be shock limited within the outboard region. In general, the transition criterion based on the dual N-factor method with N TS = N CF = 6 is reasonably successful at correlating with the measured transition fronts at Re MAC = 15 million and AoA = 1.5, 2 degrees; however, the low values of the correlating N-factors at Re MAC = 17.5 million support the hypothesis that the measured transition at the higher Reynolds number may have been strongly influenced by the merging of turbulent wedges that originate from surface imperfections near the leading edge.

Boundary layer transition↗

NASA NTRS1987

Observation of millimeter-wave oscillations from resonant tunneling diodes and some theoretical considerations of ultimate frequency limits

Recent observations of oscillation frequencies up to 56 GHz in resonant tunneling structures are discussed in relation to calculations by several authors of the ultimate frequency limits of these devices. It is found that calculations relying on the Wentzel-Kramers-Brillouin (WKB) approximation give limits well below the observed oscillation frequencies. Two other techniques for calculating the upper frequency limit were found to give more reasonable results. One method employs the solution of the time-dependent Schroedinger equation obtained by Kundrotas and Dargys (1986); the other uses the energy width of the transmission function for electrons through the double-barrier structure. This last technique is believed to be the most accurate since it is based on general results for the lifetime of any resonant state. It gives frequency limits on the order of 1 THz for two recently fabricated structures. It appears that the primary limitation of the oscillation frequency for double-barrier resonant-tunneling diodes is imposed by intrinsic device circuit parameters and by the transit time of the depletion layer rather than by time delays encountered in the double-barrier region.

Sollner, T. C. L. G.↗

NASA NTRS1974

First order latitude effects in the solar wind

The Weber-Davis model of the solar wind is generalized to include the effects of latitude. The principal assumptions of perfect electrical conductivity, rotational symmetry, a polytropic relation between pressure and density, and a flow aligned magnetic field in a system rotating with the Sun, are retained. The original three dimensional magnetohydrodynamic flow problem is reduced to a two dimensional hydrodynamic flow problem. The solution at 1 Au is most sensitive to a latitudinal dependence in the coronal boundary temperature and least sensitive to a latitudinal dependence in the magnetic field magnitude. A solution obtained for an approximate dipolar variation in the coronal magnetic field magnitude predicts that the latitudinal flow is initially toward the equator due to magnetic channeling; however, this effect is rapidly overcome and the latitudinal flow at 1 Au is toward the pole and not significantly different from the solution for constant boundary conditions.

Winge, C. R., Jr.↗

NASA NTRS1986

Fugacity and concentration gradients in a gravity field

Equations are reviewed which show that at equilibrium fugacity and concentration gradients can exist in gravitational fields. At equilibrium, the logarithm of the ratio of the fugacities of a species at two different locations in a gravitational field is proportional to the difference in the heights of the two locations and the molecular weight of the species. An analogous relation holds for the concentration ratios in a multicomponent system. The ratio is calculated for a variety of examples. The kinetics for the general process are derived, and the time required to approach equilibrium is calculated for several systems. The following special topics are discussed: ionic solutions, polymers, multiphase systems, hydrostatic pressure, osmotic pressure, and solubility gradients in a gravity field.

May, C. E.↗

NASA NTRSDate not supplied

Improvements in Iterative Convergence of FUN3D Solutions

This paper presents a hierarchical adaptive nonlinear iteration method (HANIM) implemented in the NASA computational fluid dynamics code, FUN3D, to improve robustness and computational efficiency. In contrast to the legacy FUN3D iterative solver that relies on an approximate Jacobian, a simple multicolor Gauss-Seidel point-implicit iteration scheme, and linear Courant-Friedrichs-Lewy number (CFL) ramping, HANIM is based upon a hierarchy of modules including preconditioner, generalized conjugate residual, realizability check, nonlinear control, and CFL adaption modules. HANIM performance is systematically compared with the performance of the legacy solver of FUN3D and a baseline solver based on a preconditioner alone. Iterative solutions are compared for three benchmark cases: a subsonic separated flow around a hemisphere cylinder, a supersonic flow through a long duct, and a subsonic flow over the NASA wing-fuselage juncture model. Two Reynolds-averaged Navier-Stokes turbulence models are used in these computations, namely, the negative variant of the linear one-equation Spalart-Allmaras model and its nonlinear extension based on quadratic constitutive relations.

CFD↗

NASA NTRS1990

Concerning the interaction of non-stationary cross-flow vortices in a three-dimensional boundary layer

Recently there has been much work devoted to considering some of the many and varied interaction mechanisms which may be operative in three-dimensional boundary layer flows. This paper is concerned with resonant triads of crossflow vortices. The effects of interactions upon resonant triads is examined where each member of the triad has the property of being linearly neutrally stable so that the importance of the interplay between modes can be relatively easily assessed. Modes within the boundary layer flow above a rotating disc are investigated because of the similarity between this disc flow and many important practical flows and, secondly, because the selected flow is an exact solution of the Navier-Stokes equations which makes its theoretical analysis especially attractive. It is demonstrated that the desired triads of linearly neutrally stable modes can exist within the chosen boundary layer flow. Evolution equations are obtained to describe the development of the amplitudes of these modes once the interaction mechanism is accounted for. It is found that the coefficients of the interaction terms within the evolution equations are, in general, given by quite intricate expressions although some elementary numerical work shows that the evaluation of these coefficients is practicable. The basis of the work lends itself to generalization to more complicated boundary layers, and effects of detuning or non-parallelism could be provided for within the asymptotic framework.

Bassom, Andrew P.↗

NASA NTRS2017

On High-Order Upwind Methods for Advection

In the fourth installment of the celebrated series of five papers entitled "Towards the ultimate conservative difference scheme", Van Leer (1977) introduced five schemes for advection, the first three are piecewise linear, and the last two, piecewise parabolic. Among the five, scheme I, which is the least accurate, extends with relative ease to systems of equations in multiple dimensions. As a result, it became the most popular and is widely known as the MUSCL scheme (monotone upstream-centered schemes for conservation laws). Schemes III and V have the same accuracy, are the most accurate, and are closely related to current high-order methods. Scheme III uses a piecewise linear approximation that is discontinuous across cells, and can be considered as a precursor of the discontinuous Galerkin methods. Scheme V employs a piecewise quadratic approximation that is, as opposed to the case of scheme III, continuous across cells. This method is the basis for the on-going "active flux scheme" developed by Roe and collaborators. Here, schemes III and V are shown to be equivalent in the sense that they yield identical (reconstructed) solutions, provided the initial condition for scheme III is defined from that of scheme V in a manner dependent on the CFL number. This equivalence is counter intuitive since it is generally believed that piecewise linear and piecewise parabolic methods cannot produce the same solutions due to their different degrees of approximation. The finding also shows a key connection between the approaches of discontinuous and continuous polynomial approximations. In addition to the discussed equivalence, a framework using both projection and interpolation that extends schemes III and V into a single family of high-order schemes is introduced. For these high-order extensions, it is demonstrated via Fourier analysis that schemes with the same number of degrees of freedom 𝐾 per cell, in spite of the different piecewise polynomial degrees, share the same sets of eigenvalues and thus, have the same stability and accuracy. Moreover, these schemes are accurate to order 2𝐾−1, which is higher than the expected order of 𝐾.

high-order methods↗

NASA NTRS1978

Mass loss and OH maser emission from Mira variables

The mass, radius, and luminosity of 26 Mira variables that are known OH sources of radio emission at 1612 MHz have been estimated. The time-independent solution of Salpeter's (1974) stellar-wind equation and a period-density relation are used to solve for basic stellar parameters, with the aid of the terminal expansion velocity of the OH maser cloud. Masses obtained from these calculations are consistent with other estimated values for masses of Mira variables. Good agreement is obtained when comparing the rate of mass loss as determined from Reimers's (1977) semiempirical relation with estimates of the mass-loss rate as deduced from theoretical models involving radiation pressure on grains. These calculations suggest a strong correlation between the mass-loss rate and the pulsation period. Arguments concerning the general properties of silicate grains from radiation-pressure-driven stellar-wind equations are discussed.

Michalitsianos, A. G.↗

NASA NTRS1982

Numerical grid generation; Proceedings of the Symposium on Numerical Generation of Curvilinear Coordinate Systems and Their Use in the Numerical Solution of Partial Differential Equations, Nashville, TN, April 13-16, 1982

General curvilinear coordinate systems are considered along with the error induced by coordinate systems, basic differential models for coordinate generation, elliptic grid generation, conformal grid generation, algebraic grid generation, orthogonal grid generation, patched coordinate systems, and solid mechanics applications of boundary fitted coordinate systems. Attention is given to coordinate system control and adaptive meshes, the application of body conforming curvilinear grids for finite difference solution of external flow, the use of solution adaptive grids in solving partial differential equations, adaptive gridding for finite difference solutions to heat and mass transfer problems, and the application of curvilinear coordinate generation techniques to the computation of internal flows. Other topics explored are related to the solution of nonlinear water wave problems using boundary-fitted coordinate systems, the numerical modeling of estuarine hydrodynamics on a boundary-fitted coordinate system, and conformal grid generation for multielement airfoils.

Thompson, J. F.↗

NASA NTRS1991

Applications of the conjugate gradient FFT method in scattering and radiation including simulations with impedance boundary conditions

The theoretical and computational aspects related to the application of the Conjugate Gradient FFT (CGFFT) method in computational electromagnetics are examined. The advantages of applying the CGFFT method to a class of large scale scattering and radiation problems are outlined. The main advantages of the method stem from its iterative nature which eliminates a need to form the system matrix (thus reducing the computer memory allocation requirements) and guarantees convergence to the true solution in a finite number of steps. Results are presented for various radiators and scatterers including thin cylindrical dipole antennas, thin conductive and resistive strips and plates, as well as dielectric cylinders. Solutions of integral equations derived on the basis of generalized impedance boundary conditions (GIBC) are also examined. The boundary conditions can be used to replace the profile of a material coating by an impedance sheet or insert, thus, eliminating the need to introduce unknown polarization currents within the volume of the layer. A general full wave analysis of 2-D and 3-D rectangular grooves and cavities is presented which will also serve as a reference for future work.

Barkeshli, Kasra↗

NASA NTRS1989

Direct simulation of rarefied hypersonic flows

As the capability of the space transportation vehicles (STV's) expand to meet the requirements for future space exploration and utilization, the effects of rarefied hypersonic flows will play a more significant role in defining the aerodynamic and aerothermodynamic performance of STV's. This is particularly true of the low lift/drag aeroassisted STV's where aerobraking occurs at relatively high altitudes and high velocity. Because of the limitations of the continuum description as expressed by the Navier-Stokes equations and the difficulties of solving the Boltzmann equation, the particle of molecular approach has been developed over the last three decades for modeling rarefied gas effects. The direct simulation Monte Carlo (DSMC) method of Bird is the most used method today for simulating rarefied flows. The DSMC method provides a direct physical simulation as opposed to a numerical solution of a set of model equations. This is accomplished by developing phenomenological models of the relevant physical events. The DSMC method accounts for translational, thermal, chemical, and radiative nonequilibrium effects. The general features of the DSMC method, the numerical requirements for obtaining meaningful results, the modeling used to simulate high temperature gas effects, and applications of the method to calculate the flow about an aeroassist flight experiment vehicle (AFE) are reviewed. The AFE simulates a geosynchronous return while entering the Earth's upper atmosphere at approximately 10 km/s. Results obtained using a general 3-D code are presented for the more rarefied portion of the atmospheric encounter (altitudes of 200 to 100 km) emphasizing surface, flowfield, and aerodynamic characteristics of the AFE. Finally, results obtained using axisymmetric and 1-D versions of the code are presented for lower altitude conditions.

Moss, James N.↗