Navier-Stokes Computational Fluid Dynamics Research Articles

A viscous investigation on the hydrodynamic coefficients and wave loads under the interaction of side-by-side cylinders in regular waves

The paper introduces a two-dimensional Reynolds-Averaged Navier-Stokes Computational Fluid Dynamics (RANS CFD) model to investigate the effects of adjacently floating bodies on the hydrodynamic coefficients and wave loads acting on these bodies. The analysis considers two square cylinders with a narrow gap in way of the free water surface. The physical significance of subsections of the hydrodynamic coefficient matrix becomes evident when comparing coefficients obtained from oscillating-oscillating versus oscillating-fixed models. This understanding is significant, particularly in applications where only specific portions of the hydrodynamic coefficient matrix are utilized, such as in ship-ice collision scenarios (Jiang et al., 2023a). The correlation between wave loads and wave elevations reveals that the sway force and roll moment are associated with the disparity in wave elevation at both sides of each cylinder, while heave forces correlate with the averaged wave elevation on both sides of each cylinder. It is concluded that hydrodynamic interactions are sensitive to the displacement of a floater and accordingly can dominate added mass and damping effects.

Characterization of Out-of-Plane Curved Fluidic Oscillators

A novel fluidic oscillator design, in which the oscillator is curved along the primary flow direction, was evaluated. The effects of mass flow rate, hydraulic diameter, aspect ratio, curvature radius, surface roughness, inlet orientation, and nonconstant curvatures were experimentally investigated. Measurements of the oscillation Strouhal number, discharge coefficient, spreading angle, and sweeping angle were used to characterize the jet. The resulting oscillators were found to fall into one of three categories: oscillated at the same frequency and sweep angle as conventional flat oscillators; oscillated at a slightly higher frequency and lower sweep angle than flat oscillators; or no dominant oscillation frequency detected and with no sweeping action. An unsteady Reynolds-averaged Navier–Stokes computational fluid dynamics simulation revealed fundamental differences in the internal flow mechanisms between flat and curved oscillators that drive the sweeping jet. The curvature between 37.5 and 62.5% of the total length, or the region from the inlet nozzle to halfway through the main chamber, was a primary factor influencing the response type of a design. Due to the curvature of these oscillators, they have the ability to be used in geometrically constrained spaces, such as the leading edge of wings and turbine vanes.

Hypersonic FLEET velocimetry and uncertainty characterization in a tripped boundary layer

Femtosecond laser electronic excitation tagging (FLEET) velocimetry is applied in a hypersonic boundary layer behind an array of turbulence-inducing trips. One-dimensional mean velocity and root-mean-square (RMS) of velocity fluctuation profiles are extracted from FLEET emissions oriented across a 2.75∘ wedge and through a boundary layer above a flat plate in two test campaigns spanning 21 tunnel runs. The experiment was performed in the Texas A&M University Actively Controlled Expansion tunnel that operated near Mach 6.0 with a Reynolds number near 6 × 106 m−1 and a working fluid of air at a density near 2.5 × 10−2 kg m−3. Detailed analysis of random and systematic errors was performed using synthetic curves for error in the mean velocity due to emission decay and the error in the RMS velocity fluctuation due to random error. The boundary layer behind an array of turbulence-inducing trips is documented to show the breakdown of coherent structures. FLEET velocimetry is compared to the tunnel Data Acquisition System, Vibrationally Excited Nitric Oxide Monitoring results, and Reynolds-Averaged Navier–Stokes computational fluid dynamics to verify results.

Computational Fluid Dynamics Turbulence Model and Experimental Study for a Fontan Cavopulmonary Assist Device.

Head-flow HQ curves for a Fontan cavopulmonary assist device (CPAD) were measured using a blood surrogate in a mock circulatory loop and simulated with various computational fluid dynamics (CFD) models. The tests benchmarked the CFD tools for further enhancement of the CPAD design. Recommended Reynolds-Averaged Navier-Stokes (RANS) CFD approaches for the development of conventional ventricular assist devices (VAD) were found to have shortcomings when applied to the Fontan CPAD, which is designed to neutralize off-condition obstruction risks that could contribute to a major adverse event. The no-obstruction condition is achieved with a von Karman pump, utilizing large clearances and small blade heights, which challenge conventional VAD RANS-based CFD hemodynamic simulations. High-fidelity large eddy simulation (LES) is always recommended; however, this may be cost-inhibitive for optimization studies in commercial settings, thus the reliance on RANS models. This study compares head and power predictions of various RANS turbulence models, employing experimental measurements and LES results as a basis for comparison. The models include standard k-ϵ, re-normalization group k-ϵ, realizable k-ϵ, shear stress transport (SST) k-ω, SST with transitional turbulence, and Generalized k-ω. For the pressure head predictions, it was observed that the standard k-ϵ model provided far better agreement with experiment. For the rotor torque, k-ϵ predictions were 30% lower than LES, while the SST and LES torque values were near identical. For the Fontan CPAD, the findings support using LES for the final design simulations, k-ϵ model for head and general flow simulation, and SST for power, shear stress, hemolysis, and thrombogenicity predictions.

