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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
NASA NTRS1983

On the local standard of rest

Under the influence of a spiral gravitational field, there should be differences among the mean motions of different types of objects with different dispersion velocities in a spiral galaxy. The old stars with high dispersion velocity should have essentially no mean motion normal to the galactic rotation. On the other hand, young objects and interstellar gas may be moving relative to the old stars at a velocity of a few kilometer per second in both the radial (galacto-centric), and circular directions, depending on the spiral model adopted. Such a velocity is usually referred as the systematic motion or the streaming motion. The conventionally adopted local standard of rest is indeed co-moving with the young objects of the solar vicinity. Therefore, it has a net systematic motion with respect to the circular motion of an equilibrium galactic model, defined by the old stars. Previously announced in STAR as N83-24443

Yuan, C.↗
NASA NTRS1990

A model study of the response of mesospheric ozone to short-term solar ultraviolet flux variations

An investigation is conducted in order to determine the relative importance of several modeled processes in controlling the magnitude and phase of the mesospheric ozone response. A detailed one-dimensional modeling study of the mesospheric ozone response to solar UV flux variations is conducted to remove some of the deficiencies in previous studies. This study is also used to examine specifically the importance of solar zenith angle, self-consistent calculation of water vapor abundance, and temperature feedback with a nonlocal thermodynamic equilibrium radiation model. The photochemical model is described, and the assumptions made for the purpose of comparing model results with the observed ozone response obtained from a statistical analysis of Solar Mesosphere Explorer data (Keating et al., 1987) are discussed. The numerical results for the theoretical ozone response are presented. The results of selected time-dependent calculations are considered to illustrate the degree to which a relatively simple model of the mesosphere is able to capture the major characteristics of the observed response.

Summers, M. E.↗
NASA NTRS1980

Self-gravitating gas flow in barred spiral galaxies

A series of two-dimensional numerical experiments is performed in order to test the response of an isothermal, self-gravitating gas disk to a uniformly rotating, barlike gravitational potential. The barlike potential is an equilibrium stellar model from the n-body calculations of Miller and Smith (1979). In the bar-dominated, central regions of the disk, a gas bar whose phase depends primarily on the location of principal resonances in the disk is formed. This response can be understood in terms of orbit-crowding effects. In the gas-dominated outer regions of the disk, two-armed trailing spiral waves are formed. The local pitch angle of these waves increases with increasing fractional gas mass. These self-gravitating gas waves are not self-sustaining. They are driven from the ends of equilibrium stellar bars, and their phase does not depend on the location of resonances in the disk. The relevance of these self-gravitating waves to observations and models of barred spiral galaxies is discussed. It is concluded that these waves and their associated ringlike structures may be consistent with the morphological distribution of gas features in barred spiral galaxies.

Huntley, J. M.↗
NASA NTRSDate not supplied

Trajectory Engineering with Modular Patched Conics for Entry Systems and TPS (TEMPEST)

