ALBERT

Description: 〈br /〉〈span class="detail_caption"〉Notes: 〈/span〉Kerr, Richard A -- New York, N.Y. -- Science. 2002 May 10;296(5570):1006-8.〈br /〉〈span class="detail_caption"〉Record origin:〈/span〉 〈a href="http://www.ncbi.nlm.nih.gov/pubmed/12004100" target="_blank"〉PubMed〈/a〉

Description: This viewpoint comments on recent advances in understanding the design principles of biological networks. It highlights the surprising discovery of "good-engineering" principles in biochemical circuitry that evolved by random tinkering.〈br /〉〈span class="detail_caption"〉Notes: 〈/span〉Alon, U -- New York, N.Y. -- Science. 2003 Sep 26;301(5641):1866-7.〈br /〉〈span class="detail_caption"〉Author address: 〈/span〉Department of Molecular Cell Biology, Weizmann Institute of Science, Rehovot, Israel 76100. urialon@weizmann.ac.il〈br /〉〈span class="detail_caption"〉Record origin:〈/span〉 〈a href="http://www.ncbi.nlm.nih.gov/pubmed/14512615" target="_blank"〉PubMed〈/a〉

Description: 〈br /〉〈span class="detail_caption"〉Notes: 〈/span〉Bhattacharjee, Yudhijit -- New York, N.Y. -- Science. 2003 Apr 25;300(5619):565.〈br /〉〈span class="detail_caption"〉Record origin:〈/span〉 〈a href="http://www.ncbi.nlm.nih.gov/pubmed/12714717" target="_blank"〉PubMed〈/a〉

Description: Various alternative formulations of the LES equations have been explored in which additional evolution equations for variables such as the acceleration, the subgrid-scale stress tensor, or the subgrid-scale force are explicitly carried. Statistics of the velocity field obtained from the equation for the acceleration are shown to depend strongly on the initial conditions. This feature, which is independent of LES modeling issues, seems to prove that the velocity-acceleration formulation of the Navier-Stokes is not useful for numerical simulation. Equations for the subgrid-scale quantities appear to be much more stable. However, models required by this formulation of the LES problem still require additional study.

Description: Two approaches for the identification of internal gravity waves in sheared and unsheared homogeneous stratified turbulence are investigated. First, the phase angle between the vertical velocity and density fluctuations is considered. It is found, however, that a continuous distribution of the phase angle is present in weakly and strongly stratified flow. Second, a projection onto the solution of the linearized inviscid equations of motion of unsheared stratified flow is investigated. It is found that a solution of the fully nonlinear viscous Navier-Stokes equations can be represented by the linearized inviscid solution. The projection yields a decomposition into vertical wave modes and horizontal vortical modes.

Description: The bubble coalescence and interfacial instabilities that are important to modeling critical heat flux (CHF) in reduced-gravity systems can be sensitive to even minute body forces. Understanding these complex phenomena is vital to the design and safe implementation of two-phase thermal management loops proposed for space and planetary-based thermal systems. While reduced gravity conditions cannot be accurately simulated in 1g ground-based experiments, such experiments can help isolate the effects of the various forces (body force, surface tension force and inertia) which influence flow boiling CHF. In this project, the effects of the component of body force perpendicular to a heated wall were examined by conducting 1g flow boiling experiments at different orientations. FC-72 liquid was boiled along one wall of a transparent rectangular flow channel that permitted photographic study of the vapor-liquid interface at conditions approaching CHF. High-speed video imaging was employed to capture dominant CHF mechanisms. Six different CHF regimes were identified: Wavy Vapor Layer, Pool Boiling, Stratification, Vapor Counterflow, Vapor Stagnation, and Separated Concurrent Vapor Flow. CHF showed great sensitivity to orientation for flow velocities below 0.2 m/s, where very small CHF values where measured, especially with downflow and downward-facing heated wall orientations. High flow velocities dampened the effects of orientation considerably. Figure I shows representative images for the different CHF regimes. The Wavy Vapor Layer regime was dominant for all high velocities and most orientations, while all other regimes were encountered at low velocities, in the downflow and/or downward-facing heated wall orientations. The Interfacial Lift-off model was modified to predict the effects of orientation on CHF for the dominant Wavy Vapor Layer regime. The photographic study captured a fairly continuous wavy vapor layer travelling along the heated wall while permitting liquid contact only in wetting fronts, located in the troughs of the interfacial waves. CHF commenced when wetting fronts near the outlet were lifted off the wall. The Interfacial Lift-off model is shown to be an effective tool for predicting the effects of body force on CHF at high velocities.

Description: Probabilistic CFD design is needed because we are asked to do more with less. To cost effectively accomplish the design task, we need to formally quantify the effect of uncertainties (variables) in the design. Probabilistic design is one effective method to formally quantify the effect of uncertainties. Our objective is to establish a revolutionary new early design process, by developing non-deterministic physics-based probabilistic design tools, which will include all the life cycle processes. Breakthroughs will be sought in speed, accuracy, intelligence, and usability of the system. This paper is concerned with the usefulness of parametric optimization method coupled with a Navier-Stokes analysis code for the aero-thermodynamic design of turbomachinery combustor liner. The interconnection between the CFD code and NESSUS codes facilitated the coupling between the thermal profiles and structural design. We have developed new concepts for reducing the computational cost of unsteady, three-dimensional, compressible aerodynamic analyses for multistage turbomachinery flows. The flow was modeled by the three-dimensional Favre-Reynolds-averaged Navier-Stokes equations using the k-epsilon turbulence closure, which was integrated using an implicit third-order upwind solver. The methodology developed in this paper is expected to lead to the design optimization of turbomachinery blades.

Description: Parallel Plate Plastometer (PPP) is a device commonly used for measuring the viscosity of high polymers at low rates of shear in the range 10(exp 4) to 10(exp 9) poises. This device is being validated for use in measuring the viscosity of liquid glasses at high temperatures having similar ranges for the viscosity values. PPP instrument consists of two similar parallel plates, both in the range of 1 inch in diameter with the upper plate being movable while the lower one is kept stationary. Load is applied to the upper plate by means of a beam connected to shaft attached to the upper plate. The viscosity of the fluid is deduced from measuring the variation of the plate separation, h, as a function of time when a specified fixed load is applied on the beam. Operating plate speeds measured with the PPP is usually in the range of 10.3 cm/s or lower. The flow field within the PPP can be simulated using the equations of motion of fluid flow for this configuration. With flow speeds in the range quoted above the flow field between the two plates is certainly incompressible and laminar. Such flows can be easily simulated using numerical modeling with computational fluid dynamics (CFD) codes. We present below the mathematical model used to simulate this flow field and also the solutions obtained for the flow using a commercially available finite element CFD code.

