Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “large deviations”

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,207 records · Page 67

The Completed SDSS-IV extended Baryon Oscillation Spectroscopic Survey: Cosmological Implications from two Decades of Spectroscopic Surveys at the Apache Point observatory

We present the cosmological implications from final measurements of clustering using galaxies, quasars, and Ly$\alpha$ forests from the completed Sloan Digital Sky Survey (SDSS) lineage of experiments in large-scale structure. These experiments, composed of data from SDSS, SDSS-II, BOSS, and eBOSS, offer independent measurements of baryon acoustic oscillation (BAO) measurements of angular-diameter distances and Hubble distances relative to the sound horizon, $r_d$, from eight different samples and six measurements of the growth rate parameter, $f\sigma_8$, from redshift-space distortions (RSD). This composite sample is the most constraining of its kind and allows us to perform a comprehensive assessment of the cosmological model after two decades of dedicated spectroscopic observation. We show that the BAO data alone are able to rule out dark-energy-free models at more than eight standard deviations in an extension to the flat, $\Lambda$CDM model that allows for curvature. When combined with Planck Cosmic Microwave Background (CMB) measurements of temperature and polarization the BAO data provide nearly an order of magnitude improvement on curvature constraints. The RSD measurements indicate a growth rate that is consistent with predictions from Planck primary data and with General Relativity. When combining the results of SDSS BAO and RSD with external data, all multiple-parameter extensions remain consistent with a $\Lambda$CDM model. Regardless of cosmological model, the precision on $\Omega_\Lambda$, $H_0$, and $\sigma_8$, remains at roughly 1\%, showing changes of less than 0.6\% in the central values between models. The inverse distance ladder measurement under a o$w_0w_a$CDM yields $H_0= 68.20 \pm 0.81 \, \rm km\, s^{-1} Mpc^{-1}$, remaining in tension with several direct determination methods. (abridged)

79 ASTRONOMY AND ASTROPHYSICS↗

Evidence for WW/WZ vector boson scattering in the decay channel ℓνqq produced in association with two jets in proton-proton collisions at s = 13 TeV

Evidence is reported for electroweak (EW) vector boson scattering in the decay channel $\ell\nu$qq of two weak vector bosons WV (V = W or Z), produced in association with two parton jets. The search uses a data set of proton-proton collisions at 13 TeV collected with the CMS detector during 2016-2018 with an integrated luminosity of 138 fb$^{-1}$. Events are selected requiring one lepton (electron or muon), moderate missing transverse momentum, two jets with a large pseudorapidity separation and a large dijet invariant mass, and a signature consistent with the hadronic decay of a W/Z boson. The cross section is computed in a fiducial phase space defined at parton level requiring all parton transverse momenta $p_\mathrm{T}$$\gt$ 10 GeV and at least one pair of outgoing partons with invariant mass $m_\mathrm{qq}$$\gt$ 100 GeV. The measured and expected EW WV production cross sections are 1.90 $^{+0.53}_{-0.46}$ pb and 2.23 $^{+0.08}_{-0.11}$ (scale) $\pm$ 0.05 (PDF) pb, respectively, where PDF is the parton distribution function. The observed EW signal strength is $m_\mathrm{EW}$ = 0.85 $\pm$ 0.12 (stat) $^{+0.19}_{-0.17}$ (syst), corresponding to a signal significance of 4.4 standard deviations with 5.1 expected, and it is measured keeping the quantum chromodynamics (QCD) associated diboson production fixed to the standard model prediction. This is the first evidence of vector boson scattering in the $\ell\nu$qq decay channel at LHC. The simultaneous measurement of the EW and QCD associated diboson production agrees with the standard model prediction.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Double-Crucible Vertical Bridgman Technique for Stoichiometry-Controlled Chalcogenide Crystal Growth

