Search NASA⌕ Search

SEARCH · Search NASA

Results for “TWO-DIMENSIONAL FLOW”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRS1947

Preliminary Tests in the Supersonic Sphere

This report presents preliminary data obtained in the Langley supersonic sphere. The supersonic sphere is essentially a whirling mechanism enclosed in a steel shell which can be filled with either air or Freon gas. The test models for two-dimensional study are of propeller form having the same plan form and diameter but varying only in the airfoil shape and thickness ratio. Torque coefficients for the 16-006, 65-110, and the 15 percent thick ellipse models are presented, as well as pressure distributions on a circular-arc supersonic airfoil section having a maximum thickness of 10 percent chord at the 1/3-chord position. Torque coefficients were measured in both Freon and air on the 15 percent thick ellipse, and the data obtained in air and Freon are found to be in close agreement. The torque coefficients for the three previously mentioned models showed large differences in magnitude at tip Mach numbers above 1, the model with the thickest airfoil section having the largest torque coefficient. Pressure distribution on the previously mentioned circular-arc airfoil section are presented at Mach numbers of 0.69, 1.26, and 1.42. At Mach numbers of 1.26 and 1.42 the test section is in the mixed flow region where both subsonic and supersonic speeds occur on the airfoil. No adequate theory has been developed for this condition of mixed flow, but the experimental data have been compared with values of pressure based on Ackeret's theory. The experimental data obtained at a Mach number of 1.26 on the rear portion of the airfoil section agree fairly well with the values calculated by Ackeret's theory. At a Mach number of 1.42 a larger percentage of the airfoil is in supersonic flow, and the experimental data for the entire airfoil agree fairly well with the values obtained using Ackeret's theory.

Baker, John E.↗
NASA NTRS1947

Effect of Slot-Entry Skirt Extensions on Aerodynamic Characteristics of a Wing Section of the XB-36 Airplane Equipped with a Double Slotted Flap

An investigation was made in the Langley two-dimensional low-turbulence tunnel on a wing section for the XB-36 airplane equipped with a double slotted flap to determine the effect on lift and drag of various slot-entry skirt extension. A skirt extension of 0.787 deg. was found to provide the best combination of high maximum lift with flap deflected and law drag with flap retracted. The data showed that the maximum lift at intermediate (20 deg. to 45 deg.) flap deflections was lowered considerably by the slot-entry extension; but at high flap deflections the effect was small. An increase in Reynolds number from 2.4 million to 6.0 million increased the maximum.lift coefficient at a flap deflection of 55 deg. from 3.12 to 3.30 and from 1.18 to 1.40 for the flap retracted condition, but did not greatly affect the maximum lift coefficient for intermediate flap deflections. The flap and fore flap load data indicated that the maximum lift coefficients at high flap deflections are limited by a breakdown in the flow over the .flaps.

Cahill, Jones F.↗
NASA NTRS1945

Summary of Airfoil Data

Recent airfoil data for both flight and wind-tunnel tests have been collected and correlated insofar as possible. The flight data consist largely of drag measurements made by the wake-survey method. Most of the data on airfoil section characteristics were obtained in the Langley two-dimensional low-turbulence pressure tunnel. Detail data necessary for the application of NACA 6-serles airfoils to wing design are presented in supplementary figures, together with recent data for the NACA 24-, 44-, and 230-series airfoils. The general methods used to derive the basic thickness forms for NACA 6- and 7-series airfoils and their corresponding pressure distributions are presented. Data and methods are given for rapidly obtaining the approximate pressure distributions for NACA four-digit, five-digit, 6-, and 7-series airfoils. The report includes an analysis of the lift, drag, pitching-moment, and critical-speed characteristics of the airfoils, together with a discussion of the effects of surface conditions. Available data on high-lift devices are presented. Problems associated with lateral-control devices, leading-edge air intakes, and interference are briefly discussed. The data indicate that the effects of surface condition on the lift and drag characteristics are at least as large as the effects of the airfoil shape and must be considered in airfoil selection and the prediction of wing characteristics. Airfoils permitting extensive laminar flow, such as the NACA 6-series airfoils, have much lower drag coefficients at high speed and cruising lift coefficients than earlier types-of airfoils if, and only if, the wing surfaces are sufficiently smooth and fair. The NACA 6-series airfoils also have favorable critical-speed characteristics and do not appear to present unusual problems associated with the application of high-lift and lateral-control devices. Much of the data given in the NACA Advance Confidential Report entitled "Preliminary Low-Drag-Airfoil and Flap Data from Tests at Large Reynolds Number and Low Turbulence," by Eastman N. Jacobs, Ira R. Abbott, and Milton Davidson, March 1942 has been corrected and included in the present paper, which supersedes the previously published paper.

