Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Boundary Layer Stability”

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 721 records · Page 40

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

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

Transition↗

The Flux-Differencing Discontinuous Galerkin Method Applied to an Idealized Fully Compressible Nonhydrostatic Dry Atmosphere

Dynamical cores used to study the circulation of the atmosphere employ various numerical methods ranging from finite-volume, spectral element, global spectral, and hybrid methods. In this work, we explore the use of Flux-Differencing Discontinuous Galerkin (FDDG) methods to simulate a fully compressible dry atmosphere at various resolutions. We show that the method offers a judicious compromise between high-order accuracy and stability for large-eddy simulations and simulations of the atmospheric general circulation. In particular, filters, divergence damping, diffusion, hyperdiffusion, or sponge-layers are not required to ensure stability; only the numerical dissipation naturally afforded by FDDG is necessary. We apply the method to the simulation of dry convection in an atmospheric boundary layer and in a global atmospheric dynamical core in the standard benchmark of Held and Suarez

54 ENVIRONMENTAL SCIENCES↗

Quantifying the Thermodynamic Impacts on the Atmospheric Boundary Layer due to the Sea Breeze in the Coastal Houston Region

The atmospheric boundary layer (ABL) is unique in coastal regions because of kinematic and thermodynamic influences from continental and marine environments. Sea-breeze (SB) circulations act to equilibrate the land–sea temperature gradient through advecting marine air onshore. The strength of the SB varies in terms of stability, temperature, and moisture advection and influences air quality and weather forecasts. The Tracking Aerosol Convection Interactions Experiment (TRACER) collected a wealth of data on coastal boundary layer evolution, including observations from uncrewed aerial systems (UASs). Vertical profiles of temperature, humidity, and winds were collected by the OU CopterSonde UAS from June to September in the coastal region of Houston. These profiles offer 5-m vertical resolution, on average, every 30 min through diurnal transitions, SB events, and nearby deep convection. During the campaign, CopterSonde observations were gathered through 17 SB events, six of which led to convection initiation. The UAS data can resolve the thermodynamic evolution and interactions between the SB and the preexisting convective boundary layer. Results show large variability across observed SBs and their impacts on temperature and moisture. The intensity of thermodynamic changes depends on the time of sea-breeze passage and influence from the Galveston Bay Breeze, a secondary marine circulation commonly observed in this region. In quantifying the spectrum of SB impacts, equivalent potential temperature θ e is used to contextualize its role in convection initiation and evolution. In conclusion, while all SBs tend to increase θ e from moisture advection, the rate and timing of the θ e rise can distinguish convective from nonconvective cases.

54 ENVIRONMENTAL SCIENCES↗

On the role of secondary instabilities in laminar-turbulent transition of 2D and 3D boundary layers

The accurate prediction of laminar/turbulent transition is one of the fundamental problems in engineering fluid mechanics. There is almost unanimous consent that such a transition criterion should come from stability theory. Linear primary stability theory describes the initial stage of transition, but falls short of predicting transition. Only in conjunction with empirical correlations, the widely used e(sup n) method is obtained, which, however, lacks a solid physical base. Three-dimensional secondary instabilities are known to play an important role in the transition process. However, no use has been made so far of secondary instabilities, instability interactions or wave resonances to define a 'transition location'. The paper summarizes new attempts to identify certain interaction and resonance phenomena within the laminar-turbulent transition regime in two and three-dimensional boundary layers which are associated with rapid structural and temporal changes of fluctuations beyond their exponential growths.

Dallman, Uwe↗

Secondary instability in rotating-disk flow

Primary instability of the 3D boundary layer on a rotating disk introduces periodic modulation of the mean flow in the form of stationary crossflow vortices. We study the stability of this modulated mean flow with respect to secondary disturbances. Both fundamental and subharmonic resonance cases are considered, and their corresponding results indicate that the growth rate and the frequency of the secondary instability are insensitive to the exact nature of the resonance condition. The threshold primary stationary crossflow vortex amplitude for secondary instability found in this 3D incompressible boundary layer is significantly larger than that for a 2D boundary layer which is subjected to Tollmien-Schlichting instability. The secondary instability results in a pair of travelling counter-rotating vortices, tilted up and oriented at an angle to the primary stationary crossflow vortices.

