SEARCH · Search NASA

Results for “fluid structure interaction”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 649 records · Page 36

All sources NASA NTRS DOE OSTI

NASA NTRS2002

Prediction of Drag Reduction in Supersonic and Hypersonic Flows with Counterflow Jets

Computational fluid dynamics solutions of the flowfield of a truncated cone-cylinder with and without counterflow jets have been obtained for the short penetration mode (SPM) and long penetration mode (LPM) of the freestream-counterflow jet interaction flowfield. For the case without the counterflow jet, the comparison of the normalized surface pressures showed very good agreement with experimental data. For the case with the SPM jet, the predicted surface pressures did not compare as well with the experimental data upstream of the expansion corner, while aft of the expansion corner, the comparison of the solution and the data is seen to give much better agreement. The difference in the prediction and the data could be due to the transient character of the jet penetration modes, possible effects of the plasma physics that are not accounted for here, or even the less likely effect of flow turbulence, etc. For the LPM jet computations, one-dimensional isentropic relations were used to derived the jet exit conditions in order to obtain the LPM solutions. The solution for the jet exit Mach number of 3 shows a jet penetration several times longer than that of the SPM, and therefore much weaker bow shock, with an attendant reduction in wave drag. The LPM jet is, in essence, seen to be a "pencil" of fluid, with much higher dynamic pressure, embedded in the oncoming supersonic or hypersonic freestream. The methodology for determining the conditions for the LPM jet could enable a practical approach for the design and application of counterflow LPM jets for the reduction of wave drag and heat flux, thus significantly enhancing the aerodynamic characteristics and aerothermal performance of supersonic and hypersonic vehicles. The solutions show that the qualitative flow structure is very well captured. The obtained results, therefore, suggest that counterflowing jets are viable candidate technology concepts that can be employed to give significant reductions in wave drag, heat flux, and other attendant aerodynamic benefits.

Daso, Endwell O.↗

NASA NTRS2004

Experimental, Theoretical, and Computational Investigation of Separated Nozzle Flows

A detailed experimental, theoretical, and computational study of separated nozzle flows has been conducted. Experimental testing was performed at the NASA Langley 16-Foot Transonic Tunnel Complex. As part of a comprehensive static performance investigation, force, moment, and pressure measurements were made and schlieren flow visualization was obtained for a sub-scale, non-axisymmetric, two-dimensional, convergent- divergent nozzle. In addition, two-dimensional numerical simulations were run using the computational fluid dynamics code PAB3D with two-equation turbulence closure and algebraic Reynolds stress modeling. For reference, experimental and computational results were compared with theoretical predictions based on one-dimensional gas dynamics and an approximate integral momentum boundary layer method. Experimental results from this study indicate that off-design overexpanded nozzle flow was dominated by shock induced boundary layer separation, which was divided into two distinct flow regimes; three- dimensional separation with partial reattachment, and fully detached two-dimensional separation. The test nozzle was observed to go through a marked transition in passing from one regime to the other. In all cases, separation provided a significant increase in static thrust efficiency compared to the ideal prediction. Results indicate that with controlled separation, the entire overexpanded range of nozzle performance would be within 10% of the peak thrust efficiency. By offering savings in weight and complexity over a conventional mechanical exhaust system, this may allow a fixed geometry nozzle to cover an entire flight envelope. The computational simulation was in excellent agreement with experimental data over most of the test range, and did a good job of modeling internal flow and thrust performance. An exception occurred at low nozzle pressure ratios, where the two-dimensional computational model was inconsistent with the three-dimensional separation observed in the experiment. In general, the computation captured the physics of the shock boundary layer interaction and shock induced boundary layer separation in the nozzle, though there were some differences in shock structure compared to experiment. Though minor, these differences could be important for studies involving flow control or thrust vectoring of separated nozzles. Combined with other observations, this indicates that more detailed, three-dimensional computational modeling needs to be conducted to more realistically simulate shock-separated nozzle flows.

Hunter, Craig A.↗

NASA NTRS2009

SWIFT Code Assessment for Two Similar Transonic Compressors

