Search NASA⌕ Search

SEARCH · Search NASA

Results for “equilibrium constant”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Recent Advancements in the PATO Material Response Code

Introduction: Predicting the complicated multiphysics phenomena during atmospheric entry requires high-fidelity modeling tools to refine estimates of mission risks during entry. To this end, new capabilities are being added to the Porous-material Analysis Toolbox based on OpenFOAM (PATO). PATO is an open-source software for Computational Material Response (CMR) of reactive porous materials submitted to high-temperature environments. The objective of this work is to highlight current efforts to add to and improve upon the modeling capabilities of PATO. These include efforts to loosely couple PATO with other discipline specialized codes including hypersonic Computational Fluid Dynamics (CFD), to assess the interaction effects between pyrolysis gas blowing and the boundary layer, and Computational Solid Mechanics (CSM), to address modeling of mechanical erosion. Other refinements include surface phenomena modeling capabilities to address the effects of silicone-based coatings applied to the TPS during flight preparation, and a unified multiphase solver for a mixed porous-material and plain-fluid domain. Coupling CMR with CFD (CMR/CFD): A loose coupling between PATO and the Data Parallel Line Relaxation (DPLR) CFD code has been achieved by making use of a blowing boundary condition at the heatshield surface available in DPLR. Starting with heat flux estimates with no pyrolysis gas blowing at the surface, blowing gases are computed by the CMR and passed to the CFD such that aerothermal properties of the environment can be recomputed for a new CMR computation. This leads to an iterative process which is supplemented with an estimate of the radiative heat flux using the Nonequilibrium air radiation (NEQAIR) program. The entire iterative process is illustrated in Figure 1. This coupling strategy has been utilized in computing the MSL material response. The goal is to compare the coupled CMR/CFD results with material response results obtained using traditional blowing corrections. Coupling CMS with CMR: A mechanical erosion model is currently being implemented in PATO to account for the additional mass removal induced by high shear conditions. The modeling process at each timestep consists of updating the mechanical properties as a function of temperature and computing the stress tensor and displacement fields of the material. Then, a failure criteria model determines the regions in which the stress exceeds the ultimate strength values resulting in mesh movement to account for mass removal. This model allows the material response simulation to compute the recession due to both oxidation and shear-induced erosion. The model is demonstrated by computing material response of sphere-cone arc jet samples. Surface Modeling Capabilities: NuSil, a silicone-based coating, was sprayed onto the MSL and Mars 2020 heatshields to mitigate shedding of phenolic dust. To better understand the effects of the NuSil coating on the material response, a novel model has been implemented in PATO. In this model, the equilibrium of the charred NuSil surface is modeled as pure silica, and a constant offset, inspired by the classical spallation model, is added to the the char blowing rate and wall enthalpy to reproduce HyMETS experimental results. The model has also been used to estimate the 3D material response of the MSL heatshield. Unified Solver: In addition to the iterative loose coupling approach mentioned above, a multiphase unified solver is being developed to couple the environment (plain-fluid phase) and the porous-material phase. The solver is based on the volume averaged conservation of mass, momentum, and energy for the macroscale with closure models which include microscale effects through effective physicochemical properties. The unified solver has been used to compute flow through a porous plug and solve the Beavers and Joseph problem. Since the strong coupling between phases is inherent to this solver, modeling assumptions present in other coupling methods of material response are mitigated. This strategy also makes it feasible to capture the competition between surface and volume ablation in the same computational domain, which is usually not possible with other coupling approaches.

Material Response↗
NASA NTRSDate not supplied

Recent Advancements in the PATO Material Response Code

