Subgrid scale fluid velocity timescales seen by inertial particles in large-eddy simulation of particle-laden turbulence

Large-eddy simulation (LES) has emerged as a promising tool for simulating turbulent flows in general and, in recent years,has also been applied to the particle-laden turbulence with some success (Kassinos et al., 2007). The motion of inertial particles is much more complicated than fluid elements, and therefore, LES of turbulent flow laden with inertial particles encounters new challenges. In the conventional LES, only large-scale eddies are explicitly resolved and the effects of unresolved, small or subgrid scale (SGS) eddies on the large-scale eddies are modeled. The SGS turbulent flow field is not available. The effects of SGS turbulent velocity field on particle motion have been studied by Wang and Squires (1996), Armenio et al. (1999), Yamamoto et al. (2001), Shotorban and Mashayek (2006a,b), Fede and Simonin (2006), Berrouk et al. (2007), Bini and Jones (2008), and Pozorski and Apte (2009), amongst others. One contemporary method to include the effects of SGS eddies on inertial particle motions is to introduce a stochastic differential equation (SDE), that is, a Langevin stochastic equation to model the SGS fluid velocity seen by inertial particles (Fede et al., 2006; Shotorban and Mashayek, 2006a; Shotorban and Mashayek, 2006b; Berrouk et al., 2007; Bini and Jones, 2008; Pozorski and Apte, 2009).However, the accuracy of such a Langevin equation model depends primarily on the prescription of the SGS fluid velocity autocorrelation time seen by an inertial particle or the inertial particle–SGS eddy interaction timescale (denoted by $\delt T_{Lp}$ and a second model constant in the diffusion term which controls the intensity of the random force received by an inertial particle (denoted by C_0, see Eq. (7)). From the theoretical point of view, dTLp differs significantly from the Lagrangian fluid velocity correlation time (Reeks, 1977; Wang and Stock, 1993), and this carries the essential nonlinearity in the statistical modeling of particle motion. dTLp and C0 may depend on the filter width and particle Stokes number even for a given turbulent flow. In previous studies, dTLp is modeled either by the fluid SGS Lagrangian timescale (Fede et al., 2006; Shotorban and Mashayek, 2006b; Pozorski and Apte, 2009; Bini and Jones, 2008) or by a simple extension of the timescale obtained from the full flow field (Berrouk et al., 2007). In this work, we shall study the subtle and on-monotonic dependence of $\delt T_{Lp}$ on the filter width and particle Stokes number using a flow field obtained from Direct Numerical Simulation (DNS). We then propose an empirical closure model for $\delta T_{Lp}$. Finally, the model is validated against LES of particle-laden turbulence in predicting single-particle statistics such as particle kinetic energy. As a first step, we consider the particle motion under the one-way coupling assumption in isotropic turbulent flow and neglect the gravitational settling effect. The one-way coupling assumption is only valid for low particle mass loading.

Effects of Subgrid-Scale Modeling on Time Correlations in Large Eddy Simulation

The effects of the unresolved subgrid-scale (SGS) motions on the energy balance of the resolved scales in large eddy simulation (LES) have been investigated actively because modeling the energy transfer between the resolved and unresolved scales is crucial to constructing accurate SGS models. But the subgrid scales not only modify the energy balance, they also contribute to temporal decorrelation of the resolved scales. The importance of this effect in applications including the predictability problem and the evaluation of sound radiation by turbulent flows motivates the present study of the effect of SGS modeling on turbulent time correlations. This paper compares the two-point, two-time Eulerian velocity correlation in isotropic homogeneous turbulence evaluated by direct numerical simulation (DNS) with the correlations evaluated by LES using a standard spectral eddy viscosity. It proves convenient to express the two-point correlations in terms of spatial Fourier decomposition of the velocity field. The LES fields are more coherent than the DNS fields: their time correlations decay more slowly at all resolved scales of motion and both their integral scales and microscales are larger than those of the DNS field. Filtering alone is not responsible for this effect: in the Fourier representation, the time correlations of the filtered DNS field are identical to those of the DNS field itself. The possibility of modeling the decorrelating effects of the unresolved scales of motion by including a random force in the model is briefly discussed. The results could have applications to the problem of computing sound sources in isotropic homogeneous turbulence by LES

Numerical simulation of rock failure and earthquake process on mesoscopic scale

On the basis of the lattice model of MORA and PLACE, Discrete Element Method, and Molecular Dynamics approach, another kind of numerical model is developed. The model consists of a 2-D set of particles linked by three kinds of interactions and arranged into triangular lattice. After the fracture criterion and rules of changes between linking states are given, the particle positions, velocities and accelerations at every time step are calculated using a finite-difference scheme, and the configuration of particles can be gained step by step. Using this model, realistic fracture simulations of brittle solid (especially under pressure) and simulation of earthquake dynamics are made.

Extended self-similar scaling law of multi-scale eddy structure in wall turbulence

The longitudinal fluctuating velocity of a turbulent boundary layer was measured in a water channel at a moderate Reynolds number. The extended self-similar scaling law of structure function proposed by Benzi was verified. The longitudinal fluctuating velocity, in the turbulent boundary layer was decomposed into many multi-scale eddy structures by wavelet transform. The extended self-similar scaling law of structure function for each scale eddy velocity was investigated. The conclusions are I) The statistical properties of turbulence could be self-similar not only at high Reynolds number, but also at moderate and low Reynolds number, and they could be characterized by the same set of scaling exponents xi (1)(n) = n/3 and xi (2)(n) = n/3 of the fully developed regime. 2) The range of scales where the extended self-similarity valid is much larger than the inertial range and extends far deep into the dissipation range,vith the same set of scaling exponents. 3) The extended selfsimilarity is applicable not only for homogeneous turbulence, but also for shear turbulence such as turbulent boundary layers.

