Results for “closure theory”

NASA NTRS1995

Classical closure theory and Lam's interpretation of epsilon-RNG

Lam's phenomenological epsilon-renormalization group (RNG) model is quite different from the other members of that group. It does not make use of the correspondence principle and the epsilon-expansion procedure. We demonstrate that Lam's epsilon-RNG model is essentially the physical space version of the classical closure theory in spectral space and consider the corresponding treatment of the eddy viscosity and energy backscatter.

Zhou, YE↗

NASA NTRS2011

Non-Equilibrium Turbulence and Two-Equation Modeling

Two-equation turbulence models are analyzed from the perspective of spectral closure theories. Kolmogorov theory provides useful information for models, but it is limited to equilibrium conditions in which the energy spectrum has relaxed to a steady state consistent with the forcing at large scales; it does not describe transient evolution between such states. Transient evolution is necessarily through nonequilibrium states, which can only be found from a theory of turbulence evolution, such as one provided by a spectral closure. When the departure from equilibrium is small, perturbation theory can be used to approximate the evolution by a two-equation model. The perturbation theory also gives explicit conditions under which this model can be valid, and when it will fail. Implications of the non-equilibrium corrections for the classic Tennekes-Lumley balance in the dissipation rate equation are drawn: it is possible to establish both the cancellation of the leading order Re1/2 divergent contributions to vortex stretching and enstrophy destruction, and the existence of a nonzero difference which is finite in the limit of infinite Reynolds number.

Rubinstein, Robert↗

NASA NTRSDate not supplied

Demonstration of Capability to Simulate Particle Irregular Shape and Poly-Disperse Mixtures Within Lunar Lander Plume-Surface Interaction

Plume-Surface Interaction (PSI) between lander engine plumes and regolith soil creates hazards in obscuration and contamination by particle clouds, high-energy ejecta streams, and landing area cratering damage. The MSFC Fluid Dynamics Branch is developing simulation tools to offer a predictive PSI capability to NASA customers such as the Human Lander System (HLS) and Commercial Lunar Payload Services (CLPS). The Gas-Granular Flow Solver (GGFS) is the main application tool for coupled gas-particle two-phase flow simulations to predict the range of PSI effects from onset of surface erosion to deep crater formation. GGFS features an Eulerian-Eulerian modeling approach, treating both gas and granular material as interacting continuum phases. Modeling the lunar regolith granular material fluidic characteristics poses special challenges due to complex particle shapes and mixture composition. The lunar regolith is poorly sorted with broad particle size distributions and large fines content. It has significant cohesion, due to interlocking jagged particle shapes. Eulerian granular material flow modeling requires closure formulations for the granular material constitutive models (stress, friction, collisional and kinetic energy dissipation, drag, etc.). While closure models for spherical particles are available from particle kinetic theory, closure models for realistic non-spherical particles must be extracted from unit physics Discrete Element Model (DEM) particle interaction simulations and provided in the form of tabular datasets. The effects of particle irregular shape (non-spherical shape factors, angular particle surface roughness, and interlocking features) are simulated by approximating the particle features in the form of grouped elemental spheres to form composite particles in the DEM simulations. The effects of the wide range of regolith mixture particle sizes and the strong effects of the presence of the small particle sizes results in high cohesion and low porosity of the regolith mixture. The range of particle sizes is simulated by binning the particle sizes into an appropriate finite number of particle-size species and solving the problem as a species mixture. Combining these two modeling approaches enables simulations to capture both, the contributions of the irregular particle shape and the particle size distribution. The integration and maturation of the DEM-based constitutive model database generation process and poly-disperse mixture binning approach into the GGFS simulation framework are proceeding under funding by the NASA Game Changing Development program. The status of current capabilities will be presented in comparisons of crater characteristics resulting for spherical and irregular shape particles, and for mono-, bi-, and tri-disperse mixture simulations of Apollo LM plume-surface interaction. The computational results confirm the significance of including the particle shape and mixture effects. Going forward plans for the full implementation of the general poly-disperse regolith modeling capability and maturation towards NASA project application readiness under the GCD program will be presented.