Brief Presenter Biography (35 word limit): Bohdan Wesely is an Aerospace Engineer in the Entry Systems and Technology Division at Ames. He has worked on a variety of projects for NASA including integrated TPS (thermal protection system) flight hardware deliveries and testing services for commercial partners. Introduction: TEMPEST is a new trajectory analysis framework that is designed to fill the gap between dedicated flight mechanics tools and aerothermal and TPS sizing tools. The project started as an SJSU master’s thesis and has since evolved into a general conceptual design tool capable of studying a wide variety of entry problems. Development is ongoing in the Entry Systems and Technology Division at NASA ARC. Why TEMPEST: Space missions involving entry into a planetary atmosphere involve a series of unique requirements across multiple disciplines. Whether it is traditional entry descent and landing (EDL), or aerocapture, the vehicle must navigate to its target landing location or orbit state, and the TPS must protect the payload during entry. The design process typically involves iterative handoffs between various flight mechanics, flow solver, and material response level tools. During the early conceptual phase, a wide variety of feasible trajectories are simulated in a Monte Carlo scenario which broadly satisfy the mission or landing requirements. Next, computational fluid dynamics (CFD), direct simulation Monte Carlo (DSMC), and other flow solver analyses are performed at various key trajectory points to generate an aero-database, heating and TPS design requirements also emerge at this stage. At this point, with updated aerodynamics from the various flow solvers, trajectories can be re-run, this in turn can change the required freestream conditions for the CFD tools, and as a project progresses, these analyses converge, and uncertainty is reduced. However, there is always a “hand-off” occurring between two inherently coupled phenomena. Analysis Description: One of the goals with TEMPEST is to use a variety of first principles estimation methods coupled with an atmosphere model to predict vehicle aerothermodynamics across the entire flight regime while propagating a 3 or 6 degree of freedom (DoF) trajectory. Aerodynamics methods include modified Newtonian, Maxwell and Cercignani- Lampis-Lord (CLL) for continuum, transitional, and free molecular flow regimes. Aerothermodynamics include boundary layer and reference enthalpy methods, and Mutation++ for non-equilibrium chemistry modeling. TEMPEST is also capable of stitching multiple trajectory segments together to study mission scenarios like multi-pass aerocapture and aero-gravity assists. Most of the program is implemented in MATLAB using modern system objects, it relies on several C++ shared libraries for supporting tools like Gmsh, the Global Reference Atmospheric Model (GRAM), and Mutation++. The various first principles aerothermal estimation methods are discretized across either a structured axisymmetric panel mesh or an unstructured tri-mesh generated from an open-source tool such as Gmsh, this allows solutions on the same mesh to be compared across tools such as CB-Aero. CFD Coupling. A physics-aware, gaussian process CFD anchoring scheme is proposed to adjust the various first principles methods as a CFD database is populated. One goal for this anchoring module is to inform the project where CFD should be run. Full knowledge of the entire trajectory, atmosphere, and aerothermodynamics allows for easier identification of high sensitivity areas and uncertainty quantification. While the first principles effects are well known and proven accurate in existing tools such as CB- Aero and Cart3D, a physics aware CFD anchoring scheme increases tool credibility across a project lifecycle. Material Response Modeling. Correct TPS sizing is critical for optimizing mass for science payloads and ensuring mission success. The process typically involves a thermal analysis along the trajectory with surface heating environments as a boundary condition. Several design constraints are maximum bondline temperature and maximum recession with various margining techniques. The material response tool FIAT, developed out of NASA Ames, is currently being integrated into the TEMPEST environment. TPS recession, shape change, mass loss, and mass property alteration are all factors that can perturb an entry trajectory. For missions like Mars 2020, recession was minimal and was safely handled separately as a post process. For missions such as Jupiter Galileo with a high TPS mass fraction or asteroid entries, recession plays a major role. The proposed fully coupled scheme is to use an epoch-based approach where the trajectory integration is halted after a recession threshold, the energy balance and FIAT are solved at each panel, the mesh, aerodynamics, and mass properties are updated, and the trajectory continues. Several computational tradeoffs have been made during the development of TEMPEST to limit the cost of a single trajectory and preserve its utility as a conceptual, rapid iteration tool. Conclusion: Development of TEMPEST is ongoing and the project is still in its infancy. This talk aims to showcase its unique capabilities to support future NASA entry systems missions.

Bohdan O Wesely↗
NASA NTRS2014

Stellar Laboratories: 3. New Ba 5, Ba 6, and Ba 7 Oscillator Strengths and the Barium Abundance in the Hot White Dwarfs G191-B2B and RE 0503-289

Context. For the spectral analysis of high-resolution and high-signal-to-noise (S/N) spectra of hot stars, state-of-the-art non-local thermodynamic equilibrium (NLTE) model atmospheres are mandatory. These are strongly dependent on the reliability of the atomic data that is used for their calculation. Aims. Reliable Ba 5-7 oscillator strengths are used to identify Ba lines in the spectra of the DA-type white dwarf G191−B2B and the DO-type white dwarf RE 0503−289 and to determine their photospheric Ba abundances. Methods. We newly calculated Ba v-vii oscillator strengths to consider their radiative and collisional bound-bound transitions in detail in our NLTE stellar-atmosphere models for the analysis of Ba lines exhibited in high-resolution and high-S/N UV observations of G191−B2B and RE 0503−289. Results. For the first time, we identified highly ionized Ba in the spectra of hot white dwarfs. We detected Ba vi and Ba vii lines in the Far Ultraviolet Spectroscopic Explorer (FUSE) spectrum of RE 0503−289. The Ba vi/Ba vii ionization equilibrium is well reproduced with the previously determined effective temperature of 70 000 K and surface gravity of log g=7.5. The Ba abundance is 3.5 +/- 0.5 × 10(exp−4) (mass fraction, about 23 000 times the solar value). In the FUSE spectrum of G191−B2B, we identified the strongest Ba vii line (at 993.41 Å) only, and determined a Ba abundance of 4.0 +/- 0.5 × 10(exp−6) (about 265 times solar). Conclusions. Reliable measurements and calculations of atomic data are a pre-requisite for stellar-atmosphere modeling. Observed Ba vi-vii line profiles in two white dwarfs' (G191−B2B and RE 0503−289) far-ultraviolet spectra were well reproduced with our newly calculated oscillator strengths. This allowed to determine the photospheric Ba abundance of these two stars precisely.