Precise stoichiometry control in single-crystal growth is essential for both technological applications and fundamental research. However, conventional growth methods often face challenges such as non-stoichiometry, compositional gradients, and phase impurities, particularly in non-congruent melting systems. Even in congruent melting systems like Bi₂Se₃, deviations from the ideal stoichiometric composition can lead to significant property degradation, such as excessive bulk conductivity, which limits its topological applications. In this study, we introduce the double-crucible vertical Bridgman (DCVB) method, a novel approach that enhances stoichiometry control through the combined use of continuous source material feeding, traveling-solvent growth, and liquid encapsulation, which suppresses volatile element loss under high pressure. Using Bi₂Se₃ as a model system, we demonstrate that crystals grown via DCVB exhibit enhanced stoichiometric control, significantly reducing defect density and achieving much lower carrier concentrations compared to those produced by conventional Bridgman techniques. Moreover, the continuous feeding of source material enables the growth of large crystals. As a result, this approach presents a promising strategy for synthesizing high-quality, large-scale crystals, particularly for metal chalcogenides and pnictides that exhibit challenging non-congruent melting behaviors.

Crystal structure↗

The growth of density perturbations in the last ~10 billion years from tomographic large-scale structure data

In order to investigate the origin of the ongoing tension between the amplitude of matter fluctuations measured by weak lensing experiments at low redshifts and the value inferred from the cosmic microwave background anisotropies, we reconstruct the evolution of this amplitude from z ~ 2 using existing large-scale structure data. To do so, we decouple the linear growth of density inhomogeneities from the background expansion, and constrain its redshift dependence making use of a combination of 6 different data sets, including cosmic shear, galaxy clustering and CMB lensing. We analyze these data under a consistent harmonic-space angular power spectrum-based pipeline. We show that current data constrain the amplitude of fluctuations mostly in the range 0.2 < z < 0.7, where it is lower than predicted by . This difference is mostly driven by current cosmic shear data, although the growth histories reconstructed from different data combinations are consistent with each other, and we find no evidence of systematic deviations in any particular experiment. In spite of the tension with , the data are well-described by the model, albeit with a lower value of S 8 ≡ σ 8 (Ω m /0.3) 0.5 . As part of our analysis, we find constraints on this parameter of S 8 = 0.7781 ± 0.0094 (68% confidence level), reaching almost percent-level errors comparable with CMB measurements, and 3.4σ away from the value found by .

79 ASTRONOMY AND ASTROPHYSICS↗

A machine learning approach to photometric metallicities of giant stars

ABSTRACT Despite the advances provided by large-scale photometric surveys, stellar features – such as metallicity – generally remain limited to spectroscopic observations often of bright, nearby low-extinction stars. To rectify this, we present a neural network approach for estimating the metallicities and distances of red giant stars with 8-band photometry and parallaxes from Gaia EDR3 and the 2MASS and WISE surveys. The algorithm accounts for uncertainties in the predictions arising from the range of possible outputs at each input and from the range of models compatible with the training set (through drop-out). A two-stage procedure is adopted where an initial network to estimate photoastrometric parallaxes is trained using a large sample of noisy parallax data from Gaia EDR3 and then a secondary network is trained using spectroscopic metallicities from the APOGEE and LAMOST surveys and an augmented feature space utilizing the first-stage parallax estimates. The algorithm produces metallicity predictions with an average uncertainty of $\pm 0.19\, \mathrm{dex}$. The methodology is applied to stars within the Galactic bar/bulge with particular focus on a sample of 1.69 million objects with Gaia radial velocities. We demonstrate the use and validity of our approach by inspecting both spatial and kinematic gradients with metallicity in the Galactic bar/bulge recovering previous results on the vertical metallicity gradient (−0.528 ± 0.002 dex kpc−1) and the vertex deviation of the bar ($-21.29\pm 2.74\, \mathrm{deg}$).

79 ASTRONOMY AND ASTROPHYSICS↗

In situ ligand restraints from quantum-mechanical methods

In macromolecular crystallographic structure refinement, ligands present challenges for the generation of geometric restraints due to their large chemical variability, their possible novel nature and their specific interaction with the binding pocket of the protein. Quantum-mechanical approaches are useful for providing accurate ligand geometries, but can be plagued by the number of minima in flexible molecules. In an effort to avoid these issues, the Quantum Mechanical Restraints (QMR) procedure optimizes the ligand geometry in situ , thus accounting for the influence of the macromolecule on the local energy minima of the ligand. The optimized ligand geometry is used to generate target values for geometric restraints during the crystallographic refinement. As demonstrated using a sample of >2330 ligand instances in >1700 protein–ligand models, QMR restraints generally result in lower deviations from the target stereochemistry compared with conventionally generated restraints. In particular, the QMR approach provides accurate torsion restraints for ligands and other entities.

