Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “equilibrium potential”

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 325 records · Page 18

Numerical experiments on the role of radiative processes in the development and maintenance of upper level clouds

A time-dependent, two-dimensional Eulerian model is presented whose purpose is to obtain more realistic parameterizations of extended high level cloudiness, and the results of a numerical experiment using the model are reported. The model is anelastic and the Bousinesque assumption is invoked. Unresolved subgrid scale processes are parameterized as eddy diffusion processes. Two phases of water are incorporated and equilibrium between them is assumed. The effects of infrared radiative processes are parametrically represented. Two simulations were conducted with identical initial conditions; in one of them, the radiation term was never turned on. The mean values of perturbation potential temperature at each level in the domain are plotted versus height after 15, 30, and 60 minutes of simulated time. The influence of the radiative term is seen to impose a cooling trend, leading to an increased generation of ice water and an increased generation of turbulent kinetic energy in the cloud layer.

Starr, D. O'C.↗

The stability of thick, self-gravitating disks in protostellar systems

Using numerical hydrodynamics techniques, we perform a nonlinear stability analysis of accretion disk systems that contain thick, self-gravitating disks. The systems are initially represented by a point mass M(sub c) at the center and a geometrically thick, axisymmetric disk of mass M(sub d) that supports uniform specific angular momentum and obeys an n = 3/2, polytropic equation of state. The equilibrium disk structure is uniquely defined upon the specification of two key dimensionless system parameters: M(sub d)/M(sub c) and T/absolute value of W (the ratio of rotational kinetic energy if the disk to the gravitational potential energy of the system). The focus of this work is on the identification of systems within this two-dimensional parameter space that are marginally unstable toward the development of nonaxisymmetric distortions. The geometric form and relative pattern speed of the disk's distortion as well as the likelihood of disk fragmentation as a result of such instabilities is examined, particularly in the context of protostellar systems. The value of T/absolute value of W at which thick disks first become dynamically unstable to nonaxisymmetric distortions is found to vary significantly with the mass ratio of the accertion disk system. Nonaxisymmetric eigenmodes with four distinctly different characters are identified in systems with mass ratios in the range 0.2 less than or equal to (M(sub d)/M(sub c)) less than or equal to 5.

Woodward, John W.↗

Global Magnetohydrodynamic Modeling of the Solar Corona

Under this contract, we have continued our investigations of the large scale structure of the solar corona and inner heliosphere using global magnetohydrodynamic (MHD) simulations. These computations have also formed the basis for studies of coronal mass ejections (CMES) using realistic coronal configurations. We have developed a technique for computing realistic magnetohydrodynamic (MHD) computations of the solar corona and inner heliosphere. To perform computations that can be compared with specific observations, it is necessary to incorporate solar observations into the boundary conditions. We have used the Wilcox Solar Observatory synoptic maps (collected during a solar rotation by daily measurements of the line-of-sight magnetic field at central meridian) to specify the radial magnetic field (B,) at the photosphere. For the initial condition, we use a potential magnetic field consistent with the specified distribution of B, at the lower boundary, and a wind solution consistent with the specified plasma density and temperature at the solar surface. Together this initial condition forms a (non-equilibrium) approximation of the state of the solar corona for the time-dependent MHD computation. The MHD equations are then integrated in time to steady state. Here we describe solutions relevant to a recent solar eclipse, as well as Ulysses observations. We have also developed a model configuration of solar minimum, useful for studying CME initiation and propagation.

Linker, Jon A.↗

Dynamical and Chemical Behavior of the Lower Stratosphere and Interactions with the Troposphere

