SEARCH · Search NASA

Results for “weak 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 289 records · Page 16

All sources NASA NTRS DOE OSTI

NASA NTRS1993

Inertial Currents in Isotropic Plasma

The magnetospheric convection electric field contributes to Birkeland currents. The effects of the field are to polarize the plasma by displacing the bounce paths of the ions from those of electrons, to redistribute the pressure so that it is not constant along magnetic field lines, and to enhance the pressure gradient by the gradient of the bulk speed. Changes in the polarization charge during the convection of the plasma are neutralized by electrons in the form of field-aligned currents that close through the ionosphere. The pressure drives field-aligned currents through its gradient in the same manner as in quasi-static plasma, but with modifications that are important if the bulk speed is of the order of the ion thermal speed; the variations in the pressure along field lines are maintained by a weak parallel potential drop. These effects are described in terms of the field-aligned currents in steady state, isotropic, MED plasma. Solutions are developed by taking the MHD limit of two-fluid solutions and illustrated in the special case of Maxwellian plasma for which the temperature is constant along magnetic field lines. The expression for the Birkeland current density is a generalization of Vasyliunas' expression for the field-aligned current density in quasi-static plasma and provides a unifying expression when both pressure gradients and ion inertia operate simultaneously as sources of field-aligned currents. It contains a full account of different aspects of the ion flow (parallel and perpendicular velocity and vorticity) that contribute to the currents. Contributions of ion inertia to field-aligned currents will occur in regions of strong velocity shear, electric field reversal, or large gradients in the parallel velocity or number density, and may be important in the low-latitude boundary layer, plasma sheet boundary layer, and the inner edge region of the plasma sheet.

Heinemann, M.↗

NASA NTRS2007

On the Propagation and Interaction of Spherical Blast Waves

The characteristics and the scaling laws of isolated spherical blast waves have been briefly reviewed. Both self-similar solutions and numerical solutions of isolated blast waves are discussed. Blast profiles in the near-field (strong shock region) and the far-field (weak shock region) are examined. Particular attention is directed at the blast overpressure and shock propagating speed. Consideration is also given to the interaction of spherical blast waves. Test data for the propagation and interaction of spherical blast waves emanating from explosives placed in the vicinity of a solid propellant stack are presented. These data are discussed with regard to the scaling laws concerning the decay of blast overpressure.

Kandula, Max↗

NASA NTRS2013

Influence of Coupled Radiation and Ablation on the Aerothermodynamic Environment of Planetary Entry Vehicles

A review of recently published coupled radiation and ablation capabilities involving the simulation of hypersonic flowfields relevant to Earth, Mars, or Venus entry is presented. The three fundamental mechanisms of radiation coupling are identified as radiative cooling, precursor photochemistry, and ablation-radiation interaction. The impact of these mechanisms are shown to be significant for a 3 m radius sphere entering Earth at hypothetical Mars return conditions (approximately 15 km/s). To estimate the influence precursor absorption on the radiative flux for a wide range of conditions, a simplified approach is developed that requires only the non-precursor solution. Details of a developed coupled ablation approach, which is capable of treating both massively ablating flowfields in the sublimation regime and weakly ablating diffusion Climited oxidation cases, are presented. A review of the two primary uncoupled ablation approximations, identified as the blowing correction and film coefficient approximations, is made and their impact for Earth and Mars entries is shown to be significant for recession and convective heating predictions. Fully coupled ablation and radiation simulations are presented for the Mars return sphere throughout its entire trajectory. Applying to the Mars return sphere the Pioneer- Venus heritage carbon phenolic heatshield, which has properties available in the open literature, the differences between steady state ablation and coupling to a material response code are shown to be significant.

Johnston, Christopher O.↗

NASA NTRS2021

On Synthetic Absorption Line Profiles of Thermally Driven Winds from Active Galactic Nuclei

The warm absorbers observed in more than half of all nearby active galactic nuclei are tracers of ionized outflows located at parsec-scale distances from the central engine. If the smallest inferred ionization parameters correspond to plasma at a few 104 K, then the gas undergoes a transition from being bound to unbound, provided it is further heated to ∼106 K at larger radii. Dannen et al. recently discovered that, under these circumstances, thermally driven wind solutions are unsteady and even show very dense clumps due to thermal instability. To explore the observational consequences of these new wind solutions, we compute line profiles based on the one-dimensional simulations of Dannen et al. We show how the line profiles from even a simple steady-state wind solution depend on the ionization energy (IE) of absorbing ions, which is a reflection of the wind ionization stratification. To organize the diversity of the line shapes, we group them into four categories: weak Gaussians, saturated boxy profiles with and without an extended blue wing, and broad weak profiles. The lines with profiles in the last two categories are produced by ions with the highest IE that probe the fastest regions. Their maximum blueshifts agree with the highest flow velocities in thermally unstable models, both steady-state and clumpy versions. In contrast, the maximum blueshifts of the highest-IE lines in thermally stable models can be less than half of the actual solution velocities. Clumpy solutions can additionally imprint distinguishable absorption troughs at widely separated velocities.

