Turbulent collision of inertial particles: point-particle based, hybrid simulations and beyond

Point-particle based direct numerical simulation (PPDNS) has been a productive research tool for studying both single-particle and particle-pair statistics of inertial particles suspended in a turbulent carrier flow. Here we focus on its use in addressing particle-pair statistics relevant to the quantification of turbulent collision rate of inertial particles. PPDNS is particularly useful as the interaction of particles with small-scale (dissipative) turbulent motion of the carrier flow is mostly relevant. Furthermore, since the particle size may be much smaller than the Kolmogorov length of the background fluid turbulence, a large number of particles are needed to accumulate meaningful pair statistics. Starting from the relative simple Lagrangian tracking of so-called ghost particles, PPDNS has significantly advanced our theoretical understanding of the kinematic formulation of the turbulent geometric collision kernel by providing essential data on dynamic collision kernel, radial relative velocity, and radial distribution function. A recent extension of PPDNS is a hybrid direct numerical simulation (HDNS) approach in which the effect of local hydrodynamic interactions of particles is considered, allowing quantitative assessment of the enhancement of collision efficiency by fluid turbulence. Limitations and open issues in PPDNS and HDNS are discussed. Finally, on-going studies of turbulent collision of inertial particles using large-eddy simulations and particle- resolved simulations are briefly discussed.

The effects of surface tension on the elastic properties of nano structures

In the absence of external loading, surface tension will induce a residual stress field in the bulk of nano structures. However, in the prediction of mechanical properties of nano structures, the elastic response of the bulk is usually described by classical Hooke’s law, in which the aforementioned residual stress was neglected in the existing literatures. The present paper investigates the influences of surface tension and the residual stress in the bulk induced by the surface tension on the elastic properties of nano structures. We firstly present the surface elasticity in the Lagrangian and the Eulerian descriptions and point out that even in the case of infinitesimal deformations the reference and the current configurations should be discriminated; otherwise the out-plane terms of surface displacement gradient, associated with the surface tension, may sometimes be overlooked in the Eulerian descriptions, particularly for curved and rotated surfaces. Then, the residual stress in the bulk is studied through the non-classical boundary conditions and used to construct the linear elastic constitutive relations for the bulk material. Finally, these relations are adopted to analyze the size-dependent properties of pure bending of Al nanowires. The present results show that surface tension will considerably affect the effective Young’s modulus of Al nanowires, which decrease with either the decrease of nanowires thickness or the increase of the aspect ratio.

Subgrid scale fluid velocity timescales seen by inertial particles in large-eddy simulation of particle-laden turbulence