Deformation behavior and microstructure effect in 2124Al/SiCp composite

Dynamic compression tests were performed by means of a Split Hopkinson Pressure Bar (SHPB). Test materials were 2124Al alloys reinforced with 17% volume fraction of 3, 13 and 37 μm SiC particles, respectively. Under strain rate ε = 2100 l/s, SiC particles have a strong effect on σ_0.2 of the composites and the σ_0.2 increases with different SiC particle size in the following order: 2124Al-alloy → 124Al/SiC_p (37 μm) → 2124Al/SiC_p (13 μm) → 2124Al/SiC_p (3 μm), and the strain hardening of the composites depends mainly on the strain hardening of matrix, 2124A1 alloy. The results of dimensional analysis present that the flow stress of these composites not only depends on the property of reinforcement and matrix but also relates to the microstructure scale, matrix grain size, reinforcement size, the distance between reinforcements and dislocations in matrix. The normalized flow stress here is a function of inverse power of the edge-edge particle spacing, dislocation density and matrix grain size. Close-up observation shows that, in the composite containing SiC particles (3 μm), localized deformation formed readily comparing with other materials under the same loading condition. Microscopic observations indicate that different plastic flow patterns occur within the matrix due to the presence of hard particles with different sizes.

Microdamage Evolution, Energy Dissipation and Their Trans-Scale Effects on Macroscopic Failure

The process of damage evolution concerns various scales, from micro- to macroscopic. How to characterize the trans-scale nature of the process is on the challenging frontiers of solid mechanics. In this paper, a closed trans-scale formulation of damage evolution based on statistical microdamage mechanics is presented. As a case study, the damage evolution in spallation is analyzed with the formulation. Scaling of the formulation reveals that the following dimensionless numbers: reduced Mach number M, damage number S, stress wave Fourier number P, intrinsic Deborah number D*, and the imposed Deborah number De*, govern the whole process of deformation and damage evolution. The evaluation of P and the estimation of temperature increase show that the energy equation can be ignored as the first approximation in the case of spallation. Hence, apart from the two conventional macroscopic parameters: the reduced Mach number M and damage number S, the damage evolution in spallation is mainly governed by two microdamage-relevant parameters: the Deborah numbers D* and De*. Higher nucleation and growth rates of microdamage accelerate damage evolution, and result in higher damage in the target plate. In addition, the mere variation in nucleation rate does not change the spatial distribution of damage or form localized rupture, while the increase of microdamage growth rate localizes the damage distribution in the target plate, which can be characterized by the imposed Deborah number De*.

Mesofracture mechanics: a necessary link

