Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “equilibrium modeling”

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 793 records · Page 44

Quasi-Linear Parameter Varying Representation of General Aircraft Dynamics Over Non-Trim Region

For applying linear parameter varying (LPV) control synthesis and analysis to a nonlinear system, it is required that a nonlinear system be represented in the form of an LPV model. In this paper, a new representation method is developed to construct an LPV model from a nonlinear mathematical model without the restriction that an operating point must be in the neighborhood of equilibrium points. An LPV model constructed by the new method preserves local stabilities of the original nonlinear system at "frozen" scheduling parameters and also represents the original nonlinear dynamics of a system over a non-trim region. An LPV model of the motion of FASER (Free-flying Aircraft for Subscale Experimental Research) is constructed by the new method.

Shin, Jong-Yeob↗

A survey of upwind methods for flows with equilibrium and non-equilibrium chemistry and thermodynamics

Several versions of flux-vector split and flux-difference split algorithms were compared with regard to general applicability and complexity. Test computations were performed using curve-fit equilibrium air chemistry for an M = 5 high-temperature inviscid flow over a wedge, and an M = 24.5 inviscid flow over a blunt cylinder for test computations; for these cases, little difference in accuracy was found among the versions of the same flux-split algorithm. For flows with nonequilibrium chemistry, the effects of the thermodynamic model on the development of flux-vector split and flux-difference split algorithms were investigated using an equilibrium model, a general nonequilibrium model, and a simplified model based on vibrational relaxation. Several numerical examples are presented, including nonequilibrium air chemistry in a high-temperature shock tube and nonequilibrium hydrogen-air chemistry in a supersonic diffuser.

Grossman, B.↗

X-ray emission from young supernova remnants - Nonionization equilibrium abundances and emissivities

X-ray line emission from hot, low-density plasma in young supernova remnants is strongly enhanced by departures from ionization equilibrium. The X-ray emission from a Sedov blast wave has been calculated using a nonequilibrium evolutionary treatment of the ionization structure, and the resulting spectrum has been fitted to HEAO 2 SSS data for Tycho's remnant. These models yield dramatically different elemental abundances for heavy elements (Ne, Mg, Si, S, Ar, Ca, and Fe), compared with two-temperature component equilibrium models of Becker et al. (1979, 1980). Nonequilibrium broad-band X-ray emissivities result in lower mass determinations for the supernova ejecta. Areas of further improvement of remnant X-ray modeling are suggested.

Shull, J. M.↗

A wet-geology and cold-climate Mars model: Punctuation of a slow dynamic approach to equilibrium

It was suggested that Mars may have possessed a relatively warm humid climate and a vigorous hydrological cycle involving meteoric precipitation, oceans, and continental ice sheets. Baker hypothesized that these geologically active conditions may have been repeated several times; each of these dynamic epochs was followed by a collapse of the climate and hydrologic cycle of Mars into essentially current conditions, completing what is termed a 'Baker cycle'. The purpose is to present an endmember possibility that Martian glacial landscapes, including some that were previously considered to have formed under warm climatic conditions, might be explained by processes compatible with an extremely cold surface. Two aspects of hypothesized Martian glacial terrains were cited as favoring a warm climate during Baker cycles: (1) the formation of some landscapes, including possible eskers, tunnel channels, drumlins, and outwash plains, appears to have required liquid water, and (2) a liquid-surfaced ocean was probably necessary to feed the glaciers. The requirement for liquid water, if these features were correctly interpreted, is difficult to avoid; it is entirely possible that a comparatively warm climate was involved, but it is not clear that formation of landforms by wet-based glaciers actually requires a warm climate. Even less certain is the supposed requirement for liquid oceans. Formation of glaciers only requires a source of water or ice to supply an amount of precipitation that exceeds losses due to melting and sublimation. At Martian temperatures precipitation is very low, but so are melting and sublimation, so a large body of ice that is unstable with respect to sublimation may take the role of Earth's oceans in feeding the glaciers. Recent models suggest that even current Martian polar caps, long thought to be static bodies of ice and dust, might actually be slow-moving, cryogenic continental glaciers. Is it possible that subglacial processes beneath cryogenic (but wetbased) ice sheets formed the hypothesized Martian glacial landscapes?

Kargel, J. S.↗