Description: In the design of a combustor, information is necessary for the mixing of the fuel and air in order to determine the optimum combustor length. In scramjet combustors the mixing often takes place in a shear layer that is formed between the fuel and air. This research was an experimental study of shear layers in supersonic flows aimed at determining what mechanisms affect the shear layer so that the mixing could be better predicted. A second goal was to provide sufficient instream information for use in checking existing Computational Fluid Dynamic (CFD) codes. The shear layer between a supersonic two-dimensional air stream (M = 2 or M = 3) was mixed with a near sonic two-dimensional air stream (M = 1.2). Instream measurements of pitot pressure and cone static pressure were used to determine mean velocity profiles at various axial locations. These velocity profiles were used to determine the shear layer spreading rate and are compared with various predictions. Wall measurements of static pressure, temperature and skin friction were also taken and are presented. The instream measurements were also used for comparison with an existing CFD code. The upstream velocity, pressure and temperature profiles were used as a starting profile and the code was used to calculate downstream profiles for comparison with the experimental results. Reasonable agreement between the measured and calculated results was obtained.

Description: The Physics of Colloids in Space (PCS) experiment was accommodated within International Space Station (ISS) EXpedite the PRocessing of Experiments to Space Station (EXPRESS) Rack 2 and was remotely operated from early June 2001 until February 2002 from NASA Glenn Research Center's Telescience Support Center in Cleveland, Ohio and from a remote site at Harvard University in Cambridge, Massachusetts. PCS is an experiment conceived by Professor David A. Weitz of Harvard University (the Principal Investigator), focusing on the behavior of three different classes of colloid mixtures. The sophisticated light scattering instrumentation comprising PCS is capable of color imaging, and dynamic and static light scattering from 11 to 169 degrees, Bragg scattering over the range from 10 to 60 degrees, and laser light scattering at low angles from 0.3 to 6.0 degrees. The PCS instrumentation performed remarkably well, demonstrating a flexibility that enabled experiments to be performed that had not been envisioned prior to launch. While on-orbit, PCS accomplished 2400 hours of science operations, and was declared a resounding success. Each of the eight sample cells worked well and produced interesting and important results. Crystal nucleation and growth and the resulting structures of two binary colloidal crystal alloys were studied, with the long duration microgravity environment of the ISS facilitating extended studies on the growth and coarsening characteristics of the crystals. In another experiment run, the de-mixing of the colloid-polymer critical-point sample was studied as it phase-separates into two phases, one that resembles a gas and one that resembles a liquid. This process was studied over four decades of length scale, from 1 micron to 1 centimeter, behavior that cannot be observed in this sample on Earth because sedimentation would cause the colloids to fall to the bottom of the cell faster than the de-mixing process could occur. Similarly, the study of gelation and aging of another colloid-polymer sample, the colloid-polymer gel, also provided valuable information on gelation mechanisms, as did investigations on the extremely the low concentration silica and polystyrene fractal gel samples.

Description: An explicit time filter is applied to the Navier-Stokes equation prior to a space filter. The time filter is supposed to be smooth, and an exact expansion depending on the time derivatives of the velocity is derived for the associated stress tensor. On the contrary, the effect of the space filter is treated as usual and an eddy viscosity model is introduced in the LES equation. The total stress is thus represented using a new class of mixed models combining time and space derivatives of the LES field.

Description: As pointed out by Rodi standard integral solutions for jets and plumes developed for discharge into infinite, quiescent ambient are difficult to extend to complex situations, particularly in the presence of boundaries such as the sea floor or ocean surface. In such cases the assumption of similarity breaks down and it is impossible to find a suitable entrainment coefficient. The models are also incapable of describing any but the most slowly varying unsteady motions. There is therefore a need for full time-dependent modeling of the flow field for which there are three main approaches: (1) Reynolds averaged numerical simulation (RANS), (2) large eddy simulation (LES), and (3) direct numerical simulation (DNS). Rodi applied RANS modeling to both jets and plumes with considerable success, the test being a match with experimental data for time-averaged velocity and temperature profiles as well as turbulent kinetic energy and rms axial turbulent velocity fluctuations. This model still relies on empirical constants, some eleven in the case of the buoyant jet, and so would not be applicable to a partly laminar plume, may have limited use in the presence of boundaries, and would also be unsuitable if one is after details of the unsteady component of the flow (the turbulent eddies). At the other end of the scale DNS modeling includes all motions down to the viscous scales. Boersma et al. have built such a model for the non-buoyant case which also compares well with measured data for mean and turbulent velocity components. The model demonstrates its versatility by application to a laminar flow case. As its name implies, DNS directly models the Navier-Stokes equations without recourse to subgrid modeling so for flows with a broad spectrum of motions (high Re) the cost can be prohibitive - the number of required grid points scaling with Re(exp 9/4) and the number of time steps with Re(exp 3/4). The middle road is provided by LES whereby the Navier-Stokes equations are formally filtered with the filter chosen to only exclude the smallest turbulent motions. If successful, LES should provide much of the detail available to DNS but at more bearable cost. Fatica et al. in comparing LES with DNS for a low Reynolds number jet showed that the LES could simulate the temporally evolving behavior including growth of the jet thickness. It is the intention of this report to explore the application of an LES model to jets and plumes. As always, before tackling complex situations, the model must be tested for the simplest of cases and so we address only two, a non-buoyant axisymmetric jet issuing steadily from an orifice into a semi-infinite stationary environment and a buoyant jet in the same environment. The work is a continuation of Basu and Mansour.

Description: This paper compares the filtering used in Coherent Vortex Simulation (CVS) decomposition with an orthogonal wavelet basis, with the Proper Orthogonal Decomposition (POD) or Fourier filtering. Both methods are applied to a field of Direct Numerical Simulation (DNS) data of 3D forced homogeneous isotropic turbulence at microscale Reynolds number R(sub lambda) = 168. We show that, with only 3%N retained modes, CVS filtering separates the coherent vortex tubes from the incoherent background flow. The latter is structureless, has an equipartition energy spectrum, and has a Gaussian velocity probability distribution function (PDF) and an exponential vorticity PDF. On the other hand, the Fourier basis does not extract the coherent vortex tubes cleanly and leaves organized structures in the residual high wavenumber modes whose PDFs are stretched exponentials for both the velocity and the vorticity.

