Accuracy Consideration by DRP Schemes for DNS & LES of Compressible Flow Computations
Explore the source record for details and available documents.
SEARCH · Search NASA
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.
Explore the source record for details and available documents.
As next-generation exascale-class systems arrive, existing software must be updated accordingly to effectively utilize these systems. For high concurrency and energy efficiency, many of these systems utilize GPU architectures. In this work, we present a CUDA C++ implementation of FUN3D's thermochemical nonequilibrium capability for turbulent flows. Efficiency is demonstrated at scale using the Summit system at the Oak Ridge Leadership Computing Facility which is representative of future exascale systems. This work enables faster, higher fidelity, and scale-resolving simulations of thermochemical nonequilibrium flows including reentry, hypersonics, and combustion.
Explore the source record for details and available documents.
In the present work, we use a burst-mode laser and optical parametric oscillator system to perform high-spectral resolution NO-PLIF measurements of an underexpanded jet at a repetition-rate of 100 kHz, with the motivation of multi-parameter measurements of temperature, pressure, and velocity. The laser frequency of the 1064 nm seed laser for the burst-mode laser was scanned during the burst to cover two neighboring absorption line pairs near 226 nm. The peak PLIF signal intensity varies along the axial (z) direction of the underexpanded jet as the laser frequency is scanned, which we attribute to the collisional shift induced by the flow and only revealed due to the narrow linewidth of the laser. A pseudo-Voight fit is applied to the LIF excitation spectra on a pixel-by-pixel basis to measure the spectral position of the peak intensity for the two transition pairs and their amplitude. The spectral position of the peak intensity is used to derive a frequency shift, which is separated into its collisional and Doppler components using the axisymmetric nature of the flow field. The amplitude of the measured peaks is used for two-line rotational thermometry. Challenges for quantitative measurements using such an approach are discussed, including measuring the spatial variations in the energy distribution of the laser sheet at a 100 kHz repetition rate and uncertainty/variability in the step size during the fast frequency scan.
In the present work, we use a burst-mode laser and optical parametric oscillator system to perform high-spectral resolution NO-PLIF measurements of an underexpanded jet at a repetition-rate of 100 kHz, with the motivation of multi-parameter measurements of temperature, pressure, and velocity. The laser frequency of the 1064 nm seed laser for the burst-mode laser was scanned during the burst to cover two neighboring absorption line pairs near 226 nm. The peak PLIF signal intensity varies along the axial (z) direction of the underexpanded jet as the laser frequency is scanned, which we attribute to the collisional shift induced by the flow and only revealed due to the narrow linewidth of the laser. A pseudo-Voight fit is applied to the LIF excitation spectra on a pixel-by-pixel basis to measure the spectral position of the peak intensity for the two transition pairs and their amplitude. The spectral position of the peak intensity is used to derive a frequency shift, which is separated into its collisional and Doppler components using the axisymmetric nature of the flow field. The amplitude of the measured peaks is used for two-line rotational thermometry. Challenges for quantitative measurements using such an approach are discussed, including measuring the spatial variations in the energy distribution of the laser sheet at a 100 kHz repetition rate and uncertainty/variability in the step size during the fast frequency scan.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
A study reported by Yang and Hudson (1971) was extended by Nelson (1971), who presented a method of inverse design for axisymmetric diffusers. The diffusers obtained with the aid of this method were found to function well under certain conditions, while, for other cases, a modification of the analytical design procedure was required. An outline is provided of the method of inverse solution in simple axisymmetric internal flow passages. Attention is given to governing equations, the method of solution, and examples which illustrate the feasibility of the design procedure. The discussed method, which is based on viscous compressible flow theory, has some limitations. However, it is expected to yield good designs in many practical cases in which the existing design must be modified.
The steady full-potential equation is written in the form of Poisson's equation, and the solution for the velocity field is expressed in terms of an integral equation. The integral solution consists of two surface integrals and one volume integral. The solution is obtained through successive iteration cycles. Each cycle of iteration consists of two sub-cycles, an inner cycle for wake relaxation and an out cycle for the strength of the source distribution integrals representing the flow compressibility. The density gradients in the source distribution is computed by using a type-differencing scheme of the Murman-Cole type. The method is applied to delta wings and the numerical examples show that a curved shock is captured on the wing suction side beneath the leading edge vortex sheet. Recently, a modified version of the scheme was applied to rectangular wings. In this modified scheme, the surface integral terms were computed by using a bilinear distribution of vorticity on triangular vortex panels which represent the wing and its wake. The results were compared with the available experimental data and they are in good agreement.
Compressible stability equations are solved using the spectral collocation method in an attempt to study the effects of temperature difference and compressibility on the stability of Taylor-Couette flow. It is found that the Chebyshev collocation spectral method yields highly accurate results using fewer grid points for solving stability problems. Comparisons are made between the result obtained by assuming small Mach number with a uniform temperature distribution and that based on fully incompressible analysis.
Galerkin reduced-order models (ROMs) often struggle to accurately capture multiscale fluid physics in challenging flow regimes such as flows experiencing compressibility effects. This in part stems from the global nature of both the basis construction problem and the spectral formulation itself. In this work, a multi-basis ROM is developed in an affine space based on proper orthogonal decomposition (POD) by projecting the full Navier-Stokes equations expressed in terms of the specific volume, velocity, and pressure primitive variables. The model is applied to high-fidelity numerical simulation datasets obtained for a canonical compressible flow configuration: the flow over a backward facing step at different subsonic Mach numbers. It is observed that application of an eigenvalue reassignment (ER) stabilization method is required to avoid early divergence of the ROM predictions for this configuration in the three Mach numbers tested. The sensitivity of the POD-ROM results to the choice of parameters in the stabilization algorithm is discussed.
Galerkin reduced-order models (ROMs) often struggle to accurately capture multiscale fluid physics in challenging flow regimes such as flows experiencing compressibility effects. This in part stems from the global nature of both the basis construction problem and the spectral formulation itself. In this work, a multi-basis ROM is developed in an affine space based on proper orthogonal decomposition (POD) by projecting the full Navier-Stokes equations expressed in terms of the specific volume, velocity, and pressure primitive variables. The model is applied to high-fidelity numerical simulation datasets obtained for a canonical compressible flow configuration: the flow over a backward facing step at different subsonic Mach numbers. It is observed that application of an eigenvalue reassignment (ER) stabilization method is required to avoid early divergence of the ROM predictions for this configuration in the three Mach numbers tested. The sensitivity of the POD-ROM results to the choice of parameters in the stabilization algorithm is discussed.
The purpose is to re-examine the heat transfer from a hot-wire probe in the compressible subsonic flow regime; describe the three-wire hot-wire probe calibration and data reduction techniques used to measure the velocity, density, and total temperature fluctuation; and present flow quality results obtained in the Langley 0.3 meter Transonic Cryogenic Wind Tunnel and in flight with the NASA JetStar from the same three-wire hot-wire probe.
A compressible Reynolds averaged Navier-Stokes code is applied to examine the performances of basic turbulence models for unseparated and separated flows. The turbulence models considered are a zero-equation model and several two-equation models, including a new k-epsilon model with an eddy viscosity damping function depending upon the Reynolds number and the distance from the wall. The turbulence models are assessed using the following experimental flows: a channel flow, a backward facing step, and a two-dimensional compression corner at Mach 2.8. Discussions and comparisons of the experimental and numerical results are given. Mesh refinement is shown to have a significant beneficial effect on free shear layer flow results obtained with two-equation turbulence models. The results illustrate the usefulness of several turbulence models as a design tool for fluid engineering systems.
A conference was held on Computational Fluid Dynamics (CFD) and produced related papers. Topics included CFD algorithms, transition and turbulent flow, hypersonic reacting flow, incompressible flow, two phase flow and combustion, internal flow, compressible flow, grid generation and adaption, boundary layers, environmental and industrial applications, and non-Newtonian flow.
Compressibility effects on the turbulence in homogeneous shear flow are investigated. The growth of the turbulent kinetic energy was found to decrease with increasing Mach number: a phenomenon which is similar to the reduction of turbulent velocity intensities observed in experiments on supersonic free shear layers. An examination of the turbulent energy budget shows that both the compressible dissipation and the pressure-dilatation contribute to the decrease in the growth of kinetic energy. The pressure-dilatation is predominantly negative in homogeneous shear flow, in contrast to its predominantly positive behavior in isotropic turbulence. The different signs of the pressure-dilatation are explained by theoretical consideration of the equations for the pressure variance and density variance. Previously, the following results were obtained for isotropic turbulence: (1) the normalized compressible dissipation is of O(M(sub t)(exp 2)); and (2) there is approximate equipartition between the kinetic and potential energies associated with the fluctuating compressible mode. Both of these results were substantiated in the case of homogeneous shear. The dilatation field is significantly more skewed and intermittent than the vorticity field. Strong compressions seem to be more likely than strong expansions.
Binary vortex flow fields used as gaseous-fuel- retention devices for gas-core nuclear reactors