Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model physics”

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 1,765 records · Page 98

Populations of High-Luminosity Density-Bounded HII Regions in Spiral Galaxies? Evidence and Implications

In this paper we present evidence that the H II regions of high luminosity in disk galaxies may be density bounded, so that a significant fraction of the ionizing photons emitted by their exciting OB stars escape from the regions. The key piece of evidence is the presence, in the Ha luminosity functions (LFs) of the populations of H iI regions, of glitches, local sharp peaks at an apparently invariant luminosity, defined as the Stromgren luminosity Lstr), LH(sub alpha) = Lstr = 10(sup 38.6) (+/- 10(sup 0.1)) erg/ s (no other peaks are found in any of the LFs) accompanying a steepening of slope for LH(sub alpha) greater than Lstr This behavior is readily explicable via a physical model whose basic premises are: (a) the transition at LH(sub alpha) = Lstr marks a change from essentially ionization bounding at low luminosities to density bounding at higher values, (b) for this to occur the law relating stellar mass in massive star-forming clouds to the mass of the placental cloud must be such that the ionizing photon flux produced within the cloud is a function which rises more steeply than the mass of the cloud. Supporting evidence for the hypothesis of this transition is also presented: measurements of the central surface brightnesses of H II regions for LH(sub alpha) less than Lstr are proportional to L(sup 1/3, sub H(sub alpha)), expected for ionization bounding, but show a sharp trend to a steeper dependence for LH(sub alpha) greater than Lstr, and the observed relation between the internal turbulence velocity parameter, sigma, and the luminosity, L, at high luminosities, can be well explained if these regions are density bounded. If confirmed, the density-bounding hypothesis would have a number of interesting implications. It would imply that the density-bounded regions were the main sources of the photons which ionize the diffuse gas in disk galaxies. Our estimates, based on the hypothesis, indicate that these regions emit sufficient Lyman continuum not only to ionize the diffuse medium, but to cause a typical spiral to emit significant ionizing flux into the intergalactic medium. The low scatter observed in Lstr, less than 0.1 mag rms in the still quite small sample measured to date, is an invitation to widen the data base, and to calibrate against primary standards, with the aim of obtaining a precise, approx. 10(exp 5) solar luminosity widely distributed standard candle.

Beckman, J. E.↗

Is the Narrow E-Peak Distribution of Gamma-Ray Bursts Real?

Over the performance period of the research grant, the authors conducted a study of the role that the detector response plays in the detection of gamma-ray bursts. The goal of the study was to determine whether the fact that the gamma-ray bursts observed by the BATSE instrument on the Compton Gamma-ray Observatory are characterized by approximately the same characteristic energy is a consequence of the instrument's characteristics, or whether the distribution is a physical attribute of gamma-ray bursts. The authors succeeded in showing that instrumental effects are mild, and that the observed characteristic energy is a physical attribute of bursts. In the course of this research, the authors ported the computer code for calculating the BATSE detector response matrices to the Sun Solaris platform, and created a version of the code that runs under any platform that supports a Fortran 77 compiler with DEC extensions. This code has already been used by other investigators to analyze BATSE data. The authors constructed a Monte Carlo simulation of the BATSE burst trigger, with which they determined the efficiency of detecting a burst as a function of characteristic burst spectral energy. The results were then applied to BATSE observations to determine the physical model for the distribution of burst characteristic energies.

Brainerd, Jerome J.↗

A Historical Perspective on Dynamics Testing at the Langley Research Center

The experience and advancement of Structural dynamics testing for space system applications at the Langley Research Center of the National Aeronautics and Space Administration (NASA) over the past four decades is reviewed. This experience began in the 1960's with the development of a technology base using a variety of physical models to explore dynamic phenomena and to develop reliable analytical modeling capability for space systems. It continued through the 1970's and 80's with the development of rapid, computer-aided test techniques, the testing of low-natural frequency, gravity-sensitive systems, the testing of integrated structures with active flexible motion control, and orbital flight measurements, It extended into the 1990's where advanced computerized system identification methods were developed for estimating the dynamic states of complex, lightweight, flexible aerospace systems, The scope of discussion in this paper includes ground and flight tests and summarizes lessons learned in both successes and failures.

Horta, Lucas G.↗

Genetically Engineered Microelectronic Infrared Filters