Introduction: Predicting the complicated multiphysics phenomena during atmospheric entry requires high-fidelity modeling tools to refine estimates of mission risks during entry. To this end, new capabilities are being added to the Porous-material Analysis Toolbox based on OpenFOAM (PATO) [1,2,3]. PATO is an open-source software for Computational Material Response (CMR) of reactive porous materials submitted to high-temperature environments. The objective of this work is to highlight current efforts to add to and improve upon the modeling capabilities of PATO. These include efforts to loosely couple PATO with other discipline specialized codes including hypersonic Computational Fluid Dynamics (CFD), to assess the interaction effects between pyrolysis gas blowing and the boundary layer, and Computational Solid Mechanics (CSM), to address modeling of mechanical erosion. Other refinements include surface phenomena modeling capabilities to address the effects of silicone-based coatings applied to the TPS during flight preparation, and a unified multiphase solver for a mixed porous-material and plain-fluid domain. Coupling CMR with CFD (CMR/CFD): A loose coupling between PATO and the Data Parallel Line Relaxation (DPLR) [4] CFD code has been achieved by making use of a blowing boundary condition at the heatshield surface available in DPLR. Starting with heat flux estimates with no pyrolysis gas blowing at the surface, blowing gases are computed by the CMR and passed to the CFD such that aerothermal properties of the environment can be recomputed for a new CMR computation. This leads to an iterative process which is supplemented with an estimate of the radiative heat flux using the Nonequilibrium air radiation (NEQAIR) [5] program. The entire iterative process is illustrated in Figure 1. This coupling strategy has been utilized in computing the MSL material response. The goal is to compare the coupled CMR/CFD results with material response results obtained using traditional blowing corrections. Coupling CMS with CMR: A mechanical erosion model is currently being implemented in PATO to account for the additional mass removal induced by high shear conditions. The modeling process at each timestep consists of updating the mechanical properties as a function of temperature and computing the stress tensor and displacement fields of the material. Then, a failure criteria model determines the regions in which the stress exceeds the ultimate strength values resulting in mesh movement to account for mass removal. This model allows the material response simulation to compute the recession due to both oxidation and shear-induced erosion. The model is demonstrated by computing material response of sphere-cone arc jet samples. Surface Modeling Capabilities: NuSil, a silicone-based coating, was sprayed onto the MSL and Mars 2020 heatshields to mitigate shedding of phenolic dust. To better understand the effects of the NuSil coating on the material response, a novel model has been implemented in PATO. In this model, the equilibrium of the charred NuSil surface is modeled as pure silica, and a constant offset, inspired by the classical spallation model, is added to the the char blowing rate and wall enthalpy to reproduce HyMETS experimental results. The model has also been used to estimate the 3D material response of the MSL heatshield [6]. Unified Solver: In addition to the iterative loose coupling approach mentioned above, a multiphase unified solver is being developed to couple the environment (plain-fluid phase) and the porous-material phase. The solver is based on the volume averaged conservation of mass, momentum, and energy for the macroscale with closure models which include microscale effects through effective physicochemical properties. The unified solver has been used to compute flow through a porous plug and solve the Beavers and Joseph problem [7]. Since the strong coupling between phases is inherent to this solver, modeling assumptions present in other coupling methods of material response are mitigated. This strategy also makes it feasible to capture the competition between surface and volume ablation in the same computational domain, which is usually not possible with other coupling approaches.

Thermal Protection Systems↗
NASA NTRS2003

Deviation of Long-Period Tides from Equilibrium: Kinematics and Geostrophy

New empirical estimates of the long-period fortnightly (Mf) tide obtained from TOPEX/Poseidon (T/P) altimeter data confirm significant basin-scale deviations from equilibrium. Elevations in the low-latitude Pacific have reduced amplitude and lag those in the Atlantic by 30 deg or more. These interbasin amplitude and phase variations are robust features that are reproduced by numerical solutions of the shallow-water equations, even for a constant-depth ocean with schematic interconnected rectangular basins. A simplified analytical model for cooscillating connected basins also reproduces the principal features observed in the empirical solutions. This simple model is largely kinematic. Zonally averaged elevations within a simple closed basin would be nearly in equilibrium with the gravitational potential, except for a constant offset required to conserve mass. With connected basins these offsets are mostly eliminated by interbasin mass flux. Because of rotation, this flux occurs mostly in a narrow boundary layer across the mouth and at the western edge of each basin, and geostrophic balance in this zone supports small residual offsets (and phase shifts) between basins. The simple model predicts that this effect should decrease roughly linearly with frequency, a result that is confirmed by numerical modeling and empirical T/P estimates of the monthly (Mm) tidal constituent. This model also explains some aspects of the anomalous nonisostatic response of the ocean to atmospheric pressure forcing at periods of around 5 days.

Egbert, Gary D.↗
NASA NTRS2015

Nonequilibrium Supersonic Freestream Studied Using Coherent Anti-Stokes Raman Spectroscopy

