Search NASA⌕ Search

SEARCH · Search NASA

Results for “Gradient flow”

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 487 records · Page 27

Microgravity Particle Dynamics

This research seeks to identify the experiment design parameters for future flight experiments to better resolve the effects of thermal and velocity gradients on gas-solid flows. By exploiting the reduced body forces and minimized thermal convection current of reduced gravity experiments, features of gas-solid flow normally masked by gravitationally induced effects can be studied using flow regimes unattainable under unigravity. This paper assesses the physical scales of velocity, length, time, thermal gradient magnitude, and velocity gradient magnitude likely to be involved in laminar gas-solid multiphase flight experiments for 1-100 micro-m particles.

Clark, Ivan O.↗

Calculation of turbulent boundary layers with heat transfer and pressure gradient utilizing a compressibility transformation. Part 1: Summary report

The analysis uses a compressibility transformation and utilizes higher order closure rules to complete the transformation. By requiring that the momentum equations in differential form be satisfied at the wall and at the sublayer edge, correspondence rules are obtained which relate the variable property (VP) flow to a constant property (CP) flow in which mass transfer and pressure gradient occur simultaneously. A new CP formulation is developed and numerical results for a variety of cases are presented. Comparisons with earlier forms of the transformation and with experiment are included. For the zero pressure gradient case some differences between the various predictions are observed. For the several pressure gradient cases examined, the results are found to be essentially identical to those given by first order closure rules; i.e., by a form of transformation which relates the VP flow to a CP flow with pressure gradient but zero mass transfer.

Economos, C.↗

Zonal Two Equation Kappa-Omega Turbulence Models for Aerodynamic Flows

Two new versions of the kappa-omega two-equation turbulence model will be presented. The new Baseline (BSL) model is designed to give results similar to those of the original kappa-omega model of Wilcox, but without its strong dependency on arbitrary freestream values. The BSL model is identical to the Wilcox model in the inner 50% of the boundary-layer but changes gradually to the standard kappa-epsilon model (in a kappa- omega formulation) towards the boundary-layer edge. The free shear layers. The second version of the model is called Shear-Stress Transport (SST) model. It is a variation of the BSL model with the additional ability to account for the transport of the principal turbulent shear stress in adverse pressure gradient boundary-layers. The model is based on Bradshaw's assumption that the principal shear-stress is proportional to the turbulent kinetic energy, which is introduced into the definition of the eddy-viscosity. Both models are tested for a large number of different flowfields. The results of the BSL model are similar to those of the original kappa-omega model, but without the undesirable freestream dependency. The predictions of the SST model are also independent of the freestream values but show better agreement with experimental data for adverse pressure gradient boundary-layer flows.

Menter, Florian R.↗

Pulsatile Newtonian frictional losses in a rigid tube.

Analytic expressions for the instantaneous average velocity and friction factor have been obtained for pulsatile Newtonian laminar flow with mean flow. Experimental pulsatile pressure gradients were found to be complex wave forms whose characteristics strongly influenced the experimental pulsatile average velocity and friction factor. Experimental and theoretical reduced instantaneous friction factor histories varied in a non-sinusoidal manner throughout the pulse cycle and exhibited extremum points. At high pulsation pressure ratios, substantial deviation occurred in the portion of the pulse cycle where the instantaneous average velocity approached zero. At low pulsation pressure ratios, the experimental and theoretical time average friction factor equaled unity and was independent of pulsation modulus. One potential application of this result is in the study of blood flow through small blood vessels.

Phillips, E. M.↗

Theoretical prediction of thick wing and pylon-fuselage-fanpod-nacelle aerodynamic characteristics at subcritical speeds. Part 1: Theory and results

