Search NASA⌕ Search

SEARCH · Search NASA

Results for “gradient flow dynamics”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 55 records · Page 3

On the use of potential vorticity for the diagnosis of stratospheric synoptics

According to quasi-geostrophic approximation, regions of strong cross flow gradient of geopotential vortivity support Rossby wave propagation. Using potential vorticity as a dynamical tracer, the polar vortex can be considered as a material entity. In addition to synoptic charts of geopotential heights and temperatures, maps of potential vorticity are used to describe some distinct features of the circulation of the winter stratosphere.

Rose, K.↗

Lee waves: Benign and malignant

The flow of an incompressible fluid over an obstacle will produce an oscillation in which buoyancy is the restoring force, called a gravity wave. For disturbances of this scale, the atmosphere may be treated as dynamically incompressible, even though there exists a mean static upward density gradient. Even in the linear approximation - i.e., for small disturbances - this model explains a great many of the flow phenomena observed in the lee of mountains. However, nonlinearities do arise importantly, in three ways: (1) through amplification due to the decrease of mean density with height; (2) through the large (scaled) size of the obstacle, such as a mountain range; and (3) from dynamically singular levels in the fluid field. These effects produce a complicated array of phenomena - large departure of the streamlines from their equilibrium levels, high winds, generation of small scales, turbulence, etc. - that present hazards to aircraft and to lee surface areas. The nonlinear disturbances also interact with the larger-scale flow in such a manner as to impact global weather forecasts and the climatological momentum balance. If there is no dynamic barrier, these waves can penetrate vertically into the middle atmosphere (30-100 km), where recent observations show them to be of a length scale that must involve the coriolis force in any modeling. At these altitudes, the amplitude of the waves is very large, and the phenomena associated with these wave dynamics are being studied with a view to their potential impact on high performance aircraft, including the projected National Aerospace Plane (NASP). The presentation shows the results of analysis and of state-of-the-art numerical simulations, validated where possible by observational data, and illustrated with photographs from nature.

Wurtele, M. G.↗

Numerical simulation of aerothermal loads in hypersonic engine inlets due to shock impingement

The effect of shock impingement on an axial corner simulating the inlet of a hypersonic vehicle engine is modeled using a finite-difference procedure. A three-dimensional dynamic grid adaptation procedure is utilized to move the grids to regions with strong flow gradients. The adaptation procedure uses a grid relocation stencil that is valid at both the interior and boundary points of the finite-difference grid. A linear combination of spatial derivatives of specific flow variables, calculated with finite-element interpolation functions, are used as adaptation measures. This computational procedure is used to study laminar and turbulent Mach 6 flows in the axial corner. The description of flow physics and qualitative measures of heat transfer distributions on cowl and strut surfaces obtained from the analysis are compared with experimental observations. Conclusions are drawn regarding the capability of the numerical scheme for enhanced modeling of high-speed compressible flows.

Ramakrishnan, R.↗

Convective thermocapillary instabilities in liquid bridges

A steady shear flow is generated in an axisymmetric liquid bridge surrounded by a passive gas through the imposition of a temperature gradient along the bridge and the driving of the motion by thermocapillarity. This dynamic state is susceptible to convective instabilities, which can lead to propagating hydrothermal waves that feed on the underlying temperature gradients. The convective instabilities of axisymmetric return flow are presently given as functions of the liquid's Prandtl number and the surface Biot number of the interface.

Xu, J.-J.↗

Pick-up ion pressure gradients modulating the solar wind dynamics