First principles simulation of reacting hypersonic flow over a blunt wedge

This article presents molecular-level analysis of a reactive, near-continuum, Mach 21 nitrogen flow over a blunt wedge using the direct molecular simulation (DMS) method. The flow conditions lead to internal energy excitation and dissociation in the flow field, resulting in thermal and chemical nonequilibrium in the flow. Thermal nonequilibrium in the vibrational mode is observed to extend to the molecular level, where the vibrational energy distributions at various points in the flow field are observed to be non-Boltzmann. Furthermore, this is the first reactive DMS calculation where the wall is assumed to be isothermal and full momentum accommodation of the particles is enforced, hence incorporating viscous wall effects. Since the DMS method uses a quantum mechanically generated interaction potential as its only modeling input, all thermochemical and transport properties of the flow field can directly be attributed to the ab initio potential energy surface. Using the DMS solution as a benchmark, this article assesses the performance of Navier–Stokes computational fluid dynamics solutions using lower fidelity two-temperature models. Two models are chosen as points of comparison: the well-known Park two-temperature model and the recently developed modified Marrone and Treanor model.

Aerodynamic Modeling of a Flying Wing Featuring Ludwig Prandtl’s Bell Spanload

This paper presents the aerodynamic modeling of a flying wing featuring Ludwig Prandtl’s bell spanload. The aerodynamic models are developed using a medium fidelity vortex-lattice method and using a Reynolds-Averaged Navier–Stokes computational fluid dynamics solution across a wide range of in-flow angle conditions. A methodology is developed to directly compare the spanwise force distributions from each method. Excellent agreement is seen in the prediction of spanwise inertial aerodynamic forces from each method, as well as in stability and control derivatives dominated by lift and moment contributions. The phenomenon of proverse yaw, caused by the twist distribution necessary to produce the bell spanload, is seen in the high and low-order analysis, with similar off-axis control power predicted.

CFD Analysis of Biofouling Effect on Submarine Resistance and Wake

It is well known that biofouling increases a ship’s resistance and nominal wake. For submarines, any change to the circumferential variation of the nominal wake in the propeller plane will affect the variation of the flow over the propeller blade, and hence the fluctuating forces, and noise, generated by the propeller. The ANSYS FLUENT commercial Reynolds-Averaged Navier Stokes Computational Fluid Dynamics solver was used to investigate the influence of both the longitudinal and vertical distribution of biofouling on the resistance and wake, including the circumferential variation of the nominal wake, on a submarine, using the well-known Suboff standard submarine. For the present work, the k-ε turbulence model was selected, as this is commonly used in this field and is generally considered acceptable. To handle different boundary layer thicknesses in the flow fields, the whole y+ formulation was employed, enabling automatic switching between low and high Reynolds boundary wall models. The numerical solver used for the simulations is based on the finite volume method, which discretizes the RANS equations. In this approach, a segregated model was utilized in the solver, and the convection terms were discretized using the second-order upwind scheme to enhance solution accuracy. The criteria for the near wall are between 30 and 100, and the value of y+ for the present case is 84. It is shown that fouling over only the forward third of the submarine results in a greater increase in resistance than fouling over only the aft third. Fouling over only the lower half of the submarine results in greater resistance than fouling over only the forward third, but less than fouling over the whole of the hull. Fouling over only the forward third of the hull has less influence on the circumferential variation of the wake than fouling over the aft third only of the hull. The results show the importance of keeping the forward area of the hull clean when considering resistance only, whereas keeping the aft area of the hull clean is important when considering the uniformity of the nominal wake into the propeller.

Assessment of advanced RANS turbulence models for prediction of complex flows in compressors

