Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “cascade”

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 1,171 records · Page 65

Development of an unstructured solution adaptive method for the quasi-three-dimensional Euler and Navier-Stokes equations

A general solution adaptive scheme-based on a remeshing technique is developed for solving the two-dimensional and quasi-three-dimensional Euler and Favre-averaged Navier-Stokes equations. The numerical scheme is formulated on an unstructured triangular mesh utilizing an edge-based pointer system which defines the edge connectivity of the mesh structure. Jameson's four-stage hybrid Runge-Kutta scheme is used to march the solution in time. The convergence rate is enhanced through the use of local time stepping and implicit residual averaging. As the solution evolves, the mesh is regenerated adaptively using flow field information. Mesh adaptation parameters are evaluated such that an estimated local numerical error is equally distributed over the whole domain. For inviscid flows, the present approach generates a complete unstructured triangular mesh using the advancing front method. For turbulent flows, the approach combines a local highly stretched structured triangular mesh in the boundary layer region with an unstructured mesh in the remaining regions to efficiently resolve the important flow features. One-equation and two-equation turbulence models are incorporated into the present unstructured approach. Results are presented for a wide range of flow problems including two-dimensional multi-element airfoils, two-dimensional cascades, and quasi-three-dimensional cascades. This approach is shown to gain flow resolution in the refined regions while achieving a great reduction in the computational effort and storage requirements since solution points are not wasted in regions where they are not required.

Jiang, Yi-Tsann↗

An investigation of small scales of turbulence in a boundary layer at high Reynolds numbers

The assumption that turbulence at large wave-numbers is isotropic and has universal spectral characteristics which are independent of the flow geometry, at least for high Reynolds numbers, has been a cornerstone of closure theories as well as of the most promising recent development in the effort to predict turbulent flows, viz. large eddy simulations. This hypothesis was first advanced by Kolmogorov based on the supposition that turbulent kinetic energy cascades down the scales (up the wave-numbers) of turbulence and that, if the number of these cascade steps is sufficiently large (i.e. the wave-number range is large), then the effects of anisotropies at the large scales are lost in the energy transfer process. Experimental attempts were repeatedly made to verify this fundamental assumption. However, Van Atta has recently suggested that an examination of the scalar and velocity gradient fields is necessary to definitively verify this hypothesis or prove it to be unfounded. Of course, this must be carried out in a flow with a sufficiently high Reynolds number to provide the necessary separation of scales in order unambiguously to provide the possibility of local isotropy at large wave-numbers. An opportunity to use our 12-sensor hot-wire probe to address this issue directly was made available at the 80'x120' wind tunnel at the NASA Ames Research Center, which is normally used for full-scale aircraft tests. An initial report on this high Reynolds number experiment and progress toward its evaluation is presented.

Wallace, James M.↗

Optical system storage design with diffractive optical elements

Optical data storage systems are gaining widespread acceptance due to their high areal density and the ability to remove the high capacity hard disk from the system. In magneto-optical read-write systems, a small rotation of the polarization state in the return signal from the MO media is the signal which must be sensed. A typical arrangement used for detecting these signals and correcting for errors in tracking and focusing on the disk is illustrated. The components required to achieve these functions are listed. The assembly and alignment of this complex system has a direct impact on cost, and also affects the size, weight, and corresponding data access rates. As a result, integrating these optical components and improving packaging techniques is an active area of research and development. Most designs of binary optic elements have been concerned with optimizing grating efficiency. However, rigorous coupled wave models for vector field diffraction from grating surfaces can be extended to determine the phase and polarization state of the diffracted field, and the design of polarization components. A typical grating geometry and the phase and polarization angles associated with the incident and diffracted fields are shown. In our current stage of work, we are examining system configurations which cascade several polarization functions on a single substrate. In this design, the beam returning from the MO disk illuminates a cascaded grating element which first couples light into the substrate, then introduces a quarter wave retardation, then a polarization rotation, and finally separates s- and p-polarized fields through a polarization beam splitter. The input coupler and polarization beam splitter are formed in volume gratings, and the two intermediate elements are zero-order elements.

Kostuk, Raymond K.↗

Scale disparity and spectral transfer in anisotropic numerical turbulence