The theoretical development and the comparison of results with data of a thick wing and pylon-fuselage-fanpod-nacelle analysis are presented. The analysis utilizes potential flow theory to compute the surface velocities and pressures, section lift and center of pressure, and the total configuration lift, moment, and vortex drag. The skin friction drag is also estimated in the analysis. The perturbation velocities induced by the wing and pylon, fuselage and fanpod, and nacelle are represented by source and vortex lattices, quadrilateral vortices, and source frustums, respectively. The strengths of these singularities are solved for simultaneously including all interference effects. The wing and pylon planforms, twists, cambers, and thickness distributions, and the fuselage and fanpod geometries can be arbitrary in shape, provided the surface gradients are smooth. The flow through nacelle is assumed to be axisymmetric. An axisymmetric center engine hub can also be included. The pylon and nacelle can be attached to the wing, fuselage, or fanpod.

Tulinius, J. R.↗

Improved two-equation k-omega turbulence models for aerodynamic flows

Two new versions of the k-omega two-equation turbulence model will be presented. The new Baseline (BSL) model is designed to give results similar to those of the original k-omega model of Wilcox, but without its strong dependency on arbitrary freestream values. The BSL model is identical to the Wilcox model in the inner 50 percent of the boundary-layer but changes gradually to the high Reynolds number Jones-Launder k-epsilon model (in a k-omega formulation) towards the boundary-layer edge. The new model is also virtually identical to the Jones-Lauder model for free shear layers. The second version of the model is called Shear-Stress Transport (SST) model. It is based on the BSL model, but has the additional ability to account for the transport of the principal shear stress in adverse pressure gradient boundary-layers. The model is based on Bradshaw's assumption that the principal shear stress is proportional to the turbulent kinetic energy, which is introduced into the definition of the eddy-viscosity. Both models are tested for a large number of different flowfields. The results of the BSL model are similar to those of the original k-omega model, but without the undesirable freestream dependency. The predictions of the SST model are also independent of the freestream values and show excellent agreement with experimental data for adverse pressure gradient boundary-layer flows.

Menter, Florian R.↗

Spatially Local Surrogate Modeling of Subgrid-Scale Effects in Idealized Atmospheric Flows: A Deep Learned Approach Using High-Resolution Simulation Data

Abstract We introduce a machine learned surrogate model from high-resolution simulation data to capture the subgrid-scale effects in dry, stratified atmospheric flows. We use deep neural networks (NNs) to model the spatially local state differences between a coarse-resolution simulation and a high-resolution simulation. The setup enables the capture of both dissipative and antidissipative effects in the state differences. The NN model is able to accurately capture the state differences in offline tests outside the training regime. In online tests intended for production use, the NN-coupled coarse simulation has higher accuracy over a significant period of time compared to the coarse-resolution simulation without any correction. We provide evidence of the capability of the NN model to accurately capture high-gradient regions in the flow field. With the accumulation of the errors, the NN-coupled simulation becomes computationally unstable after approximately 90 coarse simulation time steps. Insights gained from these surrogate models further pave the way for formulating stable, complex, physics-based spatially local NN models which are driven by traditional subgrid-scale turbulence closure models. Significance Statement Flows in the atmosphere are highly chaotic and turbulent, comprising flow structures of broad scales. For effective computational modeling of atmospheric flows, the effects of the small- and large-scale structures need to be captured by the simulations. Capturing the small-scale structures requires fine-resolution simulations. Even with the current state-of-the-art supercomputers, it can be prohibitively expensive to simulate these flows when computed for the entire earth over climate time scales. Thus, it is necessary to focus on the larger-scale structures using a coarse-resolution simulation while capturing the effects of the smaller-scale structures using some parameterization (approximation) scheme and incorporating it into the coarse-resolution simulation. We use machine learning to model the effects of the small-scale structures (subgrid-scale effects) in atmospheric flows. Data from a fine-resolution simulation is used to compute the missing subgrid-scale effects in coarse-resolution simulations. We then use machine learning models to approximate these differences between the coarse- and fine-resolution simulations. We see improved accuracy for the coarse-resolution simulations when corrected using these machine learned models.

