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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
NASA NTRS2009

Construction of Low Dissipative High Order Well-Balanced Filter Schemes for Non-Equilibrium Flows

The goal of this paper is to generalize the well-balanced approach for non-equilibrium flow studied by Wang et al. [26] to a class of low dissipative high order shock-capturing filter schemes and to explore more advantages of well-balanced schemes in reacting flows. The class of filter schemes developed by Yee et al. [30], Sjoegreen & Yee [24] and Yee & Sjoegreen [35] consist of two steps, a full time step of spatially high order non-dissipative base scheme and an adaptive nonlinear filter containing shock-capturing dissipation. A good property of the filter scheme is that the base scheme and the filter are stand alone modules in designing. Therefore, the idea of designing a well-balanced filter scheme is straightforward, i.e., choosing a well-balanced base scheme with a well-balanced filter (both with high order). A typical class of these schemes shown in this paper is the high order central difference schemes/predictor-corrector (PC) schemes with a high order well-balanced WENO filter. The new filter scheme with the well-balanced property will gather the features of both filter methods and well-balanced properties: it can preserve certain steady state solutions exactly; it is able to capture small perturbations, e.g., turbulence fluctuations; it adaptively controls numerical dissipation. Thus it shows high accuracy, efficiency and stability in shock/turbulence interactions. Numerical examples containing 1D and 2D smooth problems, 1D stationary contact discontinuity problem and 1D turbulence/shock interactions are included to verify the improved accuracy, in addition to the well-balanced behavior.

Wang, Wei↗
NASA NTRS2017

Model Validation for Propulsion - On the TFNS and LES Subgrid Models for a Bluff Body Stabilized Flame

With advances in computational power and availability of distributed computers, the use of even the most complex of turbulent chemical interaction models in combustors and coupled analysis of combustors and turbines is now possible and more and more affordable for realistic geometries. Recent more stringent emission standards have enticed the development of more fuel-efficient and low-emission combustion system for aircraft gas turbine applications. It is known that the NOx emissions tend to increase dramatically with increasing flame temperature. It is well known that the major difficulty, when modeling the turbulence-chemistry interaction, lies in the high non-linearity of the reaction rate expressed in terms of the temperature and species mass fractions. The transport filtered density function (FDF) model and the linear eddy model (LEM), which both use local instantaneous values of the temperature and mass fractions, have been shown to often provide more accurate results of turbulent combustion. In the present, the time-filtered Navier-Stokes (TFNS) approach capable of capturing unsteady flow structures important for turbulent mixing in the combustion chamber and two different subgrid models, LEM-like and EUPDF-like, capable of emulating the major processes occurring in the turbulence-chemistry interaction will be used to perform reacting flow simulations of a selected test case. The selected test case from the Volvo Validation Rig was documented by Sjunnesson.

bluff body↗
NASA NTRS2017

A Comparison of Three Random Number Generators for Aircraft Dynamic Modeling Applications

Three random number generators, which produce Gaussian white noise sequences, were compared to assess their suitability in aircraft dynamic modeling applications. The first generator considered was the MATLAB (registered) implementation of the Mersenne-Twister algorithm. The second generator was a website called Random.org, which processes atmospheric noise measured using radios to create the random numbers. The third generator was based on synthesis of the Fourier series, where the random number sequences are constructed from prescribed amplitude and phase spectra. A total of 200 sequences, each having 601 random numbers, for each generator were collected and analyzed in terms of the mean, variance, normality, autocorrelation, and power spectral density. These sequences were then applied to two problems in aircraft dynamic modeling, namely estimating stability and control derivatives from simulated onboard sensor data, and simulating flight in atmospheric turbulence. In general, each random number generator had good performance and is well-suited for aircraft dynamic modeling applications. Specific strengths and weaknesses of each generator are discussed. For Monte Carlo simulation, the Fourier synthesis method is recommended because it most accurately and consistently approximated Gaussian white noise and can be implemented with reasonable computational effort.

Grauer, Jared A.↗
NASA NTRS1962

STABILITY OF FILM BOILING TWO-PHASE FLOW IN CRYOGENIC SYSTEMS

