Research Article  Open Access
Innovations in Body Force Modeling of Transonic Compressor Blade Rows
Abstract
Aeroengine fans and compressors increasingly operate subject to inlet distortion in the transonic flow regime. In this paper, innovations to loworder numerical modeling of fans and compressors via volumetric source terms (body forces) are presented. The approach builds upon past work to accommodate any axial fan/compressor geometry and ensures accurate work input and efficiency prediction across a range of flow coefficients. In particular, the efficiency dropoff near choke is captured. The model for a particular blade row is calibrated using data from singlepassage bladed computations. Compared to fullwheel unsteady computations which include the fan/compressor blades, the source term model approach can reduce computational cost by at least two orders of magnitude through a combination of reducing grid resolution and, critically, eliminating the need for a timeresolved approach. The approach is applied to NASA stage 67. For uniform flow, at 90% corrected speed and peakefficiency, the body force model is able to predict the totaltototal pressure rise coefficient of the stage to within 1.43% and the isentropic efficiency to within 0.03%. With a sector of reduced inlet total pressure, distortion transfer through the machine is wellcaptured and the associated efficiency penalty predicted with less than 2.7% error.
1. Introduction
In the preliminary design of compressors, maintaining low computational cost for numerical simulation of the blade aerodynamics is paramount. Traditionally, various 1D and 2D mathematical models have been employed. 1D calculations (meanline codes) yield the shortest computational time but with a reduction in both accuracy and available flow details. 2D models (streamline curvature methods, throughflow methods) provide a greater degree of accuracy; however, the flow is typically assumed to be axisymmetric.
In this paper, a 3D distributed source term modeling approach, commonly and henceforth referred to as a body force model, is developed. The approach can be applied to any axial compressor and captures efficiency dropoff near choke. This model is used to capture distortion transfer and efficiency penalties associated with inflow distortion. A body force model replaces physical blades in a numerical simulation with a volumetric source term field. Doing so considerably reduces the grid density required to capture key phenomena such as loading distribution, overall pressure rise, and work input when compared to bladed unsteady Reynoldsaveraged NavierStokes (URANS) computations. More significantly, the source term model is a “smearedout” circumferential average of the flow field; even in the case of nonaxisymmetric inflow, source terms are timeinvariant, eliminating the need for an unsteady calculation. In circumferentially nonuniform inflow, there is no frame of reference in which the flow is steady. However, the smearedout forces are not framedependent and do not vary with time if two conditions are met regarding the flow [1]: the flow nonuniformity must have a characteristic length scale significantly larger than the blade pitch and the local rotor reduced frequency must be sufficiently small, . The local rotor reduced frequency is expressed aswhere is the axial chord, is the axial velocity, and is the compressor rotational speed. Thus, the body force model is able to capture the response of a turbomachine to nonuniform inflow using only a steady calculation.
The distributed source term approach has been extensively used in the literature to study fan and compressor response to inlet flow distortion. A few recent examples which do not form an exhaustive list include the following: Defoe and Spakovszky investigated the generation and propagation of fan rotor shock noise due to inlet swirl distortion [2, 3], Thollet et al. used body force modeling to reduce the computational cost in studying intakefan interaction [4], Defoe et al. examined lowspeed boundarylayeringesting fan and compressor performance using a nonaxisymmetric throughflow method [1], and Guo and Hu observed multistage compressor performance in the presence of inlet flow distortion [5]. Most of these papers employed different approaches for obtaining the body force model used; what follows is a brief overview of approaches relevant to those adopted and built upon in the current paper.
The methodology was first introduced by Marble [6], replacing a physical blade row with an infinite number of infinitely thin blades. In studying shortwavelength stall inception and distortion transfer in multistage compressors, Gong was able to expand upon Marble’s approach by enabling the source terms to respond to local flow properties [7]. Volumetric source terms themselves are broken down in to a normal force per unit mass, , and parallel force per unit mass, . The normal component acts perpendicular to the relative streamlines, acting to reduce the deviation of the flow from the blade camber surface. The parallel component acts against the streamwise direction, thus generating viscous losses through the blade row.
Based on an adaptation from Gong’s distributed source term model, Peters developed a body force model to investigate fan inlet and nacelle design parameters for low pressure ratio fans with increased fan and inlet coupling [8]. Peters’ model includes a radial component in the normal force, accounting for blade lean and radial streamline shifts. The normal force is dependent on, among other variables, a calibrated blade loading coefficient, . The calibration process requires nonzero blade metal angle and nonzero flow deviation. Peters implemented an offdesign formulation to capture the variation in blade losses with operating condition [8]. Peters’ formulation uses a quadratic dependence on massaveraged relative Mach number at the blade row inlet, combined with a quadratic dependence on local relative velocity inherited from Gong. This formulation produces loss buckets typical of turbomachines.
More recently, an incompressible, inviscid, normal force model was developed by Hall, Greitzer, and Tan [9]. The normal force model is a function of local flow quantities and blade camber angle. Different from previous body force methods, it can be formulated without the need of calibration data (experimental or singlepassage RANS computations). This is helpful in preliminary design as full blade profiles need not be developed. The normal force model is developed such that blade rows of low solidity produce a normal force equivalent to thin airfoil lift theory . Conversely, in high solidity blades, the bladetoblade relative flow is unable to deviate from the blade metal angle.
In this paper, an innovative methodology for calibrating a body force model of a transonic compressor is demonstrated. A compressible adaptation of the incompressibleflow normal force model developed by Hall et al. [9] is developed in conjunction with an innovation to Peters’ loss model which enhances the ability to capture differences in efficiency falloff as the flow coefficient is either increased or decreased from the maximum efficiency point. The model is shown to accurately capture overall and local machine performance in both uniform and nonuniform inflow, on and offdesign. A reduction in computational cost of at least two orders of magnitude is obtained when compared against traditional fullannulus, timeresolved computations. To assess model performance, the compressor performance metrics and flow field details are compared against previously published experimental data as well as fullannulus, bladed computations.
2. Approach
2.1. Transonic Compressor Studied
Computations in this work are carried out on a transonic compressor, NASA stage 67 [11]. This singlestage axial compressor has a mean hubtotip radius ratio of 0.427, so significant radial flow streamline shifts are possible. The spinner nose is stationary; only a portion of the hub in the vicinity of the rotor rotates. All computations are carried out at 90% of the design corrected rotational speed, yielding a tip relative Mach number of 1.20 at peakefficiency. Additional design characteristics for the rotor at 90% corrected speed are given in Table 1 and for the stator in Table 2.


