Derivation Of A Nonlinear Reynolds Stress Model Using Renormalization Group Analysis And Two-Scale Expansion Technique

Adopting Yoshizawa's two-scale expansion technique, the fluctuating field is expanded around the isotropic field. The renormalization group method is applied for calculating the covariance of the fluctuating field at the lower order expansion. A nonlinear Reynolds stress model is derived and the turbulent constants inside are evaluated analytically. Compared with the two-scale direct interaction approximation analysis for turbulent shear flows proposed by Yoshizawa, the calculation is much more simple. The analytical model presented here is close to the Speziale model, which is widely applied in the numerical simulations for the complex turbulent flows.

Sand storms: CFD analysis of Reynolds stress and collision stress of particles near sand bed

Sand storm is a serious environmental threat to humans. Sand particles are transported by saltation and suspension, causing soil erosion in one place and deposition in another. In order to prevent and predict sand storms, the causes and the manners of particle motions must be studied in detail. In this paper a standard k-epsilon model is used for the gas phase simulation and the discrete element method (DEM) is used to predict the movements of particles using an in-house procedure. The data are summarized in an Eulerian-Eulerian regime after simulation to get the statistical particle Reynolds stress and particle collision stress. The results show that for the current case the Reynolds stress and the air shear stress predominate in the region 20-250 mm above the initial sand bed surface. However, in the region below 3 mm, the collision stress must be taken into account in predicting particle movement. (C) 2010 Chinese Society of Particuology and Institute of Process Engineering, Chinese Academy of Sciences. Published by Elsevier B.V. All rights reserved.

Statistics of Turbulence in Benthal Boundary Layers

This paper deals with turbulence behavior inbenthalboundarylayers by means of large eddy simulation (LES). The flow is modeled by moving an infinite plate in an otherwise quiescent water with an oscillatory and a steady velocity components. The oscillatory one aims to simulate wave effect on the flow. A number of large-scale turbulence databases have been established, based on which we have obtained turbulencestatisticsof the boundarylayers, such as Reynolds stress, turbulence intensity, skewness and flatness ofturbulence, and temporal and spatial scales of turbulent bursts, etc. Particular attention is paid to the dependences of those statistics on two nondimensional parameters, namely the Reynolds number and the current-wave velocity ratio defined as the steady current velocity over the oscillatory velocity amplitude. It is found that the Reynolds stress and turbulence intensity profile differently from phase to phase, and exhibit two types of distributions in an oscillatory cycle. One is monotonic occurring during the time when current and wave-induced components are in the same direction, and the other inflectional occurring during the time when current and wave-induced components are in opposite directions. Current component makes an asymmetrical time series of Reynolds stress, as well as turbulence intensity, although the mean velocity series is symmetrical as a sine/cosine function. The skewness and flatness variations suggest that the turbulence distribution is not a normal function but approaches to a normal one with the increasing of Reynolds number and the current-wave velocity ratio as well. As for turbulent bursting, the dimensionless period and the mean area of all bursts per unit bed area tend to increase with Reynolds number and current-wave velocity ratio, rather than being constant as in steady channel flows.

NUMERICAL AND EXPERIMENTAL STUDY ON LIQUID-SOLID FLOW IN A HYDROCYCLONE

Hydrocyclones are widely used in industry, of which the geometrical design using CFD techniques is gaining more popularity in recent years. In this study, the Euler-Euler approach and the Reynolds stress model are applied to simulate the liquid-solid flowfield in a hydrocyclone. The methodology is validated by a good agreement between experimental data and numerical results. Within the research range, the simulation indicates that the liquid-solid separation mainly occurs in the conical segment, and increasing conical height or decreasing cylindrical height helps to improve the grade efficiencies of solid particles. Based on these results, two of the same hydrocyclones are designed and installed in series to establish a liquid-solid separation system. Many experiments are then conducted under different conditions, in which the effects of the water cut and the second hydrocyclone on the separation are investigated. The results also confirm that smaller solid particles are more susceptible to the inlet conditions, and the second hydrocyclone plays a more important role as the water cut reduces.

The mechanism of energy loss and transition in a flow with submerged vegetation

The mechanism of energy balance in an open-channel flow with submerged vegetation was investigated. The energy borrowed from the local flow, energy spending caused by vegetation drag and flow resistance, and energy transition along the water depth were calculated on the basis of the computational results of velocity and Reynolds stress. Further analysis showed that the energy spending in a cross-section was a maximum around the top of the vegetation, and its value decreased progressively until reaching zero at the flume bed or water surface. The energy borrowed from the local flow in the vegetated region could not provide for spending; therefore, surplus borrowed energy in the non-vegetated region was transmitted to the vegetated region. In addition, the total energy transition in the cross-section was zero; therefore, the total energy borrowed from the flow balanced the energy loss in the whole cross-section. At the same time, we found that there were three effects of vegetation on the flow: turbulence restriction due to vegetation, turbulence source due to vegetation and energy transference due to vegetation, where the second effect was the strongest one. Crown Copyright (C) 2010 Published by Elsevier Ltd. All rights reserved.

