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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
NASA NTRS2012

Changes in Extratropical Storm Track Cloudiness 1983-2008: Observational Support for a Poleward Shift

Climate model simulations suggest that the extratropical storm tracks will shift poleward as a consequence of global warming. In this study the northern and southern hemisphere storm tracks over the Pacific and Atlantic ocean basins are studied using observational data, primarily from the International Satellite Cloud Climatology Project, ISCCP. Potential shifts in the storm tracks are examined using the observed cloud structures as proxies for cyclone activity. Different data analysis methods are employed, with the objective to address difficulties and uncertainties in using ISCCP data for regional trend analysis. In particular, three data filtering techniques are explored; excluding specific problematic regions from the analysis, regressing out a spurious viewing geometry effect, and excluding specific cloud types from the analysis. These adjustments all, to varying degree, moderate the cloud trends in the original data but leave the qualitative aspects of those trends largely unaffected. Therefore, our analysis suggests that ISCCP data can be used to interpret regional trends in cloudiness, provided that data and instrumental artefacts are recognized and accounted for. The variation in magnitude between trends emerging from application of different data correction methods, allows us to estimate possible ranges for the observational changes. It is found that the storm tracks, here represented by the extent of the midlatitude-centered band of maximum cloud cover over the studied ocean basins, experience a poleward shift as well as a narrowing over the 25 year period covered by ISCCP. The observed magnitudes of these effects are larger than in current generation climate models (CMIP3). The magnitude of the shift is particularly large in the northern hemisphere Atlantic. This is also the one of the four regions in which imperfect data primarily prevents us from drawing firm conclusions. The shifted path and reduced extent of the storm track cloudiness is accompanied by a regional reduction in total cloud cover. This decrease in cloudiness can primarily be ascribed to low level clouds, whereas the upper level cloud fraction actually increases, according to ISCCP. Independent satellite observations of radiative fluxes at the top of the atmosphere are consistent with the changes in total cloud cover. The shift in cloudiness is also supported by a shift in central position of the mid-troposphere meridional temperature gradient. We do not find support for aerosols playing a significant role in the satellite observed changes in cloudiness. The observed changes in storm track cloudiness can be related to local cloud-induced changes in radiative forcing, using ERBE and CERES radiative fluxes. The shortwave and the longwave components are found to act together, leading to a positive (warming) net radiative effect in response to the cloud changes in the storm track regions, indicative of positive cloud feedback. Among the CMIP3 models that simulate poleward shifts in all four storm track areas, all but one show decreasing cloud amount on a global mean scale in response to increased CO2 forcing, further consistent with positive cloud feedback. Models with low equilibrium climate sensitivity to a lesser extent than higher-sensitivity models simulate a poleward shift of the storm tracks.

Bender, Frida A-M.↗
DOE OSTI2024

The Onset of Nucleate Boiling Prediction in a Mini Channel with the Eulerian Multiphase Flow

The prediction of boiling incipience is a critical issue for a reactor. Vapor lock causes operating instability and undesirable temperature rise, especially for small channels. This study uses the Eulerian multiphase flow boiling model to predict the onset of nucleate boiling in a millimeter scale rectangular channel. Two wall boiling models, namely the RPI (Rensselaer Polytechnic Institute) and the non-equilibrium sub cooled model are applied. Computational fluid dynamics is employed for this study to investigate the impact of operating conditions like liquid flow rate, operating pressure, heat flux, velocity profile, turbulence model, and inlet liquid temperature. Here, the effect of fundamental boiling parameters on the surface wall temperature and vapor volume fraction is also studied. The parameters under consideration are the bubble departure diameter, bubble departure frequency, nucleate site density, quenching time period, and interface heat transfer coefficients. A series of CFD calculations is conducted by varying the considered variables systematically in a wide range of flow conditions covering laminar, transition, and turbulent flows. The sub-component heat fluxes like the liquid convective, the vapor convective, the evaporative, and the quenching terms are monitored to infer the boiling dynamics in the transition region from the single-phase flow to the nucleate boiling zone. It is found that the CFD approach to detecting the boiling incipience point agrees reasonably with available experimental data. However, its limitations, like inaccuracy for lower flow rates and early transition to boiling flow, are also noticed due to the nature of the considered RPI model.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1985

A simple model for the interaction between vertical eddy heat fluxes and static stability

A numerical model for studying the interaction of vertical eddy heat fluxes and vertical temperature structure in midlatitude regions is described. The temperature profile for the model was derived from calculations of the equilibrium among heating rates in simplified representations of large-scale vertical eddy heat flux, moist convection and radiation. An eddy flux profile is calculated based on the quasi-geostrophic, liner baroclinic instability of a single wave. Model equilibrium states for summer and winter conditions are compared with observations, and the results are discussed in detail.