One goal of the NASA Fundamental Aeronautics Program is the assessment of computational fluid dynamic (CFD) codes used for the design and analysis of many aerospace systems. This paper describes the assessment of the SWIFT turbomachinery analysis code for two similar transonic compressors, NASA rotor 37 and stage 35. The two rotors have identical blade profiles on the front, transonic half of the blade but rotor 37 has more camber aft of the shock. Thus the two rotors have the same shock structure and choking flow but rotor 37 produces a higher pressure ratio. The two compressors and experimental data are described here briefly. Rotor 37 was also used for test cases organized by ASME, IGTI, and AGARD in 1994-1998. Most of the participating codes over predicted pressure and temperature ratios, and failed to predict certain features of the downstream flowfield. Since then the AUSM+ upwind scheme and the k- turbulence model have been added to SWIFT. In this work the new capabilities were assessed for the two compressors. Comparisons were made with overall performance maps and spanwise profiles of several aerodynamic parameters. The results for rotor 37 were in much better agreement with the experimental data than the original blind test case results although there were still some discrepancies. The results for stage 35 were in very good agreement with the data. The results for rotor 37 were very sensitive to turbulence model parameters but the results for stage 35 were not. Comparison of the rotor solutions showed that the main difference between the two rotors was not blade camber as expected, but shock/boundary layer interaction on the casing.

Chima, Rodrick V.↗

NASA NTRS2000

Inertial Effects in Suspension Dynamics