LOW-DRAG↗
NASA NTRS1944

An Experimental Investigation of the Flow of Air in a Flat Broadening Channel

The wide use of diffusers, in various fields of technology, has resulted in several experimental projects to study the action and design of diffusers. Most of the projects dealt with steam (steam turbine nozzles). But diffusers have other applications - that is, ventilators, smoke ducts, air coolers, refrigeration, drying, and so forth. At present there is another application for diffusers in wind-tunnel design. Because of higher requirements and increased power of such installations more attention must be paid to the correctness of work and the decrease in losses due to every section of the tunnel. A diffuser, being one of the component parts of a tunnel , can in the event of faulty construction introduce considerable losses. Therefore, in the design of the new CAHI wind tunnel, it was suggested that an experimental study of diffusers be made, with a view to applying the results to wind tunnels. The experiments conducted by K. K. Baulin in the laboratories of CAHI upon models of diffusers of different cross sections, lengths, and angles of divergence, were a valuable source of experimental data. They were of no help, however, in reaching any conclusion regarding the optimum shape because of the complexity and diversity of the factors which all appeared simultaneously, thereby precluding the.study of the effects of any one factor separately. On the suggestion of the director of the CAHI,Prof. B. N. Ureff, it was decided to experiment on a two-dimensional diffuser model and determine the effect, of the angle of divergence. The author is acquainted with two experimental projects of like nature: the first was conducted with water, the other with air. The first of these works, although containing a wealth of experimental data, does not indicate the nature of flow or its relation to the angle of divergence. The second work is limited to four angles - that is, 12 deg, 24 deg, 45 deg, 90 deg. The study of this diffuser did not supply any information about the effect of smaller angles which, because of their advantages, are more commonly used, The author was able to acquaint himself with the second work only after the experiments were started. For these reasons, as well as because on the basis of those works no conclusion can be reached regarding the nature of flow distribution, of eddies, and so forth, experimental work was continued. The need for determining flow patterns follows from the fact that from them are determined methods of measurement - that is, the determination of velocities by means of the pitot tube, which, as is well known, gives correct indications only when placed with its axis parallel to the axis of flow. The data contained in this report were obtained from experiments conducted by the Aerodynamical Laboratories of the CAHI. The solutions to some. of the mathematical problems connected with the experiments are due to Prof. S. A, Chapligin.

Vedernikoff, A. N.↗
DOE OSTI2026

Simulations of a hypersonic turbulent boundary layer over wavy surfaces

Here, we conduct large-eddy simulations of a Mach 5.84 cold wall turbulent boundary layer over one-dimensional wavy walls with varying amplitudes and wavelengths. Across all wall topologies, a series of alternating shock and expansion waves is shown to influence the entire boundary layer, and generate repeating wave patterns in the turbulent stresses, dispersive stresses, and turbulent kinetic energy budget. The series of alternating shocks and expansions imposes repeating adverse and favourable pressure gradients across the wavy wall, and at sufficient wall amplitude, triggers flow separation in the trough of the wave. Flow separation is demonstrated to influence the behaviour of wall pressure fluctuations over the wavy wall. In attached flows, the prominent frequencies are consistent with integral-scale boundary layer turbulence, whereas in separated flows, a two-decade frequency range is present, akin to two-dimensional shock–boundary layer interactions. Counter-rotating streamwise-oriented structures are observed on the windward side of the wave, which diminish over the wave crest. A conditional analysis demonstrates that these structures are present in the upstream boundary layer, and are amplified with increasing wall amplitude. An examination of the Görtler number and probability density function (PDF) of the fluctuating lateral wall shear stress demonstrates the strong correlation between a large Görtler number and growth of the PDF tail density, suggesting that the amplification of the counter-rotating streamwise-oriented structures are linked to centrifugal instabilities in regions of streamline concavity.

boundary layers↗
DOE OSTI2026

Studies of localized mix and asymmetry caused by the stalk in inertial fusion experiments at OMEGA