The stability of film boiling two phase flow is investigated referring in particular to saturated cryogenic liquids. This type of fluid motion is unstable, such that a range of wave numbers will always be present in which disturbances are amplified. Further the mode of operation of the energy supply system may favour temperature fluctuations of large amplitude if the circuit does not provide feed-back compensation. The theoretically predicted stabilization due to capillarity is reduced considerably in liquefied gases with low surface tension, such as hydrogen. This effect causes the transition to turbulent flow to occur at low Reynolds numbers, and increases the heat transport in the turbulent regime. A comparison with experimental data, (horizontal cylinder), shows that the "critical" film Reynolds number ūδ/𝜈, at which instabilities become effective, is lowered from approximately 200, (O2), to about 30, (H2). Similar results have been found with other configurations. Therefore application of laminar film boiling solutions is restricted to small objects, entrance regions, or stagnation points. Turbulent correlations given so far are valid only for a certain group of substances with similar behavior.

LIQUID FLOW↗
NASA NTRS1952

The excitation of unstable perturbations in a laminar friction layer

With the aid of the method of small oscillations which was used successfully in the investigation of the stability of laminar velocity distributions in the presence of two-dimensional perturbations, the excitation of the unstable perturbations for the Hartree velocity distributions occurring in plane boundary-layer flow for decreasing and increasing pressure is calculated as a supplement to a former report. The results of this investigation are to make a contribution toward calculation of the transition point on cylindrical bodies.

BOUNDARY LAYER TRANSITION - EFFECT OF SURFACE ROUG↗
DOE OSTI2022

Impact of pedestal density gradient and collisionality on ELM dynamics

BOUT++ turbulence simulations are conducted to capture the underlying physics of small ELM characteristics achieved by increasing separatrix density via controlling strike points from vertical to horizontal divertor plates for three EAST discharges. BOUT++ linear simulations show that the most unstable modes change from high-n ideal ballooning modes to intermediate-n peeling–ballooning modes and eventually to peeling–ballooning stable plasmas in the pedestal. Nonlinear simulations show that the fluctuation is saturated at a high level for the lowest separatrix density. The ELM size decreases with increasing separatrix density, until the fraction of this energy lost during the ELM crash becomes less than 1% of the pedestal stored energy, leading to small ELMs. Simulations indicate that small ELMs can be triggered either by the marginally peeling–ballooning instability near the peak pressure gradient position inside the pedestal or by a local instability in the pedestal foot with a larger separatrix density gradient. The pedestal collisionality scan for type-I ELMs with steep pedestal density gradient shows that both linear growth rate and ELM size decrease with increasing collisionality. On the contrary, the pedestal collisionality and pedestal density width scan with a weak pedestal density gradient indicate small ELMs can either be triggered by a high-n ballooning mode or by a low-n peeling mode in a low collisionality region 0.04–0.1. We report the simulations indicate the weaker the linear unstable modes near marginal stability with small linear growth rate, the lower nonlinearly saturated fluctuation intensity and the smaller turbulence spreading from the linear unstable zone to stable zone in the nonlinear saturation phase, leading to small ELMs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Deep Learning for Subgrid‐Scale Turbulence Modeling in Large‐Eddy Simulations of the Convective Atmospheric Boundary Layer

Abstract In large‐eddy simulations, subgrid‐scale (SGS) processes are parameterized as a function of filtered grid‐scale variables. First‐order, algebraic SGS models are based on the eddy‐viscosity assumption, which does not always hold for turbulence. Here we apply supervised deep neural networks (DNNs) to learn SGS stresses from a set of neighboring coarse‐grained velocity from direct numerical simulations of the convective boundary layer at friction Reynolds numbers Re τ up to 1243 without invoking the eddy‐viscosity assumption. The DNN model was found to produce higher correlation between SGS stresses compared to the Smagorinsky model and the Smagorinsky‐Bardina mixed model in the surface and mixed layers and can be applied to different grid resolutions and various stability conditions ranging from near neutral to very unstable. The DNN model can capture key statistics of turbulence in a posteriori (online) tests when applied to large‐eddy simulations of the atmospheric boundary layer.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1985

Aerodynamic properties of turbulent combustion fields

Flow fields involving turbulent flames in premixed gases under a variety of conditions are modeled by the use of a numerical technique based on the random vortex method to solve the Navier-Stokes equations and a flame propagation algorithm to trace the motion of the front and implement the Huygens principle, both due to Chorin. A successive over-relaxation hybrid method is applied to solve the Euler equation for flows in an arbitrarily shaped domain. The method of images, conformal transformation, and the integral-equation technique are also used to treat flows in special cases, according to their particular requirements. Salient features of turbulent flame propagation in premixed gases are interpreted by relating them to the aerodynamic properties of the flow field. Included among them is the well-known cellular structure of flames stabilized by bluff bodies, as well as the formation of the characteristic tulip shape of flames propagating in ducts. In its rudimentary form, the mechanism of propagation of a turbulent flame is shown to consist of: (1) rotary motion of eddies at the flame front, (2) self-advancement of the front at an appropriate normal burning speed, and (3) dynamic effects of expansion due to exothermicity of the combustion reaction. An idealized model is used to illustrate these fundamental mechanisms and to investigate basic aerodynamic features of flames in premixed gases. The case of a confined flame stabilized behind a rearward-facing step is given particular care and attention. Solutions are shown to be in satisfactory agreement with experimental results, especially with respect to global properties such as the average velocity profiles and reattachment length.