Description: The objective of this study is to investigate the use of the second generation bi-orthogonal wavelet transform for the field decomposition in the Coherent Vortex Simulation of turbulent flows. The performances of the bi-orthogonal second generation wavelet transform and the orthogonal wavelet transform using Daubechies wavelets with the same number of vanishing moments are compared in a priori tests using a spectral direct numerical simulation (DNS) database of isotropic turbulence fields: 256(exp 3) and 512(exp 3) DNS of forced homogeneous turbulence (Re(sub lambda) = 168) and 256(exp 3) and 512(exp 3) DNS of decaying homogeneous turbulence (Re(sub lambda) = 55). It is found that bi-orthogonal second generation wavelets can be used for coherent vortex extraction. The results of a priori tests indicate that second generation wavelets have better compression and the residual field is closer to Gaussian. However, it was found that the use of second generation wavelets results in an integral length scale for the incoherent part that is larger than that derived from orthogonal wavelets. A way of dealing with this difficulty is suggested.

Description: The main goal of this presentation is to give some of the objectives of the testing program. This includes: develop jet noise data base for separate flow nozzles with bypass ratio's 5 to 14; evaluate effect of pylon on noise; develop low performance impact noise suppression concepts; and evaluate potential for active control of jet noise.

Description: Coherent Vortex Simulation (CVS) filtering has been applied to Direct Numerical Simulation (DNS) data of forced and unforced time-developing turbulent mixing layers. CVS filtering splits the turbulent flow into two orthogonal parts, one corresponding to coherent vortices and the other to incoherent background flow. We have shown that the coherent vortices can be represented by few wavelet modes and that these modes are sufficient to reproduce the vorticity probability distribution function (PDF) and the energy spectrum over the entire inertial range. The remaining incoherent background flow is homogeneous, has small amplitude, and is uncorrelated. These results are compared with those obtained for the same compression rate using large eddy simulation (LES) filtering. In contrast to the incoherent background flow of CVS filtering, the LES subgrid scales have a much larger amplitude and are correlated, which makes their statistical modeling more difficult.

Description: This presentation will discuss Ir for online diagnostics with acoustics; Ptot, Ttot, Pstat rake surveys for mean flow measurements; focused Schleiren for density and some turbulence structure; Laser sheet visualization for near-nozzle diagnostics; and Two-point hotwire measurements for turbulence models.

Description: The Lagrangian-averaged Navier-Stokes (LANS) equations are numerically evaluated as a turbulence closure. They are derived from a novel Lagrangian averaging procedure on the space of all volume-preserving maps and can be viewed as a numerical algorithm which removes the energy content from the small scales (smaller than some a priori fixed spatial scale alpha) using a dispersive rather than dissipative mechanism, thus maintaining the crucial features of the large scale flow. We examine the modeling capabilities of the LANS equations for decaying homogeneous turbulence, ascertain their ability to track the energy spectrum of fully resolved direct numerical simulations (DNS), compare the relative energy decay rates, and compare LANS with well-accepted large eddy simulation (LES) models.

Description: This paper presents a set of second-order closure models for low-Reynolds-number turbulence near the wall. Existing closure models for the Reynolds-stress equations were modified to show proper near-wall behavior. A dissipation-rate equation for the turbulent kinetic energy is also reformulated. The proposed models satisfy realizability and will not produce unphysical behavior. Fully developed channel flows are used for model testing. The equations are solved for the mean velocity, the Reynolds stresses, and the dissipation rate of the turbulent kinetic energy. The calculations are compared with both direct numerical simulations and with measurements. It is shown that the present models perform well in predicting the behavior of the turbulence near a wall. Significant improvements over previous models in predicting the components of the Reynolds stress tensor are obtained in the present models.

Description: The overall objective of this work is to study nucleate boiling heat transfer under microgravity conditions in such a way that while providing basic knowledge of the phenomena, it also leads to development of simulation models and correlations that can be used as design tools for a wide range of gravity levels. In the study a building block type of approach is used and both pool and low velocity flow boiling are investigated. Starting with experiments using a single bubble, the complexity of the experiments is increased to two or three inline bubbles, to five bubbles placed on a two-dimensional grid. Finally, experiments are conducted where a large number of prescribed cavities nucleate on the heater and when a commercial surface is used. So far experiments have been conducted at earth normal gravity and in the reduced gravity environment of the KC-135 aircraft whereas experiments on the space station are planned. Modeling/complete numerical simulation of the boiling process is an integral part of the total effort. Experiments conducted with single bubbles formed on a nucleation site microfabricated on a polished silicon wafer show that for gravity levels (g) varying from 1.5g(sub e) to 0.01g(sub e), the bubble diameter at departure varies approximately as (g(sub e)/g)(exp 1/2) and the growth period as (g(sub e)/g). When bubbles merge either inline or in a plane, the bubble diameter at departure is found to be smaller than that obtained for a single bubble and shows a weaker dependence on the level of gravity. The possible reason is that as the bubbles merge they create fluid circulation around the bubbles, which in turn induces a lift force that is responsible for the earlier departure of the bubbles. The verification of this proposition is being sought through numerical simulations. There is a merger of two inline, three inline, and several bubbles in a plane in the low gravity environment of the KC-135 aircraft. After merger and before departure, a mushroom type of bubble with several stems attached to the heater surface is clearly evident. Local heat fluxes during growth and departure of a single bubble were also measured. It was found that during most of the growth period of the bubble, generally the wall heat flux decreased with time because of the increased dry area under the bubble. However, the heat flux increased rapidly just prior to departure of the bubble because of the transient conduction into the cold liquid rushing to fill the space vacated by the bubble as the bubble base shrinks. The measured heat fluxes at various radial locations are found to be in qualitative agreement with the numerical predictions. Single bubble studies at earth normal gravity have also been performed on surfaces oriented at different angles to the gravitational acceleration with flow parallel to the surface. It is found that in all cases the bubbles slide along the surface before lift-off from the surface. The lift force generated as a result of the relative motion between the sliding bubbles and the imposed flow is found to play an important role when the normal force due to buoyancy is reduced. An experimental apparatus for the study of the bubble behavior with imposed flow under reduced gravity conditions has been developed and will soon be employed for experiments in the KC-135 aircraft.