A genetic algorithm is used for design of infrared filters and in the understanding of the material structure of a resonant tunneling diode. These two components are examples of microdevices and nanodevices that can be numerically simulated using fundamental mathematical and physical models. Because the number of parameters that can be used in the design of one of these devices is large, and because experimental exploration of the design space is unfeasible, reliable software models integrated with global optimization methods are examined The genetic algorithm and engineering design codes have been implemented on massively parallel computers to exploit their high performance. Design results are presented for the infrared filter showing new and optimized device design. Results for nanodevices are presented in a companion paper at this workshop.

Cwik, Tom↗

User's Manual for LINER: FORTRAN Code for the Numerical Simulation of Plane Wave Propagation in a Lined Two-Dimensional Channel

LINER is a system of Fortran 77 codes which performs a 2D analysis of acoustic wave propagation and noise suppression in a rectangular channel with a continuous liner at the top wall. This new implementation is designed to streamline the usage of the several codes making up LINER, resulting in a useful design tool. Major input parameters are placed in two main data files, input.inc and nurn.prm. Output data appear in the form of ASCII files as well as a choice of GNUPLOT graphs. Section 2 briefly describes the physical model. Section 3 discusses the numerical methods; Section 4 gives a detailed account of program usage, including input formats and graphical options. A sample run is also provided. Finally, Section 5 briefly describes the individual program files.

Reichert, R, S.↗

Improving Global Analysis and Short-Range Forecast Using Rainfall and Moisture Observations Derived from TRMM and SSM/I Passive Microwave Instruments

The Global Precipitation Mission, a satellite project under consideration as a follow-on to the Tropical Rainfall Measuring Mission (TRMM) by the National Aeronautics and Space Agency (NASA) in the United States, the National Space Development Agency (NASDA) in Japan, and other international partners, comprises an improved TRMM-like satellite and a constellation of 8 satellites carrying passive microwave radiometers to provide global rainfall measurements at 3-hour intervals. The success of this concept relies on the merits of rainfall estimates derived from passive microwave radiometers. This article offers a proof-of-concept demonstration of the benefits of using, rainfall and total precipitable water (TPW) information derived from such instruments in global data assimilation with observations from the TRMM Microwave Imager (TMI) and 2 Special Sensor Microwave/Imager (SSM/I) instruments. Global analyses that optimally combine observations from diverse sources with physical models of atmospheric and land processes can provide a comprehensive description of the climate systems. Currently, such data analyses contain significant errors in primary hydrological fields such as precipitation and evaporation, especially in the tropics. We show that assimilating the 6-h averaged TMI and SSM/I surface rainrate and TPW retrievals improves not only the hydrological cycle but also key climate parameters such as clouds, radiation, and the upper tropospheric moisture in the analysis produced by the Goddard Earth Observing System (GEOS) Data Assimilation System, as verified against radiation measurements by the Clouds and the Earth's Radiant Energy System (CERES) instrument and brightness temperature observations by the TIROS Operational Vertical Sounder (TOVS) instruments. Typically, rainfall assimilation improves clouds and radiation in areas of active convection, as well as the latent heating and large-scale motions in the tropics, while TPW assimilation leads to reduced moisture biases and improved radiative fluxes in clear-sky regions. Ensemble forecasts initialized with analyses that incorporate TMI and SSM/I rainfall and TPW data also yield better short-range predictions of geopotential heights, winds, and precipitation in the tropics. This study offers a compelling illustration of the potential of using rainfall and TPW information derived from passive microwave instruments to significantly improve the quality of 4-dimensional global datasets for climate analysis and weather forecasting applications.

Hou, Arthur Y.↗

Radiation Physics for Space and High Altitude Air Travel

Galactic cosmic rays (GCR) are of extra-solar origin consisting of high-energy hydrogen, helium, and heavy ions. The GCR are modified by physical processes as they traverse through the solar system, spacecraft shielding, atmospheres, and tissues producing copious amounts of secondary radiation including fragmentation products, neutrons, mesons, and muons. We discuss physical models and measurements relevant for estimating biological risks in space and high-altitude air travel. Ambient and internal spacecraft computational models for the International Space Station and a Mars mission are discussed. Risk assessment is traditionally based on linear addition of components. We discuss alternative models that include stochastic treatments of columnar damage by heavy ion tracks and multi-cellular damage following nuclear fragmentation in tissue.