Recent Advancements in Modeling and Simulation of Entry Systems at NASA

This paper describes recent development of modeling and simulation technologies for entry systems in support of NASA's exploration missions. Mission-tailored research and development in modeling of entry systems occurs across the Agency (e.g., within the Orion and Mars 2020 Programs), however the aim of this paper is to discuss the broad, cross-mission research conducted by NASA's Entry Systems Modeling (ESM) Project, which serves as the Agency's only concerted effort toward advancing entry systems across a range of technical disciplines. Technology development in ESM is organized and prioritized from a system-level perspective, resulting in four broad technical areas of investment: (1) Predictive material modeling, (2) Shock layer kinetics and radiation, (3) Computational and experimental aerosciences, and (4) Guidance, navigation, and control. Investments in thermal protection material modeling are geared toward high-fidelity, predictive models capable of handling complex structures, with an eye toward optimizing design performance and quantifying thermal protection system reliability. New computational tools have been developed to characterize material properties and behavior at the microstructural level, and experimental techniques (molecular beam scattering, micro-computed tomography, among others) have been developed to measure material kinetics, morphology, and other parameters needed to inform and validate detailed simulations. Advancements have also been made in macrostructural simulation capability to enable 3-D system-scale calculations of material response with complex topological features, including differential recession of tile gaps. Research and development in the area of shock layer kinetics has focused on air and CO2-based atmospheres. Capacity and capability of the NASA Ames Electric Arc Shock Tube (EAST) have been expanded in recent years and analysis of resulting data has led to several improvements in kinetic models, while simultaneously reducing uncertainties associated with radiative heat transfer predictions. First-principles calculations of fundamental kinetic, thermodynamic, and transport data, along with state-specific models for non-equilibrium flow regimes, have also yielded new insights and have the potential to vastly improve model fidelity. Aerosciences is a very broad area of interest in entry systems, yet a number of important challenges are being addressed: Coupled fluid-structure simulations of parachute inflation and dynamics; Experimental and computational studies of vehicle dynamics; Multi-phase flow with dust particles to simulate entry environments at Mars during dust storms; Studies of roughness-induced heating augmentation relevant to tiled and woven thermal protection systems; and Advanced numerical methods to optimize computational analyses for desired accuracy versus cost. Guidance and control in the context of entry systems has focused on development of methods for multi-axis control (i.e. pitch and yaw, rather than bank angle alone) of spacecraft during entry and descent. With precision landing requirements driven by Mars human exploration goals, recent efforts have yielded 6-DOF models of multi-axis control with propulsive descent of both inflatable and rigid ellipsled-like architectures.

Barnhardt, Michael D.↗

Non-LTE line-blanketed model atmospheres of hot stars. 1: Hybrid complete linearization/accelerated lambda iteration method

A new munerical method for computing non-Local Thermodynamic Equilibrium (non-LTE) model stellar atmospheres is presented. The method, called the hybird complete linearization/accelerated lambda iretation (CL/ALI) method, combines advantages of both its constituents. Its rate of convergence is virtually as high as for the standard CL method, while the computer time per iteration is almost as low as for the standard ALI method. The method is formulated as the standard complete lineariation, the only difference being that the radiation intensity at selected frequency points is not explicity linearized; instead, it is treated by means of the ALI approach. The scheme offers a wide spectrum of options, ranging from the full CL to the full ALI method. We deonstrate that the method works optimally if the majority of frequency points are treated in the ALI mode, while the radiation intensity at a few (typically two to 30) frequency points is explicity linearized. We show how this method can be applied to calculate metal line-blanketed non-LTE model atmospheres, by using the idea of 'superlevels' and 'superlines' introduced originally by Anderson (1989). We calculate several illustrative models taking into accont several tens of thosands of lines of Fe III to Fe IV and show that the hybrid CL/ALI method provides a robust method for calculating non-LTE line-blanketed model atmospheres for a wide range of stellar parameters. The results for individual stellar types will be presented in subsequent papers in this series.

Hubeny, I.↗

Tropical Convection and Climate Processes in a Cumulus Ensemble Model

