The influence of normal stress anisotropy in predicting scalar dispersion with the v–f model

A numerical study of scalar dispersion is presented to investigate the effectiveness of pairing the v²–f turbulence model with algebraic models for the scalar flux. This approach is contrasted with utilizing a full Second Moment Closure (SMC) as the flow field input to the scalar model. Predictions of scalar transport in a turbulent channel and over a wavy wall are compared to available DNS databases. The latter case includes a scalar release from a point source and therefore detailed comparisons of the three-component turbulent scalar flux are reported. It is found that the transported variable v², representing the near wall turbulent velocity fluctuation scale, can be used to increase the level of normal stress anisotropy provided to algebraic scalar models and thereby improve mean scalar prediction over that of the Standard Gradient Diffusion Hypothesis (SGDH). Improvement is most significant in the near wall region. Three specifications of the normal stresses, derived from v², are considered to provide the link from the v²–f model to the algebraic flux models used to close the scalar transport equation. Barycentric maps are used to examine the state of turbulence anisotropy in each case. As the anisotropy in the normal stress specification becomes more accurate, improvements are realized in the prediction of the spanwise flux as well as the mean concentration.     

Four-component gas/water/oil displacements in one dimension: Part I. structure of the conservation law

This paper presents the first analysis of the mathematical structure of a system of conservation laws modeling compositional flow of four components in three phases. The phase behavior that results from assuming the equilibrium volume ratios of the components in the phases are fixed (constant K-values) when up to three phases may form, is studied. We parameterize the equations in the three-phase region and show that within the three-phase region two of the characteristic curves can be found using three-phase immiscible flow theory. The third eigenvalue can also be found analytically when the K-values are constant. We show that the eigenvalue problem given by the conservation law has a discontinuity at the boundary of the two- and three-phase regions. Finally, the loss of strict hyperbolicity is discussed.     

基于改进位移模式的一维C1有限元超收敛算法

提出了基于改进位移模式的一维C1有限元超收敛算法。利用单元内部需满足平衡方程的条件,推导了超收敛计算解析公式的显式,即将高阶有限元解的位移模式用常规有限元解的位移模式表示。用常规有限元解的位移模式与高阶有限元解的位移模式之和构造新的位移模式。采用积分形式推导了单元刚度矩阵。该算法在前处理阶段使用了超收敛计算公式,在常规试函数的基础上,增加了高阶试函数,使得单元内平衡方程的残差减少,从而达到提高精度的目标。对于Hermite单元,本文的结点和单元的位移、导数都达到了h4阶的超收敛精度。     

THE ANALYSES OF THE THREE-DIMENSIONAL STRESS STRUCTURE NEAR THE CRACK TIP OF MODE I CT SPECIMENS IN ELASTICPLASTIC STATE(Ⅰ)—THE ANALYSES OF CONSTRAINT PARAMETERS AND FRACTURE PARAMETERS

In the present paper, three dimensional analyses of some general constraint parameters and fracture parameters near the crack tip. of Mode I CT specimens in two different thicknesses are carried out by employing ADINA program. The results reveal that the constraints along the thickness direction are obviously separated into two parts: the keeping similar high constraint field (Z₁) and rapid reducing constraints one (Z₂). The two fields are experimentally confiremed to correspond to the smooth region and the shear lip on the fracture face respectively. So the three dimensional stress structure of Mode I specimens can be derived through discussing the two fields respectively. The distribution of the Crack Tip Opening Displacement (CTOD) along the thickness direction and the three dimensional distribution of the void growth ratio (V_g) near the crack tip are also obtained. The two fracture parameters are in similar trends along the thickness direction, and both of them can reflect the effect of thickness and that of the loading level to a certain degree.     

The flow past a sphere in a cylindrical tube: effects of intertia,shear-thinning and elasticity

Numerical solutions are presented for the flow past a sphere placed at the centreline of a cylindrical tube for Reynolds numbers ranging from 0 to 150, using a boundary element method. Fluids are modelled by a variety of constitutive equations including the Newtonian, the Carreau and the Phan-Thien-Tanner models. The influence of inertia, shear-thinning and fluid elasticity on the flow field, drag and the pressure drop force-drag ratio is examined. Some results are compared with available experimental data.     

Hybrid RANS/LES of flow and heat transfer in round impinging jets

Fluid flow and convective heat transfer predictions are presented of round impinging jets for several combinations of nozzle-plate distances H/D = 2, 6 and 13.5 (where D is the nozzle diameter) and Reynolds numbers Re = 5000, 23,000 and 70,000 with the newest version of the k-ω model of Wilcox (2008) and three hybrid RANS/LES models. In the RANS mode of the hybrid RANS/LES models, the k-ω model is recovered. Three formulations are considered to activate the LES mode. The first model is similar to the hybrid models of Davidson and Peng (2003) and Kok et al. (2004). The turbulent length scale is replaced by the grid size in the destruction term of the k-equation and in the definition of the RANS eddy viscosity. As grid size, a maximum measure of the hexahedral grid cell is used. The second model has the same k-equation, but the eddy viscosity is the minimum of the k-ω eddy viscosity and the Smagorinsky eddy viscosity, following a proposal by Batten et al. (2004). The Smagorinsky eddy viscosity is formed with the cube root of the cell volume. The third model has, again, the same k-equation, but has an eddy viscosity which is an intermediate between the eddy viscosities of the first and second models. This is reached by using the cube root of the cell volume in the eddy viscosity formula of the first model.The simulation results are compared with experimental data for the high Reynolds number cases Re = 23,000 and Re = 70,000 and LES data for the low-Reynolds number case Re = 5000. The Reynolds numbers are defined with the nozzle diameter and the bulk velocity at nozzle outlet. At low nozzle-plate distance (the impingement plate is in the core of the jet), turbulent kinetic energy is overpredicted by RANS in the stagnation flow region. This leads to overprediction of the heat transfer rate along the impingement plate in the impact zone. At high nozzle-plate distance (the impingement plate is in the mixed-out region of the jet), the turbulence mixing is underpredicted by RANS in the shear layer of the jet which gives a too high length of the jet core. This also results in overprediction of the heat transfer rate in the impingement zone caused by too big temperature gradients at impingement.All hybrid RANS/LES models are able to correct the heat transfer overprediction of the RANS model. For good predictions at low nozzle-plate distance, it is necessary to sufficiently resolve the formation and development of the near-wall vortices in the jet impingement region. At high nozzle-plate distance, the essence is to capture the evolution and breakup of the flow unsteadiness in the shear layer of the jet, so that accurate mean and fluctuating velocity profiles are obtained in the impingement region. Although the models have a quite different theoretical justification and generate a quite different eddy viscosity in some flow regions, their overall results are very comparable. The reason is that in zones that are crucial for the results, the models behave similarly.     

