Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parametric model dependence”

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 253 records · Page 14

Hydrodynamic Stability and Magnetic Reconnection in Disks and Stars

The purpose of this grant is to study parametric instability. The simplest example of parametric instability is a harmonic oscillator with a periodic modulation of the spring constant. If the modulation frequency is close to twice the natural frequency of the oscillator, the amplitude of oscillation tends to grow exponentially. The growth rate is proportional to the strength of the modulation, but it also depends upon the closeness to resonance of the two frequencies, and upon natural damping rate or "Q" of the oscillator. Parametric instabilities are very common in physics. A familiar example is a jogger's ponytail--normally a very strongly damped pendulum, it can be destabilized by the variation in effective gravity during the jogger's stride. Observation confirms that the period of the pendulum is half that of the jogger's vertical motion. In astrophysics, parametric instability may occur by external tidal forcing, or by interaction among eigenmodes. In the latter case, an energetic eigenmode may destabilize modes of half its frequency, provided some weak nonlinearity exists to couple them. Under a previous Astrophysical Theory grant (NAGW-2419), the PI discovered a parametric instability of tidally forced disks such as the accretion disks in cataclysmic variables and X ray binaries [2]. The destabilized modes are tightly-wound, incompressible, three-dimensional waves analogous to g-modes and r-modes in stars. Later work has confirmed our analysis [4]. It was hoped that these modes might provide a source of turbulence and angular momentum transport in accretion disks. However, a follow-up investigation of this instability by local numerical simulations, although confirming the analytically estimated growth rates, found negligible angular momentum flux [3]. Other work, partly supported by the ATP, now strongly indicates that the transport mechanism in such disks is magnetohydrodynamic turbulence [6]. Nevertheless, the parametric mechanism may truncate the outer edges of disks in close binaries [2], and it may be important in disks of very low ionization such as protostellar disks, or even cataclysmic-variable disks in quiescence where the MHD mechanism may be ineffective [5]. All analyses up to 1996 were done in a local approximation where the orbital frequency, shear rate, and tidal field were treated as constants. The locally computed growth rate turns out to depend strongly on radius, and it was unclear how to average these local rates to obtain the correct global rate. This is a critical issue for accretion disks in close binaries, because the local growth rate is comparable to the orbital frequency towards the outer edge of the disk but decreases rapidly inwards. Paper #1 examined this issue in a simplified global model where the destabilizing terms vary with position. We found that the global growth rate is essentially equal to the maximum local rate, provided that the latter is smoothed over a radial range equal to the distance that the destabilized wave propagates at its group speed in one growth time. Thus, in an accretion disk, waves would grow rapidly in the outer parts but would propagate both inwards and outwards at a maximum group speed of order the disk thickness divided by the orbital period.

Goodman, Jeremy↗

Coupled Ablation, Heat Conduction, Pyrolysis, Shape Change and Spallation of the Galileo Probe

The Galileo probe enters the atmosphere of Jupiter in December 1995. This paper presents numerical methodology and detailed results of our final pre-impact calculations for the heat shield response. The calculations are performed using a highly modified version of a viscous shock layer code with massive radiation coupled with a surface thermochemical ablation and spallation model and with the transient in-depth thermal response of the charring and ablating heat shield. The flowfield is quasi-steady along the trajectory, but the heat shield thermal response is dynamic. Each surface node of the VSL grid is coupled with a one-dimensional thermal response calculation. The thermal solver includes heat conduction, pyrolysis, and grid movement owing to surface recession. Initial conditions for the heat shield temperature and density were obtained from the high altitude rarefied-flow calculations of Haas and Milos. Galileo probe surface temperature, shape, mass flux, and element flux are all determined as functions of time along the trajectory with spallation varied parametrically. The calculations also estimate the in-depth density and temperature profiles for the heat shield. All this information is required to determine the time-dependent vehicle mass and drag coefficient which are necessary inputs for the atmospheric reconstruction experiment on board the probe.

Milos, Frank S.↗