Measurements were conducted at the University of Virginia Supersonic Combustion Facility of the flow in a constant-area duct downstream of a Mach 2 nozzle. The airflow was heated to approximately 1200 K in the facility heater upstream of the nozzle. Dual-pump coherent anti-Stokes Raman spectroscopy was used to measure the rotational and vibrational temperatures of N2 and O2 at two planes in the duct. The expectation was that the vibrational temperature would be in equilibrium, because most scramjet facilities are vitiated air facilities and are in vibrational equilibrium. However, with a flow of clean air, the vibrational temperature of N2 along a streamline remains approximately constant between the measurement plane and the facility heater, the vibrational temperature of O2 in the duct is about 1000 K, and the rotational temperature is consistent with the isentropic flow. The measurements of N2 vibrational temperature enabled cross-stream nonuniformities in the temperature exiting the facility heater to be documented. The measurements are in agreement with computational fluid dynamics models employing separate lumped vibrational and translational/rotational temperatures. Measurements and computations are also reported for a few percent steam addition to the air. The effect of the steam is to bring the flow to thermal equilibrium, also in agreement with the computational fluid dynamics.

Cutler, Andrew D.↗
NASA NTRS1973

Application of modern control theory to the design of optimum aircraft controllers

The synthesis procedure presented is based on the solution of the output regulator problem of linear optimal control theory for time-invariant systems. By this technique, solution of the matrix Riccati equation leads to a constant linear feedback control law for an output regulator which will maintain a plant in a particular equilibrium condition in the presence of impulse disturbances. Two simple algorithms are presented that can be used in an automatic synthesis procedure for the design of maneuverable output regulators requiring only selected state variables for feedback. The first algorithm is for the construction of optimal feedforward control laws that can be superimposed upon a Kalman output regulator and that will drive the output of a plant to a desired constant value on command. The second algorithm is for the construction of optimal Luenberger observers that can be used to obtain feedback control laws for the output regulator requiring measurement of only part of the state vector. This algorithm constructs observers which have minimum response time under the constraint that the magnitude of the gains in the observer filter be less than some arbitrary limit.

Power, L. J.↗
NASA NTRS1972

Equilibrium configurations of Vlasov plasmas carrying a current component along an external magnetic field.

A model of equilibrium configurations of Vlasov plasmas is considered which represents a combination of the models of Harris (1962) and Nicholson (1963). These plasma configurations carry a current component along an external magnetic field. The considered slab model contains a diamagnetic current and a field-aligned current for an arbitrary ratio of particle pressure to magnetic pressure of the applied constant field. For a fixed pressure ratio and field-aligned current, the model admits a family of equilibrium solutions in which the diamagnetic currents range from zero to a maximum value. The amount of diamagnetic current flowing in a machine depends on the width of the machine, the field-aligned current and other plasma parameters.

Kan, J. R.↗
NASA NTRS1981

Study of certain launching techniques using long orbiting tethers

A study of the basic equations governing orbital transfers using long orbiting tethers is presented. A very simple approximation to the general transfer equation is derived for the case of short tethers and low eccentricity orbits. Numerical examples are calculated for the case of injection into a circular orbit from a platform in eccentric orbit and injection into eccentric orbit from a platform in circular orbit. For the case of long tethers, a method is derived for reducing tether mass and increasing payload mass by tapering the tether to maintain constant stress per unit of tether cross section. Formulas are presented for calculating the equilibrium orbital parameters taking into account the mass of the platform, tether, and payload.

Colombo, G.↗
NASA NTRS1986

Submillimeter and millimeter observations of Uranus and Neptune

Neptune and Uranus were studied at submillimeter and millimeter wavelengths for 2 yr using the UK IR telescope and the NASA IR telescope at Mauna Kea and the 12 dish at Kitt Peak, AZ. The brightness temperatures of both planets are provided at various wavelengths over the 0.35-3.3 mm interval scanned in the study. The present data, in combination with data from the 17 microns-3 mm interval scanned in the study. The present data, in combination with data from the 17 microns-3 mm interval taken during previous campaigns, are employed to model the temperature structures in the atmospheres of the two planets over the vertical range 100 mbar-8 bar. The data support the presence of active convection processes at the 1 bar level and deeper in both atmospheres. A static equilibrium however, exists between para-H2 and ortho-H2, along with a constant CH4 fraction of 2 percent. Deficiencies in the current database are detailed.

Orton, G. S.↗
NASA NTRS2000

Some Thermodynamic Considerations on the Physical and Quantum Nature of Space and Time

