基于SPH方法粒子射流破岩数值模拟与实验研究

钻井液中加入体积分数为1%~3%的钢质粒子在钻头喷嘴处高速喷出冲击岩石,实现了粒子射流冲击和钻头机械联合破岩,有效提高了破岩效率。利用瞬态非线性动力学有限元模拟软件,基于光滑粒子流体动力学(smoothed particle hydrodynamics,SPH)方法,考虑流体对粒子射流冲击的影响,建立了粒子射流冲击破岩的物理模型,获得了粒子射流参数对破岩体积的影响规律,进行了室内实验验证,验证了SPH方法的有效性。结果表明:粒子射流冲击岩石表面形成规则的V型冲击坑;同条件下粒子射流破岩体积是水射流破岩体积的2~4倍;随着粒子射流冲蚀时间的增加,粒子射流破岩体积不断增加,但破岩效率降低;粒子射流压力大于10 MPa后,粒子射流破岩效率迅速增大;喷射角度大于6°后,破岩效率迅速减小。     

基于表面力模型的液滴冲击弹性基底SPH模拟

采用光滑粒子动力学SPH方法建立液滴冲击弹性基底的流固耦合数值模型,给出描述粘性流体和弹性固体运动的SPH离散方程和数值处理格式,引入人工耗散项来抑制标准SPH方法的数值震荡。为模拟液滴的表面张力效应,通过精确检测边界粒子,采用拉格朗日插值方法计算表面法向量和曲率,结合界面理论中的连续表面力CSF方法,建立了适用于自由表面液滴的表面力模型,方形液滴变形的模拟结果与拉普拉斯理论解吻合较好。随后,采用SPH流固耦合模型模拟1.0 mm直径水滴以不同速度(0.2 m/s～3.0 m/s)冲击两种薄板型基底,分析了基底弹性变形对液滴铺展、收缩以及回弹行为的影响。     

Development and experimental validation of a continuum micromechanics model for the elasticity of wood

This contribution covers the development and validation of a microelastic model for wood, based on a four-step homogenization scheme. At a length scale of several tens of nanometers, hemicellulose, lignin, and water are intimately mixed, and build up a polymer (polycrystal-type) network. At a length scale of around one micron, fiberlike aggregates of crystalline and amorphous cellulose are embedded in an contiguous polymer matrix, constituting the so-called cell wall material. At a length scale of about one hundred microns, the material softwood is defined, comprising cylindrical pores (lumen) in the cell wall material of the preceding homogenization step. Finally, at a length scale of several millimeters, hardwood comprises larger cylindrical pores (vessels) embedded in the softwood-type material homogenized before. Model validation rests on statistically and physically independent experiments: The macroscopic stiffness values (of hardwood or softwood) predicted by the micromechanical model on the basis of tissue-independent (‘universal’) phase stiffness properties of hemicellulose, amorphous cellulose, crystalline cellulose, lignin, and water (experimental set I) for tissue-specific composition data (experimental set IIb) are compared to corresponding experimentally determined tissue-specific stiffness values (experimental set IIa).     

考虑磁头表面高度不连续性气膜润滑的数值模拟与有效算法

随着磁头滑块的飞行高度不断降低,给气体润滑方程的数值求解带来了诸如计算时间过长、甚至计算发散等方面的问题。为了获得1Tbit/in2的存储密度,磁头滑块尾部的最小飞行高度接近1.5nm。本文基于作者提出的修正气膜润滑方程的线性流率(LFR)模型,考虑磁头滑块表面高度的不连续性,建立了基于有限体积法的气膜润滑方程离散格式,并把网格自适应技术与多重网格法应用到离散方程的迭代算法中,发展了可模拟最小飞行高度为0.5nm时磁头滑块压力分布的数值模拟方法与有效算法。文中以一个具有复杂表面形状的磁头滑块为例,检验了计算方法与算法的有效性。数值结果表明:在磁头滑块最小飞行高度较低时,必须要考虑滑块表面高度的不连续性,否则就得不到收敛的数值计算结果;与FK-Boltzmann模型相比,LFR模型具有较高的计算效率,采用网格自适应技术与多重网格法能有效地提高求解气膜润滑方程的计算效率。     

基于离散缝洞网络模型的缝洞型油藏混合多尺度有限元数值模拟

针对缝洞型油藏具有多尺度特征,基于离散缝洞网络模型,建立了Darcy/Stokes-Brinkman多尺度耦合数学模型,采用多尺度混合有限元方法,对缝洞型油藏的流体流动问题进行了研究。阐述了多尺度混合有限元方法的基本原理,推导得到了Darcy/Stokes-Brinkman多尺度模型的多尺度混合有限元计算格式。数值计算结果表明,在大尺度模型上进行计算能够捕捉到小尺度上的流动特征;相对于传统有限元,多尺度混合有限元能够捕捉小尺度上的非均质特征而具有更高的计算精度,在保证计算精度的同时能够减少计算量。     

Robust design and numerical simulation on drag reduction by a mixture film for liquid turbulence in a channel

An important way of increasing the speed and lowering the fuel consumption of ships is by decreasing the frictional drag. One of the most promising techniques for reducing drag is the use of air bubbles. The goal of this investigation is to establish a set of optimum robust parametric levels for drag reduction by a mixture (air–water) film in turbulent channel flow. Based on the conditions laid out by the Taguchi orthogonal array method, turbulent flows, with air bubbles injected into a channel, are simulated using commercial computational fluid dynamics software. The local shear stress on the upper wall is computed to evaluate the efficiency of drag reduction. Many factors can affect drag reduction. The factors investigated in this study are the rate of air injection, bubble size, area of air injection, flow speed, and measured position of the shear stress. These factors have been investigated through the analysis of variance, which has revealed that the rate of air injection and water flow speed dominate the efficiency of drag reduction by a mixture film. According to the results, the drag can be reduced by an average of 83.4%; and when the configuration of the parametric levels is optimum the maximum drag reduction of 88.5% is achieved. Copyright © 2008 John Wiley & Sons, Ltd.     

