scholarly journals A similarity theory of locally homogeneous and isotropic turbulence generated by a Smagorinsky-type LES

1996 ◽  
Vol 325 ◽  
pp. 239-260 ◽  
Author(s):  
Andreas Muschinski
Keyword(s):  
Newtonian Fluid ◽  
Material Parameter ◽  
Isotropic Turbulence ◽  
Similarity Theory ◽  
Reynolds Numbers ◽  
Finite Difference Approximation ◽  
Difference Approximation ◽  
Spectral Models ◽  
Homogeneous And Isotropic Turbulence ◽  
Locally Homogeneous

A Kolmogorov-type similarity theory of locally homogeneous and isotropic turbulence generated by a Smagorinsky-type large-eddy simulation (LES) at very large LES Reynolds numbers is developed and discussed. The underlying concept is that the LES equations may be considered equations of motion of specific hypothetical fully turbulent non-Newtonian fluids, called ‘LES fluids’. It is shown that the length scale ls = csδ, which scales the magnitude of the variable viscosity in a Smagorinsky-type LES, is the ‘Smagorinsky-fluid’ counterpart of Kolmogorov's dissipation length $\eta = v^{3/4}\epsilon^{-1/4}$ for a Newtonian fluid where ν is the kinematic viscosity and ε is the energy dissipation rate. While in a Newtonian fluid the viscosity is a material parameter and the length ν depends on ε, in a Smagorinsky fluid the length ls is a material parameter and the viscosity depends on ε. The Smagorinsky coefficient cs may be considered the reciprocal of a ‘microstructure Knudsen number’ of a Smagorinsky fluid. A combination of Lilly's (1967) cut-off model with two well-known spectral models for dissipation-range turbulence (Heisenberg 1948; Pao 1965) leads to models for the LES-generated Kolmogorov coefficient αLES as a function of cs. Both models predict an intrinsic overestimation of αLES for finite values of cs. For cs = 0.2 Heisenberg's and Pao's models provide αLES = 1.74 (16% overestimation) and αLES = 2.14 (43% overestimation), respectively, if limcs → (αLES) = 1.5 is ad hoc assumed. The predicted overestimation becomes negligible beyond about cs = 0.5. The requirement cs > 0.5 is equivalent to Δ < 2ls. A similar requirement, L < 2η where L is the wire length of hot-wire anemometers, has been recommended by experimentalists. The value of limcs → (αLES) for a Smagorinsky-type LES at very large LES Reynolds numbers is not predicted by the models and remains unknown. Two critical values of cs are identified. The first critical cs is Lilly's (1967) value, which indicates the cs below which finite-difference-approximation errors become important; the second critical cs is the value beyond which the Reynolds number similarity is violated.

The Material Parameter Identification for Functionally Graded Materials by the Isoparametric Graded Finite Element

2012 ◽  
Vol 446-449 ◽  
pp. 3609-3614 ◽  
Author(s):  
Li Xin Huang ◽  
Lin Wang ◽  
Yue Chen ◽  
Qi Yao ◽  
Xiao Jun Zhou
Keyword(s):  
Parameter Identification ◽  
Functionally Graded Materials ◽  
Material Parameter ◽  
Identification Problem ◽  
Optimization Method ◽  
Functionally Graded ◽  
Finite Difference Approximation ◽  
Difference Approximation ◽  
Material Parameter Identification ◽  
Graded Materials

A material parameter identification method is proposed for functionally graded materials (FGMs) which are modeled by the isoparametric graded finite elements (IGFE). The material parameter identification problem is formulated as the problem of minimizing the objective function defined as a square sum of differences between the measured displacement and the computed displacement by the IGFE. Levenberg-Marquardt optimization method, in which the sensitivity analysis of displacements with respect to the material parameters is based on the finite difference approximation method, is used to solve the minimization problem. Numerical example is given to illustrate the validity of the proposed method for parameter identification.

scholarly journals Design and Validation of a Probe for Spatially and Temporally Resolved Measurements of Vorticity and Strain Rates in Compressible Turbulence Interactions

10.14311/994 ◽  
2007 ◽  
Vol 47 (6) ◽  
Author(s):  
S. Xanthos ◽  
M. Gong ◽  
Y. Andreopoulos
Keyword(s):  
Isotropic Turbulence ◽  
Strain Tensor ◽  
Three Dimensional ◽  
Reynolds Numbers ◽  
Compressible Turbulence ◽  
Expansion Waves ◽  
Custom Made ◽  
Homogeneous And Isotropic Turbulence ◽  
Total Temperature ◽  
Kinetic Energy Term

