Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “thermal & statistical models”

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 469 records · Page 26

High-Reynolds Number Active Blowing Semi-Span Force Measurement System Development

Recent wind-tunnel tests at the NASA Langley Research Center National Transonic Facility utilized high-pressure bellows to route air to the model for evaluating aircraft circulation control. The introduction of these bellows within the Sidewall Model Support System significantly impacted the performance of the external sidewall mounted semi-span balance. As a result of this impact on the semi-span balance measurement performance, it became apparent that a new capability needed to be built into the National Transonic Facility s infrastructure to allow for performing pressure tare calibrations on the balance in order to properly characterize its performance under the influence of static bellows pressure tare loads and bellows thermal effects. The objective of this study was to design both mechanical calibration hardware and an experimental calibration design that can be employed at the facility in order to efficiently and precisely perform the necessary loadings in order to characterize the semi-span balance under the influence of multiple calibration factors (balance forces/moments and bellows pressure/temperature). Using statistical design of experiments, an experimental design was developed allowing for strategically characterizing the behavior of the semi-span balance for use in circulation control and propulsion-type flow control testing at the National Transonic Facility.

Lynn, Keith C.↗

Simulating cosmological evolution by quantum quench of an atomic Bose-Einstein condensate

In cosmological evolution, it is the homogeneous scalar field (inflaton) that drives the universe to expand isotropically and to generate standard model particles. However, to simulate cosmology, atomic gas research has focused on the dynamics of Bose-Einstein condensates (BEC) with continuously applied forces. In this paper we argue that a complementary approach needs also to be pursued; we thus consider the analog BEC experiments in a nondriven and naturally closed atomic system. We implement this using a BEC in an optical lattice which, after a quench, freely transitions from an unstable to a stable state. Furthermore, this dynamical evolution displays the counterpart “preheating,” “reheating,” and “thermalization” phases of cosmology. Importantly, our studies of these analog processes yield tractable analytic models. Additionally, of great utility to the cold atom community, such understanding elucidates the dynamics of nonadiabatic condensate preparation. Indeed, the dynamical processes discussed here are generic and in future cold atom reequilibration experiments it will be important to observe both the preheating stage, corresponding to a fragmented condensate, and the reheating stage, corresponding to a particle cloud.

74 ATOMIC AND MOLECULAR PHYSICS↗

Lattice QCD estimates of thermal photon production from the QGP

Thermal photons produced in heavy-ion collision experiments are an important observable for understanding quark-gluon plasma (QGP). The thermal photon rate from the QGP at a given temperature can be calculated from the spectral function of the vector current correlator. Extraction of the spectral function from the lattice correlator is known to be an ill-conditioned problem, as there is no unique solution for a spectral function for a given lattice correlator with statistical errors. The vector current correlator, on the other hand, receives a large ultraviolet contribution from the vacuum, which makes the extraction of the thermal photon rate difficult from this channel. We therefore consider the difference between the transverse and longitudinal part of the spectral function, only capturing the thermal contribution to the current correlator, simplifying the reconstruction significantly. The lattice correlator is calculated for light quarks in quenched QCD at T = 470 MeV ( ∼ 1.5 T c ), as well as in 2 + 1 flavor QCD at T = 220 MeV ( ∼ 1.2 T p c ) with m π = 320 MeV . In order to quantify the nonperturbative effects, the lattice correlator is compared with the corresponding NLO + LPM LO estimate of correlator. The reconstruction of the spectral function is performed in several different frameworks, ranging from physics-informed models of the spectral function to more general models in the Backus-Gilbert method and Gaussian process regression. We find that the resulting photon rates agree within errors. Published by the American Physical Society 2024

Astronomy & Astrophysics↗

Theory of the unmagnetized plasma.

The Vlasov mathematical model of a plasma, which has come to be thought more useful than any other in describing the dynamical behavior of the majority of plasmas of interest, is first examined. Macroscopic variables and moment equations; linear electrostatics solutions; plasma oscillations, ion acoustic waves, and linear instabilities are treated, as well as external fields, 'test' charges, and nonlinear Vlasov phenomena. Plasmas are statistically described, and attention is given to the kinetic theory of the stable, uniform plasma and the Balescu-Lenard equation; two-time ensemble averages and fluctuation spectra in stable plasmas; the kinetic theory of the unstable plasma; and ensembles of Vlasov plasmas. Some illustrative experiments are described. Four appendixes deal with the electrostatic approximation and transverse waves; solution of the linearized Vlasov equation in a magnetic field; estimates of correlation functions from thermal equilibrium; and equivalence of spatially uniform BBGKY and Klimontovich correlations.