Hsiao, C. C.↗
NASA NTRS2019

Boundary Layer Stability and Transition in a Chemically Reacting Martian Atmosphere Using LASTRAC

The Thermal Protection System (TPS) is a significant portion of the mass of reentry vehicles, planetary probes, and Martian entry vehicles. Reducing this mass has benefits in terms of decreased fuel requirements and increased payload; however, due to the high risk and uncertainty, the TPS or heat shields are designed conservatively by assuming fully turbulent flow. Laminar flow results in reduced heat flux, and improved transition prediction has the potential to reduce TPS mass and uncertainty in aerothermodynamic predictions. Limited previous research exists examining the problem of transition prediction for Martian atmospheric entry, with studies available on transition on the Mars Science Laboratory TPS. Although transition was demonstrated in wind-tunnel tests, uncertainty in the transition location resulted in a TPS designed for fully turbulent flow, and therefore greater mass than required for a partially laminar condition. In this work, we extend the boundary layer stability code LASTRAC, recently modified to include chemical and thermal nonequilibrium capabilities, to a model of the Martian atmosphere. LASTRAC provides Parabolized Stability Equations (PSE) as well as Linear Stability Theory (LST) to predict the stability of a boundary layer and transition with semi-empirical e(sup N) methods. Results included in this work compare disturbance growth characteristics between air and Martian atmosphere at similar non-dimensional freestream conditions on a simple flat plate geometry. Both chemical nonequilibrium and thermochemical nonequilibrium, as well as both PSE and LST, are used.

Kline, H. L.↗
NASA NTRS2019

Adjoint-Based Anisotropic Mesh Adaptation for a Stabilized Finite-Element Flow Solver

An adjoint solver is implemented in the FUN3D stabilized finite-element flow solver. The adjoint solution is used to generate anisotropic, adapted meshes to control error in scalar out- put functionals, such as lift or drag coefficient. The inviscid and turbulent adjoints are verified with a finite-difference derivative approximation and can be used in design optimization in addition to mesh adaptation. The adjoint capability represents an extension of previous research using the FUN3D stabilized finite-element capability with metric-based mesh adaptation and interpolation-based error estimates to generate highly anisotropic adapted meshes for turbulent flows. In the present work, a metric-based approach is again utilized, where the adjoint and the primal solutions both contribute to the generation of a metric tensor field that is subsequently used to produce the required anisotropic mesh for each adaptation cycle. Adaptive results are then shown for an inviscid supersonic flow over a diamond airfoil, inviscid transonic flow over the ONERA M6 wing, and viscous laminar flow over the NACA 0012 airfoil, all using drag force as the output functional. Adjoint-based adaptation is compared with a multiscale solution-based approach that controls the L norm of Mach number interpolation error to demonstrate the effectiveness and effciency of the adjoint-based adaptive mesh technology.

Balan, Aravind↗
NASA NTRS1995

On treatments of fetch and stability sensitivity in large-area estimates of sensible heat flux over sea ice

Bulk transfer coefficients estimated as a function of atmospheric stability and sea-ice lead width are combined with Arctic meteorological observations and ice thickness data to calculate the sensitivity of turbulent flux estimates to changes in lead width, wind speed, air temperature, and ice thickness for a high-concentration ice pack. These results are considered in terms of bulk transfer parameterizations that use a fixed transfer coefficient r that address atmospheric stability only. On the basis of the fetch-sensitive parameterizations considered here, differences in lead width for widths up to about 200 m can exert a substantial influence on sensible heat transfer coefficients and heat flux from leads under typical Arctic conditions. Fluxes from an open water lead decrease by 34% if fetch increases from 10 m to 100 m. This effect is greatest for open water leads, decreases considerably as leads refreeze, and is negligible for ice thicker than about 0.3 m. If open or newly refrozen leads make up 2% of the ice cover, than an increase in mean fetch from 10 m to 100 m yields a decrease of about 2 W/sq m in areally averaged flux from the ice pack. Calculations using observed and theoretical lead width distributions suggest that parameterizing lead widths in a sea ice model can be done effectively using a single, representative lead width rather than requiring a full distribution of widths. When coupled to the lowest atmospheric boundary layer using a bulk similarity theory model, this sensitivity of heat transfer to fetch results in substantially higher near-surface air temperatures over narrow leads, with equilibrium air temperatures decreasing by about 50% as fetch increases from 10 to 100 m.