Local convective-radiative equilibrium states of the tropical atmosphere are determined by the following external forcing: 1) Insolation, 2) Surface heat and moisture exchanges (primarily radiation and evaporation), 3) Heating and moistening induced by large-scale circulation. Understanding the equilibrium states of the tropical atmosphere in different external forcing conditions is of vital importance for studying cumulus parameterization, climate feedbacks, and climate changes. We extend our previous study using the Goddard Cumulus Ensemble (GCE) Model which resolves convective-radiative processes more explicitly than global climate models do. Several experiments are carried out under fixed insolation and sea surface temperature. The prescribed SST consists of a uniform warm pool (29C) surrounded by uniform cold SST (26C). The model produces "Walker"-type circulation with the ascending branch of the model atmosphere more humid than the descending part, but the vertically integrated temperature does not show a horizontal gradient. The results are compared with satellite measured moisture by SSM/I (Special Sensor Microwave/Imager) and temperature by MSU in the ascending and descending tropical atmosphere. The vertically integrated temperature and humidity in the two model regimes are comparable to the observed values in the tropics.

Sui, Chung-Hsiung↗

Effects of Injected Current Streams on MHD Equilibrium Reconstruction of Local Helicity Injection Plasmas in a Spherical Tokamak

Open field line currents are intrinsic to DC helicity injection plasma startup and pose a challenge for inferring the plasma equilibrium with standard reconstruction analysis. Local helicity injection (LHI) is a type of DC helicity injection which uses small, modular current sources to drive force-free current along helical field lines to produce tokamak plasmas. MHD modeling and magnetic measurements during LHI indicate the injected current streams remain coherent as helical structures on the outboard edge of a core toroidal plasma that is tokamak-like in a toroidally averaged sense. To extract core plasma equilibrium properties, external magnetic diagnostics corrected for contributions from the injected current streams are fitted by a standard Grad-Shafranov equilibrium code. An iterative approach for estimating and subtracting the stream contributions from the diagnostic signals is described and applied to a model equilibrium database to reduce systematic errors introduced by the streams. Convergence is usually attained with 2 to 4 iterations, with derived equilibrium parameters matching the prescribed axisymmetric core values to within estimated experimental uncertainties. Accurate recovery of core parameters occurs when the ratio of the net toroidal windup current from the streams to the core plasma current is less than 0.2, which is typically satisfied in most experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Inverse design of equilibrium cluster fluids applied to a physically informed model

Inverse design strategies have proven highly useful for the discovery of interaction potentials that prompt self-assembly of a variety of interesting structures. However, often the optimized particle interactions do not have a direct relationship to experimental systems. In this work, we show that Relative Entropy minimization is able to discover physically meaningful parameter sets for a model interaction built from depletion attraction and electrostatic repulsion that yield self-assembly of size-specific clusters. Here, we then explore the sensitivity of the optimized interaction potentials with respect to deviations in the underlying physical quantities, showing that clustering behavior is largely preserved even as the optimized parameters are perturbed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Statistical equilibrium in cometary C2. IV - A 10 level model including singlet-triplet transitions

Resonance fluorescence theory was used to calculate the population distribution in the energy states of the C2 molecule in comets. Ten electronic states, each with 14 vibrational states, were used in the calculations. These new calculations differ from earlier work in terms of additional electronic levels and the role of singlet-triplet transitions between the b and X levels. Since transition moments are not known, calculations are made of observable flux ratios for an array of possible values. Comparison with existing observations indicates that the a-X transition is very important, and there is marginal indication that the b-X transition is present. Swan band sequence flux ratios at large heliocentric distance are needed, as are accurate Mulliken/Swan and Phillips/Ballik-Ramsay (1963) observations.

Krishna Swamy, K. S.↗

Drying model of a high salt content cementitious waste form: Effect of capillary forces and salt solution

Highlights: • Drying model for a high salt content cementitious waste form is developed. • Water vapor diffusion and capillary liquid flow are distinguished. • Capillary and salt solution effects are considered in vapor-liquid equilibrium. A water transport model coupling capillary liquid flow with vapor diffusion is developed to describe the drying process for a cementitious waste form with high salinity porewater. Vapor-liquid equilibrium is formulated as the driving force for vapor diffusion and the model accounts for pore capillary and high salinity effects on water thermodynamic activity. Pore filling and porewater surface tension as a function of pore size distribution and water saturation have been quantified for the material. Geochemical speciation modeling is used to simulate porewater activity as a function of composition over the range of saturation. The theoretical relationship between relative humidity and water saturation generally agrees with experimental measurement, and the developed model is capable of predicting drying rates under various external relative humidity conditions. The model was developed to be incorporated into reactive transport models considering the effects of drying such as salt redistribution and efflorescence.

