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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
NASA NTRS2022

Linear Disturbance Amplification Over Blunted Flat Plates in High-Speed Flows

Modal and nonmodal instability characteristics of cylindrically blunted flat plates with varying leading edge radii are described for Mach 4 and Mach 6 freestream conditions. The selection of leading edge radii and freestream parameters is informed by experimental conditions. The investigation of this 2D problem provides a slow entropy layer swallowing which allows for an isolated development of perturbations seeded upstream within different wall-normal regions of the flow. At both Mach numbers, a decrease in modal instability amplification was seen as the leading edge radius was increased. Nonmodal analysis reveals amplifying perturbations in the boundary layer as well as in the entropy layer. The medium bluntness regime exhibits the strongest amplification of nonmodal disturbances that is nonmonotonic in character. Small amplitude boundary forcing at the plate surface or volumetric forcing at various wall-normal heights was used to account for receptivity effects. While wall forcing effectively induced modal instabilities, only an actuation above the boundary layer captured disturbances that amplify within the entropy layer. The optimal nonmodal theory’s entropy-layer disturbance evolution exhibited outstanding agreement with controlled forcing, including receptivity effects. The evolution of entropy-layer disturbances from the optimal nonmodal theory showed excellent agreement with the results of receptivity to controlled forcing. Therefore, the nonmodal optimal growth analysis may provide a useful as well as efficient technique to identify the complete disturbance spectrum in blunt hypersonic configurations, where both modal and nonmodal disturbances can amplify in the boundary-layer and entropy-layer regions.

Boundary layer transition↗
NASA NTRS1990

A review of instability and noise propagation in supersonic flows

Originally analytical and numerical models were to be developed for noise production in supersonic jets, wakes and free shear layers. While the effort was concentrated initially on these aspects, other topics were also pursued, most were of interest to the Jet Noise Group of the Aeroacoustics Branch. An overview is given of subjects reviewed and the investigations that were carried out. A significant effort was devoted to numerically predicting the flow field of a turbulent supersonic wall jet. This information is necessary for computing the pressure in the far field. The wall jet was selected because it represents a generic flow that can be associated with plug nozzle in supersonic engines. It combines the characteristic of a boundary layer with that of a free shear flow. The spatially evolving flow obtained using Dash's code would form the input for the stability analysis program. This analysis would determine the large scale instability wave within the flow. The far field pressure can be computed from the shape of the evolving large scale structure by asymptotic methods. Flow characteristics obtained from a program that analyses the turbulent downstream supersonic flow in a nozzle are described and compared with experimental results.

Daudpota, Q. Isa↗
NASA NTRS1979

Accelerated development and flight evaluation of active controls concepts for subsonic transport aircraft

The flight test of an active load alleviation/extended span for the L-1011 wide-body transport aircraft, and piloted simulation work leading to use of active stability augmentation with a small tail and aft center of gravity are reported. The extended span showed the expected cruise drag reduction of 3%. The small tail is expected to reduce cruise drag by another 3%, and eventual use of more aft center of gravity with active stability augmentation will provide further fuel savings. The active load alleviation functions included maneuver load control (MLC) and elastic mode suppression (EMS), using symmetric motions of the outboard ailerons to reduce wing bending loads in maneuvers or long-term up- or down-drafts (MLC), and to damp wing bending motions in turbulence (EMS). A gust load alleviation function using the active horizontal tail to provide airplane pitch damping in turbulence was found unnecessary. The piloted simulation tests evaluated criteria for augmentation-on and augmentation-off flying qualities. of a simple pitch control law was verified at neutral static margin. The simulation tasks established the basis for follow-on construction and flight testing of a small tail with active stability augmentation.

Source record↗
NASA NTRS1997

Studies of Premixed Laminar and Turbulent Flames at Microgravity

Several topics related to the experimental and theoretical study of the combustion limits in premixed flames at reduced gravity have been investigated. The topics studied include: (1) numerical simulation of flame ball structure and stability; (2) interferometry for micro-g combustion experiments; (3) simulation of buoyancy effects in premixed flames using aqueous autocatalytic reactions; (4) premixed flame propagation in Hele-Shaw cells; and (5) experimental characterization of flame balls.

