复杂型腔充模中纤维取向的动态模拟

基于描述短纤维增强复合材料充模过程的气-固-液三相模型及同位网格有限体积法, 实现了纤维增强复合材料沿复杂型腔水平中面充模过程的动态模拟. 不仅得到了界面位置、各物理量的信息, 而且得到了纤维在型腔中的运动情况(包括纤维的平动和取向). 结果表明, 与沿型腔厚度方向纤维取向的表层-芯层结构不同, 纤维沿型腔水平中面的取向与型腔结构有关, 入口处纤维取向环绕型腔入口, 沿水平或竖直方向纤维取向与来流方向垂直, 型腔拐角处纤维取向指向拐点.     

通道中串列旋转圆柱的格子Boltzmann-虚拟区域方法的模拟

通过结合格子Boltzmann方法(LBM)和虚拟区域(Fictitiou sDomain)思想,建立格子Boltzmann-虚拟区域(LB-DF/FD)方法.采用两套网格系统,欧拉网格用于流体,拉格朗日网格用于固体.原有的LBM在计算运动固体的受力方面存在数据振荡,LB-DF/FD方法改进了此缺陷.为验证该方法,模拟圆柱绕流、圆形颗粒在无限长通道中平动及在无限大流场中转动三种情况,结果与其他数值解及理论解符合得很好.利用该方法模拟低雷诺数下通道中串列旋转圆柱周围的流场,分析圆柱间距(g)及雷诺数(Re)对流场结构的影响.给出Re=0.001,0.1和10下,0.2≤g≤8.0的流线结构、圆柱升力、阻力以及力矩等数值结果.结果表明,g对流场的结构及圆柱的受力有显著影响,Re对圆柱阻力及Stokes单元数目的影响较大.     

水中微小波纹圆柱体声散射低频共振调控

塑料类高分子材料甲基丙烯酸甲酯-亚克力(PMMA)圆柱中亚音速Rayleigh波低频隧穿共振可引起反向散射增强,在低频标准散射体设计等领域具有重要应用价值.提出一种微弱形变的规则波纹表面结构,可实现水中PMMA圆柱反向散射低频共振频率的无源调控.利用微扰法推导了水中微弱形变规则波纹圆柱反向散射低频共振频率偏移的近似解,讨论了波纹微扰系数、周期对规则波纹圆柱共振频率偏移的影响规律.基于Rayleigh波相位匹配方法分析了低频共振频率偏移的机理.研究表明:微弱形变规则波纹圆柱中亚音速Rayleigh波沿微弱形变波纹表面传播,与光滑圆柱体相比,传播路径的改变引起Rayleigh波传播相位变化,导致了Rayleigh波低频共振频率发生偏移.最后开展了微弱形变规则波纹圆柱体声散射特性水池实验,获取了其反向散射共振频率,明显观察到了规则波纹圆柱共振频率偏移现象,与理论预报结果吻合较好.     

电磁控制两自由度涡生振荡的机理研究

电介质溶液中,电磁场产生的电磁力可以控制流体的运动.本文对两自由度圆柱涡生振荡及其电磁控制机理进行了数值研究.将坐标原点建立在振动圆柱上,推导了非惯性参考系指数极坐标下的涡量流函数方程、初始边界条件及水动力表达式.对圆柱沿法向和流向的流场、受力和位移的相互影响和瞬时对应规律进行了讨论,结果表明,圆柱的涡生振荡同时受到尾涡脱落和圆柱位移的影响.其作用方式沿法向通过影响圆柱上下两侧剪切层的强度,沿流向通过改变圆柱尾部二次涡的强度,从而改变圆柱的受力和运动.其中圆柱位移的作用效果与尾涡脱落的作用效果相反且占主导.另外,在电磁力的作用下,分离点被消除,使得圆柱的尾涡和推吸壁面的效果被抑制,从而使振动的诱因被消除,圆柱迅速达到稳定状态,并在电磁推力的作用下,圆柱的位置向上游移动.     

布朗运动和湍流扩散作用下槽流中纤维悬浮流动特性的研究