To study the effect of cancellations within long-range interactions on local isotropy at the small scales, we calculate explicitly the degree of cancellation in distant interactions in the simulations of Yeung & Brasseur and Yeung, Brasseur & Wang using the single scale disparity parameter 's' developed by Zhou. In the simulations, initially isotropic simulated turbulence was subjected to coherent anisotropic forcing at the large scales and the smallest scales were found to become anisotropic as a consequence of direct large-small scale couplings. We find that the marginally distant interactions in the simulation do not cancel out under summation and that the development of small-scale anisotropy is indeed a direct consequence of the distant triadic group, as argued by Yeung, et. al. A reduction of anisotropy at later times occurs as a result of the isotropizing influences of more local energy-cascading triadic interactions. Nevertheless, the local-to-nonlocal triadic group persists as an isotropizing influence at later times. We find that, whereas long-range interactions, in general, contribute little to net energy transfer into or out of a high wavenumber shell k, the anisotropic transfer of component energy within the shell increases with increasing scale separations. These results are consistent with results by Zhou, and Brasseur & Wei, and suggest that the anisotropizing influences of long range interactions should persist to higher Reynolds numbers. The residual effect of the forced distant group in this low-Reynolds number simulation is found to be forward cascading, on average.

Zhou, YE↗

Homogeneous turbulence subjected to mean flow with elliptic streamlines

Direct numerical simulations are performed for homogeneous turbulence with a mean flow having elliptic streamlines. This flow combines the effects of rotation and strain on the turbulence. Qualitative comparisons are made with linear theory for cases with high Rossby number. The nonlinear transfer process is monitored using a generalized skewness. In general, rotation turns off the nonlinear cascade; however, for moderate ellipticities and rotation rates the nonlinear cascade is turned off and then reestablished. Turbulence statistics of interest in turbulence modeling are calculated, including full Reynolds stress budgets.

Blaisdell, G. A.↗

N2 Lyman-Birge-Hopfield dayglow from ATLAS 1

Spectral measurements of the far ultraviolet dayglow were made from the ATLAS 1 shuttle mission using an intensified charge coupled device (CCD) imaging spectrometer array. The instrument imaged relatively large (177 A) segments of the spectrum simultaneously and had the capability to image altitude in the second dimension of the two-dimensional focal plane detector when the entrance slit of the instrument was oriented perpendicular to the limb of the Earth. The far ultraviolet (FUV) channel of the instrument used a CsTe photocathode rather than the more conventionally used and more solar blind CsI photocathode material. As a result, the spectra were also measured with good sensitivity at FUV wavelengths longer than 1500 A. In this paper, data are selected of specific interest to a study of the N2 Lyman-Birge-Hopfield system using two shuttle attitudes: In the first the field of view was held at a fixed tangent ray height of approximately 140 km while scanning in wavelength, providing repeated FUV spectra over a range of latitudes and allowing comparison with our thermospheric airglow model. In the second observing sequence the line of sight of the instrument was scanned down through the atmosphere but covered about 75% of the full wavelength range. This sequence allows comparison of height dependence,including O2 absorption effects, with the thermospheric model. In all cases, relatively good agreement is obtained with the model (run strictly in a predictive mode) in comparisions in which the line of sight must be intergrated along a path that spans a considerable range of conditions (height, local time, latitude, solar zenith angle). In the case of the data sets examined in this study the vibrational population distributions show mixed results. In one case the distribution agrees well with a theoretical model based on direct electron impact excitation without a significant cascade contribution. A second case indicated a lower population in nu' = 3 and a higher population in nu' = 6 than would be expected from direct impact excitation alone. This supports previous indications that conditions resulting in the cascade mechanism may not always be operative.

Torr, Marsha R.↗

Nonlinear evolution of a large-amplitude circularly polarized Alfven wave: Low beta

The nature of turbulent cascades arising from the parametric instabilities of a monochromatic field-aligned large-amplitude circularly polarized Alfven wave is investigated via direct numerical simulation for the case of low plasma Beta and no wave dispersion. The magnetohydrodynamic code permits nonlinear couplings in the parallel direction to the ambient magnetic field and one perpendicular direction. Compressibility is included in the form of a polytropic equation of state. Anisotropic turbulent cascades, similar to those found in early incompressible two-dimensional simulations, occur after nonlinear saturation of the parallel propagating decay instability. The turbulent spectrum can be divided into three regimes: the lowest wave numbers are dominated by lower sideband remnants of the parametric process, intermediate wave numbers display nearly incompressible dynamics, and the highest wave numbers are dominated by acoustic turbulence.

