Search NASA⌕ Search

SEARCH · Search NASA

Results for “ADIABATIC EQUATION”

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 199 records · Page 11

Evolution of dynamo-generated magnetic fields in accretion disks around compact and young stars

Geometrically thin, optically thick, turbulent accretion disks are believed to surround many stars. Some of them are the compact components of close binaries, while the others are throught to be T Tauri stars. These accretion disks must be magnetized objects because the accreted matter, whether it comes from the companion star (binaries) or from a collapsing molecular cloud core (single young stars), carries an embedded magnetic field. In addition, most accretion disks are hot and turbulent, thus meeting the condition for the MHD turbulent dynamo to maintain and amplify any seed field magnetic field. In fact, for a disk's magnetic field to persist long enough in comparison with the disk viscous time it must be contemporaneously regenerated because the characteristic diffusion time of a magnetic field is typically much shorter than a disk's viscous time. This is true for most thin accretion disks. Consequently, studying magentic fields in thin disks is usually synonymous with studying magnetic dynamos, a fact that is not commonly recognized in the literature. Progress in studying the structure of many accretion disks was achieved mainly because most disks can be regarded as two-dimensional flows in which vertical and radial structures are largely decoupled. By analogy, in a thin disk, one may expect that vertical and radial structures of the magnetic field are decoupled because the magnetic field diffuses more rapidly to the vertical boundary of the disk than along the radius. Thus, an asymptotic method, called an adiabatic approximation, can be applied to accretion disk dynamo. We can represent the solution to the dynamo equation in the form B = Q(r)b(r,z), where Q(r) describes the field distribution along the radius, while the field distribution across the disk is included in the vector function b, which parametrically depends on r and is normalized by the condition max (b(z)) = 1. The field distribution across the disk is established rapidly, while the radial distribution Q(r) evolves on a considerably longer timescale. It is this evolution that is the subject of this paper.

Stepinski, Tomasz F.↗

Self-Consistent Superthermal Electron Effects on Plasmaspheric Refilling

The effects of self-consistently including superthermal electrons in the definition of the ambipolar electric field are investigated for the case of plasmaspheric refilling after a geomagnetic storm. By using the total electron population in the hydrodynamic equations, a method for incorporating superthermal electron parameters in the electric field and electron temperature calculation is developed. Also, the ambipolar electric field is included in the kinetic equation for the superthermal electrons through a change of variables using the total energy and the first adiabatic invariant. Calculations based on these changes are performed by coupling time-dependent models of the thermal plasma and superthermal electrons. Results from this treatment of the electric field and the self-consistent development of the solution are discussed in detail. Specifically, there is a decreased thermal electron density in the plasmasphere during the first few minutes of refilling, a slightly accelerated proton shock front, and a decreased superthermal electron flux due to the deceleration by the electric field. The timescales of plasmaspheric refilling are discussed and determined to be somewhat shorter than previously calculated for the thermal plasma and superthermal electron population due to the effects of the field-aligned potential.

Liemohn, M. W.↗

New integration techniques for chemical kinetic rate equations. I - Efficiency comparison

A comparison of the efficiency of several recently developed numerical techniques for solving chemical kinetic rate equations is presented. The solution procedures examined include two general-purpose codes, EPISODE and LSODE, developed as multipurpose differential equation solvers, and three specialzed codes, CHEMEQ, CREK1D, and GCKP84, developed specifically for chemical kinetics. The efficiency comparison is made by applying these codes to two practical combustion kinetics problems. Both problems describe adiabatic, constant-pressure, gas-phase chemical reactions and include all three combustion regimes: induction, heat release, and equilibration. The comparison shows that LSODE is the fastest routine currently available for solving chemical kinetic rate equations. An important finding is that an iterative solution of the algebraic enthalpy conservation equation for temperature can be significantly faster than evaluation of the temperature by integration of its time derivative. Significant increases in computational speed are realized by updating the reaction rate constants only when the temperature change exceeds an amount Delta-T that is problem dependent. An approximate expression for the automatic evaluation of Delta-T is presented and is shown to result in increased computational speed.

Radhakrishnan, K.↗

Numerical solution of the Navier-Stokes equations for laminar, transonic flows