为了探讨槽道中纤维悬浮湍流场特性,本文对修正的雷诺平均纳维-斯托克斯方程、含纤维项的湍动能和耗散率方程、纤维动力学方程以及纤维取向概率密度函数方程进行了数值研究,同时进行了相关实验以验证计算结果.研究结果表明,流场中尤其是近壁区域,纤维取向趋向于与流动方向一致,且该现象随着雷诺数和纤维浓度的减小以及纤维长径比的增大而更为明显.纤维在槽道中分布不均匀,但随着雷诺数的增大和纤维长径比的减小趋于均匀.相对于单相流,纤维悬浮流的流向平均速度剖面更陡峭,且剖面斜率随着纤维浓度、长径比的增大以及雷诺数的减小而变大,纤维的存在使湍流场的湍动能和雷诺应力减小,且减小程度随着纤维浓度和长径比的增大以及雷诺数的减小而增加.流场中的第一法向应力差小于0.05且远小于剪切应力.从壁面到中心,剪切应力增加而第一法向应力差减小.剪切应力和第一法向应力差都随着纤维浓度和长径比的增大而增大.随着雷诺数的增大,剪切应力增大而第一法向应力差减小.纤维浓度对于剪切应力和第一法向应力差的影响比纤维长径比更显著.     

晶格中心插入体的对称性及取向对二维声子晶体带隙的影响

以二维钢/气体系声子晶体为模型,采用平面波法研究了圆柱正方及六角晶格中心添加插入体的对称性及取向与带隙的关系,给出了四方、六方、八方及圆柱插入体结构的带隙分布图及带隙随柱体取向的变化关系图.发现在低填充率条件下,插入体的截面形状与晶格类型相同时最有利于能带简并态的分离而获得带隙,但填充率较高时,采用高对称性的插入体可以获得最宽的带隙.正方晶格中心插入体取向对带隙的影响要比在六角晶格中更为显著.对四方柱正方晶格声子晶体的研究表明,仅旋转原柱体要比在其中心插入柱体后旋转更容易获得低频宽带隙,单独运用添加柱体或旋转非圆柱体来降低晶格对称性以获取低频带隙的方法要比同时使用两种方法效果更好.此外,从机理上对计算结果进行了解释.     

过渡流下管束的流动、传热特性研究—串列双圆柱流场研究

本文从换热器管束中最基本的串列布置形式出发研究处于过渡流区域的小管径换热器的性能变化。通过基于作者研究开发的复合网格计算方法模拟了在低、中雷诺数条件下,串列双圆柱随中心距的变化下的流场、温度场的变化特性,并分析了其机理。计算结果表明,随着圆柱间中心距的逐渐增大,圆柱尾流呈现出从不稳定到稳定再到不稳定的变化规律。通过研究发现,两圆柱中间流场是否稳定与这一区域的u_m是否是连续的负值有着密切的关系,而流场的变化规律同两圆柱间的u_m的变化规律也十分相似。     

用Lattice Boltzmann方法模拟1/4圆腔内的定常层流运动

分析格子Boltzmann方法中二阶精度的曲线边界处理方法.应用格子Boltzmann方法及其边界处理方法模拟1/4圆腔内的定常层流运动,引入流线图和等涡线图分析流场随Re数的变化.并且发现当Re数在10~100区间内变化时,随着Re数的增大,顺时针旋转流场的涡心位置偏离x轴的角度逐渐减小,而逆时针旋转流场的涡心位置偏离x轴的角度却越来越大.     

圆柱和直桨叶突然启动瞬态流动的数值研究

采用基于动网格方法的有限体积法对圆柱和直桨叶突然启动引起的二维非定常不可压粘性流进行了数值模拟,给出了计算方法.通过计算得到了Re=5000和9500时圆柱突然启动后的流场随时间演化的非定常过程.将计算结果与实验结果进行了对比,两者吻合较好;采用大涡模拟对高雷诺数下的直桨叶突然旋转启动的流场进行了计算,得到了完成启动后叶片内部瞬态流场的分布结果.     

Statistical property of dynamic speckles produced by a rotating cylinder object