Reynolds-Averaged Navier-Stokes (RANS) Computational Fluid Dynamics (CFD) has been widely used in compressor design and analysis. However, reasonable prediction of compressor flow and its impact on compressor performance remains challenging. In this study, Menter’s Shear Stress Transport (SST) model and its variants, as well as the ω-based Reynolds stress model (Stress-BSL) are assessed. For a single rotor (Rotor 67), under the peak efficiency operating condition, all studied turbulence models predict its performance with reasonable accuracy; under the off-design conditions, SST with Helicity correction (SST-Helicity) shows superiority in predicting the effect of flow on the spanwise distribution of aerodynamic parameters. For Darmstadt’s 1.5-stage transonic axial compressor, SST-Helicity outperforms SST, SST with the Quadratic Constitutive Relation (SST-QCR) and Stress-BSL in predicting the performance as well as the spanwise distribution of aerodynamic parameters. At the design rotating speed, the stall margin given by SST-Helicity (20.90%) is the closest to the experimental measurement (24.81%), which is more than twice that by SST (8.71%) and 1.5 times that by SST-QCR (14.14%). This paper demonstrates that SST-Helicity model, together with a good quality and sufficiently refined grid, can capture the compressor flow features with reasonable accuracy, which results in a credible prediction of compressor performance and stage matching.

Prediction of disc golf drivers’ aerodynamic characteristics using Reynolds Averaged Navier Stokes computational fluid dynamics

Prediction of disc golf drivers’ aerodynamic characteristics using Reynolds Averaged Navier Stokes computational fluid dynamics

Interaction Mechanism of Transonic Squealer Tip Cooling With the Effect of High-Speed Relative Casing Motion

Abstract The relative casing motion can significantly influence the turbine blade tip aerothermal performance. In this study, experimental investigation was conducted in a newly developed high-speed disk rotor rig which can mimic engine realistic high-speed casing relative motion while enabling full optical access to a transonic turbine blade tip surface. Spatially-resolved tip heat transfer data, including heat transfer coefficient and film cooling effectiveness, were obtained for a cooled transonic squealer tip by infrared transient thermal measurement. Combined with closely coupled Reynolds-averaged Navier–Stokes computational fluid dynamics (CFD) analysis, this paper reveals an interesting interaction mechanism between the cooling injections from the pressure side and the cavity floor with and without the effect of relative casing motion. Both experimental data and CFD results show a consistent trend in both heat transfer and cooling performance. With cavity cooling only, the cooling performance reduces with the effect of relative casing motion. However, with additional cooling injection from the pressure side, a significant improvement in the combined cooling performance with the relative casing motion can be observed. Such opposite trend highlights the importance of relative casing motion when ranking different tip cooling designs. With the consideration of relative casing motion, extra tip cooling benefit can be obtained by combining cooling injections from two different locations.

The Effect of a Linear Free Surface Boundary Condition on the Steady-State Wave-Making of Shallowly Submerged Underwater Vehicles

Near-surface simulation methods for shallowly submerged underwater vehicles are necessary for the population of a variety of free-surface-affected, coefficient-based maneuvering and seakeeping models. Simulations vary in complexity and computational costs, often sacrificing accuracy for simplicity and speed. One particular simplifying assumption, the linearization of the free surface boundary conditions, is explored in this study by comparing the steady-state wave-making characteristics of a shallowly submerged prolate spheroid using two different simulation methods at several submergence depths and forward speeds. Hydrodynamic responses are compared between a time-domain boundary element method that makes use of a linearized free surface boundary condition and an inviscid, volume of fluid Reynolds-Averaged Navier–Stokes computational fluid dynamics code that imposes no explicit free surface boundary condition. Differences of up to 22.6%, 32.5%, and 33.3% are found in the prediction of steady state surge force, heave force, and pitch moment, respectively. The largest differences between the two simulation methods arise for motions occurring at small submergences and large wave-making velocities where linear free-surface assumptions become less valid. Nonlinearities that occur in such cases are revealed through physical artifacts such as wave steepening, wave breaking, and high-energy waves. A further examination of near-surface viscous forces reveals that the viscous drag on the vessel is depth dependent due to the changing velocity profile around the body.

Numerical and Experimental Analysis of the Cold Flow Physics of a Nonpremixed Industrial Gas Burner