59 BASIC BIOLOGICAL SCIENCES↗

SCUC-DER Integration Report: Integrating Distributed Energy Resources (DER) Using Advanced Unit Commitment Models and DER Aggregation Methodologies

Distributed energy resources (DERs) are continuing to grow due to regulatory, policy, and market shifts, and it needs to be ensured that small DERs are given a level playing field with traditional resources. Legacy market processes such as unit-commitment problems were designed for a power grid consisting largely of centralized power plants. In contrast, DERs consist of many small devices with distinct operating characteristics that may or may not be connected at the transmission interconnection points, limiting their visibility to the independent system operators (ISOs) who operate wholesale electricity markets. This report details the development and initial results from a simulation platform that integrates state-of-the-art security-constrained unit commitment software, detailed feeder models, and a DER aggregator model to quantify potential DER integration issues. Quantitative results to date illustrate potential infeasible scheduling solutions from SCUC when the DERA includes aggregations of energy-limited energy storage resources. Likewise, if aggregations are not penalized for dispatch deviations, they may have incentives to deviate from the SCUC-determined resource schedules. Assumptions about the amount of aggregated demand response resources (DRRs) and the ability of DERAs to follow profit incentives have an important impact – DRRs have significant flexibility and can typically feasibly meet their SCUC schedule, but on the other hand, their profit incentives can cause unscheduled increases in load before and after DRR dispatch. Computation time results on the SCUC solver and simulation platform only show a modest increase in SCUC solution time as the number of DERAs is increased, but results to date only reflect the RTS-GMLC test system; results may show more significant solver slowdown in larger transmission systems. The largest contribution to simulation time is attributed to the DERA offer generation method, which is suggested for future improvements.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Automated Fiber Placement Through Thickness Defect Stacking Optimization

In its 2022 commercial market outlook, Boeing forecasted an 80% increase in the global fleet through 2041 compared to 2019 pre-pandemic levels. This sharp rise in demand will drive pressure onto airframe manufacturers to ramp up production and find more efficient ways to design and manufacture airplanes. Complicating this challenge is the industry’s recent transformation from traditional metal-based airframes towards hybrid composite-metal aircraft. While composites have been used in aviation for decades, aircraft manufacturers are still struggling to design and manufacture quality parts at a high rate. Automated Fiber Placement (AFP) is one of the main manufacturing techniques used to produce large-scale composite parts. After a design has been created, a manufacturing strategy has to be developed based on the working material, part geometry, and machine capabilities. This process planning stage is essential to the AFP workflow and currently requires a high level of manual input from an experienced process planner. In an effort to automate and optimize this stage, the Computer Aided Process Planning (CAPP) module was developed. CAPP assists process planners in identifying optimal starting point location and layup strategy for each ply of a laminate. This Ply-Level Optimization (PLO) phase operates on the quantification of ply quality through predictable geometry-based defects such as gaps, overlaps, angle deviation, and steering. As you move from PLO to Laminate-Level Optimization (LLO) the design space grows exponentially, emphasizing the need for automated optimization. The work presented in this thesis expands CAPP’s functionality by comparing the planned fiber paths through the thickness of the laminate to mitigate stacked area defects and achieve an optimal laminate-level manufacturing strategy. Within CAPP, predicted gap and overlap defects are imported from Vericut Composites Programming (VCP) and then discretized to streamline the through-thickness comparison. Two objective functions are used to score different combinations of ply layup strategies based on defect stacking both globally and locally. Four combinatorial optimization algorithms were coupled with these objective functions to investigate the laminate-level manufacturing strategy design space and converge on the optimal plan. These algorithms were evaluated based on accuracy and efficiency through virtual testing on a complex tool surface. A separate LLO approach was developed to achieve near-optimal laminates in significantly less time. The end result is a software package which greatly reduces the required input from process planners, shortening the design-build cycle time and improving part quality.