Neutral interstellar atoms penetrate deeply into the inner heliosphere before they become ionized by various processes. As ions they are picked-up by the frozen-in magnetic fields and are convected outwards with the solar wind plasma. Thereby the primary plasma flow is mass, momentum, and energy-loaded. The dynamics of the distant multi-constituent solar wind is, however, not solely determined by these loading processes, but is also affected by the wave-mediated pick-up ion pressure gradients derivable from the pick-up ion distribution function. The action of the radial components of these pressures essentially counter balances the decelerating effect of the solar wind momentum loading, diminishing strongly the deceleration of the distant solar wind. Furthermore the latitudinal components of the pick-up ion pressures induce latitudinal forces acting on the multiconstituent solar plasma outflow and inducing nonradial bulk flow components. The enforced nonradial outflow geometry on the upwind hemisphere may partly be responsible for the magnetic flux deficit which was claimed since several years in the PIONEER-10 magnetic flux data.

Fahr, Hans J.↗

Performance of LI-1542 reusable surface insulation system in a hypersonic stream

The thermal and structural performance of a large panel of LI-1542 reusable surface insulation tiles was determined by a series of cyclic heating tests using radiant lamps and aerothemal tests in the Langley 8-foot high-temperature structures tunnel. Aerothermal tests were conducted at a free-stream Mach number of 6.6, a total temperature of 1830 K, Reynolds numbers of 2.0 and 4,900,000 per meter, and dynamic pressures of 29 and 65 kPa. The results suggest that pressure gradients in gaps and flow impingement on the header walls at the end of longitudinal gaps are sources for increased gap heating. Temperatures higher than surface radiation equilibrium temperature were measured deep in gaps and at the header walls. Also, the damage tolerance of the LI-1542 tiles appears to be very high. Tile edge erosion rate was slow; could not be tolerated in a shuttle application. Tiles soaked with water and subjected to rapid depressurization and aerodynamic heating showed no visible evidence of damage.

Hunt, L. R.↗

Analysis and Design of Crew Sleep Station for ISS

This paper details the analysis and design of the Temporary Sleep Station (TeSS) environmental control system for International Space Station (ISS). The TeSS will provide crewmembers with a private and personal space, to accommodate sleeping, donning and doffing of clothing, personal communication and performance of recreational activities. The need for privacy to accommodate these activities requires adequate ventilation inside the TeSS. This study considers whether temperature, carbon dioxide, and humidity within the TeSS remain within crew comfort and safety levels for various expected operating scenarios. Evaluation of these scenarios required the use and integration of various simulation codes. An approach was adapted for this study, whereby results from a particular code were integrated with other codes when necessary. Computational Fluid Dynamics (CFD) methods were used to evaluate the flow field inside the TeSS, from which local gradients for temperature, velocity, and species concentration such as CO (sub 2) could be determined. A model of the TeSS, containing a human, as well as equipment such as a laptop computer, was developed in FLUENT, a finite-volume code. Other factors, such as detailed analysis of the heat transfer through the structure, radiation, and air circulation from the TeSS to the US Laboratory Aisle, where the TeSS is housed, were considered in the model. A complementary model was developed in G189A, a code which has been used by NASA/JSC for environmental control systems analyses since the Apollo program. Boundary conditions were exchanged between the FLUENT and G189A TeSS models. G189A provides human respiration rates to the FLUENT model, while the FLUENT model provides local convective heat transfer coefficients to G189A model. An additional benefit from using an approach with both a systems simulation and CFD model, is the capability to verify the results of each model by comparison to the results of the other model. The G189A and FLUENT models were used to evaluate various ventilation designs for the TeSS over a range of operating conditions with varying crew metabolic load, equipment operating modes, ventilation flow rates, and with the TeSS doors open and closed. Results from the study were instrumental in the optimization of a design for the TeSS ventilation hardware. A special case was considered where failure of the TeSS ventilation system occurred. In this case, a study was conducted in order to determine the time required for the CO (sub 2) concentration inside the TeSS to increase to ISS limit values under transient conditions. A lumped-capacitance code, SINDA-FLUINT was used in this case to provide accurate predictions of the human reaction to the TeSS cabin conditions including core and skin temperatures and body heat storage. A simple two-dimensional CFD model of a crewmember inside the TeSS was developed in FLUENT in order to determine the volume envelope of the respired air from the human, which maintained a minimum velocity profile. This volume was then used in the SINDA-FLUINT model to facilitate the calculations of CO (sub 2) concentrations, dry bulb temperatures and humidity levels inside the TeSS.