The kinematics of buoyancy induced mixing

A model problem to study mixing driven by buoyancy flow fields due to a time dependent body force which consists of a steady and an oscillatory component is considered. Flow fields generated by a time dependent body force, which consists of a steady and an oscillatory component from vibrations or g-jitter, mix two fluids inside a cavity by stretching and folding the interface into various morphological patterns. These patterns show two basic structures depending on the dominant components of the body force. These structures are whirls and tendrils, corresponding to the dominance of either the steady or oscillatory component, respectively. A combination of these structures also occurs in the interaction region. In this region the effect of the steady component opposes the oscillatory component and the result is to smooth the interface, and prevent breakup of the interface. For the parametric range considered, two basic flow regimes occur: convective and chaotic. In the convective regime, the morphological patterns are topological. In the chaotic regime, mixing occurs by a repetitive sequence of bubble formation, necking, and breakup which spreads the concentration field throughout the cavity. The length stretch of the interface increases exponentially with time. The dynamical system characteristics show that in the convective regime, at a point in space, the flow field is oscillatory with elliptical phase space trajectories, and its power spectrum indicates that the flow field responds to the corresponding input frequency. Whereas, in the chaotic regime, the flow field is aperiodic, the phase space trajectories show irregular patterns, and the power spectrum shows a broadband distribution.

Duval, Walter M. B.↗

Business Intelligence Modeling in Launch Operations

This technology project is to advance an integrated Planning and Management Simulation Model for evaluation of risks, costs, and reliability of launch systems from Earth to Orbit for Space Exploration. The approach builds on research done in the NASA ARC/KSC developed Virtual Test Bed (VTB) to integrate architectural, operations process, and mission simulations for the purpose of evaluating enterprise level strategies to reduce cost, improve systems operability, and reduce mission risks. The objectives are to understand the interdependency of architecture and process on recurring launch cost of operations, provide management a tool for assessing systems safety and dependability versus cost, and leverage lessons learned and empirical models from Shuttle and International Space Station to validate models applied to Exploration. The systems-of-systems concept is built to balance the conflicting objectives of safety, reliability, and process strategy in order to achieve long term sustainability. A planning and analysis test bed is needed for evaluation of enterprise level options and strategies for transit and launch systems as well as surface and orbital systems. This environment can also support agency simulation .based acquisition process objectives. The technology development approach is based on the collaborative effort set forth in the VTB's integrating operations. process models, systems and environment models, and cost models as a comprehensive disciplined enterprise analysis environment. Significant emphasis is being placed on adapting root cause from existing Shuttle operations to exploration. Technical challenges include cost model validation, integration of parametric models with discrete event process and systems simulations. and large-scale simulation integration. The enterprise architecture is required for coherent integration of systems models. It will also require a plan for evolution over the life of the program. The proposed technology will produce long-term benefits in support of the NASA objectives for simulation based acquisition, will improve the ability to assess architectural options verses safety/risk for future exploration systems, and will facilitate incorporation of operability as a systems design consideration, reducing overall life cycle cost for future systems. The future of business intelligence of space exploration will focus on the intelligent system-of-systems real-time enterprise. In present business intelligence, a number of technologies that are most relevant to space exploration are experiencing the greatest change. Emerging patterns of set of processes rather than organizational units leading to end-to-end automation is becoming a major objective of enterprise information technology. The cost element is a leading factor of future exploration systems.

Bardina, Jorge E.↗

Acoustic and aerodynamic performance investigation of inverted velocity profile coannular plug nozzles