Cucinotta, F. A.↗

Distribution of Dust from Kuiper Belt Objects

Using an efficient computational approach, we have reconstructed the structure of the dust cloud in the Solar system between 0.5 and 100 AU produced by the Kuiper belt objects. Our simulations offer a 3-D physical model of the 'kuiperoidal' dust cloud based on the distribution of 280 dust particle trajectories produced by 100 known Kuiper belt objects; the resulting 3-D grid consists of 1.9 x 10' cells containing 1.2 x 10" particle positions. The following processes that influence the dust particle dynamics are taken into account: 1) gravitational scattering on the eight planets (neglecting Pluto); 2) planetary resonances; 3) radiation pressure; and 4) the Poynting-Robertson (P-R) and solar wind drags. We find the dust distribution highly non-uniform: there is a minimum in the kuiperoidal dust between Mars and Jupiter, after which both the column and number densities of kuiperoidal dust sharply increase with heliocentric distance between 5 and 10 AU, and then form a plateau between 10 and 50 AU. Between 25 and 45 AU, there is an appreciable concentration of kuiperoidal dust in the form of a broad belt of mostly resonant particles associated with Neptune. In fact, each giant planet possesses its own circumsolar dust belt consisting of both resonant and gravitationally scattered particles. As with the cometary belts simulated in our related papers, we reveal a rich and sophisticated resonant structure of the dust belts containing families of resonant peaks and gaps. An important result is that both the column and number dust density are more or less flat between 10 and 50 AU, which might explain the surprising data obtained by Pioneers 10 & 11 and Voyager that the dust number density remains approximately distance-independent in this region. The simulated kuiperoidal dust, in addition to asteroidal and cometary dust, might represent a third possible source of the zodiacal light in the Solar system.

Gorkavyi, Nick N.↗

Global Geometric Properties of Martian Impact Craters: A Preliminary Assessment Using Mars Orbiter Laser Altimeter (MOLA)

Impact craters on Mars have been used to provide fundamental insights into the properties of the martian crust, the role of volatiles, the relative age of the surface, and on the physics of impact cratering in the Solar System. Before the three-dimensional information provided by the Mars Orbiter Laser Altimeter (MOLA) instrument which is currently operating in Mars orbit aboard the Mars Global Surveyor (MGS), impact features were characterized morphologically using orbital images from Mariner 9 and Viking. Fresh-appearing craters were identified and measurements of their geometric properties were derived from various image-based methods. MOLA measurements can now provide a global sample of topographic cross-sections of martian impact features as small as approx. 2 km in diameter, to basin-scale features. We have previously examined MOLA cross-sections of Northern Hemisphere and North Polar Region impact features, but were unable to consider the global characteristics of these ubiquitous landforms. Here we present our preliminary assessment of the geometric properties of a globally-distributed sample of martian impact craters, most of which were sampled during the initial stages of the MGS mapping mission (i.e., the first 600 orbits). Our aim is to develop a framework for reconsidering theories concerning impact cratering in the martian environment. This first global analysis is focused upon topographically-fresh impact craters, defined here on the basis of MOLA topographic profiles that cross the central cavities of craters that can be observed in Viking-based MDIM global image mosaics. We have considered crater depths, rim heights, ejecta topologies, cross-sectional "shapes", and simple physical models for ejecta emplacement. To date (May, 1999), we have measured the geometric properties of over 1300 impact craters in the 2 to 350 km diameter size interval. A large fraction of these measured craters were sampled with cavity-center cross-sections during the first two months of MGS mapping. Many of these craters are included in Nadine Barlow's Catalogue of Martian Impact Craters, although we have treated simple craters smaller than about 7 km in greater detail than all previous investigations. Additional information is contained in the original extended abstract.

Garvin, J. B.↗

Design and Performance of McRas in SCMs and GEOS I/II GCMs

