skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


Title: Aerodynamic Shape Optimization for Delaying Dynamic Stall of Airfoils Using Cokriging Regression
The dynamic stall phenomenon produces adverse aerodynamic loading, which negatively affects the structural strength and life of aerodynamic systems. Aerodynamic shape optimization (ASO) provides a practical approach for delaying and mitigating dynamic stall characteristics without the addition of an auxiliary system. A typical ASO investigation requires multiple evaluations of accurate but time-consuming computational fluid dynamics (CFD) simulations. In the case of dynamic stall, unsteady CFD simulations are required for airfoil shape evaluation; combining it with high-dimensions of airfoil shape parameterization renders the ASO investigation computationally costly. In this study, metamodel-based optimization (MBO) is proposed using the multifidelity modeling (MFM) technique to efficiently conduct ASO investigation for computationally expensive dynamic stall cases. MFM methods combine data from accurate high-fidelity (HF) simulations and fast low-fidelity (LF) simulations to provide accurate and fast predictions. In particular, Cokriging regression is used for approximating the objective and constraint functions. The airfoil shape is parameterized using six PARSEC parameters. The objective and constraint functions are evaluated for a sinusoidally oscillating airfoil with the unsteady Reynolds-averaged Navier-Stokes equations at a Reynolds number of 135,000, Mach number of 0.1, and reduced frequency of 0.05. The initial metamodel is generated using 220 LF and 20 HF samples. The metamodel is then sequentially refined using the expected improvement infill criteria and validated with the normalized root mean square error. The refined metamodel is utilized for finding the optimal design. The optimal airfoil shape shows higher thickness, larger leading-edge radius, and an aft camber compared to baseline (NACA 0012). The optimal shape delays the dynamic stall occurrence by 3 degrees and reduces the peak aerodynamic coefficients. The performance of the MFM method is also compared with the single-fidelity metamodeling method using HF samples. Both the approaches produced similar optimal shapes; however, the optimal shape from MFM achieved a minimum objective function value while more closely satisfying the constraint at a computational cost saving of around 41%.  more » « less
Award ID(s):
1739551
PAR ID:
10297673
Author(s) / Creator(s):
;
Date Published:
Journal Name:
AIAA SciTech 2021 Forum
Volume:
AIAA 2021
Page Range / eLocation ID:
0340
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; (, Engineering Computations)
    null (Ed.)
    Purpose The purpose of this work is to investigate the similarity requirements for the application of multifidelity modeling (MFM) for the prediction of airfoil dynamic stall using computational fluid dynamics (CFD) simulations. Design/methodology/approach Dynamic stall is modeled using the unsteady Reynolds-averaged Navier–Stokes equations and Menter's shear stress transport turbulence model. Multifidelity models are created by varying the spatial and temporal discretizations. The effectiveness of the MFM method depends on the similarity between the high- (HF) and low-fidelity (LF) models. Their similarity is tested by computing the prediction error with respect to the HF model evaluations. The proposed approach is demonstrated on three airfoil shapes under deep dynamic stall at a Mach number 0.1 and Reynolds number 135,000. Findings The results show that varying the trust-region (TR) radius (λ) significantly affects the prediction accuracy of the MFM. The HF and LF simulation models hold similarity within small (λ ≤ 0.12) to medium (0.12 ≤ λ ≤ 0.23) TR radii producing a prediction error less than 5%, whereas for large TR radii (0.23 ≤ λ ≤ 0.41), the similarity is strongly affected by the time discretization and minimally by the spatial discretization. Originality/value The findings of this work present new knowledge for the construction of accurate MFMs for dynamic stall performance prediction using LF model spatial- and temporal discretization setup and the TR radius size. The approach used in this work is general and can be used for other unsteady applications involving CFD-based MFM and optimization. 
    more » « less
  2. ; ; (, Physics of Fluids)
    In this paper, we investigate the three-dimensional nature of dynamic stall. Conducting the investigation, the flow around a harmonically pitching National Advisory Committee for Aeronautics (NACA) 0012 airfoil is numerically simulated using Unsteady-Reynolds-Averaged Navier–Stokes (URANS) and multiple detached eddy simulation (DES) solvers: the Delayed-DES (DDES) and the Improved-DDES (IDDES). Two- and three-dimensional simulations are performed for each solver, and the results are compared against experimental measurements in the literature. The results showed that three-dimensional simulations surpass two-dimensional ones in capturing the stages of dynamic stall and predicting the lift coefficient values, with a distinguished performance of the DES solvers over the URANS ones. For instance, the IDDES simulations, as an inherently three-dimensional solver, predicted the necessary cascaded amalgamation process of vortices to form the adequate strength of the dynamic stall vortex. This vortex size and timing provided accurate and sufficient suction that resulted in identical matching of the numerical and experimental lift coefficients at the peak value. Hence, the hypothesis that dynamic stall has a three-dimensional nature is supported by the superiority of the three-dimensional simulation in all aspects. In conclusion, it is found that dynamic stall is intrinsically a three-dimensional phenomenon. 
    more » « less
  3. (, Dissertation)
    The rising global trend to reduce dependence on fossil fuels has provided significant motivation toward the development of alternative energy conversion methods and new technologies to improve their efficiency. Recently, oscillating energy harvesters have shown promise as highly efficient and scalable turbines, which can be implemented in areas where traditional energy extraction and conversion are either unfeasible or cost prohibitive. Although such devices are quickly gaining popularity, there remain a number of hurdles in the understanding of their underlying fluid dynamics phenomena. The ability to achieve high efficiency power output from oscillating airfoil energy harvesters requires exploitation of the complexities of the event of dynamic stall. During dynamic stall, the oncoming flow separates at the leading edge of the airfoil to form leading ledge vortex (LEV) structures. While it is well known that LEVs play a significant role in aerodynamic force generation in unsteady animal flight (e.g. insects and birds), there is still a need to further understand their spatiotemporal evolution in order to design more effective energy harvesting enhancement mechanisms. In this work, we conduct extensive experimental investigations to shed-light on the flow physics of a heaving and pitching airfoil energy harvester operating at reduced frequencies of k = fc=U1 = 0.06-0.18, pitching amplitude of 0 = 75 and heaving amplitude of h0 = 0:6c. The experimental work involves the use of two-component particle image velocimetry (PIV) measurements conducted in a wind tunnel facility at Oregon State University. Velocity fields obtained from the PIV measurements are analyzed qualitatively and quantitatively to provide a description of the dynamics of LEVs and other flow structures that may be present during dynamic stall. Due to the difficulties of accurately measuring aerodynamic forces in highly unsteady flows in wind tunnels, a reduced-order model based on the vortex-impulse theory is proposed for estimating the aerodynamic loadings and power output using flow field data. The reduced-order model is shown to be dominated by two terms that have a clear physical interpretation: (i) the time rate of change of the impulse of vortical structures and (ii) the Kutta-Joukowski force which indirectly represents the history effect of vortex shedding in the far wake. Furthermore, the effects of a bio-inspired flow control mechanism based on deforming airfoil surfaces on the flow dynamics and energy harvesting performance are investigated. The results show that the aerodynamic loadings, and hence power output, are highly dependent on the formation, growth rate, trajectory and detachment of the LEV. It is shown that the energy harvesting efficiency increases with increasing reduced frequency, peaking at 25% when k = 0.14, agreeing very well with published numerical results. At this optimal reduced frequency, the time scales of the LEV evolution and airfoil kinematics are matched, resulting in highly correlated aerodynamic load generation and airfoil motion. When operating at k > 0:14, it is shown that the aerodynamic moment and airfoil pitching motion become negatively correlated and as a result, the energy harvesting performance is deteriorated. Furthermore, by using a deforming airfoil surface at the leading and trailing edges, the peak energy harvesting efficiency is shown to increase by approximately 17% and 25% relative to the rigid airfoil, respectively. The performance enhancement is associated with enhanced aerodynamic forces for both the deforming leading and trailing edges. In addition, The deforming trailing edge airfoil is shown to enhance the correlation between the aerodynamic moment and pitching motion at higher reduced frequencies, resulting in a peak efficiency at k = 0:18 as opposed to k = 0:14 for the rigid airfoil. 
    more » « less
  4. ; ; (, AIAA Journal)
    We evaluate two leading-edge-based dynamic stall-onset criteria (namely, the maximum magnitudes of the leading-edge suction parameter and the boundary enstrophy flux) for mixed and trailing-edge stall. These criteria have been shown to successfully predict the onset of leading-edge stall at Reynolds numbers of O(10^5), where the leading-edge suction drops abruptly. However, for mixed/trailing-edge stall, leading-edge suction tends to persist even when there is significant trailing-edge reversed flow and stall is underway, necessitating further investigation into the effectiveness of these criteria. Using wall-resolved large-eddy simulations and the unsteady Reynolds-averaged Navier–Stokes method, we simulate one leading-edge stall and three mixed/trailing-edge stall cases at Reynolds numbers of 200,000 and 300,000. We contrast the progression of flow features such as trailing-edge separation and vortex formation across different stall types and evaluate the stall-onset criteria relative to critical points in the flow. We find that the criteria nearly coincide with the instance of leading-edge suction collapse and are reached in advance of dynamic stall vortex formation and lift stall for all four cases. We conclude that the two criteria effectively signal dynamic stall onset in cases where the dynamic stall vortex plays a prominent role. 
    more » « less
  5. ; ; (, 2022 AIAA SciTech Meeting)
    We evaluate two leading-edge-based dynamic stall onset criteria, namely, the maximum magnitudes of the Leading Edge Suction Parameter and the Boundary Enstrophy Flux, for mixed and trailing-edge stall. These criteria have been shown to successfully predict the onset of leading-edge stall at Reynolds numbers >= O(10^5), where the leading-edge suction drops abruptly. However, for mixed/trailing-edge stall, leading-edge suction tends to persist even when there is significant trailing-edge reversed flow and stall is underway, necessitating further investigation of the effectiveness of these criteria. Using wall-resolved, large-eddy simulations and unsteady Reynolds-Averaged Navier-Stokes method, we simulate one leading-edge stall and three mixed/trailing-edge stall cases at Reynolds numbers 2x10^5 and 3x10^6. We contrast the progression of flow features such as trailing-edge separation and vortex formation across different stall types and evaluate the stall onset criteria relative to critical points in the flow. We find that the criteria nearly coincide with the instance of leading-edge suction collapse and are reached in advance of dynamic stall vortex formation and lift stall for all four cases. We conclude that the two criteria effectively signal dynamic stall onset in cases where the dynamic stall vortex plays a prominent role. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email