The results of model scale parametric static and wind tunnel aerodynamic performance tests on unsuppressed coannular plug nozzle configurations with inverted velocity profile are discussed. The nozzle configurations are high-radius-ratio coannular plug nozzles applicable to dual-stream exhaust systems typical of a variable cycle engine for Advanced Supersonic Transport application. In all, seven acoustic models and eight aerodynamic performance models were tested. The nozzle geometric variables included outer stream radius ratio, inner stream to outer stream ratio, and inner stream plug shape. When compared to a conical nozzle at the same specific thrust, the results of the static acoustic tests with the coannular nozzles showed noise reductions of up to 7 PNdB. Extensive data analysis showed that the overall acoustic results can be well correlated using the mixed stream velocity and the mixed stream density. Results also showed that suppression levels are geometry and flow regulation dependent with the outer stream radius ratio, inner stream-to-outer stream velocity ratio and inner stream velocity ratio and inner stream plug shape, as the primary suppression parameters. In addition, high-radius ratio coannular plug nozzles were found to yield shock associated noise level reductions relative to a conical nozzle. The wind tunnel aerodynamic tests showed that static and simulated flight thrust coefficient at typical takeoff conditions are quite good - up to 0.98 at static conditions and 0.974 at a takeoff Mach number of 0.36. At low inner stream flow conditions significant thrust loss was observed. Using an inner stream conical plug resulted in 1% to 2% higher performance levels than nozzle geometries using a bent inner plug.

Knott, P. R.↗

Interlaminar stress analysis of dropped-ply laminated plates and shells by a mixed method

A mixed method of approximation based on Reissner's variational principle is developed for the linear analysis of interlaminar stresses in laminated composites, with special interest in laminates that contain terminated internal plies (dropped-ply laminates). Two models are derived, one for problems of generalized plane deformation and the other for the axisymmetric response of shells of revolution. A layerwise approach is taken in which the stress field is assumed with an explicit dependence on the thickness coordinate in each layer. The dependence of the stress field on the thickness coordinate is determined such that the three-dimensional equilibrium equations are satisfied by the approximation. The solution domain is reduced to one dimension by integration through the thickness. Continuity of tractions and displacements between layers is imposed. The governing two-point boundary value problem is composed of a system of both differential and algebraic equations (DAE's) and their associated boundary conditions. Careful evaluation of the system of DAE's was required to arrive at a form that allowed application of a one-step finite difference approximation. A two-stage Gauss implicit Runge-Kutta finite difference scheme was used for the solution because of its relatively high degree of accuracy. Patch tests of the two models revealed problems with solution accuracy for the axisymmetric model of a cylindrical shell loaded by internal pressure. Parametric studies of dropped-ply laminate characteristics and their influence on the interlaminar stresses were performed using the generalized plane deformation model. Eccentricity of the middle surface of the laminate through the ply drop-off was found to have a minimal effect on the interlaminar stresses under longitudinal compression, transverse tension, and in-plane shear. A second study found the stiffness change across the ply termination to have a much greater influence on the interlaminar stresses.

Harrison, Peter N.↗

Effect of afterbody geometry on aerodynamic characteristics of isolated nonaxisymmetric afterbodies at transonic Mach numbers

A parametric study was conducted in the Langley 16-Foot Transonic Tunnel on an isolated nonaxisymmetic fuselage model that simulates a twin-engine fighter. The effects of aft-end closure distribution (top/bottom) nozzle-flap boattail angle versus nozzle-sidewall boattail angle) and afterbody and nozzle corner treatment (sharp or radius) were investigated. Four different closure distributions with three different corner radii were tested. Tests were conducted over a range of Mach numbers from 0.40 to 1.25 and over a range of angles of attack from -3 to 9 degrees. Solid plume simulators were used to simulate the jet exhaust. For a given closure distribution in the range of Mach numbers tested, the sharp-corner nozzles generally had the highest drag, and the 2-in. corner-radius nozzles generally had the lowest drag. The effect of closure distribution on afterbody drag was highly dependent on configuration and flight condition.

Bangert, Linda S.↗

Norton-Thevenin Receptance Coupling (NTRC) as a Payload Analysis Tool