基于夏比冲击试验的材料失效模型参数

结合夏比冲击试验和ABAQUS显式动力数值模拟,对Q370d钢进行了Johnson-Cook失效模型参数研究。首先,在不考虑材料失效的情况下,通过3种不同厚度的无缺口试件冲击实验对有限元模型参数设置和材料本构模型的准确性进行了验证,同时还讨论了试件断裂区网格的合适尺寸;在此基础上,基于正交设计,通过大量的有限元数值模拟得到失效模型参数样本,利用回归分析求得冲击功与失效模型参数的回归方程组;最后结合夏比V型缺口冲击试验,求解Q370d钢的失效模型参数,并对断裂截面的力学特性进行了分析,可为工程应用提供参考。     

Validation of a numerical model for the simulation of an electrostatic powder coating process

Numerical modeling of a complete powder coating process is carried out to understand the gas-particle two-phase flow field inside a powder coating booth and results of the numerical simulations are compared with experimental data to validate the numerical results. The flow inside the coating booth is modeled as a three-dimensional turbulent continuous gas flow with solid powder particles as a discrete phase. The continuous gas flow is predicted by solving Navier–Stokes equations using a standard k–ε turbulence model with non-equilibrium wall functions. The discrete phase is modeled based on a Lagrangian approach. In the calculation of particle propagation, a particle size distribution obtained through experiments is applied. The electrostatic field, including the effect of space charge due to free ions, is calculated with the use of the user defined scalar transport equations and user defined scalar functions in the software package, FLUENT, for the electrostatic potential and charge density.     

Partition method for impact dynamics of flexible multibody systems based on contact constraint

The impact dynamics of a flexible multibody system is investigated. By using a partition method, the system is divided into two parts, the local impact region and the region away from the impact. The two parts are connected by specific boundary conditions, and the system after partition is equivalent to the original system. According to the rigid-flexible coupling dynamic theory of multibody system, system＇s rigid-flexible coupling dynamic equations without impact are derived. A local impulse method for establishing the initial impact conditions is proposed. It satisfies the compatibility con- ditions for contact constraints and the actual physical situation of the impact process of flexible bodies. Based on the contact constraint method, system＇s impact dynamic equa- tions are derived in a differential-algebraic form. The contact/separation criterion and the algorithm are given. An impact dynamic simulation is given. The results show that system＇s dynamic behaviors including the energy, the deformations, the displacements, and the impact force during the impact process change dramatically. The impact makes great effects on the global dynamics of the system during and after impact.     

基于整体位移模式的平面裂纹结构数值模拟方法

对于平面裂纹问题,针对扩展有限元法和无网格伽辽金法的不足,从结构的整体位移模式出发,提出了一种新的数值模拟方法。在整个求解域内构造其试探函数,并引入裂纹修正项描述裂尖处的奇异性和裂纹面的强间断特性;同时,提出了一种新的强制边界条件施加方法,通过引入位移边界水平集函数,将位移边界条件包含在近似位移场的表达式中,有效地解决了位移边界条件问题,减小了刚度矩阵的阶数,非常方便地消除了刚度矩阵的奇异性,降低了线性方程组的求解难度。含裂纹矩形平板结构的数值算例验证了该方法的有效性。     

Numerical simulation of polyhedral crystal growth based on a mathematical model arising from non–local thermomechanics

In this work we present some results of the numerical simulation of the growth of a crystal from its melt, taking into account faceting. The simulation is based on a numerical solution of a three–dimensional generalized Stefan problem. That problem arises from a non–local thermomechanical theory applied to a continuous system with an interface and embodies ideas from the dislocation theory of crystal growth. In the model, the crystal surface is an isotherm and the growth velocity of a crystal face depends on the velocities of the other faces and on the whole crystal configuration as well as on the temperature gradient. A front fixing formulation of the model is considered. This is a conservative form of the Isotherm Migration Method [6, 7, 8, 9, 10, 11] in spherical coordinates. The numerical solution is based on an explicit finite difference discretization of the resulting non–linear equations. We develop a theoretical analysis of the interface equations that drive the crystal face motion. Numerical results, showing evolution of complex crystals with configuration changing during the growth, are in accord with experimental results. Furthermore, numerical experiments offer useful information on the influence of certain parameters in the model on the growth process. Received: March 21, 1996     

An efficient numerical solution for linear stability of circular jet: A combination of Petrov–Galerkin spectral method and exponential coordinate transformation based on Fornberg's treatment

This paper presents the linear stability analysis of a round jet in a radially unbounded domain using a spectral Petrov–Galerkin scheme coped with exponential coordinate transformation based on Fornberg's treatment. A Fourier–Chebyshev Petrov–Galerkin spectral method is described for the computation of the linear stability equations based on half a Gauss–Lobatto mesh. Complex basis functions presented here are exponentially mapped as Chebyshev functions, which satisfy the pole condition exactly at the origin, and can be used to expand vector functions efficiently by using the solenoidal condition. The mathematical formulation is presented in detail focusing on the solenoidal vector field used for the approximation of the flow. The scheme provides spectral accuracy in the present cases and the numerical results are in agreement with former works. Copyright © 2008 John Wiley & Sons, Ltd.