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 865 records · Page 48
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Describing the Mechanism of Instability Suppression Using a Central Pilot Flame With Coupled Experiments and Simulations

Pilot flames are commonly used to extend combustor operability limits and suppress combustion oscillations in low-emissions gas turbines. Combustion oscillations, a coupling between heat release rate oscillations and combustor acoustics, can arise at the operability limits of low-emissions combustors where the flame is more susceptible to perturbations. While the use of pilot flames is common in land-based gas turbine combustors, the mechanism by which they suppress instability is still unclear. In this study, we consider the impact of a central jet pilot on the stability of a swirl-stabilized flame in a variable-length, single-nozzle combustor. Previously, the pilot flame was found to suppress the instability for a range of equivalence ratios and combustor lengths. We hypothesize that combustion oscillation suppression by the pilot occurs because the pilot provides hot gases to the vortex breakdown region of the flow that recirculate and improve the static, and hence dynamic, stability of the main flame. This hypothesis is based on a series of experimental results that show that pilot efficacy is a strong function of pilot equivalence ratio but not pilot flow rate, which would indicate that the temperature of the pilot products as well as the combustion intensity of the pilot flame play more of a role in oscillation stabilization than the length of the pilot flame relative to the main flame. Further, the pilot-flame efficacy increases with pilot-flame equivalence ratio until it matches the main-flame equivalence ratio; at pilot equivalence ratios greater than the main equivalence ratio, the pilot-flame efficacy does not change significantly with pilot equivalence ratio. To understand these results, we use large-eddy simulation (LES) to provide a detailed analysis of the flow in the region of the pilot flame and the transport of radical species in the region between the main flame and pilot flame. Furthermore, the simulation, using a flamelet/progress variable-based chemistry tabulation approach and standard eddy viscosity/diffusivity turbulence closure models, provides detailed information that is inaccessible through experimental measurements.

42 ENGINEERING↗
DOE OSTI2024

Implementation and Validation of a Generalized Actuator Disk Parameterization for Wind Turbine Simulations Within the FastEddy Model

ABSTRACT Fast and accurate large‐eddy simulation (LES) of the atmospheric boundary layer plays a crucial role in advancing wind energy research. Long‐duration wind farm studies at turbine‐resolving scales have become increasingly important to understand the intricate interactions between large wind farms and the atmospheric boundary layer. However, the prohibitive computational cost of these turbulence‐ and turbine‐resolving simulations has precluded such modeling to be exercised on a regular basis. To that end, we implement and validate the generalized actuator disk (GAD) model in the computationally efficient, graphics processing unit (GPU)–resident, LES model FastEddy. We perform single‐turbine simulations under three atmospheric stabilities (neutral, unstable, and stable) and compare them against observations from the Scaled Wind Farm Technology (SWiFT) facility and other LES codes from the recent Wakebench turbine wake model benchmark. Our idealized LES results agree well with observed wake velocity deficit and downstream recovery across stability regimes. Turbine response in terms of rotational speed, generated power, torque, and thrust coefficient are well predicted across stability regimes and are consistent with the LES results from the benchmark. The FastEddy simulations are found to be at least two orders of magnitude more efficient than the traditional CPU‐based LES models, opening the door for realistic LES simulations of full wind plants as a viable standard practice.

17 WIND ENERGY↗
DOE OSTI2025

Optimization and active stabilization of a far-infrared laser for NSTX-U high poloidal wavenumber scattering diagnostics

The far-infrared (FIR) laser output beam power and profile are important parameters in the laser-aided diagnostics, directly influencing the spatial resolution and signal-to-noise ratio of measurements. Here, this work focuses on developing a systematic control method to enhance FIR laser beam quality through optimized mirror alignment and real-time feedback-based precision cavity length tuning. A 150 W CO 2 laser, aligned with the waveguide axis using a HeNe reference laser, serves as the pump source. The sensitivity of FIR beam intensity to pump gas pressure and thermal expansion is investigated, revealing that even a 1 µm cavity expansion can significantly degrade output power stability to about two-thirds of its original value. To address this, a feedback control module has been designed and implemented for active cavity length adjustment, stabilizing the output power at ∼30 mW. In addition, maintaining a high formic acid gas pressure ($>$190 mTorr) within the cavity ensures reliable operation. The optimized FIR laser will be deployed on the National Spherical Torus eXperiment-U high poloidal wavenumber scattering system for studying electron-scale turbulence in tokamak plasmas.