54 ENVIRONMENTAL SCIENCES↗

Viscous shock-layer flowfield analysis by an explicit-implicit method

The present work extends the recently reported implicit analogue of MacCormack's earlier widely-used explicit method to external axisymmetric laminar flows with strong entropy gradients. The details of the 'numerics' of the implicit part are provided in a body-oriented coordinate system with a moving outer (shock) boundary during the transient part of the solutions. The limiting values of the Courant number are obtained when the shock boundary is treated explicitly. The solution algorithm outlined includes the treatment of the source term associated with the equations in weak conservation form. From the results obtained for two sample problems, it becomes clear that accuracy of predictions is, indeed, very good at higher values of the Courant number. There is a significant saving in overall computing time, depending on the Courant number used and the flow Reynolds number. These properties combined with the simplicity of programming the implicit analogue may appeal to researchers for using it in the analysis of three-dimensional flow problems.

Gupta, R. N.↗

A viscous shock-layer flowfield analysis by an explicit-implicit method

An implicit analogue of a widely used explicit method to external axisymmetric laminar flows with strong entropy gradients is extended. The details of the "numerics" of the implicit part are provided in a body oriented coordinate system with a moving outer (shock) boundary during the transient part of the solutions. The limiting values of the Courant number are obtained when the shock boundary is treated explicitly. The solution algorithm outlined includes the treatment of the source term associated with the equations in weak conservation form. From the results obtained for two sample problems, it becomes clear that accuracy of predictions is, indeed, very good at higher values of the Courant number. There is a significant saving in overall computing time, depending on the Courant number used and the flow Reynolds number. These properties combined with the simplicity of programming the implicit analog may appeal to researchers for using it in the analysis of 3-D flow problems.

Gupta, R. M.↗

Flexible Micropost Arrays for Shear Stress Measurement

Increased fuel costs, heightened environmental protection requirements, and noise abatement continue to place drag reduction at the forefront of aerospace research priorities. Unfortunately, shortfalls still exist in the fundamental understanding of boundary-layer airflow over aerodynamic surfaces, especially regarding drag arising from skin friction. For example, there is insufficient availability of instrumentation to adequately characterize complex flows with strong pressure gradients, heat transfer, wall mass flux, three-dimensionality, separation, shock waves, and transient phenomena. One example is the acoustic liner efficacy on aircraft engine nacelle walls. Active measurement of shear stress in boundary layer airflow would enable a better understanding of how aircraft structure and flight dynamics affect skin friction. Current shear stress measurement techniques suffer from reliability, complexity, and airflow disruption, thereby compromising resultant shear stress data. The state-of-the-art for shear stress sensing uses indirect or direct measurement techniques. Indirect measurements (e.g., hot-wire, heat flux gages, oil interferometry, laser Doppler anemometry, small scale pressure drag surfaces, i.e., fences) require intricate knowledge of the studied flow, restrictive instrument arrangements, large surface areas, flow disruption, or seeding material; with smaller, higher bandwidth probes under development. Direct measurements involve strain displacement of a sensor element and require no prior knowledge of the flow. Unfortunately, conventional "floating" recessed components for direct measurements are mm to cm in size. Whispering gallery mode devices and Fiber Bragg Gratings are examples of recent additions to this type of sensor with much smaller (m) sensor components. Direct detection techniques are often single point measurements and difficult to calibrate and implement in wind tunnel experiments. In addition, the wiring, packaging, and installation of delicate micro-electromechanical devices impede the use of most direct shear sensors. Similarly, the cavity required for sensing element displacement is sensitive to particulate obstruction. This work was focused on developing a shear stress sensor for use in subsonic wind tunnel test facilities applicable to an array of test configurations. The non-displacement shear sensors described here have minimal packaging requirements resulting in minimal or no disturbance of boundary layer flow. Compared to previous concepts, device installation could be simple with reduced cost and down-time. The novelty lies in the creation of low profile (nanoscale to 100 μm) micropost arrays that stay within the viscous sub-layer of the airflow. Aerodynamic forces, which are related to the surface shear stress, cause post deflection and optical property changes. Ultimately, a reliable, accurate shear stress sensor that does not disrupt the airflow has the potential to provide high value data for flow physics researchers, aerodynamicists, and aircraft manufacturers leading to greater flight efficiency arising from more in-depth knowledge on how aircraft design impacts near surface properties.