Ronney, Paul D.↗
NASA NTRS2024

Intrastep, Stage-Value Predictors for Diagonally-Implicit Runge-Kutta Methods

To better identify the necessary attributes of good stage-value predictors (SVPs), numerous SVPs are designed for an existing: ESDIRK4(3)7L[2]SA [26] and a new: ESDIRK4(3)8L[2]SA scheme.1 Both are stifflyaccurate, stage-order two, explicit, singly-diagonally implicit Runge–Kutta (ESDIRK) schemes. Tradeoffs are studied in the parameter spaces enforcing the constraints on accuracy, linear stability, nonlinear stability and coefficient size to determine which objectives correlate with effective predictors. The SVPs are tested on three challenging external aerodynamics problems [107 − 108 degrees of freedom (DoFs)], each with a different level of stiffness. The problems include two 3D airfoils simulations and one canonical turbulence simulation. All simulations use the compressible Navier-Stokes equations (CNSE). An entropy stable spectral collocation formulation is used for discretizing the spatial terms in the equations. Simulations are performed at a wide variety of temporal error tolerances. Problems that are sufficiently stiff (e.g., lax temporal error tolerances) benefit from SVPs designed with second-order accuracy and stability properties: A-stability, and L-stability, rather than high accuracy constraints. Simulations with modest stiffness (e.g., strict error tolerances) are better suited for SVPs designed using high accuracy constraints. Designing SVPs with enhanced stability properties is tedious but worthwhile. Simulation times are reduced with optimal SVPs by as much as 100% on some stages, with combined stepwise improvements of between 50 − 100% for both methods. A comparative study is performed with the two aforementioned methods as well as four other ESDIRKs. The newly designed ESDIRK4(3)8L[2]SA with γ ≈ 1/10, proves to be the most efficient of the six tested ESDIRK schemes simulating the CNSE.

Diagonally-Implicit Runge-Kutta↗
NASA NTRS1986

An explicit finite-volume time-marching procedure for turbulent flow calculations

A method was developed which calculates two-dimensional, transonic, viscous flow in ducts. The finite-volume, time-marching formulation is used to obtain steady flow solutions of the Reynolds-averaged form of the Navier-Stokes equations. The entire calculation is performed in the physical domain. Control volumes are chosen so that smoothing of flow properties, typically required for stability, is not required. Different time steps are used in the different governing equations. A new pressure interpolation scheme is introduced which improves the shock capturing ability of the method. A multi-volume method for pressure changes in the boundary layer allows calculations which use very long and thin control volumes (length/height - 1000). The method is compared with two test cases. Essentially incompressible turbulent boundary layer flow in an adverse pressure gradient is calculated and the computed distributions of mean velocity and shear are in good agreement with the measurements. Transonic viscous flow in a converging diverging nozzle is calculated; the Mach number upstream of the shock is approximately 1.25. The agreement between the calculated and measured shock strength and total pressure losses is good.

Nicholson, Stephen↗
NASA NTRS1987

An explicit finite-volume time-marching procedure for turbulent flow calculations

A method was developed which calculates two-dimensional, transonic, viscous flow in ducts. The finite-volume, time-marching formulation is used to obtain steady flow solutions of the Reynolds-averaged form of the Navier-Stokes equations. The entire calculation is performed in the physical domain. Control volumes are chosen so that smoothing of flow properties, typically required for stability, is not required. Different time steps are used in the different governing equations. A new pressure interpolation scheme is introduced which improves the shock capturing ability of the method. A multi-volume method for pressure changes in the boundary layer allows calculations which use very long and thin control volumes (length/height - 1000). The method is compared with two test cases. Essentially incompressible turbulent boundary layer flow in an adverse pressure gradient is calculated and the computed distributions of mean velocity and shear are in good agreement with the measurements. Transonic viscous flow in a converging diverging nozzle is calculated; the Mach number upstream of the shock is approximately 1.25. The agreement between the calculated and measured shock strength and total pressure losses is good.