Degradation of nylon-6/clay nanocomposites in NO<Subscript><Emphasis Type="Italic">x</Emphasis></Subscript>

Nylon-6 is an important engineering polymer that, in its fully spherulitic (bulk) form, has many applications in gears, rollers, and other long life cycle components. In 1993, Toyota commercialized a nylon-6/clay nanocomposite out of which it produced the timing belt cover for the Camry. Although these hybrid nanocomposites show significant improvements in their mechanical response characteristics, including yield strength and heat distortion temperature, little is known about the degradation of these properties due to environmental pollutants like NO _x . Nylon-6 fibers are severely degraded by interaction with NO _x and other pollutants, showing a strong synergy between applied load and environmental degradation. While the nanocomposites show a significant reduction in permeability of gases and water due to the incorporation of lamellar clay, their susceptibility to nondiffusional mechano-chemical degradation is unknown. The fracture toughness of these nylon-6/clay nanocomposites increases, not as a function of clay content, but as a function of the volume of nylon-6 polymer chains influenced by the clay lamellae surfaces. Both the clay and the constrained volume offer the nanocomposites some protection from the deleterious effects of NO _x . The time-to-failure at a given stress intensity factor as a function of clay content and constrained volume is discussed along with fracture toughness of the materials.     

非惯性参考系中弹性薄板在参数激励与强迫激励联合作用下的1／2亚谐共振分岔分析

本文研究了非惯性参考系中弹性薄板在范围运动与变形运动相互耦合时的1／2亚谐共振分岔,在建立了该系统的动力学控制方程的基础上,利用多尺度法得到了参数激励与强迫激励联合作用下非惯性参考系中弹性薄板1／2亚谐共振时的分岔响应方程及其分岔集,讨论了该动力系统的稳定性,给出了它的五种分响应曲线。     

A study of the in situ transverse tensile lamina strength and the classical solution of the opening mode crack in a 0/90/0 laminate

In this paper, the classical solution of the opening mode crack in the 90° layer of 0/90/0 laminates has been determined by means of Fourier transformations and the procedure of Copson for a pair of dual integral equations. The fracture behavior and the in situ transverse strength of the 90° layer have been quantitatively studied in graphite/epoxy laminates, based on the solution obtained above. The results show that the stress intensity factor of this kind of laminates, which is different from that of a single unidirectional 90° layer, decreases with the increase in thicknessb, or modulusE _L orG _LT of the 0° layer and also decreases with the decrease in the thickness of the 90° layer. So the lamination effect manifests itself and thein situ transverse strength of the 90° layer is thereby enchanced. The theoretical calculations agree with the experimental data presented by D.L. Flaggs. Supported by National Natural Science Fundation of China.     

THE NEW CRITERIA OF ELASTIC AND FATIGUE FAILURE IN THE COMPONENT OF COMPLEX STRESS STATES

In this paper, a total criterion on elastic and fatigue failure in complex stress, that is, octahedral stress strength theory on dynamic and static states on the basis of studying modern and classic strength theories. At the same time, an analysis of an independent and fairly comprehensive theoretical system is set up. It gives generalized failure factor by 36 materials and computative theory of the 11 states of complex stresses on a point, and derives the operator equation on generalized allowable strength and a computative method for a total equation can be applied to dynamic and static states. It is illustrated that the method has a good exactness through computation of eight examples of engineering. Therefore, the author suggests applying it to engineering widely.     

The free energy of mixing for n-variant martensitic phase transformations using quasi-convex analysis

The construction of effective models for materials that undergo martensitic phase transformations requires usable and accurate functional representations for the free energy density. The general representation of this energy is known to be highly non-convex; it even lacks the property of quasi-convexity. A quasi-convex relaxation, however, does permit one to make certain estimates and powerful conclusions regarding phase transformation. The general expression for the relaxed free energy is however not known in the n-variant case. Analytic solutions are known only for up to 3 variants, whereas cases of practical interests involve 7-13 variants. In this study we examine the n-variant case utilizing relaxation theory and produce a seemingly obvious but very powerful observation regarding a lower bound to the quasi-convex relaxation that makes practical evolutionary computations possible. We also examine in detail the 4-variant case where we explicitly show the relation between three different forms of the free energy of mixing: upper bound by lamination, the Reuß lower bound, and a lower estimate of the -measure bound. A discussion of the bounds and their utility is provided; sample computations are presented for illustrative purposes.