Gutowski, W. J., Jr.↗
NASA NTRS1994

Solar flare model atmospheres

Solar flare model atmospheres computed under the assumption of energetic equilibrium in the chromosphere are presented. The models use a static, one-dimensional plane-parallel geometry and are designed within a physically self-consistent coronal loop. Assumed flare heating mechanisms include collisions from a flux of nonthermal electrons and X-ray heating of the chromosphere by the corona. The heating by energetic electrons accounts explicitly for variations of the ionized fraction with depth in the atmosphere. X-ray heating of the chromosphere by the corona incorporates a flare loop geometry by approximating distant portions of the loop with a series of point sources, while treating the loop leg closest to the chromospheric footpoint in the plane-parallel approximation. Coronal flare heating leads to increased heat conduction, chromospheric evaporation and subsequent changes in coronal pressure; these effects are included self-consistently in the models. Cooling in the chromosphere is computed in detail for the important optically thick H I, Ca II and Mg II transitions using the non-local thermodynamic equilibrium (non-LTE) prescription in the program MULTI. Hydrogen ionization rates from X-ray photoionization and collisional ionization by nonthermal electrons are included explicitly in the rate equations. The models are computed in the 'impulsive' and 'equilibrium' limits, and in a set of intermediate 'evolving' states. The impulsive atmospheres have the density distribution frozen in the pre-flare configuration, while the equilibrium models assume the entire atmosphere is in hydrostatic and energetic equilibrium. The evolving atmospheres represent intermediate stages where hydrostatic equilibrium has been established in the chromosphere and corona, but the corona is not yet in energetic equilibrium with the flare heating source. Thus, for example, chromospheric evaporation is still in the process of occurring. We have computed the chromospheric radiation that results from a range of coronal heating rates, with particular emphasis on the widely observed diagnostic H(alpha). Our conclusion is that the H(alpha) fluxes and profiles actually observed in flares can only be produced under conditions of a low-pressure corona with strong beam heating. Therefore we suggest that H(alpha) in flares is produced primarily at the footprints of newly heated loops where significant evaporation has not yet occurred. As a single loop evolves in time, no matter how strong the heating rate may become, the H(alpha) flux will diminish as the corona becomes denser and hence more effective at stopping the beam. This prediction leads to several observable consequences regarding the spatial and temporal signatures of the X-ray and H(alpha) radiation during flares.

Hawley, Suzanne L.↗
DOE OSTI2020

A Model for Turbulence Spectra in the Equilibrium Range of the Stable Atmospheric Boundary Layer

Stratification can cause turbulence spectra to deviate from Kolmogorov's isotropic -5/3 power law scaling in the universal equilibrium range at high Reynolds numbers. However, a consensus has not been reached with regard to the exact shape of the spectra. Here we propose a shape of the turbulent kinetic energy and temperature spectra in horizontal wavenumber for the equilibrium range that consists of three regimes at small Froude number: the buoyancy subrange, a transition region, and the isotropic inertial subrange through dimensional analysis and substantial revision of previous theoretical approximation. These spectral regimes are confirmed by various observations in the atmospheric boundary layer. The representation of the transition region in direct numerical simulations will require large-scale separation between the Dougherty-Ozmidov scale and the Kolmogorov scale for strongly stratified turbulence at high Reynolds numbers, which is still challenging computationally. In addition, we suggest that the failure of Monin-Obukhov similarity theory in the very stable atmospheric boundary layer is due to the fact that it does not consider the buoyancy scale that characterizes the transition region.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1977

The stochastic control of the F-8C aircraft using a multiple model adaptive control /MMAC/ method. I - Equilibrium flight

The purpose of this paper is to summarize some results obtained for the adaptive control of the F-8C aircraft using the so-called MMAC method. The discussion includes the selection of the performance criteria for both the lateral and the longitudinal dynamics, the design of the Kalman filters for different flight conditions, the 'identification' aspects of the design using hypothesis testing ideas, and the performance of the closed-loop adaptive system.

Athans, M.↗
NASA NTRS2000

The Density-Potential (N(sub e)-V(sub s/c)) Relation in the High-Latitude Prenoon Ionosphere