Balachandar, S.↗

Aerodynamic Design Optimization for Natural Laminar Flow Airfoils

Natural laminar flow technology is a passive laminar flow control (LFC) strategy that seeks to delay the onset of boundary-layer transition (BLT) through shape optimization to reduce the drag of the aerodynamic vehicle. Adjoint-based design optimization for LFC is proposed in an integrated multidisciplinary framework, which includes the computational fluid dynamics (CFD), geometry and grid deformation, and linear stability analysis (LSA) for transition prediction. In particular, the BLT location is predicted using the dual N-factor method that is based on a linear stability theory (LST) eigenvalue problem. The dual N-factor criterion accounts for the amplification of planar Tollmien-Schlichting (TS) and stationary crossflow (CF) boundary-layer instabilities to predict the transition location in three-dimensional boundary-layer flows. The adjoint-based shape optimization procedure is based on an iteratively coupled CFD and LSA methodology to converge the transition location and flow solutions, as well as to calculate the sensitivities of the aerodynamic metrics of interest with respect to the flow and shape design parameters. The RAE 2822 airfoil at 0 and 30 degrees yaw angles, an angle of attack of 0.72 degrees, and subsonic conditions (M∞ = 0.19, Rec = 5.6 × 106 ) are used as baseline configurations for design optimization. The angle of attack and the vertical displacement of free-form-deformation control points are used as design variables to reduce the drag coefficient while reaching a specified lift coefficient. The optimized unswept airfoil designs achieve a 30% drag reduction accompanied by a downstream shift of the transition locations over both suction and pressure sides of the airfoil. The initial design iterations for the swept case also show a favorable trend in the drag reduction with transition delay over both sides.

Transition↗

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↗

Toward Transition Modeling in a Hypersonic Boundary Layer at Flight Conditions

An accurate physics-based transition prediction method integrated with computational fluid dynamics (CFD) solvers is pursued for hypersonic boundary layer flows over slender hypersonic vehicles at flight conditions. The geometry and flow conditions are selected to match relevant trajectory locations from the ascent phase of the HIFiRE-1 flight experiment, namely, a 7-degree half-angle cone with 2.5 mm nose radius, freestream Mach numbers in the range of 3.8 – 5.5 and freestream unit Reynolds numbers in the range of 3.3 × 10(exp 6) – 21.4 × 10(exp 6) m(exp -1). Earlier research had shown that the onset of transition during the HIFiRE-1 flight experiment correlated with an amplification factor of N ≈ 13.5 for the planar Mack modes. However, to incorporate the N-factor correlations into a CFD code, we investigate surrogate models for disturbance amplification that avoid the direct computation of stability characteristics. A commonly used approach for low-speed flows is based on an a priori database of stability characteristics for locally similar profiles. However, the results presented in this paper demonstrate that the application of this approach to hypersonic boundary layers over blunt spherical nose-tip cones leads to large, unacceptable errors in the predictions of amplification factors, mainly due to its failure in accounting for the effects of the entropy layer on the boundary-layer profiles along the length of the model. We propose and demonstrate an alternate approach that employs the stability computations for a canonical set of blunt cone configurations to train a physics-informed convolutional neural network model that is shown to provide substantially improved transition predictions for hypersonic flow configurations with entropy-layer effects. Furthermore, the excellent performance of the neural network model is also confirmed for cone configurations with nose radius and half-angle values that do not correspond to those used to build the database. Finally, the convolutional neural network model is shown to outperform the linear stability calculations for underresolved basic states.

Pedro Paredes↗

Flow-induced morphological instabilities due to temporally-modulated stagnation-point flow