Xu, Xinhang [Univ. of California, Davis, CA (Unite↗
NASA NTRS1982

On current aspects of finite element computational fluid mechanics for turbulent flows

A set of nonlinear partial differential equations suitable for the description of a class of turbulent three-dimensional flow fields in select geometries is identified. On the basis of the concept of enforcing a penalty constraint to ensure accurate accounting of ordering effects, a finite element numerical solution algorithm is established for the equation set and the theoretical aspects of accuracy, convergence and stability are identified and quantized. Hypermatrix constructions are used to formulate the reduction of the computational aspects of the theory to practice. The robustness of the algorithm, and the computer program embodiment, have been verified for pertinent flow configurations.

Baker, A. J.↗
NASA NTRS1986

Low-disturbance wind tunnels

During the past years, there was an extensive program under way at the Langley Research Center to upgrade the flow quality in several of the large wind tunnels. This effort has resulted in significant improvements in flow quality in these tunnels and has also increased the understanding of how and where changes in existing and new wind tunnels are most likely to yield the desired improvements. As part of this ongoing program, flow disturbance levels and spectra were measured in several Langley tunnels before and after modifications were made to reduce acoustic and vorticity fluctuations. A brief description of these disturbance control features is given for the Low-Turbulence Pressure Tunnel, the 4 x 7 Meter Tunnel, and the 8 Foot Transonic Pressure Tunnel. To illustrate typical reductions in disturbance levels obtained in these tunnels, data from hot-wire or acoustic sensors are presented. A concept for a subsonic quiet tunnel designed to study boundary layer stability and transition is also presented. Techniques developed at Langley in recent years to eliminate the high intensity and high-frequency acoustic disturbances present in all previous supersonic wind tunnels are described. In conclusion, the low-disturbance levels present in atmospheric flight can now be simulated in wind tunnels over the speed range from low subsonic through high supersonic.

Beckwith, I. E.↗
NASA NTRS1996

Simulation of Three-Dimensional Symmetric and Asymmetric Instabilities in Attachment-Line Boundary Layers

On a swept wing, contamination along the leading edge, Tollmien-Schlichting waves, stationary or traveling crossflow vortices, and/or Taylor-Gortler vortices can cause the catastrophic breakdown of laminar to turbulent flow, which leads to increased skin-friction drag for the aircraft. The discussion in this Note will be limited to disturbances which evolve along the attachment line (leading edge of swept wing). If the Reynolds number of the attachment-line boundary layer is greater than some critical value, then the complete wing is inevitably engulfed in turbulent flow. Essentially, there are two critical Reynolds number points that must be considered. The first is for small-amplitude disturbances, and the second is for bypass transition. The present study will use direct numerical simulations to validate a linear 2D-eigenvalue prediction method based on parabolized stability equations by Lin and Malik. This method is considered because it suggests that a number of symmetric and asymmetric modes exist and are stable or unstable on the attachment line depending on the Reynolds number. If validated, the approach would predict a number of modes which are linearly damped in the Reynolds number regime 100 to 245; however, these modes may grow nonlinearly and provide an explanation to this region.

Joslin, Ronald D.↗
NASA NTRS2019

Effect of 3D Roughness Patch on Instability Amplification in a Supersonic Boundary Layer

Surface roughness is known to have a substantial impact on the aerothermodynamic loading of high-speed vehicles, particularly via its influence on the laminar-turbulent transition process within the boundary layer. Numerical simulations are performed to investigate the effects of a distributed region of densely packed, sinusoidal shape roughness elements on a Mach 3.5 flat plate boundary layer for flow conditions corresponding to the planned conditions of an upcoming experiment in the Mach 3.5 Supersonic Low Disturbance Tunnel at the NASA Langley Research Center. Analysis of convective instabilities in the wake of the roughness patch was reported in a previous paper and the current work extends that analysis to instability amplification across the length of the roughness patch. Quasi-parallel stability analysis of the modified boundary layer flow over the patch indicates two dominant families of unstable disturbances, namely, a group of high frequency modes that amplify in localized regions along the roughness patch and another group of low frequency modes that have smaller peak amplification rates but amplify steadily both above the roughness patch and in the wake region behind it. The results suggest that the amplification factors associated with the high-frequency modes are sufficiently low, at least for the roughness patches considered in this paper, so that these modes are unlikely to have a major influence on the transition process. The amplification of the low-frequency modes within the region of the roughness patch is further quantified via direct numerical simulations. Results confirm the strongly destabilizing influence of the roughness patch on the first mode instabilities, yielding an N-factor increment of ΔN ≈ 3.6 for a roughness patch length of eight wavelengths.

Laminar turbulent transition↗
DOE OSTI2025

Evolution of intermittent filaments in the scrape-off layer of NSTX

Filamentary structures naturally arise from background turbulence in the scrape-off layer (SOL) of plasmas, leading to significant particle and heat transport that can degrade overall plasma confinement. This enhanced transport can contribute to unacceptably high heat loads on plasma-facing components. As such, understanding the physics of SOL plasma filaments is critical for predicting and mitigating their effects in future fusion devices. On the National Spherical Torus Experiment (NSTX), plasma filaments—commonly referred to as “blobs”—were investigated using the gas-puff imaging (GPI) diagnostic in the edge and SOL regions. The analysis involved identifying, segmenting, and tracking the characteristic contours of the blobs in each frame of the GPI video sequences. Their evolution was characterized through shape descriptors, velocity, and angular velocity derived from their contour coordinates. The results indicate that as the blob area increases, their shapes become more concave and less circular, suggesting reduced structural stability in larger blobs. This result aligns with previous theoretical results where it was shown that larger blobs are more susceptible to instabilities [Krasheninnikov et al., J. Plasma Phys. 74, 679–717 (2008) and D'Ippolito et al., Phys. Plasmas 18, 060501 (2011)]. A positive correlation was observed between radial velocity and radial position, suggesting radially outward acceleration of the filaments, potentially driven by decreasing viscous drag toward the far SOL. Interestingly, blobs in background SOL turbulence exhibited minimal spinning in contrast to filaments originating from edge localized modes, which show substantial rotation during their paths [Lampert et al., Phys. Plasmas 29, 102502 (2022)]. Statistical analysis of the solidity and total curvature shape descriptors, along with their temporal evolution, revealed relatively broad, near-Gaussian distributions. This suggests that blob morphology is strongly influenced by stochastic turbulent processes in the surrounding plasma environment. Blob parameters were also compared with bulk plasma and radial profile measurements. Notable trends were found between blob rotation and poloidal velocity with collisionality and line-integrated density. These findings contribute to a deeper understanding of blob dynamics and provide valuable insights for refining SOL turbulence models.

Covariance and correlation↗
NASA NTRS1981

Instability and transition in rotating disk flow

The stability of three-dimensional rotating disk flow is investigated, including the effects of Coriolis forces and streamline curvature. The numerical results show that the critical Reynolds number for establishment of stationary vortex flow is 287. These vortices spiral outward at an angle of about 11.2 deg, and transition to turbulence occurs when their total amplification is about e to the 11th. New experimental results are also reported on the spatial growth rates of the stationary vortices. It is shown that the analysis gives growth rates that compare much better with the experimental results than do results obtained using the Orr-Sommerfeld equation. The experimental results tend to support the numerical prediction that the number of stationary vortices varies directly with the Reynolds number. The calculations also indicate the existence of weakly unstable propagating (type II) modes at low Reynolds numbers (Critical Reynolds Number being approximately equal to 49).

Malik, M. R.↗
NASA NTRS1986

Kelvin-Helmholtz instability at the magnetopause boundary

An MHD model is defined for the convective Kelvin-Helmholtz instability (CKHI) which satellite data indicates produces vortices in the near-earth region of the plasma sheet through turbulence caused by shear flow. The MHD equations are solved as an initial value problem, which reveals that the CKHI growth rate is the same as that of the periodic KHI (PKHI). However, the CKHI is not as limited in the range of amplitude as is the PKHI, which experiences stability in nonlinear growth because of the periodic characteristics of the growth. A perpendicularly moving CKHI can produce a large vortex flow and shocks off the vortex boundary. The satellite data indicate that the vortices are in the size range of 10-20 earth radii. Additional effects can include the expulsion of magnetic flux and the formation of the boundary layer. The success of the model up to the point of formation of the boundary layer is taken as an incentive for further model development using more realistic magnetospheric data.

Wu, C. C.↗
NASA NTRS1992

A fast, uncoupled, compressible, two-dimensional, unsteady boundary layer algorithm with separation for engine inlets

A finite difference boundary layer algorithm was developed to model viscous effects when an inviscid core flow solution is given. This algorithm solved each boundary layer equation separately, then iterated to find a solution. Solving the boundary layer equations sequentially was 2.4 to 4.0 times faster than solving the boundary layer equations simultaneously. This algorithm used a modified Baldwin-Lomax turbulence model, a weighted average of forward and backward differencing of the pressure gradient, and a backward sweep of the pressure. With these modifications, the boundary layer algorithm was able to model flows with and without separation. The number of grid points used in the boundary layer algorithm affected the stability of the algorithm as well as the accuracy of the predictions of friction coefficients and momentum thicknesses. Results of this boundary layer algorithm compared well with experimental observations of friction coefficients and momentum thicknesses. In addition, when used interactively with an inviscid flow algorithm, this boundary layer algorithm corrected for viscous effects to give a good match with experimental observations for pressures in a supersonic inlet.

Roach, Robert L.↗
NASA NTRS1992

A fast, uncoupled, compressible, two-dimensional, unsteady boundary layer algorithm with separation for engine inlets

A finite difference boundary layer algorithm was developed to model viscous effects when an inviscid core flow solution is given. This algorithm solved each boundary layer equation separately, then iterated to find a solution. Solving the boundary layer equations sequentially was 2.4 to 4.0 times faster than solving the boundary layer equations simultaneously. This algorithm used a modified Baldwin-Lomax turbulence model, a weighted average of forward and backward differencing of the pressure gradient, and a backward sweep of the pressure. With these modifications, the boundary layer algorithm was able to model flows with and without separation. The number of grid points used in the boundary layer algorithm affected the stability of the algorithm affected the stability of the algorithm as well as the accuracy of the predictions of friction coefficients and momentum thicknesses. Results of this boundary layer algorithm compared well with experimental observations of friction coefficients and momentum thicknesses. In addition, when used interactively with an inviscid flow algorithm, this boundary layer algorithm corrected for viscous effects to give a good match with experimental observations for pressures in a supersonic inlet.

Roach, Robert L.↗
NASA NTRS1949

Exact Calculation of Laminar Boundary Layer in Longitudinal Flow over a Flat Plate with Homogeneous Suction

Lately it has been proposed to reduce the friction drag of a body in a flow for the technically important large Reynolds numbers by the following expedient: the boundary layer, normally turbulent, is artificially kept laminar up to high Reynolds numbers by suction. The reduction in friction drag thus obtained is of the order of magnitude of 60 to 80 percent of the turbulent friction drag, since the latter, for large Reynolds numbers, is several times the laminar friction drag. In considering the idea mentioned one has first to consider whether suction is a possible means of keeping the boundary layer laminar. This question can be answered by a theoretical investigation of the stability of the laminar boundary layer with suction. A knowledge, as accurate as possible, of the velocity distribution in the laminar boundary layer with suction forms the starting point for the stability investigation. E. Schlichting recently gave a survey of the present state of calculation of the laminar boundary layer with suction.

Iglisch, Rudolf↗
NASA NTRS1989

A spectral multi-domain technique applied to stability and transition of incompressible and compressible flows

Techniques for the numerical simulation of compressible and incompressible viscous flows are described and demonstrated. The focus is on a multidomain spectral method based on the flux-balance domain-interface scheme of Macaraeg and Streett (1986). The fundamental principles of the method are reviewed, and it is shown that the multidomain spectral approach facilitates the treatment of instability and the transition to turbulence. Typical numerical results for the Taylor-Couette problem and for bounded compressible high-speed mixing layers are presented in graphs and briefly characterized.

Macaraeg, Michele↗
DOE OSTI2020

Mesh-Sequenced Realizations for Evaluation of Subgrid-Scale Models for Turbulent Combustion

This paper develops a new approach for analysis of subgrid closures for turbulent combustion as modeled using direct quadrature (finite-rate chemistry) techniques. The approach, termed multiresolution analysis through mesh-sequenced realizations (MRA-MSR), conducts simultaneous, constrained large-eddy simulations on a set of hierarchically coarsened meshes. Furthermore, the availability of underlying fine-mesh (subgrid) data corresponding to coarse-mesh locations allows a clearer assessment of the effects of unresolved fluctuations on apparent reactivity. A key to MRA-MSR is the correlation of eddy structures at coarser mesh levels, which is facilitated by the transfer of filtered fine-mesh velocity information. A seven-mesh MRA-MSR hierarchy using three resolution levels is applied to one of the Sydney bluff-body stabilized methane–hydrogen flames. Analysis of the simultaneously evolved data at different resolution levels reveals several interesting trends. First, at high Damköhler numbers, there is clear evidence of attenuation of apparent reactivity due to the effects of unresolved fluctuations. Secondly, single-point, single-time filtered density functions of a normalized subgrid Damköhler number show a characteristic beta probability density function (PDF) form and display evidence of scale similarity. Interrogation of the MRA-MSR database also shows that the recently-proposed least-squares minimization (LSM) turbulence–chemistry interaction model can account for the observed diminishment in reactivity at high Damköhler numbers but cannot reduce scatter significantly. A new form of the LSM model, which makes use of the normalized subgrid Damköhler number beta PDF distribution, performs slightly better than the original model, illustrating the potential of MRA-MSR both in assessing existing closure concepts and in developing new ones.

42 ENGINEERING↗
NASA NTRS1985

Calculation of confined swirling flows

Swirling jets are an important constituent flow of many types of combustion equipment. Flame size, shape, stability and combustion intensity are all favorably influenced by swirling some portion of the incoming air with the fuel. Fluid dynamics computer codes are being developed for the mathematical simulation of flows in practical combustors. It is important that these codes have the ability to accurately and reliably calculate swirling flows. Turbulent flow calculations with swirl can be inaccurate due to at least three major reasons: limitations of the turbulence modeling used, incorrect or inadequate specifications of inlet boundary conditions, and, error introduced through the numerics. These sources of error are described, and examples provided of each. The state-of-the-art in such calculations is reviewed. The controlling character of current numerics is demonstrated, and it is concluded that as a first priority, better numerics must be arrived at, and that an improved understanding of the discharge flow from swirl generators is essential as a second priority. When these priorities are satisfied an improved turbulence model is desirable. It is also concluded that an improved understanding of swirl-induced recirculation would be an asset.

Sturgess, G. J.↗
NASA NTRS1929

A Study of Wing Flutter

Part I describes vibration tests, in a wind tunnel, of simple airfoils and of the tail plane of an M0-1 airplane model; it also describes the air flow about this model. From these tests are drawn inferences as to the cause and cure of aerodynamic wing vibrations. Part II derives stability criteria for wing vibrations in pitch and roll, and gives design rules to obviate instability. Part III shows how to design spars to flex equally under a given wing loading and thereby economically minimize the twisting in pitch that permits cumulative flutter. Resonant flutter is not likely to ensue from turbulence of air flow along past wings and tail planes in usual flying conditions. To be flutterproof a wing must be void of reversible autorotation and not have its centroid far aft of its pitching axis, i. e., axis of pitching motion. Danger of flutter is minimized by so proportioning the wing's torsional resisting moment to the air pitching moment at high-speed angles that the torsional flexure is always small. (author)

Zahm, A F↗
NASA NTRS2004

Laser Sounder for Global Measurement of CO2 Concentrations in the Troposphere from Space: Progress

We describe progress toward developing a laser-based technique for the remote measurement of the tropospheric CO2 concentrations from orbit. Our goal is to demonstrate a lidar technique and instrument technology that will permit measurements of the CO2 column abundance in the lower troposphere from aircraft at the few ppm level, with a capability of scaling to permit global CO2 measurements from orbit. Accurate measurements of the tropospheric CO2 mixing ratio from space are challenging due to the many potential error sources. These include possible interference from other trace gas species, the effects of temperature, clouds, aerosols & turbulence in the path, changes in surface reflectivity, and variability in dry air density caused by changes in atmospheric pressure, water vapor and topographic height. Some potential instrumental errors include frequency drifts in the transmitter, small transmission and sensitivity drifts in the instrument. High signal-to-noise ratios and measurement stability are needed for mixing ratio estimates at the few ppm level. We have been developing a laser sounder approach as a candidate for a future space mission. It utilizes multiple different laser transmitters to permit simultaneous measurement of CO2 and O2 extinction, and aerosol backscatter in the same measurement path. It directs the narrow co-aligned laser beams from the instrument's fiber lasers toward nadir, and measures the energy of the strong laser echoes reflected from the Earth's land and water surfaces. During the measurement its narrow linewidth lasers are rapidly tuned on- and off- selected CO2 line near 1572 nm and an O2 absorption line near 770 nm. The receiver measures the energies of the laser echoes from the surface and any clouds and aerosols in the path with photon counting detectors. Ratioing the on- to off-line echo pulse energies for each gas permits the column extinction and column densities of CO2 and O2 to be estimated simultaneously via the differential absorption lidar technique. For the on-line wavelengths, the side of the selected absorption lines are used, which due to pressure broadening, weights the measurements to the lower troposphere, where CO2 variations caused by surface sources and sinks are largest. Simultaneous measurements of O2 column abundance are made using an identical approach using an O2 line. The laser backscatter profiles from clouds and aerosols are measured with other lidar channels, which permits identifying measurements influenced by clouds and/or aerosol scattering in the path. For space use, our lidar would continuously measure at nadir in near polar circular orbit. Using dawn and dusk measurements made over the same region will make it possible to sample the diurnal variations in CO2 mixing ratios. A 1-m diameter telescope is used for the receiver for all wavelengths. When averaging over 50 seconds, our calculations show a SNR of approximately 1500 is achievable for each gas at each on- and off-line measurement. Measurements from such a mission can be used to generate monthly global maps of the lower tropospheric CO2 column abundance. Our calculations show global coverage with an accuracy of a few ppm with a spatial resolution of approximately 50,000 sq. km are achievable each month. We have demonstrated some key elements of the laser, detector and receiver approaches in the laboratory and with measurements over a 206 m horizontal path. These include stable measurements of CO2 line shapes in an absorption cell using a fiber laser amplifier seeded by a tunable diode laser, measurement of small amplitude changes at low optical signal levels with the PMT receiver, and comparison of the horizontal path measurements of CO2 against those from an in-situ instrument.

Abshire, J. B.↗