Peter A Liever↗

NASA NTRSDate not supplied

Demonstration of Capability to Simulate Particle Irregular Shape and Poly-Disperse Mixtures Within Lunar Lander Plume-Surface Interaction

Plume-Surface Interaction (PSI) between lander engine plumes and regolith soil creates hazards in obscuration and contamination by particle clouds, high-energy ejecta streams, and landing area cratering damage. The MSFC Fluid Dynamics Branch is developing simulation tools to offer a predictive PSI capability to NASA customers such as the Human Lander System (HLS) and Commercial Lunar Payload Services (CLPS). The Gas-Granular Flow Solver (GGFS) is the main application tool for coupled gas-particle two-phase flow simulations to predict the range of PSI effects from onset of surface erosion to deep crater formation. GGFS features an Eulerian-Eulerian modeling approach, treating both gas and granular material as interacting continuum phases. Modeling the lunar regolith granular material fluidic characteristics poses special challenges due to complex particle shapes and mixture composition. The lunar regolith is poorly sorted with broad particle size distributions and large fines content. It has significant cohesion, due to interlocking jagged particle shapes. Eulerian granular material flow modeling requires closure formulations for the granular material constitutive models (stress, friction, collisional and kinetic energy dissipation, drag, etc.). While closure models for spherical particles are available from particle kinetic theory, closure models for realistic non-spherical particles must be extracted from unit physics Discrete Element Model (DEM) particle interaction simulations and provided in the form of tabular datasets. The effects of particle irregular shape (non-spherical shape factors, angular particle surface roughness, and interlocking features) are simulated by approximating the particle features in the form of grouped elemental spheres to form composite particles in the DEM simulations. The effects of the wide range of regolith mixture particle sizes and the strong effects of the presence of the small particle sizes results in high cohesion and low porosity of the regolith mixture. The range of particle sizes is simulated by binning the particle sizes into an appropriate finite number of particle-size species and solving the problem as a species mixture. Combining these two modeling approaches enables simulations to capture both, the contributions of the irregular particle shape and the particle size distribution. The integration and maturation of the DEM-based constitutive model database generation process and poly-disperse mixture binning approach into the GGFS simulation framework are proceeding under funding by the NASA Game Changing Development program. The status of current capabilities will be presented in comparisons of crater characteristics resulting for spherical and irregular shape particles, and for mono-, bi-, and tri-disperse mixture simulations of Apollo LM plume-surface interaction. The computational results confirm the significance of including the particle shape and mixture effects. Going forward plans for the full implementation of the general poly-disperse regolith modeling capability and maturation towards NASA project application readiness under the GCD program will be presented.

Peter A. Liever↗

NASA NTRS1987

Evaluation of a theory for pressure-strain rate

A theoretical expression for the slow part (the nonlinear fluctuation part) of the pressure-strain rate is compared with simulations of anisotropic homogeneous flows. The objective is to determine the quantitative accuracy of the theory and to test its prediction that the generalized Rotta coefficient, a non-dimensionalized ratio of slow term to the Reynolds stress anisotropy, varies with direction and can be negative. Comparisons are made between theoretical and simulation values of the slow term itself and of the generalized Rotta coefficients. The implications of the comparison for two-point closure theories and for Reynolds stress modeling are pointed out.

Weinstock, J.↗

NASA NTRS1990

Mathematical theory of a relaxed design problem in structural optimization