The influence of periodically-modulated planar stagnation-point flow on the morphological stability of a directionally-solidifying interface is presently considered with a view to the effect of unsteady nonparallel flows on single-crystal growth. The modeling of the system assumes that the viscous boundary layer thickness is much greater than that of the solute boundary layer, and that the modulation frequency is much smaller than the strength of plane stagnation-point flow. The solidifying interface is either stabilized or destabilized depending on the ratio of the period of modulation to the solute-diffusion time.

Merchant, G. J.↗

A lag model for turbulent boundary layers developing over rough bleed surfaces

Boundary layer mass removal (bleed) through spanwise bands of holes on a surface is used to prevent or control separation and to stabilize the normal shock in supersonic inlets. The addition of a transport equation lag relationship for eddy viscosity to the rough wall algebraic turbulence model of Cebeci and Chang was found to improve agreement between predicted and measured mean velocity distributions downstream of a bleed band. The model was demonstrated for a range of bleed configurations, bleed rates, and local free stream Mach numbers. In addition, the model was applied to the boundary layer development over acoustic lining materials for the inlets and nozzles of commercial aircraft. The model was found to yield accurate results for integral boundary layer properties unless there was a strong adverse pressure gradient.

Lee, J.↗

Linearization of the boundary-layer equations of the minimum time-to-climb problem

Ardema (1974) has formally linearized the two-point boundary value problem arising from a general optimal control problem, and has reviewed the known stability properties of such a linear system. In the present paper, Ardema's results are applied to the minimum time-to-climb problem. The linearized zeroth-order boundary layer equations of the problem are derived and solved.

Ardema, M. D.↗

Phase transitions and convection in icy satellites

The effects of solid-solid phase changes on subsolidus convection in the large icy moons of the outer solar system are considered. Phase transitions affect convection via processes that distort the phase change boundary and/or influence buoyancy through thermal expansion. Linear stability analyses are carried out for ice layers with a phase change at the midplane. Two exothermic phase transitions (ice I - ice II, ice VI - ice VIII) and two endothermic transitions (ice I - ice III, ice II - ice V) are considered. For the exothermic cases, the phase change can either impede or enhance whole-layer convection. For the endothermic cases, the phse change always inhibits whole-layer convective overturn and tends to enforce two-layer convection. These results play some constraints on possible models of icy satellite evolution and structure.

Bercovici, D.↗

Effect of surface catalysis on heating to ceramic coated thermal protection systems for transatmospheric vehicles

This paper describes the effect of surface catalysis on the heat transfer rate to the heat shield of a typical Transatmospheric Vehicle (TAV) during ascent and atmospheric entry. Surface kinetics and optical properties obtained from arc-jet tests on candidate thermal protection systems (coated metals) were used in a reacting boundary layer code to estimate the heating distribution along the surface of a TAV. Thermochemical stability of the coatings is described in terms of reduction in emittance and loss of opacifiers from the coatings during the arc-jet tests.

Stewart, David A.↗

Wind tunnel results of the low-speed NLF(1)-0414F airfoil

The large performance gains predicted for the Natural Laminar Flow (NLF)(1)-0414F airfoil were demonstrated in two-dimensional airfoil tests and in wind tunnel tests conducted with a full scale modified Cessna 210. The performance gains result from maintaining extensive areas of natural laminar flow, and were verified by flight tests conducted with the modified Cessna. The lift, stability, and control characteristics of the Cessna were found to be essentially unchanged when boundary layer transition was fixed near the wing leading edge. These characteristics are very desirable from a safety and certification view where premature boundary layer transition (due to insect contamination, etc.) must be considered. The leading edge modifications were found to enhance the roll damping of the Cessna at the stall, and were therefore considered effective in improving the stall/departure resistance. Also, the modifications were found to be responsible for only minor performance penalties.

Murri, Daniel G.↗

Propulsion/airframe integration issues for waverider aircraft