Shalini Ganguly↗

NASA NTRS1988

A stiffly-stable implicit Runge-Kutta algorithm for CFD applications

A stiffly-stable implicit Runge-Kutta integration algorithm is derived for CFD applications spanning the range of semidiscrete theories. The algorithm family contains the one-step 'theta' algorithms, including backwards Euler and the trapezoidal rule, and provides a versatile framework to identify expressions governing algorithm stability characteristics. Parameters of a Runge-Kutta optimal implicit algorithm, second-order accurate in time and stiffly-stable, are established. This algorithm is implemented within a weak statement finite element semidiscrete formulation for one- and two-dimensional conservation law systems. Numerical results are compared to theta-algorithm solutions, for unsteady quasi-one-dimensional Euler predictions with shocks, and for a specially derived two-dimensional conservation law system modeling the Euler equations.

Baker, A. J.↗

NASA NTRS1991

Relativistic motion in gamma-ray bursts

Three fundamental problems affect models of gamma-ray bursts, i.e., the energy source, the ability of high-energy photons to escape the radiation region, and the comparative weakness of X-ray emission. It is indicated that relativistic bulk motion of the gamma-ray-emitting plasma generically provides a solution to all three of these problems. Results show that, if the plasma that produces gamma-ray bursts has a bulk relativistic velocity with Lorentz factor gamma of about 10, several of the most troubling problems having to do with gamma-ray bursts are solved.

Krolik, Julian H.↗

NASA NTRS1994

Inertial currents in isotropic plasma

The magnetospheric convection electric field contributes to Birkeland currents. The effects of the field are to polarize the plasma by displacing the bounce paths of the ions from those of electrons, to redistribute the pressure so that it is not constant along magnetic field lines, and to enhance the pressure gradient by the gradient of the bulk speed. Changes in the polarization charge during the convection of the plasma are neutralized by electrons in the form of field-aligned currents that close through the ionosphere. The pressure drives field-aligned currents through its gradient in the same manner as in quasi-static plasma, but with modifications that are important if the bulk speed is of the order of the ion thermal speed; the variations in the pressure along field lines are maintained by a weak parallel potential drop. These effects are described in terms of the field-aligned currents in steady state, isotropic, magnetohyrodynamic (MHD) plasma. Solutions are developed by taking the MHD limit of two-fluid solutions and illustrated in the special case of Maxwellian plasma for which the temperature is constant along magnetic field lines. The expression for the Birkeland current density is a generalization of Vasyliunas' expression for the field-aligned current density in quasi-static plasma and provides a unifying expression when both pressure gradients and ion inertia operate simultaneously as sources of field-aligned currents. It contains a full account of different aspects of the ion flow (parallel and perpendicular velocity and vorticity) that contribute to the currents. Contributions of ion inertia to field-aligned currents will occur in regions of strong velocity shear, electric field reversal, or large gradients in the parallel velocity or number density, and may be important in the low-latitude boundary layer, plasma sheet boundary layer, and the inner edge region of the plasma sheet.

Heinemann, M.↗

NASA NTRS2008

The Aeroacoustics of Slowly Diverging Supersonic Jets

This paper is concerned with utilizing the acoustic analogy approach to predict the sound from unheated supersonic jets. Previous attempts have been unsuccessful at making such predictions over the Mach number range of practical interest. The present paper, therefore, focuses on implementing the necessary refinements needed to accomplish this objective. The important effects influencing peak supersonic noise turn out to be source convection, mean flow refraction, mean flow amplification, and source non-compactness. It appears that the last two effects have not been adequately dealt with in the literature. The first of these because the usual parallel flow models produce most of the amplification in the so called critical layer where the solution becomes singular and, therefore, causes the predicted sound field to become infinite as well. We deal with this by introducing a new weakly non parallel flow analysis that eliminates the critical layer singularity. This has a strong effect on the shape of the peak noise spectrum. The last effect places severe demands on the source models at the higher Mach numbers because the retarded time variations significantly increase the sensitivity of the radiated sound to the source structure in this case. A highly refined (non-separable) source model is, therefore, introduced in this paper.

Goldstein, M. E.↗

NASA NTRS2016

Simple Test Functions in Meshless Local Petrov-Galerkin Methods

