DEM simulations of tote blenders for enhanced axial mixing efficiency

The mixing performance of a multi-bladed tote blender is investigated using a graphics processing unit-based discrete element method program.The positioning,dimensions,and applicability of the baffles are systematically studied according to the axial mixing efficiency.The results indicate that the novel inclined multi-bladed baffles can break the symmetrical axial granular flow and introduce a more efficient convective flow into the granular mixing in the axial direction of the tote blender.Owing to the joint effects of convective mixing and asymmetrical granular flow,the axial mixing efficiency is increased by a factor of nearly 20.More importantly,the novel baffle placement approach exhibits excellent applicability to different operating conditions,particle shapes,and blender sizes.Additionally,the inclined baffles can prevent the segregation caused by shape discrepancies and improve the mixing homogeneity in the steady state.The novel baffle design is promising for applications in more complex industrial blenders for achieving a high axial mixing efficiency.     

Behavior of near-field dilution of thermal buoyant jet discharged horizontally in compound open-channel

The RNG κ-ε model considering the buoyancy effect, which is solved by the hybrid finite analytic method, is used to simulate the mixture of the horizontal round thermal buoyant jet in compound open channel flow. The mixing features near the spout and flowing characteristic of the secondary currents are studied by numerical simulation. Meanwhile, （1） the distribution of the measured isovels for stream-wise velocity, （2） secondary currents, （3） the distribution of the measured isovels for temperature of typical cross-section near the spout, were obtained by the three-dimensional Micro ADV and the temperature measuring device. Compared with experimental data, the RNG κ-ε model based on buoyancy effect can preferably simulate the jet which performs the bifurcation phenomenon, jet reattachment （Conada effect） and beach secondary currents phenomenon with the effect of ambient flow, buoyancy, and secondary currents of compound section and so on.     

Numerical investigations of asymmetric flow of a micropolar fluid between two porous disks

Numerical solution is presented for the two- dimensional flow of a micropolar fluid between two porous coaxial disks of different permeability for a range of Reynolds number Re （-300≤ Re 〈 0） and permeability parameter A （1.0≤A ≤2.0）. The main flow is superimposed by the injection at the surfaces of the two disks. Von Karman＇s similarity transformations are used to reduce the governing equations of motion to a set of non-linear coupled ordinary differential equations （ODEs） in dimensionless form. An algorithm based on the finite difference method is employed to solve these ODEs and Richardson＇s extrapolation is used to obtain higher order accuracy. The results indicate that the parameters Re and A have a strong influence on the velocity and microrotation profiles, shear stresses at the disks and the position of the viscous/shear layer. The micropolar material constants cl, c2, c3 have profound effect on microrotation as compared to their effect on streamwise and axial velocity profiles. The results of micropolar fluids are compared with the results for Newtonian fluids.     

BEHAVIOR OF OBSTRUCTED SQUARE BUOYANT VERTICAL JETS IN STATIC AMBIENT （Ⅱ）——ANALYSIS ON BEHAVIOR OF FLOW FIELD

Based on a series of numerical calculations, the behavior of flow field in obstructed square buoyant vertical jet is summarized and analyzed. Based on the axial line velocity distribution, the flow after the disc can be divided into three regions, i.e., recirculation region, transitional region and self-similar region The characteristic of self-similarity of upright velocity was validated. The three regions can also be distinguished based on the axial velocity. The axial velocity in self-similar region was found to obey the same law and the formula presented by introducing the velocity expression used by Chen and Rodi. The isolines of pressure on cross-sections of different heights were displayed and the production, expansion, breaking and disappearing of negative pressure regions were found.     

INVESTIGATION OF THE STABILITY FOR INVISCID COMPRESSIBLE SWIRLING FLOW BETWEEN CONCENTRIC CYLINDERS

The temporal stability on inviscid compressible swirling flow between two concntric cylinders is investigated. First. a linearized differential equation is derived. Two stbility criteria are derired for compressible swirling flow by an analytic method analogous to Ludwieg's method. A finite-difference numerical method is then used to solve the eigenvalue problem of this differential equation, to get temporal growth rate and to check these stability criteria derived. Finally, The effect of compressbility for stability is disscused.     

Numerical simulation of motion and deformation of ring bubble along body surface

Numerical simulation for fluid flow over an attached rigid body with a deformable ring bubble is analyzed based on the velocity potential theory together with the boundary element method （BEM）. The analysis is focused on the axisymmetric case. The bubble surface is treated as a well defined air-liquid interface and is tracked by a mixed Eulerian-Lagrangian method. The points of intersection between the bubble and body are treated, specially in the numerical procedure. The auxiliary function method is adopted to calculate the pressure on the body surface and in the flow field. The convergence study is undertaken to assess the developed numerical method and the computation code. Some case studies are undertaken in which the interactions between the bubble/body and the incoming flow field are simulated. The effects of various physical parameters on the interactions are investigated.     

