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 829 records · Page 46

Grid generation and its application to separated flows

A method for the computation of recirculating flows on a nonorthogonal boundary-fitted coordinate system is presented. The algebraic grid generation technique is used to generate grids in a flow field of arbitrary shape. Computations are performed by using the set of governing equations transformed for use on a generalized curvilinear coordinate system. Turbulent flow computations were made by using the k-epsilon model for most of the flow cases. Three pressure correction algorithms - SIMPLE, SIMPLEC, and PISO - were incorporated for use on a nonorthogonal grid. Comparisons were made on the basis of numerical stability and accuracy. The present method was found to be efficient for computations of recirculating flows.

Maruszewski, J. P.↗

Active flutter suppression using dipole filters

By using traditional control concepts of gain root locus, the active suppression of a flutter mode of a flexible wing is examined. It is shown that the attraction of the unstable mode towards a critical system zero determines the degree to which the flutter mode can be stabilized. For control situations where the critical zero is adversely placed in the complex plane, a novel compensation scheme called a 'Dipole' filter is proposed. This filter ensures that the flutter mode is stabilized with acceptable control energy. The control strategy is illustrated by designing flutter suppression laws for an active flexible wing (AFW) wind-tunnel model, where minimal control effort solutions are mandated by control rate saturation problems caused by wind-tunnel turbulence.

Srinathkumar, S.↗

Direct simulations of compressible wall-bounded turbulence

Several direct numerical simulations of high-speed turbulent Couette flow were performed with a new spectral code. Mach numbers up to three and a Reynolds number of 3000 were used. A new time-integration scheme was developed to handle Mach numbers above 1.5, which require greater accuracy and stability than lower Mach numbers. At low Mach number, the large streamwise eddies found by M. J. Lee in high incompressible Couette flow simulations were reproduced. At higher Mach numbers these structures still exist, but they become considerably less organized (although the disorganization may be a function of the spanwise box size). While the same types of vortical structures seen in the incompressible flow are observed at higher Mach numbers, a new structure involving the divergence of the velocity is also observed. This structure is generally associated with low shear areas next to the walls, but it has not been determined whether it is a cause or an effect of the low shear. A 'nonphysical' simulation was performed to determine by what mechanism the Mach number affects the flow. It appears that pressure gradient (acoustic) effects are more important than variable viscosity effects in determining the wall shear, but the size of vortical structures is determined more by the local kinematic viscosity. Low-order mean statistics are provided to help quantify these effects.

Buell, Jeffrey C.↗

Droplet Cluster Behavior in Dense and Dilute Regions of a Spray

Concern regarding the efficiency, stability and safety margins of bi-propellant combustion in rocket engines has prompted the investigation of many specific aspects of spray evaporation, ignition and combustion previously not studied. The phenomenology of bi- propellant spray evaporation, ignition and combustion in a highly turbulent environment at elevated pressure is discussed.

liquid rocket propulsion bi-propellant combustion ↗

Measured Boundary Layer Transition and Rotor Hover Performance at Model Scale

An experiment involving a Mach-scaled, 11:08 f t: diameter rotor was performed in hover during the summer of 2016 at NASA Langley Research Center. The experiment investigated the hover performance as a function of the laminar to turbulent transition state of the boundary layer, including both natural and fixed transition cases. The boundary layer transition locations were measured on both the upper and lower aerodynamic surfaces simultaneously. The measurements were enabled by recent advances in infrared sensor sensitivity and stability. The infrared thermography measurement technique was enhanced by a paintable blade surface heater, as well as a new high-sensitivity long wave infrared camera. The measured transition locations showed extensive amounts, x=c>0:90, of laminar flow on the lower surface at moderate to high thrust (CT=s > 0:068) for the full blade radius. The upper surface showed large amounts, x=c > 0:50, of laminar flow at the blade tip for low thrust (CT=s < 0:045). The objective of this paper is to provide an experimental data set for comparisons to newly developed and implemented rotor boundary layer transition models in CFD and rotor design tools. The data is expected to be used as part of the AIAA Rotorcraft SimulationWorking Group

Overmeyer, Austin D.↗

