Comparative study of two fluid model and dense discrete phase model for simulations of gas-solid hydrodynamics in circulating fluidized beds

Computational fluid dynamics(CFD)has become a valuable tool to study the complex gas-solid hydrodynamics in the circulating fluidized bed(CFB).Based on the two fluid model(TFM)under the Eulerian-Eulerian framework and the dense discrete phase model(DDPM)under the Eulerian-Lagrangian framework,this work conducts the comparative study of the gas-solid hydrodynamics in a CFB riser by these two different models.Results show that DDPM could be used to predict gas-solid hydrodynamics in the circulating fluidized bed,and there are differences between TFM and DDPM,especially in the radial distribution profiles of solid phase.Sensitivity analysis results show that the gas-solid drag model exhibits significant effects on the results for both the two models.The specularity coefficient and the restitution coefficient in the TFM,as well as the reflection coefficient and the parcel number in the DDPM,exhibit less impact on the simulated results.     

Torsional static and vibration analysis of functionally graded nanotube with bi-Helmholtz kernel based stress-driven nonlocal integral model

A torsional static and free vibration analysis of the functionally graded nanotube(FGNT)composed of two materials varying continuously according to the power-law along the radial direction is performed using the bi-Helmholtz kernel based stress-driven nonlocal integral model.The differential governing equation and boundary conditions are deduced on the basis of Hamilton’s principle,and the constitutive relationship is expressed as an integral equation with the bi-Helmholtz kernel.Several nominal variables are introduced to simplify the differential governing equation,integral constitutive equation,and boundary conditions.Rather than transforming the constitutive equation from integral to differential forms,the Laplace transformation is used directly to solve the integro-differential equations.The explicit expression for nominal torsional rotation and torque contains four unknown constants,which can be determined with the help of two boundary conditions and two extra constraints from the integral constitutive relation.A few benchmarked examples are solved to illustrate the nonlocal influence on the static torsion of a clamped-clamped(CC)FGNT under torsional constraints and a clamped-free(CF)FGNT under concentrated and uniformly distributed torques as well as the torsional free vibration of an FGNT under different boundary conditions.     

Analysis on nonlinear effect of unsteady percolation in the inhomogeneous shale gas reservoir

The nonlinear effects of unsteady multi-scale shale gas percolation,such as desorption,slippage,diffusion,pressure-dependent viscosity,and compressibility,are investigated by numerical simulation.A new general mathematical model of the problem is built,in which the Gaussian distribution is used to describe the inhomogeneous intrinsic permeability.Based on the Boltzmann transformation,an efficient semi-analytical method is proposed.The problem is then converted into a nonlinear equation in an integral form for the pressure field,and a related explicit iteration scheme is constructed by numerical discretization.The validation examples show that the proposed method has good convergence,and the simulation results also agree well with the results obtained from both numerical and actual data of two vertical fractured test wells in the literature.Desorption,slippage,and diffusion have significant influence on shale gas flows.The accuracy of the usual technique that the product of viscosity and compressibility is approximated as its value at the average formation pressure is examined.     

Active control of multiscale features in wall-bounded turbulence

This study experimentally investigates the impact of a single piezoelectric(PZT)actuator on a turbulent boundary layer from a statistical viewpoint.The working conditions of the actuator include a range of frequencies and amplitudes.The streamwise velocity signals in the turbulent boundary layer flow are measured downstream of the actuator using a hot-wire anemometer.The mean velocity profiles and other basic parameters are reported.Spectra results obtained by discrete wavelet decomposition indicate that the PZT vibration primarily influences the near-wall region.The turbulent intensities at different scales suggest that the actuator redistributes the near-wall turbulent energy.The skewness and flatness distributions show that the actuator effectively alters the sweep events and reduces intermittency at smaller scales.Moreover,under the impact of the PZT actuator,the symmetry of vibration scales’velocity signals is promoted and the structural composition appears in an orderly manner.Probability distribution function results indicate that perturbation causes the fluctuations in vibration scales and smaller scales with high intensity and low intermittency.Based on the flatness factor,the bursting process is also detected.The vibrations reduce the relative intensities of the burst events,indicating that the streamwise vortices in the buffer layer experience direct interference due to the PZT control.     

Asymptotics for filtration of polydisperse suspension with small impurities

A model for deep bed filtration of a polydisperse suspension with small impurities in a porous medium is considered.Different suspended particles move with the same velocity as the carrier water and get blocked in the pore throats due to the size-exclusion mechanism of particle retention.A solution of the model in the form of a traveling wave is obtained.The global exact solution for a multiparticle filtration with one high concentration and several low concentrations of suspended particles is obtained in an explicit form.The analytic solutions for a bidisperse suspension with large and small particles are constructed.The profiles of the retained small particles change monotony with time.The global asymptotics for the filtration of a polydisperse suspension with small kinetic rates is constructed in the whole filtration zone.     