Sounding of the Cleft Ion Fountain Energization Region (SCIFER) Thermal Electron Capped Hemisphere Spectrometer (TECHS) data are used to study the plasma density, spacecraft potential relation (N(sub e)-V(sub s/c)) in the prenoon topside auroral and cleft ionosphere during a period of low solar activity. The SCIFER TECHS data show a power-law electron temperature dependence on ambient density across the transition from positive to negative spacecraft potential values in the high-latitude ionosphere. The illuminated ionospheric N(sub e)-V(sub s/c) relation is numerically modeled by imposing equilibrium of current flowing to/from the conductive surfaces of the payload. This modeling demonstrates a strong dependence of spacecraft potential on electron temperature across a wide range of densities. This electron temperature dependence is especially significant in the negative spacecraft potential regime. SCIFER TECHS observations of the N(sub e)-V(sub s/c) relation are used to extend the previous magnetospheric observation-model results to the higher density ionospheric regime.

Adrian, M. L.↗
NASA NTRS2014

Simultaneous Laser-induced Fluorescence of Nitric Oxide and Atomic Oxygen in the Hypersonic Materials Environment Test System Arcjet Facility

Simultaneous nitric oxide (NO) and atomic oxygen (O) laser induced fluorescence (LIF) experiments were performed in the Hypersonic Materials Environmental Test System (HYMETS) facility at the NASA Langley Research Center. The data serves as an experimental database for validation for chemical and thermal nonequilibrium models used in hypersonic flows. Measurements were taken over a wide range of stagnation enthalpies (6.7 - 18.5 MJ/kg) using an Earth atmosphere simulant with a composition of 75% N2, 20% O2, and 5% Ar (by volume). These are the first simultaneous measurements of NO and O LIF to be reported in literature for the HYMETS facility. The maximum O LIF mean signal intensity was observed at a stagnation enthalpy of approximately 12 MJ/kg while the maximum NO LIF mean signal intensity was observed at a stagnation enthalpy of 6.7 MJ/kg. Experimental results were compared to simple fluorescence model that assumes equilibrium conditions in the plenum and frozen chemistry in the isentropic nozzle expansion (Mach 5). The equilibrium calculations were performed using CANTERA v2.1.1 with 16 species. The fluorescence model captured the correlation in mean O and NO LIF signal intensities over the entire range of stagnation enthalpies tested. Very weak correlations between single-shot O and NO LIF intensities were observed in the experiments at all of the stagnation enthalpy conditions.

Johansen, Craig↗
NASA NTRS2005

Implications for Core Formation of the Earth from High Pressure-Temperature Au Partitioning Experiments

Siderophile elements in the Earth.s mantle are depleted relative to chondrites. This is most pronounced for the highly siderophile elements (HSEs), which are approximately 400x lower than chondrites. Also remarkable is the relative chondritic abundances of the HSEs. This signature has been interpreted as representing their sequestration into an iron-rich core during the separation of metal from silicate liquids early in the Earth's history, followed by a late addition of chondritic material. Alternative efforts to explain this trace element signature have centered on element partitioning experiments at varying pressures, temperatures, and compositions (P-T-X). However, first results from experiments conducted at 1 bar did not match the observed mantle abundances, which motivated the model described above, a "late veneer" of chondritic material deposited on the earth and mixed into the upper mantle. Alternatively, the mantle trace element signature could be the result of equilibrium partitioning between metal and silicate in the deep mantle, under P-T-X conditions which are not yet completely identified. An earlier model determined that equilibrium between metal and silicate liquids could occur at a depth of approximately 700 km, 27(plus or minus 6) GPa and approximately 2000 (plus or minus 200) C, based on an extrapolation of partitioning data for a variety of moderately siderophile elements obtained at lower pressures and temperatures. Based on Ni-Co partitioning, the magma ocean may have been as deep as 1450 km. At present, only a small range of possible P-T-X trace element partitioning conditions has been explored, necessitating large extrapolations from experimental to mantle conditions for tests of equilibrium models. Our primary objective was to reduce or remove the additional uncertainty introduced by extrapolation by testing the equilibrium core formation hypothesis at P-T-X conditions appropriate to the mantle.

Danielson, L. R.↗
NASA NTRS1994

Aerothermodynamic methods for a Mars environmental survey Mars entry

Computational fluid dynamics models for the thermodynamics and transport properties used in an equilibrium version of the Langley aerothermodynamics upwind relaxation algorithm (LAURA) for Mars atmospheric entries are described. In addition, the physical models used in a nonequilibrium version of LAURA for Mars-entry flows are described. Uncertainties exist in defining constants used in the transport properties for the equilibrium model and in many of the physical models for the nonequilibrium version. Solutions from the two codes using the best available constants are examined at the Mars-entry conditions characteristics of the Mars environmental survey mission. While the flowfields are near thermal equilibrium, chemical nonequilibrium effects are present in the entry cases examined. Convective heating at the stagnation point for these flows (assuming fully catalytic wall boundary conditions) is approximately 100 W/cm(exp 2). Radiative heating is negligible.