Classical fracture mechanics is based on the premise that small scale features could be averaged to give a larger scale property such that the assumption of material homogeneity would hold. Involvement of the material microstructure, however, necessitates different characteristic lengths for describing different geometric features. Macroscopic parameters could not be freely exchanged with those at the microscopic scale level. Such a practice could cause misinterpretation of test data. Ambiguities arising from the lack of a more precise range of limitations for the definitions of physical parameters are discussed in connection with material length scales. Physical events overlooked between the macroscopic and microscopic scale could be the link that is needed to bridge the gap. The classical models for the creation of free surface for a liquid and solid are oversimplified. They consider only the translational motion of individual atoms. Movements of groups or clusters of molecules deserve attention. Multiscale cracking behavior also requires the distinction of material damage involving at least two different scales in a single simulation. In this connection, special attention should be given to the use of asymptotic solution in contrast to the full field solution when applying fracture criteria. The former may leave out detail features that would have otherwise been included by the latter. Illustrations are provided for predicting the crack initiation sites of piezoceramics. No definite conclusions can be drawn from the atomistic simulation models such as those used in molecular dynamics until the non-equilibrium boundary conditions can be better understood. The specification of strain rates and temperatures should be synchronized as the specimen size is reduced to microns. Many of the results obtained at the atomic scale should be first identified with those at the mesoscale before they are assumed to be connected with macroscopic observations. Hopefully, "mesofracture mechanics" could serve as the link to bring macrofracture mechanics closer to microfracture mechanics.

Effects Of Subgrid-Scale Modeling On Lagrangian Statistics In Large-Eddy Simulation

The application of large-eddy simulation (LES) to turbulent transport processes requires accurate prediction of the Lagrangian statistics of flow fields. However, in most existing SGS models, no explicit consideration is given to Lagrangian statistics. In this paper, we focus on the effects of SGS modeling on Lagrangian statistics in LES ranging from statistics determining single-particle dispersion to those of pair dispersion and multiparticle dispersion. Lagrangian statistics in homogeneous isotropic turbulence are extracted from direct numerical simulation (DNS) and the LES with a spectral eddy-viscosity model. For the case of longtime single-particle dispersion, it is shown that, compared to DNS, LES overpredicts the time scale of the Lagrangian velocity correlation but underpredicts the Lagrangian velocity fluctuation. These two effects tend to cancel one another leading to an accurate prediction of the longtime turbulent dispersion coefficient. Unlike the single-particle dispersion, LES tends to underestimate significantly the rate of relative dispersion of particle pairs and multiple-particles, when initial separation distances are less than the minimum resolved scale due to the lack of subgrid fluctuations. The overprediction of LES on the time scale of the Lagrangian velocity correlation is further confirmed by a theoretical analysis using a turbulence closure theory.

Short-and resonant-range interactions between scales in turbulence and their applications

Interactions between different scales in turbulence were studied starting from the incompressible Navier-Stokes equations. The integral and differential formulae of the short-range viscous stresses, which express the short-range interactions between contiguous scales in turbulence, were given. A concept of the resonant-range interactions between extreme contiguous scales was introduced and the differential formula of the resonant-range viscous stresses was obtained. The short- and resonant-range viscous stresses were applied to deduce the large-eddy simulation (LES) equations as well as the multiscale equations, which are approximately closed and do not contain any empirical constants or relations. The properties and advantages of using the multiscale equations to compute turbulent flows were discussed. The short-range character of the interactions between the scales in turbulence means that the multiscale simulation is a very valuable technique for the calculation of turbulent flows. A few numerical examples were also given.

Small scale, grain size and substrate effects in nano-indentation experiment of film-substrate systems

The mechanical properties of film-substrate systems have been investigated through nano-indentation experiments in our former paper (Chen, S.H., Liu, L., Wang, T.C., 2005. Investigation of the mechanical properties of thin films by nano-indentation, considering the effects of thickness and different coating-substrate combinations. Surf. Coat. Technol., 191, 25-32), in which Al-Glass with three different film thicknesses are adopted and it is found that the relation between the hardness H and normalized indentation depth h/t, where t denotes the film thickness, exhibits three different regimes: (i) the hardness decreases obviously with increasing indentation depth; (ii) then, the hardness keeps an almost constant value in the range of 0.1-0.7 of the normalized indentation depth h/t; (iii) after that, the hardness increases with increasing indentation depth. In this paper, the indentation image is further investigated and finite element method is used to analyze the nano-indentation phenomena with both classical plasticity and strain gradient plasticity theories. Not only the case with an ideal sharp indenter tip but also that with a round one is considered in both theories. Finally, we find that the classical plasticity theory can not predict the experimental results, even considering the indenter tip curvature. However, the strain gradient plasticity theory can describe the experimental data very well not only at a shallow indentation depth but also at a deep depth. Strain gradient and substrate effects are proved to coexist in film-substrate nano-indentation experiments. (c) 2006 Elsevier Ltd. All rights reserved.