Exploring data-driven modeling of boundary layer transition

Prediction of laminar-turbulent transition in boundary layer flows is an important component of predicting the aerodynamic performance of a number of aerospace configurations. According to the CFD Vision 2030 [1], transition modeling represents acriticalarea in CFD simulation capability that will remain a pacing item for the foreseeable future. The fact thattransition can take placevia either one of a myriad possible paths adds to the challenges inreliable transition predictions, despite a limited knowledge of the relevant input parameters. In the low disturbance environments typical of flight applications, transition is often initiated by small amplitude disturbances in the form of linear instability waves of the laminar boundary layer. These disturbances amplify linearly at first and eventually undergo a sequence of nonlinear interactions that result in transition to turbulence. Because the nonlinear phase is rather rapid, the amplification of boundary layer instabilities is governed by the linearstability theory over a majority of the distance leading up to the onset of transition. Semi-empirical transition correlations based on the linear stability theory have been successful in explaining the observed trends in transition location within a broad class of flows. However, the application of stability theory is highly non-robust and often requires a significant domain expertise. Recent work at the NASA Langley Research Center has beenaimed at bridging the gap between physics based transition analyses such as those based on linear stability theory and practical applications that require transition prediction by users that may not be well versed in transition physics. The applications of deep learning have been at the center of these efforts. This presentation will focus on the progress achieved thus far, highlighting the applications of neural networks to selectedtransition scenarios across a range of Mach numbers and flow configuration, as well as the lessons learnedand remaining challengeswithrespect to the selection of training data and neural networks architectures, hyperparameter tuning, and the physical insights distilled from the otherwise black-box models.

M. R. Malik↗

Boundary Layer Climatology at ARM Southern Great Plains

Operational since 1992, the Atmospheric Radiation Measurement (ARM) southern great plains (SGP) site at Oklahoma, USA has become a reference research site for meteorological studies. Due to an open data policy the ARM data are used by researchers all over the world. In this report, we review the long-term climatology of the atmospheric boundary layer, SGP instrumentations, the site and some site-specific atmospheric conditions which potentially effect wind turbines within the region. As the atmospheric boundary layer is bounded and influenced by the surface, observations of surface radiation components and heat fluxes are crucial in understanding land-atmosphere interactions. The entrainment of air, updrafts, downdrafts and boundary layer height characteristics is needed for understanding the structure and growth of the atmospheric boundary layer. Therefore, measurements from both surface in-situ and remote sensing observations at SGP provide an overall climatology and their interactions from surface up to the boundary layer. Measurements from a 60 m meteorological tower, surface flux stations, disdrometers, soil temperature and moisture flux plates, coherent Doppler lidar, Raman lidar, radiosondes, and satellite data at SGP central facility were analyzed. All the measurements were generally within a few square kilometers of each other at the central facility. This report focuses on data from January 2010 to June 2020 at SGP central facility. The various sections describe the ARM SGP site and surrounding wind turbines; in-situ and remote sensing instrumentation used in the report; provides mathematical equations to analyze fluxes, turbulence and other boundary layer parameters; a climatological analysis of surface winds, fluxes and thermodynamic parameters for several years; an analysis of observed winds in the framework of Monin-Obukhov Similarity Theory; an analysis of the boundary layer winds and direction from a Doppler lidar; multi-year turbulence estimates through the boundary layer from a Doppler lidar; atmospheric boundary layer water vapor and relative humidity profiles from Raman lidar; cloud base height and boundary layer height from multiple sensors and satellite data; and finally site specific atmospheric conditions, such as nocturnal low-level jets. Diurnal, seasonal and yearly variations of surface, sub-surface and boundary layer quantities, such as wind speed, direction, temperature, atmospheric stability, soil temperature, and various atmospheric fluxes at SGP showed distinct trends useful for focused modeling studies. The applicability of surface similarity theory on ARM SGP data is also evaluated, which showed northerly flows are aligned with MO theory estimates compared to southerly flows. Boundary layer winds and direction profiles for several years from a Doppler lidar shows a consistent presence of a nocturnal low-level jet and predominant southerly wind directions through the boundary layer at SGP. The inter-annual variability at SGP is low (<3.5%), with a mean annual wind speed of approximately 7 m s -1 at 100 m above ground level. Boundary layer turbulence and moisture transport from Doppler and Raman lidars are evaluated, which provides evidence of increased water vapor mass flux into the great plains during nocturnal low-level jets. The moisture flux from nocturnal low-level jets is observed to be maximum during summer periods. A novel machine learning algorithm is implemented to accurately estimate the planetary boundary layer height, providing further insights into growth and destruction of the convective boundary layer height during various seasons and land-atmosphere conditions. The high frequency of low-level clouds during winter, spring and fall seasons is validated using the multi-sensor array and satellite estimates of cloud top height. Satellite vegetative fraction data provides insight into seasonal surface roughness and vegetation variability around SGP site.