The NASA Engineering and Safety Center (NESC) has funded a study of a new method formulated by NASA Engineers called Norton-Thevenin Receptance Coupling (NTRC) to perform coupled loads analysis (CLA). The problem that NTRC attempts to solve is the dependency of the payload organization to high CLA costs, long schedules, lack of standard capabilities to evaluate multiple configurations and unavailability of loads when needed. NTRC solves the problem by providing a tool that payload developers can use to obtain loads at a fraction of the cost of a CLA at any time that it is required. While NTRC is not intended to replace the formal load cycles performed by the launch vehicle (LV) provider, it will provide the ability to reduce the conservatism in defining preliminary design loads, assess the impact of design changes between formal load cycles, perform trade studies and perform parametric loads analysis where many different design configurations can be evaluated with a minimum amount of data required from the LV provider. NTRC condenses all the necessary information into the launch vehicle to payload/s connection points or boundary degrees of freedom (BD). The launch vehicle model is represented by its impedance at its BDs; its forcing functions are represented by the acceleration at those BDs when the payload is absent and the latter is represented by its impedance at the same BDs. Payload responses are represented by transfer functions of selected response to interface BDs. The methodology has contributed to the Loads and Dynamics discipline advancement and successfully passed Peer Reviews. NTRC is exact in the frequency domain. Time domain replication and accuracy is outstanding. A second phase is envisioned to benchmark the whole set of CLA events for the Agency's most utilized Launch Vehicles, and ready it for operational deployment at NASA.

Kaufman, Daniel↗

Creation of a Rapid High-Fidelity Aerodynamics Module for a Multidisciplinary Design Environment

In the traditional aerospace vehicle design process, each successive design phase is accompanied by an increment in the modeling fidelity of the disciplinary analyses being performed. This trend follows a corresponding shrinking of the design space as more and more design decisions are locked in. The correlated increase in knowledge about the design and decrease in design freedom occurs partly because increases in modeling fidelity are usually accompanied by significant increases in the computational expense of performing the analyses. When running high fidelity analyses, it is not usually feasible to explore a large number of variations, and so design space exploration is reserved for conceptual design, and higher fidelity analyses are run only once a specific point design has been selected to carry forward. The designs produced by this traditional process have been recognized as being limited by the uncertainty that is present early on due to the use of lower fidelity analyses. For example, uncertainty in aerodynamics predictions produces uncertainty in trajectory optimization, which can impact overall vehicle sizing. This effect can become more significant when trajectories are being shaped by active constraints. For example, if an optimal trajectory is running up against a normal load factor constraint, inaccuracies in the aerodynamic coefficient predictions can cause a feasible trajectory to be considered infeasible, or vice versa. For this reason, a trade must always be performed between the desired fidelity and the resources available. Apart from this trade between fidelity and computational expense, it is very desirable to use higher fidelity analyses earlier in the design process. A large body of work has been performed to this end, led by efforts in the area of surrogate modeling. In surrogate modeling, an up-front investment is made by running a high fidelity code over a Design of Experiments (DOE); once completed, the DOE data is used to create a surrogate model, which captures the relationships between input variables and responses into regression equations. Depending on the dimensionality of the problem and the fidelity of the code for which a surrogate model is being created, the initial DOE can itself be computationally prohibitive to run. Cokriging, a modeling approach from the field of geostatistics, provides a desirable compromise between computational expense and fidelity. To do this, cokriging leverages a large body of data generated by a low fidelity analysis, combines it with a smaller set of data from a higher fidelity analysis, and creates a kriging surrogate model with prediction fidelity approaching that of the higher fidelity analysis. When integrated into a multidisciplinary environment, a disciplinary analysis module employing cokriging can raise the analysis fidelity without drastically impacting the expense of design iterations. This is demonstrated through the creation of an aerodynamics analysis module in NASA’s OpenMDAO framework. Aerodynamic analyses including Missile DATCOM, APAS, and USM3D are leveraged to create high fidelity aerodynamics decks for parametric vehicle geometries, which are created in NASA’s Vehicle Sketch Pad (VSP). Several trade studies are performed to examine the achieved level of model fidelity, and the overall impact to vehicle design is quantified.

Aerodynamic Properties↗

The importance of properties in modeling