低温超导核磁共振陀螺仪模型

建立了单工作物质的三自由度3He低温超导核磁共振陀螺仪结构;利用了量子力学和经典动力学,经过严密的力学分析和数学演算,给出了三轴陀螺仪的工作原理和结构示意图;针对陀螺仪的交叉轴角速率耦合问题,给出了附加磁场线圈解耦法和冗余设计结构,最后建立了基于超导量子干涉仪探测磁矩,并采用最小二乘估计法来推导陀螺进动频率.陀螺仪测速范围可以达到10-9～103 rad/s,漂移为10-4 (°)/h.该陀螺仪结合低温超导技术具有高精度的前景.需要进一步对超导量子干涉磁矩检测仪的精度与陀螺性能进行研究.     

一种Cartesian气动网格的自适应划分算法

提出了针对空间Cartesian结构化网格的自适应划分算法,以满足飞机总体设计过程中进行多学科设计优化时根据方案变化生成CFD计算网格的需要,使气动特性能够得以自动化分析,并且保证优化迭代过程的顺利进行.分析了Cartesian网格划分过程中的两个关键问题,并给出了解决方法.在此基础之上,从理论上推导了能够自动适应飞机布局形式以及几何尺寸任意变化的Cartesian网格生成算法和实际划分步骤,并给出了其中部分主要的计算公式,同时论述了自适应划分过程中的若干重要问题.阐述了自适应划分所需的飞机及其部件特征尺寸的获取途径.最后展示了几个由自适应划分算法生成的Cartesian网格实例以说明其实际应用.     

DEM investigation of macro-and micro-mechanical properties of rigid-grain and soft-chip mixtures

We investigated the macro-and micro-mechanical properties of rigid-grain and soft-chip mixtures(GCMs)through numerical simulations using the discrete element method.We present a novel framework for the discrete modeling of soft chips and rigid grains in conjunction with calibration processes.Several numerical triaxial tests were also performed on GCMs with 0%,10%,20%,and 30%volumetric chip contents,P.The simulation results demonstrate that increasing P leads to higher GCM toughness,higher deviatoric peak stress,and higher corresponding shear strain.Higher P also contributes to more volume contraction and less dilation.The friction angles at both the peak and residual state significantly increase with increasing P.In view of the micro-mechanical features,strong contact force chains develop along the loading direction,which results in considerable anisotropy in the peak and residual states.Both the formation of strong force chains and rotation of grains decrease with increasing P,whereas the grain sliding percentage increases.The tensile force is mobilized with shearing and higher P leads to less mobilization of the tensile force.These findings are useful for better understanding the internal structure of GCMs with different soft-chip contents,especially in granular mixture mechanics and geomechanics.     

Interaction effects on heat conduction and thermal stress in an infinite elastic plane with two circular holes

Numerous researches have focused on the physical behavior of an elastic material in the vicinity of a single hole under the assumption that the interaction effects arising from the introduction of multiple holes remain negligible if the holes are placed sufficiently far from each other.In an effort to understand hole interaction effects on heat conduction and thermal stress,we consider the case when two circular holes are embedded in an infinite elastic material and use complex variable methods together with numerical analysis to obtain solutions describing temperature and elastic fields in the vicinity of the two circular holes.The results indicate that the interaction effects on temperature distribution and stress strongly depend on the relative size of the two holes and the distance placed between them but not on the actual size of the holes.     

向蔡树棠先生学习什么——写在送别蔡先生这一天

题记蔡树棠教授离开我们已近半年,他的音容笑貌还时不时地浮现在我的脑海里.我刚从中国科学技术大学回来,蔡先生曾在该校工作多年,而那里的几位年青学子居然不知道蔡树棠这个名字,这使我有点悲哀.这里公布我写于送别他那天的短文<向蔡树棠先生学习什么>,希望大家记住这位曾在力学园地耕耘了半个多世纪的知名学者.     

Quasi-static tensile deformation and fracture behavior of a highly particle-filled composite using digital image correlation method

Polymer bonded explosives (PBXs) are highly particle-filled composite materials. This paper experimentally studies the tensile deformation and fracture behavior of a PBX simulation by using the semi-circular bending (SCB) test. The deformation and fracture process of a pre-notched SCB sample with a random speckle pattern is recorded by a charge coupled device camera. The displacement and strain fields on the observed surface during the loading process are obtained by using the digital image correlation method. The crack opening displacement is calculated from the displacement fields, the initiation and propagation of the crack are analyzed. In addition, the damage evolution and fracture mechanisms of the SCB sample are analyzed according to the strain fields and the correlation coefficient fields at different loading steps.     