AFP↗

Large-scale spectral structure with a gap in the stably stratified atmosphere

The analysis of the large-scale structure of turbulence in a stratified atmosphere is important for the prediction and modelling of turbulence. The Monin-Obukhov similarity consideration is well suited for the modelling of small-scale turbulence. But the large-scale turbulence possesses new features and parameters, as connected with the appearance of a spectral gap in the stably stratified atmosphere, for the understanding of which an analytical investigation is needed. By a group-kinetic method, the anisotropic form of the transfer function and its governing eddy viscosity are derived. On this basis, the spectral laws k to the -5/3rd, k to the -1st, k to the -3rd are obtained for the subranges of inertia, shear, and buoyancy, respectively. A gap connects the last two subranges by the variable anisotropy, as a result of the fall of the spectrum of the vertical velocity fluctuations. The variable anisotropy causes the deviations from the spectral laws k to the -1st and k to the -3rd. An experiment supports the analytical prediction.

Larsen, S. E.↗

A millimeter-wave radiometer for detecting microbursts

This paper describes a millimeter-wave radiometer for the detection of wind shear from airborne platforms or at airport terminals. This proposed instrument will operate near the group of atmospheric oxygen absorptions centered near 60 GHz, which it will use to sense temperature from a distance. The instrument will use two channels to provide two different temperature measurements, providing the basis for solution of two equations in two unknowns, which are range to the wind shear plume and its temperature. A third channel will measure ambient atmospheric temperature. Depending on the temperature difference between the wind-shear plume and ambient, the standard deviation of range measurement accuracy is expected to be about 1 km at 5 km range, while the temperature measurement standard deviation will be about one-fourth the temperature difference between plume and ambient at this range. The instrument is expected to perform usefully at ranges up to 10 km, giving adequate warning of the presence of wind shear even for high performance jet aircraft. Other atmospheric hazards which might be detected by this radiometer include aircraft wakes and vortices, clear-air turbulence, and wind rotors, although the latter two phenomena would be detected by an airborne version of the instrument. A separate radiometer channel will be provided in the proposed instrument to detect aircraft wakes and vortices based on perturbation of the spectrum of microscopic atmospheric temperature fluctuations caused by the passage of large aircraft.

Mcmillan, Robert↗

The importance of properties in modeling

Casting and welding of superalloys, stainless steel and titanium alloys are processes which can be improved through modeling of heat flow, fluid flow, residual stress development, and microstructural evolution. These simulations require inputs of thermophysical data, some of which involves the partially or totally liquid state. In particular, these processes involve melting, flow in the liquid, and solidification. Modeling of such processes can lead to an improved understanding of defects such as shrinkage, inclusions, cracks, incomplete filling (or penetration), macrosegregation, improper grain structure, and deviations from dimensional specifications. Effective modeling can shorten process development time and improve quality. An approach to these problems is to develop efficient models; validate through correlations with thermal, distortion, and microstructural data; run parametric studies; extract knowledge based rules; and apply to adaptive closed loop control systems. With the appropriate pre- and post-processing, such analyses can be made 'user friendly'. This would include graphical user interfaces as well as realistic images and color maps. In such form, these models can be used for sensitivity analyses, which are useful in defining appropriate sensors and in the development of control strategies. Such modeling can be done at several levels, e.g., the MARO level, modeling large scale phenomena such as heat and fluid flow or material deformation; the MICRO level, modeling the development of dendrites, grains or precipitates; or at the NANO level, modeling point defects, dislocations, stacking faults, etc. There are many computational issues associated with these simulations, e.g. computational efficiency and accuracy. In addition, there are many materials issues, not the least of which is the availability of accurate high temperature thermophysical data for complex alloys. This would include latent heat of fusion, temperature dependent heat capacity and thermal conductivity (for liquid and solid), viscosity, surface tension, thermal expansion, mechanical properties, etc. Preliminary data is frequently gathered from the literature; however, this is often not available for modern alloys. If additional data are required, measurements can be used; however, these are costly, time consuming and can be erroneous due to a lack of testing standards or impure materials. Microstructural predictors can be extracted from thermal information, e.g. cooling rate and thermal gradient; the prediction of microstructure is dependent on solidus and liquidus temperature, mushy zone permeability, the solidification curve, volume changes, phase transformations, alloying effects (such as surface tension or viscosity), mold/metal reactions, metal/environment reactions, etc. Defect maps may be needed to predict the onset of shrinkage, hot cracking or 'freckling'. Constants may be needed for stress relaxation, dendrite coarsening, vaporization, etc. Visualization was used as a tool to better comprehend complex data sets associated with the analysis of directional solidification (including crystal growth) and welding. Examples include not only isotherms, but also cooling rate, growth rate and thermal gradient. The latter two are not single valued scalars, but rather time and space dependent vector fields. Efficient models were developed for both casting and welding to predict heat flow and the relationship to dendrite and grain growth. These codes include many of the non-linear effects, e.g. radiation, which dominate these processes. The home-built FDM code(s) were designed to be useful not only to the scientist, but also to the process engineer. Special output can be requested to compare directly to experimental data. Visualization procedures were developed to visualize critical results, e.g. fusion zone width at the surface opposite that where the arc is applied ('penetration'). Both elaborate and simplified distortion analyses were carried out. It is clear that extensive mechanical property data are critical in order to accurately predict residual stress patterns. A scheme is currently being developed to integrate these modeling tools into a set of control algorithms; however, the success of this approach is critically dependent on the availability of accurate high temperature thermophysical data.