2.2. Body Force Modeling Approach
The normal force model in this work is based on the incompressible, inviscid turning force model developed by Hall et al. [9]:where is the flow deviation angle, is the number of blades, is the velocity magnitude in the rotational frame of reference moving with the blade, is the streamtube radius, and is the circumferential projection of the camber surface normal. To adapt this model for transonic flow, a compressibility correction field, , is added to the local flow deviation angle:where is a function of axial and radial location within each blade row’s swept volume. This field is calibrated to enforce the correct relative flow angle throughout the blade rows; the calibration is machinespecific and not meant to be an empirical constant for all turbomachines. The calibration procedure for imposes no restrictions on the type of turbomachine that can be modeled, as will be detailed in Section 2.4. Within a blade row, it is not possible to produce both the correct swirl velocity and relative flow angle simultaneously due to the absence of blade metal blockage. This phenomenon is captured in Figure 1; a mismatched chordwise axial velocity profile dictates that the body force model can only produce one of the two aforementioned properties accurately. The consequence of this is incorrect rotor work input since the trailing edge deviation is underpredicted. This follows directly from the Euler turbine equation:where is the total enthalpy and is the swirl velocity in a stationary frame of reference. For adiabatic flow, the total enthalpy change is governed purely by work done on the fluid. To produce the correct total enthalpy rise in the body force model, the change in swirl velocity must match actual machine behavior. This is accomplished through a twopart calibration: the compressibility correction factor tuning described earlier and effective camber reduction. Based on work by Patel et al. [12], the camber is reduced linearly from leading edge to trailing edge:where is the effective camber reduction of the rotor, is a blade row specific constant, is the axial coordinate, and is the local relative flow angle. To further simplify the model, the effects of blade lean are neglected. Whereas, in Hall’s model, the normal force is perpendicular to the relative streamwise direction in a plane shared by the blade normal and the local relative velocity , the force is instead taken to be normal to the relative streamwise direction, purely in the plane (). Within the normal force calculation, is replaced by to accomplish this. The implications of this simplification are later discussed. The normal force takes the final form
The formulation allows for the correct flow angle and thus correct loading at the blade leading edge, as well as correct swirl velocity and hence total enthalpy rise at the blade trailing edge. Ensuring the correct blade loading at the leading edge is critical for the model to be able to drive the upstream flow redistribution in nonuniform flow.
Peters developed a loss model for studying the effect of ultrashort nacelles on low pressure ratio fans [8], which included a quadratic dependence on bladeinlet massaveraged relative Mach number, . This formulation is used to capture increased entropy generation at offdesign flow conditions. The low pressure fans studied operate at flow coefficients far from choke where Peters’ model was sufficient to produce the desired isentropic efficiency versus flow coefficient characteristic. The flow coefficient is defined:where is the mass flow rate, is the massaveraged stagnation density, is the reference area at domain inlet, and is the midspan blade speed. Peters’ loss model is is calculated on a plane of constant axial coordinate at the leading edge of the blade; this is denoted as station 1 for the rotor and station 2 for the stator. The coefficients , , and are used to set the peakefficiency and the slope of the efficiency curve. In Peters’ work, is set to be at maximum efficiency. Finally, is the staggered blade spacing,where is the blade solidity. This staggered blade spacing is the distance between two discrete blades in the direction normal to the camber surface. Observation of (8) clearly illustrates that these coefficients set the value of peakefficiency or minimum value of . The deficiency of this model is its inability to capture the severe dropoff in performance as the machine approaches choke. The efficiency characteristic slope is significantly different at flow coefficients above and below peakefficiency, in particular for transonic machines where peakefficiency and choking values of are not far apart. To account for this, two innovations are implemented: is shifted away from peakefficiency point, and a piecewise loss function is established:where is a constant used to alter the efficiency at flow coefficients where . and serve as local loss minima within their respective multiplicative terms, allowing for two independent “shaping” variables. The values for , , , , and are calculated iteratively by comparing the resultant efficiency characteristic to the desired efficiency characteristic from highfidelity computations or experimental data. Each variable is adjusted independently until the desired output is obtained.
2.3. Calibration Data
For the calibration of the body force model, singlepassage bladed RANS calculations as well as previously published experimental data by Fidalgo et al. [10] are used. The geometry for the rotor blades was obtained from opensource rotor data [11] and the stator blade geometry was provided by Ewan J. Gunn (Whittle Laboratory, University of Cambridge, 1 JJ Thomson Avenue, Cambridge CB3 0DY, UK) and Vincente J. Fidalgo (Turbomachinery Laboratory, University of Notre Dame, 1165 Franklin Street, Suite 200, South Bend, IN 46601, USA). The duct is extended approximately five rotor axial chords upstream of the rotor and five stator axial chords downstream of the stator. The computations are carried out using ANSYS CFX 17.0 [13]. A combination of Pointwise [14] (upstream and downstream ducts) and TurboGrid [15] (rotor and stator blade passages) was used to create the grid. At the inlet, uniform total temperature, total pressure, and axial flow are specified; the static pressure is set at the domain outlet to yield the desired flow coefficient. The shearstress transport (SST) turbulence model [16] is used and a wall value less than 40 is maintained throughout the domain. This includes a maximum value of 10 along the upstream and downstream endwalls. along the stator blade is kept under 20 and for the rotor blade is largely kept under 20 other than in the tip gap region due to the locally high velocities. At 90% rotational speed, the flow is generally attached, such that the SST model behaves as the model, where is far below the recommended maximum of 300 [13]. Additional cell refinement is deemed unnecessary; increased grid density would minimally change the resultant viscous losses, whereas the importance of this work is body force modeling approach development. The details of the nearwall flow features are therefore not critical. The approach ensures that the efficiency and its variation with flow coefficient are well matched between body force and singlepassage computations, regardless of the accuracy of the singlepassage calibration data. Total grid count for this singlepassage calculation is 3.6 million cells. Initial grid sizing was based on a study performed by Fidalgo et al. [10] on the same NASA 67 machine and a grid independence study was carried out to assess the impact of further refinement. The final grid selected was the second of four grids tested. Grid independence was assumed once the rotor pressure ratio and isentropic efficiency changes were each less than 1%.
As discussed in Section 2.2, only peakefficiency singlepassage results are required for calibration of the normal force model. For the calibration of the parallel force coefficients, data is required over a range of flow coefficients. In this paper, experimental data provided by Fidalgo et al. [10] is used as a calibration source to reduce the number of singlepassage computations required. Within the rotor, the experimental total pressure ratio and isentropic efficiency at various flow coefficients are inputs to the parallel force coefficient calibration process. The body force loss model coefficients are calculated with the objective of matching these properties as closely as possible. In the stator, singlepassage computations are needed to tune the loss model as no statoralone experimental data is available. Thus, the stator body force model coefficients are calculated on the basis of matching the entropybased loss coefficient to the singlepassage results.
2.4. Calibration of the Body Force Model
A highlevel overview of the body force model calibration process is illustrated in Figure 2. The machine and operating speed of interest must first be selected, from which the singlepassage RANS speedline is calculated. Using both the singlepassage results and, in the case of the machine studied here, experimental data, the normal force model is calibrated and are determined). Following this, the parallel force model is calibrated simultaneously above and below by tuning each of the loss model coefficients.
Calculating the field requires calibrating against singlepassage RANS data. Initially, this field constant is set to zero. An iterative process is required to calculate the final field. The calibration constant iswhere is the compressibility calibration field obtained at each iteration and is the total number of iterations. Within each iteration, is calculated:where SP and BF represent single passage and body force, respectively. Each calculation is done concurrently in the rotor and stator. Calibration is complete once (11) produces a change less than 1% at each location in the entire field between iterations, for both blade rows. This iterative formulation enforces the body force to accurately predict the turning no matter the machine. The number of iterations required is expected to rise as flow compressibility and blade metal blockage increase.
The next step is to calculate the camber reduction, . To do so, the value is altered until the rotorexit massaveraged total temperature ratio is best matched with the reference values at peakefficiency flow coefficient.
The last step is to calibrate the parallel force model. To reduce the computational cost, educated guesses for the initial values of the model constants are important. Given that the flow coefficient at peakefficiency is known from experimental data, a simplification to (8) is made by combining the coefficients into :Using this simplified form of the viscous model, is determined using the same method employed by Peters [8] at a single flow coefficient (peakefficiency). From here, the viscous force field at peakefficiency remains invariant. Thus there is a constraint on the model constants such that the total value of at peakefficiency must be conserved. Previously used values for (0.0336) and (0.6321) in the work by Patel [12] served as initial guesses. is initially set to be equal to at peakefficiency as in Peters [8]. The revised parallel force formulation is used in the rotor, and Peters’ parallel force formulation is used unaltered in the stator. In the stator, the flow is further from choke, so that Peters’ loss model is sufficient. Without computational resources as a limiting factor, the revised parallel force model could be used in both blade rows but significant improvements in accuracy would not be expected.
2.5. Body Force Computational Domain
The body force grid in this work is a slice of the fullannulus machine, created using Pointwise. The axial extents of the domain are identical to the singlepassage grid. The rotor and stator blade regions encompass the swept volume of their respective blade rows. The structured grid is constructed on the meridional plane and subsequently revolved. The result is an axisymmetric grid, and uniform inflow calculation results are independent of the circumferential extent. The blade row models are implemented through the use of source terms in the momentum and energy equations within CFX. This is accomplished through a series of expressions and user functions. The forces are frameindependent, allowing for the use of a continuous grid with no mixing planes or interfaces of any kind. A meridional view of the body force domain and axial measurement locations is given in Figure 3. The boundary conditions are set in an identical manner to those in the singlepassage computations.
NASA stage 67 contains sharp hub radius changes in the rotor and stator. Within the stator, these area contractions coincide with the adverse pressure gradient produced by flow straightening and endwall boundary layer development. Due to these conditions, nonphysical severe flow separation was observed. The separation is not consistent with the results observed in the singlepassage model. To avoid this flow separation, the noslip condition is not enforced on endwalls in the body force computations. Within the rotor and stator blade regions, the parallel force model accounts for the endwall losses. Using slipwall conditions leads to the absence of boundary layer development upstream of the rotor. While it is desired to be able to capture this effect, the use of slipwalls only affects the flow in the bottom 5% span and upper 5% span. For the purpose of this study, the overall performance trends are still adequately captured.
Without the physical presence of blades in the body force model, the pressure surface to suction surface flow over the rotor blade tips is not captured. Having no source terms in the tip gap region in the distributed source term model was attempted, but, without physical blades, nonphysical backflow occurs. Thus, in this work, blade forces are extended into the tip gap, which for this machine is less than 0.4% of the blade span, yielding a tip region swept volume that is only 0.64% of the entire swept volume of the rotor. A consequence of this approach is that work is done on the flow in the tip gap, meaning the total enthalpy rise is overestimated, and the tip leakage flow details cannot be observed. However, the relatively small tip gap regions mean that the model is still able to predict the overall stage work to an acceptable degree, as will be shown.
A series of five grids were tested for the body force model, with increasing refinement in the axial, radial, and circumferential directions. Grid independence was assumed as the rotor work at peakefficiency changed significantly less than 1% between each of the successive grids. Rotor efficiency is not a concern in determining grid independence as the parallel force model is calibrated based on the selected grid. Instead, the additional constraint placed on grid selection is the tradeoff between refinement and computational cost. As grid density increases, work input by the rotor changes at a rate that is acceptable, <0.5% between all cases. The efficiency, however, is largely dependent on the grid discretization. The selected grid has a fullwheel cell count of cells, corresponding to the second level of the five grid densities tested. For comparison, this grid contains a cell count reduction of over 95% relative to the bladed RANS computational domain. The grid contains 60 cells in the spanwise direction, 50 cells in the axial direction for each blade row, and 0.71 cells/degree in the circumferential direction.
3. Results
3.1. Resultant Body Force Model
The field converged after 19 iterations and the resultant spatial distribution for NASA stage 67 is depicted in Figure 4. It is evident that a larger compressibility correction is required at increasing span values in the rotor and along the stator leading edge. This is expected, as in the rotor, at increasing radius, there is an associated increase in relative Mach number while in the stator the swirl and thus the Mach number are highest around the leading edge. Additionally, the compressibility correction is largest where the blade thickness gradient is largest. In the rotor, the correction factor field closely resembles the relative velocity field. In the stator near the hub, the correction is large to match the flow angles in the noslip endwall singlepassage calculations.
After seven iterations, the camber reduction constant converged to a value of . The resultant for the rotor is shown in Figure 5 and the effects of camber reduction are shown in Figure 6; the work input by the rotor is shifted downwards to match machine performance. Important to note is that the process of camber reduction does not appreciably alter the characteristic slope. In incompressible flow with axial inflow and neglecting deviation, the rate of change of work coefficient with flow coefficient isso that altering the trailing edge metal angle alters the slope of the rotor work characteristic. What is found, however, is that the camber reduction term produces what appears to be a vertical shift rather than a change in slope. This is due to two separate factors. Only a small range of flow coefficients are considered, all of which are far away from . Therefore, although the slope is changing, the apparent impact is more the updown shift of the characteristic over the flow coefficient region of interest. The machine studied in this paper is a transonic compressor and the experimental data shown in Figure 6 (open circles) clearly shows that the characteristic is nonlinear due to compressibility effects as the flow approaches choke; this effect is not captured in the body force model (solid line) because of the lack of metal blockage.
The parallel force coefficients found from calibration are shown in Table 3 and the resultant rotor isentropic efficiency characteristic is shown in Figure 7. At flow coefficients above peakefficiency, the model is able to capture the efficiency reduction as the flow approaches choke. At flow coefficients below peakefficiency, the reduction in isentropic efficiency is less severe, matching the performance predicted by the experimental results.