Two meshless local Petrov-Galerkin (MLPG) methods based on two different trial functions but that use a simple linear test function were developed for beam and column problems. These methods used generalized moving least squares (GMLS) and radial basis (RB) interpolation functions as trial functions. These two methods were tested on various patch test problems. Both methods passed the patch tests successfully. Then the methods were applied to various beam vibration problems and problems involving Euler and Beck's columns. Both methods yielded accurate solutions for all problems studied. The simple linear test function offers considerable savings in computing efforts as the domain integrals involved in the weak form are avoided. The two methods based on this simple linear test function method produced accurate results for frequencies and buckling loads. Of the two methods studied, the method with radial basis trial functions is very attractive as the method is simple, accurate, and robust.

Raju, Ivatury S.↗

NASA NTRS1989

Theoretical investigation of EM wave generation and radiation in the ULF, ELF, and VLF bands by the electrodynamic orbiting tether

The problem of electromagnetic wave generation by an electrodynamic tethered satellite system is important both for the ordinary operation of such systems and for their possible application as orbiting transmitters. The tether's ionospheric circuit closure problem is closely linked with the propagation of charge-carrying electromagnetic wave packets away from the tethered system. Work is reported which represents a step towards a solution to the problem that takes into account the effects of boundaries and of vertical variations in plasma density, collision frequencies, and ion species. The theory of Alfen wave packet generation by an electrodynamic tethered system in an infinite plasma medium is reviewed, and brief summary of previous work on the problem is given. The consequences of the presence of the boundaries and the vertical nonuniformity are then examined. One of the most significant new features to emerge when ion-neutral collisions are taken into account is the coupling of the Alfven waves to the fast magnetosonic wave. This latter wave is important, as it may be confined by vertical variations in the Alfven speed to a sort of leaky ionospheric wave guide, the resonances of which could be of great importance to the signal received on the Earth's surface. The infinite medium solution for this case where the (uniform) geomagnetic field makes an arbitrary angle with the vertical is taken as the incident wave-packet. Even without a full solution, a number of conclusions can be drawn, the most important of which may be that the electromagnetic field associated with the operation of a steady-current tethered system will probably be too weak to detect on the Earth's surface, even for large tethered currents. This is due to the total reflection of the incident wave at the atmospheric boundary and the inability of a steady-current tethered system to excite the ionospheric wave-guide. An outline of the approach to the numerical problem is given. The use of numerical integrations and boundary conditions consistent with a conducting Earth is proposed to obtain the solution for the horizontal electromagnetic field components at the boundary of the ionosphere with the atmospheric cavity.

Estes, Robert D.↗

NASA NTRS1985

On the relation between 'mixing length' and 'direct interaction approximation' theories of turbulence

The capabilities of direct interaction approximations (DIA) and weak-coupling approximations (WCA) in modeling turbulence in plasma physics studies are examined. DIA and WCA formalisms are reviewed and a DIA expression for homogeneous incompressible hydrodynamics is simplified to derive an energy spectrum in an arbitrary wavenumber space. A solution is defined for the expression when applied to long wavelength phenomena generated by an instability. The equation is further generalized to yield a transport equation for modeling the spectrum generated by homogeneous turbulence in a linearly unstable system undergoing nonlinear interaction. Finally, the technique is applied to two-dimensional turbulence in a low pressure, weakly ionized plasma embedded in a homogeneous magnetic field and experiencing instabilities caused by electron density and electric potential gradients, such as found in the ionospheric E region.

Sudan, R. N.↗

NASA NTRS1980

Gravitational waves from rotating and precessing rigid bodies. II - General solutions and computationally useful formulas

A rigid, freely precessing Newtonian body emits gravitational radiation. In this paper the classical-mechanics results for free precession which are needed in order to calculate the weak-field slow-motion quadrupole-moment gravitational waves are reviewed. Within that formalism, algorithms for computing the exact gravitational power radiated and waveforms produced by arbitrary rigid-body freely precessing source are given. Also presented are the dominant terms in series expansions of the waveforms for the case of an almost-spherical object precessing with a small wobble angle. These series expansions, which retain the precise frequency dependence of the waves, may be useful for gravitational astronomers when freely precessing sources begin to be observed.

Zimmermann, M.↗

NASA NTRS1980

The approximation of radiative effects in relativistic gravity - Gravitational radiation reaction and energy loss in nearly Newtonian systems

An argument is presented to determine the accuracy with which a solution of Einstein's field equations of gravitation must be approximated in order to describe the dominant effects of gravitational radiation emission from weak-field systems. Several previous calculations are compared in the light of this argument, and some apparent discrepancies among them are resolved. The majority of these calculations support the 'quadrupole formulae' for gravitational radiation energy loss and radiation reaction.

Walker, M.↗