Equivalent-barotropic calculations, in tandem with Lagrangian analyses, reveal how changes of total ozone follow from vertical and horizontal transport by planetary waves. Those calculations also throw light on how diabatic motions comprising the Brewer-Dobson circulation develop from quasi-horizontal advection by planetary waves. Potential temperature along a material surface indicates organized subsidence inside the polar-night vortex, resembling tracer observations from UARS. Lagrangian histories illustrate that this sinking motion follows in large part from parcels being driven out of thermodynamic equilibrium by planetary waves, especially at high latitudes. Irreversible heat transfer then produces a net drift of air across isentropic surfaces as parcels orbit about the displaced vortex. By driving mean-meridional overturning in the stratosphere, this downward drift is ultimately responsible for transferring ozone from the tropics to the extratropical lower stratosphere. It also introduces horizontal structure into the distribution of total ozone, which surfaces clearly in ozone trends. High-resolution global cloud imagery constructed from 6 satellites simultaneously observing the Earth was used to investigate the spectrum of equatorial waves generated by tropical convection and propagating vertically into the stratosphere. The results indicate that temperature variability is dominated by planetary-scale equatorial waves like the Kelvin mode, which agrees with satellite observations of the tropical stratosphere. However, the Kelvin mode accounts for only about 30 - 50% of the eastward momentum flux radiating into the stratosphere, the remainder coming from gravity waves. An algorithm was developed to determine 3-dimensional atmospheric motion from satellite tracer measurements. Based on Lagrangian constraints, the algorithm circumvents limitations of the traditional scheme for inferring motion from temperature measurements and determines the circulation in the tropics as reliably as elsewhere. A study of deep convection revealed that the highest towers (those penetrating into stratospheric air and controlling tropopause height and composition through convective mixing) occur in close association with the diurnal cycle of convection. Clouds colder than 220 K develop almost entirely in association with the diurnal cycle of convection over tropical landmasses and substantially in association with it even over maritime regions.

Salby, Murry L.↗

Ion Exchange Technology Development in Support of the Urine Processor Assembly

The urine processor assembly (UPA) on the International Space Station (ISS) recovers water from urine via a vacuum distillation process. The distillation occurs in a rotating distillation assembly (DA) where the urine is heated and subjected to sub-ambient pressure. As water is removed, the original organics, salts, and minerals in the urine become more concentrated and result in urine brine. Eventually, water removal will concentrate the urine brine to super saturation of individual constituents, and precipitation occurs. Under typical UPA DA operating conditions, calcium sulfate or gypsum is the first chemical to precipitate in substantial quantity. During preflight testing with ground urine, the UPA achieved 85% water recovery without precipitation. However, on ISS, it is possible that crewmember urine can be significantly more concentrated relative to urine from ground donors. As a result, gypsum precipitated in the DA when operating at water recovery rates at or near 85%, causing the failure and subsequent re14 NASA Tech Briefs, September 2013 placement of the DA. Later investigations have demonstrated that an excess of calcium and sulfate will cause precipitation at water recovery rates greater than 70%. The source of the excess calcium is likely physiological in nature, via crewmembers' bone loss, while the excess sulfate is primarily due to the sulfuric acid component of the urine pretreatment. To prevent gypsum precipitation in the UPA, the Precipitation Prevention Project (PPP) team has focused on removing the calcium ion from pretreated urine, using ion exchange resins as calcium removal agents. The selectivity and effectiveness of ion exchange resins are determined by such factors as the mobility of the liquid phase through the polymer matrix, the density of functional groups, type of functional groups bound to the matrix, and the chemical characteristics of the liquid phase (pH, oxidation potential, and ionic strength). Previous experience with ion exchange resins has demonstrated that the most effective implementation for an ion exchange resin is a cartridge, or column, in which the resin is contained. Based on the results of equilibrium and sub-scale dynamic column testing, a possible solution for mitigating the calcium precipitation issue on the ISS has been identified. From an original pool of 13 ion exchange resins, two candidates have been identified that demonstrate substantial calcium removal on the sub-scale. The dramatic reduction in resin performance from published calcium uptake demonstrates the need for thorough evaluation of resins at the low pH and strong oxidizing environment present in the UPA. Chemical variations in the influent (calcium concentrations and pretreatment dosing) appear to have a noticeable impact on the calcium capacity of the resin. Low calcium concentrations and high pretreatment dosing will likely result in a decrease in calcium capacity. Conversely, low pre trea t - ment dosing will likely result in an increase in calcium capacity. In contrast, investigations at a variety of flow rates, length-to-diameter ratios, resin volumes, and flow regimes (continuous versus pulsed) show that changes in physical parameters do not have substantial impacts on resin performance in the very low specific velocity ranges of interest. This result is particularly useful because most commercial applications at higher specific velocities do show a relatively strong relationship between flow and capacity. The lack of a strong relationship will allow more flexibility in the implementation of an ion exchange bed for flight. Verification of subscale tests with flight-scale resin beds is recommended prior to implementation in the on-orbit UPA.

Mitchell, Julie↗

Modeling and Simulation of Radiative Compressible Flows in Aerodynamic Heating Arc-Jet Facility