Abstract The flow field of a nonpremixed industrial gas burner is analyzed with Reynolds-averaged Navier–Stokes computational fluid dynamics validated against velocity and pressure measurements. Combustion is not modeled because the aim is optimizing the predictive capabilities of the cold flow before including chemistry. The system's complex flow physics, affected by a 90 deg turn, backward and forward facing steps, and transversal jets in the mainstream is investigated at full and partial load. The sensitivity of the computed flow field to inflow boundary condition setup, approach for resolving/modeling wall-bounded flows, and turbulence closure is assessed. In the first sensitivity analysis, the inflow boundary condition is prescribed using measured total pressure or measured velocity field. In the second, boundary layers are resolved down to the wall or modeled with wall functions. In the third sensitivity analysis, the turbulence closure uses the k−ω shear stress transport eddy viscosity model or two variants of the Reynolds stress model. The agreement between the predictions of most simulation setups among themselves and with the measurements is good. For a given type of inflow condition and wall flow treatment, the ω-based Reynolds stress model gives the best agreement with measurements among the considered turbulence models at full load. At partial load, the comparison with measured data highlights some scatter in the predictions of different patterns of the flow measurements. Overall, the findings of this study provide insight into the fluid dynamics of industrial gas burners and guidelines for their simulation-based analysis.

High-Fidelity Aerostructural Optimization with a Geometrically Nonlinear Flutter Constraint

When designing aircraft, avoiding dynamic aeroelastic instabilities such as flutter is a key requirement. One way to meet this requirement is to use a multidisciplinary design optimization subject to a flutter constraint. Flutter-constrained design optimizations have used geometrically linear detailed models, which do not accurately predict flutter for very flexible aircraft, or geometrically nonlinear low-order models, which do not accurately trade off cruise range and structural mass. This paper presents a framework for integrating a geometrically nonlinear, subsonic flutter constraint that captures the large in-flight deflections of very flexible aircraft into high-fidelity gradient-based aerostructural optimization. The cruise range and stress constraints are computed accurately with detailed aerostructural analyses, which use a built-up finite element model coupled to Reynolds-averaged Navier–Stokes computational fluid dynamics. The detailed model is condensed to a low-order aeroelastic model to compute the geometrically nonlinear flutter constraint and its adjoint derivatives with computational cost and robustness suitable for optimization. The framework is demonstrated by maximizing the cruise range of a subsonic high-aspect-ratio wing with respect to panel thicknesses, sweep, and span. The impact of the geometrically nonlinear flutter constraint highly depends on the in-flight deformation level. At low deflection levels, sweeping the wing backward is the most effective flutter prevention method, and the flutter constraint has little impact on the achievable cruise range. At high deflection levels, shortening the wing span is necessary to suppress flutter, reducing the achievable cruise range by up to 5.9%. This work is a step toward making multidisciplinary design optimization a practical tool for designing the energy-efficient, very flexible aircraft of the future, which require geometrically nonlinear flutter analyses early in the design cycle to prevent flutter.

A Review of Solution Stabilization Techniques for RANS CFD Solvers

Nonlinear, time-linearized and adjoint Reynolds-averaged Navier-Stokes (RANS) computational fluid dynamics (CFD) solvers are widely used to assess and improve the aerodynamic and aeroelastic performance of aircrafts and turbomachines. While RANS CFD solver technologies are relatively mature for applications at design conditions where the flow is benign, their use in off-design conditions, featuring flow instabilities, such as separations and shock wave/boundary layer interactions, still faces many challenges, with tight residual convergence being a major difficulty. To cope with this, several solver stabilization techniques have been proposed. However, a systematic and comparative study of these techniques has not been reported, to some extent hindering the wide deployment of these methods for industrial applications. In this paper, we critically review the existing methods for solver convergence stabilization, with the main purpose of explaining the rationale behind the algorithms and providing a systematic view of the seemingly different methods. Specifically, mathematical formulations and implementation details of these methods, example applications, and the pros and cons of the methods are discussed in detail, along with suggestions for further improvements. This review is expected to give CFD method developers an overview of the various solution stabilization methods and application engineers an idea how to choose a suitable method for their respective applications.

CFD study to predict the effect of a passing ship on moored ships in a confined waterway

This study aims to present a practical practice to obtain a reliable numerical solution based on an Unsteady Reynolds Averaged Navier-Stokes (URANS) Computational Fluid Dynamics (CFD) to estimate the forces acting on moored ships during a passage of a ship in a confined waterway. The motion of the passing ship is implemented using an overset-grid method. The convergence of the numerical solutions depending on the spatial and temporal discretization are systematically analyzed, and the computational results are validated with a recent benchmark test data. In particular, the effect of the initial state of the passing ship in the unsteady simulation is dealt with in-depth. The simulation results when the passing ship departs at a constant speed and when it accelerates from a standstill are compared, and the influence of undesirable wave radiation generated with the departure of the passing ship at a constant speed is discussed.