New Ba V↗
NASA NTRS2014

Stellar Laboratories II. New Zn Iv and Zn v Oscillator Strengths and Their Validation in the Hot White Dwarfs G191-B2B and RE0503-289

Context. For the spectral analysis of high-resolution and high-signal-to-noise (SN) spectra of hot stars, state-of-the-art non-local thermodynamic equilibrium (NLTE) model atmospheres are mandatory. These are strongly dependent on the reliability of the atomic data that is used for their calculation. In a recent analysis of the ultraviolet (UV) spectrum of the DA-type white dwarf G191B2B,21 Zn iv lines were newly identified. Because of the lack of Zn iv data, transition probabilities of the isoelectronic Ge vi were adapted for a first, coarse determination of the photospheric Zn abundance.Aims. Reliable Zn iv and Zn v oscillator strengths are used to improve the Zn abundance determination and to identify more Zn lines in the spectra of G191B2B and the DO-type white dwarf RE 0503289. Methods. We performed new calculations of Zn iv and Zn v oscillator strengths to consider their radiative and collisional bound-bound transitions in detail in our NLTE stellar-atmosphere models for the analysis of the Zn iv v spectrum exhibited in high-resolution and high-SN UV observations of G191B2B and RE 0503289. Results. In the UV spectrum of G191B2B, we identify 31 Zn iv and 16 Zn v lines. Most of these are identified for the first time in any star. We can reproduce well almost all of them at log Zn 5.52 0.2 (mass fraction, about 1.7 times solar). In particular, the Zn iv Zn v ionization equilibrium, which is a very sensitive Teff indicator, is well reproduced with the previously determined Teff 60 000 2000 K and log g 7.60 0.05. In the spectrum of RE 0503289, we identified 128 Zn v lines for the first time and determined log Zn 3.57 0.2 (155 times solar). Conclusions. Reliable measurements and calculations of atomic data are a pre-requisite for stellar-atmosphere modeling. Observed Zn iv and Zn v line profiles in two white dwarf (G191B2B and RE 0503289) ultraviolet spectra were well reproduced with our newly calculated oscillator strengths. This allowed us to determine the photospheric Zn abundance of these two stars precisely.

strengths↗
NASA NTRS2003

Volcanic Origin of Alkali Halides on Io

The recent observation of NaCl (gas) on Io confirms our earlier prediction that NaCl is produced volcanically. Here we extend our calculations by modeling thermochemical equilibrium of O, S, Li, Na, K, Rb, Cs, F, Cl, Br, and I as a function of temperature and pressure in a Pele-like volcanic gas with O/S/Na/Cl/K = 1.518/1/0.05/0.04/0.005 and CI chondritic ratios of the other (as yet unobserved) alkalis and halogens. For reference, the nominal temperature and pressure for Pele is 1760 plus or minus 210 K and 0.01 bars based on Galileo data and modeling.

Schaefer, L.↗
NASA NTRS2018

Large-Eddy Simulation of Axisymmetric Compression Corner Flow

The Wall-Modeled Large Eddy Simulation (WMLES) approach is used to study the interaction of a shock wave with a high Reynolds number turbulent boundary layer. Since the near wall region is modeled, high Reynolds number turbulent flows can be simulated at a moderate computational cost. The case considered is that of an axisymmetric Mach 2.85 turbulent boundary layer over a 30 compression corner. The Reynolds number of the boundary layer upstream of the interaction based on momentum thickness (Re theta = u sub infinity theta/v sub infinity) is ~12,000. The geometry and flow conditions match the experiments of Dunagan et al. (NASA TM 88227, 1986). The simulations were performed using equilibrium and non-equilibrium wall models. The agreement with experiment is encouraging for the finest grid with respect to the separation bubble length, unsteady shock structure and wall pressure distribution. Sensitivity ofWMLES results to grid, wall model, and blockage effects in the tunnel are reported.

Iyer, Prahladh S.↗
NASA NTRSDate not supplied

Thermodynamic Modeling of the Vapor in Equilibrium With Apollo 17 Basalts.