Numerical simulations of an arc heated flow inside NASA's 20 [MW] Aerodynamics heating facility (AHF) are performed in order to investigate the three-dimensional swirling flow and the current distribution inside the wind tunnel. The plasma is considered in Local Thermodynamics Equilibrium(LTE) and is composed of Air-Argon gas mixture. The governing equations are the Navier-Stokes equations that include source terms corresponding to Joule heating and radiative cooling. The former is obtained by solving an electric potential equation, while the latter is calculated using an innovative massively parallel ray-tracing algorithm. The fully coupled system is closed by the thermodynamics relations and transport properties which are obtained from Chapman-Enskog method. A novel strategy was developed in order to enable the flow solver and the radiation calculation to be preformed independently and simultaneously using a different number of processors. Drastic reduction in the computational cost was achieved using this strategy. Details on the numerical methods used for space discretization, time integration and ray-tracing algorithm will be presented. The effect of the radiative cooling on the dynamics of the flow will be investigated. The complete set of equations were implemented within the COOLFluiD Framework. Fig. 1 shows the geometry of the Anode and part of the constrictor of the Aerodynamics heating facility (AHF). Fig. 2 shows the velocity field distribution along (x-y) plane and the streamline in (z-y) plane.

Pasma flows↗

Degeneracy effects of neutrino mass ejection in supernovae

A neutrino mechanism is discussed in order to explain supernovae in massive stars. An argument is presented for supernova mass ejection through leptonic neutrino transport characteristics suppressed by the arbitrary zero chemical potential condition. Results show that lepton conservation effects may be important in supernova neutrino transport. At low temperature and density the diffusion approximation becomes less precise because of the long mean free paths of low energy neutrinos. The amount of equilibrium neutrino spectrum affected here is small over most of the collapsing supernova structure.

Mazurek, T. J.↗

Adsorption and excess fission xenon

The adsorption of Xe and Kr on lunar soil 10084 was measured by a method that employs only very low fractions of monolayer coverage. Results are presented as parameters for calculation of the Henry constant for adsorption as a function of temperature. The adsorption potentials are about 3 kcal/mole for Kr and 5 kcal/mole for Xe; heating the sample in vacuum increased the Xe potential to nearly 7 kcal/mole. Henry constants at the characteristic lunar temperature are about 0.3 cu cm STP/g-atm. These data were applied to consider whether adsorption is important in producing the excess fission Xe effect characteristic of highland breccias. Sorption equilibrium with a transient lunar atmosphere vented fission Xe produces concentrations seven orders of magnitude lower than observed concentrations. Higher concentrations result because of the resistance of the regolith to upward diffusion of Xe. A diffusion coefficient of 0.26 sq cm/sec is estimated for this process.

Podosek, F. A.↗

Complex time dependent wave packet technique for thermal equilibrium systems - Electronic spectra

A time dependent wave packet method is presented for the rapid calculation of the properties of systems in thermal equilibrium and is applied, as an illustration, to electronic spectra. The thawed Gaussian approximation to quantum wave packet dynamics combined with evaluation of the density matrix operator by imaginary time propagation is shown to give exact electronic spectra for harmonic potentials and excellent results for both a Morse potential and for the band contours of the three transitions of the visible electronic absorption spectrum of the iodine molecule. The method, in principle, can be extended to many atoms (e.g., condensed phases) and to other properties (e.g., infrared and Raman spectra and thermodynamic variables).

Reimers, J. R.↗

Numerical simulations of non-equilibrium shock layers with efficient implicit schemes

Current and future calculations of nonequilibrium shock layers require the use of a very large number of equations, due to a multiplicity of chemical species, excited states, and internal energy modes. The computational cost associated with the use of standard implicit methods becomes prohibitive; it is therefore desirable to examine the potential of several methods and determine if any can be projected to be more efficient and accurate for large systems of equations. Here, the performance of several implicit schemes on several simple practical examples of reacting flows is examined. The Euler equations are solved by three different implicit methods, and two methods of coupling between the fluid dynamics and the chemistry are studied. Several cases of stiffness are considered and both 1D and 2D examples are computed.

Cambier, Jean-Luc↗

Dynamics of Galaxy Clusters and Expectations from Astro-H