In an inertially confined fusion experiment, a capsule of fusion fuel is compressed to high temperatures and densities using lasers, creating a hot plasma in which thermonuclear fusion reactions take place. To achieve efficient conversion of laser energy to fuel heating, the implosion must be as spherically symmetric as possible. One common source of asymmetry that degrades performance is the stalk and glue spot that hold the capsule in place before the experiment. The exact nature and magnitude of this perturbation had not been previously quantified. In this work, inertially confined fusion experiments were executed with an additional stalk on the capsule to measure and quantify the effect of the stalk and glue spot on the symmetry and performance. It is found that the stalk drives low-mode asymmetry by shadowing the laser drive on the part of the shell near the glue spot, enhancing the hot-spot velocity toward the stalk by 40–80 km/s and elongating the hot-spot along the stalk axis. It is also found that localized mix due to the additional stalk has a negligible effect on implosion performance. Two-dimensional radiation-hydrodynamic simulations using the xRAGE code reproduce the direction of the experimentally observed change to hot-spot velocity and implosion performance metrics, but not the observed elongation, suggesting that this feature results from an interaction between stalk perturbations and other preexisting asymmetry seeds, such as laser drive asymmetry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Mode transitions and spoke structures in E × B Penning discharge

Two-dimensional particle-in-cell simulations in the (radial-azimuthal) plane perpendicular to the axial direction of a cylindrical $E$ x $B$ Penning discharge are presented. The low-pressure discharge is self-consistently supported by plasma ionization from the electron beam injected axially, along the direction of the external magnetic field. It is shown that with the increasing strength of the external magnetic field, the discharge undergoes a sequence of transitions between several azimuthal modes. Azimuthal m > 1 spiral arm structures are excited at low magnetic field values as plasma confinement improves and the radial density profile becomes peaked. With a larger field, spiral arms with m > 1 are replaced by the m = 1 spoke mode, most clearly seen in plasma density. A transition from spiral arms to the spoke regime occurs when the plasma potential in the center changes from weakly positive (or zero) to negative. Further increase in the magnetic field results in a well-developed m = 1 spoke mode with additional small-scale higher-frequency m > 1 structures inside and around the spoke. It is shown that while ionization and collisions affect some characteristics of the observed fluctuations, the basic features of the spoke and m > 1 spiral structure remained similar without ionization. The role of energy conservation in small-scale high-frequency modes and spoke dynamics is discussed. It is demonstrated that in regimes with the m = 1 spoke mode, additional m = 4 harmonics of the ion and electron fluxes to the wall appear due to the square boundary. The frequency of the m = 1 mode is weakly affected by the geometry of the boundary.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

A discontinuous Galerkin spectral element method for compressible reacting flows

High-order methods have recently been shown to be an effective tool for high-fidelity flow computations like direct numerical simulations and large-eddy simulations because of their strong balance between accuracy and computational cost. In this work, a high-order discontinuous Galerkin spectral element method (DGSEM) is developed to solve the chemically reacting Navier-Stokes equations. To handle the disparate length and time scales associated with these equations, we develop a novel method which combines the spectral accuracy of the SEM with the flexibility of the DG approach. The framework, implemented in the spectral element code Nek5000, is well suited to capture turbulence in smooth regions of the flow, while maintaining numerical stability in the presence of shocks. An entropy-residual based artificial viscosity is added to smooth shocked regions of flow, and a positivity-preserving limiter is implemented to suppress non-physical oscillations. These enhancements support the numerical stability of the hydrodynamic sub-step, which is decoupled from the chemistry integration through a second-order operator splitting method. Here, a series of smooth and discontinuous validation cases are presented in increasing physical and computational complexity for both inviscid and viscous flows. In particular, simulations of canonical one-dimensional and two-dimensional detonations are performed, and the high-order numerical results are validated against available literature data. Additional validation studies are carried out for classical three-dimensional numerical simulations of incompressible and compressible turbulent flows.

Compressible reacting flows↗
DOE OSTI2026

Time-dependent-bases with local CUR decomposition method for accelerating turbulent combustion simulations

