Search NASA⌕ Search

SEARCH · Search NASA

Results for “energy 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 271 records · Page 15

Modeling the effects of latitudinal gradients in stellar winds, with application to the solar wind

A steady, axisymmetric, quasi-radial, global model previously developed for stellar winds with embedded magnetic fields has been extended to include latitudinal gradient effects on the azimuthal velocity and magnetic field. The linear results at large radii are presented for large-amplitude latitudinal variations in the radial magnetic field, mass loss rate, and radial velocity of the wind. The magnetohydrodynamic (MHD) equations predict meridional flows that develop naturally from internal magnetic stresses. The flows open flux tubes in the star's equatorial plane, redistributing mass and magnetic flux as a function of stellar latitude. The plasma spins up to conserve angular momentum in fields and plasma. The results are generally applicable to stellar winds (including radiatively driven winds), provided that the internal structure is not dominated by rotation. The asymptotic solutions do not explicitly depend on the form of the energy equation, although the assumed O(1) state which drives these solutions depends on the deposition of energy and momentum throughout the wind.

Nerney, S.↗

Finite-element solutions for geothermal systems

Vector potential and scalar potential are used to formulate the governing equations for a single-component and single-phase geothermal system. By assuming an initial temperature field, the fluid velocity can be determined which, in turn, is used to calculate the convective heat transfer. The energy equation is then solved by considering convected heat as a distributed source. Using the resulting temperature to compute new source terms, the final results are obtained by iterations of the procedure. Finite-element methods are proposed for modeling of realistic geothermal systems; the advantages of such methods are discussed. The developed methodology is then applied to a sample problem. Favorable agreement is obtained by comparisons with a previous study.

Chen, J. C.↗

Radiative transfer in spherical atmospheres

A method for defining spherical model atmospheres in radiative/convective and hydrostatic equilibrium is presented. A finite difference form is found for the transfer equation and a matrix operator is developed as the discrete space analog (in curvilinear coordinates) of a formal integral in plane geometry. Pressure is treated as a function of temperature. Flux conservation is maintained within the energy equation, although the correct luminosity transport must be assigned for any given level of the atmosphere. A perturbed integral operator is used in a complete linearization of the transfer and constraint equations. Finally, techniques for generating stable solutions in economical computer time are discussed.

Kalkofen, W.↗

A model of concurrent flow flame spread over a thin solid fuel

A numerical model is developed to examine laminar flame spread and extinction over a thin solid fuel in lowspeed concurrent flows. The model provides a more precise fluid-mechanical description of the flame by incorporating an elliptic treatment of the upstream flame stabilization zone near the fuel burnout point. Parabolic equations are used to treat the downstream flame, which has a higher flow Reynolds number. The parabolic and elliptic regions are coupled smoothly by an appropriate matching of boundary conditions. The solid phase consists of an energy equation with surface radiative loss and a surface pyrolysis relation. Steady spread with constant flame and pyrolysis lengths is found possible for thin fuels and this facilitates the adoption of a moving coordinate system attached to the flame with the flame spread rate being an eigen value. Calculations are performed in purely forced flow in a range of velocities which are lower than those induced in a normal gravity buoyant environment. Both quenching and blowoff extinction are observed. The results show that as flow velocity or oxygen percentage is reduced, the flame spread rate, the pyrolysis length, and the flame length all decrease, as expected. The flame standoff distance from the solid and the reaction zone thickness, however, first increase with decreasing flow velocity, but eventually decrease very near the quenching extinction limit. The short, diffuse flames observed at low flow velocities and oxygen levels are consistent with available experimental data. The maximum flame temperature decreases slowly at first as flow velocity is reduced, then falls more steeply close to the quenching extinction limit. Low velocity quenching occurs as a result of heat loss. At low velocities, surface radiative loss becomes a significant fraction of the total combustion heat release. In addition, the shorter flame length causes an increase in the fraction of conduction downstream compared to conduction to the fuel. These heat losses lead to lower flame temperatures, and ultimately, extinction. This extinction mechanism differs from that of blowoff, where the flame is unable to be stabilized due to the high flow velocity.

Ferkul, Paul V.↗

Distributed Relaxation Multigrid and Defect Correction Applied to the Compressible Navier-Stokes Equations