Montgomery, D. C.↗

Thermal N+ in the Inner Magnetosphere

There has been much interest in O+ in the magnetosphere since it was first reported by Shelley et al. [1971]. However, relatively little has been done with N+ even though it is the second most abundant ion in the ionosphere at the F2 peak. What observations there are show that there is a large range in the magnitude of the ratio of the N+ density to the O+ density, that N+ is a significant ion in the ionosphere, and that the concentration of this ion vanes with season, time of day, solar cycle, latitude, and geomagnetic conditions. Most observational studies have not been used with a large database, nor have they been of a statistical nature. We use the DE1 RIMS observations to survey where N+ is found, in what concentrations relative to O+, and how this concentration varies with changes in the solar input, season, and time of day. We also gauge our understanding of N+ in the ionosphere and plasmasphere by comparing model (FLIP) results with the observations.

Craven, P. D.↗

Air separation and N 2 purification with Ba 0.15 Sr 0.85 FeO 3-$δ$ via a two-step thermochemical process

Thermochemical air separation to produce high-purity N 2 was demonstrated in a vertical tube reactor via a two-step reduction–oxidation cycle with an A-site substituted perovskite Ba 0.15 Sr 0.85 FeO 3–δ (BSF1585). BSF1585 particles were synthesized and characterized in terms of their chemical, morphological, and thermophysical properties. A thermodynamic cycle model and sensitivity analysis using computational heat and mass transfer models of the reactor were used to select the system operating parameters for a concentrating solar thermal-driven process. Thermal reduction up to 800 °C in air and temperature-swing air separation from 800 °C to minimum temperatures between 400 and 600 °C were performed in the reactor containing a 35 g packed bed of BSF1585. The reactor was characterized for dispersion, and air separation was characterized via mass spectrometry. Gas measurements indicated that the reactor produced N 2 with O 2 impurity concentrations as low as 0.02 % for > 30 min of operation. Additionally, a parametric study of air flow rates suggested that differences in observed and thermodynamically predicted O 2 impurities were due to imperfect gas transport in the bed. Temperature swing reduction/oxidation cycling experiments between 800 and 400 °C in air were conducted with no statistically significant degradation in N 2 purity over 50 cycles.

14 SOLAR ENERGY↗

Design and Construction of a High-Resolution Hodoscope for the GlueX Experiment with High-Statistics Analysis of the p0, ¿, and ¿1 Photoproduction Cross Sections from the RadPhi Experiment

Differential cross sections for forward-angle photoproduction of p0, ¿, and ¿ 1 pseudoscalar mesons were measured using data from the RadPhi experiment conducted in Hall B at Jef ferson Lab. RadPhi utilized a tagged bremsstrahlung photon beam incident on a stationary 9Be target, with a detector system configured to trigger on a recoil proton in coincidence with multiple neutral showers in the calorimeter. Events were reconstructed and subjected to kinematic constraints, with background suppressed via sideband subtraction guided by Monte Carlo modeling of background contributions. Cross sections were extracted over the photon energy range 4.4– 5.4 GeV and binned in invariant momentum transfer t, providing measurements from one of the first high-statistics experiments of forward ¿ and ¿1 pro duction from a nuclear target at these energies. Acceptance corrections were applied using a detailed GEANT-based simulation of the detector geometry and response. The resulting cross sections are consistent with 2020 CLAS results, when scaled by the number of protons in beryllium, and show broad agreement with other data and theoretical models. In parallel, a high-resolution photon tagger detector, the Tagger Microscope (TAGM), was designed, constructed, and commissioned for the GlueX experiment in Hall D at Jefferson Lab. The TAGM was developed to provide high-rate tagging capability in the coherent bremsstrahlung peak by detecting post-bremsstrahlung electrons across a one GeV range along the focal plane of the tagging spectrometer. The detector consists of a 5ˆ102 array of 2ˆ2 mm2 square BCF-20 plastic scintillating fibers thermally fused to BCF-98 light guide fibers optically coupled to silicon photomultipliers. These fibers are mounted in a precision machined framework enabling fine positional adjustments to maintain precise alignment with post-bremsstrahlung electron trajectories, while ensuring mechanical rigidity, thermal stability, optical isolation, minimal inactive area, and radiation shielding for electronics. The construction effort involved extensive testing of fiber quality, light transmission, thermal fusing, radiation hardness, and defect analysis using SEM and EDX techniques. Following its installation and commissioning, the TAGM became a critical component of the GlueX beamline, enabling high-rate tagging essential for studies of hybrid mesons and gluonic ex citations.