It is suggested that the Planck h = m(sub k)c Lambda(sub k) and the Boltzmann k = m(sub k)c nu(sub k)Constants have stochastic foundation. It is further suggested that a body of fluid at equilibrium is composed of a spectrum of molecular clusters (energy levels) the size of which are governed by the Maxwell-Boltzmann distribution function. Brownian motions are attributed to equilibrium between suspensions and molecular clusters. Atomic (molecular) transition between different size atomic- (molecular-) clusters (energy levels) is shown to result in emission/absorption of energy in accordance with Bohr's theory of atomic spectra. Physical space is identified as a tachyonic fluid that is Dirac's stochastic ether or de Broglie's hidden thermostat. Compressibility of physical space, in accordance with Planck's compressible ether, is shown to result in the Lorentz-Fitzgerald contraction, thus providing a causal explanation of relativistic effect in accordance with the perceptions of Poincare and Lorentz. The invariant Schrodinger equation is derived from the invariant Bernoulli equation for incompressible potential flow. Following Heisenberg a temporal uncertainty relation is introduced as Delta(nu(sub Beta)) Delta(Rho(sub Beta)) > = k.

Sohrab, Siavash H.↗
NASA NTRS1975

Experimental results for the transpired turbulent boundary layer in an adverse pressure gradient

The fluid mechanics of transpired incompressible turbulent boundary layers under zero and adverse pressure gradient conditions is investigated using an open-ended wind tunnel with a porous floor in the test section and a secondary air system for supply and metering of the transpiration air. All velocity profiles and turbulence profiles are obtained by linearized constant-temperature hot-wire anemometry. The wall shear stress is determined by measuring the shear stress away from the wall and extrapolating to the wall by integrating the boundary layer equations for the shear-stress profile. Equilibrium boundary layers are obtained when the transpiration velocity is varied such that the blowing parameter and the Clauser pressure gradient parameter are held constant. The experimental results obtained are presented in tabular and graphical forms.

Andersen, P. S.↗
NASA NTRS1972

Nutational stability of a dual-spin satellite under the influence of applied reaction torques.

The nutational stability of a dual-spin satellite under the influence of an applied reaction torque is analyzed. It is assumed that since the solar paddle attachments to the hub of the spacecraft are not rigidly locked, the effect of the solar panels can be replaced by a constant reaction torque acting on the hub of the spacecraft. This could result in the satellite having an equilibrium motion about an axis displaced from the nominal axis of symmetry. The variational equations of motion are developed about such an equilibrium position using the SAS-A spacecraft as a model. Energy dissipation on the rotor as well as the main body is included. This nonautonomous set of differential equations are linearized and transformed to an autonomous set using the Lyapunov Reducibility Theorem. The stability of the kinematically similar system is examined numerically using representative SAS-A parameters for the case when either pair of solar panels is assumed to be loosely attached. Stability is verified for small system nutation angles (0.1 degree) although the time constant associated with the least damped mode is approximately one order of magnitude larger than that for the SAS-A system without the influence of the reaction torques.

Bainum, P. M.↗
NASA NTRSDate not supplied

A Strategic Approach for Dense, Integrated, Vehicle Navigation

Drone usage has been on the rise in recent years with applications that include parcel delivery, wildlife protection, precision farming, law enforcement, and industrial inspection, just to name a few. Once regulations and safety policies are put in place to allow for the widespread use of unmanned drones, the number of aircraft in the National Airspace System (NAS) is expected to skyrocket to millions, potentially congesting the airspace which increases the likelihood of separation violations and possibly incidents. Currently, flight infrastructure can only support a few thousand aircraft flying over the United States National Airspace System (NAS) at any given time. A delay at one airport can send ripple effects throughout the system, causing more delays and missed connections. In air traffic control, separation is the concept of keeping an “ownship” aircraft outside a minimum distance from “intruder” aircraft to reduce the risk of the aircraft colliding, as well as preventing accidents due to secondary factors, such as wake turbulence. Maintaining proper separation is often a safety critical property for fixed-wing drones in the airspace. This paper addresses drone separation in time and distance for high volume corridors (en route) and lanes (on ground), merging as well as crossing intersections of multiple corridors/lanes. In this paper, the term drone is applied to both Unmanned Aerial Vehicle (UAV) and small Unmanned Aircraft System (UAS) vehicles operating autonomously. There exists a gamut of approaches to the merging and crossing problems. At one end of the extreme are the conservative yet low cost and verifiable solutions of today that deal with two drones at a time. At the other end are complex Machine Learning-based solutions with high computing requirements for fully autonomous drones of the future that are expected to handle all contentions. This paper presents a feasible and verifiable strategic approach to these problems that is based on distributed cooperation between the drones and the infrastructure. Three phases of the strategic approach (Prepare, Adjust, Commit) are presented. Simulation results are presented that show the proposed approach is stable and resilient to induced perturbations and guarantees a set of fixed-wing drones to merge and cross intersections by adjusting their speed based on their distance to the aircraft in front of them while remaining in the equilibrium state. The equilibrium state is defined as the state when a set of n aircraft move at a relatively constant speed and uniform spacing from each other in a congested system. A congested system is defined as the state when at least one aircraft cannot move at its maximum allowed speed. Unlike existing centralized and pre-planned approaches, the proposed solution is fully distributed and enables autonomous aircraft to decide to adjust their speed and distance with respect to the preceding aircraft, dynamically. Simulation results are presented that assess the feasibility of the approach.