The distributed-relaxation multigrid and defect- correction methods are applied to the two- dimensional compressible Navier-Stokes equations. The formulation is intended for high Reynolds number applications and several applications are made at a laminar Reynolds number of 10,000. A staggered- grid arrangement of variables is used; the coupled pressure and internal energy equations are solved together with multigrid, requiring a block 2x2 matrix solution. Textbook multigrid efficiencies are attained for incompressible and slightly compressible simulations of the boundary layer on a flat plate. Textbook efficiencies are obtained for compressible simulations up to Mach numbers of 0.7 for a viscous wake simulation.

Thomas, J. L.↗

Semi-classical approximation description of static properties of nuclei

In this report, we present recent development in semi-classical description of static properties of nuclei. By employing the Wigner transform, we derive simple semiclassical approximations for evaluating properties of nuclear system. We present results of calculations of static properties, using the energy density functional associated with Skyrme effective nucleon–nucleon interaction. In particular, we consider properties of the Wigner distribution function, line of beta stability, Coulomb and symmetry energies, deformation energies, equation of state, nuclear radii, and single-particle level density.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

One-dimensional hybrid satellite track model for the Dynamics Explorer 2 (DE 2) satellite

A one-dimensional hybrid satellite track model has been developed to calculate the high-latitude thermospheric/ionospheric structure below the satellite altitude using Dynamics Explorer 2 (DE 2) satellite measurements and theory. This model is based on Emery et al. (1985) satellite track code but also includes elements of Roble et al. (1987b) global mean thermosphere/ionosphere model. A number of parameterizations and data handling techniques are used to input satellite data from several DE 2 instruments into this model. Profiles of neutral atmospheric densities are determined from the Mass Spectrometer Incoherent Scatter 1990 (MSIS-90) model and measured neutral temperatures. Measured electron precipitation spectra are used in an auroral model to calculate particle impact ionization rates below the satellite. These rates are combined with a solar ionization rate profile and used to solve the O(+) diffusion equation, with the measured electron density as an upper boundary condition. The calculated O(+) density distribution, as well as the ionization profiles, are then used in a photochemical equilibrium model to calculate the electron and molecular ion densities. The electron temperature is also calculated by solving the electron energy equation with an upper boundary condition determined by the DE 2 measurement. The model enables calculations of altitude profiles of conductivity and Joule heating rate along and below the satellite track. In a first application of the new model, a study is made of thermospheric and ionospheric structure below the DE 2 satellite for a single orbit which occurred on October 25, 1981. The field-aligned Poynting flux, which is independently obtained for this orbit, is compared with the model predictions of the height-integrated energy conversion rate. Good quantitative agreement between these two estimates has been reached. In addition, measurements taken at the incoherent scatter radar site at Chatanika (65.1 deg N, 147.4 deg W) during a DE 2 overflight are compared with the model calculations. A good agreement was found in lower thermospheric conductivities and Joule heating rate.

Deng, Wei↗

One-Dimensional Hybrid Satellite Track Model for the Dynamics Explorer 2 (DE 2) Satellite

A one-dimensional hybrid satellite track model has been developed to calculate the high-latitude thermospheric/ionospheric structure below the satellite altitude using Dynamics Explorer 2 (DE 2) satellite measurements and theory. This model is based on Emery et al. satellite track code but also includes elements of Roble et al. global mean thermosphere/ionosphere model. A number of parameterizations and data handling techniques are used to input satellite data from several DE 2 instruments into this model. Profiles of neutral atmospheric densities are determined from the MSIS-90 model and measured neutral temperatures. Measured electron precipitation spectra are used in an auroral model to calculate particle impact ionization rates below the satellite. These rates are combined with a solar ionization rate profile and used to solve the O(+) diffusion equation, with the measured electron density as an upper boundary condition. The calculated O(+) density distribution, as well as the ionization profiles, are then used in a photochemical equilibrium model to calculate the electron and molecular ion densities. The electron temperature is also calculated by solving the electron energy equation with an upper boundary condition determined by the DE 2 measurement. The model enables calculations of altitude profiles of conductivity and Joule beating rate along and below the satellite track. In a first application of the new model, a study is made of thermospheric and ionospheric structure below the DE 2 satellite for a single orbit which occurred on October 25, 1981. The field-aligned Poynting flux, which is independently obtained for this orbit, is compared with the model predictions of the height-integrated energy conversion rate. Good quantitative agreement between these two estimates has been reached. In addition, measurements taken at the incoherent scatter radar site at Chatanika (65.1 deg N, 147.4 deg W) during a DE 2 overflight are compared with the model calculations. A good agreement was found in lower thermospheric conductivities and Joule heating rate.

