Search NASA⌕ Search

Engineering topics

Rubinstein, R.

Publications and source records attributed to Rubinstein, R..

13 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS2001

Effects of Eddy Viscosity on Time Correlations in Large Eddy Simulation

Subgrid-scale (SGS) models for large. eddy simulation (LES) have generally been evaluated by their ability to predict single-time statistics of turbulent flows such as kinetic energy and Reynolds stresses. Recent application- of large eddy simulation to the evaluation of sound sources in turbulent flows, a problem in which time, correlations determine the frequency distribution of acoustic radiation, suggest that subgrid models should also be evaluated by their ability to predict time correlations in turbulent flows. This paper compares the two-point, two-time Eulerian velocity correlation evaluated from direct numerical simulation (DNS) with that evaluated from LES, using a spectral eddy viscosity, for isotropic homogeneous turbulence. It is found that the LES fields are too coherent, in the sense that their time correlations decay more slowly than the corresponding time. correlations in the DNS fields. This observation is confirmed by theoretical estimates of time correlations using the Taylor expansion technique. Tile reason for the slower decay is that the eddy viscosity does not include the random backscatter, which decorrelates fluid motion at large scales. An effective eddy viscosity associated with time correlations is formulated, to which the eddy viscosity associated with energy transfer is a leading order approximation.

He, Guowei↗
NASA NTRS2001

Turbulence Modeling Workshop

Advances in turbulence modeling are needed in order to calculate high Reynolds number flows near the onset of separation and beyond. To this end, the participants in this workshop made the following recommendations. (1) A national/international database and standards for turbulence modeling assessment should be established. Existing experimental data sets should be reviewed and categorized. Advantage should be taken of other efforts already under-way, such as that of the European Research Community on Flow, Turbulence, and Combustion (ERCOFTAC) consortium. Carefully selected "unit" experiments will be needed, as well as advances in instrumentation, to fill the gaps in existing datasets. A high priority should be given to document existing turbulence model capabilities in a standard form, including numerical implementation issues such as grid quality and resolution. (2) NASA should support long-term research on Algebraic Stress Models and Reynolds Stress Models. The emphasis should be placed on improving the length-scale equation, since it is the least understood and is a key component of two-equation and higher models. Second priority should be given to the development of improved near-wall models. Direct Numerical Simulations (DNS) and Large Eddy Simulations (LES) would provide valuable guidance in developing and validating new Reynolds-averaged Navier-Stokes (RANS) models. Although not the focus of this workshop, DNS, LES, and hybrid methods currently represent viable approaches for analysis on a limited basis. Therefore, although computer limitations require the use of RANS methods for realistic configurations at high Reynolds number in the foreseeable future, a balanced effort in turbulence modeling development, validation, and implementation should include these approaches as well.

Rubinstein, R.↗
NASA NTRS2000

Mapping Closure Approximation to Conditional Dissipation Rate for Turbulent Scalar Mixing

A novel mapping closure approximation (MCA) technique is developed to construct a model for the conditional dissipation rate (CDR) of a scalar in homogeneous turbulence. It is shown that the CDR model from amplitude mapping closure is incorrect in asymptotic behavior for unsymmetric binary mixing. The correct asymptotic behavior can be described by the CDR model formulated by the MCA technique. The MCA approach is outlined for constructing successive approximation to probability density function (PDF) and conditional moment.

He, Gouwei↗
NASA NTRS1993

Time dependent turbulence modeling and analytical theories of turbulence

By simplifying the direct interaction approximation (DIA) for turbulent shear flow, time dependent formulas are derived for the Reynolds stresses which can be included in two equation models. The Green's function is treated phenomenologically, however, following Smith and Yakhot, we insist on the short and long time limits required by DIA. For small strain rates, perturbative evaluation of the correlation function yields a time dependent theory which includes normal stress effects in simple shear flows. From this standpoint, the phenomenological Launder-Reece-Rodi model is obtained by replacing the Green's function by its long time limit. Eddy damping corrections to short time behavior initiate too quickly in this model; in contrast, the present theory exhibits strong suppression of eddy damping at short times. A time dependent theory for large strain rates is proposed in which large scales are governed by rapid distortion theory while small scales are governed by Kolmogorov inertial range dynamics. At short times and large strain rates, the theory closely matches rapid distortion theory, but at long times it relaxes to an eddy damping model.

Rubinstein, R.↗
NASA NTRS1992

Renormalization group methods for the Reynolds stress transport equations

