垂直入射时无限长分层柱电磁散射的改进算法及应用

结合非均匀球粒子对平面波散射的散射场计算的改进算法,提出了平面波垂直入射无限长分层圆柱散射场快速稳定而有效的改进电磁散射算法。与已有算法相比,改进算法所能计算的无限长非均匀介质圆柱的尺寸参量突破10000,计算层数达到106,并且计算时间很短,最多仅为几秒。该算法可以用于不同的波段以及不同领域的任意无耗或吸收无限长圆柱体散射场的计算。最后将该算法应用于梯度折射率聚合物光纤(GI-POF)散射特性的研究,为非接触、在线测量聚合物光纤折射率分布提供了理论依据。     

空气喷射对旋转圆筒表面传质特性的影响

通过对旋转雷诺数R,喷射雷诺数Rej,喷嘴宽度和圆筒直径比B／d,喷射距离与喷嘴宽度之比L／B等影响旋转圆筒表面传质特性参数的实验研究,发现：在空气喷射作用下,由旋转产生的第一临界点消失,只存在第二临界点;平均舍伍德数Sh与Rej,B／d和L／B之间的关联式为Sh=0．515Re^2/3（B／d）^0.3（L／B）^0.1;空气喷射对圆筒表面传质的强化作用大于旋转的作用．     

Experimental study on wake flow structures of screen cylinders using PIV

Experiments were conducted in a water flume using Particle Image Velocimetry (PIV) to study the evolution of the vortical structures in the wakes of four types of screen cylinders at a Reynolds number of about 3200. The results were compared with that of a bare cylinder. The screen cylinders were made of stainless steel screen meshes of various porosities (37%, 48%, 61% and 67%) rolled into cylindrical shapes. Smoke wire flow visualisations in a wind tunnel were also conducted in support of the PIV tests. Depending on the porosity of the screen mesh, two vortex formation mechanisms for the screen cylinder wakes were identified. One was associated with wake instability and the other was associated with shear-layer (Kelvin-Helmholtz) convective instability which involved merging through pairing and tripling of small-scale vortices within the shear layers. The former was responsible for the formation of large-scale vortices in the bare cylinder and the screen cylinder wakes with 37% and 48% porosities, while the latter was responsible for the screen cylinder wakes with 61% and 67% porosities. The results also showed that with increasing porosity, the vortex formation region was extended farther downstream and the Reynolds shear stress, the Turbulent Kinetic Energy (TKE) and vortex intensity were decreased constantly.     

Finite deformation swelling of a temperature-sensitive hydrogel cylinder under combined extension-torsion

The swelling behavior of a temperature-sensitive poly-N-isopropylacrylamide(PNIPAM) hydrogel circular cylinder is studied subjected to combined extension-torsion and varied temperature. In this regard, a semi-analytical solution is proposed for general combined loading. A finite element(FE) analysis is conducted, subjecting a hydrogel cylinder to the combined extension-torsion and the varied temperature to evaluate the validity and accuracy of the solution. A user-defined UHYPER subroutine is developed and verified under free and constrained swelling conditions. The FE results illustrate excellent agreement with the semi-analytical solution. Due to the complexity of the problem, some compositions and applied loading factors are analyzed. It is revealed that for larger cross-linked density and larger ending temperature, the cylinder yields higher stresses and smaller radial swelling deformation. Besides, the radial and hoop stresses increase by applying larger twist and axial stretch. The hoop stresses intersect at approximately R/Rout = 0.58, where the hoop stress vanishes. Besides, the axial force has direct and inverse relationships with the axial stretch and the twist, respectively. However, the resultant torsional moment behavior is complex, and the position of the maximum point varies significantly by altering the axial stretch and the twist.     

Acoustic radiation force on an elastic cylinder in a Gaussian beam near an impedance boundary

Acoustic radiation force (ARF) is studied by considering an infinite elastic cylinder near an impedance boundary when the cylinder is illuminated by a Gaussian beam. The surrounding fluid is an ideal fluid. Using the method of images and the translation-addition theorem for the cylindrical Bessel function, the resulting sound field including the incident wave, its reflection from the boundary, the scattered wave from the elastic cylinder, and its image are expressed in terms of the cylindrical wave function. Then, we deduce the exact equations of the axial and transverse ARFs. The solutions depend on the cylinder position, cylinder material, beam waist, reflection coefficient, distance from the impedance boundary, and absorption in the cylinder. To analyze the effects of the various factors intuitively, we simulate the radiation force for non-absorbing elastic cylinders made of stainless steel, gold, and beryllium as well as for an absorbing elastic cylinder made of polyethylene, which is a well-known biomedical polymer. The results show that the impedance boundary, cylinder material, absorption in the cylinder, and cylinder position in the Gaussian beam significantly affect the magnitude and direction of the force. Both stable and unstable equilibrium regions are found. Moreover, a larger beam waist broadens the beam domain, corresponding to non-zero axial and transverse ARFs. More importantly, negative ARFs are produced depending on the choice of the various factors. These results are particularly important for designing acoustic manipulation devices operating with Gaussian beams.     