Wohl, Christopher J.↗

Nonideal stability analysis of differentially rotating plasmas with global curvature effects

The linear stability of global nonaxisymmetric modes in differentially rotating, magnetized, nonideal plasma is critical to classifying turbulence and transport phenomena. We investigate the competition between the local magneto-rotational instability (MRI) and the magneto-curvature instability (MCI)—a distinct nonaxisymmetric low-frequency curvature-driven global branch that appears alongside MRI. Here, to accomplish this, we developed a nonideal global spectral method, which is validated against NIMROD code simulations. This spectral approach allows for the direct derivation of an extended effective potential formalism and a resistive Alfvénic resonance condition, providing a framework for direct analysis of energy contributions and confinement mechanisms. Our study reveals that the global, low-frequency MCI persists at low magnetic Reynolds numbers (Rm), whereas the localized, high-frequency MRI is stabilized by diffusive broadening of its structure around its Alfvénic resonances. Consequently, we identify the global MCI branch as the primary onset mechanism for nonaxisymmetric magnetohydrodynamic instability in systems with finite curvature, e.g., astrophysical rotators. We establish distinct parameter regimes for mode dominance: MCI prevails in geometrically moderate-thickness disks with intermediate curvature and radial gaps, while MRI dominates in thin, low-curvature disks with large radial gaps. Mode competition is also highly sensitive to the flow profile, particularly vorticity and its gradient, with nonuniform shear profiles exhibiting more robust instability due to flow curvature (i.e., the second derivative of the flow profile) and shear contributions. A key outcome is the development of spectral diagrams derived from the global spectral method. These diagrams comprehensively map dominant instabilities and their characteristics, offering a predictive tool for critical onset parameters (i.e., flow curvature, magnetic field, and Rm) and facilitating the interpretation of experimental and simulation results. Notably, these diagrams demonstrate that the global MCI is generally the sole unstable mode at the initial onset of nonaxisymmetric instability.