ДОСЛІДЖЕННЯ МОРЕХІДНОСТІ РИБАЛЬСЬКИХ СУДЕН ЗА ДОПОМОГОЮ CFD МЕТОДІВ

Seakeeping performance is quite important for certain ship types. Fishing vessels often operate in areas with frequent storms, and sometimes it is necessary to continue fishing despite the worsening sea state. The development of hull forms providing good seakeeping performance of fishing vessels is the problem of current interest in ship design. The conventional approach to seakeeping studies is testing models in ship model basins. However, it is time-consuming and expensive, especially when many hull form variants are studied. For this reason, computer calculations based on the theory of ship motions and strip theory were developed and introduced. Today the more advanced methods of computational fluid dynamics (CFD) can be applied to the problem. The study of ship motions with the help of Reynolds-averaged Navier-Stokes (RANS) CFD method is considered in this paper. A suggested numerical model implies the generation of waves through initial and boundary conditions, which express fully developed waves with preset parameters. An object of research is the seagoing trawler with an overall length of 44.6 m. Three versions of hull form are used in simulations. All versions have similar transom afterbodies, but different forebodies: one similar to Axe Bow, a bulbous bow and one similar to X-bow. Ship movement with headings 180° and 150° at speed of 3.5 knots was studied to reveal differences in added resistance and ship motions. While added resistance and characteristics of heaving and pitching have shown no clear advantages among the considered shapes of forebodies, the version similar to X-bow has demonstrated a significant decrease in rolling at heading 150°. The results of the study have shown that the suggested numerical setup in combination with the CFD methods described can be used for quite realistic simulations of ship behaviour in rough seas. Keywords: seakeeping simulations, seakeeping of fishing vessels, CFD for ship design, improved seakeeping performance, simulation of ship motions.

Surrogate Model Development for Optimized Blended-Wing-Body Aerodynamics

In the conceptual design phase of conventional-configuration aircraft, calibrated low-fidelity methods provide sufficiently accurate estimates of aerodynamic coefficients. It has been observed, however, that for blended-wing–body aircraft, important flow effects are not captured adequately with low-fidelity aerodynamic tools. Consequently, high-fidelity methods become necessary to study blended-wing–body aerodynamics. Since repeated function calls are needed in an optimization loop, high-fidelity analysis is prohibitively expensive in the conceptual design phase, where several optimization scenarios are considered. In this paper, the integration of high-fidelity data for blended-wing–body aircraft for a mission calculation module is presented. A surrogate model based on Gaussian processes (GPs) with acceptably low prediction error is sought as an alternative to Reynolds-averaged Navier–Stokes computational fluid dynamics. Three adaptations are considered: sparse GPs, mixtures of GP approximators, and need-based filtering for GP. The results provide benchmark values for this case and show that the combination of subsonic and transonic behaviors in the training set is problematic and that, for the considered datasets, sparse GP models suffer from oversmoothing, while mixtures of GP models suffer from overfitting. From the error levels, it is observed that a GP with an infinitely differentiable squared exponential kernel based on reduced data pertinent to mission analysis is the most effective option.

Atmospheric rotating rig testing of a swept blade tip and comparison with multi-fidelity aeroelastic simulations

Abstract. One promising design solution for increasing the energy production of modern horizontal axis wind turbines is the installation of curved tip extensions. However, since the aeroelastic response of such geometrical add-ons has not been characterized yet, there are currently uncertainties in the application of traditional aerodynamic numerical models. The objective of the present work is twofold. On the one hand, it represents the first effort in the experimental characterization of curved tip extensions in atmospheric flow. On the other hand, it includes a comprehensive validation exercise, accounting for different numerical models for aerodynamic load prediction. The experiments consist of controlled field tests at the outdoor rotating rig at the Risø campus of the Technical University of Denmark (DTU), and consider a swept tip shape. This geometry is the result of an optimized design, focusing on locally maximizing power performance within load constraints compared to an optimal straight tip. The tip model is instrumented with spanwise bands of pressure sensors and is tested in atmospheric inflow conditions. A range of fidelities of aerodynamic models are then used to aeroelastically simulate the test cases and to compare with the measurement data. These aerodynamic codes include a blade element momentum (BEM) method, a vortex-based method coupling a near-wake model with a far-wake model (NW), a lifting-line hybrid wake model (LL), and fully resolved Navier–Stokes computational fluid dynamics (CFD) simulations. Results show that the measured mean normal loading can be captured well with the vortex-based codes and the CFD solver. The observed trends in mean loading are in good agreement with previous wind tunnel tests of a scaled and stiff model of the tip extension. The CFD solution shows a highly three-dimensional flow at the very outboard part of the curved tip that leads to large changes of the angle of the resultant force with respect to the chord. Turbulent simulations using the BEM code and the vortex codes resulted in a good match with the measured standard deviation of the normal force, with some deviations of the BEM results due to the missing root vortex effect.

