Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Wall-normal vortices”

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 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS2018

The Effect of Forward-Facing Steps on Stationary Crossflow Instability Growth and Breakdown

The effect of a forward-facing step on stationary crossflow transition was studied using standard stereo particle image velocimetry (PIV) and time-resolved PIV. Step heights ranging from 53 to 71% of the boundary-layer thickness were studied in detail. The steps above a critical step height of approximately 60% of the boundary-layer thickness had a significant impact on the stationary crossflow growth downstream of the step. For the critical cases, the stationary crossflow amplitude grew suddenly downstream of the step, decayed for a short region, then grew again. The adverse pressure gradient upstream of the step resulted in a region of crossflow reversal. A secondary set of vortices, rotating in the opposite direction to the primary vortices, developed underneath the uplifted primary vortices. The wall-normal velocity disturbance (V' ) created by these secondary vortices impacted the step, and is believed to feed into the strong vortex that developed downstream of the step. A large but very short negative crossflow region formed for a short region downstream of the step due to a sharp inboard curvature of the streamlines near the wall. For the larger step height cases, a crossflow-reversal region formed just downstream of the strong negative crossflow region. This crossflow reversal region is believed to play an important role in the growth of the stationary crossflow vortices downstream of the step, and may be a good indication of the critical forward-facing step height.

Eppink, Jenna L.↗
DOE OSTI2025

Effect of nozzle translation on a bounded vortex flow

A bounded vortex flow is generated by a nozzle that combines azimuthally tilted downward jets in a circular array and a central suction port. When the nozzle is directed toward a bottom surface (called the impingement surface), the flow develops a strong intake vortex spanning from the impingement surface to the suction port, which acts to generate strong shear stress on the impingement surface. This strong shear stress can aid in removing particles attached to the impingement surface or in enhancement of heat or mass transfer, where in typical applications, the nozzle would be translated along the impingement surface. This paper reports on a computational study of the effect of nozzle translation on the bounded vortex flow field. The nozzle translation induces a bending of the bounded vortex in the downstream flow direction, which becomes larger with increase in the nozzle translation speed. We identify an interesting oscillatory phenomenon by which the bounded vortex can resist the bending induced by the translational flow and rotate about a mean position downstream of the suction port. At sufficiently high translation speeds, the oscillatory process can no longer resist the downwash effect from nozzle translation and the vortex eventually detaches from the impingement surface, limiting the effectiveness of this flow device.

42 ENGINEERING↗
DOE OSTI2023

Transition to turbulence in viscoelastic channel flow of dilute polymer solutions

The transition to turbulence in a plane Poiseuille flow of dilute polymer solutions is studied by direct numerical simulations of a finitely extensible nonlinear elastic fluid with the Peterlin closure. The range of Reynolds number ($Re$)$2000 \le Re \le 5000$is studied but with the same level of elasticity in viscoelastic flows. The evolution of a finite-amplitude perturbation and its effects on the transition dynamics are investigated. A viscoelastic flow begins transition at an earlier time than its Newtonian counterparts, but the transition time appears to be insensitive to polymer concentration in the dilute and semi-dilute regimes studied. Increasing polymer concentration, however, decreases the maximum attainable energy growth during the transition process. The critical or minimum perturbation amplitude required to trigger transition is computed. Interestingly, both Newtonian and viscoelastic flows follow almost the same power-law scaling of$Re^\gamma$with the critical exponent$\gamma \approx -1.25$, which is in close agreement with previous studies. However, a shift downward is observed for viscoelastic flow, suggesting that smaller perturbation amplitudes are required for the transition. A mechanism of the early transition is investigated by the evolution of wall-normal and spanwise velocity fluctuations and flow structure. The early growth of these fluctuations and the formation of quasi-streamwise vortices around low-speed streaks are promoted by polymers, hence causing an early transition. These vortical structures are found to support the critical exponent$\gamma \approx -1.25$. Once the transition process is completed, polymers play a role in dampening the wall-normal and spanwise velocity fluctuations and vortices to attain a drag-reduced state in viscoelastic turbulent flows.

Mechanics↗
NASA NTRS1993

Vortex generation and wave-vortex interaction over a concave plate with roughness and suction

The generation and amplification of vortices by surface homogeneities, both in the form of surface waviness and of wall-normal velocity, is investigated using the nonlinear parabolic stability equations. Transients and issues of algebraic growth are avoided through the use of a similarity solution as initial condition for the vortex. In the absence of curvature, the vortex decays as the square root of 1/x when flowing over streamwise aligned riblets of constant height, and grows as the square root of x when flowing over a corresponding streamwise aligned variation of blowing/suction transpiration velocity. However, in the presence of wall inhomogeneities having both streamwise and spanwise periodicity, the growth of the vortex can be much larger. In the presence of curvature, the vortex develops into a Gortler vortex. The 'direct' and 'indirect' interaction mechanisms possible in wave-vortex interaction are presented. The 'direct' interaction does not lead to strong resonance with the flow conditions investigated. The 'indirect' interaction leads to K-type transition.

Bertolotti, Fabio↗
NASA NTRS2019

Nonlinear Görtler Vortices and Their Secondary Instability in a Hypersonic Boundary Layer

Nonlinear development of the Görtler instability over a concave surface gives rise to a highly distorted inflectional flow field in the boundary layer that leads to both wall-normal and spanwise gradients in the flow. Such nonlinear structures are susceptible to strong, high-frequency secondary instabilities that may lead to the onset of laminar-turbulent transition. The present numerical study uses direct numerical simulations and linear secondary instability theory to investigate finite amplitude Görtler vortices and their secondary instability characteristics, respectively, in the hypersonic flow over an axisymmetric cone with a concave aft body. To complement previous studies in the literature wherein the Görtler instability was usually studied for a flat plate and initiated at some upstream location by imposing an eigenfunction as the inflow condition or by blowing and suction at the wall, the present investigation is focused on fully realizable Görtler instability that is excited by an azimuthally periodic array of surface protuberances. Furthermore, while the previous work had mostly focused on the secondary instability of Görtler vortices with cross-plane velocity contours that resembled bell-shaped structures, the present results confirm that fully developed mushroom structures also exist in the hypersonic regime when the Görtler vortex amplitude is sufficiently large. Computations further reveal that the dominant modes of secondary instability in these mushroom-shaped structures correspond to an antisymmetic (i.e., sinuous) “stem” mode that concentrates within the strong, nearly wall-normal internal shear layers surrounding the stem regions underneath the caps of the mushroom structures. Additionally, there exist a multitude of other significantly unstable secondary instability modes of both symmetric and antisymmetric types. Analogous to the secondary instability of crossflow vortices in hypersonic flows, secondary instability modes of both symmetric and antisymmetric types. Analogous to the secondary instability of crossflow vortices in hypersonic flows, secondary instability modes originating from the Mack mode instability play an important role during the nonlinear breakdown process.

Laminar turbulent transition↗
NASA NTRS2016

Influence of Stationary Crossflow Modulation on Secondary Instability

A likely scenario for swept wing transition on subsonic aircraft with natural laminar flow involves the breakdown of stationary crossflow vortices via high frequency secondary instability. A majority of the prior research on this secondary instability has focused on crossflow vortices with a single dominant spanwise wavelength. This paper investigates the effects of the spanwise modulation of stationary crossflow vortices at a specified wavelength by a subharmonic stationary mode. Secondary instability of the modulated crossflow pattern is studied using planar, partial-differential-equation based eigenvalue analysis. Computations reveal that weak modulation by the first subharmonic of the input stationary mode leads to mode splitting that is particularly obvious for Y-type secondary modes that are driven by the wall-normal shear of the basic state. Thus, for each Y mode corresponding to the fundamental wavelength of results in unmodulated train of crossflow vortices, the modulated flow supports a pair of secondary modes with somewhat different amplification rates. The mode splitting phenomenon suggests that a more complex stationary modulation such as that induced by natural surface roughness would yield a considerably richer spectrum of secondary instability modes. Even modest levels of subharmonic modulation are shown to have a strong effect on the overall amplification of secondary disturbances, particularly the Z-modes driven by the spanwise shear of the basic state. Preliminary computations related to the nonlinear breakdown of these secondary disturbances provide interesting insights into the process of crossflow transition in the presence of the first subharmonic of the dominant stationary vortex.

Vortex structure↗
DOE OSTI2020

Some observations concerning "laminarization" in heated vertical tubes

For low "turbulent" Reynolds numbers in strongly-heated vertical gas flow, fundamental results of existing direct numerical simulation (DNS) databases have been examined to deduce differences between cases which "revert" to turbulent and those that yield integral parameters which correspond to laminar flow (henceforth laminarizing or laminarized). The objectives of this paper are (1) to examine the streamwise evolution of near-wall flow structures and (2) to determine which, if any, proposed laminarization parameters could be used to discriminate between turbulent and laminarizing flows resulting. Views of streamlines in r-θ cross sections showed evidence of a ring of irregular-shaped streamwise vortices near the wall at all locations. The transverse advection of these vortices is hypothesized to provide a path of least resistance to the transport of streamwise momentum in the wall-normal direction. It is demonstrated that a steady streamwise laminar vortex can provide apparent turbulent momentum transport (as a consequence of the Reynolds decomposition) such that a reasonable mean velocity profile is predicted near the wall. For the present cases, the values of the non-dimensional radius (y c,w + = Re τ,w ) and of the modified wall Reynolds number did discriminate whether turbulent or laminarizing flow resulted.

99 GENERAL AND MISCELLANEOUS↗
NASA NTRS2010

Development and Breakdown of Goertler Vortices in High Speed Boundary Layers

The nonlinear development of G rtler instability over a concave surface gives rise to a highly distorted stationary flow in the boundary layer that has strong velocity gradients in both spanwise and wall-normal directions. This distorted flow is susceptible to strong, high frequency secondary instability that leads to the onset of transition. For high Mach number flows, the boundary layer is also subject to the second mode instability. The nonlinear development of G rtler vortices and the ensuing growth and breakdown of secondary instability, the G rtler vortex interactions with second mode instabilities as well as oblique second mode interactions are examined in the context of both internal and external hypersonic configurations using nonlinear parabolized stability equations, 2-D eigenvalue analysis and direct numerical simulation. For G rtler vortex development inside the Purdue Mach 6 Ludwieg tube wind tunnel, multiple families of unstable secondary eigenmodes are identified and their linear and nonlinear evolution is examined. The computation of secondary instability is continued past the onset of transition to elucidate the physical mechanisms underlying the laminar breakdown process. Nonlinear breakdown scenarios associated with transition over a Mach 6 compression cone configuration are also explored.

Li, Fei↗
NASA NTRS2020

Mechanisms of Stationary Crossflow Instability Growth and Breakdown Induced by Forward-Facing Steps

An experimental study is performed to determine the mechanisms by which a forward-facing step impacts the growth and breakdown to turbulence of the stationary cross-flow instability. Particle image velocimetry measurements are obtained in the boundary layer of a 30◦swept flat plate with a pressure body. Step heights range from 53 % to 71 % of the boundary-layer thickness. The critical step height is approximately 60 % of the boundary-layer thickness for the current study, although it is also shown that the critical step height depends on the initial amplitude of the stationary cross-flow vortices. For the critical cases, the stationary cross-flow amplitude grows sharply downstream of the step, decays for a short region and then grows again. The initial growth region is linear, and can be explained primarily through the impact of the step on the mean flow. Namely, the step causes abrupt changes to the mean flow, resulting in large values of wall-normal shear, as well as highly inflectional profiles, due to either cross-flow reversal, separation or both. These inflectional profiles are highly unstable for the stationary cross-flow. Additionally,the reversed flow regions are significantly modulated by the stationary cross-flow vortices. The second region of growth occurs due to the stationary-cross-flow-induced modulation of the shear layer, which leads to multiple smaller wavelength stream wise vortices. High-frequency fluctuations indicate that the unsteady transition mechanism for the critical cases relates to the shedding of vortices downstream of reattachment of the modulated separated regions.

boundary layer stability↗
NASA NTRS2011

Effects of Cavities and Protuberances on Transition over Hypersonic Vehicles

Surface protuberances and cavities on a hypersonic vehicle are known to cause several aerodynamic or aerothermodynamic issues. Most important of all, premature transition due to these surface irregularities can lead to a significant rise in surface heating. To help understand laminar-turbulent transition induced by protuberances or cavities on a Crew Exploration Vehicle (CEV) surface, high-fidelity numerical simulations are carried out for both types of trips on a CEV wind tunnel model. Due to the large bluntness, these surface irregularities reside in an accelerating subsonic boundary layer. For the Mach 6 wind tunnel conditions with a roughness Reynolds number Re(sub kk) of 800, it was found that a protuberance with a height to boundary layer thickness ratio of 0.73 leads to strong wake instability and spontaneous vortex shedding, while a cavity with identical geometry only causes a rather weak flow unsteadiness. The same cavity with a larger Reynolds number also leads to similar spontaneous vortex shedding and wake instability. The wake development and the formation of hairpin vortices for both protuberance and cavity were found to be qualitatively similar to that observed for an isolated hemisphere submerged in a subsonic, low speed flat-plate boundary layer. However, the shed vortices and their accompanying instability waves were found to be slightly stabilized downstream by the accelerating boundary layer along the CEV surface. Despite this stabilizing influence, it was found that the wake instability spreads substantially in both wall-normal and azimuthal directions as the flow is evolving towards a transitional state. Similarities and differences between the wake instability behind a protuberance and a cavity are investigated. Computations for the Mach 6 boundary layer over a slender cylindrical roughness element with a height to the boundary layer thickness of about 1.1 also shows spontaneous vortex shedding and strong wake instability. Comparisons of detailed flow structures associated with protuberances at subsonic and supersonic edge Mach numbers indicate distinctively different instability mechanisms.

Chang, Chau-Lyan↗
NASA NTRS2014

Direct Numerical Simulation and Theories of Wall Turbulence with a Range of Pressure Gradients

A new Direct Numerical Simulation (DNS) of Couette-Poiseuille flow at a higher Reynolds number is presented and compared with DNS of other wall-bounded flows. It is analyzed in terms of testing semi-theoretical proposals for universal behavior of the velocity, mixing length, or eddy viscosity in pressure gradients, and in terms of assessing the accuracy of two turbulence models. These models are used in two modes, the traditional one with only a dependence on the wall-normal coordinate y, and a newer one in which a lateral dependence on z is added. For pure Couette flow and the Couette-Poiseuille case considered here, this z-dependence allows some models to generate steady streamwise vortices, which generally improves the agreement with DNS and experiment. On the other hand, it complicates the comparison between DNS and models.

Coleman, G. N.↗
DOE OSTI2023

Uniform momentum and temperature zones in unstably stratified turbulent flows

Wall-bounded turbulent flows exhibit a zonal arrangement, in which streamwise velocity organizes into uniform momentum zones (UMZs), separated by thin layers of elevated interfacial shear. While significant research efforts have focused on these structural features in neutrally stratified flows, the effects of unstable thermal stratification on UMZs and on analogous uniform temperature zones (UTZs) have not been considered previously. In this article, statistical properties of UMZs and UTZs are investigated using a suite of large eddy simulations of unstably stratified turbulent channel flow spanning weakly to highly convective conditions. When normalized by the friction velocity and stability-dependent mixing length, the mean velocity gradient based on UMZ interfacial velocity jumps and the vorticity thickness exhibits good collapse for all stabilities, establishing a link between UMZ properties and scaling predictions from Monin–Obukhov similarity theory. A similar relationship is found between UTZ properties and surface-layer scaling of the mean temperature gradient. In the mixed layer, mean UMZ depth is quasi-constant with wall-normal distance, while the deepest UTZs are found in the centre of the boundary layer. These instantaneous structures are found to be linked to the well-mixed velocity and temperature profiles in the convective mixed layer. Conditional averaging indicates that both UMZ and UTZ interfaces are associated with ejections of momentum and warm updrafts below the interface and sweeps of momentum and cool downdrafts above the interface. These results demonstrate a tangible connection between instantaneous structural features, mean properties and scaling laws in unstably stratified flows.

Mechanics↗
NASA NTRS2014

DNS of Laminar-Turbulent Transition in Swept-Wing Boundary Layers

Direct numerical simulation (DNS) is performed to examine laminar to turbulent transition due to high-frequency secondary instability of stationary crossflow vortices in a subsonic swept-wing boundary layer for a realistic natural-laminar-flow airfoil configuration. The secondary instability is introduced via inflow forcing and the mode selected for forcing corresponds to the most amplified secondary instability mode that, in this case, derives a majority of its growth from energy production mechanisms associated with the wall-normal shear of the stationary basic state. An inlet boundary condition is carefully designed to allow for accurate injection of instability wave modes and minimize acoustic reflections at numerical boundaries. Nonlinear parabolized stability equation (PSE) predictions compare well with the DNS in terms of modal amplitudes and modal shape during the strongly nonlinear phase of the secondary instability mode. During the transition process, the skin friction coefficient rises rather rapidly and the wall-shear distribution shows a sawtooth pattern that is analogous to the previously documented surface flow visualizations of transition due to stationary crossflow instability. Fully turbulent features are observed in the downstream region of the flow.

Duan, L.↗
NASA NTRS2013

Direct Numerical Simulation of Transition in a Swept-Wing Boundary Layer

Direct numerical simulation (DNS) is performed to examine laminar to turbulent transition due to high-frequency secondary instability of stationary crossflow vortices in a subsonic swept-wing boundary layer for a realistic natural-laminar-flow airfoil configuration. The secondary instability is introduced via inflow forcing derived from a two-dimensional, partial-differential-equation based eigenvalue computation; and the mode selected for forcing corresponds to the most amplified secondary instability mode which, in this case, derives a majority of its growth from energy production mechanisms associated with the wall-normal shear of the stationary basic state. Both the growth of the secondary instability wave and the resulting onset of laminar-turbulent transition are captured within the DNS computations. The growth of the secondary instability wave in the DNS solution compares well with linear secondary instability theory when the amplitude is small; the linear growth is followed by a region of reduced growth resulting from nonlinear effects before an explosive onset of laminar breakdown to turbulence. The peak fluctuations are concentrated near the boundary layer edge during the initial stage of transition, but rapidly propagates towards the surface during the process of laminar breakdown. Both time-averaged statistics and flow visualization based on the DNS reveal a sawtooth transition pattern that is analogous to previously documented surface flow visualizations of transition due to stationary crossflow instability. The memory of the stationary crossflow vortex is found to persist through the transition zone and well beyond the location of the maximum skin friction.

Duan, Lian↗
DOE OSTI2022

Unstable turbulent channel flow response to spanwise-heterogeneous heat fluxes: Prandtl's secondary flow of the third kind

Turbulent secondary flows are defined as Prandtl's secondary flow of the first or second kind, the former produced by stretching and/or tilting of vorticity, the latter produced via spatial heterogeneity of Reynolds stresses. Both mechanisms are instantaneously active within inertia-dominated wall turbulence; Reynolds stress spatial heterogeneity is required for Reynolds-averaged secondary flows. Spanwise-variable surface roughness can induce turbulent stress spatial heterogeneity in the spanwise–wall-normal plane and provide sustenance for streamwise-aligned mean secondary flows. Herein, we demonstrate that turbulent secondary flows can also be sustained by spanwise variability in the surface heat flux in unstably stratified turbulent channels, defined hereafter as Prandtl's secondary flow of the third kind. Support for this mechanism is established with scaling arguments, while large-eddy simulation is used to model inertia-dominated channel turbulence responding to a lower boundary with uniform aerodynamic/hydrodynamic roughness but spanwise-variable surface heat flux. Transport equations for streamwise vorticity and turbulent kinetic energy, k, outline the conditions needed for third-kind production: shear and buoyancy production over the elevated heat flux regions necessitates lateral entrainment of low- k fluid, inducing mean counter-rotating secondary cells aligned such that upwelling and downwelling occur over the high and low heat flux regions, respectively. Buoyancy-driven production of k alters aggregate flow response and thus is a distinctly different mechanism responsible for sustenance of secondary flows.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1999

Swept-Wing Receptivity Studies Using Distributed Roughness

Transition to turbulence in swept-wing flows has resisted correlation with linear theory because of its sensitivity to freestream conditions and 3-D roughness and because one of the principal instability modes quickly 'becomes nonlinear. In the face of such a formidable problem, two rather long-term fundamental efforts have been underway at DLR Gottinberg and Arizona State University that address swept-wing transition. These efforts have been recently reviewed by Bippes (1997) and Reibert and Saric (1997). Thus, the present work is a continuation of a series of studies on swept-wing boundary layers which have led to a better understanding of the transition process. In particular, we have taken advantage of the sensitivity to 3-D roughness and the modal nature of the instability in order to propose a particular control strategy. Complementing the two aforementioned reviews, general reviews of the swept-wing transition problem are found in Arnal (1997) and Kachanov (1996). Other recent reviews include Reshotko (t997), Crouch (1997), and Herbert (1997a,b). The failure of linear theory is discussed in Reed et al. (1996). The historical work is found in Reed and Sar-ic (1989). The basic idea is that the combination of sweep and chordwise pressure gradient within the boundary layer creates a velocity component perpendicular to the inviscid streamline. This crossflow profile is inflectional and exhibits both traveling and stationary unstable waves called crossflow vortices that are (approximately) aligned along the inviscid streamlines. Under conditions of low freestream turbulence levels, the dominant crossflow wave is stationary (Reibert and Saric t997) while moderate to high turbulence levels initiate dominant traveling waves (Dehle and Bippes 1996; Bippes 1997). 'Me mechanism is relatively insensitive to sound and 2-D surface roughness (Radeztsky et al. 1993) but very sensitive to 3-D roughness near the attachment line. We concentrate our work on low-turbulence freestream flows and stationary crossflow waves. Although the v' and w' components of the disturbances are very small, by convecting streamwise momentum in the wall-normal direction, they produce 0(l) changes inu'.Thus the mean flow is highly distorted with localized inflection points. Transition is then triggered by a high-frequency secondary instability of the distorted mean profile.

Saric, William S.↗
DOE OSTI2023

Airfoil Computational Fluid Dynamics - 2k shapes, 25 AoA's, 3 Re numbers

This dataset contains aerodynamic quantities - including flow field values (momentum, energy, and vorticity) and summary values (coefficients of lift, drag, and momentum) - for 1,830 airfoil shapes computed using the HAM2D CFD (computational fluid dynamics) model. The airfoil shapes were designed using the separable shape tensor parameterization that encodes two-dimensional shapes as elements of the Grassmann manifold. This data-driven approach learns two independent spaces of parameter from a collection of sample airfoils. The first captures large-scale, linear perturbations, and the second defines small-scale, higher-order perturbations. For this dataset, we used the G2Aero database of over 19,000 airfoil shapes to learn a parameter space that captured a wide array of shape characteristics. We sampled airfoil designs over both parameter spaces to explore the full range of possible shape variations. The aerodynamic quantities for the generated airfoil were obtained using the HAM2D code, which is a finite-volume Reynolds-averaged Navier-Stokes (RANS) flow solver. We employ a fifth-order WENO scheme for spatial reconstruction with Roe's flux difference scheme for inviscid flux and second-order central differencing for viscous flux. A preconditioned GMRES method is applied for implicit integration. The Spalart-Allmaras 1-eq turbulence model is used for the turbulence closure, and the Medida-Baeder 2-eq transition model is applied to account for the effects of laminar turbulent transition. The airfoil grid is generated with a total of 400 points on the airfoil surface, the initial wall-normal spacing of y+ = 1, and an outer boundary located at 300 chord lengths away from the wall. The CFD simulations are performed at a freestream Mach number of 0.1, for or three different Reynolds' numbers (3M, 6M, and 9M), and for 25 angles of attack from -4 deg. to 20 deg. with 1 degree increments. Across all these various parameters, this dataset includes the results from over 250,000 CFD simulations. The simulations were performed using the Bridges-2 system at the Pittsburgh Supercomputing Center in February 2023 as part of the INTEGRATE project funded by the Advanced Research Projects Agency - Energy, in the U.S. Department of Energy. The data was collected, reformatted, and preprocessed for this OEDI submission in July 2023 under the Foundational AI for Wind Energy project funded by the U.S. Department of Energy Wind Energy Technologies Office. This dataset is intended to serve as a benchmark against which new artificial intelligence (AI) or machine learning (ML) tools may be tested. Baseline AI/ML methods for analyzing this dataset have been implemented, and a link to their repository containing those models has been provided. The .h5 data file structure can be found in the GitHub Repository resource under explore_airfoil_2k_data.ipynb.

2k↗
DOE OSTI2023

Airfoil Computational Fluid Dynamics - 9k shapes, 2 AoA's

This dataset contains aerodynamic quantities - including flow field values (momentum, energy, and vorticity) and summary values (coefficients of lift, drag, and momentum) - for 8,996 airfoil shapes, computed using the HAM2D CFD (computational fluid dynamics) model. The airfoil shapes were designed using the separable shape tensor parameterization that encodes two-dimensional shapes as elements of the Grassmann manifold. This data-driven approach learns two independent spaces of parameter from a collection of sample airfoils. The first captures large-scale, linear perturbations, and the second defines small-scale, higher-order perturbations. For this data, we used the G2Aero database of over 19,000 airfoil shapes to learn a parameter space that captured a wide array of shape characteristics. We fixed the linear deformations to be the mean over the database and sampled new shapes over a four-dimensional parameter space of higher-order perturbation. This sampling approaches allows for isolated analysis of non-linear airfoil shape deformations while holding other aspects (e.g., airfoil thickness) approximately constant. The aerodynamic quantities for the generated airfoil were obtained using the HAM2D code, which is a finite-volume Reynolds-averaged Navier-Stokes (RANS) flow solver. We employ a fifth-order WENO scheme for spatial reconstruction with Roe's flux difference scheme for inviscid flux and second-order central differencing for viscous flux. A preconditioned GMRES method is applied for implicit integration. The Spalart-Allmaras 1-eq turbulence model is used for the turbulence closure, and the Medida-Baeder 2-eq transition model is applied to account for the effects of laminar turbulent transition. The airfoil grid is generated with a total of 400 points on the airfoil surface, the initial wall-normal spacing of y+ = 1, and an outer boundary located at 300 chord lengths away from the wall. The CFD simulations are performed at a freestream Mach number of 0.1, Reynolds number of 9M, and at two angles of attack, 4 deg. and 12 deg. The simulations were performed using the Bridges-2 system at the Pittsburgh Supercomputing Center in February 2023 as part of the INTEGRATE project funded by the Advanced Research Projects Agency - Energy in the U.S. Department of Energy. The data was collected, reformatted, and preprocessed for this OEDI submission in July 2023 under the Foundational AI for Wind Energy project funded by the U.S. Department of Energy Wind Energy Technologies Office. This dataset is intended to serve as a benchmark against which new artificial intelligence (AI) or machine learning (ML) tools may be tested. Baseline AI/ML methods for analyzing this dataset have been implemented, and a link to their repository containing those models has been provided. The .h5 data file structure can be found in the GitHub Repository resource under explore_airfoil_9k_data.ipynb.

9k↗