Haywood, Alexander [Princeton Univ., NJ (United St↗

Development of new flux splitting schemes

Maximizing both accuracy and efficiency has been the primary objective in designing a numerical algorithm for computational fluid dynamics (CFD). This is especially important for solutions of complex three dimensional systems of Navier-Stokes equations which often include turbulence modeling and chemistry effects. Recently, upwind schemes have been well received for their capability in resolving discontinuities. With this in mind, presented are two new flux splitting techniques for upwind differencing. The first method is based on High-Order Polynomial Expansions (HOPE) of the mass flux vector. The second new flux splitting is based on the Advection Upwind Splitting Method (AUSM). The calculation of the hypersonic conical flow demonstrates the accuracy of the splitting in resolving the flow in the presence of strong gradients. A second series of tests involving the two dimensional inviscid flow over a NACA 0012 airfoil demonstrates the ability of the AUSM to resolve the shock discontinuity at transonic speed. A third case calculates a series of supersonic flows over a circular cylinder. Finally, the fourth case deals with tests of a two dimensional shock wave/boundary layer interaction.

Liou, Meng-Sing↗

Passivation of Flexible YBCO Superconducting Current Lead With Amorphous SiO2 Layer

Adiabatic demagnetization refrigerators (ADR) are operated in space to cool detectors of cosmic radiation to a few 10s of mK. A key element of the ADR is a superconducting magnet operating at about 0.3 K that is continually energized and de-energized in synchronism with a thermal switch, such that a piece of paramagnetic salt is alternately warm in a high magnetic field and cold in zero magnetic field. This causes the salt pill or refrigerant to cool, and it is able to suck heat from an object, e.g., the sensor, to be cooled. Current has to be fed into and out of the magnets from a dissipative power supply at the ambient temperature of the spacecraft. The current leads that link the magnets to the power supply inevitably conduct a significant amount of heat into the colder regions of the supporting cryostat, resulting in the need for larger, heavier, and more powerful supporting refrigerators. The aim of this project was to design and construct high-temperature superconductor (HTS) leads from YBCO (yttrium barium copper oxide) composite conductors to reduce the heat load significantly in the temperature regime below the critical temperature of YBCO. The magnet lead does not have to support current in the event that the YBCO ceases to be superconducting. Cus - tomarily, a normal metal conductor in parallel with the YBCO is a necessary part of the lead structure to allow for this upset condition; however, for this application, the normal metal can be dispensed with. Amorphous silicon dioxide is deposited directly onto the surface of YBCO, which resides on a flexible substrate. The silicon dioxide protects the YBCO from chemically reacting with atmospheric water and carbon dioxide, thus preserving the superconducting properties of the YBCO. The customary protective coating for flexible YBCO conductors is silver or a silver/gold alloy, which conducts heat many orders of magnitude better than SiO2 and so limits the use of such a composite conductor for passing current across a thermal gradient with as little flow of heat as possible to make an efficient current lead. By protecting YBCO on a flexible substrate of low thermal conductivity with SiO2, a thermally efficient and flexible current lead can be fabricated. The technology is also applicable to current leads for 4 K superconducting electronics current biasing. A commercially available thin-film YBCO composite tape conductor is first stripped of its protective silver coating. It is then mounted on a jig that holds the sample flat and acts as a heat sink. Silicon dioxide is then deposited onto the YBCO to a thickness of about 1 micron using PECVD (plasma-enhanced chemical vapor deposition), without heating the YBCO to the point where degradation occurs. Since SiO2 can have good high-frequency electrical properties, it can be used to coat YBCO cable structures used to feed RF signals across temperature gradients. The prime embodiment concerns the conduction of DC current across the cryogenic temperature gradient. The coating is hard and electrically insulating, but flexible.

Johannes, Daniel↗

The lagRST Model: A Turbulence Model for Non-Equilibrium Flows

This study presents a new class of turbulence model designed for wall bounded, high Reynolds number flows with separation. The model addresses deficiencies seen in the modeling of nonequilibrium turbulent flows. These flows generally have variable adverse pressure gradients which cause the turbulent quantities to react at a finite rate to changes in the mean flow quantities. This "lag" in the response of the turbulent quantities can t be modeled by most standard turbulence models, which are designed to model equilibrium turbulent boundary layers. The model presented uses a standard 2-equation model as the baseline for turbulent equilibrium calculations, but adds transport equations to account directly for non-equilibrium effects in the Reynolds Stress Tensor (RST) that are seen in large pressure gradients involving shock waves and separation. Comparisons are made to several standard turbulence modeling validation cases, including an incompressible boundary layer (both neutral and adverse pressure gradients), an incompressible mixing layer and a transonic bump flow. In addition, a hypersonic Shock Wave Turbulent Boundary Layer Interaction with separation is assessed along with a transonic capsule flow. Results show a substantial improvement over the baseline models for transonic separated flows. The results are mixed for the SWTBLI flows assessed. Separation predictions are not as good as the baseline models, but the over prediction of the peak heat flux downstream of the reattachment shock that plagues many models is reduced.

Lillard, Randolph P.↗

Laminar flow control is maturing

Recent research demonstrates that laminar flow (LF) can be reliable in flight and that the support system need not be complex. Shaping produces favorable pressure gradients for maintaining natural laminar flow (NLF), and laminar flow control (LFC) techniques such as full chord suction promise a fuel-saving payoff of up to 30 percent on long-range missions. For large aircraft, current research is concentrated on hybrid LFC concepts which combine suction and pressure-gradient control. At NASA Ames, an F-14 with variable wing sweep has been flight tested with smooth surface gloves on the wings; preliminary results indicate high transition Reynolds numbers to sweep angles as large as 25 deg. In addition, a 757 was flight tested with an NLF glove on the right wing just outboard of the engine pylon; and the LF was found to be suprisingly robust.

Wagner, Richard D.↗

The upper-branch stability of compressible boundary layer flows

The upper-branch linear and nonlinear stability of compressible boundary layer flows is studied using the approach of Smith and Bodonyi (1982) for a similar incompressible problem. Both pressure gradient boundary layers and Blasius flow are considered with and without heat transfer, and the neutral eigenrelations incorporating compressibility effects are obtained explicitly. The compressible nonlinear viscous critical layer equations are derived and solved numerically and the results indicate some solutions with positive phase shift across the critical layer. Various limiting cases are investigated including the case of much larger disturbance amplitudes and this indicates the structure for the strongly nonlinear critical layer of the Benney-Bergeon (1969) type. It is also shown how a match with the inviscid neutral inflexional modes arising from the generalized inflexion point criterion, is achieved.

Gajjar, J. S. B.↗

Aero-acoustics of Drag Generating Swirling Exhaust Flows

Aircraft on approach in high-drag and high-lift configuration create unsteady flow structures which inherently generate noise. For devices such as flaps, spoilers and the undercarriage there is a strong correlation between overall noise and drag such that, in the quest for quieter aircraft, one challenge is to generate drag at low noise levels. This paper presents a rigorous aero-acoustic assessment of a novel drag concept. The idea is that a swirling exhaust flow can yield a steady, and thus relatively quiet, streamwise vortex which is supported by a radial pressure gradient responsible for pressure drag. Flows with swirl are naturally limited by instabilities such as vortex breakdown. The paper presents a first aero-acoustic assessment of ram pressure driven swirling exhaust flows and their associated instabilities. The technical approach combines an in-depth aerodynamic analysis, plausibility arguments to qualitatively describe the nature of acoustic sources, and detailed, quantitative acoustic measurements using a medium aperture directional microphone array in combination with a previously established Deconvolution Approach for Mapping of Acoustic Sources (DAMAS). A model scale engine nacelle with stationary swirl vanes was designed and tested in the NASA Langley Quiet Flow Facility at a full-scale approach Mach number of 0.17. The analysis shows that the acoustic signature is comprised of quadrupole-type turbulent mixing noise of the swirling core flow and scattering noise from vane boundary layers and turbulent eddies of the burst vortex structure near sharp edges. The exposed edges are the nacelle and pylon trailing edge and the centerbody supporting the vanes. For the highest stable swirl angle setting a nacelle area based drag coefficient of 0.8 was achieved with a full-scale Overall Sound Pressure Level (OASPL) of about 40dBA at the ICAO approach certification point.

Shah, P. N.↗

Implicit Large-Eddy Simulations of Zero-Pressure Gradient, Turbulent Boundary Layer

A set of direct simulations of zero-pressure gradient, turbulent boundary layer flows are conducted using various span widths (62-630 wall units), to document their influence on the generated turbulence. The FDL3DI code that solves compressible Navier-Stokes equations using high-order compact-difference scheme and filter, with the standard recycling/rescaling method of turbulence generation, is used. Results are analyzed at two different Re values (500 and 1,400), and compared with spectral DNS data. They show that a minimum span width is required for the mere initiation of numerical turbulence. Narrower domains ((is) less than 100 w.u.) result in relaminarization. Wider spans ((is) greater than 600 w.u.) are required for the turbulent statistics to match reference DNS. The upper-wall boundary condition for this setup spawns marginal deviations in the mean velocity and Reynolds stress profiles, particularly in the buffer region.

Turblent boundary layer↗