Tailoring particle shape for enhancing the homogeneity of powder mixtures:Experimental study and DEM modelling

The effect of particle shape modification on the segregation reduction of enzyme granules in laundry detergent powder mixtures was investigated,both experimentally and computationally using Deseret Element Method(DEM).The shape of modified enzyme particles was in such a way that the large and dense enzyme particles were layered by other fine particles in the detergent powder,by means of a process known in the literature as“seeded granulation”.It is found that the homogeneity of modified enzyme particles could be improved significantly comparing to the original spherical enzyme particles in powder mixtures.Overall,the results of this research demonstrated that the segregation-induced properties of the dense/spherical enzyme particles could be lowered by altering their shape,which could enable the enzyme particles to behave almost similar to other ingredients during the pile formation process.     

Improving the hydrodynamic performance of the SUBOFF bare hull model: a CFD approach

The main aims of this study are to investigate the hydrodynamic performance of an autonomous underwater vehicle(AUV),calculate its hydrodynamic coefficients,and consider the flow characteristics of underwater bodies.In addition,three important parts of the SUBOFF bare hull,namely the main body,nose,and tail,are modified and redesigned to improve its hydrodynamic performance.A three-dimensional(3D)simulation is carried out using the computational fluid dynamics(CFD)method.To simulate turbulence,the k-ωshear stress transport(SST)model is employed,due to its good prediction capability at reasonable computational cost.Considering the effects of the length-to-diameter ratio(LTDR)and the nose and tail shapes on the hydrodynamic coefficients,it is concluded that a hull shape with bullet nose and sharp tail with LTDR equal to 7.14 performs better than the SUBOFF model.The final proposed model shows lower drag by about 14.9%at u=1.5 m·s^-1.Moreover,it produces 8 times more lift than the SUBOFF model at u=6.1 m·s^-1.These effects are due to the attachment of the fluid flow at the tail area of the hull,which weakens the wake region.     

Analysis of a two-grid method for semiconductor device problem

The mathematical model of a semiconductor device is governed by a system of quasi-linear partial differential equations.The electric potential equation is approximated by a mixed finite element method,and the concentration equations are approximated by a standard Galerkin method.We estimate the error of the numerical solutions in the sense of the Lqnorm.To linearize the full discrete scheme of the problem,we present an efficient two-grid method based on the idea of Newton iteration.The main procedures are to solve the small scaled nonlinear equations on the coarse grid and then deal with the linear equations on the fine grid.Error estimation for the two-grid solutions is analyzed in detail.It is shown that this method still achieves asymptotically optimal approximations as long as a mesh size satisfies H=O(h^1/2).Numerical experiments are given to illustrate the efficiency of the two-grid method.     

Shape of free-fall arch in quasi-2D silo

In this study,the optical flow method is used to measure the velocity distribution of a granular flow in a rectangular quasi-two-dimensional silo.Using the velocity gradient,a free-fall arch(FFA)is obtained and its geometric characteristics are calculated.A parabola-shaped FFA structure is discovered above the orifice in the steady flow state.The shape of the FFA affects the flow rate through the orifice.Furthermore,as jamming begins to occur,the geometry of the FFA disappears gradually from both sides and then from the middle;finally,the FFA disappears completely in the state of jamming.As the boundary between finite-stress and stress-free regions,the FFA facilitates further studies regarding the discontinuity of the stress area above the orifice.     

A drag force formula for heterogeneous granular flow systems based on finite average statistical method

The existing drag models are mostly based on the assumption of homogenous fluidization.However,the use of a homogeneous drag model to predict a heterogeneous granular flow system will cause a deviation.In this study,we developed a drag force model based on the assumption of heterogeneous fluidization.To prevent weakening of the heterogeneous characteristics in the drag force formula,we propose a finite average statistical method to filter the information of the heterogeneous granular cluster.The filtered information was used to fit the modified drag formula,which can reflect the heterogeneity of the granular cluster considering different configurations.A comparison shows that the new proposed drag formula filtered by the finite average statistical method fits well with energy minimization multi-scale simulation results.     