Lunar mare basalts sampled during the Apollo 17 mission provide insight to the chemical evolution of mare basalt magmas and their associated vapor as these samples erupted and cooled rapidly at the lunar surface. We focus on minerals present on vesicle and vug surfaces in Apollo 17 basalts, such as 71036, which was opened as part of the Apollo Next Generation Sample Analysis (ANGSA) program. The mineral assemblage observed includes native Fe, SiO2, and a Mg-P mineral, possibly merrillite [1, 2]. The pet-rographic context of these phases suggests that they may have been deposited directly from the vapor phase after eruption. Here we use thermodynamic modeling to assess the plausibility of the vapor-phase deposition hypothesis.

D. L. Thompson↗
NASA NTRS1997

The role of reafference in recalibration of limb movement control and locomotion

The reafference model has frequently been used to explain spatial constancy during eye and head movements. We have found that its basic concepts also form part of the information processing necessary for the control and recalibration of reaching movements. Reaching was studied in a novel force environment--a rotating room that creates centripetal forces of the type that could someday substitute for gravity in space flight, and Coriolis forces which are side effects of rotation. We found that inertial, noncontacting Coriolis forces deviate the path and endpoint of reaching movements, a finding that shows the inadequacy of equilibrium position models of movement control. Repeated movements in the rotating room quickly lead to normal movement patterns and to a failure to perceive the perturbing forces. The first movements made after rotation stops, without Coriolis forces present, show mirror-image deviations and evoke perception of a perturbing force even though none is present. These patterns of sensorimotor control and adaptation can largely be explained on the basis of comparisons of efference copy, reafferent muscle spindle, and cutaneous mechanoreceptor signals. We also describe experiments on human locomotion using an apparatus similar to that which Mittelstaedt used to study the optomotor response of the Eristalis fly. These results show that the reafference principle relates as well to the perception of the forces acting on and exerted by the body during voluntary locomotion.

NASA Discipline Neuroscience↗
NASA NTRS1994

Nonequilibrium radiation and chemistry models for aerocapture vehicle flowfields

The primary accomplishments of the project were as follows: (1) From an overall standpoint, the primary accomplishment of this research was the development of a complete gasdynamic-radiatively coupled nonequilibrium viscous shock layer solution method for axisymmetric blunt bodies. This method can be used for rapid engineering modeling of nonequilibrium re-entry flowfields over a wide range of conditions. (2) Another significant accomplishment was the development of an air radiation model that included local thermodynamic nonequilibrium (LTNE) phenomena. (3) As part of this research, three electron-electronic energy models were developed. The first was a quasi-equilibrium electron (QEE) model which determined an effective free electron temperature and assumed that the electronic states were in equilibrium with the free electrons. The second was a quasi-equilibrium electron-electronic (QEEE) model which computed an effective electron-electronic temperature. The third model was a full electron-electronic (FEE) differential equation model which included convective, collisional, viscous, conductive, vibrational coupling, and chemical effects on electron-electronic energy. (4) Since vibration-dissociation coupling phenomena as well as vibrational thermal nonequilibrium phenomena are important in the nonequilibrium zone behind a shock front, a vibrational energy and vibration-dissociation coupling model was developed and included in the flowfield model. This model was a modified coupled vibrational dissociation vibrational (MCVDV) model and also included electron-vibrational coupling. (5) Another accomplishment of the project was the usage of the developed models to investigate radiative heating. (6) A multi-component diffusion model which properly models the multi-component nature of diffusion in complex gas mixtures such as air, was developed and incorporated into the blunt body model. (7) A model was developed to predict the magnitude and characteristics of the shock wave precursor ahead of vehicles entering the Earth's atmosphere. (8) Since considerable data exists for radiating nonequilibrium flow behind normal shock waves, a normal shock wave version of the blunt body code was developed. (9) By comparing predictions from the models and codes with available normal shock data and the flight data of Fire II, it is believed that the developed flowfield and nonequilibrium radiation models have been essentially validated for engineering applications.

Carlson, Leland A.↗
NASA NTRS2011

Turbulence Modeling for Shock Wave/Turbulent Boundary Layer Interactions