36 MATERIALS SCIENCE↗

LGM Paleoclimate Constraints Inform Cloud Parameterizations and Equilibrium Climate Sensitivity in CESM2

The Community Earth System Model version 2 (CESM2) simulates a high equilibrium climate sensitivity (ECS > 5°C) and a Last Glacial Maximum (LGM) that is substantially colder than proxy temperatures. In this study, we examine the role of cloud parameterizations in simulating the LGM cooling in CESM2. Through substituting different versions of cloud schemes in the atmosphere model, we attribute the excessive LGM cooling to the new CESM2 schemes of cloud microphysics and ice nucleation. Further exploration suggests that removing an inappropriate limiter on cloud ice number (NoNimax) and decreasing the time-step size (substepping) in cloud microphysics largely eliminate the excessive LGM cooling. NoNimax produces a more physically consistent treatment of mixed-phase clouds, which leads to an increase in cloud ice content and a weaker shortwave cloud feedback over mid-to-high latitudes and the Southern Hemisphere subtropics. Microphysical substepping further weakens the shortwave cloud feedback. Based on NoNimax and microphysical substepping, we have developed a paleoclimate-calibrated CESM2 (PaleoCalibr), which simulates well the observed twentieth century warming and spatial characteristics of key cloud and climate variables. PaleoCalibr has a lower ECS (∼4°C) and a 20% weaker aerosol-cloud interaction than CESM2. PaleoCalibr represents a physically more consistent treatment of cloud microphysics than CESM2 and is a valuable tool in climate change studies, especially when a large climate forcing is involved. Our study highlights the unique value of paleoclimate constraints in informing the cloud parameterizations and ultimately the future climate projection.

equilibrium climate sensitivity↗

Transient cryogenic liquid discharge in normal and micro-gravity

Transient discharge of liquid cryogens from partially filled tanks under microgravity is considered. The discharge is made to a low-pressure site through a straight tube. The ultimate objective of the study is to predict the transient mass flow rate. In addition, the pressure history in the tank is of interest, since the pressure falls below the triple point toward the end of the dump, causing solidification in the tank and/or the discharge line. As first steps, experiments have been performed in normal gravity and attempts are made to predict the transient rate and pressure history under this condition. Experiments are performed with liquid nitrogen dumped to sites with pressures below the triple point. Homogeneous equilibrium and nonequilibrium two-phase critical flow models are compared to the experimental results. Equilibrium and nonequilibrium tank pressure history models are also compared to experimental data.

Purwin, T. P.↗

Dispersive Charge Transfer State Electroluminescence in Organic Solar Cells

Abstract The notion of quasi‐equilibrium is central to most solar cells; however, it has been questioned in organic photovoltaics (OPVs) owing to strong energetic disorder that frustrates efficient relaxation of electrons and holes within their respective density of states (DOS). Here, modulation electroluminescence (EL) spectroscopy is applied to OPVs and it is found that the frequency response of charge transfer (CT) state EL on the high energy side of the spectrum differs from that of the low energy side. This observation confirms that static disorder contributes substantially to the linewidth of the steady‐state EL spectrum and is unambiguous proof that the distribution of CT states formed by electrical injection in the dark is not in quasi‐equilibrium. These results emphasize the need for caution when analyzing OPV cells on the basis of reciprocity models that assume quasi‐equilibrium holds, and highlight a new method to study this unusual aspect of OPV operation.

14 SOLAR ENERGY↗

Thermal and dynamic partition of dumbbell interstitials in complex concentrated alloys

Complex concentrated alloys (CCAs) are promising candidates for applications in extreme conditions, such as irradiation where interstitial mediated diffusion is important. In CCAs with N principal elements, $\frac{N(N+1)}{2}$ types of dumbbell interstitials exist. Currently, there is no way to predict the thermal partition (fractional concentration at equilibrium) and the dynamic partition (fractional time an interstitial spends during diffusion) of each type of dumbbell interstitial. To mitigate this issue, this work proposes a theoretical model for computing the equilibrium concentrations and thermal partition of dumbbell interstitials in CCAs and validates the model using grand canonical Monte Carlo simulations. Lattice kinetic Monte Carlo simulations show that the thermal partition is equivalent to the dynamic partition, and both are governed by composition and formation energies of dumbbells. As a result, the model proposed provides a foundation for understanding radiation enhanced diffusion and induced segregation in CCAs under irradiation.