Nicholson, Stephen↗
DOE OSTI2026

Impact of impurities on peeling–ballooning modes and turbulence in tokamak plasmas

This study investigates the impact of various impurity species on peeling–ballooning (PB) modes and microturbulence in tokamak plasmas through the extension of traditional two-fluid and gyro-landau-fluid (GLF) models. By incorporating finite Larmor radius (FLR) effects, the analysis provides a comprehensive understanding of impurity-driven impact and its interaction with plasma turbulence. Depending on charge state and local plasma conditions, heavy impurities may exhibit gyro-radii larger than those of main ions, which are captured in the extended GLF model presented. Following the presentation of modified two-fluid equations incorporating impurity effects, we systematically analyze the distinctions between impurity and main ion dynamics and their resultant feedback mechanisms on plasma behavior. Derivation of the linear dispersion relation enables quantification of impurity-mediated modifications to: plasma vorticity, diamagnetic drift and gyroviscous effects, electron Hall physics, and FLR effects. BOUT++ – based linear simulations corroborate this formalism, demonstrating systematic stabilization of PB modes upon impurity seeding. And then operational implications for practical impurity control strategies in tokamak devices are proposed. The results underscore the necessity of impurity management to maintain stability and optimize plasma confinement, with specific focus on how FLR effects contribute to transport dynamics. This work paves the way for enhanced modeling and simulation efforts, supporting the development of strategies to control impurity-induced turbulence and improve overall reactor performance.

BOUT++ simulation↗
DOE OSTI2022

A direct numerical simulation study of the dilution tolerance of propane combustion under spark-ignition engine conditions

Modern spark ignition internal combustion (IC) engines rely on highly diluted fuel-air mixtures to achieve high brake thermal efficiencies. To support this, new engine designs have introduced high stroke-to-bore ratios and cylinder head designs that promote high tumble flow and turbulence intensities. However, mixture dilution through exhaust gas recirculation (EGR) is limited by combustion instabilities manifested in the form of cycle-to-cycle variability. Propane has been observed to have superior EGR dilution tolerance than gasoline, which makes it a very competitive low-carbon fuel for the new IC engines without sacrificing efficiency. Two-dimensional direct numerical simulations (DNS) are performed with detailed chemistry to study and contrast the effect of turbulence intensity and dilution on propane and iso-octane premixed flames at high pressure conditions similar to those in-cylinder. A new reduced mechanism for propane consisting of 53 transported species and 17 quasi-steady state species is developed based on a previously published mechanism and used in these simulations. Three levels of turbulence intensity and two levels of exhaust gas dilution are chosen based on conditions relevant to IC engine operation. The DNS results are analyzed based on the evolution of the flame surface area and the statistics of its driving terms, which are found to be similar for both fuels when there is no dilution but considerably different under high dilution. The analysis of the DNS data provides fundamental insights into the underlying mechanisms for improved stability under dilution.

42 ENGINEERING↗
NASA NTRS2014

Conceptual Design Optimization of an Augmented Stability Aircraft Incorporating Dynamic Response Performance Constraints

This research focused on incorporating stability and control into a multidisciplinary de- sign optimization on a Boeing 737-class advanced concept called the D8.2b. A new method of evaluating the aircraft handling performance using quantitative evaluation of the sys- tem to disturbances, including perturbations, continuous turbulence, and discrete gusts, is presented. A multidisciplinary design optimization was performed using the D8.2b transport air- craft concept. The con guration was optimized for minimum fuel burn using a design range of 3,000 nautical miles. Optimization cases were run using xed tail volume coecients, static trim constraints, and static trim and dynamic response constraints. A Cessna 182T model was used to test the various dynamic analysis components, ensuring the analysis was behaving as expected. Results of the optimizations show that including stability and con- trol in the design process drastically alters the optimal design, indicating that stability and control should be included in conceptual design to avoid system level penalties later in the design process.

Welstead, Jason↗
DOE OSTI2023

Gaussian process hydrodynamics