An implicit finite difference solution of the Navier-Stokes equations yielded time histories of the transonic laminar flow development about a circular cylinder and NACA-0018 airfoil. Reynolds numbers ranged from those corresponding to purely laminar flow to those corresponding to significant turbulence in the boundary layer. Body thermal conditions of an adiabatic wall and a specified body temperature were considered. Versatility in treating arbitrary bodies was incorporated by using numerically generated, body-fitted coordinate transformations. Solution of the simultaneous difference equations for the dependent variables was obtained using an accelerated Gauss-Seidel iterative scheme. Computational results are presented in the form of velocity vector fields, Mach number contours, aerodynamic coefficients, heat transfer rates at the body surface, and body temperature distributions. Truncation analyses of first and second derivative difference approximations resulted in general criteria for numerically generated coordinates so that flow near a body is more accurately represented.

Turner, L., III↗

A critical analysis of the accuracy of several numerical techniques for combustion kinetic rate equations

A detailed analysis of the accuracy of several techniques recently developed for integrating stiff ordinary differential equations is presented. The techniques include two general-purpose codes EPISODE and LSODE developed for an arbitrary system of ordinary differential equations, and three specialized codes CHEMEQ, CREK1D, and GCKP4 developed specifically to solve chemical kinetic rate equations. The accuracy study is made by application of these codes to two practical combustion kinetics problems. Both problems describe adiabatic, homogeneous, gas-phase chemical reactions at constant pressure, and include all three combustion regimes: induction, heat release, and equilibration. To illustrate the error variation in the different combustion regimes the species are divided into three types (reactants, intermediates, and products), and error versus time plots are presented for each species type and the temperature. These plots show that CHEMEQ is the most accurate code during induction and early heat release. During late heat release and equilibration, however, the other codes are more accurate. A single global quantity, a mean integrated root-mean-square error, that measures the average error incurred in solving the complete problem is used to compare the accuracy of the codes. Among the codes examined, LSODE is the most accurate for solving chemical kinetics problems. It is also the most efficient code, in the sense that it requires the least computational work to attain a specified accuracy level. An important finding is that use of the algebraic enthalpy conservation equation to compute the temperature can be more accurate and efficient than integrating the temperature differential equation.

Radhadrishnan, Krishnan↗

Calculation of flow distribution in large radius ratio stages of axial flow turbines and comparison of theory and experiment

A method of calculating stage parameters and flow distribution of axial turbines is described. The governing equations apply to space between the blade rows and are based on the assumption of rotationally symmetrical, compressible, adiabatic flow conditions. Results are presented for stage design and flow analysis calculations. Theoretical results from the calculation system are compared with experimental data from low pressure steam turbine tests.

Herzog, J.↗

Effect of wall temperature on a supersonic turbulent boundary layer

Measurement of mean flow profiles in a fully developed Mach 3 turbulent boundary layer with negligible pressure gradient is reported. Data were acquired at several streamwise locations for wall-to-total temperature ratios of 0.94, 0.71, and 0.54. The results demonstrate that the velocity defect formulation of the law-of-the-wake, which successfully correlates compressible, adiabatic boundary layers, is also valid for nonadiabatic flows. It is also shown that for adiabatic walls, the linear Crocco relation between total temperature and velocity does not provide a valid test of the nature of the boundary-layer flow for practical cases where the Prandtl number departs from unity. Finally, the turbulent shear stress, mixing length, and eddy viscosity were extracted from the 'time-averaged' conservation equations using the measured mean flow profiles and found to be insensitive to wall temperature. In particular, the latter properties are in good agreement with earlier compressible, adiabatic correlations of turbulent transport properties.

Laderman, A. J.↗

The effect of aerodynamic and structural detuning on turbomachine supersonic unstalled torsional flutter

The effects of alternate-blade structural detuning and adjacent-blade alternate-circumferential-spacing aerodynamic detuning on the supersonic unstalled torsional flutter stability of a turbomachine rotor are investigated analytically. An unsteady aerodynamic model employing influence coefficients is constructed for the case of a flat-plate-airfoil cascade in torsion-mode harmonic oscillation in a supersonic inviscid isentropic adiabatic irrotational perfect-gas inlet flow with a subsonic leading-edge locus. The influence coefficients and equations of motion are derived; the model is verified by applying it to the 12-blade cascade-B flow geometry studied by Verdon and McCune (1975); and the results are presented graphically. It is found that the rotor can be stabilized over the entire reduced frequency range by applying a combination of structural and aerodynamic detuning as the passive flutter-control mechanism.

Hoyniak, D.↗