Description: Dynamic wetting, the displacement of one fluid by another immiscible fluid on a surface, controls many natural and technological phenomena, such as coating, printing, spray painting and lubricating. Particularly in coating and spraying applications, contact lines advance across pre-existing fluid films. Most previous work has focused on contact lines advancing across films sufficiently thick that they behave as simple Newtonian fluids. Ultrathin films, where the film thickness may impinge on fundamental length scales in the fluid, have received less attention. In this talk, we will discuss the effects of ultrathin polymer films on dynamic wetting. We measure the interface shape within microns of moving contact lines advancing across preexisting films and compare the measurements to existing models of viscous bending for interfaces advancing across dry surfaces and 'thick' (in the sense that they behave as liquids) films. In the experiments, we advance a contact line of 10-poise and 1-poise polydimethylsiloxane (silicone oil) across pre-coated films of the same fluid with thickness from a single chain thickness (approx. 10 A) through a couple of radii of gyration (100-200 A) to films so thick they are likely bulk in behavior (10(exp 3) A). All films are physisorbed, i.e. they readily rinse from the surface. Thus, molecules in the film are not anchored to the surface and can move within the film if the hydrodynamics dictate such motion. For films of the thickness of a single chain (approx. 10 A), our experiments indicate that the advancing fluid behaves just as it would if it advanced over a dry surface. For the thicker films (10(exp 3) A), we find behavior indicating that the molecules in the film are acting as a fluid with the bulk properties. In this regime, results for the two different fluids are identical when the experiments are performed at the same pre-existing film thickness and advancing capillary number, Ca. For film of thickness of a few radii of gyration (approx. 100-200 A), the behavior depends on Ca of the advancing meniscus. At low Ca, the viscous bending of the interface near the contact line does not behave as it would on a dry surface. It has a lower curvature than expected. However, at higher Ca, the viscous bending is described by the model for spreading over a dry surface. These results show that the fluid flow in the film does behave differently than bulk as the film thickness becomes comparable to molecular length scale. But even more intriguing is the unusual velocity dependence of that behavior where the film behaves more solid-like at higher contact line speeds. We will discuss these results in terms of the properties of confined polymer melts.

Description: This paper presents turbulence modeling from NASA's perspective. The topics include: 1) Hierarchy of Solution Methods; 2) Turbulence Modeling Focus; 3) Linear Eddy Viscosity Models; and 4) Nonlinear Eddy Viscosity Algebraic Stress Models.

Description: This paper presents the purpose and expectations of turbulence modeling, in viewgraph form.

Description: An important component in biotechnology, particularly in the area of protein engineering and rational drug design is the knowledge of the precise three-dimensional molecular structure of proteins. The quality of structural information obtained from X-ray diffraction methods is directly dependent on the degree of perfection of the protein crystals. As a consequence, the growth of high quality macromolecular crystals for diffraction analyses has been the central focus for biochemists, biologists, and bioengineers. Macromolecular crystals are obtained from solutions that contain the crystallizing species in equilibrium with higher aggregates, ions, precipitants, other possible phases of the protein, foreign particles, the walls of the container, and a likely host of other impurities. By changing transport modes in general, i.e., reduction of convection and sedimentation, as is achieved in "microgravity", researchers have been able to dramatically affect the movement and distribution of macromolecules in the fluid, and thus their transport, formation of crystal nuclei, and adsorption to the crystal surface. While a limited number of high quality crystals from space flights have been obtained, as the recent National Research Council (NRC) review of the NASA microgravity crystallization program pointed out, the scientific approach and research in crystallization of proteins has been mainly empirical yielding inconclusive results. We postulate that we can reduce convection in ground-based experiments and we can understand the different aspects of convection control through the use of strong magnetic fields and field gradients. Whether this limited convection in a magnetic field will provide the environment for the growth of high quality crystals is still a matter of conjecture that our research will address. The approach exploits the variation of fluid magnetic susceptibility with concentration for this purpose and the convective damping is realized by appropriately positioning the crystal growth cell so that the magnetic susceptibility force counteracts terrestrial gravity. The general objective is to test the hypothesis of convective control using a strong magnetic field and magnetic field gradient and to understand the nature of the various forces that come into play. Specifically we aim to delineate causative factors and to quantify them through experiments, analysis and numerical modeling. Once the basic understanding is obtained, the study will focus on testing the hypothesis on proteins of pyruvate dehydrogenase complex (PDC), proteins E1 and E3. Obtaining high crystal quality of these proteins is of great importance to structural biologists since their structures need to be determined. Specific goals for the investigation are: 1. To develop an understanding of convection control in diamagnetic fluids with concentration gradients through experimentation and numerical modeling. Specifically solutal buoyancy driven convection due to crystal growth will be considered. 2. To develop predictive measures for successful crystallization in a magnetic field using analyses and numerical modeling for use in future protein crystal growth experiments. This will establish criteria that can be used to estimate the efficacy of magnetic field flow damping on crystallization of candidate proteins. 3. To demonstrate the understanding of convection damping by high magnetic fields to a class of proteins that is of interest and whose structure is as yet not determined. 4. To compare quantitatively, the quality of the grown crystals with and without a magnetic field. X-ray diffraction techniques will be used for the comparative studies. In a preliminary set of experiments, we studied crystal dissolution effects in a 5 Tesla magnet available at NASA Marshall Space Flight Center (MSFC). Using a Schlieren setup, a 1mm crystal of Alum (Aluminum-Potassium Sulfate) was introduced in a 75% saturated solution and the resulting dissolution plume was observed. The experiment was conducted both in the presence and absence of a magnetic field gradient. The magnet produces a gradient field of approx. 1 Tesla2/cm. Image analysis of the recorded images indicated an enhanced plume velocity that was of the order of the measurement limit. For this experiment, both the gradient and gravity fields are in the same direction resulting in an enhanced effective gravity that tends to accelerate the observed plume velocity. While the results are not conclusive, pending further tests, it clearly points out the inadequacy of the MSFC magnet for conducting protein crystallization experiments and the need for a stronger magnet. In spacebased experiments, however, where the gravitational effects are small, only a weak magnetic field will be required to control or mitigate the effects of convective contamination.