Galaxy clusters span a range of dynamical states, from violent mergers -- the most energetic events in the Universe -- to systems near hydrostatic equilibrium that allow us to map their dark matter distribution using X-ray observations of the intracluster gas. Accurate knowledge of the cluster physics, and in particular, the physics of the hot intracluster gas, is required to realize the full potential of clusters as cosmological probes. So far, we have been studying the cluster dynamics indirectly, deducing merger geometries, cluster masses, etc., using X-ray brightness and gas temperature mapping. For the first time, the calorimeter onboard Astro-H will provide direct measurements of line-of-sight velocities and turbulent broadening in the intracluster gas, testing many of our key assumptions about clusters. This talk will summarize expectations for cluster dynamic studies with this new instrument.

Markevitch, Maxim↗

Thermal Conductivity and Specific Heat Measurements of Candidate Structural Materials for the JWST Optical Bench

The James Webb Space Telescope will include an optical bench known as the integrated science instrument module (ISIM). Candidate structural materials for the ISIM must have low density, high stiffness, high thermal conductivity, and low thermal expansion coefficient at the operating temperature of 30 Kelvin. The specific heat is also important in modeling the on-orbit cooldown. We built two different systems for measuring the thermal conductivity and specific heat of samples between 4 Kelvin and 290 Kelvin. Both experiments were carefully designed to minimize potential errors due to radiative heat transfer. We chose the cooling system and instrumentation to allow long-term unattended operation. Software was developed to automate each experiment. It used an algorithm designed to ensure that each system was in thermal equilibrium before a measurement was taken. We describe the two experiments and present the data.

Canavan, E. R.↗

Direct Numerical Simulation of a Weakly Stratified Turbulent Wake

Direct numerical simulation (DNS) is used to investigate a time-dependent turbulent wake evolving in a stably stratified background. A large initial Froude number is chosen to allow the wake to become fully turbulent and axisymmetric before stratification affects the spreading rate of the mean defect. The uncertainty introduced by the finite sample size associated with gathering statistics from a simulation of a time-dependent flow is reduced, compared to earlier simulations of this flow. The DNS reveals the buoyancy-induced changes to the turbulence structure, as well as to the mean-defect history and the terms in the mean-momentum and turbulence-kinetic-energy budgets, that characterize the various states of this flow - namely the three-dimensional (essentially unstratified), non-equilibrium (or 'wake-collapse') and quasi-two-dimensional (or 'two-component') regimes observed elsewhere for wakes embedded in both weakly and strongly stratified backgrounds. The wake-collapse regime is not accompanied by transfer (or 'reconversion') of the potential energy of the turbulence to the kinetic energy of the turbulence, implying that this is not an essential feature of stratified-wake dynamics. The dependence upon Reynolds number of the duration of the wake-collapse period is demonstrated, and the effect of the details of the initial/near-field conditions of the wake on its subsequent development is examined.

Redford, J. A.↗

Calculation of State Specific Rate Coefficients for Non-Equilibrium Hypersonics Applications: from H(Psi) = E(Psi) to k(T) = A *exp(-E(sub a)/RT)

Development of High-Fidelity Physics-Based Models to describe hypersonic flight through the atmospheres of Earth and Mars is underway at NASA Ames Research Center. The goal is to construct chemistry models of the collisional and radiative processes that occur in the bow shock and boundary layers of spacecraft during atmospheric entry that are free of empiricism. In this talk I will discuss our philosophy and describe some of our progress. Topics to be covered include thermochemistry, internal energy relaxation, collisional dissociation and radiative emission and absorption. For this work we start by solving the Schrodinger equation to obtain accurate interaction potentials and radiative properties. Then we invoke classical mechanics to compute state-specific heavy particle collision cross sections and reaction rate coefficients. Finally, phenomenological rate coefficients and relaxation times are determined from master equation solutions.

aerothermodynamics↗

Potentials of fissioning plasmas

Successful experiments with the nuclear pumping of lasers have demonstrated that in a gaseous medium the kinetic energy of fission fragments can be converted directly into nonequilibrium optical radiation. This confirms the concept that the fissioning medium in a gas-phase nuclear reactor shows an internal structure such as a plasma in near thermal equilibrium varying up to a state of extreme nonequilibrium. During 20 years of research under NASA support major elements of the fissioning plasma reactor were demonstrated in theory and experiment, culminating in a proof-of-principle reactor test conducted at the Los Alamos Scientific Laboratory. It is concluded that the construction of a gaseous fuel reactor power plant is within the reach of present technology.