3.2. Example Applications of the Body Force Model
This section examines the body force model’s global and local performance metrics in both uniform and nonuniform flow, compared against fullannulus URANS computations.
3.2.1. Uniform Flow
To evaluate the model accuracy, a comparison is made against the uniform inflow singlepassage RANS calculation. At the peakefficiency flow coefficient, 0.505, the key performance metrics are shown in Table 4. The rotor isentropic efficiency is correct to within 0.03%. This degree of accuracy is expected, as the parallel force is calibrated to match the singlepassage isentropic efficiency. In the stator, the entropybased loss coefficient is used to assess accuracy:Here is the outlet static temperature, is the entropy change across the stator blade row, and represents the flow of kinetic energy per unit mass at the blade outlet. To calculate the massweighted average , streamlines are traced through the stator. Using these streamlines, the local value of is calculated for each streamline at a plane of constant axial coordinate near the stator outlet, plane 3 in Figure 3. The error in the loss coefficient is found to be 3.37% relative to singlepassage RANS results.

At flow coefficients away from peakefficiency, the model is compared against data from Fidalgo et al. [10]. The authors assessed NASA stage 67 with uniform inflow at 90% of design rotational speed both experimentally and using fullannulus URANS calculations. In Figure 8, the body force model results are plotted against these experimental and CFD results for the rotor isentropic efficiency, total temperature ratio, and total pressure ratio at various flow coefficients. The isentropic efficiency is well matched at flow coefficients both below and above design condition, even outperforming Fidalgo et al. [10] bladed URANS results at flow coefficients below peakefficiency. At the highest flow coefficients investigated, the body force model is unable to capture choking effects. This is expected, as blade metal blockage is not modeled. The total pressure ratio trend is also wellpredicted at all flow coefficients, aside from choke conditions, again due to the absence of blockage. There is a systematic increase in total pressure ratio relative to the previously published results (between 1.9% and 4.1%). This error stems from various sources: overprediction of rotor work input due to the artificial tip gap work, additional overprediction of rotor work input near choke, and a small error in the rotor isentropic efficiency, especially in its inability to match the severe dropoff near choking. The largest error is expected to be found in the total pressure versus flow coefficient characteristic as it is strongly dependent on the accuracy of both the normal and parallel force models.
In observing the spanwise mass flux profile at rotor trailing edge, shown in Figure 9, the effect of neglecting tip gap leakage is more evident. In the outer 8% span, the body force model does a poor job in reproducing the desired mass flux distribution. Additionally, slipwall conditions create another discrepancy in the results in the bottom 5% span, where the mass flux is overpredicted. At spanwise locations between the two endwalls, the model does an excellent job in capturing the mass flux of the singlepassage calculation. Finally, the spanwise work input from the rotor is shown in Figure 10. The error stems from the nature of the normal force modeling approach. Firstly, neglecting the blade lean in the normal force calculation may have led to the underprediction near the hub, the area with most significant lean. In addition, the relative flow angle is firstly constrained to match the singlepassage results, after which, camber reduction is performed. The method used in this study places emphasis on producing the correct stage work input. As a result, the massweighted average of total temperature from the body force model compares favorably to that of the singlepassage computations. However, the magnitude at each spanwise location leaves room for improvement.
3.2.2. Nonuniform Inflow: Low Total Pressure Sector
Recall that the intended application of the body force model approach is prediction of compressor response to nonuniform inflow. To assess the model in such conditions, a comparison was made against a nonuniform inflow study performed by Fidalgo et al. [10]. In their work, they subjected NASA stage 67 to a total pressure distortion. Figure 11 shows the far upstream inlet profile, at station 0. An operating speed of 90% of the design rotational speed is again used, with an inletcorrected flow coefficient based on massaveraged, total quantities of . The distorted region is set to a total pressure of 0.89 at the domain inlet. The experimental results at inlet involve significant variation from the ideal flow which exists in the CFD. This is due to limitations in how the distortion is produced experimentally given the discontinuous distribution desired. In this sense, the distortion studied is more severe than can be expected to be found in practical applications, such as boundary layer ingestion. For this reason, comparisons in this section will focus on the distributed source term model performance versus URANS results.
The rotor work input and trailing edge absolute swirl angle of the body force model are well matched with URANS results, as shown in the middle part of Figure 12. Both the rotor work input and flow angle follow the URANS trend line near perfectly, while the total pressure ratio trends are captured. The same behavior was observed in uniform flow as discussed earlier. Additionally, the circumferential smearingout of the flow field in the computations with the model means that the model is unable to capture the swirl variations due to stator upstream influence.
In addition to the circumferential flow characteristics, the body force model also captures blade trailing edge spanwise flow properties. As shown in Figure 13, the total pressure, total temperature, and absolute swirl angles are well matched at a circumferential coordinate location of , a spanwise location just outside of the distortion region; this location is chosen due to the availability of URANS data there. Two notable exceptions are near the hub and tip, for the same reasons discussed earlier (for uniform flow).
Overall performance metrics for the rotor are listed in Table 5. Comparing the results of the URANS calculation by Fidalgo et al. to the performance of the body force model, the change in rotor pressure ratio is perfectly matched to three significant figures (no change between the clean and distorted cases). The change in rotor isentropic efficiency, or the penalty associated with the distortion, is wellcaptured by the distributed source term model: 1.50% in Fidalgo et al. URANS and 1.46% in the body force model (2.67% error). Given that this is an exaggerated inlet distortion, this degree of accuracy emphasizes the effectiveness of the model for capturing distortion transfer and the response of the turbomachinery to the inflow distortion.