Shear flow induced vibrations of long slender cylinders with a wake oscillator model

A time domain model is presented to study the vibrations of long slender cylinders placed in shear flow. Long slender cylinders such as risers and tension legs are widely used in the field of ocean engineering. They are subjected to vortex-induced vibrations(VIV) when placed within a transverse incident flow. A three dimensional model coupled with wake oscillators is formulated to describe the response of the slender cylinder in cross-flow and in-line directions. The wake oscillators are distributed along the cylinder and the vortex-shedding frequency is derived from the local current velocity. A non-linear fiuid force model is accounted for the coupled effect between cross-flow and in-line vibrations. The comparisons with the published experimental data show that the dynamic features of VIV of long slender cylinder placed in shear flow can be obtained by the proposed model,such as the spanwise average displacement,vibration frequency,dominant mode and the combination of standing and traveling waves. The simulation in a uniform flow is also conducted and the result is compared with the case of nonuniform flow. It is concluded that the flow shear characteristic has significantly changed the cylinder vibration behavior.     

Periodical streaming potential and electro-viscous effects in microchannel flow

This paper presents an analytical solution to periodical streaming potential, flow-induced electric field and velocity of periodical pressure-driven flows in twodimensional uniform microchannel based on the Poisson-Boltzmann equations for electric double layer and Navier-Stokes equation for liquid flow. Dimensional analysis indicates that electric-viscous force depends on three factors： （1） Electric-viscous number representing a ratio between maximum of electric-viscous force and pressure gradient in a steady state, （2） profile function describing the distribution profile of electro-viscous force in channel section, and （3） coupling coefficient reflecting behavior of arnplitude damping and phase offset of electro-viscous force. Analytical results indicate that flow-induced electric field and flow velocity depend on frequency Reynolds number （Re = wh^2/v）. Flow-induced electric field varies very slowly with Re when Re 〈 1, and rapidly decreases when Re 〉 1. Electro-viscous effect on flow-induced electric field and flow velocity are very significant when the rate of the channel width to the thickness of electric double layer is small.     

Magnetohydrodynamics instability of interfacial waves between two immiscible incompressible cylindrical fluids

The problem of nonlinear instability of interfacial waves between two immiscible conducting cylindrical fluids of a weak Oldroyd 3-constant kind is studied. The system is assumed to be influenced by an axial magnetic field, where the effect of surface tension is taken into account. The analysis, based on the method of multiple scale in both space and time, includes the linear as well as the nonlinear effects. This scheme leads to imposing of two levels of the solvability conditions, which are used to construct like-nonlinear Schr6dinger equations （1-NLS） with complex coefficients. These equations generally describe the competition between nonlinearity and dispersion. The stability criteria are theoret- ically discussed and thereby stability diagrams are obtained for different sets of physical parameters. Proceeding to the nonlinear step of the problem, the results show the appearance of dual role of some physical parameters. Moreover, these effects depend on the wave kind, short or long, except for the ordinary viscosity parameter. The effect of the field on the system stability depends on the existence of viscosity and differs in the linear case of the problem from the nonlinear one. There is an obvious difference between the effect of the three Oldroyd constants on the system stability. New instability regions in the parameter space, which appear due to nonlinear effects, are shown.     

Numerical study on the performance of nozzle flow for supersonic chemical oxygen–iodine lasers

Laser performance is greatly dependent on its operating conditions due to the strong coupling among multi- physics such as gas-dynamics, chemical reaction kinetics and optics in the mixing nozzle of COIL. In this paper, 3D CFD technology is used to simulate the mixing and reactive flow of subsonic cross jet scheme at different conditions. Results obtained show that the jet penetration depth plays a dominant role in the spatial distribution of small signal gains. In the case of over-penetration, unsteady flow structures are induced by impinging between the opposing jets. The optimum spatial distribution of the chemical performance cannot be obtained even if the full penetration condition is achieved through the subsonic transverse jet mixing scheme in the COIL nozzle flow.     

一种基于直接计算高阶奇异积分的断裂力学双边界积分方程分析法