Experimental analysis of crack tip fields in rubber materials under large deformation

A three-nested-deformation model is proposed to describe crack-tip fields in rubber-like materials with large deformation.The model is inspired by the distribution of the measured in-plane and out-of-plane deformation.The inplane displacement of crack-tip fields under both Mode I and mixed-mode(Mode I-II) fracture conditions is measured by using the digital Moire’ method.The deformation characteristics and experimental sector division mode are investigated by comparing the measured displacement fields under different fracture modes.The out-of-plane displacement field near the crack tip is measured using the three-dimensional digital speckle correlation method.     

基于数字图像相关性方法的脆性材料拉氏反分析实验技术

利用分离式Hopkinson压杆作为加载系统,借助超高速数字图像相关性分析方法,发展了长杆试件拉氏反分析实验技术,并用于研究脆性材料在小变形条件下的动态本构特性。通过超高速相机实时拍摄冲击加载下长杆试件变形的散斑图像,再对散斑图像进行数字图像相关性（digital image correlation,DIC）分析,获得长杆试件表面速度场和应变场。随后,以脆性材料PMMA（polymethyl methacrylate）为例,从DIC分析得到的速度场中提取出不同拉格朗日位置上质点速度时程曲线,构建路径线连接整个速度场,再结合零初始条件,数值求解得到了试件中的应力时程曲线,消去时间参数后,获得了脆性材料PMMA的动态应力-应变曲线,并与Hopkinson压杆实验和准静态压缩实验的结果进行了对比,揭示了PMMA材料在小应变条件下的黏弹性本构响应特征。     

基于显微CT图像的近场动力学建模与复合材料微结构破坏模拟

随着纤维增强复合材料应用领域的不断扩展且用量激增,亟需理清复合材料微观结构损伤对宏观力学性能影响的内在机制。因此,发展针对纤维增强复合材料微结构破坏过程的建模与高效模拟方法就显得十分重要。本文借助显微CT(Micro-computed Tomography)扫描技术,提出了一种基于显微CT图像中像素点离散的近场动力学建模与模拟方法。一方面,近场动力学作为一种由积分方程建模的非局部理论,便于采用基于空间点离散的数值计算方法,相比传统的连续介质力学能够更有效地模拟材料从连续变形到裂纹萌生与扩展(非连续变形)的全过程。另一方面,对显微CT图像使用像素点灰度阈值分割处理技术,能够快速建立含有复合材料原位微结构信息的空间点离散模型。该离散模型可以直接用于微结构破坏过程的近场动力学模拟,从而避免了传统的数值模拟技术需要依据像素点先建立光滑的几何模型、再划分成有限单元网格的复杂前处理过程,并且极大地保留了复合材料的原位组分分布信息。数值模拟结果表明,基于显微CT图像的近场动力学建模方法能够精确捕捉到复合材料微结构信息,并能准确模拟纤维增强复合材料的微结构破坏过程。     

Analysis of a multiaxial test on a C/C composite by using digital image correlation and a damage model

The “planar” digital image correlation technique needs a single CCD camera to acquire the surface patterns of a zone of a specimen in the underformed and deformed states. With these two images, one can determine in-plane displacement and strain fields. The digital image correlation technique used herein is based on Fast Fourier Transforms, which are very effective in reducing the computation cost. Its performance is assessed and discussed on artificial signals and in a real experimental situation. The technique is utilized to analyze experimental results of a plane shear experiment and validate a damage meso-model describing different degradations in a C/C composite material.     

<Emphasis Type="Italic">In-situ</Emphasis> AFM Experiments with Discontinuous DIC Applied to Damage Identification in Biomaterials

Natural materials (e.g. nacre, bone, and spider silk) exhibit unique and outstanding mechanical properties. This performance is due to highly evolved hierarchical designs. Building a comprehensive understanding of the multi-scale mechanisms that enable this performance represents a critical step toward realizing strong and tough bio-inspired materials. This paper details a multi-scale experimental investigation into the toughening mechanisms in natural nacre. By applying extended digital image correlation and other image processing techniques, quantitative information is extracted from otherwise prodominantly qualitative experiments. In situ three point bending fracture tests are performed to identify and quantify the toughening mechanisms involved during the fracture of natural nacre across multiple length scales. At the macro and micro scales, fracture tests performed in situ with a macro lens and optical microscope enable observation of spreading of damage outward from the crack tip. This spreading is quantified using an iso-contour technique to assess material toughness. At the nanoscale, fracture tests are performed in situ an atomic force microscope to link the larger-scale damage spreading to sliding within the tablet-based microstructure. To quantify the magnitude of sliding and its distribution, images from the in situ AFM fracture tests are analyzed using new algorithms based on digital image correlation techniques which allow for discontinuous displacement fields. Ultimately, this comprehensive methodology provides a framework for broad experimental investigations into the failure mechanisms of bio- and bio-inspired materials.     