A second-order dynamic subgrid-scale stress model

A second-order dynamic model based on the general relation between the subgrid-scale stress and the velocity gradient tensors was proposed. A priori test of the second-order model was made using moderate resolution direct numerical simulation date at high Reynolds number ( Taylor microscale Reynolds number R-lambda = 102 similar to 216) for homogeneous, isotropic forced flow, decaying flow, and homogeneous rotating flow. Numerical testing shows that the second-order dynamic model significantly improves the correlation coefficient when compared to the first-order dynamic models.

Stress-strain fields and the effectiveness shear properties for three-phase composites with imperfect interface

Based on the 'average stress in the matrix' concept of Mori and Tanaka (:Mori, T., Tanaka, K., 1973. Average stress in matrix and average elastic energy of materials with misfitting inclusion. Acta Metall. 21, 571-580) a micromechanical model is presented for the prediction of the elastic fields in coated inclusion composites with imperfect interfaces. The solutions of the effective elastic moduli for this kind of composite are also obtained. In two kinds of composites with coated particulates and fibers, respectively, the interface imperfections are takes to the assumption that the interface displacement discontinues are linearly related to interface tractions like a spring layer of vanishing thickness. The resulting effective shear modulus for each material and the stress fields in the composite are presented under a transverse shear loading situation.

Flow Stress of Biomorphous Metal-Matrix Composites

For metal-matrix composites (MMCs), interfacial debonding between the ductile matrix and the reinforcing hard inclusions is an important failure mode. A fundamental approach to improving the properties of MMCs is to optimize their microstructure to achieve maximum strength and toughness. Here, we investigate the flow stress of a MMC with a nanoscale microstructure similar to that of bone. Such a 'biomorphous' MMC would be made of staggered hard and slender nanoparticles embedded in a ductile matrix. We show that the large aspect ratio and the nanometer size of inclusions in the biomorphous MMC lead to significantly improved properties with increased tolerance of interfacial damage. In this case, the partially debonded inclusions continue to carry mechanical load transferred via longitudinal shearing of the matrix material between neighboring inclusions. The larger the inclusion aspect ratio, the larger is the flow stress and work hardening rate for the composite. Increasing the volume concentration of inclusion also makes the biomorphous MMC more tolerant of interfacial damage.

Control of vortex shedding at moderate Reynolds numbers

The suppression method of vortex shedding from a circular cylinder has been studied experimentally in the Reynolds number range from 300 to 1600. The test is performed in a water channel. The model cylinder is 1 cm in diameter and 38 cm in length. A row of small rods of 0.18 cm in diameter and 1.5 cm in length are perpendicularly connected to the surface of the model cylinder and distributed along the meridian, The distance between the neighboring rods and the angle of attack of the rods can be changed so that the suppression effect on vortex shedding can be adjusted. The results show that vortex shedding can be suppressed effectively if the distance between the neighboring rods is smaller than 3 times and the cylinder diameter and the angle of attack is in the range of 30degreesless than or equal tobeta<90&DEG;.

Burst Detection in Turbulent Channel Flows Based on Large Eddy Simulation Databases

Reliable turbulent channel flow databases at several Reynolds numbers have been established by large eddy simulation (LES), with two of them validated by comparing with typical direct numerical simulation (DNS) results. Furthermore, the statistics, such as velocity profile, turbulent intensities and shear stress, were obtained as well as the temporal and spatial structure of turbulent bursts. Based on the LES databases available, the conditional sampling methods are used to detect the structures of burst events. A method to deterimine the grouping parameter from the probability distribution function (pdf) curve of the time separation between ejection events is proposed to avoid the errors in detected results. And thus, the dependence of average burst period on thresholds is considerably weakened. Meanwhile, the average burst-to-bed area ratios are detected. It is found that the Reynolds number exhibits little effect on the burst period and burst-to-bed area ratio.

Increasing Critical Sensitivity of the Load/Unload Response Ratio before Large Earthquakes with Identified Stress Acclumulation Pattern