Here, this study presents a novel reduced-order modeling framework, Time-Dependent Bases with Local CUR decomposition (TDB-L-CUR), designed to efficiently and accurately approximate the species transport equations in reacting flow simulations. The method extends the existing TDB-CUR approach for chemically reacting flows (Jung et al. Comput. Methods Appl. Mech. Engrg. 437 (2025) 117758), which leverages matrix decomposition techniques to form a global-in-space, time-dependent low-dimensional manifold. While TDB-CUR performs well in homogeneous systems, it may be less well-suited to spatially heterogeneous systems such as turbulent flames, where higher-rank approximations are typically required. The proposed TDB-L-CUR framework introduces two methodological extensions to the baseline approach. First, it applies unsupervised clustering to partition the physical domain into distinct regions, enabling spatially localized manifold construction, thereby reducing the rank required for the reduced-order representation. Second, it incorporates a computational singular perturbation (CSP)-based scheme for identifying and penalizing fast species, allowing for spatio-temporally adaptive mitigation of chemical stiffness. The proposed framework is validated on a hierarchy of test cases, including a one-dimensional premixed flame, a two-dimensional nonpremixed ignition case with vortex interaction, and a three-dimensional turbulent premixed flame. TDB-L-CUR significantly improves accuracy over TDB-CUR while further reducing computational cost. The fully on-the-fly formulation of TDB-L-CUR (i.e., requiring no offline training or prior knowledge) makes it a robust and scalable tool for reduced-order modeling of reactive flows.

Local manifold↗
DOE OSTI2025

Generative AI models for learning flow maps of stochastic dynamical systems in bounded domains

Simulating stochastic differential equations (SDEs) in bounded domains, presents significant computational challenges due to particle exit phenomena, which requires accurate modeling of interior stochastic dynamics and boundary interactions. Despite the success of machine learning-based methods in learning SDEs, existing learning methods are not applicable to SDEs in bounded domains because they cannot accurately capture the particle exit dynamics. We present a unified hybrid data-driven approach that combines a conditional diffusion model with an exit prediction neural network to capture both interior stochastic dynamics and boundary exit phenomena. Our ML model consists of two major components: a neural network that learns exit probabilities using binary cross-entropy loss with rigorous convergence guarantees, and a training-free diffusion model that generates state transitions for non-exiting particles using closed-form score functions. The two components are integrated through a probabilistic sampling algorithm that determines particle exit at each time step and generates appropriate state transitions. Here, the performance of the proposed approach is demonstrated via three test cases: a one-dimensional simplified problem for theoretical verification, a two-dimensional advection-diffusion problem in a bounded domain, and a three-dimensional problem of interest to magnetically confined fusion plasmas.

Bounded domains↗
DOE OSTI2025

Simulation of an inductively coupled plasma with a two-dimensional Darwin particle-in-cell code

A two-dimensional particle-in-cell code for the simulation of low-frequency electromagnetic processes in laboratory plasmas has been developed. The code uses the Darwin method omitting the electromagnetic wave propagation. The Darwin method separates the electric field into solenoidal and irrotational parts. The solenoidal electric field is calculated with a new algorithm based on the equation for the electric field vorticity. The system of linear equations in the new algorithm is readily solved using a standard iterative method. The irrotational electric field is the electrostatic field calculated with the direct implicit algorithm. The code is verified by reproducing the two-stream instability, electron electromagnetic waves, and shear Alfvén waves. The code is applied to simulate an inductively coupled plasma with the driving current flowing around the plasma region. In this simulation, a ring of dense plasma forms at the initial stage but then the density becomes maximal in the center and decays monotonically toward the walls. The skin effect is in the transitional mode between local and non-local, and the electron velocity distribution function is non-Maxwellian.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

A two-dimensional numerical study of the magneto-Rayleigh–Taylor instability with FLASH: Application to the staged Z-pinch concept

Magnetically driven implosions involving a liner collapsing onto a target are inherently vulnerable to the magneto-Rayleigh–Taylor instability (MRTI). Among the various approaches proposed to achieve fusion conditions within the target, the staged Z-pinch (SZP) concept employs a high-Z liner, the advantages of which remain an active area of investigation. Consequently, ongoing design optimization efforts are essential, while critical physical processes such as magnetic field (B-field) diffusion and radiation transport increase the complexity of required simulations. In this study, we utilize the new capabilities of the FLASH code to simulate in 2D, for the first time using this code, staged Z-pinch configurations (designated SZP1 and SZP1*), focusing on their stability. First, a comparison of simulation results with theoretical predictions of MRTI growth provide new insights into the often-overlooked influences of high mesh resolution and initial perturbation seeding on instability dynamics. These findings then lead to a discussion on potential improvements for future SZP experiments. These include the use of an axial B-field stabilization and the optimization of radiation transport processes during the implosion. The results presented herein establish the framework for simulating multi-dimensional Z pinches using the FLASH code and pave the way for the development of innovative experimental configurations leveraging its advanced simulation capabilities.