Description: The Earth's atmosphere from 90 km to 200 km provides the last aerothermodynamics frontier. Present NASA programs which require but also can provide an understanding of the aerodynamics and aerothermodynamics of the free molecule and transition flows that exist at these altitudes are the Aeroassisted OTV, Entry Research Vehicle and the Tethered Satellite. Each of these programs provides a unique opportunity to do flight research in the rarefied upper atmosphere. However, the Tethered Satellite Program provides, because of its capability to obtain global, in-situ, steady state data, the greatest potential to: (1)define the performance of aerodynamic shapes as a function of environmental characteristics (free molecule, transition, slip flow regimes); (2)define the characteristics of the upper atmosphere and the global variability of properties such as composition temperature, pressure and density. Such data are required to accomplish the systematic development and verification of analytical prediction techniques required to support advance configuration designs.

Description: The wall-matching methodology of Wilcox is modified to include a solid-wall, thermal-conduction model. This coupled fluid-thermal-structure model is derived assuming that the wall thermal-structure behavior is locally one-dimensional and that structural deformations, due to thermally induced stresses, are not significant. The one-dimensional coupled fluid-thermal-structure model is derived such that the wall temperature is removed as an independent boundary condition variable. The one-dimensional coupled fluid-thermal-structure model is also derived for the general case of an arbitrary mixture of thermally prefect gases and a wall of arbitrary thickness and conductivity by using a compressible, streamwise-pressure-gradient-corrected, wall-matching function and Fourier's law of heat conduction. The resulting model was implemented in the VULCAN CFD code as a new boundary condition type. VULCAN was then used to simulate a two-dimensional Mach 6 wind tunnel facility nozzle flow to demonstrate/validate the one-dimensional coupled fluid-thermal-structure model. The nozzle internal-wall surface temperature and heat transfer distributions computed using the one-dimensional coupled fluid-thermal-structure model are compared to wall temperature and heat transfer distributions from an iterative multi-dimensional analysis obtained by coupling the VULCAN CFD code and the MSC/NASTRAN-thermal code. The one-dimensional coupled fluid-thermal-structure model analysis is shown to be very robust and in excellent agreement with the multi-dimensional iteratively coupled analysis. It is also shown that the one-dimensional analysis can be used as an initial guess for the multi-dimensional iteratively coupled analysis.

Description: Results from Direct Numerical Simulations of temporal, supercritical mixing layers for two species systems are analyzed to elucidate species-specific turbulence aspects. The two species systems, O2/H2 and C7HG16/N2, have different thermodynamic characteristics; thus, although the simulations are performed at similar reduced pressure (ratio of the pressure to the critical pressure), the former system is dose to mixture ideality and has a relatively high solubility with respect to the latter, which exhibits strong departures from mixture ideality Due to the specified, smaller initial density stratification, the C7H16/N2 layers display higher growth and increased global molecular mixing as well as larger turbulence levels. However, smaller density gradients at the transitional state for the O2/H2 system indicate that on a local basis, the layer exhibits an enhanced mixing, this being attributed to the increased solubility and to mixture ideality. These thermodynamic features are shown to affect the irreversible entropy production (i.e. the dissipation), which is larger for the O2/H2 layer and is primarily concentrated in high density-gradient magnitude regions that are distortions of the initial density stratification boundary. In contrast, the regions of largest dissipation in the C7H16/N2 layer are located in high density-gradient magnitude regions resulting from the mixing of the two fluids.

Description: High temperature composite heat exchangers are an enabling technology for a number of aeropropulsion applications. They offer the potential for mass reductions of greater than fifty percent over traditional metallics designs and enable vehicle and engine designs. Since they offer the ability to operate at significantly higher operating temperatures, they facilitate operation at reduced coolant flows and make possible temporary uncooled operation in temperature regimes, such as experienced during vehicle reentry, where traditional heat exchangers require coolant flow. This reduction in coolant requirements can translate into enhanced range or system payload. A brief review of the approaches and challengers to exploiting this important technology are presented, along with a status of recent government-funded projects.

Description: Nonequilibrium phenomena in hypersonic flows are examined on the basis of theoretical models and selected experimental data, in an introduction intended for second-year graduate students of aerospace engineering. Chapters are devoted to the physical nature of gas atoms and molecules, transitions of internal states, the formulation of the master equation of aerothermodynamics, the conservation equations, chemical reactions in CFD, the behavior of air flows in nonequilibrium, experimental aspects of nonequilibrium flow, a review of experimental results, and gas-solid interaction. Diagrams, graphs, and tables of numerical data are provided.

Description: A modeling of the vortex-airfoil interaction is presented in which the finite-area of the real vortices is taken into consideration. Two vortex models are used. In the first, a disturbed piece of vorticity layer is simulated by four rows of discrete vortices of small strength. In the second, a number of discrete vortices is arranged within a circle. The first model may simulate a shear layer or a wake, while the second, a well-formed vortex. The method was applied to the calculation of the pressure induced on the surface of the airfoil by the interacting vortex. Both models give similar results. It was found that for large distances of the vortex from the surface of the airfoil, the consideration or not of the finite-area of the vortex is not a significant factor in determining the induced pressure field. However, when the distance of the vortex from the surface is reduced, its shape is distorted and the induced pressure pulses have lower amplitude than the ones induced by an equivalent point vortex. In the limit, where the vortex impinges on the leading edge of the airfoil, it is split into two and the time dependent pressure coefficient takes even negative values at some time intervals.

Description: The flows around highly sweptback wings and bodies of revolution at high angle of attack are described, and inviscid model approximations and mathematical formulation of the problem are given to steady and unsteady incompressible flows. A general presentation of the methods of solution is given, with emphasis on current computational techniques. Detailed descriptions of the nonlinear vortex-lattice and vortex-panel techniques are presented to show how the boundary conditions are enforced using iteration. Typical numerical results are compared with the available experimental data.

Description: A brief review is presented of various problems which are confronted in the development of an unsteady finite difference potential code. This review is conducted mainly in the context of what is done for a typical small disturbance and full potential methods. The issues discussed include choice of equation, linearization and conservation, differencing schemes, and algorithm development. A number of applications including unsteady three-dimensional rotor calculation, are demonstrated.