McIntyre, James [Univ. of Connecticut, Storrs, CT ↗

Effusion rate, length, and area relationships for some lava flows on Hawaii and Mount Etna with planetary implications

A model for the radiative cooling of thermally well mixed lava flows is presented and the relationship between effusion rate and length and area is analyzed. If radiative cooling is the prime mode of heat loss for a lava flow, one should expect to see a stronger correlation between the effusion rate and the plan area of the flow, than between effusion rate and just flow length. Different flows on a single volcano with differing initial temperatures, volatile content, and gross compositions should yield different areas for a given effusion rate. Likewise, a range of slopes for the relationship between effusion rate and flow area should result from comparisons between different volcanoes. As a test of these ideas, available data on the effusion rates, lengths, and areas of Hawaiian and Etnean flow is studied. It was found that: (1) the effusion rate/area correlation was statistically more significant than the correlation between effusion rate and length for four out of the five eruption episodes which met the necessary criteria of more than three individual flows with area, length, and effusion rate independently measured; (2) that there exists a minimum length and area for a given effusion rate, reflecting competition between overall characteristic proportionality between effusion rate and flow length, width, and area.

Pieri, D.↗

Effects of non-thermal termolecular reactions on wedge-induced oblique detonation waves

The shock-induced combustion ramjet (Shcramjet) based on oblique detonation waves (ODWs) is among the promising choices for hypersonic propulsion systems. An understanding of the ignition, propagation, and stability of ODWs is critical to harnessing their propulsive potential. In such high speed reacting flows, there is a high probability of occurrence of non-thermal reactions due to the presence of non-trivial amounts of highly reactive radicals including H, O and OH apart from O 2 as demonstrated recently [M. P. Burke, S. J. Klippenstein, Nat. Chem. 9 (2017) 1078–1082, Y. Tao, A. W. Jasper, Y. Georgievskii, S. J. Klippenstein, R. Sivaramakrishnan, Proc. Combust. Inst. 38 (2021) 515–522]. The present work focuses on examining the initiation, propagation and structure of oblique detonation waves in stoichiometric H 2 -air mixtures through numerical simulations with and without non-thermal reactivity on a two-dimensional adaptive grid. Non-thermal reactions were included in the macroscopic kinetic model as chemically termolecular reactions facilitated by the H + OH radical-radical recombination and the H + O 2 radical-molecule association reactions. Since, the non-thermal reactions result in a corresponding decrease in the reaction fluxes of the incipient recombination/association reactions, an additional simulation was performed by applying corrections to the respective incipient recombination/association rate constants using the methodology demonstrated by Tao et al. [Y. Tao, A. W. Jasper, Y. Georgievskii, S. J. Klippenstein, R. Sivaramakrishnan, Proc. Combust. Inst. 38 (2021) 515–522]. Results show that, under ODWE relevant conditions, non-thermal reactivity fundamentally alters the induction length, intensity as well as the structure of the ODW. Specifically, it is found that non-thermal reactivity leads to a noticeable reduction in initiation length and a simultaneous increase in instantaneous peak heat release rate and the degree of unsteadiness of the ODW. Finally, statistical analysis of key thermo-chemical variables is performed to elucidate the important species as well as reactions responsible for the observed variations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Evidence for the Thermal Sunyaev Zeldovich Effect Associated with Quasar Feedback

Using a radio-quiet subsample of the Sloan Digital Sky Survey spectroscopic quasar catalogue, spanning redshifts 0.5-3.5, we derive the mean millimetre and far-infrared quasar spectral energy distributions (SEDs) via a stacking analysis of Atacama Cosmology Telescope and Herschel-Spectral and Photometric Imaging REceiver data. We constrain the form of the far-infrared emission and find 3 sigma-4 sigma evidence for the thermal Sunyaev-Zel'dovich (SZ) effect, characteristic of a hot ionized gas component with thermal energy (6.2 plus or minus 1.7) × 10 (exp 60) erg. This amount of thermal energy is greater than expected assuming only hot gas in virial equilibrium with the dark matter haloes of (1-5) × 10(exp 12) h(exp −1) solar mass that these systems are expected to occupy, though the highest quasar mass estimates found in the literature could explain a large fraction of this energy. Our measurements are consistent with quasars depositing up to (14.5 +/- 3.3)tau (sub 8)(exp -1) per cent of their radiative energy into their circumgalactic environment if their typical period of quasar activity is tau(sub 8) x 108 yr. For high quasar host masses, approximately 10(exp 13) h(exp −1) solar mass, this percentage will be reduced. Furthermore, the uncertainty on this percentage is only statistical and additional systematic uncertainties enter at the 40 per cent level. The SEDs are dust dominated in all bands and we consider various models for dust emission. While sufficiently complex dust models can obviate the SZ effect, the SZ interpretation remains favoured at the 3 sigma-4 sigma level for most models.

Sunyaev-Zeldovich (SZ)↗

Material Interactions with Solar Wind Ion Environments

Solar wind composition is dominated by hydrogen (approx.96%) and helium (approx.3 to 4%) with a minor fraction (less than or equal to 1%) of heavy ions. Hydrogen (helium) ions impact spacecraft surfaces with energies from 0.5 to 5 keV (1.8 to 21 keV) due to variations in solar wind velocity from 300 km/s to 1000 km/sec with extremes of a few 10 s keV during periods of extremely high solar wind velocity exceeding 1000 km/sec. Mean impact energies are typically on the order of approximately 1 keV and 4 keV for hydrogen ions and helium ions, respectively. These energies are typically of the peak of the energy dependent light ion sputter yields for hydrogen and helium on many metals. In addition, light ions with kilovolt energies have been shown to produce blister (or exfoliation) damage to metal surfaces due to formation of high pressure gas bubbles within the materials when exposed to ion fluences on the order of 10(exp 16 to (10(exp 17 ions/sq cm. A number of spacecraft designs for current and future missions include gossamer polymer structures with thin metallic reflection coatings to shield instruments from the Sun or solar sail propulsion systems for use in a variety of locations in the inner solar system from 0.5 to 1 AU. In addition, there is interest in designing spacecraft for solar physics missions requiring operations as close to the Sun as 0.16 to 0.2 AU. Integrity of the metallic coatings is critical in many of these applications since degradation will result in modification of material thermal properties or exposure of polymers to solar UV photons which can compromise mission requirements. This paper will evaluate the relative contributions of sputtering and blister formation to material degradation in solar wind environments over a range of radial distances from the Sun to demonstrate where solar wind environments become important for materials selection. We will first review the physics and results from laboratory measurements of light ion sputtering, blistering, and exfoliation of metallic surfaces to establish the order of magnitude ion fluence required for significant surface damage. Solar wind ion fluence environments will then be evaluated due to variations in solar wind conditions as a function of solar cycle for varying distances from the Sun using models for radial variations in solar wind ion number density, temperature, and velocity to determine where sputtering and blistering is most likely to be an issue. Finally, ion fluence statistics for varying radial distances from the Sun will be shown to establish the mission duration and radial distances from the Sun where missions will encounter sufficient ion fluence to exhibit damage to metallic surfaces.

Minow, Joseph I.↗

DefectTrack: a deep learning-based multi-object tracking algorithm for quantitative defect analysis of in-situ TEM videos in real-time

Abstract In-situ irradiation transmission electron microscopy (TEM) offers unique insights into the millisecond-timescale post-cascade process, such as the lifetime and thermal stability of defect clusters, vital to the mechanistic understanding of irradiation damage in nuclear materials. Converting in-situ irradiation TEM video data into meaningful information on defect cluster dynamic properties (e.g., lifetime) has become the major technical bottleneck. Here, we present a solution called the DefectTrack , the first dedicated deep learning-based one-shot multi-object tracking (MOT) model capable of tracking cascade-induced defect clusters in in-situ TEM videos in real-time. DefectTrack has achieved a Multi-Object Tracking Accuracy (MOTA) of 66.43% and a Mostly Tracked (MT) of 67.81% on the test set, which are comparable to state-of-the-art MOT algorithms. We discuss the MOT framework, model selection, training, and evaluation strategies for in-situ TEM applications. Further, we compare the DefectTrack with four human experts in quantifying defect cluster lifetime distributions using statistical tests and discuss the relationship between the material science domain metrics and MOT metrics. Our statistical evaluations on the defect lifetime distribution suggest that the DefectTrack outperforms human experts in accuracy and speed.

42 ENGINEERING↗

The 2025 Evaluation of Experimental Thermonuclear Reaction Rates (ETR25)

This work describes the formalism for estimating thermonuclear reaction rates for astrophysical applications, emphasizing modern statistical approaches such as Monte Carlo sampling and Bayesian models. We discuss related topics including the calculation of resonance energies from nuclear Q values, indirect estimates of particle partial widths, and matching of reaction rates at elevated temperatures to statistical model results. We have evaluated available experimental data on cross sections, resonance energies and strengths, partial widths, lifetimes, spin-parities, and spectroscopic factors. Based on these results, we have estimated numerical values of 78 experimental charged-particle thermonuclear reaction rates for target nuclei in the A = 2–40 mass region, for temperatures ranging from 1 MK to 10 GK. For each reaction, three rate values are provided: low, median, and high, corresponding to the 16th, 50th, and 84th percentiles, respectively, of the cumulative reaction rate probability density distribution. Additionally, we present the factor uncertainty of each rate at each temperature grid point. These results enable users to sample the reaction rate probability density in nucleosynthesis calculations, facilitating uncertainty estimates of nuclidic abundances. The rates presented here refer to their laboratory values. For use in stellar model simulations, these values need to be corrected for the effects of thermal excitations of the interacting nuclei. For each reaction, we include graphs that illustrate the fractional contributions to the overall reaction rate along with the associated uncertainty. These visuals are designed to assist both stellar modelers and nuclear experimentalists by identifying the primary sources of rate uncertain=^texttx);ty at specific stellar temperatures. A graphical comparison with earlier Monte Carlo rates is also provided.

Nuclear astrophysics↗

Scenario Planning Management Actions to Restore Cold Water Stream Habitat: Comparing Mechanistic and Statistical Modeling Approaches

ABSTRACT Under the United States Clean Water Act, states are required to periodically assess state waters to determine compliance with water quality criteria (including temperature) and then to develop total maximum daily loads (TMDLs) for impaired waters as necessary to bring them into compliance. We compared the performance of mechanistic stream temperature models (HeatSource, QUAL2K, and QUAL2Kw) applied to the mainstem of three TMDL watersheds (Middle Fork John Day, OR; Wind River, WA; South Fork Nooksack, WA) with that of spatial stream network (SSN) models applied to the full watersheds and used these to evaluate the potential effectiveness of restoration strategies. SSN models performed well with slightly lesser accuracy (RMSE = 0.47–0.87) for mainstem predictions than mechanistic models (RMSE = 0.4) but provided additional benefits to inform management, including information on spatial and temporal heterogeneity of restoration effectiveness throughout the watershed. Of the four scenarios considered (restoration of riparian zones to potential natural vegetation, channel narrowing, increasing flow by restricting irrigation withdrawals, and combined applications), riparian zone restoration was consistently the most effective in reducing temperatures at the outlet, mainstem, and throughout the watersheds. Predicted restoration effectiveness for thermal regimes varied significantly both within and among watersheds. A focus on water quality criteria exceedance only at the watershed outlet or along the mainstem reach can obscure knowledge of restoration potential for fish habitat in tributaries and headwaters, potential for creation of thermal refuge areas along the mainstem critical for maintaining migration corridors, and thermal regime heterogeneity across space and time.

Fuller, M. R.↗

Superclustering with the Atacama Cosmology Telescope and Dark Energy Survey. II. Anisotropic Large-scale Coherence in Hot Gas, Galaxies, and Dark Matter

Statistics that capture the directional dependence of the baryon distribution in the cosmic web enable unique tests of cosmology and astrophysical feedback. We use constrained oriented stacking of thermal Sunyaev–Zel’dovich (tSZ) maps to measure the anisotropic distribution of hot gas 2.5–40 Mpc away from galaxy clusters embedded in massive filaments and superclusters. The cluster selection and orientation (at a scale of ∼15 Mpc) use Dark Energy Survey (DES) Year 3 data, while expanded tSZ maps from the Atacama Cosmology Telescope Data Release 6 enable a ∼3× more significant measurement of the extended gas compared to the technique’s proof-of-concept. Decomposing stacks into cosine multipoles of order m , we detect a dipole ( m = 1) and quadrupole ( m = 2) at 8σ–10σ, as well as evidence for m = 4 signal at up to 6σ, indicating sensitivity to late-time non-Gaussianity. We compare to Cardinal simulations with spherical gas models pasted onto dark matter halos. The fiducial tSZ data can discriminate between two models that deplete pressure differently in low-mass halos (mimicking astrophysical feedback), preferring higher average pressure in extended structures. However, uncertainty in the amount of cosmic infrared background contamination reduces the constraining power. Additionally, we apply the technique to DES galaxy density and weak lensing to study for the first time their oriented relationships with tSZ. In the tSZ-to-lensing relation, averaged on 7.5 Mpc (transverse) scales, we observe dependence on redshift but not shape or radial distance. Thus, on large scales, the superclustering of gas pressure, galaxies, and total matter is coherent in shape and extent.

79 ASTRONOMY AND ASTROPHYSICS↗

Vibrational Power Flow Analysis of Rods and Beams

A new method to model vibrational power flow and predict the resulting energy density levels in uniform rods and beams is investigated. This method models the flow of vibrational power in a manner analogous to the flow of thermal power in a heat conduction problem. The classical displacement solutions for harmonically excited, hysteretically damped rods and beams are used to derive expressions for the vibrational power flow and energy density in the rod and beam. Under certain conditions, the power flow in these two structural elements will be shown to be proportional to the energy density gradient. Using the relationship between power flow and energy density, an energy balance on differential control volumes in the rod and beam leads to a Poisson's equation which models the energy density distribution in the rod and beam. Coupling the energy density and power flow solutions for rods and beams is also discussed. It is shown that the resonant behavior of finite structures complicates the coupling of solutions, especially when the excitations are single frequency inputs. Two coupling formulations are discussed, the first based on the receptance method, and the second on the travelling wave approach used in Statistical Energy Analysis. The receptance method is the more computationally intensive but is capable of analyzing single frequency excitation cases. The traveling wave approach gives a good approximation of the frequency average of energy density and power flow in coupled systems, and thus, is an efficient technique for use with broadband frequency excitation.

Wohlever, James Christopher↗

Stellar fibril magnetic systems. I - Reduced energy state

The remarkable fibril structure of the magnetic fields at the surface of the sun (with fibrils compressed to 1,000-2,000 gauss) lies outside existing statistical theories of magnetohydrodynamic turbulence. The total energy of the fibril field is enhanced by a factor of more than 100 above the energy for the mean field in a continuum state. The magnetic energy density within a fibril is of the order of 100 times the local kinetic energy density, so that no simple application of equipartition principles is possible. It is pointed out that the total energy of the atmosphere (thermal + gravitational + magnetic) is reduced by the fibril state of the field by avoiding the magnetic inhibition of the convective overturning, suggesting that the formation of the observed intense fibril state may be in response to the associated energy reduction. Calculation of the minimum total energy of a polytropic atmosphere permeated by magnetic fibrils yields theoretical fibril fields of the order of 1-5 kilogauss when characteristics appropriate to the solar convective zone are introduced, in rough agreement with the actual fields of 1-2 kilogauss. The polytrope model, although crude, establishes that a large reduction in total energy is made possible by the fibril state.

Parker, E. N.↗

Combined Experimental and Numerical Simulations of Thermal Barrier Coated Turbine Blades Erosion

A combined experimental and computational study was conducted to investigate the erosion of thermal barrier coated (TBC) blade surfaces by alumina particles ingestion in a single stage turbine. In the experimental investigation, tests of particle surface interactions were performed in specially designed tunnels to determine the erosion rates and particle restitution characteristics under different impact conditions. The experimental results show that the erosion rates increase with increased impingement angle, impact velocity and temperature. In the computational simulations, an Euler-Lagrangian two stage approach is used in obtaining numerical solutions to the three-dimensional compressible Reynolds Averaged Navier-Stokes equations and the particles equations of motion in each blade passage reference frame. User defined functions (UDF) were developed to represent experimentally-based correlations for particle surface interaction models which were employed in the three-dimensional particle trajectory simulations to determine the particle rebound characteristics after each surface impact. The experimentally based erosion UDF model was used to predict the TBC erosion rates on the turbine blade surfaces based on the computed statistical data of the particles impact locations, velocities and angles relative to the blade surface. Computational results are presented for the predicted TBC blade erosion in a single stage commercial APU turbine, for a NASA designed automotive turbine, and for the NASA turbine scaled for modern rotorcraft operating conditions. The erosion patterns in the turbines are discussed for uniform particle ingestion and for particle ingestion concentrated in the inner and outer 5 percent of the stator blade span representing the flow cooling the combustor liner.

Hamed, Awate↗