A numerical investigation of CO2 dilution on the thermochemical characteristics of a swirl stabilized diffusion flame

The turbulent combustion flow modeling is performed to study the effects of CO_2 addition to the fuel and oxidizer streams on the thermochemical characteristics of a swirl stabilized diffusion flame. A flamelet approach along with three well-known turbulence models is utilized to model the turbulent combustion flow field. The k-ω shear stress transport(SST) model shows the best agreement with the experimental measurements compared with other models. Therefore, the k-ω SST model is used to study the effects of CO_2 dilution on the flame structure and strength, temperature distribution, and CO concentration. To determine the chemical effects of CO_2 dilution, a fictitious species is replaced with the regular CO_2 in both the fuel stream and the oxidizer stream. The results indicate that the flame temperature decreases when CO_2 is added to either the fuel or the oxidizer stream. The flame length reduction is observed at all levels of CO_2 dilution. The H radical concentration indicating the flame strength decreases, following by the thermochemical effects of CO_2 dilution processes. In comparison with the fictitious species dilution, the chemical effects of CO_2 addition enhance the CO mass fraction. The numerical simulations show that when the dilution level is higher, the rate of the flame length reduction is more significant at low swirl numbers.     

Static output feedback stabilization for second-order singular systems using model reduction methods

In this paper,the static output feedback stabilization for large-scale unstable second-order singular systems is investigated.First,the upper bound of all unstable eigenvalues of second-order singular systems is derived.Then,by using the argument principle,a computable stability criterion is proposed to check the stability of secondorder singular systems.Furthermore,by applying model reduction methods to original systems,a static output feedback design algorithm for stabilizing second-order singular systems is presented.A simulation example is provided to illustrate the effectiveness of the design algorithm.     

Geometric and material nonlinearities of sandwich beams under static loads

A new theory developed from extended high-order sandwich panel theory(EHSAPT)is set up to assess the static response of sandwich panels by considering the geometrical and material nonlinearities simultaneously.The geometrical nonlinearity is considered by adopting the Green-Lagrange-type strain for the face sheets and core.The material nonlinearity is included as a piecewise function matched to the experimental stress-strain curve using a polynomial fitting technique.A Ritz technique is applied to solve the governing equations.The results show that the stress stiffening feature is well captured in the geometric nonlinear analysis.The effect of the geometric nonlinearity in the face sheets on the displacement response is more significant when the stiffness ratio of the face sheets to the core is large.The geometric nonlinearity decreases the shear stress and increases the normal stress in the sandwich core.By comparison with open literature and finite element simulations,the present nonlinear EHSAPT is shown to be sufficiently precise for estimating the nonlinear static response of sandwich beams by considering the geometric and material nonlinearities simultaneously.     

第七届全国周培源大学生力学竞赛试题

一、小球在高脚玻璃杯中的运动(20分) 一半球形高脚玻璃杯,半径r=5cm,其质量m1=0.3kg,杯底座半径R=5cm,厚度不计,杯脚高度h=10cm.如果有一个质量m2=0.1kg的光滑小球自杯子的边缘由静止释放后沿杯的内侧滑下,小球的半径忽略不计.     

罐内爆炸成型技术在修复结晶器中的实验研究

随着连铸生产的发展,易损件结晶器铜管的消耗越来越大,使用报废的铜管结晶器越积越多.为重复利用结晶器,采用罐内爆炸成型技术做修复实验.针对200×200规格的结晶器铜管修复选择参数为:导爆索间距t=20mm;导爆索距离铜管外表20mm;炸药筒长度依据铜管长度加100mm确定.依据方案给出具体实验案例,成型后的结晶器经过测量表明修复后的结晶器完全可以再次利用.尤其,为实现厂内、连续、方便的修复结晶器,提出在1kgTNT当量爆炸罐体内实现.实验表明结晶器可以实现厂内便捷修复.     

Distortional hardening plasticity model for paperboard