Giamei, A. F.↗

Compositional and phase dependence of elastic modulus of crystalline and amorphous Hf 1- x Zr x O 2 thin films

The elastic moduli of amorphous and crystalline atomic layer-deposited Hf 1-x Zr x O 2 (HZO, x = 0, 0.31, 0.46, 0.79, 1) films prepared with TaN electrodes on silicon substrates were investigated using picosecond acoustic measurements. The moduli of the amorphous films were observed to increase between 211 ± 6 GPa for pure HfO 2 and 302 ± 9 GPa for pure ZrO 2 . In the crystalline films, it was found that the moduli increased upon increasing the zirconium composition from 248 ± 6 GPa for monoclinic HfO 2 to 267 ± 9 GPa for tetragonal ZrO 2 . Positive deviations from this increase were observed for the Hf 0.69 Zr 0.31 O 2 and Hf 0.54 Zr 0.46 O 2 compositions, which were measured to have moduli of 264 ± 8 GPa and 274 ± 8 GPa, respectively. These two compositions contained the largest fractions of the ferroelectric orthorhombic phase, as assessed from polarization and diffraction data. The biaxial stress states of the crystalline films were characterized through sin 2 (ψ) x-ray diffraction analysis. The in-plane stresses were all found to be tensile and observed to increase with the increasing zirconium composition, between 2.54 ± 0.6 GPa for pure HfO 2 and 5.22 ± 0.5 GPa for pure ZrO 2 . The stresses are consistent with large thermal expansion mismatches between the HZO films and silicon substrates. Furthermore, these results demonstrate a device-scale means to quantify biaxial stress for investigation on its effect on the ferroelectric properties of hafnia-based materials.

36 MATERIALS SCIENCE↗

A Multi-Objective Bayesian Optimized Human Assessed Multi-Target Generated Spectral Recommender System for Rapid Pareto Discoveries of Material Properties

Optimization for different tasks like material characterization, synthesis, and functional properties for desired applications over multi-dimensional control parameter and function spaces need a rapid strategic search through active learning. However, in all cases prior to optimization, the target material properties are assumed known and fixed, which mostly deviates from real-world scenarios in material synthesis. This can be critical for running expensive experiments on new materials, when the experimental results are fuzzy for any scientific outcomes due to improper target setting, ultimately wasting time and cost. The failure rate and cost are even higher over exploring on multi-target space, where we want to learn the pareto among multiple properties, to jointly optimize during material synthesis for desired applications. To address the challenge, here we introduce the human-operator attempt flexibility in the active learning based automated experiment framework, with generating multiple human assessed targets through a voting-based recommender system during real-time microscope measurements over the large material image space, sequentially learn/update multiple desired targets through a weighting system, and adaptively search in multiple material properties functional space for non-dominated pareto discoveries to maximize the custom structural similarity based acquisition function. We term this a multi-objective Bayesian optimized human assessed multi-target generated spectral recommender systems (MOBO-HAM-SRS). The approach has been demonstrated to peizoresponse force spectroscopy of a ferroelectric thin film, exploring with different kernels and acquisition functions. This work shows an advancement towards human-AI collaborated automated experiments, steering optimization trajectories through human overpowering AI at the early stage when uncertainty is high and AI overpowering human at the later stage with rapid exploration towards optimal goal, following human-assessed multiple targets properties.