Various attempts have been made to construct a rigorous mathematical theory of optimization for size, shape, and topology (i.e. layout) of an elastic structure. If these are represented by a finite number of parametric functions, as Armand described, it is possible to construct an existence theory of the optimum design using compactness argument in a finite dimensional design space or a closed admissible set of a finite dimensional design space. However, if the admissible design set is a subset of non-reflexive Banach space such as L(sup infinity)(Omega), construction of the existence theory of the optimum design becomes suddenly difficult and requires to extend (i.e. generalize) the design problem to much more wider class of design that is compatible to mechanics of structures in the sense of variational principle. Starting from the study by Cheng and Olhoff, Lurie, Cherkaev, and Fedorov introduced a new concept of convergence of design variables in a generalized sense and construct the 'G-Closure' theory of an extended (relaxed) optimum design problem. A similar attempt, but independent in large extent, can also be found in Kohn and Strang in which the shape and topology optimization problem is relaxed to allow to use of perforated composites rather than restricting it to usual solid structures. An identical idea is also stated in Murat and Tartar using the notion of the homogenization theory. That is, introducing possibility of micro-scale perforation together with the theory of homogenization, the optimum design problem is relaxed to construct its mathematical theory. It is also noted that this type of relaxed design problem is perfectly matched to the variational principle in structural mechanics.

Kikuchi, Noboru↗

NASA NTRS2006

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.

Rubinstein, Robert↗

NASA NTRS2001

Discontinuous Galerkin Methods for NonLinear Differential Systems

This talk considers simplified finite element discretization techniques for first-order systems of conservation laws equipped with a convex (entropy) extension. Using newly developed techniques in entropy symmetrization theory, simplified forms of the discontinuous Galerkin (DG) finite element method have been developed and analyzed. The use of symmetrization variables yields numerical schemes which inherit global entropy stability properties of the PDE (partial differential equation) system. Central to the development of the simplified DG methods is the Eigenvalue Scaling Theorem which characterizes right symmetrizers of an arbitrary first-order hyperbolic system in terms of scaled eigenvectors of the corresponding flux Jacobian matrices. A constructive proof is provided for the Eigenvalue Scaling Theorem with detailed consideration given to the Euler equations of gas dynamics and extended conservation law systems derivable as moments of the Boltzmann equation. Using results from kinetic Boltzmann moment closure theory, we then derive and prove energy stability for several approximate DG fluxes which have practical and theoretical merit.

Barth, Timothy↗

NASA NTRS1995

Analytical and phenomenological studies of rotating turbulence

A framework, which combines mathematical analysis, closure theory, and phenomenological treatment, is developed to study the spectral transfer process and reduction of dimensionality in turbulent flows that are subject to rotation. First, we outline a mathematical procedure that is particularly appropriate for problems with two disparate time scales. The approach which is based on the Green's method leads to the Poincare velocity variables and the Poincare transformation when applied to rotating turbulence. The effects of the rotation are now reflected in the modifications to the convolution of a nonlinear term. The Poincare transformed equations are used to obtain a time-dependent analog of the Taylor-Proudman theorem valid in the asymptotic limit when the non-dimensional parameter mu is identical to Omega(t) approaches infinity (Omega is the rotation rate and t is the time). The 'split' of the energy transfer in both direct and inverse directions is established. Secondly, we apply the Eddy-Damped-Quasinormal-Markovian (EDQNM) closure to the Poincare transformed Euler/Navier-Stokes equations. This closure leads to expressions for the spectral energy transfer. In particular, an unique triple velocity decorrelation time is derived with an explicit dependence on the rotation rate. This provides an important input for applying the phenomenological treatment of Zhou. In order to characterize the relative strength of rotation, another non-dimensional number, a spectral Rossby number, which is defined as the ratio of rotation and turbulence time scales, is introduced. Finally, the energy spectrum and the spectral eddy viscosity are deduced.

Mahalov, Alex↗

NASA NTRS1977

A model for dispersion from area sources in convective turbulence

Four independent estimates of the vertical distribution of the eddy coefficient for dispersion of a passive contaminant from an extensive area source in a convective layer have been presented. The estimates were based on the following methods: (1) a second-order closure prediction, (2) field data of pollutant concentrations over Los Angeles, (3) lab measurements of particle dispersion, and (4) assumption of equality between momentum and mass transfer coefficients in the free convective limit. It is suggested that K-values estimated both from second-order closure theory and from Los Angeles measurements are systematically underestimated.