Ghosh, S.↗

Intermittent turbulence in solar wind from the south polar hole

The magnetic fields measured by the Ulysses spacecraft are used to study solar wind turbulence in the fast solar wind from the south polar hole. The spacecraft was at about 46 deg south latitude and 3.9 AU. For a magnetic field with a Gaussian distribution the power spectrum (second-order structure function) is sufficient to completely characterize the turbulence. However, the actual distribution is non-Gaussian so that the effects of intermittency must be taken into account. The observed spectral exponents include effects of intermittency and cannot be directly compared with the standard second-order spectral theories such as the Kolmogorov and Kraichnan theories. To permit a better comparison of the observations with the theoretical models, we study the structure characteristics of the data. We find the exponents of the second-order structure functions (power spectra) and the higher-order normalized structure functions for the components of the magnetic fields. We show that these sets of exponents can be approximately described by two basic numbers: the spectral exponent and the intermittency exponent. The intermittency exponent characterizes correlation properties of the energy cascade from large to small scales. Before comparing the observations to the theoretically expected values, a reduction must be made to the observed spectral exponent. The amount of the reduction depends on both the intermittency exponent and the model of the energy cascade assumed in the turbulence theory. We reduce the measured spectral indices according to a simple model for Alfven turbulence that is described here. We then compare our reduced spectral indices with second-order spectral theory. The reduced spectral indices for the period range of 1 min to about a half hour are remarkably constant and in good agreement with the value of 3/2. Thus our treatment is self-consistent. Our tentative conclusion is that the high-frequency turbulence appears to agree with the model of random-phased Alfven waves. This tentative conclusion must be tested by further theoretical and observational work.

Ruzmaikin, A. A.↗

Multifractal scaling of the kinetic energy flux in solar wind turbulence

The geometrical and scaling properties of the energy flux of the turbulent kinetic energy in the solar wind have been studied. By present experimental technology in solar wind measurements, we cannot directly measure the real volumetric dissipation rate, epsilon(t), but are constrained to represent it by surrogating the energy flux near the dissipation range at the proton gyro scales. There is evidence for the multifractal nature of the so defined dissipation field epsilon(t), a result derived from the scaling exponents of its statistical q-th order moments. The related generalized dimension D(q) has been determined and reveals that the dissipation field has a multifractal structure. which is not compatible with a scale-invariant cascade. The associated multifractal spectrum f(alpha) has been estimated for the first time for MHD turbulence in the solar wind. Its features resemble those obtained for turbulent fluids and other nonlinear multifractal systems. The generalized dimension D(q) can, for turbulence in high-speed streams, be fitted well by the functional dependence of the p-model with a comparatively large parameter, p = 0.87. indicating a strongly intermittent multifractal energy cascade. The experimental value for D(p)/3, if used in the scaling exponent s(p) of the velocity structure function, gives an exponent that can describe some of the observations. The scaling exponent mu of the auto correlation function of epsilon(t) has also been directly evaluated. It has the value of 0.37. Finally. the mean dissipation rate was determined, which could be used in solar wind heating models.

Marsch, E.↗

Analysis of Subgrid Eddy Diffusivity in Turbulence Using Direct Numerical Simulation Databases

Direct numerical simulation data in isotropic turbulence at Taylor-scale Reynolds number 90 are used to compute subgrid transfer and spectral eddy diffusivities for passive scalars at Schmidt numbers 1 and 1/8. The results are consistent with spectral transfer processes dominated by a local forward cascade in wavenumber space which is a consequence of moderately nonlocal interactions in the scalar fields. Data decompositions show that resolvable-subgrid interactions make (in general) substantially stronger contributions than fully-subgrid interactions to the eddy diffusivity, which has a non-negligible inverse cascade, or backscatter, component. The eddy diffusivity has a strong cusp at wavenumbers close to the cutoff value separating resolvable and subgrid ranges. At low wavenumbers the eddy diffusivity is generally small, but there the details depend on Schmidt number. The behavior of the subgrid turbulent Schmidt number (the ratio of effective eddy viscosity to eddy diffusivity) is observed to depend strongly on Schmidt number. In interpreting the results at a given spectral cutoff it is important to note that the spectrum of less diffusive scalars (of higher Schmidt number) has more high-wavenumber content.