Thom, K.↗

Plasma Chemical Conversion and Resource Generation Beyond Low-Earth Orbit

As humanity gears up for its return to the moon after more than half a century, collaborative efforts between NASA, Artemis Accords Partners, and private industry are underway to establish the necessary infrastructure and technologies for lunar habitation and eventual Mars exploration. However, the traditional ISS resupply and waste management model is impractical and economically infeasible for prolonged missions to the Moon and Mars. Advanced chemical conversion technologies are needed to generate vital consumable products from local planetary resources (ISRU) and recycled gasses and waste within semi-closed loop life-support systems. Low-temperature plasma reactors are emerging power-to-gas technologies with the potential to facilitate various chemical synthesis processes with hardware commonality and redundancy. In plasma-based systems, electrical power is used to ionize a feedstock gas, creating a highly reactive environment that leverages electron excitation chemistry to break stable molecular bonds and form value added products. Unlike thermal chemical processes, plasma reactors operate at non-equilibrium conditions, allowing for lower-temperature operation and instantaneous start-up, making them adaptable to intermittent power availability. Moreover, their scalability permits deployment in both portable astronaut systems and large-scale industrial setups for colonies. One promising application of plasmas is for CO 2 conversion. Carbon dioxide comprises 96% of the Martian atmosphere and is a byproduct of human respiration, which typically must be scrubbed and vented from space habitats. A plasma source integrated with membrane separation technology could generate a stream of oxygen for life support and rocket propellant. CO 2 splitting may also be beneficial as a precursor to manufacture carbon-based products and fuels in situ, like methane, methanol, and polyethylene. Plasma-assisted CO 2 conversion is a simpler case to study without the concern of selectivity and is the first step toward complex chemical synthesis. This work presents preliminary experimental case study from a plasma reactor for CO 2 conversion and casts a vision for the potential of plasma technologies in a sub-architecture for resource production to enable the next generation of human spaceflight activities.

Plasma↗

Static current-sheet models of quiescent prominences

A particular class of theoretical models idealize the prominence to be a discrete flat electric-current sheet suspended vertically in a potential magnetic field. The weight of the prominence is supported by the Lorentz force in the current sheet. These models can be extended to have curved electric-current sheets and to vary three-dimensionally. The equation for force balance is 1 over 4 pi (del times B) times Bdel p- p9 z=zero. Using Cartesian coordinates we take, for simplicity, a uniform gravity with constant acceleration g in the direction -z. If we are interested not in the detailed internal structure of the prominence, but in the global magnetic configuration around the prominence, we may take prominence plasma to be cold. Consideration is given to how such equilibrium states can be constructed. To simplify the mathematical problem, suppose there is no electric current in the atmosphere except for the discrete currents in the cold prominence sheet. Let us take the plane z =0 to be the base of the atmosphere and restrict our attention to the domain z greater than 0. The task we have is to solve for a magnetic field which is everywhere potential except on some free surface S, subject to suit able to boundary conditions. The surface S is determined by requiring that it possesses a discrete electric current density such that the Lorentz force on it is everywhere vertically upward to balance the weight of the material m(S). Since the magnetic field is potential in the external atmosphere, the latter is decoupled from the magnetic field and its plane parallel hydrostatic pressure and density can be prescribed.

Wu, F.↗

The thermal environment of transatmospheric vehicles

A transatmospheric vehicle using primarily air-breathing propulsion must fly in the denser part of the atmosphere to achieve adequate acceleration to reach orbital speed. The potentially long ascent times, combined with the need for a low-drag configuration, result in a severe aerothermodynamics environment. To achieve low drag, the vehicle must have a relatively sharp nose and wing leading edges. The ascent peak stagnation point and equilibrium wall temperatures for the wing leading edge can reach values of 4000 K and 3000 K, respectively, for high dynamic pressure trajectories, making some form of mass addition cooling mandatory. The corresponding temperatures during entry are about 1500 K lower. The vehicle windward centerline temperatures are more moderate, however, with values peaking around 1500 K. Therefore, radiative cooling should be effective over large areas of the vehicle. The windward, centerline heat loads are relatively insensitive to the dynamic pressure of the ascent trajectory, in contrast to the stagnation point and wing leading edge. The windward surface entry heat loads are much lower, but depend strongly on the flightpath.

Tauber, Michael E.↗