Application of wake oscillators to two-dimensional vortex-induced vibrations of circular cylinders in oscillatory flows

A nonlinear time-domain simulation model for predicting two-dimensional vortex-induced vibration (VIV) of a flexibly mounted circular cylinder in planar and oscillatory flow is presented. This model is based on the utilization of van der Pol wake oscillators, being unconventional since wake oscillators have typically been applied to steady flow VIV predictions. The time-varying relative flow–cylinder velocities and accelerations are accounted for in deriving the coupled hydrodynamic lift, drag and inertia forces leading to the cylinder cross-flow and in-line oscillations. The system fluid–structure interaction equations explicitly contain the time-dependent and hybrid trigonometric terms. Depending on the Keulegan–Carpenter number (KC) incorporating the flow maximum velocity and excitation frequency, the model calibration is performed, entailing a set of empirical coefficients and expressions as a function of KC and mass ratio. Parametric investigations in cases of varying KC, reduced flow velocity, cylinder-to-flow frequency ratio and mass ratio are carried out, capturing some qualitative features of oscillatory flow VIV and exploring the effects of system parameters on response prediction characteristics. The model dependence of hydrodynamic coefficients on the Reynolds number is studied. Discrepancies and limitations versus advantages of the present model with different feasible solution scenarios are illuminated to inform the implementation of wake oscillators as a computationally efficient prediction model for VIV in oscillatory flows.     

超声共振法测定圆柱体材料弹性常数的影响因素分析

为明确超声共振频谱分析法(RUS)测定圆柱体试样弹性常数的影响因素, 由此针对各向同性圆柱体开展实验, 观察圆柱体试样测量结果的分散性和稳定性, 且主要研究横纵比和有效模态数量对测量结果的影响, 并与拉伸试验的结果进行比较. 以测量结果标准差、变异系数和95%置信度下均值的区间估计3个指标来衡量结果的稳定性和分散性, 且从模态简并角度对实验结果进行分析. 最终结果表明, 当超声共振频谱分析法的有效模态数量在17左右、试样横纵比接近1时, 测量结果的稳定性和分散性较小.     

不可压粘性流体浸入边界直接力法及软件实现

浸入边界法通过在N-S方程中施加体积力模拟不可滑移固壁边界及动边界,避免生成复杂贴体网格及动网格,极大地节省了网格建模时间及动网格计算消耗。本文提出一种新型附加体积力简化计算方法,将简化附加体积力以源项形式嵌入动量方程迭代中,通过用户自定义函数对CFD软件FLUENT二次开发,实现了浸入边界法和通用流体力学求解器的耦合计算。通过静止圆柱和动圆柱绕流数值模拟进行了验证,并探讨了插值函数对计算精度的影响。研究表明,通过引入浸入边界模型,能够提高计算效率,并实现结构网格背景下复杂边界和动边界的高效建模。     

近壁面柱体涡激振动的迟滞效应

柱体涡激振动是典型的流固耦合问题,其响应规律标识码在升速流动和远离壁面条件下获得的. 而自然环境流动通常不断经历升速和降速过程,近壁面柱体的涡激振动可呈现与远离标识码体不同的响应特征. 本研究结合大型波流水槽,设计了具有微结构阻尼的柱体涡激振动装置. 基于量纲分析,开展系列水槽标识码验,通过同步测量柱体涡激振动位移时程和绕流流场变化,研究了升降流速作用下柱体涡激振动触发和停振的临界速度(即上临标识码临界速度)变化规律,探究了近壁面柱体涡激振动迟滞效应. 采用自下向上激光扫射的 PIV 流场测量系统,对比分析了固定柱体标识码振动柱体的绕流特征. 实验观测表明,近壁面柱体涡激振动触发的临界速度呈现随壁面间距比减小而逐渐减小的变化趋势;但标识码速条件下的涡激振动停振所对应的下临界速度却明显小于升速时的涡激振动触发所对应的上临界速度. 采用上临界与下临界约标识码差值可定量表征涡激振动迟滞程度,研究发现该值随着柱体间距比减小呈线性增大趋势. 涡激振动迟滞现象通常伴随振幅阶跃标识码阶跃值则随着间距比减小而非线性减小.      

Effect of rotation in magneto-thermoelastic transversely isotropic hollow cylinder with three-phase-lag model

This article deals with the various heat source responses in a transversely isotropic hollow cylinder under the purview of three-phase-lag (TPL) generalized thermoelasticity theory. In presence of magnetic field and due to the rotating behavior of the cylinder, the governing equations are redefined for generalized thermoelasticity with thermal time delay. In order to obtain the stress, displacement and temperature field, the field functions are expressed in terms of modified Bessel functions in Laplace transformed domain. When the outer radius of hollow cylinder tends to infinity, the corresponding results are discussed. Finally an appropriate Laplace transform inversion technique is adopted.