A distortional hardening elasto-plastic model at finite strains suitable for modeling of orthotropic materials is presented. As a prototype material, paperboard is considered. An in-plane model is established. The model developed is motivated from non-proportional loading tests on paperboard where the paperboard is pre-strained in one direction and then loaded in the perpendicular direction. A softening effect is revealed in the pre-strained samples. The observed experimental findings cannot be accurately predicted by current models for paperboard. To be able to model the softening effects, a yield surface based on multiple hardening variables is introduced. It is shown that the model parameters can be obtained from simple uniaxial experiments. The model is implemented in a finite element framework which is used to illustrate the behavior of the model at some specific loading situations and is compared with strain fields obtained from Digital Image Correlation experiments.     

Momentum and heat transfer in the magnetohydrodynamic stagnation-point flow of a viscoelastic fluid toward a stretching surface

An analysis is made of the steady two-dimensional stagnation-point flow of an incompressible viscoelastic fluid over a flat deformable surface when the surface is stretched in its own plane with a velocity proportional to the distance from the stagnation-point. It is shown that for a viscoelastic conducting fluid of short memory (obeying Walters’ Bʹ model), a boundary layer is formed when the stretching velocity of the surface is less than the inviscid free-stream velocity and velocity at a point increases with increase in the Hartmann number. On the other hand an inverted boundary layer is formed when the surface stretching velocity exceeds the velocity of the free stream and the velocity decreases with increase in the Hartmann number. A novel result of the analysis is that the flow near the stretching surface is that corresponding to an inviscid stagnation-point flow when the surface stretching velocity is equal to the velocity of the free stream. Temperature distribution in the boundary layer is found when the surface is held at constant temperature and surface heat flux is determined. It is found that in the absence of viscous and Ohmic dissipation and strain energy in the flow, temperature at a point decreases with increase in the Hartmann number.     

Steady sliding frictional contact problem for a 2d elastic half-space with a discontinuous friction coefficient and related stress singularities

The steady sliding frictional contact problem between a moving rigid indentor of arbitrary shape and an isotropic homogeneous elastic half-space in plane strain is extensively analysed. The case where the friction coefficient is a step function (with respect to the space variable), that is, where there are jumps in the friction coefficient, is considered. The problem is put under the form of a variational inequality which is proved to always have a solution which, in addition, is unique in some cases. The solutions exhibit different kinds of universal singularities that are explicitly given. In particular, it is shown that the nature of the universal stress singularity at a jump of the friction coefficient is different depending on the sign of the jump.     

Interface damage modeled by spring boundary conditions for in-plane elastic waves

In-plane elastic wave propagation in the presence of a damaged interface is investigated. The damage is modeled as a distribution of small cracks and this is transformed into a spring boundary condition. First the scattering by a single interface crack is determined explicitly in the low frequency limit for the case of a plane wave normally incident to the interface. The transmission at an interface with a random distribution of small cracks is then determined and is compared to periodically distributed cracks. The cracked interface is then described by a distributed spring boundary condition. As an illustration the dispersion relation of the first modes in a thick plate with a damaged interface in the middle is given.     

Numerical study on bifurcations in the wake of a circular disk

Using Large-eddy simulation (LES), the dynamics in the wake of a circular disk with an aspect ratio of d/w = 5 is numerically studied. The circular disk is normal to the main flow, and Reynolds number ranges from 115 to 300. The first bifurcation is confirmed for Re = 120, leading to the steady state mode with a reflectional symmetry and a double-thread wake extending to the downstream. The Hopf bifurcation is found for Re = 152, and the planar symmetry is lost, which is different from that observed in the sphere wake; it is called the “reflectional-symmetry-breaking (RSB)” mode and the hairpin vortices in this mode are always shedding in a fixed orientation. The third bifurcation is captured for Re = 166, which is named the “standing wave (SW)” mode; the planar symmetry lost in RSB mode is recovered and the hairpin vortices are shedding in the oppositely sided orientations, unlike the ones observed in the sphere wake. The fourth bifurcation, referred to as “zigzag (ZZ)” mode, is observed for Re = 265 and the planar symmetry is lost again; the hairpin vortices are shedding in an irregular orientation and propagating in a zigzagged way; and a few small-scale structures begin to appear. Three different vortex shedding regimes are found in RSB, SW and ZZ modes, respectively. Results show that the recirculation region plays a significant role in the mode transitions, and the stagnation point of recirculation zone is conjectured to be the initial region causing the wake instability.     

