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 883 records · Page 49

Plasma diffusion at the magnetopause - The case of lower hybrid drift waves

The diffusion expected from the quasi-linear theory of the lower hybrid drift instability at the earth's magnetopause is recalculated. The resulting diffusion coefficient is marginally large enough to explain the thickness of the boundary layer under quiet conditions, based on observational upper limits for the wave intensities. Thus, one possible model for the boundary layer could involve equilibrium between the diffusion arising from lower hybrid waves and various loss processes.

Treumann, R. A.↗

Estimation of the initial equilibrium constants in the formation of tetragonal lysozyme nuclei

Results are presented from a study of the equilibria, kinetic rates, and the aggregation pathway which leads from a lysozyme monomer crystal to a tetragonal crystal, using dialyzed and recrystallized commercial hen eggwhite lysozyme. Relative light scattering intensity measurements were used to estimate the initial equilibrium constants for undersaturated lysozyme solutions in the tetragonal regime. The K1 value was estimated to be (1-3) x 10 exp 4 L/mol. Estimates of subsequent equilibrium constants depend on the crystal aggregation model chosen or determined. Experimental data suggest that tetragonal lysozyme crystal grows by addition of aggregates preformed in the bulk solution, rather than by monomer addition.

Pusey, Marc L.↗

X-ray emission-line spectra of photoionized plasmas - Density sensitivity of the Fe L-shell series

The circumsource environments of accretion-powered X-ray sources are likely to support relatively dense (greater than 10 exp 11/cu cm) photoionized X-ray emission-line regions. The Fe L-shell ions provide a versatile class of plasma diagnostics in this regime, their multielectron structures resulting in diverse spectral phenomena. Attention is given to the spectral response of Fe L-shell ions to variations in electron density over the range 10 exp 11 to 10 exp 16/cu cm. It is found that density-sensitive line ratios exist in the wavelength interval 12-17 A for the ions Fe XVII-XXI. The prominent role of radiative recombination in the population kinetics distinguishes the density-sensitive Fe lines in photoionized plasmas from those which operate in coronal equilibrium plasmas. The results of detailed atomic modeling of these ions are presented and applications to spectroscopic observations of accretion-driven X-ray sources are discussed.

Liedahl, Duane A.↗

The role of large-scale eddies in the climate equilibrium. I - Fixed static stability

An efficient two-level model on a sphere that is based on the balance equations with fixed static stability is developed and used to study how eddies arising from baroclinic instability interact with the temperature structure. The model gives a much better simulation of the eddy momentum flux and of the total eddy forcing of the zonal-mean temperature and zonal wind fields than do quasi-geostrophic beta-plane models. Nonetheless, the results are qualitatively similar. The midlatitude eddy regimes range between two extreme cases. In one, the eddies have no effect on the temperature and zonal wind fields, and in the other (similar to the observed atmosphere), the eddy forcing of the temperature and zonal wind fields is dominated by the eddy heat flux. Quantitatively, some of the model's results differ significantly from those based on the quasi-geostrophic beta-plane. For example, the temperature structure is much more sensitive to the external forcing, and the eddy heat flux is less sensitive to the temperature structure.

Zhou, Shuntai↗

Experimental investigation of nozzle/plume aerodynamics at hypersonic speeds

Much of the work involved the Ames 16-Inch Shock Tunnel facility. The facility was reactivated and upgraded, a data acquisition system was configured and upgraded several times, several facility calibrations were performed and test entries with a wedge model with hydrogen injection and a full scramjet combustor model, with hydrogen injection, were performed. Extensive CFD modeling of the flow in the facility was done. This includes modeling of the unsteady flow in the driver and driven tubes and steady flow modeling of the nozzle flow. Other modeling efforts include simulations of non-equilibrium flows and turbulence, plasmas, light gas guns and the use of non-ideal gas equations of state. New experimental techniques to improve the performance of gas guns, shock tubes and tunnels and scramjet combustors were conceived and studied computationally. Ways to improve scramjet engine performance using steady and pulsed detonation waves were also studied computationally. A number of studies were performed on the operation of the ram accelerator, including investigations of in-tube gasdynamic heating and the use of high explosives to raise the velocity capability of the device.