Adaptive mesh refinement↗
DOE OSTI2025

Multimode turbulent flow measurements using magnetic resonance imaging- and laser-based techniques and computational fluid dynamics simulations

We studied the flow field characteristics of a turbulent flow over a regularized cube array with a perpendicular injection flow through the floor between the second and third cubical elements, representing the complex flow interactions of a 3D jet and the wake flows behind cubical obstacles. Four different experimental measurements were performed: two magnetic resonance imaging-based measurements for three-dimensional three-component velocity (MRV) and concentration (MRC) and two laser-based techniques, particle image velocimetry (PIV) and planar laser-induced fluorescence (PLIF), for two-dimensional two-component velocity and concentration measurement, respectively. The mainstream Reynolds number is Re = 15 000⁠, based on the primary inlet velocity U m and channel height D H ⁠, whereas the injector Reynolds number is Re j = 3400⁠, based on the injector velocity U j and the injector's exit width D j ⁠. Numerical simulations were performed for the studied flow configuration of turbulent flow over a regularized cube array using Reynolds-averaged Navier–Stokes (RANS) and large-eddy simulation (LES) approaches. Results obtained from experimental measurements—including MRV, MRC, PIV, and PLIF—as well as RANS and LES simulations are discussed and compared along several horizontal and vertical planes of the studied configuration. In addition, 3D turbulent flow structures, such as leading-edge vortex, horseshoe vortex, and jet shear-layer vortex, and the isosurfaces of scalar concentration successfully revealed by MRV and MRC techniques were found to be in very good agreement with those 3D features extracted from RANS and LES simulations. In conclusion, the high-resolution experimental and numerical database obtained from this study could be useful for validation and verification of numerical codes.

Computational fluid dynamics↗
DOE OSTI2026

Data-based filtered dissipation rate modelling for multi-modal turbulent combustion: evaluating a priori model generalizability

Manifold-based models offer a computationally efficient alternative to directly transporting the thermochemical state in computational simulations of turbulent reacting flows, projecting the high-dimensional thermochemical state-space onto a low-dimensional manifold. Recent efforts have yielded a manifold-based model applicable to multi-modal combustion, enabling reconstruction of the thermochemical state from solutions to two-dimensional manifold equations in mixture fraction and generalized progress variable that are parameterised by three scalar dissipation rates. In coarse-grained simulations such as Large Eddy Simulation (LES), closure of the multi-modal manifold equations and subfilter variances/covariance requires closure of three filtered scalar dissipation rates. Here, the present work adopts a data-based approach, providing closure for the three filtered scalar dissipation rates via deep neural networks (DNNs). High-fidelity datasets corresponding to an autoigniting n-dodecane jet flame and a bluff body swirl-stabilized confined lifted spray flame of two aviation fuels (Jet-A and C1) with different ignition propensities are leveraged to generate training data that spans a diverse range of thermodynamic conditions and combustion modes, including low- and high-temperature ignition regimes in addition to premixed and nonpremixed behaviour. A final DNN model is trained to enforce inherent physical constraints by learning nonlinear functional transformations of the three filtered scalar dissipation rates. The generalizability of this constrained DNN model is demonstrated a priori via conditional statistics evaluated on the lifted spray flame with C1–a configuration that had not been included in the training data. Excellent DNN agreement with conditional DNS statistics is observed, and integrated gradients are computed to identify the most sensitive input variables. The similarity of the marginal PDFs of the most informative input variables and outputs across configurations are quantified via the Wasserstein metric, demonstrating that data-based models may successfully generalize to unseen parametric conditions so long as the most informative input variables share similar distributions across training and testing datasets.

Data-based modelling↗
DOE OSTI2025

Depinning, melting, and sliding of generalized Wigner crystals in Moiré systems