Casting and welding of superalloys, stainless steel and titanium alloys are processes which can be improved through modeling of heat flow, fluid flow, residual stress development, and microstructural evolution. These simulations require inputs of thermophysical data, some of which involves the partially or totally liquid state. In particular, these processes involve melting, flow in the liquid, and solidification. Modeling of such processes can lead to an improved understanding of defects such as shrinkage, inclusions, cracks, incomplete filling (or penetration), macrosegregation, improper grain structure, and deviations from dimensional specifications. Effective modeling can shorten process development time and improve quality. An approach to these problems is to develop efficient models; validate through correlations with thermal, distortion, and microstructural data; run parametric studies; extract knowledge based rules; and apply to adaptive closed loop control systems. With the appropriate pre- and post-processing, such analyses can be made 'user friendly'. This would include graphical user interfaces as well as realistic images and color maps. In such form, these models can be used for sensitivity analyses, which are useful in defining appropriate sensors and in the development of control strategies. Such modeling can be done at several levels, e.g., the MARO level, modeling large scale phenomena such as heat and fluid flow or material deformation; the MICRO level, modeling the development of dendrites, grains or precipitates; or at the NANO level, modeling point defects, dislocations, stacking faults, etc. There are many computational issues associated with these simulations, e.g. computational efficiency and accuracy. In addition, there are many materials issues, not the least of which is the availability of accurate high temperature thermophysical data for complex alloys. This would include latent heat of fusion, temperature dependent heat capacity and thermal conductivity (for liquid and solid), viscosity, surface tension, thermal expansion, mechanical properties, etc. Preliminary data is frequently gathered from the literature; however, this is often not available for modern alloys. If additional data are required, measurements can be used; however, these are costly, time consuming and can be erroneous due to a lack of testing standards or impure materials. Microstructural predictors can be extracted from thermal information, e.g. cooling rate and thermal gradient; the prediction of microstructure is dependent on solidus and liquidus temperature, mushy zone permeability, the solidification curve, volume changes, phase transformations, alloying effects (such as surface tension or viscosity), mold/metal reactions, metal/environment reactions, etc. Defect maps may be needed to predict the onset of shrinkage, hot cracking or 'freckling'. Constants may be needed for stress relaxation, dendrite coarsening, vaporization, etc. Visualization was used as a tool to better comprehend complex data sets associated with the analysis of directional solidification (including crystal growth) and welding. Examples include not only isotherms, but also cooling rate, growth rate and thermal gradient. The latter two are not single valued scalars, but rather time and space dependent vector fields. Efficient models were developed for both casting and welding to predict heat flow and the relationship to dendrite and grain growth. These codes include many of the non-linear effects, e.g. radiation, which dominate these processes. The home-built FDM code(s) were designed to be useful not only to the scientist, but also to the process engineer. Special output can be requested to compare directly to experimental data. Visualization procedures were developed to visualize critical results, e.g. fusion zone width at the surface opposite that where the arc is applied ('penetration'). Both elaborate and simplified distortion analyses were carried out. It is clear that extensive mechanical property data are critical in order to accurately predict residual stress patterns. A scheme is currently being developed to integrate these modeling tools into a set of control algorithms; however, the success of this approach is critically dependent on the availability of accurate high temperature thermophysical data.

Giamei, A. F.↗

GridTool: A surface modeling and grid generation tool

GridTool is designed around the concept that the surface grids are generated on a set of bi-linear patches. This type of grid generation is quite easy to implement, and it avoids the problems associated with complex CAD surface representations and associated surface parameterizations. However, the resulting surface grids are close to but not on the original CAD surfaces. This problem can be alleviated by projecting the resulting surface grids onto the original CAD surfaces. GridTool is designed primary for unstructured grid generation systems. Currently, GridTool supports VGRID and FELISA systems, and it can be easily extended to support other unstructured grid generation systems. The data in GridTool is stored parametrically so that once the problem is set up, one can modify the surfaces and the entire set of points, curves and patches will be updated automatically. This is very useful in a multidisciplinary design and optimization process. GridTool is written entirely in ANSI 'C', the interface is based on the FORMS library, and the graphics is based on the GL library. The code has been tested successfully on IRIS workstations running IRIX4.0 and above. The memory is allocated dynamically, therefore, memory size will depend on the complexity of geometry/grid. GridTool data structure is based on a link-list structure which allows the required memory to expand and contract dynamically according to the user's data size and action. Data structure contains several types of objects such as points, curves, patches, sources and surfaces. At any given time, there is always an active object which is drawn in magenta, or in their highlighted colors as defined by the resource file which will be discussed later.