Abstract We present a Gaussian process (GP) approach, called Gaussian process hydrodynamics (GPH) for approximating the solution to the Euler and Navier-Stokes (NS) equations. Similar to smoothed particle hydrodynamics (SPH), GPH is a Lagrangian particle-based approach that involves the tracking of a finite number of particles transported by a flow. However, these particles do not represent mollified particles of matter but carry discrete/partial information about the continuous flow. Closure is achieved by placing a divergence-free GP prior ξ on the velocity field and conditioning it on the vorticity at the particle locations. Known physics (e.g., the Richardson cascade and velocity increment power laws) is incorporated into the GP prior by using physics-informed additive kernels. This is equivalent to expressing ξ as a sum of independent GPs ξ l , which we call modes, acting at different scales (each mode ξ l self-activates to represent the formation of eddies at the corresponding scales). This approach enables a quantitative analysis of the Richardson cascade through the analysis of the activation of these modes, and enables us to analyze coarse-grain turbulence statistically rather than deterministically. Because GPH is formulated by using the vorticity equations, it does not require solving a pressure equation. By enforcing incompressibility and fluid-structure boundary conditions through the selection of a kernel, GPH requires significantly fewer particles than SPH. Because GPH has a natural probabilistic interpretation, the numerical results come with uncertainty estimates, enabling their incorporation into an uncertainty quantification (UQ) pipeline and adding/removing particles (quanta of information) in an adapted manner. The proposed approach is suitable for analysis because it inherits the complexity of state-of-the-art solvers for dense kernel matrices and results in a natural definition of turbulence as information loss. Numerical experiments support the importance of selecting physics-informed kernels and illustrate the major impact of such kernels on the accuracy and stability. Because the proposed approach uses a Bayesian interpretation, it naturally enables data assimilation and predictions and estimations by mixing simulation data and experimental data.

Mathematics↗
DOE OSTI2023

Simulation of radioactive plume transport in the atmosphere including dynamics of particle aggregation and breakup

Accurate prediction of the atmospheric transport of debris particles relies heavily on our knowledge of the size distribution of the particles within a debris cloud. Assuming a fixed particle size during simulations is not always viable since the size distribution of the debris can change during transport. Various microphysical processes, such as aggregation and breakup, influence debris particles and dictate any changes to the size distribution. To track those changes that can occur, a population balance model can be adopted and instituted within a model framework. Nonetheless, many of the models that simulate the transport of radioactive debris following a device-driven fission incident have historically neglected to consider these processes. As such, this work describes our effort to develop a modeling framework capable of simulating the transport and deposition of a radioactive plume generated from a fission incident with a dynamic population balance including particle aggregation and breakup. Herein the impact of aggregation and breakup, individually and collectively, on the particle size distribution is explored using the developed framework. When simulating aggregation, for example, six mechanisms, including Brownian coagulation, the convective enhancement to Brownian coagulation, van der Waals-viscous force correction for Brownian coagulation, gravitational collection, turbulent inertial motion, and turbulent shear, are considered. Brownian coagulation and its corrections have, as one would expect, a large impact on relatively small aggregates. Aggregates with a diameter that is less than or equal to 1.0 μm, for instance, comprise 50.6 vol % of all aggregates in the absence of aggregation and 31.2 vol % when Brownian coagulation and its corrections are considered. Gravitational collection and, to a much lesser extent, turbulent shear and turbulent inertial motion are, conversely, of great importance to relatively large aggregates (i.e., diameter greater than 3.0 μm). Additionally, the individual effects of atmospheric and particle parameters, such as wind speed and particle density, are examined. Of the parameters examined, turbulent energy dissipation and aggregate fractal dimension (i.e., aggregate shape with lower values representing more irregular particles) were of substantial importance since both terms directly impact aggregate stability and, by extension, the breakup rate. Large-scale transport and deposition simulations in a dry atmosphere are also presented and discussed as a proof of concept.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
NASA NTRS1982

ISEE-1 and -2 observations of laminar bow shocks - Velocity and thickness