Bogdanoff, David W.↗

The Diffuse Extreme Ultraviolet Background

Observations of the diffuse EUV background towards 138 different directions using the spectrometers aboard the Extreme Ultraviolet Explorer satellite (EUVE) have been combined into a spectrum from 150A to 730A and represent an effective exposure of 18 million seconds. There is no significant evidence of any non-local line flux in the resultant spectrum such as that from a hot coronal plasma. These results are inconsistent with the Wisconsin C and B broad-band surveys assuming the source is a logT = 5.8 - 6.1 hot plasma in ionization equilibrium with solar abundances, confirming the previous result of Jelinksy, Vallerga and Edelstein) (hereafter Paper 1) using an observation along the ecliptic with the same instrument. To make these results consistent with the previous broad-band surveys, the plasma responsible for the emission must either be depleted in Fe by a factor of approximately 6, be behind an absorbing slab of neutral H with a column of 2 x 10(exp 19)/sq cm, or not be in collisional ionization equilibrium (CIE). One such non-CIE model (Breitswerdt and Schmutzier) that explains the soft x-ray results is also inconsistent with this EUV data.

Vallerga, John↗

Aerothermodynamic Analysis of Stardust Sample Return Capsule with Coupled Radiation and Ablation

An aerothermodynamic analysis of the forebody aeroshell of the Stardust Sample Return Capsule is carried out by using the axisymmetric viscous shock-layer equations with and without fully coupled radiation and ablation. Formulation of the viscous shock-layer equations with shoulder radius as the length scale and implementation of the Vigneron pressure condition allow resolution of the flowfield over the shoulder. With a predominantly supersonic outflow over the shoulder, a globally iterated solution or viscous shock-layer equations can be obtained. The stagnation-point results are obtained along a specified trajectory, whereas detailed calculations along the body are provided at the peak-heating point. The equilibrium calculations with ablation injection are the focus of the present study because of the lack of a general chemical nonequilibrium analysis that accounts for both surface and flowfield effect. The equilibrium calculations also provide a simple way to conserve surface (and flowfield) elemental composition for the current small ablation injection rates, where the surface elemental composition is a mixture of freestream and ablator elements. Therefore, the coupled laminar and turbulent flow solutions with radiation and ablation are obtained by using the equilibrium flow chemistry, whereas a nonequilibrium chemistry model is used for solutions without ablation and turbulence. Various computed results are compared with those obtained by the other researchers.

Gupta, Roop N.↗

X-38 Seal Development

An X-38 Crew Return Vehicle Seal Development is presented. The contents include: 1) X-38 Crew Return Vehicle; 2) X-38 TPS Configuration; 3) X-38 Seal Locations; 4) X-38 Rudder/Fin Seal Assembly; 5) Baseline X-38 Rudder/Fin Seal Design; 6) Rudder/Fin Seal to Bracket Assembly; 7) X-38 Rudder/Fin Vertical Rub Surface Inconel-0.10 inches; 8) X-38 Rudder/Fin Seal Analysis; 9) Seal Analysis Model; and 10) Governing Differential Equations for Equilibrium Thermal Assumption. The X-38 Rudder/Fin Seal temperature and pressure properties are also given.

Curry, Donald M.↗

Afterbody Heating Predictions for a Mars Science Laboratory Entry Vehicle

The Mars Science Laboratory mission intends to deliver a large rover to the Martian surface within 10 km of its target site. One candidate entry vehicle aeroshell consists of a 3.75-m diameter, 70-deg sphere-cone forebody and a biconic afterbody similar to that of Viking. This paper presents computational fluid dynamics predictions of laminar afterbody heating rates for this configuration and a 2010 arrival at Mars. Computational solutions at flight conditions used an 8-species Mars gas model in chemical and thermal non-equilibrium. A grid resolution study examined the effects of mesh spacing on afterbody heating rates and resulted in grids used for heating predictions on a reference entry trajectory. Afterbody heating rate reaches its maximum value near 0.6 W/sq cm on the first windward afterbody cone at the time of peak freestream dynamic pressure. Predicted afterbody heating rates generally are below 3% of the forebody laminar nose cap heating rate throughout the design trajectory. The heating rates integrated over time provide total heat load during entry, which drives thermal protection material thickness.