Distributed↗
NASA NTRS1973

Application of modern control theory to the design of optimum aircraft controllers

A procedure is described for synthesis of optimal aircraft control systems by application of the concepts of optimal control theory to time-invariant linear systems with quadratic performance criteria. Essential in this synthesis procedure is the solution of the Riccati matrix equation which results in a constant linear feedback control law for an output regulator which maintains a plant in an equilibrium in the presence of impulse disturbances. An algorithm is derived for designing maneuverable output regulators with selected state variables for feedback.

Power, L. J.↗
NASA NTRS1980

Interaction of a cumulus cloud ensemble with the large-scale environment. II

The closure assumption of the Arakawa-Schubert (1974) cumulus parameterization is examined with attention to the kinetic energy budget of a cumulus subensemble. It is concluded that dissipation should depend primarily on cloud type, and that cloud-work functions calculated under different synoptic conditions should consequently be quasi-constant for each cloud type. Observational evidence for the cloud-work function quasi-equilibrium is presented. Cloud-work functions are calculated from a variety of data sets in the tropics and subtropics. The results demonstrate that the cloud-work functions fall into a well-defined narrow range for each subensemble, although the thermodynamical vertical structures for each data set are quite different.

Lord, S. J.↗
NASA NTRS1986

Meridional distribution of stratospheric trace constituents

Vertically stratified tracers such as methane and nitrous oxide tend to have constant mixing ratio surfaces that slope downward toward the poles in the meridional plane. The equilibrium tracer slope results from the competition between the slope steepening effects of advection by the diabatic circulation and the slope flattening effects of quasi-isentropic eddy transport and photochemical loss. The diabatic circulation iself is, however, driven primarily by eddy transports, which maintain the departure of stratospheric temperatures from radiative equilibrium. If the eddy transports are weak, the diabatic circulation is also weak and the slope is small. Using a simple beta-plane channel model and an eddy diffusion parameterization for the eddy potential vorticity and tracer transports, it is shown that the slope is a maximum for a value of eddy diffusion such that the dynamical time scale lies between the radiative and chemical time scales.

Holton, J. R.↗
NASA NTRS1989

The mean ozone profile and its temperature sensitivity in the upper stratosphere and lower mesosphere - An analysis of LIMS observations

Multiple simultaneous LIMS (Limb Infrared Monitor of the Stratosphere) observations of both O3 and temperature T in the upper stratosphere and lower mesosphere are analyzed to derive O3-T correlations for comparisons with model calculations. The zonally averaged O3 profile is compared with results from a simplified photochemical model that assumes O3 to be in a photochemical steady state. The model profile is systematically lower than the observed profile. Key parameters are identified in which changes can systematically increase the O3 profile in both the stratosphere and the mesosphere. The LIMS-derived values for the sensitivity of O3 to changes in T are compared with equilibrium model calculations which include the T-driven opacity feedback effect on photodissociation rate constants. The theoretical O3 response to temperature perturbations is investigated. It is shown how the O3-T sensitivity coefficient is affected by zonal and vertical advection terms as well as by the photochemical coupling between O3 and T.

Froidevaux, Lucien↗
NASA NTRS1993

Motion of the heliospheric termination shock - A gas dynamic model

A simple quantitative model is presented for the heliospheric termination shock's anticipated movement in response to upstream solar wind condition variations, under the assumption that the termination shock is initially a strong gasdynamic shock that is at rest relative to the sun, and that there is a discontinuous increase or decrease in the dynamical pressure upstream of the shock. The model suggests that the termination shock is constantly in motion, and that the mean position of the shock lies near the mean equilibrium position which corresponds to the balance between the mean solar wind dynamical pressure and the mean interstellar pressure.

Barnes, Aaron↗