We numerically examine the depinning, sliding, and melting of commensurate and incommensurate Wigner crystals on two-dimensional hexagonal periodic substrates near fillings of 1 / 3 , 1 / 2 , and 2 / 3 to model the dynamics of generalized Wigner crystals in moiré heterostructures. At low temperatures where thermal fluctuations are irrelevant, for commensurate fillings we find a strongly pinned state that depins into a sliding crystal, while at incommensurate fillings the depinning threshold is strongly reduced. For fillings below 1 / 2 , the depinning occurs in a two-step process. Above the first depinning threshold, there is an extended range of drives where the conduction occurs via the sliding of antikinks along the charge stripes, followed by a second threshold where all the charges begin to slide. At finite temperatures, the external driving reduces the effective melting temperature at commensurate fillings and enhances creep at incommensurate fillings. We show that depinning into different sliding states, such as moving fluids or moving crystals, produces nonlinear features in the transport curves. We also show that transport is asymmetric on either side of a commensurate filling due to the different dynamics for interstitials versus holes in the commensurate structure. A variety of sliding states should be accessible for generalized Wigner crystals at finite temperatures even for low drives. If experiments can realize sliding dynamics in moiré systems, it would open another class of commensurate and incommensurate phases for study. Additionally, the sliding dynamics and antikink flow could provide unique functionalities for charge transport in moiré heterostructures. Published by the American Physical Society 2025

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Measurement of Local Impedance Characteristics of In-Plane PEM Electrolyzer Component Features via Segmented Cell

Segmented cells are a class of advanced diagnostic devices that enable the measurement of current distribution within the active area of an electrochemical cell, allowing for the evaluation of localized impacts from operating conditions, flow fields and component inhomogeneities. Segmented cell devices have been applied to support R&D efforts for the advancement of polymer electrolyte membrane fuel cell technologies and more recently proton exchange membrane water electrolysis (PEMWE). Beyond current density distributions, segmented cells have been successfully deployed to collect the distribution of other relevant parameters such as high frequency resistance and electrode potentials. Very recently, a promising impedance method has been reported for water electrolysis where a one-dimensional segmented cell is coupled to a multichannel potentiostat, to enable the measurement of local cell impedances, further expanding the diagnostic capabilities of the segmented cell. In this study, we will present recent efforts to accurately measure local impedance characteristics for a PEMWE cell using a two-dimensional segmented cell device coupled with a multichannel potentiostat. The distribution of local resistances, i.e. high frequency, charge transfer, mass transport and catalyst layer, and iR-free overpotentials were successfully measured after eliminating sources of random and systematic error. We will discuss the application of this diagnostic on a cell containing a 1 cm feature simulating an ionomer skin. The results from these tests give insights if such an irregularity needs to be classified as a defect, and highlight the usefulness of coupling impedance spectroscopy with segmented cells for holistic spatial diagnostics of electrochemical devices.

08 HYDROGEN↗
DOE OSTI2025

A High-Order Discontinuous Galerkin Spectral Element Method for Compressible Reacting Flows

High-order methods have recently been shown to be an effective tool for high-fidelity flow computations like direct numerical simulations and large eddy simulations due to their strong balance between accuracy and computational cost. In this work, a high-order discontinuous Galerkin spectral element method (DGSEM) is developed to solve the chemically reactive Euler equations encountered in high-speed combustion. To handle the disparate length and time scales associated with these equations, we develop a novel method which combines the spectral accuracy of the SEM with the flexibility of DG approach. Thus, the framework is well suited to capture turbulence in smooth regions of the flow, while maintaining numerical stability in the presence of shocks. The numerical method is implemented within the spectral element solver Nek5000. Validation cases are conducted for both non-reactive and reactive discontinuous flows to demonstrate the solver capability. In particular, canonical one-dimensional and two-dimensional detonation simulations are performed and the high-order numerical results are validated against available literature data.

computational fluid dynamics (CFD)↗
NASA NTRS1947

Investigation of Axial-flow Fan and Compressor Rotors Designed for Three-dimensional Flow

An investigation has been conducted to determine whether three-dimensional flows may be utilized in axial-flow fan and compressor rotors so that the spanwise load distribution may be varied to obtain high pressure rise. Two rotors, one with approximately uniform and one with solid-body downstream tangential-velocity distributions, were designed and tested at the design blade angle. Radial surveys of total pressure, static pressure, and flow angle were made upstream and downstream of the test rotors through a quantity-coefficient range. Tests of the solid-body rotor were also conducted at a large value of tip clearance. The results indicated that the three-dimensional flows may be utilized with high efficiency and that the three-dimensional theory used in conjunction with two-dimensional cascade data is sufficiently accurate for design purposes. The tests also showed that the tip-clearance losses of rotors highly loaded at the tips are not excessive. The existing three-dimensional theory in simplified for and an illustrative rotor design are presented in appendixes.

Kahane, A↗