高斯光束照射以固定角速度旋转的粗糙圆柱体表面,反射空间形成随着被测表面运动变化的动态散斑.研究了照射面积大于圆柱曲面条件下在远场衍射区形成的动态散斑统计特性,得到了不同入射条件及不同圆柱时动态散斑强度起伏空间-时间归一化相关函数.结果表明:随着圆柱半径、旋转速度以及照射光斑有效面积的增大,动态散斑强度起伏相关性快速减小;当圆柱旋转速度恒定,动态散斑光强起伏相关函数的峰值随着探测点之间距离的增大而减小,但峰值的位置随之而增大;在近似点照射情况下随着圆柱半径的增大空间相关长度基本不变,而相关时间性明显增大.     

The Force Acting on a Cylinder in a Ring Flow of a Viscous Fluid with a Small Eccentric Displacement

A system consisting of two circular cylinders one inside the other with parallel axes is considered. The outer cylinder of radius R₂ is fixed, and the inner cylinder of radius R₁ rotates with a sufficiently large angular velocity. The region between the cylinders is filled with an incompressible viscous fluid and, in the case of coaxial cylinders, Couette flow along circular trajectories arises. Upon an eccentric small displacement of the axis of the inner cylinder, the symmetry of the flow is disturbed and a force exerted on the inner cylinder by the fluid is created. Within the ideal fluid model, the force depends linearly on the transverse velocities and accelerations of the cylinder. In a viscous fluid, the force depends on the previous motion of the cylinder. It is expressed in terms of the velocity functional by analogy with the Basset force acting on a ball moving in a viscous fluid with a variable velocity.     

Three-mode optomechanical system for angular velocity detection

We propose a scheme for measuring the angular velocity of absolute rotation using a three-mode optomechanical system in which one mode of the two-dimensional(2 D) mechanical resonator is coupled to an optical cavity. When the total system rotates, the Coriolis force acting on the 2 D mechanical resonator due to the absolute rotation will affect the mechanical motion and thus change the phase of the output field from the cavity. The angular velocity of the absolute rotation can be estimated by monitoring the spectrum of the output field from the cavity via homodyne measurement. The minimum measurable angular velocity, which is determined by the noise spectrum, is calculated. The working range of the gyroscope for measuring angular velocity is discussed.     

基于浸入边界-多松弛时间格子玻尔兹曼通量求解法的流固耦合算法研究

针对流固耦合问题,发展了基于浸入边界-多松弛时间格子玻尔兹曼通量求解法(immersed boundary method multi-relaxation-time lattice Boltzmann flux solver,IB-MRT-LBFS)的弱耦合算法.依据多尺度Chapman-Enskog展开,建立不可压宏观方程状态变量和通量与格子玻尔兹曼方程中粒子密度分布函数之间的关系;采用强制浸入边界法处理流固界面使固壁表面满足无滑移边界条件,根据修正的速度求解动量方程力源项;结构运动方程采用四阶龙格-库塔法求解.格子模型与浸入边界法的引入使流固耦合计算可以在笛卡尔网格下进行,无需生成贴体网格及运用动网格技术,简化了计算过程.数值模拟了单圆柱横向涡激振动、单圆柱及串列双圆柱双自由度涡激振动问题.结果表明,IB-MRT-LBFS能够准确预测圆柱涡激振动的锁定区间、振动响应、受力情况以及捕捉尾流场结构形态,验证了该算法在求解流固耦合问题的有效性和可行性.     

Closed-form solution for a pair of touching cylindrical conductors in an external electric field

We give expressions for the electric potential and field around two conducting cylinders in contact, in closed form and in terms of elementary functions. The polarizability tensor, and hence the torque acting on the cylinder pair, follow from these results. The surface charge density, and the total charge per unit length on each cylinder, can also be evaluated in closed form. The equal and opposite forces on the two cylinders are found in the case of equal cylinder radii.     

Lattice Boltzmann Simulation of Sedimentation of a Single Elastic Dumbbell in a Newtonian Fluid