基于PVM的网络并行有限元面向对象的设计

子结构是有限元并行计算常用的一种方法,本文采用面向对象的方法,首先对子结构进行了面向对象的设计,得到了其类层次结构图;然后针对工作站网络有限元并行计算环境。提出了基于PVM消息传递平台上的Shadow—Mirror数据传输模型,该模型在有限元并行计算数据传输时,充分发挥数据面向对象的特性,采用设置数据缓冲区、短消息合并等方法以缩短数据通信时间,并据此编制了相应的程序。计算结果表明,使用文中提出的面向对象的Shadow—Mirror数据传输模型可以得到较为理想的并行加速比,而且随着问题规模增大,并行加速比增高。本文研究内容为进一步开展基于工作站网络的并行有限元研究提供了一个可参考的基础。     

On the Effect of Anisotropy on the Stability of Convection in Rotating Porous Media

We investigate natural convection in an anisotropic porous layer subjected to centrifugal body forces. The Darcy model (including centrifugal and permeability anisotropy effects) is used to describe the flow and a modified energy equation (including the effects of thermal anisotropy) is used in the current analysis. The linear stability theory is used to evaluate the critical Rayleigh number for the onset of convection in the presence of thermal and mechanical anisotropy. It is found that the convection is stabilized when the thermal anisotropy ratio (which is a function of the thermal and mechanical anisotropy parameters) is increased in magnitude.     

Deformation of an ellipsoidal ferrogel sample in a uniform magnetic field

The elongation of a ferroelastic material sample (whose initial shape is a sphere or an ellipsoid of revolution) under the action of an external magnetic field is studied in an in approximation of small strains. For a sphere, there is a classical estimate obtained under the assumption that elongating in the direction of the field, it becomes a spheroid and the stress and strain fields remain uniform. In the present calculation, it is assumed that the body is an ellipsoid (a sphere in a particular case) only in the absence of an external field; the shape of the sample in the presence of a field is not specified in advance but is found from the condition of balance of surface forces (elastic and magnetic). For the spherical case, the problem is solved exactly: it is shown, that the contour of the deformed body is described by a third-order algebraic equation. The case where the initial configuration is an ellipsoid of revolution is studied numerically. It is shown that in all versions, the refined solution leads to an appreciable increase in the elongation of the sample compared to the classical estimate. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 46, No. 3, pp. 153–164, May–June, 2005.     

Experiments on the Onset of Impacting Motion Using a Pipe Conveying Fluid

The transition from stable periodic nonimpacting motion to impacting motion, due to variations of parameters, is observable in a wide range of vibro-impact systems. Recent theoretical studies suggest a possible scenario for this type of transition. A key element in the proposed scenario is fulfilled if the oscillatory motion involved in the transition is born in a supercritical Hopf bifurcation. If the onset of impacting motion is close to the Hopf bifurcation, the impacting motion is likely to be chaotic. A numerical simulation of a system of articulated pipes conveying fluid is used to illuminate the theory. An experimental setup is presented, where a cantilevered pipe conveying fluid is unilaterally constrained. Results from experiments are found to be in good qualitative agreement with the theory.     

Acoustical scattering identification with local impedance through a spectral approach

Acoustical scattering in waveguides is studied in this paper. The Wave Finite Element (WFE) approach is mainly used, since it allows the reduction of problems dealing with periodic waveguides. The paper deals with guided acoustical propagation, that is, propagation in a main direction is privileged. The scattering by a locally reacting lining is first studied. The liner can be characterised by its local impedance in this case. The equivalent surface impedance is therefore calculated. Then, scattering by a porous layer is considered. A full three-dimensional modelling of the lining is preferred since porous materials are bulk reacting. The scattering matrix of the lined part is computed, and acoustical scattering of high-order modes and conversion between modes are highlighted. The acoustic power attenuation is further evaluated. The response of ducts subjected to constraining boundary conditions is also calculated. Numerical results are presented and compared to those obtained with conventional approaches.