Yeung, P. K.↗

Optimization of Air-Breathing Engine Concept

The design optimization of air-breathing propulsion engine concepts has been accomplished by soft-coupling the NASA Engine Performance Program (NEPP) analyzer with the NASA Lewis multidisciplinary optimization tool COMETBOARDS. Engine problems, with their associated design variables and constraints, were cast as nonlinear optimization problems with thrust as the merit function. Because of the large number of mission points in the flight envelope, the diversity of constraint types, and the overall distortion of the design space; the most reliable optimization algorithm available in COMETBOARDS, when used by itself, could not produce satisfactory, feasible, optimum solutions. However, COMETBOARDS' unique features-which include a cascade strategy, variable and constraint formulations, and scaling devised especially for difficult multidisciplinary applications-successfully optimized the performance of subsonic and supersonic engine concepts. Even when started from different design points, the combined COMETBOARDS and NEPP results converged to the same global optimum solution. This reliable and robust design tool eliminates manual intervention in the design of air-breathing propulsion engines and eases the cycle analysis procedures. It is also much easier to use than other codes, which is an added benefit. This paper describes COMETBOARDS and its cascade strategy and illustrates the capabilities of the combined design tool through the optimization of a high-bypass- turbofan wave-rotor-topped subsonic engine and a mixed-flow-turbofan supersonic engine.

Patnaik, Surya N.↗

Application of the kappa-omega Turbulence Model to Quasi-Three-Dimensional Turbomachinery Flows

Many computational fluid dynamics codes for turbomachinery use the Baldwin-Lomax (B-L) turbulence model. It is easy to implement in two dimensions and works well for predicting overall turbomachinery performance. However, it is awkward to implement in three dimensions, often has difficulty finding the length scale, has a crude transition model, and neglects freestream turbulence, surface roughness, and mass injection. The kappa-omega model developed by Wilcox is an appealing alternative for several reasons. First, it is the only two-equation model that can be integrated to the wall without requiring damping functions or the distance to the wall, and hence, should behave well numerically. Second, the effects of freestream turbulence, surface roughness, and mass injection are easily included in the model. Finally, transition can be simulated using the low Reynolds number version of the model. Menter applied the kappa-model to external flows and showed very good results for flows with adverse pressure gradients. Liu and Zheng described their implementation of the kappa-model in a cascade code that included an area change term to account for endwall convergence. They validated the model for a flat plate, and compared the B-L and kappa-models to measured surface pressures for a low-pressure turbine cascade. Since they did not use the low Reynolds number version of the model, their results showed problems resulting from early transition. In this Note the low Reynolds number kappa-model was incorporated in the author's quasi-three-dimensional turbomachinery analysis code. The code includes the effects of rotation, radius change, and stream-surface thickness variation, and also includes the B-L turbulence model. The kappa-omega model was implemented using many of Menter's recommendations and an implicit approximate-factorization scheme described by Baldwin and Barth. The model was tested for a transonic compressor with rotation and variable stream-surface radius and height, and for a transonic turbine vane with transition and heat transfer. Results were compared to the B-L model and to experimental data.

Chima, Rodrick V.↗

Low Speed, 2-D Rotor/Stator Active Noise Control at the Source Demonstration