Samareh-Abolhassani, Jamshid↗

Research on Orbital Plasma-Electrodynamics (ROPE)

Since the development of probe theory by Langmuir and Blodgett, the problem of current collection by a charged spherically or cylindrically symmetric body has been investigated by a number of authors. This paper overviews the development of a fully three-dimensional particle simulation code which can be used to understand the physics of current collection in three dimensions and can be used to analyze data resulting from the future tethered satellite system (TSS). According to the TSS configurations, two types of particle simulation models were constructed: a simple particle simulation (SIPS) and a super particle simulation (SUPS). The models study the electron transient response and its asymptotic behavior around a three dimensional, highly biased satellite. The potential distribution surrounding the satellite is determined by solving Laplace's equation in the SIPS model and by solving Poisson's equation in the SUPS model. Thus, the potential distribution in space is independent of the density distribution of the particles in the SUPS model but it does depend on the density distribution of the particles in the SUPS model. The evolution of the potential distribution in the SUPS model is described. When the spherical satellite is charged to a highly positive potential and immersed in a plasma with a uniform magnetic field, the formation of an electron torus in the equatorial plane (the plane in perpendicular to the magnetic field) and elongation of the torus along the magnetic field are found in both the SIPS and the SUPS models but the shape of the torus is different. The areas of high potential that exist in the polar regions in the SUPS model exaggerate the elongation of the electron torus along the magnetic field. The current collected by the satellite for different magentic field strengths is investigated in both models. Due to the nonlinear effects present in SUPS, the oscillating phenomenon of the current collection curve during the first 10 plasma periods can be seen (this does not appear in SIPS). From the parametric studies, it appears that the oscillating phenomenon of the current collection curve occurs only when the magnetic field strength is less than 0.2 gauss for the present model.

Wu, S. T.↗

Radiative Forcing by Contrails

A parametric study of the instantaneous radiative impact of contrails is presented using three different radiative transfer models for a series of model atmospheres and cloud parameters. Contrails are treated as geometrically and optically thin plane parallel homogeneous cirrus layers in a static atmospheres The ice water content is varied as a function of ambient temperature. The model atmospheres include tropical, mid-latitude, and subarctic summer and winter atmospheres Optically thin contrails cause a positive net forcing at top of the atmosphere. At the surface the radiative forcing is negative during daytime. The forcing increases with the optical depth and the amount of contrail cover. At the top of the atmosphere a mean contrail cover of 0.1% with average optical depth of 0.2 to 0.5 causes about 0.01 to 0.03 W/m(exp 2)a daily mean instantaneous radiative forcing. Contrails cool the surface during the day and heat the surface during the night, and hence reduce the daily temperature amplitude The net effect depends strongly on the daily variation of contrail cloud cover. The indirect radiative forcing due to particle changes in natural cirrus clouds may be of the same magnitude as the direct one due to additional cover.

Meerkoetter, R.↗

Influence of Cooling Channel Geometry on the Thermal Response in Silicon Nitride Plates Studied

