Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Turbulence Stabilization”

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 271 records · Page 15

On the energy principle and ion tearing in the magnetotail

We re-examine the use of the energy principle as applied to the tearing instability in the magnetotail. We demonstrate that when a magnetic field component normal to the current sheet is present, electron pitch-angle diffusion (PAD) either by micro-turbulence or by chaotic orbits cannot remove the strong stabilization of the tearing mode caused by electron compressibility. We find that our conclusions are in agreement with those of Pellat et al. (1991), who argued on the basis of canonical P(sub y) conservation that the stabilization of the ion tearing mode cannot be removed by the introduction of PAD. Our results are at variance with those of Kuznetsova and Zelenyi (1991), who argued that the application of the energy principle used by Pellat et al. (1991) is incorrect, and that tearing is in fact unstable in the limit of strong PAD. We show that the disagreement between these two studies can be traced to an incorrect orbit evaluation first introduced by Coroniti (1980) and subnsequently used by Kuznetsova and Zelenyi (1991).

Brittnacher, M.↗

Evaluation of the Antarctic Boundary Layer Thermodynamic Structure In MERRA2 Using Dropsonde Observations from the Concordiasi Campaign

Recent high-resolution dropsonde observations from the 2010 Concordiasi field campaign in austral spring season show that surface-based inversions (SBIs) over Antarctica are frequently eroded, with well-mixed boundary layers occurring 33% and 18% of the time in West and East Antarctica, respectively. In this study, using the dropsonde observations, we evaluate the performance of the Modern-Era Retrospective analysis for Research and Applications, version 2 (MERRA2) in representing the Antarctic boundary layer thermodynamic structure. Results show that MERRA2 has a good overall representation of the Antarctic surface stability and correctly predicts 82% of the SBIs. However, an under prediction of less stable boundary layer occurrence, especially over the elevated East Antarctic plateau, is favored during conditions of increased lower tropospheric stability associated with model dynamics, indicating difficulty in parameterizing turbulence in very stable boundary layers. In addition, a lower tropospheric cool bias (first model level and above) is observed in the MERRA2 reanalysis, especially over West Antarctica, which amplifies in the boundary layer during mixed conditions. The near-surface cold bias is most pronounced when the model fails to predict mixed layers over West Antarctica, and is expected to negatively impact the representation of surface energy budget and melt processes.

Manisha Ganeshan↗

Distinct Surface Response to Black Carbon Aerosols