Biswas, Arpan↗

A Survey of Large-Scale Variations in Thermospheric Oxygen Column Density with Magnetic Activity as Inferred from Observations of the FUV Dayglow

Brightness of the terrestrial far-ultraviolet (FUV) dayglow is dominated by the 130.4-nm emission of neutral atomic oxygen, OI, and variations in the brightness observed from altitudes high above the emitting region reflect variations in thermospheric oxygen density. This paper summarizes the results of an initial survey of the Dynamics Explorer 1 observations of the FUV dayglow through a presentation of 13 representative events selected to demonstrate the spatial extent and short-term temporal stability of the brightness perturbations. The emphasis here is on the morning sector of local time and the polar cap for observations obtained in the time interval from September 23, 1981, through January 19, 1982. An analytic expression is derived for the average response of the FUV photometer to the dayglow during periods of high-latitude magnetic quiescence. The remaining observations in this time interval are then analyzed for their deviations from the established quiet time values. Deviations of -40% to +30% are found following intervals of increased magnetic activity. The most significant decreases (-30% to -40%) are observed equatorward of the instantaneous aurora] oval only after sustained periods (approx. 6 hours) of intense magnetic activity (average AE greater than approx. 700 nT). Decreases extend equatorward from the aurora to geographic latitudes as low as approx. 30 deg N. Decreases of lesser magnitude that do not extend as far equatorward are associated with sustained periods of more moderate activity in which the average value of AE is smaller (approx. 300-400 nT). Also, the spatial extent and magnitude of the decreases in the morning sector appear greater when the IMF B(sub y) component is positive. In both cases, decreases are readily observed within the polar cap. Localized enhancements of +20% to +30% occur much less frequently and are detected at the middle latitudes, well equator-ward of the auroral oval.

Nicholas, A. C.↗

Least-squares/parabolized Navier-Stokes procedure for optimizing hypersonic wind tunnel nozzles

A new procedure is demonstrated for optimizing hypersonic wind-tunnel-nozzle contours. The procedure couples a CFD computer code to an optimization algorithm, and is applied to both conical and contoured hypersonic nozzles for the purpose of determining an optimal set of parameters to describe the surface geometry. A design-objective function is specified based on the deviation from the desired test-section flow-field conditions. The objective function is minimized by optimizing the parameters used to describe the nozzle contour based on the solution to a nonlinear least-squares problem. The effect of the changes in the nozzle wall parameters are evaluated by computing the nozzle flow using the parabolized Navier-Stokes equations. The advantage of the new procedure is that it directly takes into account the displacement effect of the boundary layer on the wall contour. The new procedure provides a method for optimizing hypersonic nozzles of high Mach numbers which have been designed by classical procedures, but are shown to produce poor flow quality due to the large boundary layers present in the test section. The procedure is demonstrated by finding the optimum design parameters for a Mach 10 conical nozzle and a Mach 6 and a Mach 15 contoured nozzle.

Korte, John J.↗

A GCM case study on the maintenance of short-term subtropical wind maxima in the summer hemisphere during SOP-1, FGGE

In order to substantiate earlier findings of pronounced short-term deviations in Southern Hemisphere tropical convection during the first Special Observing Period (SOP-1) of FGGE, the hypothesis that the observed summertime subtropical wind maxima were strengthened by episodes of enhanced tropical divergence is numerically tested using a fourth-order GCM. FGGE observational data for January 6-20, 1979 suggested that strong convection over the South Pacific helped to maintain a nearby subtropical jet. The dominant term in the zonal momentum budget was the Coriolis force applied to the diabatically driven meridional circulation. The present GCM study simulates these processes to a reasonable degree of accuracy. Experimental reduction of the latent energy supplied to the GCM is used to substantiate that summertime subtropical westerly maxima are largely forced by divergent outflow from the tropics.