Edquist, Karl T.↗

Viking Afterbody Heating Computations and Comparisons to Flight Data

Computational fluid dynamics predictions of Viking Lander 1 entry vehicle afterbody heating are compared to flight data. The analysis includes a derivation of heat flux from temperature data at two base cover locations, as well as a discussion of available reconstructed entry trajectories. Based on the raw temperature-time history data, convective heat flux is derived to be 0.63-1.10 W/sq cm for the aluminum base cover at the time of thermocouple failure. Peak heat flux at the fiberglass base cover thermocouple is estimated to be 0.54-0.76 W/sq cm, occurring 16 seconds after peak stagnation point heat flux. Navier-Stokes computational solutions are obtained with two separate codes using an 8-species Mars gas model in chemical and thermal non-equilibrium. Flowfield solutions using local time-stepping did not result in converged heating at either thermocouple location. A global time-stepping approach improved the computational stability, but steady state heat flux was not reached for either base cover location. Both thermocouple locations lie within a separated flow region of the base cover that is likely unsteady. Heat flux computations averaged over the solution history are generally below the flight data and do not vary smoothly over time for both base cover locations. Possible reasons for the mismatch between flight data and flowfield solutions include underestimated conduction effects and limitations of the computational methods.

Edquist, Karl T.↗

Viking Afterbody Heating Computations and Comparisons to Flight Data

Computational fluid dynamics predictions of Viking Lander 1 entry vehicle afterbody heating are compared to flight data. The analysis includes a derivation of heat flux from temperature data at two base cover locations, as well as a discussion of available reconstructed entry trajectories. Based on the raw temperature-time history data, convective heat flux is derived to be 0.63-1.10 W/cm2 for the aluminum base cover at the time of thermocouple failure. Peak heat flux at the fiberglass base cover thermocouple is estimated to be 0.54-0.76 W/cm2, occurring 16 seconds after peak stagnation point heat flux. Navier-Stokes computational solutions are obtained with two separate codes using an 8- species Mars gas model in chemical and thermal non-equilibrium. Flowfield solutions using local time-stepping did not result in converged heating at either thermocouple location. A global time-stepping approach improved the computational stability, but steady state heat flux was not reached for either base cover location. Both thermocouple locations lie within a separated flow region of the base cover that is likely unsteady. Heat flux computations averaged over the solution history are generally below the flight data and do not vary smoothly over time for both base cover locations. Possible reasons for the mismatch between flight data and flowfield solutions include underestimated conduction effects and limitations of the computational methods.

Edquist, Karl T.↗

Role of Climate Change in Global Predictions of Future Tropospheric Ozone and Aerosols

A unified tropospheric chemistry-aerosol model within the Goddard Institute for Space Studies general circulation model II is applied to simulate an equilibrium CO2-forced climate in the year 2100 to examine the effects of climate change on global distributions of tropospheric ozone and sulfate, nitrate, ammonium, black carbon, primary organic carbon, secondary organic carbon, sea salt, and mineral dust aerosols. The year 2100 CO2 concentration as well as the anthropogenic emissions of ozone precursors and aerosols/aerosol precursors are based on the Intergovernmental Panel on Climate Change Special Report on Emissions Scenarios (SRES) A2. Year 2100 global O3 and aerosol burdens predicted with changes in both climate and emissions are generally 5-20% lower than those simulated with changes in emissions alone; as exceptions, the nitrate burden is 38% lower, and the secondary organic aerosol burden is 17% higher. Although the CO2-driven climate change alone is predicted to reduce the global O3 concentrations over or near populated and biomass burning areas because of slower transport, enhanced biogenic hydrocarbon emissions, decomposition of peroxyacetyl nitrate at higher temperatures, and the increase of O3 production by increased water vapor at high NOx levels. The warmer climate influences aerosol burdens by increasing aerosol wet deposition, altering climate-sensitive emissions, and shifting aerosol thermodynamic equilibrium. Climate change affects the estimates of the year 2100 direct radiative forcing as a result of the climate-induced changes in burdens and different climatological conditions; with full gas-aerosol coupling and accounting for ozone and direct radiative forcings by the O2, sulfate, nitrate, black carbon, and organic carbon are predicted to be +0.93, -0.72, -1.0, +1.26, and -0.56 W m(exp -2), respectively, using present-day climate and year 2100 emissions, while they are predicted to be +0.76, -0.72, 0.74, +0.97, and -0.58 W m(exp -2), respectively, with year 2100 climate and emissions.