Based on the lattice Boltzmann method (LBM), the sedimentations of a single elastic dumbbell in a Newtonian fluid under different initial positions and orientations, and also that of the elastic dumbbells with different free lengths of the spring under the same initial conditions have been simulated. All of the numerical results show that the final orientations of the elastic dumbbells are in the same horizontal direction, and the final positions of their centroids are all on the centerline of the tube no matter what the initial positions and orientations of the elastic dumbbell or the free lengths of the spring are. When the elastic dumbbell finally falls down vertically, the two circular cylinders of the elastic dumbbell rotate around their own symmetry-axis respectively, and their angular velocities are equal but opposite to each other. For the sedimentations of the elastic dumbbells with different free lengths of the spring, the shorter of the free length is, the faster the final angular velocity and vertical velocity of the circular cylinder will be.     

串列双圆柱时均流场特性分析

对Re=12 000,间距比L/D=1.167,2.333,3.500和4.667的串列双圆柱后方速度场进行了实验测量,分析串列双圆柱后方不同剖面处的速度分布规律和湍流度分布规律.并通过流函数理论模型对小间距比串列双圆柱后方流场进行了分析,得出如下结论:与单圆柱相比,当串列双圆柱间距比较小时,后方圆柱的自由剪切层明显向内偏移.随着间距比的变大,由于前方圆柱的尾流对后方圆柱的干扰,后方圆柱的自由剪切层变得越来越模糊.对于间距比较大的串列双圆柱,其对后方流场的扰动较强,致使后方流场湍流强度和最大速度衰减量较大.通过流函数理论模型分析发现,在小间距比条件下,串列双圆柱由于两个圆柱的相互干扰,使得圆柱后方涡相互靠近,并且后方涡向外倾斜的角度也减小,从而导致了自由剪切层向内侧偏移.      

Lattice Boltzmann Simulation of Sedimentation of a Single Elastic Dumbbell in a Newtonian Fluid

Based on the lattice Boltzmann method (LBM), the sedimentations of a single elastic dumbbell in a Newtonian fluid under different initial positions and orientations, and also that of the elastic dumbbells with different free lengths of the spring under the same initial conditions have been simulated. All of the numerical results show that the final orientations of the elastic dumbbells are in the same horizontal direction, and the final positions of their centroids are all on the centerline of the tube no matter what the initial positions and orientations of the elastic dumbbell or the free lengths of the spring are. When the elastic dumbbell finally falls down vertically, the two circular cylinders of the elastic dumbbell rotate around their own symmetry-axis respectively, and their angular velocities are equal but opposite to each other. For the sedimentations of the elastic dumbbells with different free lengths of the spring, the shorter of the free length is, the faster the final angular velocity and vertical velocity of the circular cylinder will be.     

Exact Solutions for the Flow of Fractional Maxwell Fluid in Pipe-Like Domains

This paper presents an analysis of unsteady flow of incompressible fractional Maxwell fluid filled in the annular region between two infinite coaxial circular cylinders. The fluid motion is created by the inner cylinder that applies a longitudinal time-dependent shear stress and the outer cylinder that is moving at a constant velocity. The velocity field and shear stress are determined using the Laplace and finite Hankel transforms. Obtained solutions are presented in terms of the generalized G and R functions. We also obtain the solutions for ordinary Maxwell fluid and Newtonian fluid as special cases of generalized solutions. The influence of different parameters on the velocity field and shear stress is also presented using graphical illustration. Finally, a comparison is drawn between motions of fractional Maxwell fluid, ordinary Maxwell fluid and Newtonian fluid.     

Anisotropy of ultrasonic velocity and elastic properties in normal human myocardium.

Measurements of ultrasonic quasilongitudinal velocity were made in the muscle fiber plane of excised human myocardium. Multiple adjacent planes across the left ventricular wall were interrogated to assess the transmural dependence of velocity. For each measurement plane, data were obtained in 2-deg increments through the full 360 deg relative to the myofibers. An approximate 1.3% magnitude of anisotropy was observed with maximum velocity along the muscle fibers and minimum velocity perpendicular to the muscle fibers. The known transmural shift in myofiber orientation was evidenced in the anisotropy of velocity as angular shifts between plots obtained from adjacent transmural planes within the same specimen. Measured values of velocity and density were used to estimate the effective C33 and C11 elastic constants of a thin layer of normal myocardium.