Maslanik, J. A.↗
NASA NTRS1994

The effects of profiles on supersonic jet noise

The effect of velocity profiles on supersonic jet noise are studied by using stability calculations made for a shock-free coannular jet, with both the inner and outer flows supersonic. The Mach wave emission process is modeled as the noise generated by the large scale turbulent structures or the instability waves in the mixing region. Both the vortex-sheet and the realistic finite thickness shear layer models are considered. The stability calculations were performed for both inverted and normal velocity profiles. Comparisons are made with the results for an equivalent single jet, based on equal thrust, mass flow rate and exit area to that of the coannular jet. The advantages and disadvantages of these velocity profiles as far as noise radiation is concerned are discussed. It is shown that the Rayleigh's model prediction of the merits and demerits of different velocity profiles are in good agreement with the experimental data.

Tiwari, S. N.↗
NASA NTRS1987

On a finite element CFD algorithm for compressible, viscous and turbulent aerodynamic flows

This paper develops and analyses individual construction aspects of an efficient and accurate finite element algorithm for prediction of viscous and turbulent flow fields of impact in aerodynamics. The theoretical construction employs a Taylor weak statement (TWS) for coincident embedding of stability mechanisms within a classic Galerkin finite element formulation of semidiscrete approximation error orthogonalization. A wide variety of the stabilizing mechanisms of independently derived CFD algorithms are contained within the TWS theory. An implicit construction that meets the requirement of efficient convergence to steady state is developed. The theoretical asymptotic error estimates of the TWS finite element algorithm for supersonic and viscous boundary layer flows are verified. Application to a three-dimensional turbulent flow is cited.

Baker, A. J.↗
NASA NTRS1978

Nonlinear stability of solar type 3 radio bursts. 2: Application to observations near 1 AU

A set of rate equations including strong turbulence effects and anomalous resitivity are solved using parmeters which model several solar type 3 bursts. Exciter distributions observed at 1 AU are excitation of the linear bump-in-tail instability, amplifying Langmuir waves above the threshold for the oscillating two stream instability (OTSI). The OTSI, and the attendant anomalous resistivity produce a rapid spectral transfer of Langmuir waves to short wavelengths, out of resonance with the electron exciter. Further energy loss of the beam is thus precluded. The various parameters needed to model the bursts are extrapolated inside 1 AU with similar results. Again, the OTSI is excited and decouples the electron beam from the Langmuir radiation. Reabsorption of the Langmuir waves by the beam is shown to be unimportant in all cases, even at 0.1 AU. The theory provides a natural explanation for the observed realationship between radio flux, and the electron flux.

Goldstein, M. L.↗
NASA NTRS1979

Nonlinear stability of solar type III radio bursts. II - Application to observations near 1 AU

A set of rate equations including strong turbulence effects and anomalous resistivity are solved using parameters which model several solar type III bursts. Analysis of these bursts has led to quantitative comparisons between several of the observed phenomena and the theory. Through use of an analytic model for the time evolution of the energetic electron exciter, it is found that the exciter distributions observed at 1 AU are unstable to the excitation of the linear bump-in-tail instability, amplifying Langmuir waves above the threshold for the oscillating two-stream instability (OTSI). The OTSI and the attendant anomalous resistivity produce a rapid spectral transfer of Langmuir waves to short wavelengths, out of resonance with the electron exciter. In addition, the various parameters needed to model the bursts are extrapolated inside 1 AU with similar results. Finally, reabsorption of the Langmuir waves by the beam is shown to be unimportant in all cases, even at 0.1 AU.

Goldstein, M. L.↗
NASA NTRSDate not supplied

Shock-Wave/Boundary-Layer Interaction (SWBLI) Experiments in the Presence of Transition-to-Turbulence on a Flat Plate Model in the NASA LaRC 31-Inch Mach 10 Air Tunnel