A comparison of adiabatic and diabatic forcing in an intense extratropical cyclone system

The roles of adiabatic and diabatic forcings in the development of a cyclone are examined. The two types of height tendency equations (classical quasi-geostrophic and extended forms) employed to solve for height tendency are described. The extended equation, in which geostrophic wind is replaced with observed wind, the effects of diabatic heating are added, and three-dimensional varying static stability is permitted, is utilized to analyze the evolution of the January 9-11, 1975 cyclone event. The analysis reveals that vorticity advection has a primary effect on the movement and propagation of the wave system; the differential thermal advection has a secondary, but significant role in wave development; latent heating effects are minimal; and the vertical advection of static stability slows the propagation and development of the wave.

Tsou, Chih-Hua↗

Sideslip in a viscous compressible gas

By an analysis of the Navier-Stokes equations it is shown that the aerodynamic coefficients of an infinite rectangular swept wing in an isothermal or adiabatic flow of a compressible gas can be determined from the aerodynamic coefficients of the unswept wing. When the flow is neither isothermal nor adiabatic, a three-dimensional boundary layer theory is developed and applied to the special case of a swept flat plate.

THEORIES - WINGS, SWEEPBACK & SWEEPFORWARD↗

A diagnosis of the development of a winter anticyclone over North America

This paper examines the 48-h life cycle of a winter anticyclone occurring over North America from 18 to 20 January 1979 using Goddard Laboratory for Atmospheres FGGE level 3b (SOP 1) global analyses on a 4 deg latitude by 5 deg longitude grid. Applying the relatively new methodology of the Zwack-Okossi equation, results show that anticyclonic vorticity advection and cold-air advection acted to develop the anticyclone, while adiabatic warming in the descending air opposed development. Other forcing processes made only small contributions to anticyclone changes. Vertical profiles of the development quantities reveal that vorticity and temperature advections, as well as the adiabatic warming, maximized in the 200-300-mb layer.

King, Melinda L.↗

Heat Capacity Anomaly Near the Lower Critical Consolute Point of Triethylamine-Water

The heat capacity of the binary liquid mixture triethylamine-water has been measured near its lower critical consolute point using a scanning, adiabatic calorimeter. Two data runs are analyzed to provide heat capacity and enthalpy data that are fitted by equations with background terms and a critical term that includes correction to scaling. The critical exponent a was determined to be 0.107 +/- 0.006, consistent with theoretical predictions. When alpha was fixed at 0.11 to determine various amplitudes consistently, our values of A(+) and A(-) agreed with a previous heat capacity measurement, but the value of A(-) was inconsistent with values determined by density or refractive index measurements. While our value for the amplitude ratio A(+)/ A(-) = 0.56 +/- 0.02 was consistent with other recent experimental determinations in binary liquid mixtures, it was slightly larger than either theoretical predictions or recent experimental values in liquid-vapor systems. The correction to scaling amplitude ratio D(+)/D(-) = 0.5 +/- 0.1 was half of that predicted. As a result of several more precise theoretical calculations and experimental determinations, the two-scale-factor universality ratio X, which we found to be 0.019 +/- 0.003, now is consistent among experiments and theories. A new 'universal' amplitude ratio R(sup +/-)(sub Bcr) involving the amplitudes for the specific heat was tested. Our determination of R(sup +/-)(sub Bcr) = -0.5 +/- 0.1 and R(sup -)(sub Bcr) = 1.1 +/- 0.1 is smaller in magnitude than predicted and is the first such determination in a binary fluid mixture.

Flewelling, Anne C.↗

A generalized version of the Rankine-Hugoniot relations including ionization, dissociation and related phenomena

For purposes of computing shocks in stellars atmospheres and winds we have developed a generalized version of the Rankine-Hugoniot relations including ionization, dissociation, radiation and related phenomena such as excitation, rotation and vibration of molecules. The new equations are given in analytical form. They are valid as long as the internal energy E, the total pressure P, and the first adiabatic coefficient gamma(sub 1) can be evaluated. However, we have not treated shock structures. In the case of non-LTE we have to employ an approximation for gamma(sub 1) because in that case no definition exists. Our new version of the Rankine-Hugoniot relations can easily be used for many purposes including ab-initio modeling. In our derivation we introduce a parameter gamma(sub H), which is definded as the ratio of the enthalpy H (sometimes called heat function w) to the internal energy E (sometimes called U). Using this parameter we solve the equations for changing mu and (d(natural log P)/d(natural log rho))(sub ad) identically equal to gamma(sub 1) on both sides of the shock. Both gamma(sub H) and gamma(sub 1), and also mu are functions of pressure P and temperature T. We present: (1) the derivation, (2) examples of gamma(sub 1) (P,T) and gamma(sub H) (P,T) which include/exclude ionization and radiation, and (3) as an example the differences in post-shock parameters as function of the pre-shock temperature for the case with ionization and without radiation.