Within the stator blade passage, the qualitative flow response to the distortion is in excellent agreement with the results by Fidalgo et al. [10]. Midspan circumferential traverses of important flow properties downstream of the stator are shown in Figure 14. As expected the stator exhibits an absence of wakes in the source term model. Nevertheless, it is apparent that the flow angles would be in agreement if the blade wakes are filtered out. The total temperature is wellpredicted, meaning that the body force model is able to capture the distortion transfer through the stator. The total pressure profile at stator exit is nearly identical to the rotor exit, aside from a vertical shift down, due to excessive stator losses. Circumferential flow redistribution within the stator is minimal as illustrated by the similarity between rotorexit flow properties and stator exit flow properties. This is expected in a URANS computation, as it has been demonstrated that physical blades block circumferential flow redistribution over lengthscales greater than a blade pitch [17].
To quantify the overall stage performance of the body force model in nonuniform inflow, the stator trailing edge total temperature RMS error is 0.31%. The flow field at this location is a result of distortion transfer through both blade rows, showing that distortion transfer and blade response to distortion are wellcaptured by model. In comparing the total pressure at stator exit of the body force model to the URANS results, an RMS error of 3.1% exists. This is largely due to excessive loss forces in both the rotor and stator. To isolate the model from this loss force overestimation, an alternative RMS error is calculated. Artificially shifting the body force model data such that the mean value is equal to that of the URANS data, an RMS error of 0.57% is found. This supports the postulation that the model is able to capture distortion transfer as well as blade row losses; the accuracy of the calibration of the loss model is the limiting factor.
4. Summary and Conclusions
In this paper, a transonic body force modeling approach is developed by expanding upon existing flow turning and loss models. The model requires calibration of a compressibility correction factor that is both machine and operating speed dependent. Additionally, a camber reduction term is implemented to account for the absence of mass flow blockage. A twosided loss model is implemented to allow for an enhanced ability to capture choking effects. The resulting model obtained for NASA stage 67 is compared against singlepassage uniform inflow and fullannulus nonuniform inflow.
The key outcomes from this work are the following: a normal force compressibility correction in combination with camber modification is successful in producing a model which captures spanwise work input; the body force model is able to capture global performance metrics (rotor work and isentropic efficiency) in both uniform and nonuniform inflow; the body force model is able to capture localized flow properties and flow redistribution in both uniform and nonuniform inflow; and the modeling approach used is expected to be robust, since there are no limitations to blade geometry, flow regime, or incoming flow properties.
At peakefficiency, in uniform flow, the body force model is able to predict the rotor total pressure ratio within 1.43% as well as the isentropic efficiency within 0.03%. In the stator domain, the entropybased loss coefficient is found to have an error of 3.37%. Considering the rotor total temperature ratio, total pressure ratio, and isentropic efficiency away from peakefficiency, the largest error is found within the total pressure ratio. Even still, the overestimation in total pressure ratio is limited to 4.1% compared to URANS results. In the case of a severe nonuniform inflow, the stator trailing edge total temperature ratio at midspan has an RMS error of 0.31%, showing that the distortion transfer and blade work input are wellpredicted throughout the entire stage.
Future work will involve an automated, optimized scheme to calibrate the normal and parallel force models with greatly reduced user effort. Additionally, a mass flow blockage source term could be implemented to allow for the removal of the term and for more accurate blade loading. The model obtained is useful for studying the various effects of nonuniform flows on transonic compressors. Most importantly, the reduced computational cost when compared to fullannulus URANS allows for many different distortions to be studied. The intended application for this model is the study of efficiency scaling associated with varying distortion patterns and intensities.
Nomenclature
Reference area  
Staggered blade spacing  
Number of blades  
Lift coefficient per unit span  
Blade axial chord  
Volumetric source term per unit mass  
Local rotor reduced frequency  
Enthalpy  
Blade loading coefficient  
Combined loss coefficient  
Design loss coefficient  
Offdesign loss coefficient  
Corrected mass flow rate  
Relative Mach number  
Reference Mach number  
Camber surface normal  
Pressure  
Radial coordinate  
Entropy  
Temperature  
Velocity component in stationary frame  
Velocity magnitude in stationary frame  
Velocity magnitude in rotating frame  
Cartesian coordinates  
Absolute swirl angle  
Relative swirl angle  
Flow deviation  
Compressibility correction  
Entropybased loss coefficient  
Isentropic efficiency  
Circumferential coordinate  
Local blade camber angle  
Camber reduction constant  
Density  
Blade solidity  
Flow coefficient  
Work coefficient  
Corrected rotational speed. 
04  Axial measurement planes 
Inlet  
Isentropic  
Leading edge  
Normal to streamwise direction  
Parallel to streamwise direction  
Total/stagnation quantity  
Trailing edge  
Freestream. 
Massaveraged quantity. 
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that they have no conflicts of interest regarding this publication of this article.
Acknowledgments
This work is funded by the NSERC Discovery Grants program. It is based upon the work completed in Hill’s thesis [18]. Computational resources are provided by the facilities of the Shared Hierarchical Academic Research Computing Network (SHARCNET: www.sharcnet.ca) and Compute/Calcul Canada. The authors would like to thank Dr. Ewan Gunn and Dr. Vincente Fidalgo for the stage 67 stator model used in this work.
References
 J. J. Defoe, M. Etemadi, and D. K. Hall, “Fan performance scaling with inlet distortions,” Journal of Turbomachinery, vol. 140, no. 7, 11 pages, 2018. View at: Publisher Site  Google Scholar
 J. J. Defoe and Z. S. Spakovszky, “Effects of boundarylayer ingestion on the aeroacoustics of transonic fan rotors,” Journal of Turbomachinery, vol. 135, no. 3, 2013. View at: Google Scholar
 J. J. Defoe and Z. S. Spakovszky, “Shock Propagation and MPT Noise From a Transonic Rotor in Nonuniform Flow,” Journal of Turbomachinery, vol. 135, no. 1, 2012. View at: Google Scholar
 W. Thollet, G. Dufour, X. Carbonneau, and F. Blanc, “Bodyforce modeling for aerodynamic analysis of air intake  fan interactions,” International Journal of Numerical Methods for Heat & Fluid Flow, vol. 26, no. 7, pp. 2048–2065, 2016. View at: Publisher Site  Google Scholar
 J. Guo and J. Hu, “Development of body force model for steady inlet distortions in highspeed multistage compressor,” Proceedings of the Institution of Mechanical Engineers, Part G: Journal of Aerospace Engineering, vol. 231, no. 9, pp. 1650–1659, 2017. View at: Publisher Site  Google Scholar
 F. Marble, ThreeDimensional Flow in Turbomachines. High Speed Aerodynamics and Jet Propulsion, 10, 83–166, 1964.
 Y. Gong, C. S. Tan, K. A. Gordon, and E. M. Greitzer, “A computational model for shortwavelength stall inception and development in multistage compressors,” Journal of Turbomachinery, vol. 121, no. 4, pp. 726–734, 1999. View at: Publisher Site  Google Scholar
 A. Peters, Z. S. Spakovszky, W. K. Lord, and B. Rose, “Ultrashort nacelles for low fan pressure ratio propulsors,” Journal of Turbomachinery, vol. 137, no. 2, 2015. View at: Google Scholar
 D. K. Hall, E. M. Greitzer, and C. S. Tan, “Analysis of fan stage conceptual design attributes for boundary layer ingestion,” Journal of Turbomachinery, vol. 139, no. 7, 2017. View at: Google Scholar
 V. J. Fidalgo, C. A. Hall, and Y. Colin, “A Study of FanDistortion Interaction Within the NASA Rotor 67 Transonic Stage,” Journal of Turbomachinery, vol. 134, no. 5, 2012. View at: Google Scholar
 A. Strazisar, J. Wood, M. Hathaway, and K. Suder, Laser Anemometer Measurements in a Transonic AxialFlow Fan Rotor, 1989.
 K. Patel, C. Novak, and J. Defoe, “Novel numerical approach for generation and propagation of rotorstator interaction noise,” in Proceedings of the 22nd AIAA/CEAS Aeroacoustics Conference, 2016, France, June 2016. View at: Google Scholar
 ANSYS CFX Solver Theory Guide, Canonsburg, PA, USA, 2013.
 Pointwise User’s Manual, 4th edition, 2016.
 ANSYS TurboGrid User’s Guide, Canonsburg, PA, USA, 2013.
 F. Menter, “Zonal Two Equation kw Turbulence Models For Aerodynamic Flows,” in Proceedings of the 23rd Fluid Dynamics, Plasmadynamics, and Lasers Conference, Orlando, FL, U.S.A.. View at: Publisher Site  Google Scholar
 E. J. Gunn and C. A. Hall, “Aerodynamics of Boundary Layer Ingesting Fans,” in Proceedings of the ASME Turbo Expo 2014: Turbine Technical Conference and Exposition, p. V01AT01A024, Düsseldorf, Germany. View at: Publisher Site  Google Scholar
 D. Hill, Compressor Performance Scaling in the Presence of NonUniform Flow, University of Windsor, Canada, 2017.
Copyright
Copyright © 2018 David J. Hill and Jeffrey J. Defoe. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.