The Yakhot-Orszag renormalization group is used to analyze the pressure gradient-velocity correlation and return to isotropy terms in the Reynolds stress transport equations. The perturbation series for the relevant correlations, evaluated to lowest order in the epsilon-expansion of the Yakhot-Orszag theory, are infinite series in tensor product powers of the mean velocity gradient and its transpose. Formal lowest order Pade approximations to the sums of these series produce a rapid pressure strain model of the form proposed by Launder, Reece, and Rodi, and a return to isotropy model of the form proposed by Rotta. In both cases, the model constants are computed theoretically. The predicted Reynolds stress ratios in simple shear flows are evaluated and compared with experimental data. The possibility is discussed of deriving higher order nonlinear models by approximating the sums more accurately. The Yakhot-Orszag renormalization group provides a systematic procedure for deriving turbulence models. Typical applications have included theoretical derivation of the universal constants of isotropic turbulence theory, such as the Kolmogorov constant, and derivation of two equation models, again with theoretically computed constants and low Reynolds number forms of the equations. Recent work has applied this formalism to Reynolds stress modeling, previously in the form of a nonlinear eddy viscosity representation of the Reynolds stresses, which can be used to model the simplest normal stress effects. The present work attempts to apply the Yakhot-Orszag formalism to Reynolds stress transport modeling.

Rubinstein, R.↗
NASA NTRS1992

Renormalization group analysis of the Reynolds stress transport equation

The pressure velocity correlation and return to isotropy term in the Reynolds stress transport equation are analyzed using the Yakhot-Orszag renormalization group. The perturbation series for the relevant correlations, evaluated to lowest order in the epsilon-expansion of the Yakhot-Orszag theory, are infinite series in tensor product powers of the mean velocity gradient and its transpose. Formal lowest order Pade approximations to the sums of these series produce a fast pressure strain model of the form proposed by Launder, Reece, and Rodi, and a return to isotropy model of the form proposed by Rotta. In both cases, the model constant are computed theoretically. The predicted Reynolds stress ratios in simple shear flows are evaluated and compared with experimental data. The possibility is discussed of driving higher order nonlinear models by approximating the sums more accurately.

Rubinstein, R.↗
NASA NTRS1991

Nonlinear Reynolds stress model for turbulent shear flows

A nonlinear algebraic Reynolds stress model, derived using the renormalization group, is applied to equilibrium homogeneous shear flow and fully developed flow in a square duct. The model, which is quadratically nonlinear in the velocity gradients, successfully captures the large-scale inhomogeneity and anisotropy of the flows studied. The ratios of normal stresses, as well as the actual magnitudes of the stresses are correctly predicted for equilibrium homogeneous shear flow. Reynolds normal stress anisotropy and attendant turbulence driven secondary flow are predicted for a square duct. Profiles of mean velocity and normal stresses are in good agreement with measurements. Very close to walls, agreement with measurements diminishes. The model has the benefit of containing no arbitrary constants; all values are determined directly from the theory. It seems that near wall behavior is influenced by more than the large scale anisotropy accommodated in the current model. More accurate near wall calculations may well require a model for anisotropic dissipation.

Barton, J. Michael↗
NASA NTRS1991

Young's fringe particle image velocimetry in the presence of random particle motions

In the Young's fringe approach to particle image velocimetry, random particle motions cause loss of fringe visibility by decorrelating the coherent interexposure particle separations responsible for fringe formation. Since the visibility reduction is determined by the random motion through a Fourier transform relation analogous to the Van Cittert Zernike Theorem, it has been proposed that the random motion can be characterized statistically by analyzing fringe visibility distributions in the transform plane. This paper assesses the accuracy of such measurements. In particular, the effects of finite particle population and of correlated random motion are evaluated. The theory applies to diffusive motion, turbulence, and to random motion caused by mean velocity inhomogeneities.

Rubinstein, R.↗
NASA NTRS1988

STAEBL/General composites with hygrothermal effects (STAEBL/GENCOM)

A computer code has been developed to perform structural optimization of turbine blades made from angle ply fiber composite laminates. Design variables available for optimization include geometric parameters such as blade thickness disribution and root chord, and composite material parameters such as ply angles and numbers of plies of each constituent material. Design constraints include resonance margins, forced response margins, maximum stress, and maximum ply combined stress. A general description of this code is given. Design optimization studies for typical blades are presented.

Rubinstein, R.↗
NASA NTRS1986

Structural tailoring of SSME turbopump blades (SSME/STAEBL)

Computer structural optimization is applied to the design of Space Shuttle main engine turbopump blades. The optimization is implemented by the program SSME/STAEBL. A general description of this program is given. Design optimization studies for typical blade designs are presented.

Rubinstein, R.↗
NASA NTRS1985

Structural Tailoring of SSME Blades (vanes)

The engine blade design optimization program STAEBL (Structural Tailoring of Engine Blades) is available at the NASA Lewis computer facility. The analysis capabilities of this program were extended to typical loading conditions for SSME turbopump blades including thermal and pressure loading. Input files for representative SSME blade designs were developed and sample optimization studies for these blades completed. The structural tailoring program combines a general optimization package and a finite element blade analysis package. The analysis package's capabilities include natural frequency, maximum stress, and forced response computation, and fatigue life and flutter analysis. Optimization is performed using the feasible directions method. The current design is modified by perturbing the design variables so that the design constraints are satisfied while the objective function, such as blade weight, is reduced at the maximum rate. The program's geometric design variables include blade thickness distribution, thickness to chord ratios, and root chord.

Rubinstein, R.↗