A custom-made hot-wire vorticity probe was designed and developed capable of measuring the time-dependent highly fluctuating three dimensional velocity and vorticity vectors, and associated total temperature, in non-isothermal and inhomogeneous flows with reasonable spatial and temporal resolution. These measurements allowed computation of the vorticity stretching/tilting terms, vorticity generation through dilatation terms, full dissipation rate of the kinetic energy term and full rate-of-strain tensor. The probe has been validated experimentally in low-speed boundary layers and used in the CCNY Shock Tube Research Facility, where interactions of planar expansion waves or shock waves with homogeneous and isotropic turbulence have been investigated at several Reynolds numbers. 

scholarly journals Stable and Efficient Modeling of Anelastic Attenuation in Seismic Wave Propagation

2012 ◽  
Vol 12 (1) ◽  
pp. 193-225 ◽  
Author(s):  
N. Anders Petersson ◽  
Björn Sjögreen
Keyword(s):  
Wave Equation ◽  
Finite Difference ◽  
Half Space ◽  
Material Model ◽  
Second Order ◽  
Finite Difference Approximation ◽  
Difference Approximation ◽  
Space Problem ◽  
Standard Linear Solid ◽  
Standard Linear

AbstractWe develop a stable finite difference approximation of the three-dimensional viscoelastic wave equation. The material model is a super-imposition of N standard linear solid mechanisms, which commonly is used in seismology to model a material with constant quality factor Q. The proposed scheme discretizes the governing equations in second order displacement formulation using 3N memory variables, making it significantly more memory efficient than the commonly used first order velocity-stress formulation. The new scheme is a generalization of our energy conserving finite difference scheme for the elastic wave equation in second order formulation [SIAM J. Numer. Anal., 45 (2007), pp. 1902-1936]. Our main result is a proof that the proposed discretization is energy stable, even in the case of variable material properties. The proof relies on the summation-by-parts property of the discretization. The new scheme is implemented with grid refinement with hanging nodes on the interface. Numerical experiments verify the accuracy and stability of the new scheme. Semi-analytical solutions for a half-space problem and the LOH.3 layer over half-space problem are used to demonstrate how the number of viscoelastic mechanisms and the grid resolution influence the accuracy. We find that three standard linear solid mechanisms usually are sufficient to make the modeling error smaller than the discretization error.

On the time scales and structure of Lagrangian intermittency in homogeneous isotropic turbulence

2019 ◽  
Vol 867 ◽  
pp. 438-481 ◽  
Author(s):  
R. Watteaux ◽  
G. Sardina ◽  
L. Brandt ◽  
D. Iudicone
Keyword(s):  
Time Scales ◽  
Isotropic Turbulence ◽  
Flow Characteristics ◽  
Residence Times ◽  
Particle Trajectories ◽  
Local Flow ◽  
Available Energy ◽  
Homogeneous And Isotropic Turbulence ◽  
History Of ◽  
Optimum Manner

We present a study of Lagrangian intermittency and its characteristic time scales. Using the concepts of flying and diving residence times above and below a given threshold in the magnitude of turbulence quantities, we infer the time spectra of the Lagrangian temporal fluctuations of dissipation, acceleration and enstrophy by means of a direct numerical simulation in homogeneous and isotropic turbulence. We then relate these time scales, first, to the presence of extreme events in turbulence and, second, to the local flow characteristics. Analyses confirm the existence in turbulent quantities of holes mirroring bursts, both of which are at the core of what constitutes Lagrangian intermittency. It is shown that holes are associated with quiescent laminar regions of the flow. Moreover, Lagrangian holes occur over few Kolmogorov time scales while Lagrangian bursts happen over longer periods scaling with the global decorrelation time scale, hence showing that loss of the history of the turbulence quantities along particle trajectories in turbulence is not continuous. Such a characteristic partially explains why current Lagrangian stochastic models fail at reproducing our results. More generally, the Lagrangian dataset of residence times shown here represents another manner for qualifying the accuracy of models. We also deliver a theoretical approximation of mean residence times, which highlights the importance of the correlation between turbulence quantities and their time derivatives in setting temporal statistics. Finally, whether in a hole or a burst, the straining structure along particle trajectories always evolves self-similarly (in a statistical sense) from shearless two-dimensional to shear bi-axial configurations. We speculate that this latter configuration represents the optimum manner to dissipate locally the available energy.

Unit Mobility Ratio Displacement Calculations for Pattern Floods in Homogeneous Medium