Keener, John F.↗

Dynamic localization and second-order subgrid-scale models in large eddy simulations of channel flow

The objective here is to test the Dynamic Localization (DL) model in a wall-bounded channel flow for numerical stability and accuracy of results. Algebraic stress models suggest that the model for the residual subgrid-scale (SGS) Reynolds stress and scalar flux should generally have terms comprising most of the unique products of the resolved strain (S) and rotation (R) tensors with S and the resolved scalar gradient. The standard dynamic SGS model uses a simple (Smagorinsky) base model for the residual Reynolds stress, which is made proportional to S, and down-gradient base models for residual scalar fluxes; these correspond to the lowest, 'first-order' terms in algebraic stress models. Temporal scaling terms in these base models are formed from the magnitude of the resolved strain rate. While this is appropriate for simple shear flows, it may not be appropriate for more complicated flows (relevant to geophysical and astrophysical problems) that include any combination of shear, rotation, buoyancy, etc. On the other hand, the coefficient in the dynamic SGS model readily adjusts itself to different flow conditions and may adequately take account of these effects without the need for more complicated base models. Cabot (1993) has begun to test the dynamic SGS model in buoyant flows (Rayleigh-Benard and internally heated convection) with and without buoyancy terms explicitly included in the scaling terms of the base model; no great differences were found in large eddy simulation (LES) results for the different base model scalings. The second objective in this work is to test base models with additional, 'second-order' terms (e.g., S(sup 2) and RS for the residual Reynolds stress). These terms have been found to improve large-scale flow predictions by kappa-epsilon models in the presence of rotation and shear. Second-order base models will be tested here in the LES of channel flow with and without solid-body rotation and compared with results from the standard first-order base models to determine if there are significant differences or improvements in results that would warrant the added complexity of the second-order base models.

Cabot, William H.↗

Flowfield-Dependent Mixed Explicit-Implicit (FDMEL) Algorithm for Computational Fluid Dynamics

Despite significant achievements in computational fluid dynamics, there still remain many fluid flow phenomena not well understood. For example, the prediction of temperature distributions is inaccurate when temperature gradients are high, particularly in shock wave turbulent boundary layer interactions close to the wall. Complexities of fluid flow phenomena include transition to turbulence, relaminarization separated flows, transition between viscous and inviscid incompressible and compressible flows, among others, in all speed regimes. The purpose of this paper is to introduce a new approach, called the Flowfield-Dependent Mixed Explicit-Implicit (FDMEI) method, in an attempt to resolve these difficult issues in Computational Fluid Dynamics (CFD). In this process, a total of six implicitness parameters characteristic of the current flowfield are introduced. They are calculated from the current flowfield or changes of Mach numbers, Reynolds numbers, Peclet numbers, and Damkoehler numbers (if reacting) at each nodal point and time step. This implies that every nodal point or element is provided with different or unique numerical scheme according to their current flowfield situations, whether compressible, incompressible, viscous, inviscid, laminar, turbulent, reacting, or nonreacting. In this procedure, discontinuities or fluctuations of an variables between adjacent nodal points are determined accurately. If these implicitness parameters are fixed to certain numbers instead of being calculated from the flowfield information, then practically all currently available schemes of finite differences or finite elements arise as special cases. Some benchmark problems to be presented in this paper will show the validity, accuracy, and efficiency of the proposed methodology.

Garcia, S. M.↗

Upper-Ocean Heat Balance Processes and the Walker Circulation in CMIP5 Model Projections