Historically, there have been a limited number of studies involving Transitional Shock-Wave/Boundary-Layer Interactions (XSWBLI) in hypersonic flows, the majority of which have been performed at Mach Numbers of 7.5 and below. Additionally, there have been even fewer such experiments which incorporate non-intrusive and optical diagnostics. NASA Langley Research Center (NASA LaRC), in collaboration with The University of Texas at San Antonio (UTSA) and The University of Tennessee Space Institute (UTSI), is currently planning and preparing for a series of experimental aerodynamic tests using the NASA LaRC 31-inch Mach 10 Air Tunnel. The experiments will primarily focus on studying the dynamics of Shock-Wave/Boundary-Layer Interactions (SWBLI) in the presence of transition (XSWBLI) and turbulence (SWTBLI) on a large, flat plate model at a freestream Mach number of 10. Besides continued interest in understanding the flow at high Mach numbers, executing such tests at Mach 10 avoids some known concerns and will aid in solving new problems. XSWBLI are a highly unsteady phenomena and the generation of XSWBLI at lower Mach numbers can prove to be a significant challenge. This concern will be somewhat mitigated in the Mach 10 flow, as boundary layers become increasingly stabilized, but some difficulty in achieving transition is still expected. Modelling SWBLI in the presence of transitional and turbulent boundary layers has also proven to be difficult, so a high-Mach number experiment which applies non-intrusive and optical diagnostics will aid in solving a unique problem as well as advancing the understanding and characterization of an aerodynamic surface at Mach 10.The NASA LaRC 31-inch Wind Tunnel incorporates a test section having a cross-section of 31-inches × 31-inches and provides optical access to the test section via three (3) ultra-violet (UV) transmitting windows on the top, side, and bottom. The Mach 10 operating conditions will consist of pressures (P0) ranging from 2.4 to 10MPa (348 to 1,450 psi) at a temperature (T0) of 1,000 K (1,800 °R). The NASA LaRC 31-inch wind tunnel will deliver 1-minute blow-down runs, of which 30-45 seconds will be consumed to reach the aforementioned test conditions. The pump down time between runs will be approximately 45-60 minutes, which will provide for 6-10 test runs per day over the course of a weeklong entry. The proposed flat plate model will be constructed of stainless steel with geometry consisting of a10-inch × 30-inch upper surface and a 2-inch thickness. The model will be designed to support multiple leading-edge inserts in order to examine blunt, round, and sharp leading edges at Mach 10. Adjustments will be made to manage boundary layer thickness as well as the strength and size of the shock interaction region for each test. Measurements and diagnostics will be performed though high-speed Schlieren, on-body high-speed pressure transducers, IR thermography, and oil flow visualization. Hypersonic vehicle applications have and will continue to emerge at the forefront of aerospace. The NASA-UTSA-UTSI team realizes the value of characterizing SWBLI, XSWBLI, and SWTBLI behaviors well beyond single-digit hypersonic Mach numbers. As hypersonic vehicle speed capabilities continue to increase, measurement and diagnostic methods for Mach numbers of 10 and beyond will provide the foundation for a firm understanding of flow field behavior which will directly influence the advancement of technologies towards the design and manufacturing of high-speed aerodynamic surfaces, controls, thermal protection systems, acoustic treatments, and structural components. The use of non-intrusive and optical diagnostic methods in such experimentation is pivotal in developing the visualization and empirical data necessary to advance the aforementioned technology areas.

Shockwave/Boundary-Layer Interactions↗
NASA NTRSDate not supplied

Boundary Layer Stability and Transition in a Chemically Reacting Martian Atmosphere using LASTRAC

The Thermal Protection System (TPS) is a significant portion of the mass of reentry vehicles, planetary probes, and Martian entry vehicles. Reducing this mass has benefits in terms of decreased fuel requirements and increased payload; however, due to the high risk and un-certainty, the TPS or heat shields are designed conservatively by assuming fully turbulent flow. Laminar flow results in reduced heat flux, and improved transition prediction has the potential to reduce TPS mass and uncertainty in aerothermodynamic predictions. Limited previous research exists examining the problem of transition prediction for Martian atmospheric entry, with studies available on transition on the Mars Science Laboratory TPS. Although transition was demonstrated in wind-tunnel tests, uncertainty in the transition location resulted in a TPS designed for fully turbulent flow, and therefore greater mass than required for a partially-laminar condition. In this work, we extend the boundary layer stability code LASTRAC, recently modified to include chemical and thermal nonequilibrium capabilities, to a model of the Martian atmosphere. LASTRAC provides Parabolized Stability Equations (PSE) as well as Linear Stability Theory (LST) to predict the stability of a boundary layer and transition with semi-empirical eN methods. Results included in this work compare disturbance growth characteristics between air and Martian atmosphere at similar nondimensional freestream conditions on a simple flat plate geometry. Both chemical nonequilibrium and thermochemical nonequilibrium, as well as both PSE and LST, are used.

Heather L Kline↗