Liao, Hong↗

Chandra X-ray Grating Spectrometry of Eta Carinae near X-ray Minimum: I. Variability of the Sulfur and Silicon Emission Lines

We report on variations in important X-ray emission lines in a series of Chandra grating spectra of the supermassive colliding wind binary star eta Car, including key phases around the X-ray minimum/periastron passage in 2003.5. The X-rays arise from the collision of the slow, dense wind of eta Car with the fast, low-density wind of an otherwise hidden companion star. The X-ray emission lines provide the only direct measure of the flow dynamics of the companion's wind along the wind-wind collision zone. We concentrate here on the silicon and sulfur lines, which are the strongest and best resolved lines in the X-ray spectra. Most of the line profiles can be adequately fit with symmetric Gaussians with little significant skewness. Both the silicon and sulfur lines show significant velocity shifts and correlated increases in line widths through the observations. The R = forbidden-to-intercombination ratio from the Si XIII and S XV triplets is near or above the low-density limit in all observations, suggesting that the line-forming region is > 1.6 stellar radii from the companion star, and that the emitting plasma may be in a non-equilibrium state. We show that simple geometrical models cannot simultaneously fit both the observed centroid variations and changes in line width as a function of phase. We show that the observed profiles can be fitted with synthetic profiles with a reasonable model of the emissivity along the wind-wind collision boundary. We use this analysis to help constrain the line formation region as a function of orbital phase, and the orbital geometry. Subject headings: X-rays: stars -stars: early-type-stars: individual (q Car)

Henley, D. B.↗

Molybdenum Valence in Basaltic Silicate Melts

The moderately siderophile element molybdenum has been used as an indicator in planetary differentiation processes, and is particularly relevant to core formation [for example, 1-6]. However, models that apply experimental data to an equilibrium differentiation scenario infer the oxidation state of molybdenum from solubility data or from multivariable coefficients from metal-silicate partitioning data [1,3,7]. Partitioning behavior of molybdenum, a multivalent element with a transition near the J02 of interest for core formation (~IW-2) will be sensitive to changes in JO2 of the system and silicate melt structure. In a silicate melt, Mo can occur in either 4+ or 6+ valence state, and Mo6+ can be either octahedrally or tetrahedrally coordinated. Here we present first XANES measurements of Mo valence in basaltic run products at a range of P, T, and JO2 and further quantify the valence transition of Mo.

Danielson, L. R.↗

Ringberg15: Earth's Climate Sensitivities

To assess gaps in understanding of Earth's climate sensitivities a workshop was organised under the auspices of the WCRP (World Climate Research Programme) Grand Science Challenge on Clouds, Circulation and Climate Sensitivity (Ringberg15). The workshop took place in March 2015 and gathered together over thirty experts from around the world for one week. Attendees each gave short presentations and participated in moderated discussions of specific questions related to understanding Earth's climate sensitivities. Most of the time was focused on understanding of the equilibrium climate sensitivity, defined as the equilibrium near-surface warming associated with a doubling of atmospheric carbon dioxide. The workshop produced nine recommendations, many of them focusing on specific research avenues that could be exploited to advance understanding of climate sensitivity. Many of these dealt, in one fashion or another, with the need to more sharply focus research on identifying and testing story lines for a high (larger than 4 degrees Kelvin) or low (less than 2 degrees Kelvin) equilibrium climate sensitivity. Additionally, a subset of model intercomparison projects (CFMIP (Cloud Feedback Model Intercomparison Project), PMIP (Palaeoclimate Modelling Intercomparison Project), PDRMIP (Precipitation Driver and Response Model Intercomparison Project), RFMIP (Radiative Forcing Model Intercomparison Project) and VolMIP (Volcanic Forcings Model Intercomparison Project)) that have been proposed for inclusion within CMIP were identified as being central to resolving important issues raised at the workshop; for this reason modelling groups were strongly encouraged to participate in these projects. Finally the workshop participants encouraged the WCRP to initiate and support an assessment process lead by the Grand Science Challenge on Clouds, Circulation and Climate Sensitivity on the topic of Earth's Climate Sensitivities, culminating in a report that will be published in 2019, forty years after the seminal report by Charney and co-authors.