Considerable uncertainty remains as to the importance of mechanisms governing decadal and longer variability of the Walker Circulation, its connection to the tropical climate system, and prospects for tropical climate change in the face of anthropogenic forcing. Most contemporary climate models suggest that in response to elevated CO2 and a warmer but more stratified atmosphere, the required upward mass flux in tropical convection will diminish along with the Walker component of the tropical mean circulation as well. Alternatively, there is also evidence to suggest that the shoaling and increased vertical stratification of the thermocline in the eastern Pacific will enable a muted SST increase there-- preserving or even enhancing some of the dynamical forcing for the Walker cell flow. Over the past decade there have been observational indications of an acceleration in near-surface easterlies, a strengthened Pacific zonal SST gradient, and globally-teleconnected dislocations in precipitation. But is this evidence in support of an ocean dynamical thermostat process posited to accompany anthropogenic forcing, or just residual decadal fluctuations associated with variations in warm and cold ENSO events and other stochastic forcing? From a modeling perspective we try to make headway on this question by examining zonal variations in surface energy fluxes and dynamics governing tropical upper ocean heat content evolution in the WCRP CMIP5 model projections. There is some diversity among model simulations; for example, the CCSM4 indicates net ocean warming over the IndoPacific region while the CSIRO model concentrates separate warming responses over the central Pacific and Indian Ocean regions. The models, as with observations, demonstrate strong local coupling between variations in column water vapor, downward surface longwave radiation and SST; but the spatial patterns of changes in the sign of this relationship differ among models and, for models as a whole, with observations. Our analysis focuses initially on probing the inter-model differences in energy fluxes / transports and Walker Circulation response to forcing. We then attempt to identify statistically the El Nino- / La Nina-related ocean heat content variability unique to each model and regress out the associated energy flux, ocean heat transport and Walker response on these shorter time scales for comparison to that of the anthropogenic signals.

Robertson, F. R.↗

Large eddy simulations of time-dependent and buoyancy-driven channel flows

The primary goal of this work has been to assess the performance of the dynamic SGS model in the large eddy simulation (LES) of channel flows in a variety of situations, viz., in temporal development of channel flow turned by a transverse pressure gradient and especially in buoyancy-driven turbulent flows such as Rayleigh-Benard and internally heated channel convection. For buoyancy-driven flows, there are additional buoyant terms that are possible in the base models, and one objective has been to determine if the dynamic SGS model results are sensitive to such terms. The ultimate goal is to determine the minimal base model needed in the dynamic SGS model to provide accurate results in flows with more complicated physical features. In addition, a program of direct numerical simulation (DNS) of fully compressible channel convection has been undertaken to determine stratification and compressibility effects. These simulations are intended to provide a comparative base for performing the LES of compressible (or highly stratified, pseudo-compressible) convection at high Reynolds number in the future.

Cabot, William H.↗

Aerobrake plasmadynamics - Macroscopic effects

The flow around an aerobraking spacecraft (such as the Aeroassist Flight Experiment reentry vehicle) will contain a region of partially ionized gas, that is, a plasma. It is shown here by numerical simulation that macroscopic plasmadynamic effects (which are not included in standard aerothermodynamic simulations) will have a noticeable effect on the reentry flow field. In particular, there are thermoelectric phenomena which can have a major influence in flow dynamics at the front of an ionizing bowshock. These thermoelectric phenomena arise because of the presence of large density and temperature gradients at the front of a reentry bowshock, and they include strong local magnetic fields, electric currents, and ohmic heating. One important result is the dramatic increase in temperature (over the case where plasma effects are neglected) at a reentry shock front; the implication is that macroscopic plasmadynamic effects can no longer be neglected in simulations of hypersonic reentry flow fields.

Shebalin, John V.↗

Treating convection in sequential solvers

The treatment of the convection terms in the sequential solver, a standard procedure found in virtually all pressure based algorithms, to compute the flow problems with sharp gradients and source terms is investigated. Both scalar model problems and one-dimensional gas dynamics equations have been used to study the various issues involved. Different approaches including the use of nonlinear filtering techniques and adoption of TVD type schemes have been investigated. Special treatments of the source terms such as pressure gradients and heat release have also been devised, yielding insight and improved accuracy of the numerical procedure adopted.

Shyy, Wei↗