Accurate aerodynamic computational predictions are essential for the safety of space vehicles, but these computations are of limited accuracy when large pressure gradients are present in the flow. The goal of the current project is to improve the state of compressible turbulence modeling for high speed flows with shock wave / turbulent boundary layer interactions (SWTBLI). Emphasis will be placed on models that can accurately predict the separated region caused by the SWTBLI. These flows are classified as nonequilibrium boundary layers because of the very large and variable adverse pressure gradients caused by the shock waves. The lag model was designed to model these nonequilibrium flows by incorporating history effects. Standard one- and two-equation models (Spalart Allmaras and SST) and the lag model will be run and compared to a new lag model. This new model, the Reynolds stress tensor lag model (lagRST), will be assessed against multiple wind tunnel tests and correlations. The basis of the lag and lagRST models are to preserve the accuracy of the standard turbulence models in equilibrium turbulence, when the Reynolds stresses are linearly related to the mean strain rates, but create a lag between mean strain rate effects and turbulence when nonequilibrium effects become important, such as in large pressure gradients. The affect this lag has on the results for SWBLI and massively separated flows will be determined. These computations will be done with a modified version of the OVERFLOW code. This code solves the RANS equations on overset grids. It was used for this study for its ability to input very complex geometries into the flow solver, such as the Space Shuttle in the full stack configuration. The model was successfully implemented within two versions of the OVERFLOW code. Results show a substantial improvement over the baseline models for transonic separated flows. The results are mixed for the SWBLI assessed. Separation predictions are not as good as the baseline models, but the over prediction of the peak heat flux downstream of the reattachment shock that plagues many models is reduced.

Lillard, Randolph P.↗
NASA NTRSDate not supplied

Wall-modeled LES of Flow Over a Gaussian Bump

We perform wall-modeled large eddy simulations (WMLES) of turbulent flow over a Gaussian-shaped bump geometry, to assess its performance in the accelerating and separation regions of the flow. The flow conditions are based on the ongoing CFD validation experiments of Slotnick [NATO STO-MP-AVT-307, 2019]. The oncoming flow Mach number is 0.176, and two Reynolds numbers are simulated that are about 10000 and 36000 based on boundary layer thickness upstream of the bump. Preliminary Reynolds-averaged Navier-Stokes simulations are first performed to assess the effects of Mach number, Reynolds number, tunnel top and sidewall effects. Finally, WMLES results with an equilibrium wall model will be presented at two Reynolds numbers to assess their performance for this flow by making detailed comparisons with available experimental and higher-fidelity numerical data.

turbulence, wall model, large eddy simulation↗
NASA NTRS2002

Chandra Discovery of Luminous Supersoft X-Ray Sources in M81

A Chandra ACIS-S imaging observation of the nearby galaxy M81 (NGC 3031) reveals 9 luminous soft X-ray sources. The local environments, X-ray spectral properties, and X-ray light curves of the sources are presented and discussed in the context of prevailing physical models for supersoft sources. It is shown that the sample falls within expectations based on population synthesis models taken from the literature though the high observed luminosities (L approx.2e36 to approx.3e38 ergs in the 0.2--2.0-keV band) and equivalent blackbody temperatures (T approx.40 to 80 eV) place the brightest detected M81 objects at the high luminosity end of the class of supersoft sources defined by previous ROSAT and Einstein studies of nearby galaxies. This is interpreted as a natural consequence of the higher sensitivity of Chandra to hotter and more luminous systems. Most of the sources can be explained as canonical supersoft sources, secreting white dwarfs powered by steady surface nuclear burning, with X-ray spectra well-fit by hot white dwarf local thermodynamic equilibrium atmosphere models. An exceptionally bright source is scrutinized in greater detail as its estimated barometric luminosity, L approx. 1.5e39 ergs, greatly exceeds theoretical estimates for supersoft sources. This source may be beyond the stability limit and undergoing a phase of mass outflow under extreme conditions. Alternatively, a model in which the observed X-ray spectrum arises from an accretion disk around a blacklists of mass approx.1200/sqrt(cosi) solar masses (viewed at an inclination angle 1) cannot be excluded.

Swartz, Douglas A.↗
NASA NTRS1995

Wind Variability of B Supergiants: The Rapid Rotator HD 64760 (B0.5 Ib) - No. 1