For the radiative impact of individual climate forcings, most previous studies focused on the global mean values at the top of the atmosphere (TOA), and less attention has been paid to surface processes, especially for black carbon (BC) aerosols. In this study, the surface radiative responses to five different forcing agents were analyzed by using idealized model simulations. Our analyses reveal that for greenhouse gases, solar irradiance, and scattering aerosols, the surface temperature changes are mainly dictated by the changes of surface radiative heating, but for BC, surface energy redistribution between different components plays a more crucial role. Globally, when a unit BC forcing is imposed at TOA, the net shortwave radiation at the surface decreases by −5.87±0.67 W m −2 (W m −2 ) −1 (averaged over global land without Antarctica), which is partially offset by increased downward longwave radiation (2.32±0.38 W m −2 (W m −2 ) −1 from the warmer atmosphere, causing a net decrease in the incoming downward surface radiation of −3.56±0.60 W m −2 (W m −2 ) −1 . Despite a reduction in the downward radiation energy, the surface air temperature still increases by 0.25±0.08 K because of less efficient energy dissipation, manifested by reduced surface sensible (−2.88±0.43 W m −2 (W m −2 ) −1 ) and latent heat flux (−1.54±0.27 W m −2 (W m −2 ) −1 ), as well as a decrease in Bowen ratio (−0.20±0.07 (W m −2 ) −1 ). Such reductions of turbulent fluxes can be largely explained by enhanced air stability (0.07±0.02 K (W m −2 ) −1 ), measured as the difference of the potential temperature between 925 hPa and surface, and reduced surface wind speed (−0.05±0.01 m s −1 (W m −2 ) −1 ). The enhanced stability is due to the faster atmospheric warming relative to the surface, whereas the reduced wind speed can be partially explained by enhanced stability and reduced Equator-to-pole atmospheric temperature gradient. These rapid adjustments under BC forcing occur in the lower atmosphere and propagate downward to influence the surface energy redistribution and thus surface temperature response, which is not observed under greenhouse gases or scattering aerosols. Our study provides new insights into the impact of absorbing aerosols on surface energy balance and surface temperature response.

radiative impact↗

Time analysis of polydisperse sprays in complex turbulent environments

The two-phase flow field of a small swirl-stabilized burner has been investigated, using a phase Doppler particle analyzer to obtain the velocity, drop size, and time of arrival components. Water/kerosene comparisons were made for the case of nonreacting flow. It is found that the spray flow field is strongly influenced by the aerodynamic flow field and by the reaction and fluid properties. It is noted that vortex shedding produces clusters of drops.

Rudoff, R. C.↗

High-Order Stabilized Finite Elements on Dynamic Meshes

The development of dynamic mesh capability for turbulent flow simulations using the Streamlined Upwind Petrov-Galerkin (SUPG) discretization is described. The current work extends previous research to include high-order spatial accuracy, including the satisfaction of the discrete geometric conservation law (GCL) on curved elements. Two closely-related schemes are described and the ability of these schemes to satisfy the GCL, while also maintaining temporal accuracy and conservation is assessed. Studies indicate that although one scheme discretizes the time derivative in conservative form, both schemes exhibit temporal conservation errors that decrease according to the expected design order of accuracy. The source of the temporal conservation errors is examined, and it is demonstrated that many finite-volume and finite-element schemes can also be expected to have difficulty strictly satisfying conservation in time. The effects on conservation are examined and, while present in the simulations, are seen to be negligible for the problems considered.

Anderson, W. Kyle↗

Analogies between transitional and turbulent boundary layers

One of the interesting aspects of transitional and turbulent boundary layers is the development of counter-rotating streamwise vortices near the wall. The most regular pattern is found in a boundary layer on a concave wall where the generation mechanism is known to be the Goertler instability. The origin of these vortices in other translational and turbulent boundary layers is presently unknown. Since the counter-rotating vortices are located in a region of strong shear, low-speed fluid is pumped away from the wall which coalesces into regions of low momentum lying between the vortices. As this pumping action continues, localized inflectional velocity profiles become apparent in the transitional and turbulent boundary layers. The oscillations which develop upon these profiles scale with the local thickness and velocity difference in the same manner as the two-dimensional steady free shear layer stability problems. The oscillations grow to large amplitude and break down into new turbulence in both the transitional and turbulent boundary layers.

Blackwelder, R. F.↗

Aerodynamic Parameter Estimation Using Reconstructed Turbulence Measurements

A classical method for reconstructing atmospheric turbulence from onboard measurements of airdata and inertial sensors was improved. The reconstructed turbulence was included in a system identification analysis to estimate nondimensional stability and control derivatives in a longitudinal short period model using the maximum likelihood equation-error method with Fourier transform data. Practical aspects of the approach are discussed, such as reconstruction accuracy, data collinearity, model structure, and real-time estimation. Results using simulation and flight test data for a subscale airplane indicated that accurate turbulence reconstructions and parameter estimates could be obtained from maneuvering flight in high levels of turbulence.

Jared A Grauer↗

Transition Prediction for the 30P30N Airfoil and the CRM-HL using FUN3D

The ability to accurately model transition over high-lift configurations can lead to improved agreement between stationary Reynolds-Averaged Navier-Stokes (RANS) computations and experimental measurements. In this work, selected transport-equation-based RANS models for transition prediction in the NASA FUN3D flow solver are used to simulate flow over the 30P30N three-element airfoil and the High-Lift Common Research Model (CRM-HL). The primary RANS-based transition model is the Langtry-Menter (LM) 𝜸-Re𝜽𝒕 model coupled with the Menter shear-stress transport model. Additionally, the existing capability for transition prediction based on linear stability correlations is used to allow the modeling of transition over multiple elements of the 30P30N configuration. This approach pairs an algebraic transition model with the Spalart-Allmaras (SA) turbulence model that relies on an imposed transition location in FUN3D and the parabolized stability equations as implemented in the LASTRAC code. Predictions based on these different transition prediction approaches and models are compared with each other, experimental data, and previous simulations by evaluating the lift and drag coefficients, the surface-pressure/skin-friction distributions, and the locations of transition onset. Results are presented for multiple angles of attack for the 30P30N airfoil and CRM-HL. For the CRM-HL, a grid resolution study is performed for the LM transition model and the SA turbulence model.

High Lift↗

The internal structure of the convective envelope of the sun

An axisymmetric hydrodynamic model of the sun consisting of a core at radiant equilibrium and a convective envelope is studied with consideration of gravitational forces, rotation, turbulence effects, and electromagnetic forces. The dimensionless parameters that control the equations of plasma motion within the sun and the boundary conditions are investigated; estimates are given for these parameters on the basis of experimental data and the existing theoretical models. Stability conditions are derived for convective equilibrium, developed cellular convection, and developed turbulent convection with consideration of rotation and turbulent Reynolds stresses. These conditions characterize the possible physical state of the matter in the convective envelope. Using estimates of the dimensionless parameters and models of some of them, a method of successive approximations is proposed for solution of the problem as a whole.

Yavorskaya, I. M.↗

Analysis of Instabilities in Non-Axisymmetric Hypersonic Boundary Layers Over Cones

Hypersonic flows over circular cones constitute one of the most important generic configurations for fundamental aerodynamic and aerothermodynamic studies. In this paper, numerical computations are carried out for Mach 6 flows over a 7-degree half-angle cone with two different flow incidence angles and a compression cone with a large concave curvature. Instability wave and transition-related flow physics are investigated using a series of advanced stability methods ranging from conventional linear stability theory (LST) and a higher-fidelity linear and nonlinear parabolized stability equations (PSE), to the 2D eigenvalue analysis based on partial differential equations. Computed N factor distribution pertinent to various instability mechanisms over the cone surface provides initial assessments of possible transition fronts and a guide to corresponding disturbance characteristics such as frequency and azimuthal wave numbers. It is also shown that strong secondary instability that eventually leads to transition to turbulence can be simulated very efficiently using a combination of advanced stability methods described above.

Li, Fei↗

Critical Evaluation of the Conference on Low Reynolds Number Aerodynamics at the University of Notre Dame

As a result of the continued interest in designing efficient low Reynolds number systems, the University of Notre Dame decided to organize a Conference on Low Reynolds Number Aerodynamics in June 1989. This Conference followed the 1986 International Conference in London by about three years and the first Conference on Low Reynolds Number Airfoil Aerodynamics at Notre Dame in 1985 by four years. The emphasis of the 1989 Conference was to assess the state-of-the-art in the chord Reynolds number range from about 10,000 to about 700,000. Applications of current interest include high altitude remotely or robotically piloted vehicles, ultra-light, and human powered vehicles as well as mini-RPVs at low altitudes. Other examples include small axial-flow fans used to cool electronic equipment in the unpressurized sections of high-altitude aircraft and gas turbine blades. High Reynolds number airfoil design strategies attempt to control the onset and development of turbulent boundary layers. This is difficult at low Reynolds numbers because of the increased stability of attached laminar boundary layers. Therefore, laminar separation is common even at small angles of attack at low Reynolds numbers. Under these conditions, the development of a turbulent boundary layer usually depends on the formation of a transitional separation bubble. The purpose of this Conference on Low Reynolds Number Aerodynamics was to bring together those researchers who have been active in areas closely related to this subject. It is clear from the papers presented that a great deal of progress has been made in understanding the occurrence and behavior of laminar separation and transition as well as their overall effect on the performance of airfoils at low chord Reynolds numbers. This progress has brought us closer to our goal of improving analytical methods for the design and evaluation of a variety of practical applications.

Mueller, Thomas J.↗

Laminar to Turbulence Transition in Boundary Layers due to Tripping Devices

Direct numerical simulations (DNS)and linear stability analysis were performed to identify the mechanisms by which two-and three-dimensional trips promote early transition in subsonic boundary layers over a flat plate. Investigations were conducted for two boundary layer transition trip configurations: one is a two-dimensional rectangular-shaped rod and the other is a three-dimensional zigzag tape. Reynolds numbers based on the trip heights are 683 and 300 for the two-and three-dimensional trips, respectively. Simulations with the two-dimensional trip showed that the flow separates at the top of the trip and forms a long recirculation zone downstream. However, the flow remains steady in the simulation without any external disturbances. Linear stability analysis showed that the inflectional profile in the recirculation zone is strongly unstable to frequencies in the range of 100-300 Hz. Simulations with the external disturbances forced at the frequency of 200 Hz showed that disturbances grow strongly inside the recirculation zone as observed in the Klebanoff1 experiment.

Transition↗

Crossflow Instability on a Wedge-Cone at Mach 3.5

As a follow-on activity to the HyBoLT flight experiment, a six degree half angle wedge-cone model at zero angle of attack has been employed to experimentally and computationally study the boundary layer crossflow instability at Mach 3.5 under low disturbance freestream conditions. Computed meanflow and linear stability analysis results are presented along with corresponding experimental Pitot probe data. Using a model-mounted probe survey apparatus, data acquired to date show a well defined stationary crossflow vortex pattern on the flat wedge surface. This effort paves the way for additional detailed, calibrated flow field measurements of the crossflow instability, both stationary and traveling modes, and transition-to-turbulence under quiet flow conditions as a means of validating existing stability theory and providing a foundation for dynamic flight instrumentation development.

Beeler, George B.↗

Intermittency Model for Coupled CFD-Stability Transition Analysis in Hypersonic Flow

Accurate prediction of aerothermodynamic loads on hypersonic vehicles requires precise modeling of surface quantities across the boundary layer transition zone. The peak heating loads and total heat transfer are determined by parameters such as the transition zone length and the magnitude of potential overshoots in heat flux and skin friction beyond their respective values in fully turbulent flows. While previous studies on CFD integrated transition modeling using linear stability correlations have shown promise in modeling these features, they did not develop an intermittency model for hypersonic flows. This paper presents a data-driven approach to develop a model of this type by correlating experimental transition data with the relevant flow parameters. The resulting model demonstrates significant improvements over previous low-speed models in terms of predicting heat transfer distributions during the transition process associated with first and second mode instabilities in axisymmetric high-speed flows. Separate correlations for flight and ground test conditions are developed, and the potential to combine these correlations is discussed.

Transition↗

Non-Turbulent Liquid-Bearing Polar Clouds: Observed Frequency of Occurrence and Simulated Sensitivity to Gravity Waves

A common feature of polar liquid-bearing clouds (LBCs) is radiatively-driven turbulence, which may variously alter cloud lifecycle via vertical mixing, droplet activation, and subsequent feedbacks. However, polar LBCs are commonly initiated under stable, non-turbulent conditions. Using long-term data from the North Slope of Alaska and McMurdo, Antarctica, we show that non-turbulent conditions prevail in ~25% of detected LBCs, surmised to be preferentially early in their lifecycle. We conclude that non-turbulent LBCs are likely common over the polar regions owing primarily to atmospheric temperature and stability. Such stable environments are known to support gravity wave activity. Using large-eddy simulations we find that short to intermediate period gravity waves may catalyze turbulence formation when aerosol particles available for activation are sufficiently small. We posit that the frequent occurrence of non-turbulent LBCs over the polar regions has implications for polar aerosol-cloud interactions and their parameterization in large-scale models.

Clouds↗

Hi-alpha forebody design. Part 1: Methodology base and initial parametrics

The use of Computational Fluid Dynamics (CFD) has been investigated for the analysis and design of aircraft forebodies at high angle of attack combined with sideslip. The results of the investigation show that CFD has reached a level of development where computational methods can be used for high angle of attack aerodynamic design. The classic wind tunnel experiment for the F-5A forebody directional stability has been reproduced computationally over an angle of attack range from 10 degrees to 45 degrees, and good agreement with experimental data was obtained. Computations have also been made at combined angle of attack and sideslip over a chine forebody, demonstrating the qualitative features of the flow, although not producing good agreement with measured experimental pressure distributions. The computations were performed using the code known as cfl3D for both the Euler equations and the Reynolds equations using a form of the Baldwin-Lomax turbulence model. To study the relation between forebody shape and directional stability characteristics, a generic parametric forebody model has been defined which provides a simple analytic math model with flexibility to capture the key shape characteristics of the entire range of forebodies of interest, including chines.

Mason, William H.↗

A compressible Navier-Stokes code for turbulent flow modeling

An implicit, finite volume code for solving two dimensional, compressible turbulent flows is described. Second order upwind differencing of the inviscid terms of the equations is used to enhance stability and accuracy. A diagonal form of the implicit algorithm is used to improve efficiency. Several zero and two equation turbulence models are incorporated to study their impact on overall flow modeling accuracy. Applications to external and internal flows are discussed.

Coakley, T. J.↗

The stability of laminar flow past a sphere

As a contribution to the problem of turbulence on a surface of rotation, the method of small oscillation is applied to the flow past a sphere. It was found that the method developed for two-dimensional flow is applicable without modifications. The frictional layer in the vicinity of the stagnation point of a surface of rotation is less stable against small two-dimensional disturbances than in the stagnation point itself, as proved from an analysis of the velocity distribution made by Homann.

Pretsch, J↗