NASA NTRS1983

A coupled ice-ocean model of upwelling in the marginal ice zone

A dynamical coupled ice-ocean numerical model for the marginal ice zone (MIZ) is suggested and used to study upwelling dynamics in the MIZ. The nonlinear sea ice model has a variable ice concentration and includes internal ice stress. The model is forced by stresses on the air/ocean and air/ice surfaces. The main coupling between the ice and the ocean is in the form of an interfacial stress on the ice/ocean interface. The ocean model is a linear reduced gravity model. The wind stress exerted by the atmosphere on the ocean is proportional to the fraction of open water, while the interfacial stress ice/ocean is proportional to the concentration of ice. A new mechanism for ice edge upwelling is suggested based on a geostrophic equilibrium solution for the sea ice medium. The upwelling reported in previous models invoking a stationary ice cover is shown to be replaced by a weak downwelling due to the ice motion. Most of the upwelling dynamics can be understood by analysis of the divergence of the across ice edge upper ocean transport. On the basis of numerical model, an analytical model is suggested that reproduces most of the upwelling dynamics of the more complex numerical model.

Roed, L. P.↗

NASA NTRS1990

Reynolds stress closure in jet flows using wave models

A collection of papers is presented. The outline of this report is as follows. Chapter three contains a description of a weakly nonlinear turbulence model that was developed. An essential part of the application of such a closure scheme to general geometry jets is the solution of the local hydrodynamic stability equation for a given jet cross-section. Chapter four describes the conformal mapping schemes used to map such geometries onto a simple computational domain. Chapter five describes a solution of a stability problem for circular, elliptic, and rectangular geometries. In chapter six linear models for the shock shell structure in non-circular jets is given. The appendices contain reprints of papers also published during this study including the following topics: (1) instability of elliptic jets; (2) a technique for predicting the shock cell structure in non-circular jets using a vortex sheet model; and (3) the resonant interaction between twin supersonic jets.

Morris, Philip J.↗

NASA NTRS1991

Time-domain finite elements in optimal control with application to launch-vehicle guidance

A time-domain finite element method is developed for optimal control problems. The theory derived is general enough to handle a large class of problems including optimal control problems that are continuous in the states and controls, problems with discontinuities in the states and/or system equations, problems with control inequality constraints, problems with state inequality constraints, or problems involving any combination of the above. The theory is developed in such a way that no numerical quadrature is necessary regardless of the degree of nonlinearity in the equations. Also, the same shape functions may be employed for every problem because all strong boundary conditions are transformed into natural or weak boundary conditions. In addition, the resulting nonlinear algebraic equations are very sparse. Use of sparse matrix solvers allows for the rapid and accurate solution of very difficult optimization problems. The formulation is applied to launch-vehicle trajectory optimization problems, and results show that real-time optimal guidance is realizable with this method. Finally, a general problem solving environment is created for solving a large class of optimal control problems. The algorithm uses both FORTRAN and a symbolic computation program to solve problems with a minimum of user interaction. The use of symbolic computation eliminates the need for user-written subroutines which greatly reduces the setup time for solving problems.

Bless, Robert R.↗

NASA NTRS1991

Fundamental parameters of brightest cluster galaxies

Photometry and spectroscopy of a sample of the brightest cluster galaxies (BCGs) are presented, and the relationships between scale length, surface brightness, velocity dispersion, and line strength are determined and compared with the relations determined for lower luminosity ellipticals (LLEs). It is found that the relation between log r(e), log sigma, and mean surface brightness reported by Dressler et al. for LLEs also describes the relation for BCGs fairly well. Fundamental plane solutions for BCGs are derived which can be used as distance indicators with a predicted error of about 21 percent per galaxy. Only a weak correlation is found between the Na D 5893 line strength and velocity dispersion for BCGs, in contrast to the significant correlation between these quantities reported for LLEs. It is found that the mean mass-to-light ratio for BCGs is a factor of about two greater than for the LLEs in the Virgo and Coma sample of Dressler et al.

Oegerle, William R.↗

NASA NTRS1996

Linear Instability of a Uni-Directional Transversely Sheared Mean Flow

The effect of spanwise-periodic mean-flow distortions (i.e. streamwise-vortex structures) on the evolution of small-amplitude, single-frequency instability waves in an otherwise two-dimensional shear flow is investigated. The streamwise-vortex structures are taken to be just weak enough so that the spatially growing instability waves behave (locally) like linear perturbations about a uni-directional transversely sheared mean flow. Numerical solutions are computed and discussed for both the mean flow and the instability waves. The influence of the streamwise-vortex wavelength on the properties of the most rapidly growing instability wave is also discussed.

Wundrow, David W.↗