54 ENVIRONMENTAL SCIENCES↗

Large Eddy simulation of turbulence: A subgrid scale model including shear, vorticity, rotation, and buoyancy

The Reynolds numbers that characterize geophysical and astrophysical turbulence (Re approximately equals 10(exp 8) for the planetary boundary layer and Re approximately equals 10(exp 14) for the Sun's interior) are too large to allow a direct numerical simulation (DNS) of the fundamental Navier-Stokes and temperature equations. In fact, the spatial number of grid points N approximately Re(exp 9/4) exceeds the computational capability of today's supercomputers. Alternative treatments are the ensemble-time average approach, and/or the volume average approach. Since the first method (Reynolds stress approach) is largely analytical, the resulting turbulence equations entail manageable computational requirements and can thus be linked to a stellar evolutionary code or, in the geophysical case, to general circulation models. In the volume average approach, one carries out a large eddy simulation (LES) which resolves numerically the largest scales, while the unresolved scales must be treated theoretically with a subgrid scale model (SGS). Contrary to the ensemble average approach, the LES+SGS approach has considerable computational requirements. Even if this prevents (for the time being) a LES+SGS model to be linked to stellar or geophysical codes, it is still of the greatest relevance as an 'experimental tool' to be used, inter alia, to improve the parameterizations needed in the ensemble average approach. Such a methodology has been successfully adopted in studies of the convective planetary boundary layer. Experienc e with the LES+SGS approach from different fields has shown that its reliability depends on the healthiness of the SGS model for numerical stability as well as for physical completeness. At present, the most widely used SGS model, the Smagorinsky model, accounts for the effect of the shear induced by the large resolved scales on the unresolved scales but does not account for the effects of buoyancy, anisotropy, rotation, and stable stratification. The latter phenomenon, which affects both geophysical and astrophysical turbulence (e.g., oceanic structure and convective overshooting in stars), has been singularly difficult to account for in turbulence modeling. For example, the widely used model of Deardorff has not been confirmed by recent LES results. As of today, there is no SGS model capable of incorporating buoyancy, rotation, shear, anistropy, and stable stratification (gravity waves). In this paper, we construct such a model which we call CM (complete model). We also present a hierarchy of simpler algebraic models (called AM) of varying complexity. Finally, we present a set of models which are simplified even further (called SM), the simplest of which is the Smagorinsky-Lilly model. The incorporation of these models into the presently available LES codes should begin with the SM, to be followed by the AM and finally by the CM.

Canuto, V. M.↗

Spray combustion stability

The central purpose of this project is the improvement of liquid-fueled rocket motor design technology in order to assist the establishment of economical commercial access to space through the development of engines with enhanced performance and reliability. Specific research effort is focused on spray physics and associated combustion instability phenomena. Results concerning high pressure droplet gasification model, droplet turbulent dispersion model, and spray atomization model will contribute to the development of new computational tools for design of stable liquid propellant rocket engines.

Liang, Pak-Yan↗

Stably trapped proton fluxes in the Jovian magnetosphere

A model of the energetic proton fluxes in the Jovian magnetosphere is constructed based on the inward radial diffusion of protons from the solar wind and the plasma turbulent precipitation loss of protons from the radiation belts. Outside 12 Jovian radii the proton fluxes follow a loss-free radial diffusion profile. Inside 12 Jovian radii the proton fluxes should be near the stably trapped limit flux set by convective marginal stability to the electromagnetic ion cyclotron wave and the quasi-electrostatic ion loss-cone wave.

Coroniti, F. V.↗

Finite element analysis of low speed viscous and inviscid aerodynamic flows

A weak interaction solution algorithm was established for aerodynamic flow about an isolated airfoil. Finite element numerical methodology was applied to solution of each of differential equations governing potential flow, and viscous and turbulent boundary layer and wake flow downstream of the sharp trailing edge. The algorithm accounts for computed viscous displacement effects on the potential flow. Closure for turbulence was accomplished using both first and second order models. The COMOC finite element fluid mechanics computer program was modified to solve the identified equation systems for two dimensional flows. A numerical program was completed to determine factors affecting solution accuracy, convergence and stability for the combined potential, boundary layer, and parabolic Navier-Stokes equation systems. Good accuracy and convergence are demonstrated. Each solution is obtained within the identical finite element framework of COMOC.

Baker, A. J.↗

Implicit treatment of diffusion terms in lower-upper algorithms

A method is presented which allows diffusion terms to be treated implicitly in the lower-upper (LU) algorithm (which is a commonly used method for solving 'compressible' Euler and Navier-Stokes equations) so that the algorithm's good stability properties will not be impaired. The new method generalizes the concept of LU factorization from that associated with the sign of eigenvalues to that associated with backward- and forward-difference operators without regard to eigenvalues. The method is verified in a turbulent boundary layer study.

Shih, T. I.-P.↗

Vibroacoustic Tailoring of a Rod-Stiffened Composite Fuselage Panel with Multidisciplinary Considerations

An efficient multi-objective design tailoring procedure seeking to improve the vibroacoustic performance of a fuselage panel while maintaining or reducing weight is presented. The structure considered is the pultruded rod stitched efficient unitized structure, a highly integrated composite structure concept designed for a noncylindrical, next-generation flight vehicle fuselage. Modifications to a baseline design are evaluated within a six-parameter design space including spacing, flange width, and web height for both frame and stringer substructure components. The change in sound power radiation attributed to a design change is predicted using finite-element models sized and meshed for analyses in the 500 Hz, 1 kHz, and 2 kHz octave bands. Three design studies are carried out in parallel while considering a diffuse acoustic field excitation and two types of turbulent boundary-layer excitation. Kriging surrogate models are used to reduce the computational costs of resolving the vibroacoustic and weight objective Pareto fronts. The resulting Pareto optimal designs are then evaluated under a static pressurization ultimate load to assess structural strength and stability. Results suggest that choosing alternative configurations within the considered design space can reduce weight and improve vibroacoustic performance without compromising strength and stability of the structure under the static load condition considered, but the tradeoffs are significantly influenced by the spatial characteristics of the assumed excitation field.

Allen, Albert R.↗

Electron Jet of Asymmetric Reconnection

We present Magnetospheric Multiscale observations of an electron-scale current sheet and electron outflow jet for asymmetric reconnection with guide field at the subsolar magnetopause. The electron jet observed within the reconnection region has an electron Mach number of 0.35 and is associated with electron agyrotropy. The jet is unstable to an electrostatic instability which generates intense waves with E(sub parallel lines) amplitudes reaching up to 300 mV/m and potentials up to 20% of the electron thermal energy. We see evidence of interaction between the waves and the electron beam, leading to quick thermalization of the beam and stabilization of the instability. The wave phase speed is comparable to the ion thermal speed, suggesting that the instability is of Buneman type, and therefore introduces electron-ion drag and leads to braking of the electron flow. Our observations demonstrate that electrostatic turbulence plays an important role in the electron-scale physics of asymmetric reconnection.

Khotyaintsev, Yu. V.↗

Design and Modeling of an 80 bar Oxy-Combustor for Direct Fired Supercritical CO2 Applications

As supercritical CO2 power cycles for fossil energy power generation continue to generate interest, there are significant issues and unanswered questions regarding injector design, flame stabilization, wall heat transfer, CO emissions, combustion dynamics and other combustion phenomenon. For natural gas, direct-fired cycles with carbon capture it is believed that Computational Fluid Dynamics (CFD) modeling will play an essential role in the combustor design process. To accurately model turbulent reacting flows at these unique conditions, experimental data is needed to validate CFD codes and sub-models and is currently lacking at conditions relevant for these cycles. This paper presents the conceptual design and CFD simulations of an experimental 80 bar oxy-combustion facility and test article currently under construction at NETL. The facility is targeted towards the testing of a single injector, direct-fired sCO2 combustor at the 100 kW thermal output level. While these conditions do not reflect the actual Allam cycle operating conditions (300 bar) they are viewed as a stepping stone in the model validation process at supercritical conditions. Reynolds Averaged Navier-Stokes as well as Large Eddy Simulation are used to model the turbulent combustion process and aid in the design of the combustor and injector. Process parameters including oxidizer preheat temperature and combustor flowrate are investigated.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Editorial: Resolving atmospheric flow in complex environments: recent experiments in terrain and forest canopies

The characterization of atmospheric flows in complex environments, which may include steep terrain slopes and heterogeneous vegetation and/or forest cover, is a long-standing challenge in boundary-layer meteorology. Atmospheric observations are complicated by the presence of transient, terrain-induced flow features, forest-canopy-atmosphere interactions, and atmospheric stability effects, not to mention the logistical hurdles involved with instrument deployment, data analysis, and quality control. Furthermore, challenges in atmospheric modeling arise due to numerical errors associated with complex terrain flows, as well as reliance on simplified parameterizations for unresolved processes such as turbulent mixing and land-surface or forest-canopy-atmosphere interactions. These modeling challenges are exacerbated in the so-called “gray zone,” wherein features of interest have length scales that are similar to the model grid spacing, or when the principal flow layer is smaller than the grid spacing (e.g., slope flows).

54 ENVIRONMENTAL SCIENCES↗

SSME LOX post flow analysis/fluid structure interaction. Volume 1: Flow analysis

The realization of measures to improve the performance of the Space Shuttle is, to a large extent, dependent on an improved understanding of the fluid flow phenomena occurring in the main engine. The overall arrangement of the primary components of the Space Shuttle Main Engine (SSME) are presented. The impingement of the hot gases from the transfer ducts onto the LOX posts causes them to be subjected to severe thermal and gas dynamic loads, which in the past have resulted in the occasional breakage of some elements of the outer row of posts during test firings of the engine, particularly at higher power levels. Large velocities in the gap between the LOX posts may also be a contributing factor in causing dynamic stability problems. The deforming structural response of the posts to the pressure loading may likely affect the gas flowfield by producing a moving flowfield boundary, thereby creating a dynamically coupled unsteady fluid-structure system. The objective was to investigate the three-dimensional, turbulent flow around a simplified SSME LOX post array using an existing Reynolds averaged Navier-Stokes flow solver and a suitable turbulence model to parameterize the turbulent shear stresses. Numerical computations were performed to analyze the effect on the flowfield of varying the spacing between the LOX posts, which were modeled as rigid, three-dimensional circular cylinders. The methodology used in the computations is described. Results are presented.

Burke, Roger W.↗

Direct Numerical Simulation of Transition Due to Traveling Crossflow Vortices

Previous simulations of laminar breakdown mechanisms associated with stationary crossflow instability over a realistic swept-wing configuration are extended to investigate the alternate scenario of transition due to secondary instability of traveling crossflow modes. Earlier analyses based on secondary instability theory and parabolized stability equations have shown that this alternate scenario is viable when the initial amplitude of the most amplified mode of the traveling crossflow instability is greater than approximately 0.03 times the initial amplitude of the most amplified stationary mode. The linear growth predictions based on the secondary instability theory and parabolized stability equations agree well with the direct numerical simulation. Nonlinear effects are initially stabilizing but subsequently lead to a rapid growth followed by the onset of transition when the amplitude of the secondary disturbance exceeds a threshold value. Similar to the breakdown of stationary vortices, the transition zone is rather short and the boundary layer becomes completely turbulent across a distance of less than 15 times the boundary layer thickness at the completion of transition.

Li, Fei↗