相对于有限元法,边界单元法在求解断裂问题上有着独特的优势,现有的边界单元法中主要有子区域法和双边界积分方程法.采用一种改进的双边界积分方程法求解二维、三维断裂问题的应力强度因子,对非裂纹边界采用传统的位移边界积分方程,只需对裂纹面中的一面采用面力边界积分方程,并以裂纹间断位移为未知量直接用于计算应力强度因子.采用一种高阶奇异积分的直接法计算面力边界积分方程中的超强奇异积分;对于裂纹尖端单元,提供了三种不同形式的间断位移插值函数,采用两点公式计算应力强度因子.给出了多个具体的算例,与现存的精确解或参考解对比,可得到高精度的计算结果.      

汶川地震诱发文家沟巨型滑坡 碎屑流基本特征及成因机制初步分析


2008年5月12日,汶川M80地震在四川省绵竹市清平乡文家沟内诱发一巨型滑坡。通过现场调查得知,滑坡前后缘高差455m,厚度20～30m,滑面为基岩层面,初始方量2750×107m3。滑体在运动中转化为碎屑流。滑坡-碎屑流总的水平运动距离为4022m,垂直运动距离为1443m,遗留的堆积物体积达5×107m3。滑坡距映秀—北川断裂仅36km,位于其下盘,地震烈度达XI度。滑坡导致文家沟中48人遇害,并形成一条完整的地震次生地质灾害链。初步分析表明滑坡启动速度快,滑坡向碎屑流转化过程明显、地点明确。碎屑流运动过程复杂,伴有强烈的“气垫效应”和“前缘气浪冲击效应”。作者认为,文家沟滑坡的高启动速度是长持时强烈地震动作用的结果,与山体的猛烈碰撞是导致滑体解体并转化为碎屑流的原因。     

Preface: Orbits around asteroid

One of the core issues in modern celestial mechanics is the orbital dynamics in the near-regime gravitational field of as- teroids, which provides deep insights into the mathematical nature of a class of nonlinear systems, and plays as a critical basis for in situ explorations of different science goals. Lots of efforts have been made to reveal the characteristics of orbital motion in the vicinity of asteroids, and to improve the skills of asteroid research in methodology.     

Preface: Current research progress on mechanics of high speed rail

The construction and operation of high-speed rail（HSR）grid within the past two decades in China,in terms of the scale,may be possibly comparable to any national-wide construction in the history of China,even the Great Wall.By counting railways with commercial train service at the speed of200 km/h,China has the world’s longest HSR network with over 19 369.8 km of track in service today and this number     

ACTA MECHANICA SINICA 2014 ANNUAL SUBJECT INDEX

正~~     

Effet de l'anisotropie élastique cristalline sur la distribution des facteurs de Schmid à la surface des polycristauxInfluence of crystalline elasticity anisotropy on Schmid factor distribution at the free surface of polycrystals

Experimental studies of the plasticity mechanisms of polycrystals are usually based on the Schmid factor distribution supposing crystalline elasticity isotropy. A numerical evaluation of the effect of crystalline elasticity anisotropy on the apparent Schmid factor distribution at the free surface of polycrystals is presented. Cubic elasticity is considered. Order II stresses (averaged on all grains with the same crystallographic orientation) as well as variations between averages computed on grains with the same crystallographic orientation but with different neighbour grains are computed. The Finite Element Method is used. Commonly studied metals presenting an increasing anisotropy degree are considered (aluminium, nickel, austenite, copper). Concerning order II stresses in strongly anisotropic metals, the apparent Schmid factor distribution is drifted towards small Schmid factor values (the maximum Schmid factor is equal to 0.43 instead of 0.5) and the slip activation order between characteristic orientations of the crystallographic standard triangle is modified. The computed square deviations of the stresses averaged on grains with the same crystallographic orientation but with different neighbour grains are a bit higher than the second order ones (inter-orientation scatter). Our numerical evaluations agree quantitatively with several observations and measures of the literature concerning stress and strain distribution in copper and austenite polycrystals submitted to low amplitude loadings. Hopefully, the given apparent Schmid factor distributions could help to better understand the observations of the plasticity mechanisms taking place at the free surface of polycrystals. To cite this article: M. Sauzay, C. R. Mecanique 334 (2006).     

Nonlinear behavior of matrix-inclusion composites under high confining pressure: application to concrete and mortar

This paper is devoted to a micromechanics-based simulation of the response of concrete to hydrostatic and oedometric compressions. Concrete is described as a composite made up of a cement matrix in which rigid inclusions are embedded. The focus is put on the role of the interface between matrix and inclusion which represent the interfacial transition zone (ITZ). A plastic behavior is considered for both the matrix and the interfaces. The effective response of the composite is derived from the modified secant method adapted to the situation of imperfect interfaces. To cite this article: T.H. Le et al., C. R. Mecanique 336 (2008).