Wake/blade-row interaction noise produced by the Annular Cascade Facility at Purdue University has been modeled using the LINFLO analysis. Actuator displacements needed for complete cancellation of the propagating acoustic response modes have been determined, along with the associated actuator power requirements. As an alternative, weighted least squares minimization of the total far-field sound power using individual actuators has also been examined. Attempts were made to translate the two-dimensional aerodynamic results into three-dimensional actuator requirements. The results lie near the limit of present actuator technology. In order to investigate the concept of noise control at the source for active rotor/stator noise control at the source, various techniques for embedding miniature actuators into vanes were examined. Numerous miniature speaker arrangements were tested and analyzed to determine their suitability as actuators for a demonstration test in the Annular Cascade Facility at Purdue. The best candidates demonstrated marginal performance. An alternative concept to using vane mounted speakers as control actuators was developed and tested. The concept uses compression drivers which are mounted externally to the stator vanes. Each compression driver is connected via a tube to an air cavity in the stator vane, from which the driver signal radiates into the working section of the experimental rig. The actual locations and dimensions of the actuators were used as input parameters for a LINFLO computational analysis of the actuator displacements required for complete cancellation of tones in the Purdue experimental rig. The actuators were designed and an arrangement determined which is compatible with the Purdue experimental rig and instrumentation. Experimental tests indicate that the actuators are capable of producing equivalent displacements greater than the requirements predicted by the LINFLO analysis. The acoustic output of the actuators was also found to be unaffected by the presence of air flow representative of the Purdue experimental rig. A test of the active noise control at the source concept for rotor/stator active noise control was demonstrated. This 2-D test demonstrated conclusively the simultaneous reduction of two acoustic modes. Reductions of over 10 dB were obtained over a wide operating range.

Simonich, John C.↗

Prediction of Unsteady Transitional Layers in Turbomachinery Using Navier Stokes Equations

The objective of the research reported in this presentation is to develop computational techniques for the prediction of unsteady transitional flows associated with the rotor stator interaction in turbomachinery. Three low-Reynolds number turbulence models are incorporated in two unsteady Navier-Stokes codes (one is pressure based and the other is time marching with Runge-Kutta time stepping) and evaluated for accuracy in predicting the onset and the end of unsteady transitional patches due to wake passing. The best model is then used for modification and improvement for the leading edge effect. An existing steady Navier-Stokes code was modified to include pseudo-time stepping, which provided acceleration from 5 to 25 times that of the original code. A systematic validation procedure was implemented to assess the effects of the grid, artificial dissipation, physical, and the pseudo-time step for an accurate prediction of transitional flows resulting from the rotor-stator interaction. The ability of the Navier-Stokes code to predict the unsteady transitional flow on a turbomachinery blade is demonstrated. The unsteady pressure and velocity fields are in good agreement with the experimental data and the prediction from the Euler/boundary layer approach. The numerical solver was able to capture all zones (wake induced transitional strip, wake induced turbulent strip, calmed region, etc.) associated with wake induced transition in a compressor cascade. Another significant step is the assessment of k-epsilon turbulence models, including the leading edge modifications. Best results were obtained from the FLB model. The LB model predicted earlier inception of the transition and shorter transition length. Modification of the k-epsilon model was found to be essential for an accurate prediction of the unsteady transitional flow in a compressor cascade. The CH model failed to predict the unsteady transitional flow. Predicted boundary layer was turbulent from the leading edge, even with the modification of the k-epsilon model near the stagnation point. A comparison of the instantaneous shape factor, the skin friction coefficient, and the momentum thickness indicates that the Navier-Stokes predictions are reasonably good. Interaction between the upstream wake and stator wake results in shedding of unsteady vortices from the trailing edge and increased dissipation in the stator wake and, as a consequence, increased rate of decay of the stator wake. The procedure developed at Penn State should aid the designers of turbomachinery in allowing for the unsteady transitional flows due to rotor-stator interaction.

Lakshminarayana, B.↗

The Development of Kolmogoroff-Like Power Spectra in the Expanding Solar Wind

Power spectra of the solar wind fluctuations consistently exhibit a -5/3 power-law slope consistent with the idea that the medium is undergoing a turbulent cascade as seen in ordinary fluids. This is surprising both because the radial streams and the magnetic field threading the plasma will induce anisotropies and because the expansion of the wind will tend to lead to the suppression of nonlinear cascades. These conditions violate the assumptions used by Kolmogoroff to derive the -5/3 law. We have studied this issue using a compressible IMHD code in spherical coordinates and have shown that a -5/3 spectrum results from a broad-band flat-spectrum input condition that is sheared and distorted by a current sheet as the wind expands. We determine spectra from time series taken at selected points in the domain as is done with observational spacecraft data. The spectra are very like those we have seen in nonexpanding runs and exhibit evolution and compressive characteristics very similar to those seen in observations. We will report on these results in addition to a new set of runs intended to constrain the necessary and sufficient conditions for the spectra to have this form. The simulation also allows us to examine the anisotropy for the spectra to attempt to determine why the result of an isotropic magnetofluid is obtained in a highly anisotropic situation.