We present the results of a 6 day time series of observations of the rapidly rotating B0.5 Ib star HD 64760. We point out several reasons why such intermediate luminosity B supergiants are ideal targets for wind variability studies and then present our results that show the following: continuous wind activity throughout the 6 day run with the wind never in steady state for more than a few hr; wind variability very near nu = 0 km sec(exp -1) in the resonance lines from the lower ionization stages (Al III and C II); a distinct correlation between variability in the Si III ; lambda(lambda)1300 triplets, the strong C III (lambda)1247 singlet, and the onset of extremely strong wind activity, suggesting a connection between photospheric and wind activity; long temporal coherence in the behavior of the strong absorption events; evidence for large-scale spatial coherence, implied by a whole scale, simultaneous weakening in the wind absorption over a wide range in velocities; and ionization variability in the wind accompanying the largest changes in the absorption strengths of the wind lines. In addition, modeling of the wind lines provides the following information about the state the wind in HD 64760. The number of structures on the portion of a constant velocity surface occulting the stellar disk at a particular time must be quite small, while the number on the entire constant velocity surface throughout the wind must be large. The escape probability at low velocity is overestimated by a normal beta approx. 1 velocity law, perhaps due to the presence of low-velocity shocks deep in the wind or a shallow velocity gradient at low velocity. Estimates of the ionization structure in the wind indicate that the ionization ratios are not those expected from thermal equilibrium wind models or from an extrapolation of previous O star results. The large observed q(N V)/q(Si IV) ratio is almost certainly due to distributed X-rays, but the level of ionization predicted by distributed X-ray wind models is inconsistent with the predicted mass-loss rate. Thus, it is impossible to reconcile the observed ionization ratios and the predicted mass-loss rate within the framework of the available models.

Massa, Derck↗
NASA NTRS2005

Aqueous Alteration and Hydrogen Generation on Parent Bodies of Unequilibrated Ordinary Chondrites: Thermodynamic Modeling for the Semarkona Composition

Ordinary chondrites are the most abundant class of meteorites that could represent rocky parts of solar system bodies. However, even the most primitive unequilibrated ordinary chondrites (UOC) reveal signs of mild alteration that affected the matrix and peripheral zones of chondrules. Major chemical changes include oxidation of kamacite, alteration of glass, removal of alkalis, Al, and Si from chondrules, and formation of phases enriched in halogens, alkalis, and hydrogen. Secondary mineralogical changes include formation of magnetite, ferrous olivine, fayalite, pentlandite, awaruite, smectites, phosphates, carbonates, and carbides. Aqueous alteration is consistent with the oxygen isotope data for magnetite. The presence of secondary magnetite, Ni-rich metal alloys, and ferrous silicates in UOC implies that H2O was the oxidizing agent. However, oxidation by H2O means that H2 is produced in each oxidative pathway. In turn, production of H2, and its redistribution and possible escape should have affected total pressure, as well as the oxidation state of gas, aqueous and mineral phases in the parent body. Here we use equilibrium thermodynamic modeling to explore water-rock reactions in UOC. The chemical composition of gas, aqueous, and mineral phases is considered.

Zolotov, M. Y.↗
NASA NTRS1985

Three-dimensional viscous-shock-layer applications for the Space Shuttle orbiter

A nonorthogonal three-dimensional viscous-shock-layer code was modified and used to compute the windward surface flowfield over a Shuttle-like configuration. Code modifications were made in the methods used to evaluate shock boundary conditions and in the implementation of the equilibrium air models. Additional viscous terms were also retained in the governing equations. Effects of these modifications are demonstrated. Three-dimensional flowfield solutions obtained with the modified code showed good comparison of surface heating and pressure with actual flight data. An apparently three-dimensional effect, characterized by a drop in surface heating on the Shuttle windward centerline in the wing root area, was predicted by the present method. A preliminary review of Shuttle ground-based experimental and flight data confirmed this trend. The present three-dimensional solutions were also compared with the predictions of several approximate methods. Those comparisons showed the approximate methods to be valid for predicting heating along the Shuttle windward centerline uner both perfect gas and equilibrium air conditions.

Thompson, R. A.↗
NASA NTRS2005

Fast Whole-Engine Stirling Analysis

An experimentally validated approach is described for fast axisymmetric Stirling engine simulations. These simulations include the entire displacer interior and demonstrate it is possible to model a complete engine cycle in less than an hour. The focus of this effort was to demonstrate it is possible to produce useful Stirling engine performance results in a time-frame short enough to impact design decisions. The combination of utilizing the latest 64-bit Opteron computer processors, fiber-optical Myrinet communications, dynamic meshing, and across zone partitioning has enabled solution times at least 240 times faster than previous attempts at simulating the axisymmetric Stirling engine. A comparison of the multidimensional results, calibrated one-dimensional results, and known experimental results is shown. This preliminary comparison demonstrates that axisymmetric simulations can be very accurate, but more work remains to improve the simulations through such means as modifying the thermal equilibrium regenerator models, adding fluid-structure interactions, including radiation effects, and incorporating mechanodynamics.

Dyson, Rodger W.↗