While many propulsion concepts and technologies developed for nonwaverider-type hypersonic vehicles may apply to waveriders, some aspects of these configurations require unique technological approaches. An evaluation is made of such distinctive opportunities in the cases of engine cycle selection, inlets, nozzle designs and integration, longitudinal stability, and thermal management. Also discussed are waverider requirements for control surface effectiveness, inlet boundary layer ingestion effects, and structural/configurational optimization, giving attention to trades in volumetric/structural efficiency and vehicle L/D.

Blankson, Isaiah M.↗

Stability Analysis of the Flow Over a Swept Forward-Facing Step Using PIV Base Flows

Step excrescences are a type of surface imperfection encountered on swept wings of commercial aircraft, commonly because accessibility requirements prevent creating the wing’s surface from a single panel. If the step height is too large, super-critical, the flow can undergo an early transition to turbulence, which can be highly detrimental for the performance of the wing. In the case of a forward-facing step on a swept wing in a low-disturbance environment, stationary crossflow vortices can develop a significant amplitude upstream of the step and hence dominate the structure of the boundary-layer flow over the step. The main goal of the present investigation is to illuminate the path to transition supported by the fascinatingly complex flow field in the direct downstream vicinity of a step with a super critical height. The high-resolution, stereographic Particle Image Velocimetry (PIV) measurement dataset presently available for this flow field provides a complete description of the laminar flow for the execution of BiGlobal stability analysis in a plane parallel to the step. Although the notorious sensitivity of stability results to the description of the base flow demands a very careful uncertainty analysis of those results, it is argued that this very fact can be leveraged to produce new insight into the supported perturbation dynamics. In performing the analysis, several unsteady mode families are discovered that display the explosive perturbation expected for early transition to be induced. In considering domain widths equal to an integer-multiple of the incident crossflow-vortex wavelength and analyzing an extent of 5 crossflow-vortex wavelengths parallel to the step, it is found that the stability results converge while increasing the domain width. It is demonstrated, moreover, that the results for the wider domains can be approximated by appropriately averaging the results on neighboring single-crossflow-vortex-wavelength domains covering the same region. Besides being useful for computational purposes, this observed property suggests interpreting the instability mechanism as a distorted primary mechanism rather than a “proper” secondary mechanism. This follows in the context of the secondary in-stability analysis of three-dimensional boundary layers, because the secondary mechanism is usually characterized by being localized in a pocket of strong shear, while the distorted primary mechanism typically has an infinite support in the direction parallel to the step. Even though the growth rates are found to be sensitive to the interrogation-window size inherent to the PIV post-processing procedure, the spatial structure of the eigenfunctions is found to be relatively insensitive. Lastly, the spatial structure of the eigenfunctions corresponding to all velocity components are matched with the shape functions determined by computing the Spectral Proper Orthogonal Decomposition (SPOD) of a time-resolved measurement of the perturbation content.

forward-facing step↗

Logarithmic profile of temperature in sheared and unstably stratified atmospheric boundary layers

The impact of buoyancy on the mean velocity, temperature, and scalar concentration profiles in the lower atmosphere is typically investigated within the framework of Monin-Obukhov similarity theory (MOST). MOST is the theoretical foundation for parametrizing surface-atmosphere exchanges in nearly all weather, climate, and hydrological models.According to MOST, the classic logarithmic profiles of mean velocity and temperature break down as the buoyancy effects become important. However, recent studies on turbulent Rayleigh-Bénard convection and natural convection along vertical walls suggest that the mean temperature in the near-surface region still follows a logarithmic profile.Motivated by these new results, we study the mean potential temperature profile in sheared and unstably stratified atmospheric boundary layers using direct numerical simulations and field observations. We find that the mean potential temperature profile remains logarithmic across a wide range of stability parameters, which characterizes the relative importance of buoyancy versus shear effects. Compared to MOST, our results suggest that the buoyancy force does not modify the logarithmic nature of the mean potential temperature profile, but instead modulates its slope, which is no longer universal and differs from 1/κ,whereκisthe von Kármán constant. Furthermore, this study provides another perspective on scalar turbulence in the atmospheric boundary layer

58 GEOSCIENCES↗