Description: The computational treatment of unsteady transonic flows is discussed, reviewing the historical development and current techniques. The fundamental physical principles are outlined; the governing equations are introduced; three-dimensional linearized and two-dimensional linear-perturbation theories in frequency domain are described in detail; and consideration is given to frequency-domain FEMs and time-domain finite-difference and integral-equation methods. Extensive graphs and diagrams are included.

Description: This lecture is introductory to the subject of unsteady subsonic and supersonic flows. The primary objective is to present fundamental concepts in order to promote an understanding of the relations between the basic physical problems and their mathematical formulation as well as to establish a common foundation for the more detailed presentations of subsequent lectures in this session. Linearized (small-perturbation) potential flow is emphasized, although needs beyond that limit are indicated. The basic equations, concepts, and procedures common to all the methods are reviewed first, followed by the development, discussion, and status of methods for creating two-dimensional incompressible flow, strip theory, subsonic lifting-surface theory, subsonic/supersonic surface-panel methods, and supersonic lifting-surface theory.

Description: The prediction of missile aerodynamic characteristics is presently undertaken through the application of supersonic paneling methods and nonlinear corrections to the prediction of missile aerodynamic characteristics. Attention is given to supersonic panel methods and line-singularity methods for the modeling of axisymmetric bodies, in combination with corrections for nonlinear flow phenomena, which are applied to complete missile, inlets, and wing-body combinations. The LRCDM2 computer program is used as an example of the methods presented.

Description: The discrete vortex cloud approach models a missile airframe's vortex wake by means that are capable of treating a variety of configurations over a range of flow conditions. Attention is given to the sheets of vorticity formed on the lee side of a missile at moderate angles-of-attack. While three-dimensional attached flow models are used to represent the missile body, two-dimensional, incompressible, separated flow models are used to represent the separated vortex wake. The predicted pressure distribution of the body under the influence of the freestream and the separation vortex wake are used to calculate aerodynamic loads on the body. The separation vortex wake is represented by clouds of discrete vortices in cross flow planes normal to the body axis.

Description: Applications of laser velocimetry to the measurement of turbulent flow properties of strong transonic viscous-inviscid interactions are reviewed. The data resulting from these studies are then discussed in relation to their importance in the development of improved viscous-flow calculation methods. Also considered are the current limitations of laser velocimetry, the need for further improvements in the method, and potential future applications.

Description: The results of Reynolds-averaged time-dependent inviscid and turbulent compressible Navier-Stokes computations using the implicit finite-difference approach of Steger (1978), modified by incorporating a pressure boundary condition, (PBC) to account for wall interference are compared with experimental data on a NACA 64A010 airfoil (Johnson and Bachalo, 1980) in graphs and briefly characterized. The computational approach is the same as that used by King and Johnson (1980), but a 137 x 50 mesh is used instead of a 97 x 35 mesh, and special care is taken in resolving the nose, shock, and trailing-edge regions. Imposition of PBC is shown to improve significantly the accuracy of the computations for the flowfield on the upper surface of the airfoil, shifting the shock forward to its experimentally measured position in the case of turbulent flow. The failure of the method, even with PBC, to match the experimental shock location in the case of a flow with a separation bubble is attributed to inadequacies in the algebraic turbulence model employed (Baldwin and Lomax, 1978).

Description: Current progress in the computational analysis of rotary-wing flowfields is surveyed, and some typical results are presented in graphs. Topics examined include potential theory, rotating coordinate systems, lifting-surface theory (moving singularity, fixed wing, and rotary wing), panel methods (surface singularity representations, integral equations, and compressible flows), transonic theory (the small-disturbance equation), wake analysis (hovering rotor-wake models and transonic blade-vortex interaction), limitations on computational aerodynamics, and viscous-flow methods (dynamic-stall theories and lifting-line theory). It is suggested that the present algorithms and advanced computers make it possible to begin working toward the ultimate goal of turbulent Navier-Stokes calculations for an entire rotorcraft.

Description: Unsteady aerodynamic and aeroelastic stability calculations based upon transonic small disturbance (TSD) potential theory are presented. Results from the two-dimensional XTRAN2L code and the three-dimensional XTRAN3S code are compared with experiment to demonstrate the ability of TSD codes to treat transonic effects. The necessity of nonisentropic corrections to transonic potential theory is demonstrated. Dynamic computational effects resulting from the choice of grid and boundary conditions are illustrated. Unsteady airloads for a number of parameter variations including airfoil shape and thickness, Mach number, frequency, and amplitude are given. Finally, samples of transonic aeroelastic calculations are given. A key observation is the extent to which unsteady transonic airloads calculated by inviscid potential theory may be treated in a locally linear manner.

Description: An algebraic procedure for the generation of boundary-fitted grids about wing-fuselage configurations is presented. A wing-fuselage configuration is specified by cross sections and mathematically represented by Coons' patches. A configuration is divided into sections so that several grid blocks that either adjoin each other or partially overlap each other can be generated. Each grid has six exterior surfaces that map into a computational cube. Grids are first determined on the six boundary surfaces and then in the interior. Grid curves that are on the surface of the configuration are derived from the intersection of planes with the Coons' patch definition. Single-valued functions relating approximate arc lengths along the grid curves to a computational coordinate define the distribution of grid points. The two-boundary technique and transfinite interpolation are used to determine the boundary surface grids that are not on the configuration, and transfinite interpolation with linear blending functions is used to determine the interior grid.

Description: An overview of helicopter aerodynamics technology is presented with emphasis on rotor wake and airloads methodology developed at the United Technologies Research Center (UTRC). The evolution over the past twenty years of various levels of computerized wake geometry models at UTRC, such as undistorted wake, prescribed empirical wake, predicted distorted wake, and generalized wake models for the hover and forward flight regimes, is reviewed. The requirement for accurate wake modeling for flow field and airload prediction is demonstrated by comparisons of theoretical and experimental results. These results include blade pressure distributions predicted from a recently developed procedure for including the rotor wake influence in a full potential flow analysis. Predictions of the interactional aerodynamics of various helicopter components (rotor, fuselage, and tail) are also presented. It is concluded that, with advanced computers and the rapidly progressing computational aerodynamics technology, significant progress toward reliable prediction of helicopter airloads is forseeable in the near future.