Engine manufacturers are continually attempting to improve the performance and efficiency of internal combustion engines. Usually they raise the operating temperature or reduce the cooling air requirement for the hot section turbine components. However, the success of these attempts depends on finding materials that are lightweight, are strong, and can withstand high temperatures. Ceramics are among the top candidate materials considered for such harsh applications. They hold low-density, high-temperature strength, and thermal conductivity, and they are undergoing investigation as potential materials for replacing nickel-base alloys and superalloys that are currently used for engine hot-section components. Ceramic structures can withstand higher operating temperatures and a harsh combustion environment. In addition, their low densities relative to metals help reduce component mass. The long-term objectives of the High Temperature Propulsion Components (HOTPC) Project are to develop manufacturing technology, thermal and environmental barrier coatings (TBC/EBC), and the analytical modeling capability to predict thermomechanical stresses in minimally cooled silicon nitride turbine nozzle vanes under simulated engine conditions. Two- and three-dimensional finite element analyses with TBC were conducted at the NASA Glenn Research Center. Nondestructive evaluation was used to determine processing defects. The study included conducting preliminary parametric analytical runs of heat transfer and stress analyses under steady-state conditions to demonstrate the feasibility of using cooled Si3N4 parts for turbine applications. The influence of cooling-channel shapes (such as circular, square, and ascending-order cooling channels) on cooling efficiency and thermal stresses was investigated. Temperature distributions were generated for all cases considered under both cooling and no-cooling conditions, with air being the cooling medium. The table shows the magnitude of the maximum and minimum temperature obtained for the plates under air cooling. Each channel's cross-sectional shape delivered a different temperature; however, the two-dimensional analyses for circular and square or equal-side rectangular holes produced close results. Moreover, the model of the panel with ascending order cooling channels experienced the lowest temperature. A difference of near 260 C was found among the three cooling-hole configurations investigated. The ascending-order cooling channels arrangement showed superior performance by attaining the lowest temperature (1077 C) in comparison to the circular (1379 C) and square (1343 C) channels for the same cooling-hole size. This indicates that the panel with ascending-order cooling channels is the most suitable configuration regardless of the complexity involved in its manufacture. More details pertaining to this study are reported.

Abdul-Aziz, Ali↗

Polar Heat Flow on Io

Recently, Galileo spacecraft data have revealed Io's polar regions to be much warmer than previously expected. This unexpected development came from Photo-Polarimeter Radiometer (PPR) data which show that the minimum night temperatures are in the range of 90-95 K virtually everywhere on Io. The minimum night temperatures show no dependence upon latitude and, when away from the sunset terminator, they show no dependence upon time of night. This is indeed bizarre behavior for surface units which generally had been assumed to be passive with respect to Io's pervasive volcanism. Night temperatures of 90-95 K at high, polar latitudes are particularly hard to explain. Even assuming infinite thermal inertia, at these latitudes there is insufficient sunlight to support these warm night temperatures. Thus, through the process of elimination of other possibilities, we come to the conclusion that these surfaces are volcanically heated. Taking previously passive units and turning them into new sources of heat flow is a radical departure from previous thermophysical model paradigms. However, the geological interpretation is straight forward. We are simply seeing the effect of old, cool lava flows which cover most of the surface of Io but yet have some heat to radiate. Under these new constraints, we have taken on the challenge of formulating a physical model which quantitatively reproduces all of the observations of Io's thermal emission. In the following we introduce a new parametric model which suffices to identify a previously unrecognized polar component of Io's heat flow.

Veeder, G. J.↗

Estimating Temperature Retrieval Accuracy Associated With Thermal Band Spatial Resolution Requirements for Center Pivot Irrigation Monitoring and Management

