Reduction of Large Dynamical Systems by Minimization of Evolution Rate

Reduction of a large system of equations to a lower-dimensional system of similar dynamics is investigated. For dynamical systems with disparate timescales, a criterion for determining redundant dimensions and a general reduction method based on the minimization of evolution rate are proposed

Dynamical System Analysis of Reynolds Stress Closure Equations

In this paper, we establish the causality between the model coefficients in the standard pressure-strain correlation model and the predicted equilibrium states for homogeneous turbulence. We accomplish this by performing a comprehensive fixed point analysis of the modeled Reynolds stress and dissipation rate equations. The results from this analysis will be very useful for developing improved pressure-strain correlation models to yield observed equilibrium behavior

Towards understanding turbulent scalar mixing

In an effort towards understanding turbulent scalar mixing, we study the effect of molecular mixing, first in isolation and then by accounting for the effects of the velocity field. The chief motivation for this approach stems from the strong resemblance of the scalar probability density function (PDF) obtained from the scalar field evolving from the heat conduction equation that arises in a turbulent velocity field. However, the evolution of the scalar dissipation is different for the two cases. We attempt to account for these differences, which are due to the velocity field, using a Lagrangian frame analysis. After establishing the usefulness of this approach, we use the heat-conduction simulations (HCS), in lieu of the more expensive direct numerical simulations (DNS), to study many of the less understood aspects of turbulent mixing. Comparison between the HCS data and available models are made whenever possible. It is established that the beta PDF characterizes the evolution of the scalar PDF during mixing from all types of non-premixed initial conditions

Second Moment Closure Near the Two-component Limit

The purpose of this paper is to explore some wider implications of the two-component limit for both single point turbulence models and spectral closure theories. Although the two-component limit arises most naturally in inhomogeneous problems like wall-bounded turbulence, the analysis will be restricted to homogeneous turbulence. But since homogeneous turbulence is the crucial case for realizability, the conclusions will nevertheless be applicable to modeling. Th essential point of our argument is that whereas the evolution of the stochastic velocity field is Markovian because it is governed by the Navier-Stokes equations, the exact stress evolution equation is not Markovian because it is unclosed. This property of moment evolution has been stressed by Kraichnan (1959). We will show that modeling stress evolution at the two-component limit with a closure that is Markovian in the stresses alone leads to basic inconsistencies in single-point modeling and, perhaps surprisingly, in spectral modes as well

Prandtl number effects on the hydrodynamic stability of compressible boundary layers: flow-thermodynamic interactions

Hydrodynamic stability of compressible boundary layers is strongly influenced by Mach number ($M$), Prandtl number ($Pr$) and thermal wall boundary condition. These effects manifest on the flow stability via the flow-thermodynamic interactions. Comprehensive understanding of stability flow physics is of fundamental interest and important for developing predictive tools and closure models for integrated transition-to-turbulence computations. The flow-thermodynamic interactions are examined using linear analysis and direct numerical simulations (DNS) in the following parameter regime: $0.5 \leq M \leq 8$; and, $0.5 \leq Pr \leq 1.3$. For adiabatic wall boundary condition, increasing Prandtl number has a destabilizing effect. In this work, we characterize the behavior of production, pressure-strain correlation and pressure-dilatation as functions of Mach and Prandtl numbers. First and second instability modes exhibit similar stability trends but the underlying flow physics is shown to be diametrically opposite. The Prandtl-number influence on instability is explicated in terms of the base flow profile with respect to the different perturbation mode shapes

Spectrum and energy transfer in steady Burgers turbulence

The spectrum, energy transfer, and spectral interactions in steady Burgers turbulence are studied using numerically generated data. The velocity field is initially random and the turbulence is maintained steady by forcing the amplitude of a band of low wavenumbers to be invariant in time, while permitting the phase to change as dictated by the equation. The spectrum, as expected, is very different from that of Navier-Stokes turbulence. It is demonstrated that the far range of the spectrum scales as predicted by Burgers. Despite the difference in their spectra, in matters of the spectral energy transfer and triadic interactions Burgers turbulence is similar to Navier-Stokes turbulence

Partially-Averaged Navier Stokes Model for Turbulence: Implementation and Validation

Partially-averaged Navier Stokes (PANS) is a suite of turbulence closure models of various modeled-to-resolved scale ratios ranging from Reynolds-averaged Navier Stokes (RANS) to Navier-Stokes (direct numerical simulations). The objective of PANS, like hybrid models, is to resolve large scale structures at reasonable computational expense. The modeled-to-resolved scale ratio or the level of physical resolution in PANS is quantified by two parameters: the unresolved-to-total ratios of kinetic energy (f(sub k)) and dissipation (f(sub epsilon)). The unresolved-scale stress is modeled with the Boussinesq approximation and modeled transport equations are solved for the unresolved kinetic energy and dissipation. In this paper, we first present a brief discussion of the PANS philosophy followed by a description of the implementation procedure and finally perform preliminary evaluation in benchmark problems

A modified restricted Euler equation for turbulent flows with mean velocity gradients

The restricted Euler equation captures many important features of the behavior of the velocity gradient tensor observed in direct numerical simulations (DNS) of isotropic turbulence. However, in slightly more complex flows the agreement is not good, especially in regions of low dissipation. In this paper, it is demonstrated that the Reynolds-averaged restricted Euler equation violates the balance of mean momentum for virtually all homogeneous turbulent flows with only two major exceptions: isotropic and homogeneously-sheared turbulence. A new model equation which overcomes this shortcoming and is more widely applicable is suggested. This modele is derived from the Navier-Stokes equation with a restricted Euler type approximation made on the fluctuating velocity gradient field. Analytical solutions of the proposed modified restricted Euler equation appear to be difficult to obtain. Hence, a strategy for numerically calculating the velocity gradient tensor is developed. Preliminary calculations tend to indicate that the modified restricted Euler equation captures many important aspects of the behavior of the fluctuating velocity gradients in anisotropic homogeneous turbulence