Goldstein, M. L.↗

Clouds in ECMWF's 30 KM Resolution Global Atmospheric Forecast Model (TL639)

Global models of the general circulation of the atmosphere resolve a wide range of length scales, and in particular cloud structures extend from planetary scales to the smallest scales resolvable, now down to 30 km in state-of-the-art models. Even the highest resolution models do not resolve small-scale cloud phenomena seen, for example, in Landsat and other high-resolution satellite images of clouds. Unresolved small-scale disturbances often grow into larger ones through non-linear processes that transfer energy upscale. Understanding upscale cascades is of crucial importance in predicting current weather, and in parameterizing cloud-radiative processes that control long term climate. Several movie animations provide examples of the temporal and spatial variation of cloud fields produced in 4-day runs of the forecast model at the European Centre for Medium-Range Weather Forecasts (ECMWF) in Reading, England, at particular times and locations of simultaneous measurement field campaigns. model resolution is approximately 30 km horizontally (triangular truncation TL639) with 31 vertical levels from surface to stratosphere. Timestep of the model is about 10 minutes, but animation frames are 3 hours apart, at timesteps when the radiation is computed. The animations were prepared from an archive of several 4-day runs at the highest available model resolution, and archived at ECMWF. Cloud, wind and temperature fields in an approximately 1000 km X 1000 km box were retrieved from the archive, then approximately 60 Mb Vis5d files were prepared with the help of Graeme Kelly of ECMWF, and were compressed into MPEG files each less than 3 Mb. We discuss the interaction of clouds and radiation in the model, and compare the variability of cloud liquid as a function of scale to that seen in cloud observations made in intensive field campaigns. Comparison of high-resolution global runs to cloud-resolving models, and to lower resolution climate models is leading to better understanding of the upscale cascade and suggesting new cloud-radiation parameterizations for climate models.

Cahalan, R. F.↗

Energy Calibration of the Scintillating Optical Fiber Calorimeter Chamber (SOFCAL)

The Scintillating Optical Fiber Calorimeter (SOFCAL) detector is designed to make direct measures of the primary cosmic ray spectrum from -200 GeV/amu - 20 TeV/amu. The primary particles are resolved into groups according to their charge (p, He, CNO, Medium Z, Heavy Z) using both active and passive components integrated into the detector. The principal part of SOFCAL is a thin ionization calorimeter that measures the electromagnetic cascades that result from these energetic particles interacting in the detector. The calorimeter is divided into two sections: a thin passive emulsion/x-ray film calorimeter, and a fiber calorimeter that uses crossing layers of small scintillating optical fibers to sample the energy deposition of the cascades. The energy determination is made by fitting the fiber data to transition curves generated by Monte Carlo simulations. The fiber data must first be calibrated using the electron counts from the emulsion plates in the calorimeter for a small number of events. The technique and results of this calibration will be presented together with samples of the data from a balloon flight.

Christl, M. C.↗

Effect of Film-Hole Shape on Turbine Blade Film Cooling Performance

The detailed heat transfer coefficient and film cooling effectiveness distributions as well as tile detailed coolant jet temperature profiles on the suction side of a gas turbine blade A,ere measured using a transient liquid crystal image method and a traversing cold wire and a traversing thermocouple probe, respectively. The blade has only one row of film holes near the gill hole portion on the suction side of the blade. The hole geometries studied include standard cylindrical holes and holes with diffuser shaped exit portion (i.e. fanshaped holes and laidback fanshaped holes). Tests were performed on a five-blade linear cascade in a low-speed wind tunnel. The mainstream Reynolds number based on cascade exit velocity was 5.3 x 10(exp 5). Upstream unsteady wakes were simulated using a spoke-wheel type wake generator. The wake Strouhal number was kept at 0 or 0.1. Coolant blowing ratio was varied from 0.4 to 1.2. Results show that both expanded holes have significantly improved thermal protection over the surface downstream of the ejection location, particularly at high blowing ratios. However, the expanded hole injections induce earlier boundary layer transition to turbulence and enhance heat transfer coefficients at the latter part of the blade suction surface. In general, the unsteady wake tends to reduce film cooling effectiveness.

Han, J. C.↗