Large-eddy simulation (LES) has emerged as a promising tool for simulating turbulent flows in general and, in recent years,has also been applied to the particle-laden turbulence with some success (Kassinos et al., 2007). The motion of inertial particles is much more complicated than fluid elements, and therefore, LES of turbulent flow laden with inertial particles encounters new challenges. In the conventional LES, only large-scale eddies are explicitly resolved and the effects of unresolved, small or subgrid scale (SGS) eddies on the large-scale eddies are modeled. The SGS turbulent flow field is not available. The effects of SGS turbulent velocity field on particle motion have been studied by Wang and Squires (1996), Armenio et al. (1999), Yamamoto et al. (2001), Shotorban and Mashayek (2006a,b), Fede and Simonin (2006), Berrouk et al. (2007), Bini and Jones (2008), and Pozorski and Apte (2009), amongst others. One contemporary method to include the effects of SGS eddies on inertial particle motions is to introduce a stochastic differential equation (SDE), that is, a Langevin stochastic equation to model the SGS fluid velocity seen by inertial particles (Fede et al., 2006; Shotorban and Mashayek, 2006a; Shotorban and Mashayek, 2006b; Berrouk et al., 2007; Bini and Jones, 2008; Pozorski and Apte, 2009).However, the accuracy of such a Langevin equation model depends primarily on the prescription of the SGS fluid velocity autocorrelation time seen by an inertial particle or the inertial particle–SGS eddy interaction timescale (denoted by $\delt T_{Lp}$ and a second model constant in the diffusion term which controls the intensity of the random force received by an inertial particle (denoted by C_0, see Eq. (7)). From the theoretical point of view, dTLp differs significantly from the Lagrangian fluid velocity correlation time (Reeks, 1977; Wang and Stock, 1993), and this carries the essential nonlinearity in the statistical modeling of particle motion. dTLp and C0 may depend on the filter width and particle Stokes number even for a given turbulent flow. In previous studies, dTLp is modeled either by the fluid SGS Lagrangian timescale (Fede et al., 2006; Shotorban and Mashayek, 2006b; Pozorski and Apte, 2009; Bini and Jones, 2008) or by a simple extension of the timescale obtained from the full flow field (Berrouk et al., 2007). In this work, we shall study the subtle and on-monotonic dependence of $\delt T_{Lp}$ on the filter width and particle Stokes number using a flow field obtained from Direct Numerical Simulation (DNS). We then propose an empirical closure model for $\delta T_{Lp}$. Finally, the model is validated against LES of particle-laden turbulence in predicting single-particle statistics such as particle kinetic energy. As a first step, we consider the particle motion under the one-way coupling assumption in isotropic turbulent flow and neglect the gravitational settling effect. The one-way coupling assumption is only valid for low particle mass loading.

Smoothed Particle Hydrodynamics (SPH): an Overview and Recent Developments

Smoothed particle hydrodynamics (SPH) is a meshfree particle method based on Lagrangian formulation, and has been widely applied to different areas in engineering and science. This paper presents an overview on the SPH method and its recent developments, including (1) the need for meshfree particle methods, and advantages of SPH, (2) approximation schemes of the conventional SPH method and numerical techniques for deriving SPH formulations for partial differential equations such as the Navier-Stokes (N-S) equations, (3) the role of the smoothing kernel functions and a general approach to construct smoothing kernel functions, (4) kernel and particle consistency for the SPH method, and approaches for restoring particle consistency, (5) several important numerical aspects, and (6) some recent applications of SPH. The paper ends with some concluding remarks.

光滑粒子动力学方法中粒子分布与数值稳定性分析

光滑粒子动力学(SPH)作为一种拉格朗日型无网格粒子方法,已经成功地应用于包括含多相流动界面以及移动边界的可压缩和不可压缩流体运动的研究中.通过对Poiseuille流动的深入研究,探索了SPH方法中粒子分布对计算精度的影响,揭示了一种因为粒子不规则分布而导致的数值不稳定现象.研究显示,这种数值不稳定性起源于SPH方法粒子近似过程中的不连续性.使用了一种新的粒子近似格式以确保SPH方法中粒子近似的连续性.计算结果表明,这种新的粒子近似格式对于规则和不规则的粒子分布都能得到稳定精度的结果.

节流器内液-固两相流固体颗粒冲蚀数值模拟

建立了考虑颗粒碰撞的颗粒冲蚀计算模型,该数学模型包括:在Eulerian坐标系下求解连续相流场;在Lagrangian坐标系下运用离散颗粒硬球模型求解颗粒碰撞;应用半实验关联式求解颗粒冲蚀速率。对水力加砂压裂施工中节流器内液-固两相流的固体颗粒运动和冲蚀特性进行了数值模拟。计算结果表明,固体颗粒密集于节流器入口到出口的一段狭长区域内,冲蚀速率随流体速度呈指数性变化。颗粒直径越大,冲蚀速率也越大。节流器内冲蚀最严重的位置发生在距离节流器出口上边缘10mm以内的局部区域。

Localization and tracking of phonating finless porpoises using towed stereo acoustic data-loggers

Cetaceans produce sound signals frequently. Usually, acoustic localization of cetaceans was made by cable hydrophone arrays and multichannel recording systems. In this study, a simple and relatively inexpensive towed acoustic system consisting of two miniature stereo acoustic data-loggers is described for localization and tracking of finless porpoises in a mobile survey. Among 204 porpoises detected acoustically, 34 individuals (similar to 17%) were localized, and 4 of the 34 localized individuals were tracked. The accuracy of the localization is considered to be fairly high, as the upper bounds of relative distance errors were less than 41% within 173 m. With the location information, source levels of finless porpoise clicks were estimated to range from 180 to 209 dB re 1 mu Pa pp at 1 m with an average of 197 dB (N=34), which is over 20 dB higher than that estimated previously from animals in enclosed waters. For the four tracked porpoises, two-dimensional swimming trajectories relative to the moving survey boat, absolute swimming speed, and absolute heading direction are deduced by assuming the animal movements are straight and at constant speed in the segment between two consecutive locations.

小波变换和神经网络在化学中的应用

本论文分为两部分：小波变换和神经网络在化学中的应用。小波变换是新近出现的数学方法。近年来，在化学中得到广泛的应用，本论文介绍小波变换和我分辨分析的原理和方法，并将其应用到信号的压缩和滤噪。在研究中提出了常用小波变换数据压缩的三种方法，将紧支集小波和正交三次B－样条小波压缩4－苯乙基邻苯二甲酸酐的红外光谱数据进行了对比，计算表明正交三次B－样条小波变换方法效果较好，而在全部保留模糊信号及只保留锐化信号中数值较大的系数时，压缩比大而重建光谱数据与原始光谱数据间的均方差较小。应用小波变换对表面等离子体子共振（Surface Plasmon Resonance，SPR）仪的信号进行滤噪处理，利用SPR仪器信号和噪音的频率特性而将其分离，取得良好效果。本文对神经网络在化学中应用进行了较深入的研究，并对影响神经网络的诸多因素进行了探讨。在神经网络和多元回归等在化学应用中过多的变量会导致数学模型的预测结果变差，因而选择合适变量是很重要的。本文对比了传统的统计方法（前进选择法，后退剔除法，逐步回归法），Leaps-and-Bounds回归法，正交变换法，主成分分析以及最新的优化技术遗传算法，得到了一些有意义的结果。同时提出了组合算法和前进选择法的得合算法，结果表明这种算法在一定程度上避免了局部最优且减少了计算量。本论文还利用上述方法进行了一些定量结构活性相关性研究，主要内容：1）环境中有毒有害有机物（苯酚，苯胺，硝基苯）的定量结构活性相性研究。2）抗艾滋病类药物（HEPT）定量结构活性相关性研究。3）抗肿瘤类药物（氮芥子气类化合物和2－甲酸吡啶缩氨基硫脲类化合物定量结构活性相关性研究。4）苯酚和苯胺类衍生物色谱比移值预测。5）将神经网络用于茶叶的分类。

Java程序中数组越界和空指针错误的静态分析

介绍了静态分析的相关背景知识,对Java程序中数组越界和空指针错误的常见表现形式进行了归类,并通过实验评估了现有的Java静态分析工具功能,指出了其在跨过程分析中的不足,并对如何改进进行了讨论.

A four-phase confocal elliptical cylinder model for predicting the effective thermal conductivity of coated fibre composites

A four-phase confocal elliptical cylinder model is proposed from which a generalised self-consistent method is developed for predicting the thermal conductivity of coated fibre reinforced composites. The method can account for the influence of the fibre section shape ratio on conductivity, and the physical reasonableness of the model is demonstrated by using the fibre distribution function. An exact solution is obtained for thermal conductivity by applying conformal mapping and Laurent series expansion techniques of the analytic function. The solution to the three-phase confocal elliptical model, which simulates composites with idealised fibre-matrix interfaces, is arrived at as the degenerated case. A comparison with other available micromechanics methods, Hashin and Shtrikman's bounds and experimental data shows that the present method provides convergent and reasonable results for a full range of variations in fibre section shapes and for a complete spectrum of the fibre volume fraction. Numerical results show the dependence of the effective conductivities of composites on the aspect ratio of coated fibres and demonstrate that a coating is effective in enhancing the thermal transport property of a composite. The present solutions are helpful to analysis and design of composites.

Application of the characteristic CIP method to a shallow water model on the sphere

Semi-implicit algorithms are popularly used to deal with the gravitational term in numerical models. In this paper, we adopt the method of characteristics to compute the solutions for gravity waves on a sphere directly using a semi-Lagrangian advection scheme instead of the semi-implicit method in a shallow water model, to avoid expensive matrix inversions. Adoption of the semi-Lagrangian scheme renders the numerical model always stable for any Courant number, and which saves CPU time. To illustrate the efficiency of the characteristic constrained interpolation profile (CIP) method, some numerical results are shown for idealized test cases on a sphere in the Yin-Yang grid system.

耗散粒子动力学处理复杂固体壁面的一种有效方法

耗散粒子动力学(dissipative particle dynamics,DPD)作为一种介观尺度拉格朗日型粒子方法,已经成功地应用于微纳米流动和生化科技的研究中. 复杂固体壁面的处理和壁面边界条件的实施一直是DPD方法发展及应用的一个障碍. 提出了处理复杂固体壁面的一种新的方法. 复杂固体区域通过冻结随机分布并且达到平衡状态的DPD粒子代表;所冻结的DPD粒子位于临近流动区域的一个截距内;在靠近固体壁面的流动区域中设置流动反弹层,当流动DPD粒子进入此流动层后反弹回流动区域. 应用这种固体壁面处理方法.

Role of the N

New measurement by CELSIUS-WASA Collaboration on the pp →pnπ+ reaction reveals clear evidence for the presence of the Roper resonance N∗(1440) which has been ignored in previous theoretical calculations. In this article, based on an effective Lagrangian approach and available knowledge on the Roper resonance, we investigate the role of the Roper resonance for the pp→pnπ+ reaction. It is found that the contribution from the Roper resonance N∗(1440) becomes significant for kinetic energy above 1.1 GeV, consistent with the new experimental observation. The t -channel σ-meson exchange is dominant for the production of the Roper resonance.

Evidence of N-*(1535) resonance contribution in the pn -> d phi reaction

The N ∗(1535) resonance contributions to the pn → dφ reaction are evaluated in an effective Lagrangian model. The π-, η-, and ρ-meson exchange are considered. It is shown that the contributions from π- and ρ-meson exchange are dominant, while the contribution from η-meson exchange is negligibly small. Our theoretical results reproduce the experimental data of both total cross section and angular distribution well. This is more evidence that the N ∗(1535) resonance has a large s ¯s component leading to a large coupling to Nφ, which may be the real origin of the Okubo-Zweig-Iizuka rule violation in the πN and pN reactions.

N*(1535) contribution in pp -> pp eta ' and pn -> d phi

In an effective Lagrangian model we find that the N*(1535) resonance contribution might be important to the interpretation of the present data of the pp -> pp eta' and pn -> d phi reactions. The strong coupling strength of N*(1.535) to eta' and phi are indicated, and the possible implication to the intrinsic component of N*(1535) is explored. These results may provide hints to the real origin of the OZI rule violation. It is stressed that further measurements could be performed at the Cooling Storage Ring (CSR) at Lanzhou of China.