Description: In the present treatment of the calculation of forces on a wing that is suddenly brought into motion at a constant speed, attention is given to the unsteady potential's contribution to the force balance. Total bound vorticity is produced at the initial impulse. The results obtained are independent of wing aspect ratio; as time increases, this effect on the drag force becomes smaller as the vortex emanating from the trailing edge is left behind. The second contributor to induced drag is the spanwise vorticity shedding that results from the spanwise load distribution of three-dimensional wings. This contribution grows with time as the length of the wake grows.

Description: Laminar heating distributions have been measured on a 1.9 percent scale model of a generic aeroassisted vehicle taking the shape of a spherically blunted, 13-deg/7-deg biconic whose forecone section is bent upward (by 7 deg) to furnish self-trim capability at a 20-deg angle-of-attack. The results thus obtained were compared with data gathered for a straight biconic. While no Reynolds number effect on heating was noted on the windward side of the forecone, the opposite was true of the leeward side, where a Reynolds number increase caused circumferential flow separation at lower angles of attack. Generally, windward heating was predicted to within 10 percent with a computer code solving the steady, three-dimensional parabolized Navier-Stokes equations.

Description: Interferometry methods were applied to the investigation of steady and unsteady flows in large scale transonic wind tunnels. Holographic interferometry was demonstrated to provide reliable flow visualization and quantitative results for a number of two-dimensional flows. These conclusions were based on extensive comparisons with results obtained by other means. Data obtained on a NACA 64A010 airfoil with an oscillating flap installed in the Ames 11-foot transonic tunnel are presented. Interferograms were recorded at a free stream Mach number of 0.8, flap frequency of 30 Hertz and chord Reynolds numbers of 6.6 x 10 to the 6th and 12.3 x 10 to the 6th. The interferometric results were reduced to dynamic surface pressures, Mach contours and wake flow profiles. A new interferometry method that is capable of providing real-time interferometry data is also discussed.

Description: Flow and sound field data are presented for a 2.54 cm diameter air jet at a Mach number of 0.50 and a Reynolds number of 300,000. Distributions of mean velocity, turbulence intensities, Reynolds stress, spectral components of turbulence as well as of the near field pressure, together with essential characteristics of the far field sound are reported. This detailed set of data for one particular flow, erstwhile unavailable in the literature, is expected to help promoote and calibrate subsonic jet noise theories. 'Source locations' in terms of the turbulence maxima, coupling between the entrainment dynamics and the near pressure field, the sound radiation paths, and the balance in mass, momentum and sound energy fluxes are discussed. The results suggest that the large scale coherent structures of the jet govern the 'source locations' by controlling the turbulence and also strongly influence the near field pressure fluctuations.

Description: Computational methods solving the thin shear layer formulation of the compressible, Reynolds-averaged Navier-Stokes equations are presently used to investigate the strongly interactive flow field about aircraft afterbodies. Solutions for a variety of axisymmetric afterbody and nozzle geometries are solved by means of a time-dependent implicit numerical algorithm for both subsonic and supersonic external flows, and the results obtained are compared with experimental data. A novel adaptive-grid technique is used to resolve flow regimes having large gradients, as well as to improve the accuracy and efficiency of the computational scheme.

Description: The aerodynamic performance and controllability of advanced, highly maneuverable supersonic aircraft can be enhanced by means of 'vortex management', which refers to the purposeful manipulation and reordering of stable and concentrated vortical structures due to flow separations from highly swept leading edges and slender forebodies at moderate-to-high angles-of-attack. Attention is presently given to a variety of results obtained in the course of experiments on generic research models at NASA Langley, clarifying their underlying aerodynamics and evaluating their performance-improvement potential. The vortex-management concepts discussed encompass aerodynamic compartmentation of highly swept leading edges, vortex lift augmentation and modulation, and forebody vortex manipulation.

Description: The CRAY multitasking system was developed in order to utilize all four processors and sharply reduce the wall clock run time. This paper describes the techniques used to modify the computational fluid dynamics code ARC3D for this run and analyzes the achieved speedup. The ARC3D code solves either the Euler or thin-layer N-S equations using an implicit approximate factorization scheme. Results indicate that multitask processing can be used to achieve wall clock speedup factors of over three times, depending on the nature of the program code being used. Multitasking appears to be particularly advantageous for large-memory problems running on multiple CPU computers.

Description: (Previously cited in issue 20, p. 2915, Accession no. A86-42687)

Description: The NCOREL full-potential method with an entropy correction is presently applied to supersonic missile flowfield problems. After defining the salient characteristics of the method, a combination of linear theory with NCOREL and experimental data is used to isolate the nonlinear features of the supersonic flow so that the influence of geometry and flow conditions on the development of such flow nonlinearities can be appreciated. Comparisons of experimental longitudinal force and moment data with NCOREL and various linear theory predictions are presented for several generic missile airframe configurations of circular and elliptic cross section. The NCOREL code solves the nonconservative full potential equation in a spherical coordinate system; exact boundary conditions are defined on the missile surface.

Description: A comprehensive evaluation is made of experimental data compiled to date for the flowfields and aerodynamic forces that occur at high angles of attack for low aspect ratio wings with delta, rectangular, clipped delta, and strake/wing planform geometries. Attention is given to wing leading edge-generated vortex breakdown, aspect ratio and compressibility effects, and strake vortex effects on main wing areas. Although the nonlinear effects created by a wing-body combination significantly alter wing-alone aerodynamics, the wing-alone data presented are vital to the development of prediction methodologies for large angle of attack aerodynamics.

Description: A solution procedure is presented for the lifting transonic flow past modern rotor configurations in forward flight. In this procedure, the three-dimensional, unsteady Euler equations are solved in strong conservation form on a body-fitted moving coordinate system. A hybrid procedure of second order spatial accuracy and first order temporal accuracy is used to integrate the governing equations. In lifting flows, the effect of the elements of wake not captured by the computational procedure, and other aeroelastic effects are accounted for as local angle of attack corrections. Detailed comparisons with experimental data are presented for a 1/7 scale model of the Cobra OLS rotor, and for a three-bladed rotor tested in France. Some preliminary results are also presented for a three-dimensional blade vortex interaction problem.