36 MATERIALS SCIENCE↗

On Kinetics Modeling of Vibrational Energy Transfer

Two models of vibrational energy exchange are compared at equilibrium to the elementary vibrational exchange reaction for a binary mixture. The first model, non-linear in the species vibrational energies, was derived by Schwartz, Slawsky, and Herzfeld (SSH) by considering the detailed kinetics of vibrational energy levels. This model recovers the result demanded at equilibrium by the elementary reaction. The second model is more recent, and is gaining use in certain areas of computational fluid dynamics. This model, linear in the species vibrational energies, is shown not to recover the required equilibrium result. Further, this more recent model is inconsistent with its suggested rate constants in that those rate constants were inferred from measurements by using the SSH model to reduce the data. The non-linear versus linear nature of these two models can lead to significant differences in vibrational energy coupling. Use of the contemporary model may lead to significant misconceptions, especially when integrated in computer codes considering multiple energy coupling mechanisms.

Gilmore, John O.↗

Evaluation of a distributed catchment scale water balance model

The validity of some of the simplifying assumptions in a conceptual water balance model is investigated by comparing simulation results from the conceptual model with simulation results from a three-dimensional physically based numerical model and with field observations. We examine, in particular, assumptions and simplifications related to water table dynamics, vertical soil moisture and pressure head distributions, and subsurface flow contributions to stream discharge. The conceptual model relies on a topographic index to predict saturation excess runoff and on Philip's infiltration equation to predict infiltration excess runoff. The numerical model solves the three-dimensional Richards equation describing flow in variably saturated porous media, and handles seepage face boundaries, infiltration excess and saturation excess runoff production, and soil driven and atmosphere driven surface fluxes. The study catchments (a 7.2 sq km catchment and a 0.64 sq km subcatchment) are located in the North Appalachian ridge and valley region of eastern Pennsylvania. Hydrologic data collected during the MACHYDRO 90 field experiment are used to calibrate the models and to evaluate simulation results. It is found that water table dynamics as predicted by the conceptual model are close to the observations in a shallow water well and therefore, that a linear relationship between a topographic index and the local water table depth is found to be a reasonable assumption for catchment scale modeling. However, the hydraulic equilibrium assumption is not valid for the upper 100 cm layer of the unsaturated zone and a conceptual model that incorporates a root zone is suggested. Furthermore, theoretical subsurface flow characteristics from the conceptual model are found to be different from field observations, numerical simulation results, and theoretical baseflow recession characteristics based on Boussinesq's groundwater equation.

Troch, Peter A.↗

Atmospheric reentry of the in-core thermionic SP-100 reactor system

Presumed end-of-life atmospheric reentry of the GA SP-100 system was studied to assess dispersal feasibility and associated hazards. Reentry was studied by sequential use of an orbital trajectory and a heat analysis computer program. Two heating models were used. The first model assumed a thermal equilibrium condition between the stagnation point aerodynamic heating and the radiative cooling of the skin material surface. The second model allowed for infinite conductivity of the skin material. Four reentering configurations were studied representing stages of increased SP-100 breakup: (1) radiator, shield and reactor, (2) shield and reactor, (3) reactor with control drums, and (4) reactor without control drums. Each reentering configuration was started from a circular orbit at 116 km having an inertial velocity near Mach 25. The assumed failing criterion was the attainment of melting temperature of a critical system component. The reentry analysis revealed breakup of the vessel in the neighborhood of 61 km altitude and scattering of the fuel elements. Subsequent breakup of the fuel elements was not predicted. Oxidation of the niobium skin material was calculated to cause an increase in surface temperature of less than ten percent. The concept of thermite analogs for enhancing reactor reentry dispersal was assessed and found to be feasible in principle. A conservative worst-case hazards analysis was performed for radioactive and nonradioactive toxic SP-100 materials assumed to be dispersed during end-of-life reentry. The hazards associated with this phase of the SP-100 mission were calculated to be insignificant.

Stamatelatos, M. G.↗