A parametric study on the hydrodynamics of tidal turbine blade deformation

The blades of axial-flow rotors, which are typically made of flexible composite materials, can experience significant deformations through their operation. However, the impact of these deformations on rotor hydrodynamics is not well understood. Blade deformations can be separated into three main components: flapwise, a thrust-driven deformation in the direction of the flow; edgewise, relatively small deformations in the plane of rotor motion and driven by rotor torque; and twist deformations that affect each blade’s angle of attack distribution and are generated by the torsional moment acting about the blade’s spanwise axis as well as by structural interactions. This work evaluates the hydrodynamic effects of the decoupled flapwise and twist deformations on a tidal turbine blade, using blade-resolved Computational Fluid Dynamics (CFD) simulations with the objective of identifying and quantifying the associated flow phenomena. The deformation cases were generated by scaling the static deformation shapes from a structural model of the turbine blade. The dataset used in this work, part of a larger research project, consists of 48 blade-resolved Reynolds-Averaged Navier–Stokes CFD steady-state simulations. It was found that twist deformation effects are significant, and can be adequately described using blade element theory. Flapwise deformations, on the other hand, produce different phenomena affecting the rotor loading and performance in ways that cannot be explained on a blade sectional basis through blade element theory. It is found that the hydrodynamic impact of flapwise deformations can be explained through two different mechanisms: a pressure drop on the suction side of the blade that generates inboard load augmentation, and an increase in near-tip radial flow effects that moves the onset of tip-loss effects, and shedding of thrust and power, further inboard. The studied cases show a significant impact of flapwise deformation on integrated power (up to a 20% drop), while the impact on integrated thrust remains limited, with variations between -4% to +2% in thrust coefficient observed depending on blade bending.

Optimization of Turbine Blade Aerodynamic Designs Using CFD and Neural Network Models

Optimization methods have been widely applied to the aerodynamic design of gas turbine blades. While applying optimization to high-fidelity computational fluid dynamics (CFD) simulations has proven capable of improving engineering design performance, a challenge has been overcoming the prolonged run-time due to the computationally expensive CFD runs. Reduced-order models and, more recently, machine learning methods have been increasingly used in gas turbine studies to predict performance metrics and operational characteristics, model turbulence, and optimize designs. The application of machine learning methods allows for utilizing existing knowledge and datasets from different sources, such as previous experiments, CFD, low-fidelity simulations, 1D or system-level studies. The present study investigates inserting a machine learning model that utilizes such data into a high-fidelity CFD driven optimization process, and hence effectively reduces the number of required evaluations of the CFD model. Artificial Neural Network (ANN) models were trained on data from over three thousand two-dimensional (2D) CFD analyses of turbine blade cross-sections. The trained ANN models were then used as surrogates in a nested optimization process alongside a full three-dimensional Navier–Stokes CFD simulation. The much lower evaluation cost of the ANN model allows for tens of thousands of design evaluations to guide the search of the best blade profiles to be used in the more expensive, high-fidelity CFD runs, improving the progress of the optimization while reducing the required computation time. It is estimated that the current workflow achieves a five-fold reduction in computational time in comparison to an optimization process that is based on three-dimensional (3D) CFD simulations alone. The methodology is demonstrated on the NASA/General Electric Energy Efficient Engine (E3) high pressure turbine blade and found Pareto front designs with improved blade efficiency and power over the baseline. Quantitative analysis of the optimization data reveals that some design parameters in the present study are more influential than others, such as the lean angle and tip scaling factor. Examining the optimized designs also provides insight into the physics, showing that the optimized designs have a lower amount of pressure drop near the trailing edge, but have an earlier onset of pressure drop on the suction side surface when compared to the baseline design, contributing to the observed improvements in efficiency and power.