Deng, Wei↗

The physical origin of dark energy constraints from rubin observatory and CMB-S4 lensing tomography

ABSTRACT We seek to clarify the origin of constraints on the dark energy equation of state parameter from CMB lensing tomography, that is the combination of galaxy clustering and the cross-correlation of galaxies with CMB lensing in a number of redshift bins. We focus on the analytic understanding of the origin of the constraints. Dark energy information in these data arises from the influence of three primary relationships: distance as a function of redshift (geometry), the amplitude of the power spectrum as a function of redshift (growth), and the power spectrum as a function of wavenumber (shape). We find that the effects from geometry and growth play a significant role and partially cancel each other out, while the shape effect is unimportant. We also show that Dark Energy Task Force figure of merit forecasts from the combination of LSST galaxies and CMB-S4 lensing are comparable to the forecasts from cosmic shear in the absence of the CMB lensing map, thus providing an important independent check. Compared to the forecasts with the LSST galaxies alone, combining CMB lensing and LSST clustering information increases the FoM by roughly a factor of 3–4 in the optimistic scenario where systematics are fully under control. We caution that achieving these forecasts will likely require a full analysis of higher-order biasing, photometric redshift uncertainties, and stringent control of other systematic limitations, which are outside the scope of this work, whose primary purpose is to elucidate the physical origin of the constraints.

79 ASTRONOMY AND ASTROPHYSICS↗

Near-wall modeling of compressible turbulent flow

A near-wall two-equation model for compressible flows is proposed. The model is formulated by relaxing the assumption of dynamic field similarity between compressible and incompressible flows. A postulate is made to justify the extension of incompressible models to ammount for compressibility effects. This requires formulation the turbulent kinetic energy equation in a form similar to its incompressible counterpart. As a result, the compressible dissipation function has to be split into a solenoidal part, which is not sensitive to changes of compressibility indicators, and a dilatational part, which is directly affected by these changes. A model with an explicit dependence on the turbulent Mach number is proposed for the dilatational dissipation rate.

So, Ronald M. C.↗

Analysis of laminar flow between stationary and rotating disks with inflow

The laminar flow between a rotating and a stationary disk with inflow was analyzed. Solutions to the dimensionless governing equations are sought by expanding each of the velocity components in powers of inverse radius. The equations to leading order are those for the configuration with no inflow. The subsequent orders yield sets of linear ordinary differential equations. Solutions are obtained for the first two of these subsequent orders. The solutions indicate that inflow tends to increase the magnitude of the azimuthal velocity in the flow between the two disks and to decrease the torque on the rotating disk. For Prandtl number one, an energy integral is obtained which relates the temperature distribution to the velocity distribution for all Reynolds numbers and therefore eliminates the needs for separate solution of the energy equation.

Rohatgi, U.↗

A compressible solution of the Navier-Stokes equations for turbulent flow about an airfoil

A compressible time dependent solution of the Navier-Stokes equations including a transition turbulence model is obtained for the isolated airfoil flow field problem. The equations are solved by a consistently split linearized block implicit scheme. A nonorthogonal body-fitted coordinate system is used which has maximum resolution near the airfoil surface and in the region of the airfoil leading edge. The transition turbulence model is based upon the turbulence kinetic energy equation and predicts regions of laminar, transitional, and turbulent flow. Mean flow field and turbulence field results are presented for an NACA 0012 airfoil at zero and nonzero incidence angles of Reynolds number up to one million and low subsonic Mach numbers.

Shamroth, S. J.↗

Analysis and design of an ultrahigh temperature hydrogen-fueled MHD generator

A coupled gas dynamics/radiative heat transfer analysis of partially ionized hydrogen, in local thermodynamic equilibrium, flowing through an ultrahigh temperature (10,000-20,000 K) magnetohydrodynamic (MHD) generator is performed. Gas dynamics are modeled by a set of quasi-one-dimensional, nonlinear differential equations which account for friction, convective and radiative heat transfer, and the interaction between the ionized gas and applied magnetic field. Radiative heat transfer is modeled using nongray, absorbing-emitting 2D and 3D P-1 approximations which permit an arbitrary variation of the spectral absorption coefficient with frequency. Gas dynamics and radiative heat transfer are coupled through the energy equation and through the temperature- and density-dependent absorption coefficient. The resulting nonlinear elliptic problem is solved by iterative methods. Design of such MHD generators as onboard, open-cycle, electric power supplies for a particular advanced airbreathing propulsion concept produced an efficient and compact 128-MWe generator characterized by an extraction ratio of 35.5 percent, a power density of 10,500 MWe/cu m, and a specific (extracted) energy of 324 MJe/kg of hydrogen. The maximum wall heat flux and total wall heat load were 453 MW/sq m and 62 MW, respectively.