This study explores the use of synthetic thermal center pivot irrigation scenes to estimate temperature retrieval accuracy for thermal remote sensed data, such as data acquired from current and proposed Landsat-like thermal systems. Center pivot irrigation is a common practice in the western United States and in other parts of the world where water resources are scarce. Wide-area ET (evapotranspiration) estimates and reliable water management decisions depend on accurate temperature information retrieval from remotely sensed data. Spatial resolution, sensor noise, and the temperature step between a field and its surrounding area impose limits on the ability to retrieve temperature information. Spatial resolution is an interrelationship between GSD (ground sample distance) and a measure of image sharpness, such as edge response or edge slope. Edge response and edge slope are intuitive, and direct measures of spatial resolution are easier to visualize and estimate than the more common Modulation Transfer Function or Point Spread Function. For these reasons, recent data specifications, such as those for the LDCM (Landsat Data Continuity Mission), have used GSD and edge response to specify spatial resolution. For this study, we have defined a 400-800 m diameter center pivot irrigation area with a large 25 K temperature step associated with a 300 K well-watered field surrounded by an infinite 325 K dry area. In this context, we defined the benchmark problem as an easily modeled, highly common stressing case. By parametrically varying GSD (30-240 m) and edge slope, we determined the number of pixels and field area fraction that meet a given temperature accuracy estimate for 400-m, 600-m, and 800-m diameter field sizes. Results of this project will help assess the utility of proposed specifications for the LDCM and other future thermal remote sensing missions and for water resource management.

Ryan, Robert E.↗

Preliminary Sizing Completed for Single- Stage-To-Orbit Launch Vehicles Powered By Rocket-Based Combined Cycle Technology

Single-stage-to-orbit (SSTO) propulsion remains an elusive goal for launch vehicles. The physics of the problem is leading developers to a search for higher propulsion performance than is available with all-rocket power. Rocket-based combined cycle (RBCC) technology provides additional propulsion performance that may enable SSTO flight. Structural efficiency is also a major driving force in enabling SSTO flight. Increases in performance with RBCC propulsion are offset with the added size of the propulsion system. Geometrical considerations must be exploited to minimize the weight. Integration of the propulsion system with the vehicle must be carefully planned such that aeroperformance is not degraded and the air-breathing performance is enhanced. Consequently, the vehicle's structural architecture becomes one with the propulsion system architecture. Geometrical considerations applied to the integrated vehicle lead to low drag and high structural and volumetric efficiency. Sizing of the SSTO launch vehicle (GTX) is itself an elusive task. The weight of the vehicle depends strongly on the propellant required to meet the mission requirements. Changes in propellant requirements result in changes in the size of the vehicle, which in turn, affect the weight of the vehicle and change the propellant requirements. An iterative approach is necessary to size the vehicle to meet the flight requirements. GTX Sizer was developed to do exactly this. The governing geometry was built into a spreadsheet model along with scaling relationships. The scaling laws attempt to maintain structural integrity as the vehicle size is changed. Key aerodynamic relationships are maintained as the vehicle size is changed. The closed weight and center of gravity are displayed graphically on a plot of the synthesized vehicle. In addition, comprehensive tabular data of the subsystem weights and centers of gravity are generated. The model has been verified for accuracy with finite element analysis. The final trajectory was rerun using OTIS (Boeing Corporation's trajectory optimization software package), and the sizing output was incorporated into a solid model of the vehicle using PRO/Engineer computer-aided design software (Parametric Technology Corporation, Waltham, MA).

Roche, Joseph M.↗

Thermal balance of the atmospheres of Jupiter and Uranus

Two-dimensional, radiative-convective-dynamical models of the visible atmospheres of Jupiter and Uranus are presented. Zonally-averaged temperatures and heat fluxes are calculated numerically as functions of pressure and latitude. In addition to radiative heat fluxes, the dynamical heat flux due to large-scale baroclinic eddies is included and is parametrized using a mixing length theory which gives heat fluxes similar to those of Stone. The results for Jupiter indicate that the internal heat flow is non-uniform in latitude and nearly balances the net radiative flux leaving the atmosphere. The thermal emission is found to be uniform in latitude in agreement with Pioneer and Voyager observations. Baroclinic eddies are calculated to transport only a small amount of the meridional heat flow necessary to account for the uniformity of thermal emission with latitude. The bulk of the meridional heat transfer is found to occur very deep in the stable interior of Jupiter as originally proposed by Ingersoll and Porco. The relative importance of baroclinic eddies vs. internal heat flow in the thermal balance of Uranus depends on the ratio of emitted thermal power to absorbed solar power. The thermal balance of Uranus is compared to that of Jupiter for different values of this ratio.

Friedson, A. J.↗