Crane, G.↗

NASA NTRS1991

Modeling turbulent compressible flows - The mass fluctuating velocity and squared density

This paper deals with single-point closure theory for compressible turbulent flow, including the effects of compressibility on the turbulence. In particular, the combination of the pressure dilatation and the dilatation dissipation, terms which appear on the turbulent kinetic energy equation, are modeled. Model parameters in these transport equations are determined by comparing predictions with boundary layer measurements. Finally, predictions with a k-epsilon model, including the new formulations, are presented for the compressible shear layer.

Taulbee, D.↗

NASA NTRS1993

Renormalization group estimates of transport coefficients in the advection of a passive scalar by incompressible turbulence

The advection of a passive scalar by incompressible turbulence is considered using recursive renormalization group procedures in the differential sub grid shell thickness limit. It is shown explicitly that the higher order nonlinearities induced by the recursive renormalization group procedure preserve Galilean invariance. Differential equations, valid for the entire resolvable wave number k range, are determined for the eddy viscosity and eddy diffusivity coefficients, and it is shown that higher order nonlinearities do not contribute as k goes to 0, but have an essential role as k goes to k(sub c) the cutoff wave number separating the resolvable scales from the sub grid scales. The recursive renormalization transport coefficients and the associated eddy Prandtl number are in good agreement with the k-dependent transport coefficients derived from closure theories and experiments.

Zhou, YE↗

NASA NTRS1993

A boundary-layer model for Mars - Comparison with Viking lander and entry data

A 1D boundary-layer model of Mars based on a momentum equation that describes friction, pressure gradient, and Coriolis forces is presented. Frictional forces and convective heating are computed using the level-2 turbulence closure theory of Mellor and Yamada (1974). The model takes into account the radiative effects of CO2 gas and suspended dust particles. Both radiation and convection depend on surface temperatures which are computed from a surface heat budget. Model predictions are compared with available observations from Viking landers. It is concluded that, in general, the model reproduces the basic features of the temperature data. The agreement is particularly good at entry time for the V L-2 site, where the model and observations are within several degrees at all levels for which data are available.

Haberle, Robert M.↗

NASA NTRS1993

An investigation of small scales of turbulence in a boundary layer at high Reynolds numbers

The assumption that turbulence at large wave-numbers is isotropic and has universal spectral characteristics which are independent of the flow geometry, at least for high Reynolds numbers, has been a cornerstone of closure theories as well as of the most promising recent development in the effort to predict turbulent flows, viz. large eddy simulations. This hypothesis was first advanced by Kolmogorov based on the supposition that turbulent kinetic energy cascades down the scales (up the wave-numbers) of turbulence and that, if the number of these cascade steps is sufficiently large (i.e. the wave-number range is large), then the effects of anisotropies at the large scales are lost in the energy transfer process. Experimental attempts were repeatedly made to verify this fundamental assumption. However, Van Atta has recently suggested that an examination of the scalar and velocity gradient fields is necessary to definitively verify this hypothesis or prove it to be unfounded. Of course, this must be carried out in a flow with a sufficiently high Reynolds number to provide the necessary separation of scales in order unambiguously to provide the possibility of local isotropy at large wave-numbers. An opportunity to use our 12-sensor hot-wire probe to address this issue directly was made available at the 80'x120' wind tunnel at the NASA Ames Research Center, which is normally used for full-scale aircraft tests. An initial report on this high Reynolds number experiment and progress toward its evaluation is presented.

Wallace, James M.↗

NASA NTRS1994

Development of the NASA/FLAGRO computer program for analysis of airframe structures