AEROTHERMODYNAMIC METHODS↗
NASA NTRS2004

The Influence of Tungsten on the Chemical Composition of a Temporally Evolving Nanostructure of a Model Ni-Al-Cr Superalloy

The influence of W on the temporal evolution of gamma' precipitation toward equilibrium in a model Ni-Al-Cr alloy is investigated by three-dimensional atom-probe (3DAP) microscopy and transmission electron microscopy (TEM). We report on the alloys Ni-10 Al-8.5 Cr (at.%) and Ni-10 Al-8.5 Cr-2 W (at.%), which were aged isothermally in the gamma + gamma' two-phase field at 1073 K, for times ranging from 0.25 to 264 h. Spheroidal-shaped gamma' precipitates, 5-15 nm diameter, form during quenching from above the solvus temperature in both alloys at a high number density (approx. 10(exp 23/cu m). As gamma' precipitates grow with aging at 1073 K, a transition from spheriodal- to cuboidal-shaped precipitates is observed in both alloys. The elemental partitioning and spatially resolved concentration profiles across the gamma' precipitates are obtained as a function of aging time from three-dimensional atom-by-atom reconstructions. Proximity histogram concentration profiles of the quaternary alloy demonstrate that W concentration gradients exist in gamma' precipitates in the as-quenched and 0.25-h aging states, which disappear after 1 h of aging. The diffusion coefficient of W in gamma' is estimated to be 6.2 x 10(exp -20) sq m/s at 1073 K. The W addition decreases the coarsening rate constant, and leads to stronger partitioning of Al to gamma' and Cr to gamma.

Sudbrack, Chantal K.↗
NASA NTRS1993

Aerothermodynamics of a MESUR Martian entry

Computational fluid dynamics models for the thermodynamics and transport properties used in an equilibrium version of the Langley Aerothermodynamics Upwind Relaxation Algorithm (LAURA) for Martian Atmospheric entries are described. In addition, the physical models used in a nonequilibrium version of LAURA for Martian-entry flows are described. Uncertainties exist in defining constants used in the transport properties for the equilibrium model and in many of the physical models for the nonequilibrium version. Solutions from the two codes using the best available constants are examined at the Martian-entry conditions characteristic of the Mars Environmental Survey (MESUR) mission. While the flowfields are near thermal equilibrium, chemical nonequilibrium effects are present in the entry cases examined. Convective heating at the stagnation point for these flows (assuming fully catalytic wall boundary conditions) is approximately 100 W/sq cm. Radiative heating is negligible.

Mitcheltree, R. A.↗
NASA NTRS1995

Modeling and simulation of a turbulent far wake

Work continued on two projects which had been started during previous years. Both projects involve calculations of the subsonic, turbulent far wake of a two-dimensional object at a Reynolds number of 1000 (based on wake momentum thickness). This flow was used as a test case for direct comparison of various turbulence models and a direct numerical simulations (DNS) of this flow were undertaken. In the turbulence model comparison studies, for any particular model tested, a unique self-similar solution was obtained far enough downstream, regardless of inlet conditions. Furthermore, different turbulence models led to different far-wake equilibrium solutions. No turbulence model could correctly predict all features of the turbulent far wake. For example, the spreading rate and turbulent shear stresses were underpredicted by all the standard models (both two-equation and full Reynolds stress models). In cases where a more correct spreading rate was achieved, it was at the expense of the turbulent kinetic energy, which was overpredicted. In general, the Algebraic Dissipation Rate Model of Gatski and Speziale, 1992, when added to any of the standard models, improved the results dramatically. Also, full Reynolds stress closure models did a much better job at predicting the shapes of both the mean and turbulence profiles, but the spreading rate was not significantly improved over that predicted by the simpler two-equation models. There are two main conclusions from these studies: First, in a comparison such as this, it is not enough to compare just one parameter, like the spreading rate. A good prediction for one parameter does not necessarily imply good predictions for all parameters in a flow. Second, since no turbulence model could correctly predict the turbulent far wake, much of the important physics of turbulent free shear flows is apparently lost by the assumptions inherent in today's methods of turbulence modeling; turbulence models must be improved. Direct simulations of this flow were begun last year in order to provide a data base through which some of the deficiencies of the existing turbulence models could be identified. Quantities such as the pressure-strain correlation, turbulent diffusion, and the dissipation rate tensor can be easily calculated from the DNS results, whereas these quantities are nearly impossible to measure experimentally. Improvements to existing turbulence models (and development of new models) require knowledge about flow quantities such as these. During this summer, diagnostics codes were written which will calculate the parameters mentioned above, along with other single-point and multi-point statistics. The DNS calculations are still in progress at the time of this writing. When these calculations are complete, the diagnostics codes will be applied so that the results can aid turbulence modelers. In addition, the results will show whether or not there exists a universal equilibrium turbulent far wake, independent of initial conditions.