Description: Rotary-wing computational fluid dynamics is reaching a point where many three-dimensional, unsteady, finite-difference codes are becoming available. This paper gives a brief review of five such codes, which treat the small disturbance, conservative and nonconservative full-potential, and Euler flow models. A discussion of the methods of applying these codes to the rotor environment (including wake and trim considerations) is followed by a comparison with various available data. These data include tests of advancing lifting and nonlifting, and hovering model rotors with significant supercritical flow regions. The codes are also compared for computational efficiency.

Description: The propagation characteristics of several helicopter airfoil profiles have been investigated using the transonic small disturbance equation. A test case was performed to generate a moving shock that propagated off the airfoil. Various grids were then examined to determine their ability to accurately capture these propagating shock waves. Finally, the case of airfoil-vortex interactions was thoroughly studied over a wide range of Mach numbers and airfoil shapes with particular emphasis on the transonic regime; this results in a highly conplicated fluctuation of lift, drag, and pitching moment. The calculated acoustic intensity levels, along with the details of the computational flow field, provide new insights into the understanding of transonic airfoil-vortex interactions.

Description: An experimental and theoretical investigation of rotor/wing aerodynamic interactions in hover is described. The experimental investigation consisted of both a large-scale and small-scale test. A 0.658-scale, V-22 rotor and wing was used in the large-scale test. Wind download, wing surface pressure, rotor performance, and rotor downwash data from the large-scale test are presented. A small-scale experiment was conducted to determine how changes in the rotor/wing geometry affected the aerodynamic interactions. These geometry variations included the distance between the rotor and wing, wing incidence angle, and configurations both with the rotor axis at the tip of the wing (tilt rotor configuration) and with the rotor axis at the center of the wing (compound helicopter oonfiguration). A wing with boundary-layer control was also tested to evaluate the effect of leading and trailing edge upper surface blowing on the wing download. A computationally efficient, semi-empirical theory was developed to predict the download on the wing. Finally, correlations between the theoretical predictions and test data are presented.

Description: It is proposed that the study of Rusak et al. (1985), which reports numerical modeling sensitivities on longitudinal force/moment properties for a vortex-lattice method incorporating free vortex filaments to represent the leading-edge vortex separation, employs a formula that is strongly affected by the particular points of analysis chosen. This results in a narrowly applicable curve fit, where numerical sensitivities of the theory are inappropriately traded off against physical effects that are not modeled in that theory. Attention is also given to questionable drag estimate computations.

Description: Steady, high speed, compressible separated flows modeled through numerical simulations resulting from solutions of the mass-averaged Navier-Stokes equations are reviewed. Emphasis is placed on benchmark flows that represent simplified (but realistic) aerodynamic phenomena. These include impinging shock waves, compression corners, glancing shock waves, trailing edge regions, and supersonic high angle of attack flows. A critical assessment of modeling capabilities is provided by comparing the numerical simulations with experiment. The importance of combining experiment, numerical algorithm, grid, and turbulence model to effectively develop this potentially powerful simulation technique is stressed.

Description: The analysis and the incorporation into a multigrid scheme of several vectorizable algorithms are discussed. von Neumann analyses of vertical-line, horizontal-line, and alternating-direction ZEBRA algorithms were performed; and the results were used to predict their multigrid damping rates. The algorithms were then successfully implemented in a transonic conservative full-potential computer program. The convergence acceleration effect of multiple grids is shown, and the convergence rates of the vectorizable algorithms are compared with those of standard successive-line overrelaxation (SLOR) algorithms.

Description: Numerical procedures for solving the thin-shear-layer Navier-Stokes equations and for the interaction of solutions to inviscid and boundary-layer equations are described and evaluated. To allow appraisal of the numerical and fluid dynamic abilities of the two schemes, they have been applied to one airfoil as a function of angle of attack at two slightly different Reynolds numbers. The NACA 0012 airfoil has been chosen because it allows comparison with measured lift, drag, and moment and with surface-pressure distributions. Calculations have been performed with algebraic eddy-viscosity formulations, and they include consideration of transition. The results are presented in a form that allows easy appraisal of the accuracy of both procedures and of the relative costs. The interactive procedure is computationally efficient but restrictive relative to the thin-layer Navier-Stokes procedure. The latter procedure does a better job of predicting drag than does the former. In both procedures, the location of transition is crucial for accurate or detailed computations, particularly at high angles of attack. When the upstream influence of pressure field through the shear layer is important, the thin-layer Navier-Stokes procedure has an edge over the interactive procedure.

Description: Recent progress in the development of finite element methodology for the prediction of aerothermal loads is described. Three dimensional, inviscid computations are presented, but emphasis is placed on development of an approach extendable to three dimensional viscous flows. Progress in key research areas is described. Initial 30 results from the computational procedure are described.

Description: The numerical method developed by Schiff and Sturek (1980) on the basis of the thin-layer parabolized Navier-Stokes equations of Schiff and Steger (1980) is extended to the case of turbulent supersonic flows on pointed bodies at high angles of attack. The governing equations, the numerical scheme, and modifications to the algebraic eddy-viscosity turbulence model are described; and results for three cones and one ogive-cylinder body (obtained using grids of 50 nonuniformly spaced points in the radial direction between the body and the outer boundary) are presented graphically and compared with published experimental data. The grids employed are found to provide sufficient spatial resolution of the leeward-side vortices; when combined with the modified turbulence model, they are shown to permit accurate treatment of flows with large regions of crossflow separation.

Description: The present use of a parabolized Navier-Stokes solver to accurately simulate the flowfield in a supersonic inlet yields good agreement between numerical analysis and experiment for a Mach 7.4 inlet under cruise conditions, with an internal compression ratio of 8. The significance of real gas effects on the performance calculation of a hypersonic inlet is demonstrated, with small changes in the ratio of specific heats resulting in a substantial change in the calculated pitot pressure ratio.

Description: Through a series of flights in artificial clouds, ice accretions on the main rotor of a UH-1H helicopter were documented in detail upon landing by silicone-rubber molds for both hover and level flights. Full scale reproductions of typical accretions in hover were fabricated by means of epoxy castings and used for a wind-tunnel test program. Surface static pressure distributions were recorded and used to evaluate lift and pitching moment increments while drag was determined by wake surveys. For comparison, accreted ice shapes are presented for two level flight cases as well as preliminary analytical predictions.

Description: The mechanism that locates a shock wave in a transonic flow in one and two dimensions is examined. It is found that in one dimension the shock is located by specifying the downstream pressure whereas in two dimensions the shock is located by the application of an entropy condition at the sonic line.