The ISEE-1 and -2 spacecraft have provided the opportunity for studying the thickness and velocity of the bow shock on a regular basis. There are two models for the thickness of quasi-perpendicular laminar shocks. Galeev (1976) has calculated the thickness of the nearly perpendicular laminar bow shock by means of a weak ion sound turbulence (WIST) approach. Morse and Greenstadt (1976), on the other hand, assume that the shock thickness adjust itself to stay close to the drift velocity for marginal stability of current driven electrostatic waves. The present investigation is concerned with an evaluation of the competing theories, taking into account ISEE-1 and -2 observations of laminar bow shocks. It is found that the formula used by Morse and Greenstadt does not predict the observed thicknesses as well as the formula given by Galeev.

Russell, C. T.↗
NASA NTRS1990

Reynolds stress closure in jet flows using wave models

A collection of papers is presented. The outline of this report is as follows. Chapter three contains a description of a weakly nonlinear turbulence model that was developed. An essential part of the application of such a closure scheme to general geometry jets is the solution of the local hydrodynamic stability equation for a given jet cross-section. Chapter four describes the conformal mapping schemes used to map such geometries onto a simple computational domain. Chapter five describes a solution of a stability problem for circular, elliptic, and rectangular geometries. In chapter six linear models for the shock shell structure in non-circular jets is given. The appendices contain reprints of papers also published during this study including the following topics: (1) instability of elliptic jets; (2) a technique for predicting the shock cell structure in non-circular jets using a vortex sheet model; and (3) the resonant interaction between twin supersonic jets.

Morris, Philip J.↗
NASA NTRS1992

Boundary layer structure over the ocean observed by LEANDRE 1 during a tramontane event

A new airborne backscatter lidar, has been developed by CNRS (Service d'Aeronomie, (SA) Laboratoire de Meteorologie Dynamique (LMD) and the Institut des Sciences de 1'Univers) in the frame of the LEANDRE research program. It has been qualified on board the ARAT in autumn 1989 and spring 1990 and was involved in its first cooperative campaign during PYREX in October and November 1990. During this campaign, lidar observations of the perturbations induced on tropospheric flow and boundary layer structure were performed, and results are presented. A large number of experiments were performed, for synoptic situation description (meteorological radiosoundings, constant level balloons) and local flow analysis (aircrafts, radars, sodars). For the first time in such an experiment, a lidar has been flown on a research aircraft to perform altitude resolved observations of these perturbations, and we will present here results obtained for deflected flow structure. In the presence of a synoptic northerly flow, part of it is deflected to the east by the Pyrenees, and accelerated over the Mediterranean by the mountain surroundings. In this case, a low level wind is generated (the Tramontane) bringing cold and dry air over the Mediterranean Sea. As the sea is still at a warm temperature in November (around 17 degrees), an Internal Marine Boundary Layer rapidly grows over the first tens of kilometers and stabilizes at about 1 km depth, corresponding to an altitude just below the Lifting Condensation Level. The whole Marine Atmospheric Boundary Layer (MABL) is characterized by highly turbulent motions bringing large humid particles from the surface up to its top. The lidar signal due to scattering by these particles is then representative of the turbulent kinetic energy in this layer.

Flamant, C.↗
NASA NTRS1993

Effect of micron-sized roughness on transition in swept-wing flows

Boundary-layer transition-to-turbulence studies are conducted in the Arizona State University Unsteady Wind Tunnel on a 45-degree swept airfoil. The pressure gradient is designed so that the initial stability characteristics are purely crossflow-dominated. Flow visualization and hot-wire measurements show that the development of the crossflow vortices is influenced by roughness near the attachment-line. Comparisons of transition location are made between a painted surface, a machine-polished surface, and a hand-polished surface. Then, isolated 6 micron roughness elements are placed near the attachment line on the airfoil surface under conditions of the final polish (0.25 micron rms). These elements amplify a centered stationary crossflow vortex and its neighbors, resulting in localized early transition. The diameter, height, and location of these roughness elements are varied in a systematic manner. Spanwise hot-wire measurements are taken behind the roughness element to document the enhanced vortices. These scans are made at several different chord locations to examine vortex growth.

Radeztsky, Ronald H., Jr.↗