The NASA/FLAGRO (NASGRO) computer program was developed for fracture control analysis of space hardware and is currently the standard computer code in NASA, the U.S. Air Force, and the European Agency (ESA) for this purpose. The significant attributes of the NASGRO program are the numerous crack case solutions, the large materials file, the improved growth rate equation based on crack closure theory, and the user-friendly promptive input features. In support of the National Aging Aircraft Research Program (NAARP); NASGRO is being further developed to provide advanced state-of-the-art capability for damage tolerance and crack growth analysis of aircraft structural problems, including mechanical systems and engines. The project currently involves a cooperative development effort by NASA, FAA, and ESA. The primary tasks underway are the incorporation of advanced methodology for crack growth rate retardation resulting from spectrum loading and improved analysis for determining crack instability. Also, the current weight function solutions in NASGRO or nonlinear stress gradient problems are being extended to more crack cases, and the 2-d boundary integral routine for stress analysis and stress-intensity factor solutions is being extended to 3-d problems. Lastly, effort is underway to enhance the program to operate on personal computers and work stations in a Windows environment. Because of the increasing and already wide usage of NASGRO, the code offers an excellent mechanism for technology transfer for new fatigue and fracture mechanics capabilities developed within NAARP.

Forman, R. G.↗

NASA NTRS1996

Crack Development in Cross-Ply Laminates Under Uniaxial Tension

This study addresses matrix-dominated failures in carbon fiber/polymer matrix composite laminates in a cross-ply lay-up. The events of interest are interlaminar fracture in the form of transverse cracks in the 90' plies and longitudinal splitting in the 0 deg plies and interlaminar fracture in the form of 0/90 delamination. These events were observed using various nondestructive evaluation (NDE) techniques during static tensile tests. Acoustic emission (AE) x radiography, and edge view microscopy were the principal ones utilized in a real-time environment. A comparison of the NDE results with an analytical model based on the classical linear fracture mechanics concept of strain energy release rate as a criterion for crack growth was performed. The virtual crack closure theory was incorporated with a finite element model to generate strain energy release rate curves for the analytical case. Celion carbon fiber/polyimide matrix (G30-500/PMR-15) was the material tested with cross-ply lay-ups of (0(2)/90(6))s and (0(4)/90(4))s. The test specimens contained thermally induced cracks caused by the high-temperature processing. The analytical model was updated to compensate for the initial damage and to study further accumulation by taking into account the crack interactions. By correlating the experimental and analytical data, the critical energy release rates were found for the observable events of interest.

Gyekenyesi, Andrew L.↗

NASA NTRS1997

Energy Transfer and Triadic Interactions in Compressible Turbulence

Using a two-point closure theory, the Eddy-Damped-Quasi-Normal-Markovian (EDQNM) approximation, we have investigated the energy transfer process and triadic interactions of compressible turbulence. In order to analyze the compressible mode directly, the Helmholtz decomposition is used. The following issues were addressed: (1) What is the mechanism of energy exchange between the solenoidal and compressible modes, and (2) Is there an energy cascade in the compressible energy transfer process? It is concluded that the compressible energy is transferred locally from the solenoidal part to the compressible part. It is also found that there is an energy cascade of the compressible mode for high turbulent Mach number (M(sub t) greater than or equal to 0.5). Since we assume that the compressibility is weak, the magnitude of the compressible (radiative or cascade) transfer is much smaller than that of solenoidal cascade. These results are further confirmed by studying the triadic energy transfer function, the most fundamental building block of the energy transfer.

Bataille, F.↗

NASA NTRS2011

Two-Point Turbulence Closure Applied to Variable Resolution Modeling

Variable resolution methods have become frontline CFD tools, but in order to take full advantage of this promising new technology, more formal theoretical development is desirable. Two general classes of variable resolution methods can be identified: hybrid or zonal methods in which RANS and LES models are solved in different flow regions, and bridging or seamless models which interpolate smoothly between RANS and LES. This paper considers the formulation of bridging methods using methods of two-point closure theory. The fundamental problem is to derive a subgrid two-equation model. We compare and reconcile two different approaches to this goal: the Partially Integrated Transport Model, and the Partially Averaged Navier-Stokes method.

Girimaji, Sharath S.↗