The Load/Unload Response Ratio (LURR) method is proposed for short-to-intermediate-term earthquake prediction [Yin, X.C., Chen, X.Z., Song, Z.P., Yin, C., 1995. A New Approach to Earthquake Prediction — The Load/Unload Response Ratio (LURR) Theory, Pure Appl. Geophys., 145, 701–715]. This method is based on measuring the ratio between Benioff strains released during the time periods of loading and unloading, corresponding to the Coulomb Failure Stress change induced by Earth tides on optimally oriented faults. According to the method, the LURR time series usually climb to an anomalously high peak prior to occurrence of a large earthquake. Previous studies have indicated that the size of critical seismogenic region selected for LURR measurements has great influence on the evaluation of LURR. In this study, we replace the circular region usually adopted in LURR practice with an area within which the tectonic stress change would mostly affect the Coulomb stress on a potential seismogenic fault of a future event. The Coulomb stress change before a hypothetical earthquake is calculated based on a simple back-slip dislocation model of the event. This new algorithm, by combining the LURR method with our choice of identified area with increased Coulomb stress, is devised to improve the sensitivity of LURR to measure criticality of stress accumulation before a large earthquake. Retrospective tests of this algorithm on four large earthquakes occurred in California over the last two decades show remarkable enhancement of the LURR precursory anomalies. For some strong events of lesser magnitudes occurred in the same neighborhoods and during the same time periods, significant anomalies are found if circular areas are used, and are not found if increased Coulomb stress areas are used for LURR data selection. The unique feature of this algorithm may provide stronger constraints on forecasts of the size and location of future large events.

Can stress-strain relationships be obtained from indentation curves using conical and pyramidal indenters?

Applying the scaling relationships developed recently for conical indentation in elastic-plastic solids with work-hardening, we examine the question of whether stress-strain relationships of such solids can be uniquely determined by matching the calculated loading and unloading curves with that measured experimentally. We show that there can be multiple stress-strain curves for a given set of loading and unloading curves. Consequently, stress-strain relationships may not be uniquely determined from loading and unloading curves alone using a conical or pyramidal indenter.

Strain gradient theory with couple stress for crystalline solids

A new phenomenological strain gradient theory for crystalline solid is proposed. It fits within the framework of general couple stress theory and involves a single material length scale Ics. In the present theory three rotational degrees of freedom omega (i) are introduced, which denote part of the material angular displacement theta (i) and are induced accompanying the plastic deformation. omega (i) has no direct dependence upon u(i) while theta = (1 /2) curl u. The strain energy density omega is assumed to consist of two parts: one is a function of the strain tensor epsilon (ij) and the curvature tensor chi (ij), where chi (ij) = omega (i,j); the other is a function of the relative rotation tensor alpha (ij). alpha (ij) = e(ijk) (omega (k) - theta (k)) plays the role of elastic rotation reason The anti-symmetric part of Cauchy stress tau (ij) is only the function of alpha (ij) and alpha (ij) has no effect on the symmetric part of Cauchy stress sigma (ij) and the couple stress m(ij). A minimum potential principle is developed for the strain gradient deformation theory. In the limit of vanishing l(cs), it reduces to the conventional counterparts: J(2) deformation theory. Equilibrium equations, constitutive relations and boundary conditions are given in detail. For simplicity, the elastic relation between the anti-symmetric part of Cauchy stress tau (ij), and alpha (ij) is established and only one elastic constant exists between the two tensors. Combining the same hardening law as that used in previously by other groups, the present theory is used to investigate two typical examples, i.e., thin metallic wire torsion and ultra-thin metallic beam bend, the analytical results agree well with the experiment results. While considering the, stretching gradient, a new hardening law is presented and used to analyze the two typical problems. The flow theory version of the present theory is also given.

Investigation of Dynamic Fracture Toughness under Compressive-Shear Combinede Stress Wave Loading

提出亚微秒单脉冲应力波载荷作用下II型裂纹的平板冲击实验技术。加载率为dK/dt-10~8MPa·m~{”/d}·s~{-1}。实验中由锰铜应力片和弹性波理论分别测定和计算了压应力；通过微观分析确定了动态裂纹的平均扩展长度；引进等效应力强度因子,用动态断裂理论确定了60号钢的动态断裂韧性K_{Id}和K_{IId}；建立了亚微秒冲击载荷作用下确定材料动态断裂韧性的方法。

Numerical Simulation of Transient Thermal Field and Residual Stress in Laser Melting Process

以激光熔凝表面强韧化处理为背景,应用空间弹塑性有限单元和高精度数值算法同时考虑材料组织性能的变化模拟工件的温度场及残余应力,研究激光熔凝加工中瞬时温度场及残余应力数值模拟,同时考虑相变潜热及相变塑性的影响,用算例验证了模型的正确性,给出了不同时刻温度场分布及残余应力分布。