The present work analyses the dynamics of a suspension of heavy particles in shear flow. The magnitude of the particle inertia is given by the Stokes number St = m(gamma/6(pi)a, which is the ratio of the viscous relaxation time of a particle tau(sub p) = m=6pi(eta)a to the flow time gamma(sup -1). Here, m is the mass of the particle, a is its size, eta is the viscosity of the suspending fluid and gamma is the shear rate. The ratio of the Stokes number to the Reynolds number, Re = (rho)f(gamma)a(exp 2)/eta, is the density ratio rho(sub p)/rho(sub f). Of interest is to understand the separate roles of particle (St) and fluid (Re) inertia in the dynamics of suspensions. In this study we focus on heavy particles, rho(sub p)/rho(sub f) much greater than 1, for which the Stokes number is finite, but the Reynolds number is sufficiently small for inertial forces in the fluid to be neglected; thus, the fluid motion is governed by the Stokes equations. On the other hand, the probability density governing the statistics of the suspended particles satisfies a Fokker-Planck equation that accounts for both configuration and momentum coordinates, the latter being essential for finite St. The solution of the Fokker-Planck equation is obtained to O(St) via a Chapman-Enskog type-procedure, and the conditional velocity distribution so obtained is used to derive a configuration-space Smoluchowski equation with inertial corrections. The inertial effects are responsible for asymmetry in the relative trajectories of two spheres in shear flow, in contrast to the well known symmetric structure in the absence of inertia. Finite St open trajectories in the plane of shear suffer a downward lateral displacement resulting from the inability of a particle of finite mass to follow the curvature of the zero-Stokes-number pathlines. In addition to the induced asymmetry, the O(St) inertial perturbation dramatically alters the nature of the near-field trajectories. The stable closed orbits (for St = 0) in the plane of shear now spiral in, approaching particle-particle contact in the limit. All trajectories starting from an initial offset of O(St(sup 1/2) or less (which remain open for St = 0) also spiral in. The asymmetry of the trajectories leads to a non-Newtonian rheology and diffusive behavior. The latter because a given particle (moving along a finite St open trajectory) suffers a net displacement in the transverse direction after a single interaction. A sequence of such uncorrelated displacements leads to the particle executing a random walk. The inertial diffusivity tensor is anisotropic on account of differing strengths of interaction in the gradient and vorticity directions. Since the entire region (constituting an in finite area) of closed orbits in the plane of shear spirals onto contact for #finite St, the latter represents a singular surface for the pair-distribution function. The exact form of the pair-distribution function at contact is still, however, indeterminate in the absence of non-hydrodynamic effects. It should also be noted that finite St non-rectilinear flows do not support a spatially uniform number density owing to the cross-streamline inertial migration of particles.

J. F. Brady↗

NASA NTRS1999

The Effects of Buoyancy and Dilution on the Structure and Lift-off of Coflow Laminar Diffusion Flames

The ability to predict the coupled effects of complex transport phenomena with detailed chemical kinetics in diffusion flames is critical in the modeling of turbulent reacting flows and in understanding the processes by which soot formation and radiative transfer take place. In addition, an understanding of the factors that affect flame extinction in diffusion flames is critical in the suppression of fires and in improving engine efficiency. The goal of our characterizations of coflow laminar diffusion flames is to bring to microgravity the multidimensional diagnostic tools available in normal gravity, and in so doing provide a broader understanding of the successes and limitations of current combustion models. This will lead to a more detailed understanding of the interaction of convection, diffusion and chemistry in both buoyant and nonbuoyant environments. As a sensitive marker of changes in the flame shape, the number densities of excited-state CH (A(exp 2)delta, denoted CH*), and excited-state OH (A(exp 2)Sigma, denoted OH*) are measured in mu-g and normal gravity. Two-dimensional CH* and OH* number densities are deconvoluted from line-of-sight chemiluminescence measurements made on the NASA KC-135 reduced-gravity aircraft. Measured signal levels are calibrated, post-flight, with Rayleigh scattering. Although CH* and OH* kinetics are not well understood, the CH*, OH*, and ground-state CH distributions are spatially coincident in the flame anchoring region. Therefore, the ground-state CH distribution, which is easily computed, and the readily measured CH*/OH* distributions can be used to provide a consistent and convenient way of measuring lift-off height and flame shape in the diffusion flame under investigation. Given that the fuel composition affects flame chemistry and that buoyancy influences the velocity profile of the flow, we have the opportunity to computationally and experimentally study the roles of fluids and chemistry. In performing this microgravity study, improvements to the computational model have been made and new calculations performed for a range of gravity and flow conditions. Furthermore, modifications to the experimental approach were required as a consequence of the constraints imposed by existing microgravity facilities. Results from the computations and experiments are presented in the following sections.

Walsh, Kevin T.↗

NASA NTRS2010

Analysis of Nozzle Jet Plume Effects on Sonic Boom Signature

An axisymmetric full Navier-Stokes computational fluid dynamics (CFD) study was conducted to examine nozzle exhaust jet plume effects on the sonic boom signature of a supersonic aircraft. A simplified axisymmetric nozzle geometry, representative of the nozzle on the NASA Dryden NF-15B Lift and Nozzle Change Effects on Tail Shock (LaNCETS) research airplane, was considered. The highly underexpanded nozzle flow is found to provide significantly more reduction in the tail shock strength in the sonic boom N-wave pressure signature than perfectly expanded and overexpanded nozzle flows. A tail shock train in the sonic boom signature, similar to what was observed in the LaNCETS flight data, is observed for the highly underexpanded nozzle flow. The CFD results provide a detailed description of the nozzle flow physics involved in the LaNCETS nozzle at different nozzle expansion conditions and help in interpreting LaNCETS flight data as well as in the eventual CFD analysis of a full LaNCETS aircraft. The current study also provided important information on proper modeling of the LaNCETS aircraft nozzle. The primary objective of the current CFD research effort was to support the LaNCETS flight research data analysis effort by studying the detailed nozzle exhaust jet plume s imperfect expansion effects on the sonic boom signature of a supersonic aircraft. Figure 1 illustrates the primary flow physics present in the interaction between the exhaust jet plume shock and the sonic boom coming off of an axisymmetric body in supersonic flight. The steeper tail shock from highly expanded jet plume reduces the dip of the sonic boom N-wave signature. A structured finite-volume compressible full Navier-Stokes CFD code was used in the current study. This approach is not limited by the simplifying assumptions inherent in previous sonic boom analysis efforts. Also, this study was the first known jet plume sonic boom CFD study in which the full viscous nozzle flow field was modeled, without coupling to a sonic boom propagation analysis code, from the stagnation chamber of the nozzle to the far field external flow, taking into account all nonisentropic effects in the shocks, boundary layers, and free shear layers, and their interactions at distances up to 30 times the nozzle exit diameter from the jet centerline. A CFD solution is shown in Figure 2. The flow field is very complicated and multi-dimensional, with shock-shock and shockplume interactions. At the time of this reporting, a full three-dimensional CFD study was being conducted to evaluate the effects of nozzle vectoring on the aircraft tail shock strength.

Bui, Trong↗

NASA NTRS2009

Design of Experiments for Both Experimental and Analytical Study of Exhaust Plume Effects on Sonic Boom

Computational fluid dynamics (CFD) analysis has been performed to study the plume effects on sonic boom signature for isolated nozzle configurations. The objectives of these analyses were to provide comparison to past work using modern CFD analysis tools, to investigate the differences of high aspect ratio nozzles to circular (axisymmetric) nozzles, and to report the effects of under expanded nozzle operation on boom signature. CFD analysis was used to address the plume effects on sonic boom signature from a baseline exhaust nozzle. Nearfield pressure signatures were collected for nozzle pressure ratios (NPRs) between 6 and 10. A computer code was used to extrapolate these signatures to a ground-observed sonic boom N-wave. Trends show that there is a reduction in sonic boom N-wave signature as NPR is increased from 6 to 10. As low boom designs are developed and improved, there will be a need for understanding the interaction between the aircraft boat tail shocks and the exhaust nozzle plume. These CFD analyses will provide a baseline study for future analysis efforts. For further study, a design of experiments has been conducted to develop a hybrid method where both CFD and small scale wind tunnel testing will validate the observed trends. The CFD and testing will be used to screen a number of factors which are important to low boom propulsion integration, including boat tail angle, nozzle geometry, and the effect of spacing and stagger on nozzle pairs. To design the wind tunnel experiment, CFD was instrumental in developing a model which would provide adequate space to observe the nozzle and boat tail shock structure without interference from the wind tunnel walls.

Castner, Raymond S.↗

NASA NTRS1995

Large HAWT wake measurement and analysis

From the theoretical fluid dynamics point of view, the wake region of a large horizontal-axis wind turbine has been defined and described, and numerical models of wake behavior have been developed. Wind tunnel studies of single turbine wakes and turbine array wakes have been used to verify the theory and further refine the numerical models. However, the effects of scaling, rotor solidity, and topography on wake behavior are questions that remain unanswered. In the wind tunnel studies, turbines were represented by anything from scaled models to tea strainers or wire mesh disks whose solidity was equivalent to that of a typical wind turbine. The scale factor compensation for the difference in Reynolds number between the scale model and an actual turbine is complex, and not typically accounted for. Though it is wise to study the simpler case of wakes in flat topography, which can be easily duplicated in the wind tunnel, current indications are that wind turbine farm development is actually occurring in somewhat more complex terrain. Empirical wake studies using large horizontal-axis wind turbines have not been thoroughly composited, and, therefore, the results have not been applied to the well-developed theory of wake structure. The measurement programs have made use of both in situ sensor systems, such as instrumented towers, and remote sensors, such as kites and tethered, balloonborne anemometers. We present a concise overview of the work that has been performed, including our own, which is based on the philosophy that the MOD-2 turbines are probably their own best detector of both the momentum deficit and the induced turbulence effect downwind. Only the momentum deficit aspects of the wake/machine interactions have been addressed. Both turbine power output deficits and wind energy deficits as measured by the onsite meteorological towers have been analyzed from a composite data set. The analysis has also evidenced certain topographic influences on the operation of spatially diverse wind turbines.

Miller, A. H.↗

NASA NTRS2000

Physics of Hard Spheres Experiment: Significant and Quantitative Findings Made

Direct examination of atomic interactions is difficult. One powerful approach to visualizing atomic interactions is to study near-index-matched colloidal dispersions of microscopic plastic spheres, which can be probed by visible light. Such spheres interact through hydrodynamic and Brownian forces, but they feel no direct force before an infinite repulsion at contact. Through the microgravity flight of the Physics of Hard Spheres Experiment (PHaSE), researchers have sought a more complete understanding of the entropically driven disorder-order transition in hard-sphere colloidal dispersions. The experiment was conceived by Professors Paul M. Chaikin and William B. Russel of Princeton University. Microgravity was required because, on Earth, index-matched colloidal dispersions often cannot be density matched, resulting in significant settling over the crystallization period. This settling makes them a poor model of the equilibrium atomic system, where the effect of gravity is truly negligible. For this purpose, a customized light-scattering instrument was designed, built, and flown by the NASA Glenn Research Center at Lewis Field on the space shuttle (shuttle missions STS 83 and STS 94). This instrument performed both static and dynamic light scattering, with sample oscillation for determining rheological properties. Scattered light from a 532- nm laser was recorded either by a 10-bit charge-coupled discharge (CCD) camera from a concentric screen covering angles of 0 to 60 or by sensitive avalanche photodiode detectors, which convert the photons into binary data from which two correlators compute autocorrelation functions. The sample cell was driven by a direct-current servomotor to allow sinusoidal oscillation for the measurement of rheological properties. Significant microgravity research findings include the observation of beautiful dendritic crystals, the crystallization of a "glassy phase" sample in microgravity that did not crystallize for over 1 year in 1g (Earth's gravity), and the emergence of face-centered-cubic (FCC) crystals late in the coarsening process (as small crystallites lost particles to the slow ripening of large crystallites). Significant quantitative findings from the microgravity experiments have been developed describing complex interactions among crystallites during the growth process, as concentration fields overlap in the surrounding disordered phase. Time-resolved Bragg scattering under microgravity captures one effect of these interactions quite conclusively for the sample at a volume fraction of 0.528. From the earliest time until the sample is almost fully crystalline, the size and overall crystallinity grow monotonically, but the number of crystallites per unit volume (number density) falls. Apparently nucleation is slower than the loss of crystallites because of the transfer of particles from small to large crystals. Thus, coarsening occurs simultaneously with growth, rather than following the completion of nucleation and growth as is generally assumed. In the same sample, an interesting signature appears in the apparent number density of crystallites and the volume fraction within the crystallites shortly before full crystallinity is reached. A brief upturn in both indicates the creation of more domains of the size of the average crystallite simultaneous with the compression of the crystallites. Only the emergence of dendritic arms offers a reasonable explanation. The arms would be "seen" by the light scattering as separate domains whose smaller radii of curvature would compress the interior phase. In fiscal year 1999, numerous papers, a doctoral dissertation, and the PHaSE final report were produced. Although this flight project has been completed, plans are in place for a follow-on colloid experiment by Chaikin and Russel that employs a light microscope within Glenn's Fluids and Combustion Facility on the International Space Station. PHaSE is providing us with a deeper understanding of the nure of phase transitions. The knowledge derived has added to the understanding of condensed matter. In addition, the burgeoning study of the dynamics of colloidal self-assembly may lead to the development of a range of photonic materials that control the desirable properties of light. Thus, applications of ordered colloidal structures include not only ultrastructure ceramics, but also photonic crystals and photothermal nanosecond light-switching devices. Industries dealing with semiconductors, electro-optics, ceramics, and composites stand to benefit from such advancements.

Doherty, Michael P.↗

DOE OSTI2024

Thermodynamic and Kinetic Modulation of Methylammonium Lead Bromide Crystallization Revealed by In Situ Monitoring

Hybrid organic–inorganic perovskite (HOIP) crystals are promising optoelectronic materials, but little is known about either the thermodynamic and kinetic controls on crystal growth or the underlying growth mechanism(s). Herein, we use fluid cell atomic force microscopy (AFM) and solution nuclear magnetic resonance (NMR) spectroscopy to investigate growth of the model HOIP crystal CH 3 NH 3 PbBr 3 (MAPbBr 3 ) and to determine how formic acid (HCOOH) modulates the thermodynamics and kinetics of growth. The results show that growth of MAPbBr 3 in dimethylformamide (DMF) proceeds through the classical pathway by the spreading of atomic crystal steps generated at screw dislocations on the {100} surface. Temperature dependent step velocity measurements demonstrate that with increasing concentration, HCOOH decreases both the solubility of MAPbBr 3 and the kinetic coefficient (b) of step movement. 1 H-NMR measurements indicate that HCOOH increases the lifetime of the methylammonium (MA + ) ions and promotes the association of MAPbBr 3 , thus tuning the solubility of the perovskite. HCOOH also alters the molecular tumbling motion and bulk diffusion of the MA + ions, possibly via H-bonding. Further, our findings establish a direct correlation between the mesoscale crystal growth kinetics and the molecular-scale interactions between organic additives and constituent ions, providing unprecedented insights for developing predictive syntheses of HOIP crystals with defined size, crystal habit and shape, and defect distribution.

36 MATERIALS SCIENCE↗

NASA NTRS1999

Novel Approaches to Cellular Transplantation from the US Space Program

Research in the treatment of type I diabetes is entering a new era that takes advantage of our knowledge in an ever increasing variety of scientific disciplines. Some may originate from very diverse sources, one of which is the Space Program at National Aeronautics and Space Administration (NASA). The Space Program contributes to diabetes-related research in several treatment modalities. As an ongoing effort for medical monitoring of personnel involved in space exploration activities NASA and the extramural scientific community investigate strategies for noninvasive estimation of blood glucose levels. Part of the effort in the space protein crystal growth program is high-resolution structural analysis insulin as a means to better understand the interaction with its receptor and with host immune components and as a basis for rational design of a "better" insulin molecule. The Space Program is also developing laser technology for potential early cataract detection as well as a noninvasive analyses for addressing preclinical diabetic retinopathy. Finally, NASA developed an exciting cell culture system that affords some unique advantages in the propagation and maintenance of mammalian cells in vitro. The cell culture system was originally designed to maintain cell suspensions with a minimum of hydrodynamic and mechanical sheer while awaiting launch into microgravity. Currently the commercially available NASA bioreactor (Synthecon, Inc., Houston, TX) is used as a research tool in basic and applied cell biology. In recent years there is continued strong interest in cellular transplantation as treatment for type I diabetes. The advantages are the potential for successful long-term amelioration and a minimum risk for morbidity in the event of rejection of the transplanted cells. The pathway to successful application of this strategy is accompanied by several substantial hurdles: (1) isolation and propagation of a suitable uniform donor cell population; (2) management of host immune rejection; (3) protection from the autoimmune component of the disease; and (4) anatomic placement of the engrafted cells that permits timely response to blood sugar levels as well as effective release and deployment of insulin. Bioreactor technology may provide some critical advances for surmounting some of these scientific hurdles. The NASA bioreactor is a horizontally rotating cylinder that is completely filled with culture medium. Gaseous exchange is maintained by a concentric cylinder of permeable silicon. In slow rotation (15-25 rpm) particles of small mass such as cells and tissue aggregates remain suspended in the rotating body of fluid. This novel approach suspends cells without stirring thus, allowing objects of different size mass to colocate and interact in a very low shear environment. In fact, analysis of the forces acting on individual cells in the rotating bioreactor reveals that the cells are continuously falling through the fluid medium. The conditions in the bioreactor permit assembly of cells into aggregates. three-dimensional tissue growth, synthesis of intercellular matrix,10 differentiation, and some sinusoid formations that may serve as a surrogate vasculature for larger tissue segments.

Pellis, Neal R.↗

NASA NTRS1994

Towards Understanding the Fluid Dynamic Phenomenon of Interest to Rocket Base Heating: A Review

The significance of the base heating problem for rockets during ascent is due to the complex interaction between the rocket nozzle plumes and the external-flow which can change the flow characteristics in the base region dramatically. At lower altitudes the external-flow merges with the plume-flow, without the formation of a large separated flow region, and the cooler external-flow promotes convective cooling of the base wall. Under these conditions the majority of the base heating is due to radiative heating from the shock heated plume gases. At higher altitudes, however, the process of base heating is not so straightforward. The plume and the base flow expands dramatically and separated flow regions occur in the base area. Hot exhaust gases from the rocket nozzle will be entrained into the separated flow regions and produce a convective component to the base wall heating. Further, if the rocket exhaust-gas contains soot, the soot can increase the emission from the gas and dramatically increase the wall absorption coefficient for radiative heating if it is deposited on the walls . In addition, if the rocket exhaust gas is fuel rich, the fuel can bum in the separated flow regions and further increase the base heating. The base burning phenomenon, and the increased base heating caused by it at higher altitudes, have been observed for the Space Shuttle and Saturn Rocket. Under these conditions, the total heating is significantly higher than the heating without separated flow in the base region, and the increase in heating is directly attributable to the fluid dynamic complexity of the base region. Realistic simulation of the base heating requires that the calculated flow environment reproduce the fluid dynamic flow features accurately. Thus, it will be necessary to introduce into the CFD codes the capability for the flow to respond to the complex vehicle geometry, the effect of turbulence, the ability to accurately reproduce the plume shock/shear layer structures and capture such phenomenon as shock induced base separation and base-burning phenomenon. A survey of experimental, theoretical and computational work that details the fluid dynamics of the base flow environment will be presented in the proposed paper. CFD simulations of rocket base flows using standard CFD codes such as OVERFLOW or GASP will be explored to capture these phenomenon accurately. Merits and limitations of these codes for base flow environment predictions will be explored.

Venkatapathy, E.↗

NASA NTRS1999

Double Diffusive Convection in Materials Processing

A great number of crystals grown in space are plagued by convective motions which contribute to structural flaws. The character of these instabilities is not well understood but is associated with density variations in the presence of residual gravity (g-jitter). As a specific example, past HgCdTe crystal growth space experiments by Lehoczky and co-workers indicate radial compositional asymmetry in the grown crystals. In the case of HgCdTe the rejected component into the melt upon solidification is HgTe which is denser than the melt. The space grown crystals indicate the presence of three dimensional flow with the heavier HgTe-rich material clearly aligned with the residual gravity (0.55-1.55 micro g) vector. This flow stems from double-diffusive convection, namely, thermal and solutal buoyancy driven flow in the melt. The study of double-diffusive convection is multi-faceted and rather vast. In our investigation, we seek to focus on one specific aspect of this discipline that is of direct relevance to materials processing especially crystal growth, namely, the side ways heating regime. This problem has been widely studied, both experimentally and numerically, in the context of solar ponds wherein the system is characterized by a linear salt (solutal) gradient with an imposed lateral temperature gradient. The induced flow instabilities arise from the wide disparity between the fluid thermal diffusivity and the solute diffusivity. The extension of the analysis to practical crystal growth applications has however not been rigorously made and understood. One subtle but important difference in crystal growth systems is the fact that die system solute gradient is non-linear (typically exponential). Besides, the crystal growth problem has the added complexities of solidification, both lateral and longitudinal thermal gradients and segregation phenomena in systems where binary and ternary compounds are being grown. This paper treats the side ways heating problem alone in a model fluid system. Results from detailed numerical calculations, mainly two dimensional are provided. The interactions between a non-linear solute gradient and an imposed transverse thermal gradient are investigated. The buoyancy effects are treated in the traditional Boussinesq approximation and also in a more complete density formulation to address recent concerns of the first approach especially in simulations of the system response in a reduced gravity environment. Detailed flow, temperature and solute field plots along with heat and mass transfer results are presented in the paper. Implications to practical crystal growth systems as discerned from the modeling results are also explored and reported.

Ramachandra, Narayanan↗

NASA NTRS2014

Statistical Mechanics of Turbulent Dynamos

Incompressible magnetohydrodynamic (MHD) turbulence and magnetic dynamos, which occur in magnetofluids with large fluid and magnetic Reynolds numbers, will be discussed. When Reynolds numbers are large and energy decays slowly, the distribution of energy with respect to length scale becomes quasi-stationary and MHD turbulence can be described statistically. In the limit of infinite Reynolds numbers, viscosity and resistivity become zero and if these values are used in the MHD equations ab initio, a model system called ideal MHD turbulence results. This model system is typically confined in simple geometries with some form of homogeneous boundary conditions, allowing for velocity and magnetic field to be represented by orthogonal function expansions. One advantage to this is that the coefficients of the expansions form a set of nonlinearly interacting variables whose behavior can be described by equilibrium statistical mechanics, i.e., by a canonical ensemble theory based on the global invariants (energy, cross helicity and magnetic helicity) of ideal MHD turbulence. Another advantage is that truncated expansions provide a finite dynamical system whose time evolution can be numerically simulated to test the predictions of the associated statistical mechanics. If ensemble predictions are the same as time averages, then the system is said to be ergodic; if not, the system is nonergodic. Although it had been implicitly assumed in the early days of ideal MHD statistical theory development that these finite dynamical systems were ergodic, numerical simulations provided sufficient evidence that they were, in fact, nonergodic. Specifically, while canonical ensemble theory predicted that expansion coefficients would be (i) zero-mean random variables with (ii) energy that decreased with length scale, it was found that although (ii) was correct, (i) was not and the expected ergodicity was broken. The exact cause of this broken ergodicity was explained, after much investigation, by greatly extending the statistical theory of ideal MHD turbulence. The mathematical details of broken ergodicity, in fact, give a quantitative explanation of how coherent structure, dynamic alignment and force-free states appear in turbulent magnetofluids. The relevance of these ideal results to real MHD turbulence occurs because broken ergodicity is most manifest in the ideal case at the largest length scales and it is in these largest scales that a real magnetofluid has the least dissipation, i.e., most closely approaches the behavior of an ideal magnetofluid. Furthermore, the effects grow stronger when cross and magnetic helicities grow large with respect to energy, and this is exactly what occurs with time in a real magnetofluid, where it is called selective decay. The relevance of these results found in ideal MHD turbulence theory to the real world is that they provide at least a qualitative explanation of why confined turbulent magnetofluids, such as the liquid iron that fills the Earth's outer core, produce stationary, large-scale magnetic fields, i.e., the geomagnetic field. These results should also apply to other planets as well as to plasma confinement devices on Earth and in space, and the effects should be manifest if Reynolds numbers are high enough and there is enough time for stationarity to occur, at least approximately. In the presentation, details will be given for both theoretical and numerical results, and references will be provided.

Shebalin, John V.↗

NASA NTRS2014

MHD Turbulence and Magnetic Dynamos

Shebalin, John V↗

NASA NTRS1990

On the instability of hypersonic flow past a flat plate

The instability of hypersonic boundary-layer flows over flat plates is considered. The viscosity of the fluid is taken to be governed by Sutherland's law, which gives a much more accurate representation of the temperature dependence of fluid viscosity at hypersonic speeds than Chapman's approximate linear law; although at lower speeds the temperature variation of the mean state is less pronounced so that the Chapman law can be used with some confidence. Attention is focussed on the so-called (vorticity) mode of instability of the viscous hypersonic boundary layer. This is thought to be the fastest growing inviscid disturbance at hypersonic speeds; it is also believed to have an asymptotically larger growth rate than any viscous or centrifugal instability. As a starting point the instability of the hypersonic boundary layer which exists far downstream from the leading edge of the plate is investigated. In this regime the shock that is attached to the leading edge of the plate plays no role, so that the basic boundary layer is non-interactive. It is shown that the vorticity mode of instability of this flow operates on a significantly different lengthscale than that obtained if a Chapman viscosity law is assumed. In particular, it is found that the growth rate predicted by a linear viscosity law overestimates the size of the growth rate by O(M(exp 2). Next, the development of the vorticity mode as the wavenumber decreases is described, and it is shown that acoustic modes emerge when the wavenumber has decreased from it's O(1) initial value to O(M (exp -3/2). Finally, the inviscid instability of the boundary layer near the leading edge in the interaction zone is discussed and particular attention is focussed on the strong interaction region which occurs sufficiently close to the leading edge. It is found that the vorticity mode in this regime is again unstable, and that it is concentrated in the transition layer at the edge of the boundary layer where the temperature adjusts from its large, O(M(exp 2), value in the viscous boundary layer, to its O(1) free stream value. The existence of the shock indirectly, but significantly, influences the instability problem by modifying the basic flow structure in this layer.

Blackaby, Nicholas↗

NASA NTRS2017

Vortex Generators in a Streamline-Traced, External-Compression Supersonic Inlet

Vortex generators within a streamline-traced, external-compression supersonic inlet for Mach 1.66 were investigated to determine their ability to increase total pressure recovery and reduce total pressure distortion. The vortex generators studied were rectangular vanes arranged in counter-rotating and co-rotating arrays. The vane geometric factors of interest included height, length, spacing, angle-of-incidence, and positions upstream and downstream of the inlet terminal shock. The flow through the inlet was simulated numerically through the solution of the steady-state, Reynolds-averaged Navier-Stokes equations on multi-block, structured grids using the Wind-US flow solver. The vanes were simulated using a vortex generator model. The inlet performance was characterized by the inlet total pressure recovery and the radial and circumferential total pressure distortion indices at the engine face. Design of experiments and statistical analysis methods were applied to quantify the effect of the geometric factors of the vanes and search for optimal vane arrays. Co-rotating vane arrays with negative angles-of-incidence positioned on the supersonic diffuser were effective in sweeping low-momentum flow from the top toward the sides of the subsonic diffuser. This distributed the low-momentum flow more evenly about the circumference of the subsonic diffuser and reduced distortion. Co-rotating vane arrays with negative angles-of-incidence or counter-rotating vane arrays positioned downstream of the terminal shock were effective in mixing higher-momentum flow with lower-momentum flow to increase recovery and decrease distortion. A strategy of combining a co-rotating vane array on the supersonic diffuser with a counter-rotating vane array on the subsonic diffuser was effective in increasing recovery and reducing distortion.

computational fluid dynamics↗

NASA NTRS2012

Characterization of Deficiencies in the Frequency Domain Forced Response Analysis Technique for Supersonic Turbine Bladed Disks

Turbine blades in rocket and jet engine turbomachinery experience enormous harmonic loading conditions. These loads result from the integer number of upstream and downstream stator vanes as well as the other turbine stages. Assessing the blade structural integrity is a complex task requiring an initial characterization of whether resonance is possible and then performing a forced response analysis if that condition is met. The standard technique for forced response analysis in rocket engine turbines is to decompose a computational fluid dynamics (CFD).generated flow field into its harmonic components, and to then perform a frequency response analysis at the problematic natural frequencies using cyclically symmetric structural dynamic models. Recent CFD analysis and water-flow testing at NASA/MSFC, though, indicates that this technique may miss substantial harmonic and non ]harmonic excitation sources that become present in complex flows. This complex content can only be captured by a CFD flow field encompassing at least an entire revolution. A substantial development effort to create a series of software programs to enable application of the 360 degree forcing function in a frequency response analysis on cyclic symmetric models has been completed (to be described in a future paper), but the question still remains whether the frequency response analysis itself is capable of capturing the excitation content sufficiently. Two studies comparing frequency response analysis with transient response analysis, therefore, of bladed-disks undergoing this complex flow environment have been performed. The first is of a bladed disk with each blade modeled by simple beam elements and the disk modeled with plates (using the finite element code MSC/NASTRAN). The focus of this model is to be representative of response of realistic bladed disks, and so the dimensions are roughly equivalent to the new J2X rocket engine 1st stage fuel pump turbine. The simplicity of the model allows the CFD load to be able to be readily applied, along with analytical and experimental variations in both the temporal and spatial fourier components of the excitation. In addition, this model is a first step in identifying response differences between transient and frequency forced response analysis techniques. The second phase assesses this difference for a much more realistic solid model of a bladed-disk in order to evaluate the effect of the spatial variation in loading on blade dominated modes. Neither research on the accuracy of the frequency response method when used in this context or a comprehensive study of the effect of test-observed variation on blade forced response have been found in the literature, so this research is a new contribution to practical structural dynamic analysis of gas turbines. The primary excitation of the upstream nozzles interacts with the blades on fuel pump of the J2X causes the 5th Nodal diameter modes to be excited, as explained by Tyler and Sofrin1, so a modal analysis was first performed on the beam/plate model and the 5ND bladed-disk mode at 40167 hz was identified and chosen to be the one excited at resonance (see figure 1). The first forced response analysis with this model focuses on identifying differences between frequency and transient response analyses. A hypothesis going into the analysis was that perhaps the frequency response was enforcing a temporal periodicity that did not really exist, and so therefore it would overestimate the response. As high dynamic response was a considerable source of stress in the J2X, examining this concept could potentially be beneficial for the program.

Brown, Andrew M.↗