1966 ◽  
Vol 6 (03) ◽  
pp. 217-227 ◽  
Author(s):  
Hubert J. Morel-Seytoux
Keyword(s):  
Oil Recovery ◽  
Numerical Solutions ◽  
Mobility Ratio ◽  
Finite Difference Approximation ◽  
Difference Approximation ◽  
Numerical Techniques ◽  
Sweep Efficiency ◽  
Closed Form Solutions ◽  
Regular Patterns ◽  
Quantities Of Interest

Abstract The influence of pattern geometry on assisted oil recovery for a particular displacement mechanism is the object of investigation in this paper. The displacement is assumed to be of unit mobility ratio and piston-like. Fluids are assumed incompressible and gravity and capillary effects are neglected. With these assumptions it is possible to calculate by analytical methods the quantities of interest to the reservoir engineer for a great variety of patterns. Specifically, this paper presentsvery briefly, the methods and mathematical derivations required to obtain the results of engineering concern, andtypical results in the form of graphs or formulae that can be used readily without prior study of the methods. Results of this work provide checks for solutions obtained from programmed numerical techniques. They also reveal the effect of pattern geometry and, even though the assumptions of piston-like displacement and of unit mobility ratio are restrictive, they can nevertheless be used for rather crude but quick, cheap estimates. These estimates can be refined to account for non-unit mobility ratio and two-phase flow by correlating analytical results in the case M=1 and the numerical results for non-Piston, non-unit mobility ratio displacements. In an earlier paper1 it was also shown that from the knowledge of closed form solutions for unit mobility ratio, quantities called "scale factors" could be readily calculated, increasing considerably the flexibility of the numerical techniques. Many new closed form solutions are given in this paper. INTRODUCTION BACKGROUND Pattern geometry is a major factor in making water-flood recovery predictions. For this reason many numerical schemes have been devised to predict oil recovery in either regular patterns or arbitrary configurations. The numerical solutions, based on the method of finite difference approximation, are subject to errors often difficult to evaluate. An estimate of the error is possible by comparison with exact solutions. To provide a variety of checks on numerical solutions, a thorough study of the unit mobility ratio displacement process was undertaken. To calculate quantities of interest to the reservoir engineer (oil recovery, sweep efficiency, etc.), it is necessary to first know the pressure distribution in the pattern. Then analytical procedures are used to calculate, in order of increasing difficulty: injectivity, breakthrough areal sweep efficiency, normalized oil recovery and water-oil ratio as a function of normalized PV injected. BACKGROUND Pattern geometry is a major factor in making water-flood recovery predictions. For this reason many numerical schemes have been devised to predict oil recovery in either regular patterns or arbitrary configurations. The numerical solutions, based on the method of finite difference approximation, are subject to errors often difficult to evaluate. An estimate of the error is possible by comparison with exact solutions. To provide a variety of checks on numerical solutions, a thorough study of the unit mobility ratio displacement process was undertaken. To calculate quantities of interest to the reservoir engineer (oil recovery, sweep efficiency, etc.), it is necessary to first know the pressure distribution in the pattern. Then analytical procedures are used to calculate, in order of increasing difficulty: injectivity, breakthrough areal sweep efficiency, normalized oil recovery and water-oil ratio as a function of normalized PV injected.

Experiments of Compressible Homogeneous and Isotropic Turbulence Interacting With Expansion Waves

2003 ◽  
Author(s):  
Savvas S. Xanthos ◽  
Yiannis Andreopoulos
Keyword(s):  
Mach Number ◽  
Shock Tube ◽  
Total Pressure ◽  
Large Scale ◽  
Isotropic Turbulence ◽  
Pressure Fluctuations ◽  
Incident Shock ◽  
Expansion Waves ◽  
Curvature Effects ◽  
Homogeneous And Isotropic Turbulence

The interaction of traveling expansion waves with grid-generated turbulence was investigated in a large-scale shock tube research facility. The incident shock and the induced flow behind it passed through a rectangular grid, which generated a nearly homogeneous and nearly isotropic turbulent flow. As the shock wave exited the open end of the shock tube, a system of expansion waves was generated which traveled upstream and interacted with the grid-generated turbulence; a type of interaction free from streamline curvature effects, which cause additional effects on turbulence. In this experiment, wall pressure, total pressure and velocity were measured indicating a clear reduction in fluctuations. The incoming flow at Mach number 0.46 was expanded to a flow with Mach number 0.77 by an applied mean shear of 100 s−1. Although the strength of the generated expansion waves was mild, the effect on damping fluctuations on turbulence was clear. A reduction of in the level of total pressure fluctuations by 20 per cent was detected in the present experiments.

Sign in / Sign up

Export Citation Format

Share Document