Cimbala, John M.↗
NASA NTRS1982

Development and application of the GIM code for the Cyber 203 computer

The GIM computer code for fluid dynamics research was developed. Enhancement of the computer code, implicit algorithm development, turbulence model implementation, chemistry model development, interactive input module coding and wing/body flowfield computation are described. The GIM quasi-parabolic code development was completed, and the code used to compute a number of example cases. Turbulence models, algebraic and differential equations, were added to the basic viscous code. An equilibrium reacting chemistry model and implicit finite difference scheme were also added. Development was completed on the interactive module for generating the input data for GIM. Solutions for inviscid hypersonic flow over a wing/body configuration are also presented.

Stainaker, J. F.↗
NASA NTRS1983

Stochastic self-propagating star formation in three-dimensional disk galaxy simulations

The behavior of stochastic self-propagating star formation (SSPSF) in three dimensions is examined in simulation for the first time, emphasizing the effect of the added dimension on the sensitivity of spiral structure to the probability of star formation. The model produces global equilibrium spiral structure over a much more restricted range of star formation probabilities and relaxation times than in two dimensions. Spiral structure also occurs as a transient phenomenon in the runs which eventually fill or evolve to structured nonspiral states. The equilibrium spirals are not as distinctive as those produced by two-dimensional models. However, there are refinements which may modify these results, such as allowing for the depletion of gas in the interstellar medium due to conversion to long-lived low mass stars.

Statler, T.↗
NASA NTRS1988

Development of the general interpolants method for the CYBER 200 series of supercomputers

The General Interpolants Method (GIM) is a 3-D, time-dependent, hybrid procedure for generating numerical analogs of the conservation laws. This study is directed toward the development and application of the GIM computer code for fluid dynamic research applications as implemented for the Cyber 200 series of supercomputers. An elliptic and quasi-parabolic version of the GIM code are discussed. Turbulence models, algebraic and differential equations, were added to the basic viscous code. An equilibrium reacting chemistry model and an implicit finite difference scheme are also included.

Stalnaker, J. F.↗
NASA NTRS2001

Quasi-Equilibrium States in the Tropics Simulated by a Cloud-Resolving Model: Specific Features and Budget Analysis - Part 1

A series of long-term integrations using the two-dimensional Goddard Cumulus Ensemble (GCE) model were performed by altering imposed environmental components to produce various quasi-equilibrium thermodynamic states. Model results show that the genesis of a warm/wet quasi-equilibrium state is mainly due to either strong vertical wind shear (from nudging) or large surface fluxes (from strong surface winds), while a cold/dry quasi-equilibrium state is attributed to a remarkably weakened mixed-wind shear (from vertical mixing due to deep convection) along with weak surface winds. In general, latent heat flux and net large-scale temperature forcing, the two dominant physical processes, dominate in the beginning stage of the simulated convective systems, then considerably weaken in the final stage, which leads to quasi-equilibrium states. A higher thermodynamic regime is found to produce a larger rainfall amount, as convective clouds are the leading source of rainfall over stratiform clouds even though the former occupy much less area. Moreover, convective clouds are more likely to occur in the presence of strong surface winds (latent heat flux), while stratiform clouds (especially the well-organized type) are favored in conditions with strong wind shear (large-scale forcing). The convective systems, which consist of distinct cloud types due to the variation in horizontal winds, are also found to propagate differently. Accordingly, convective systems with mixed-wind shear generally propagate in the direction of shear, while the system with strong (multidirectional) wind shear propagates in a more complex way. Based on the results from the temperature (Q1) and moisture (Q2) budgets, cloud-scale eddies are found to act as a hydrodynamic 'vehicle' that cascades the heat and moisture vertically. Several other specific features such as atmospheric stability, CAPE, and mass fluxes are also investigated and found to be significantly different between diverse quasi-equilibrium states. Detailed comparisons between the various states are presented.

Shie, C.-L.↗