The design of a prognostic cloud scheme named McRAS (Microphysics of clouds with Relaxed Arakawa-Schubert Scheme) for general circulation models (GCMs) will be discussed. McRAS distinguishes three types of clouds: (1) convective, (2) stratiform, and (3) boundary-layer types. The convective clouds transform and merge into stratiform clouds on an hourly time-scale, while the boundary-layer clouds merge into the stratiform clouds instantly. The cloud condensate converts into precipitation following the auto-conversion equations of Sundqvist that contain a parametric adaptation for the Bergeron-Findeisen process of ice crystal growth and collection of cloud condensate by precipitation. All clouds convect, advect, as well as diffuse both horizontally and vertically with a fully interactive cloud-microphysics throughout the life-cycle of the cloud, while the optical properties of clouds are derived from the statistical distribution of hydrometeors and idealized cloud geometry. An evaluation of McRAS in a single column model (SCM) with the GATE Phase III and 5-ARN CART datasets has shown that together with the rest of the model physics, McRAS can simulate the observed temperature, humidity, and precipitation without many systematic errors. The time history and time mean incloud water and ice distribution, fractional cloudiness, cloud optical thickness, origin of precipitation in the convective anvil and towers, and the convective updraft and downdraft velocities and mass fluxes all show a realistic behavior. Performance of McRAS in GEOS 11 GCM shows several satisfactory features but some of the remaining deficiencies suggest need for additional research involving convective triggers and inhibitors, provision for continuously detraining updraft, a realistic scheme for cumulus gravity wave drag, and refinements to physical conditions for ascertaining cloud detrainment level.

Sud, Yogesh C.↗

The Coronae of Moderate-Mass Giants in the Hertzsprung Gap and the Clump