Sensitivity↗

Modeling Mg II During Solar Flares. II. Nonequilibrium Effects

To extract the information that the Mg II NUV spectra (observed by the Interface Region Imaging Spectrograph) carry about the chromosphere during solar flares, and to validate models of energy transport via model–data comparison, forward modeling is required. The assumption of statistical equilibrium (SE) is typically used to obtain the atomic level populations from snapshots of flare atmospheres, due to computational necessity. However, it is possible that relying on SE could lead to spurious results. We compare solving the atomic level populations via SE versus a non-equilibrium (NEQ) time-dependent approach. This was achieved using flare simulations from RADYN alongside the minority species version MS_RADYN from which the time-dependent Mg II atomic level populations and radiation transfer were computed in complete frequency redistribution. The impacts on the emergent profiles, light curves, line ratios, and formation heights are discussed. In summary we note that NEQ effects during flares are typically important only in the initial stages and for a short period following the cessation of the energy injection. An analysis of the timescales of ionization equilibrium reveals that for most of the duration of the flare, when the temperatures and densities are sufficiently enhanced, the relaxation timescales are short (τrelax < 0.1 s), so that the equilibrium solution is an adequate approximation. These effects vary with the size of the flare, however. In weaker flares, effects can be more pronounced. We recommend that NEQ effects be considered when possible but that SE is sufficient at most stages of the flare.

Kerr, Graham S.↗

Condensation Front Migration in a Protoplanetary Nebula

Condensation front dynamics are investigated in the mid-solar nebula region. A quasi-steady model of the evolving nebula is combined with equilibrium vapor pressure curves to determine evolutionary condensation fronts for selected species. These fronts are found to migrate inwards from the far-nebula to final positions during a period of 10(exp 7) years. The physical process governing this movement is a combination of local viscous heating and luminescent heating from the central star. Two luminescent heating models are used and their effects on the ultimate radial position of the condensation front are discussed. At first the fronts move much faster than the nebular accretion velocity, but after a time the accreting gas and dust overtakes the slowing condensation front.

Davis, Sanford S.↗

Predictive Modeling of Carbon Ablators

Efforts to build a Predictive Material Modeling (PMM) framework from the micro-scale to the macro-scale are presented in this abstract. To reduce the need for extensive testing, accelerate the design cycle process, and reduce uncertainty margins applied to final designs, NASA is developing simulation and modeling tools that enable characterization of material properties and response to high-enthalpy environments. The Porous Microstructure Analysis (PuMA) code has been developed for computing macroscale (volume averaged) properties of porous materials using microscale images from micro-computed tomography (micro-CT). Microscale modeling requires a realistic representation of a material microstructure; these are obtained either synthetically during the design of the material or through X-ray micro-CT. Volume averaged properties are then used to inform macroscale material response models, such as those implemented in the Porous-material Analysis Toolbox based on OpenFOAM (PATO) software, also actively developed by NASA. The computational model in PATO is a generic heat and mass transfer model for porous reactive materials containing several solid phases and a single gas phase. The detailed chemical interactions occurring between the solid phases and the gas phase are modeled at the pore scale assuming local thermal equilibrium. These tools were developed to efficiently interface with other pre-existing codes such as SPARTA (direct simulation Monte Carlo), DPLR (hypersonic CFD), NEQAIR (radiative transport) and DAKOTA (uncertainty quantification and optimization). Detailed flight data (Mars Science Laboratory [MSL] Entry Descent and Landing Instrument [MEDLI]) is critical for validating these computational tools for NASA applications. Examples of modeling ablative material response using these codes will be presented including 3D simulations of the full-scale heatshield of the MSL capsule. The simulations demonstrate the ability of the modern material response code, PATO, to handle the material response of geometrically complex and large domains, through the use of massively parallel computations.

Thermal Protection Systems↗