Moder, Jeffrey P.↗

Testing cosmology with double source lensing

Double source lensing provides a dimensionless ratio of distance ratios, a "remote viewing" of cosmology through distances relative to the gravitational lens, beyond the observer. We use this to test the cosmological framework, particularly with respect to spatial curvature and the distance duality relation. We derive a consistency equation for constant spatial curvature, allowing not only the investigation of flat vs curved but of the Friedmann-Lemaître-Robertson-Walker framework itself. For distance duality, we demonstrate that the evolution of the lens mass profile slope must be controlled to ≳ 5 times tighter fractional precision than a claimed distance duality violation. Furthermore, using LensPop forecasts of double source lensing systems in Euclid and LSST surveys we also explore constraints on dark energy equation of state parameters and any evolution of the lens mass profile slope.

79 ASTRONOMY AND ASTROPHYSICS↗

The prediction of the turbulent flow field about an isolated airfoil

A compressible time-dependent solution of the Navier-Stokes equations including a transition-turbulence model is obtained for the isolated airfoil flow field problem. The equations are solved by a consistently split linearized block implicit scheme due to Briley and McDonald. A nonorthogonal body fitted coordinate system is used which has maximum resolution near the airfoil surface and in the region of the airfoil leading edge. The transition-turbulence model is based upon the turbulence kinetic energy equation and predicts regions of laminar, transitional and turbulent flow. Mean flow field and turbulence field results are presented for an NACA 0012 airfoil at zero and nonzero incidence angles at Reynolds number up to one million and low subsonic Mach numbers.

Shamroth, S. J.↗

An analysis of temperature effect in a finite journal bearing with spatial tilt and viscous dissipation

The analysis presented herein deals with the evaluation of the pressure, velocity, and temperature profiles in a finite-length plane journal bearing. The geometry of the case under study consists of a spatially tilted shaft. The two-dimensional Reynolds equation accounts for the variation of the clearance gap h with x and z and is used to model the pressure field. The latter is solved for a variety of shaft tilt angles and then used to calculate the two-dimensional flow field. Finally, the flow field is used in the energy equation to solve for the film temperature profile, when the effect of viscous dissipation is taken into account.

Braun, M. J.↗

Minimal upper mantle temperature variations consistent with observed heat flow and plate velocities

Heat-flow and plate-velocity measurements are used to model upper-mantle temperature variations; the results are intended to form the basis for a gravity-field-variation model and are also applicable to the interpretation of petrological and seismological data. A 5-deg grid is used, a depth of 280 km is chosen as the fully convecting level, and it is assumed that 85 percent of the global heat production (4.0 x 10 to the 13th W) comes from below this level. The velocity field is calculated by integrating the momentum equations down to 280 km (assuming laterally homogeneous viscosity and density) and then used to determine the temperature fields. The results are presented in graphs, tables, and maps. The largest lateral temperature differences are about 1500 C and occur in the top 20 km, where the largest (about 200 C/Myr) nonlinear terms of the energy equation are also found. Below 50 km, heat transfer becomes more convective than conductive, and the most significant temperature variation appears in the form of negative 'tongues' as cold as -825 C at about 100 km. Temperature variations of at least + or - 180 C are calculated at the fully convecting level.

Kaula, W. M.↗

Studying Dark Energy, Black Holes and Cosmic Feedback at X-ray Wavelengths: NASA's Constellation-X Mission

Among the most important topics in modern astrophysics are the nature of the dark energy equation of state, the formation and evolution of supermassive black holes in concert with galaxy bulges, and the self-regulating symmetry imposed by both stellar and AGN feedback. All of these topics are readily addressed with observations at X-ray wavelengths. For instance, theoretical models predict that the majority (98%) of the energy and metal content in starburst superwinds exists in the hot million-degree gas. The Constellation-X observatory is being developed to perform spatially resolved high-resolution X-ray spectroscopy so that we may directly measure the absolute element abundances and velocities of this hot gas. This talk focuses on the driving science behind this mission, which is one of two flagship missions in NASA's Beyond Einstein program. A general overview of the observatory's capabilities and basic technology will also be given.

Hornschemeier, A.↗