We have used the Roentgensatellit (ROSAT), the Extreme Ultraviolet Explorer (EUVE), and the Hubble Space Telescope (HST) to measure X-ray and ultraviolet emissions of moderate-mass (APPROX. 2-3 solar mass) giants in the Hertzsprung gap (spectral types early F to mid-G) and the post-helium flash "clump" (approx. G8-K0). Our motivation was to document the evolution of hot coronae (T greater than 10(exp 6)K) along the post-main-sequence trajectories traveled by such stars in order to gain insight concerning the "X-ray deficiency" of the F-GO giants and the strong braking of stellar rotation at the red edge of the Hertzsprung gap. With few exceptions, Hertzsprung gap and clump giants observed by ROSAT PSPC show hot (T approx. 10(exp 7)K) coronal energy distributions, regardless of any X-ray deficiency, EUVE spectra of gap star 31 Com (G0 111) indicate a broad coronal emission measure hump at approx. 10(exp 7.2)K, while the active clump giant beta Ceti (K0 III) displays a sharp peak at approx. 10(exp 6.8)K, as seen previously in the mixed clump/gap binary Capella (alpha Aur: G8 III + G0 III). The gap giants upsilon Peg (F8 III) and 24 UMa (G4 III) have EUV emissions of intermediate temperature (approx. 10(exp 7.0)K). The stars 31 Com, psi(sup 3) Psc (G0 III), and beta Cet exhibit redshifted transition zone (TZ: approx. 10(exp 5)K) lines in HST GHRS spectra, as reported earlier in Procyon (alpha CMi: F5 IV-V) and Capella G0. Such redshifts on the Sun are thought to signify flows in magnetic loops. beta Cas (F2 III)-a rare soft coronal source among the gap stars-displays blueshifts of C iv and 0 iv], although emissions at cooler and hotter temperatures are near the photospheric velocity. The remarkably broad line profiles of the fastest rotating gap giants suggest that the 10(exp 5) K "subcoronal " emission zones extend to h approx. R(sub *) above the photosphere, about 50 scale heights. In contrast to the TZ line redshifts, the upper chromospheric emissions (e.g., Mg II and Si III) of 31 Com and upsilon(sup 3) Psc have blueshifted cores. Blue-asymmetric peaks in the solar Mg a lines are thought to indicate dynamical heating in the chromosphere. Observations of the H(sub I) Ly(alpha) feature of 31 Com taken 9 months apart reveal striking profile changes, reminiscent of those noted previously in the Ly(alpha) blue peak of the Capella G0 star. We used the far-ultraviolet diagnostics, in combination with ROSAT X-ray photometry and EUVE high-excitation line strengths, to constrain physical models of the stellar outer atmospheres. Quasi-static magnetic loops can simulate the empirical coronal emission measures of the giant stars, but the inferred pressures for sensible loop lengths conflict with direct measurements of subcoronal densities. Furthermore, the high rate of emission at approx. 10(exp 5) K cannot be explained by thermal conduction down the legs of hot quasi-static loops. On the other hand, the possible existence of elongated (l - R(sub *) emission structures on the gap giants leads to a speculative scenario to explain the X-ray deficiency. It is based on the increased importance of the dynamical filling phase ("explosive evaporation") of the loop life cycle; conductive cooling, yielding TZ emissions at the footpoints, when the heating is interrupted; and the possibility for transitions between " hot " and " cool " energy balance solutions owing to dynamical suspension and centrifugal trapping of the cooling gas. The long loops might represent a vestigial global " magnetosphere " inherited from the main-sequence phase, which ultimately is disrupted near approx. G0 by the deepening convective envelope and growth of a more solar-like dynamo. Coronal emissions might be boosted temporarily as the X-ray deficiency is removed but soon would be quenched by wind braking previously inhibited by the magnetospheric "dead zone."

Ayres, Thomas R.↗

Computation of Reacting Flows in Combustion Processes

The objective of this research is to develop an efficient numerical algorithm with unstructured grids for the computation of three-dimensional chemical reacting flows that are known to occur in combustion components of propulsion systems. During the grant period (1996 to 1999), two companion codes have been developed and various numerical and physical models were implemented into the two codes.

Keith, Theo G., Jr.↗

Accuracy Assessment of Response Surface Approximations for Supersonic Turbine Design

There is a growing trend to employ CFD tools to supply the necessary information for design optimization of fluid dynamics components/systems. Such results are prone to uncertainties due to reasons including discretization. errors, incomplete convergence of computational procedures, and errors associated with physical models such as turbulence closures. Based on this type of information, gradient-based optimization algorithms often suffer from the noisy calculations, which can seriously compromise the outcome. Similar problems arise from the experimental measurements. Global optimization techniques, such as those based on the response surface (RS) concept are becoming popular in part because they can overcome some of these barriers. However, there are also fundamental issues related to such global optimization technique such as RS. For example, in high dimensional design spaces, typically only a small number of function evaluations are available due to computational and experimental costs. On the other hand, complex features of the design variables do not allow one to model the global characteristics of the design space with simple quadratic polynomials. Consequently a main challenge is to reduce the size of the region where we fit the RS, or make it more accurate in the regions where the optimum is likely to reside. Response Surface techniques using either polynomials or and Neural Network (NN) methods offer designers alternatives to conduct design optimization. The RS technique employs statistical and numerical techniques to establish the relationship between design variables and objective/constraint functions, typically using polynomials. In this study, we aim at addressing issues related to the following questions: (1) How to identify outliers associated with a given RS representation and improve the RS model via appropriate treatments? (2) How to focus on selected design data so that RS can give better performance in regions critical to design optimization? (3) How to combine NN and polynomial techniques for improving the accuracy of the RS model?

Papila, Nilay↗

Information Retrieval from SAGE II and MFRSR Multi-Spectral Extinction Measurements

Direct beam spectral extinction measurements of solar radiation contain important information on atmospheric composition in a form that is essentially free from multiple scattering contributions that otherwise tend to complicate the data analysis and information retrieval. Such direct beam extinction measurements are available from the solar occultation satellite-based measurements made by the Stratospheric and Aerosol Gas Experiment (SAGE II) instrument and by ground-based Multi-Filter Shadowband Radiometers (MFRSRs). The SAGE II data provide cross-sectional slices of the atmosphere twice per orbit at seven wavelengths between 385 and 1020 nm with approximately 1 km vertical resolution, while the MFRSR data provide atmospheric column measurements at six wavelengths between 415 and 940 nm but at one minute time intervals. We apply the same retrieval technique of simultaneous least-squares fit to the observed spectral extinctions to retrieve aerosol optical depth, effective radius and variance, and ozone, nitrogen dioxide, and water vapor amounts from the SAGE II and MFRSR measurements. The retrieval technique utilizes a physical model approach based on laboratory measurements of ozone and nitrogen dioxide extinction, line-by-line and numerical k-distribution calculations for water vapor absorption, and Mie scattering constraints on aerosol spectral extinction properties. The SAGE II measurements have the advantage of being self-calibrating in that deep space provides an effective zero point for the relative spectral extinctions. The MFRSR measurements require periodic clear-day Langley regression calibration events to maintain accurate knowledge of instrument calibration.

Lacis, Andrew A.↗

Charge Generation and Propagation in Igneous Rocks

Resistivity changes, ground potentials, electromagnetic (EM) and luminous signals prior to or during earthquakes have been reported, in addition to ground uplift and tilt, and to changes in the seismic wave propagation parameters. However, no physical model exists that ties these diverse phenomena together. Through time-resolved impacts experiments it has been observed that, when igneous rocks (gabbro, diorite, granite) are impacted at low velocities (approx. 100 m/sec), highly mobile electronic charge carriers are generated, spreading from a small volume near the impact point, causing electric potentials, EM and light emission. The rock becomes momentarily conductive. When impacted at higher velocities (approx. 1.5 km/sec), the propagation of the P and S waves is registered through the transient piezoelectric response of quartz. At the same time, the rock volume is filled with mobile charge carriers, and a positive surface potential is registered. During the next 1-2 msec the surface potential oscillates, due to electron injection from ground. These observations are consistent with positive holes, e.g. defect electrons in the O(2-) sublattice, that can travel via the O 2p-dominated valence band of the silicate minerals at the speed of a phonon-mediated charge transfer. Before activation, the positive hole charge carriers lay dormant in form of positive hole pairs, PHP, electrically inactive, chemically equivalent to peroxy links in the structures of constituent minerals. PHPs are introduced by way of hydroxyl (O3Si-OH) incorporated into nominally anhydrous minerals when they crystallize in water-laden environments. Given that sound waves of even relatively low intensity appear to cause PHPs dissociation, thus generating mobile positive holes, it is proposed that microfracturing during rock deformation cause PHP dissociation. Depending on where and how much the rock volume is stressed, the positive holes are expected to form fluctuating charge clouds in the earthquake source region that may account for earthquake-related electrical phenomena and the reported low frequency EM signals.

Freund, Friedemann↗

On the Dynamics of Some Discretizations of Convection-Diffusion Equations

Numerical discretizations of differential equations which model physical processes can possess dynamics quite different from that of the equations themselves. Recently the emphasis has been on the the dynamics of numerical discretizations for Ordinary Differential Equations (ODEs). For Partial Differential Equations (PDEs) using a method of lines approach the situation is more complex. First, the spatial discretisation may introduce dynamics not present in the original equations; second, the solution of the resulting system of ODEs is open to the modified dynamics of the ODE solver used. These two effects may interact in a complex manner. In this talk we present some results of our recent work on the dynamics of discretizations of convection-diffusion equations, including those produced using Total Variation Diminishing (TVD) schemes and adaptive grid techniques. A more general overview of the area may be found on our accompanying poster presentation.

Sweby, Peter K.↗

Numerical Simulation of Two Phase Flows

Two phase flows can be found in broad situations in nature, biology, and industry devices and can involve diverse and complex mechanisms. While the physical models may be specific for certain situations, the mathematical formulation and numerical treatment for solving the governing equations can be general. Hence, we will require information concerning each individual phase as needed in a single phase. but also the interactions between them. These interaction terms, however, pose additional numerical challenges because they are beyond the basis that we use to construct modern numerical schemes, namely the hyperbolicity of equations. Moreover, due to disparate differences in time scales, fluid compressibility and nonlinearity become acute, further complicating the numerical procedures. In this paper, we will show the ideas and procedure how the AUSM-family schemes are extended for solving two phase flows problems. Specifically, both phases are assumed in thermodynamic equilibrium, namely, the time scales involved in phase interactions are extremely short in comparison with those in fluid speeds and pressure fluctuations. Details of the numerical formulation and issues involved are discussed and the effectiveness of the method are demonstrated for several industrial examples.

Liou, Meng-Sing↗

Status Of The UPS Space-Marching Flow Solver

The status of the three-dimensional parabolized Navier-Stokes solver UPS is described. The UPS code, initiated at NASA Ames Research Center in 1986, continues to develop and evolve through application to supersonic and hypersonic flow fields. Hypersonic applications have motivated enhancement of the physical modeling capabilities of the code, specifically real gas modeling, boundary conditions, and turbulence and transition modeling. The UPS code has also been modified to enhance robustness and efficiency in order to be practically used in concert with an optimization code for supersonic transport design. These developments are briefly described along with some relevant results for generic test problems obtained during verification of the enhancements. Included developments and results have previously been published and widely disseminated domestically.

Lawerence, Scott L.↗