韧性金属裂纹尖端损伤区内应变场测量——数字散斑相关法应用

本文用数字散斑相关方法测量了五种不同幂硬化指数韧性金属材料(铝和铜),双边裂纹尖端细观区域内应变场.对所得结果用韧性损伤模型进行了分析.在此法中以金属自然表面结构为散斑场,不同加载状态的散斑场进行比较,得到相对变形与应变.图象之间相关性 C 是变形参数或是位移及其导数的泛函.使其相关性 C 取最大值的试凑变形即为其真实变形场.这一方法在细观测量中应用得到满意的结果.     

Research on synchronous measurement technique of temperature and deformation fields using multispectral camera with bilateral telecentric lens

The hot-section parts easily occur the creep-fatigued interaction under the condition of mechanical-thermal coupled load during the period of service, which may lead to the damage of the parts, and therefore, the measurement and characterization of thermal-deformed fields of the parts are important to understand its damage process. Aiming at relevant demand, the bilateral telecentric-multispectral imaging system was established, the research of synchronous measurement technique of the temperature and deformation fields was developed. On the one hand, the measurement technology for surface temperature of the object was developed using the two-color images captured by the multispectral camera with bilateral telecentric lens and combined with colorimetric method. On the other hand, the 2D-DIC measurement technique of the multispectral camera was developed by conducting digital image correlation analysis using the blue light images before and after deformation, which can measure the high temperature deformation field of the object (the blue light images were filtered by multispectral camera). Results showed that the bilateral telecentric lens is used to replace the ordinary optical lens for imaging, which can effectively eliminate the distortion of the multispectral imaging system. Since the temperature measurement process of this measurement system is little affected by the emissivity of the object, therefore, it has excellent robustness. The thermal expansion coefficients of the nickel alloys are evaluated at the temperature ranges of 700–1000 ℃, indicating this system can achieve the synchronous and precise measurement of the temperature and deformation fields of the object.     

Modeling the effects of material non-linearity using moving window micromechanics

In the analysis of materials with random heterogeneous microstructure the assumption is often made that material behavior can be represented by homogenized or effective properties. While this assumption yields accurate results for the bulk behavior of composite materials, it ignores the effects of the random microstructure. The spatial variations in these microstructures can focus, initiate and propagate localized non-linear behavior, subsequent damage and failure. In previous work a computational method, moving window micromechanics (MW), was used to capture microstructural detail and characterize the variability of the local and global elastic response. Digital images of material microstructure described the microstructure and a local micromechanical analysis was used to generate spatially varying material property fields. The strengths of this approach are that the material property fields can be consistently developed from digital images of real microstructures, they are easy to import into finite element models (FE) using regular grids, and their statistical characterizations can provide the basis for simulations further characterizing stochastic response. In this work, the moving window micromechanics technique was used to generate material property fields characterizing the non-linear behavior of random materials under plastic yielding; specifically yield stress and hardening slope, post yield. The complete set of material property fields were input into FE models of uniaxial loading. Global stress strain curves from the FE–MW model were compared to a more traditional micromechanics model, the generalized method of cells. Local plastic strain and local stress fields were produced which correlate well to the microstructure. The FE–MW method qualitatively captures the inelastic behavior, based on a non-linear flow rule, of the sample continuous fiber composites in transverse uniaxial loading.     

Study of a Crack at a Fastener Hole by Digital Image Correlation

In this work the efficacy of using digital image correlation to determine stress intensity factors for a crack emanating from a fastener hole has been investigated. To this end a fatigue crack was grown in pure mode I from a 50 mm diameter hole in an Al 7010 alloy plate test-piece. This crack was then loaded elastically under several combinations of mixed mode (I + II) displacements. In each case, images of the sample surface before and after the deformation were recorded using a high resolution digital camera. The surface preparation consisted only of scratching the surface lightly with silicon carbide abrasive paper. The crack location and resulting displacements were then calculated using digital image correlation. The analytical displacement fields for a traction free crack under arbitrary loading conditions based on the Muskhelishvili’s complex function approach were fitted to the experimentally measured displacement fields and the mixed mode stress intensity factor was determined in each case. Good agreement with the nominal applied values was obtained. The uncertainty of the crack tip position has a major influence on the accuracy of the stress intensity factors and so the Sobel edge finding filter was successfully applied to experimental displacement fields to establish precisely the crack tip location. This paper was presented at the 2007 SEM Congress held in Springfield, Massachusetts, USA