Multi-scale analysis of AFM tip and surface interactions

Thoroughly understanding AFM tip-surface interactions is crucial for many experimental studies and applications. It is important to realize that despite its simple appearance, the system of tip and sample surface involves multiscale interactions. In fact, the system is governed by a combination of molecular force (like the van der Waals force), its macroscopic representations (such as surface force) and gravitational force (a macroscopic force). Hence, in the system, various length scales are operative, from sub-nanoscale (at the molecular level) to the macroscopic scale. By integrating molecular forces into continuum equations, we performed a multiscale analysis and revealed the nonlocality effect between a tip and a rough solid surface and the mechanism governing liquid surface deformation and jumping. The results have several significant implications for practical applications. For instance, nonlocality may affect the measurement accuracy of surface morphology. At the critical state of liquid surface jump, the ratio of the gap between a tip and a liquid dome (delta) over the dome height (y(o)) is approximately (n-4) (for a large tip), which depends on the power law exponent n of the molecular interaction energy. These findings demonstrate that the multiscale analysis is not only useful but also necessary in the understanding of practical phenomena involving molecular forces. (c) 2007 Elsevier Ltd. All rights reserved.

Crack growth by grain boundary cavitation in the transient and extensive creep regimes

Crack growth due to cavity growth and coalescence along grain boundaries is analyzed under transient and extensive creep conditions in a compact tension specimen. Account is taken of the finite geometry changes accompanying crack tip blunting. The material is characterized as an elastic-power law creeping solid with an additional contribution to the creep rate arising from a given density of cavitating grain boundary facets. All voids are assumed present from the outset and distributed on a given density of cavitating grain boundary facets. The evolution of the stress fields with crack growth under three load histories is described in some detail for a relatively ductile material. The full-field plane strain finite element calculations show the competing effects of stress relaxation due to constrained creep, diffusion and crack tip blunting. and of stress increase due to the instantaneous elastic response to crack growth. At very high crack growth rates the Hui-Riedel fields dominate the crack tip region. However. the high growth rates are not sustained for any length of time in the compact tension geometry analyzed. The region of dominance of the Hui-Riedel field shrinks rapidly so that the near-tip fields are controlled by the HRR-type field shortly after the onset of crack growth. Crack growth rates under various conditions of loading and spanning the range of times from small scale creep to extensive creep are obtained. We show that there is a strong similarity between crack growth history and the behaviour of the C(t) and C(t) parameters. so that crack growth rates correlate rather well with C(t) and C(t). A relatively brittle material is also considered that has a very different near-tip stress field and crack growth history.

A Simplified Evaluation Procedure of the Scale of Fluctuation Combined Recurrence Averaging and Weighted Linear Regression

The divergence of properties from one location to another within a soil mass is termed spatial variability, which traditionally includes three parameters the mean, the standard deviation, and the scale of fluctuation, in order to stochastically describe a soil property. Among them, determining the scale of fluctuation in the evaluation of spatial variability of soil profiles is not easy due to soil condition complexity. A simplified procedure is presented in the paper to determine the scale of fluctuation combined recurrence averaging and weighted linear regression. The alternative approach utilizes widely usable spreadsheet to solve the problem more directly and efficiently.

The Relationship between Decision Trees and the Scale of Train Data Set

Decision Trees need train samples in the train data set to get classification rules. If the number of train data was too small, the important information might be missed and thus the model could not explain the classification rules of data. While it is not affirmative that large scale of train data set can get well model. This Paper analysis the relationship between decision trees and the train data scale. We use nine decision tree algorithms to experiment the accuracy, complexity and robustness of decision tree algorithms. Some results are demonstrated.

Monte Carlo Modeling of Micro Scale Gas Flows

An information preservation (IP) method has been used to simulate many micro scale gas flows. It may efficiently reduce the statistical scatter inherent in conventional particle approaches such as the direct simulation Monte Carlo (DSMC) method. This paper reviews applications of IP to some benchmark problems. Comparison of the IP results with those given by experiment, DSMC, and the linearized Boltzmann equation, as well as the Navier-Stokes equations with a slip boundary condition, and the lattice Boltzmann equation, shows that the IP method is applicable to micro scale gas flows over the entire flow regime from continuum to free molecular.