Nieuwenhuijzen, H.↗

Adiabatic Shock Capturing in Perfect Gas Hypersonic Flows

This paper considers the streamline-upwind Petrov/Galerkin (SUPG) method applied to the compressible Euler and Navier-Stokes equations in conservation-variable form. The spatial discretization, including a modified approach for interpolating the inviscid flux terms in the SUPG finite element formulation, is briefly reviewed. Of particular interest is the behavior of the shock capturing operator, which is required to regularize the scheme in the presence of strong, shock-induced gradients. A standard shock capturing operator which has been widely used in previous studies by several authors is presented and discussed. Specific modifications are then made to this standard operator which are designed to produce a more physically consistent discretization in the presence of strong shock waves. The actual implementation of the term in a finite dimensional approximation is also discussed. The behavior of the standard and modified scheme is then compared for several supersonic/hypersonic flows. The modified shock capturing operator is found to preserve enthalpy in the inviscid portion of the flowfield substantially better than the standard operator.

Kirk, Benjamin S.↗

Galactic cosmic ray modulation from 1965-1970.

Numerical solutions of the cosmic-ray equation of transport within the solar cavity and including the effects of diffusion, convection, and energy losses due to adiabatic deceleration, have been used to reproduce the modulation of galactic electrons, protons, and helium nuclei observed during the period from 1965 to 1970. Kinetic energies between 10 and 10,000 MeV/nucleon are considered. Computed and observed spectra are given for the years 1965, 1968, 1969, and 1970 together with the diffusion coefficients. These diffusion coefficients are assumed to be of separable form in rigidity and radial dependence, and are consistent with the available magnetic-field power spectra. The force-field solutions are given for these diffusion coefficients and galactic spectra and are compared with the numerical solutions. It is shown that the energy losses and convection lead to near-earth nuclei spectra at kinetic energies less than or equal to 100 MeV/nucleon in which the differential intensity is proportional to the kinetic energy with little dependence on the form of the galactic spectrum. This dependence is in agreement with the observed spectra of all species of atomic nuclei and it is argued that this provides strong observational evidence for the presence of energy losses in the propagation process, and for the exclusion of low-energy galactic nuclei from near earth.

Urch, I. H.↗

Asymptotic theory of intermediate- and high-degree solar acoustic oscillations

A second-order asymptotic approximation is developed for adiabatic nonradial p-modes of a spherically symmetric star. The exact solutions of adiabatic oscillations are assumed in the outermost layers, where the asymptotic description becomes invalid, which results in a eigenfrequency equation with model-dependent surface phase shift. For lower degree modes, the phase shift is a function of frequency alone; for high-degree modes, its dependence on the degree is explicitly taken into account.

Brodsky, M.↗

The computation of isentropic atmospheric trajectories using a 'discrete model' formulation

An explicit technique for calculating atmospheric trajectories is presented as an alternative method to the standard implicit scheme of Danielsen (1961). The technique uses the inviscid equations of motion and the discrete model formulation derived by Greenspan (1972, 1973) to compute trajectories on isentropic surfaces, assuming adiabatic flow. The discrete model formulation is designed specifically for a Lagrangian system and objectively accounts for the geostrophic departures, local psi-tendencies, and the subsequent accelerations along the entire length of the trajectory. Application of the discrete formulation to a diagnostic case study yielded favorable results.

Petersen, R. A.↗

Exact solutions and low-frequency instability of the adiabatic auroral arc model

The adiabatic auroral arc model couples a kinetic theory parallel current driven by mirror forces to horizontal ionospheric currents; the resulting equations are nonlinear. Some exact stationary solutions to these equations, some of them based on the Liouville equation, are developed, with both latitudinal and longitudinal spatial variations. These Liouville equation exact solutions are related to stability boundaries of low-frequency instabilities such as Kelvin-Helmholtz, as shown by a study of a simplified model.

Cornwall, John M.↗