Hurrell, James W.↗

A time-accurate finite volume method valid at all flow velocities

A finite volume method to solve the Navier-Stokes equations at all flow velocities (e.g., incompressible, subsonic, transonic, supersonic and hypersonic flows) is presented. The numerical method is based on a finite volume method that incorporates a pressure-staggered mesh and an incremental pressure equation for the conservation of mass. Comparison of three generally accepted time-advancing schemes, i.e., Simplified Marker-and-Cell (SMAC), Pressure-Implicit-Splitting of Operators (PISO), and Iterative-Time-Advancing (ITA) scheme, are made by solving a lid-driven polar cavity flow and self-sustained oscillatory flows over circular and square cylinders. Calculated results show that the ITA is the most stable numerically and yields the most accurate results. The SMAC is the most efficient computationally and is as stable as the ITA. It is shown that the PISO is the most weakly convergent and it exhibits an undesirable strong dependence on the time-step size. The degenerated numerical results obtained using the PISO are attributed to its second corrector step that cause the numerical results to deviate further from a divergence free velocity field. The accurate numerical results obtained using the ITA is attributed to its capability to resolve the nonlinearity of the Navier-Stokes equations. The present numerical method that incorporates the ITA is used to solve an unsteady transitional flow over an oscillating airfoil and a chemically reacting flow of hydrogen in a vitiated supersonic airstream. The turbulence fields in these flow cases are described using multiple-time-scale turbulence equations. For the unsteady transitional over an oscillating airfoil, the fluid flow is described using ensemble-averaged Navier-Stokes equations defined on the Lagrangian-Eulerian coordinates. It is shown that the numerical method successfully predicts the large dynamic stall vortex (DSV) and the trailing edge vortex (TEV) that are periodically generated by the oscillating airfoil. The calculated streaklines are in very good comparison with the experimentally obtained smoke picture. The calculated turbulent viscosity contours show that the transition from laminar to turbulent state and the relaminarization occur widely in space as well as in time. The ensemble-averaged velocity profiles are also in good agreement with the measured data and the good comparison indicates that the numerical method as well as the multipletime-scale turbulence equations successfully predict the unsteady transitional turbulence field. The chemical reactions for the hydrogen in the vitiated supersonic airstream are described using 9 chemical species and 48 reaction-steps. Consider that a fast chemistry can not be used to describe the fine details (such as the instability) of chemically reacting flows while a reduced chemical kinetics can not be used confidently due to the uncertainty contained in the reaction mechanisms. However, the use of a detailed finite rate chemistry may make it difficult to obtain a fully converged solution due to the coupling between the large number of flow, turbulence, and chemical equations. The numerical results obtained in the present study are in good agreement with the measured data. The good comparison is attributed to the numerical method that can yield strongly converged results for the reacting flow and to the use of the multiple-time-scale turbulence equations that can accurately describe the mixing of the fuel and the oxidant.

Kim, S.-W.↗

Oscillating cascade aerodynamics at large mean incidence

The aerodynamics of a cascade of airfoils oscillating in torsion about the midchord is investigated experimentally at a large mean incidence angle and, for reference, at a low mean incidence angle. The airfoil section is representative of a modern, low aspect ratio, fan blade tip section. Time-dependent airfoil surface pressure measurements were made for reduced frequencies of up to 1.2 for out-of-phase oscillations at a Mach number of 0.5 and chordal incidence angles of 0 deg and 10 deg; the Reynolds number was 0.9 x l0(exp 6). For the 10 deg chordal incidence angle, a separation bubble formed at the leading edge of the suction surface. The separated flow field was found to have a dramatic effect on the chordwise distribution of the unsteady pressure. In this region, substantial deviations from the attached flow data were found with the deviations becoming less apparent in the aft region of the airfoil for all reduced frequencies. In particular, near the leading edge the separated flow had a strong destabilizing influence while the attached flow had a strong stabilizing influence.

Buffum, Daniel H.↗