Formulation for Simultaneous Aerodynamic Analysis and Design Optimization

An efficient approach for simultaneous aerodynamic analysis and design optimization is presented. This approach does not require the performance of many flow analyses at each design optimization step, which can be an expensive procedure. Thus, this approach brings us one step closer to meeting the challenge of incorporating computational fluid dynamic codes into gradient-based optimization techniques for aerodynamic design. An adjoint-variable method is introduced to nullify the effect of the increased number of design variables in the problem formulation. The method has been successfully tested on one-dimensional nozzle flow problems, including a sample problem with a normal shock. Implementations of the above algorithm are also presented that incorporate Newton iterations to secure a high-quality flow solution at the end of the design process. Implementations with iterative flow solvers are possible and will be required for large, multidimensional flow problems.

Hou, G. W.↗

(abstract) Three Dimensional Ice-Flow Velocity Estimation Using Satellite Radar Interferometry

An understanding of the flow dynamics of an ice sheet's outlet glaciers and ice streams requires knowledge of their flow velocity and strain rates (i.e., velocity gradients). Prior to the recent advent of satellite radar interferometry, it was not possible to measure detailed ice-flow velocity over the vast featureless areas that comprise most of the ice sheets. Since the launch of ERS-1, the use of satellite radar interferometry data for making densly sampled ice-flow velocity measurements has been firmly established by several studies. We have combined data from nonparallel orbits with surface slope information to make vector ice-flow measurements for the Ryder Glacier, Greenland. Our results for the Ryder are promising and indicate that repeat-pass interferometric data can be used to make vector measurements of ice velocity.

ice sheets radar interferometry glaciers flow dyna↗

The symmetric turbulent plane wake downstream of a sharp trailing edge

The analysis and modeling of the symmetric turbulent plane wake downstream of a sharp trailing edge is addressed. A compact description of the flow near the trailing edge is formulated using the results of a previous asymptotic analysis. The new description retains the two-layered structure identified in the previous work, and it clarifies the principal dynamics of the flow in the near-wake outer layer, away from the wake centerline. For zero-pressure-gradient flow, the near-wake outer layer is shown to be represented to leading order by the similarity solution that governs the outer region of the surface boundary layer. The leading perturbation in the outer layer due to the developing near-wake inner-layer flow is identified, and this is shown to be asymptotically smaller than undetermined higher-order terms associated with the surface boundary-layer flow. Results of the new near-wake analysis are used to formulate an algebraic eddy viscosity model for wake flow predictions at arbitrary distances from the trailing edge. The model is used in a numerical solution of the boundary layer equations, and computed velocity and Reynolds stress profiles are shown to compare well with experimental data.

Bogucz, E. A.↗

Second coefficient of viscosity in air

Acoustic attenuation measurements in air were analyzed in order to estimate the second coefficient of viscosity. Data over a temperature range of 11 C to 50 C and at relative humidities between 6 percent and 91 percent were used. This analysis showed that the second coefficient of viscosity varied between 1900 and 20,000 times larger than the dynamic or first coefficient of viscosity over the temperature and humidity range of the data. In addition, the data showed that the molecular relaxation effects, which are responsible for the magnitude of the second coefficient of viscosity, place severe limits on the use of time-independent, thermodynamic equations of state. Compressible flows containing large streamwise velocity gradients, like shock waves, which cause significant changes in particle properties to occur during time intervals shorter than hundredths of seconds, must be modeled using dynamic equations of state. The dynamic model approach is described briefly.

Ash, Robert L.↗

Electrohydrodynamically coupled Kelvin-Helmholtz instability

An inviscid model of the surface dynamics at the interface between two parallel-flowing insulating dielectric fluids, subject to a tangential electric field with a weak gradient, is developed. The two flows are oppositely directed, resulting in an electrohydrodynamically coupled Kelvin-Hemholtz instability. A simple stability criterion is used successfully to correlate experimental data obtained with a test configuration.

Perry, M. P.↗