基于格子Boltzmann大密度比模型的多孔介质内气泡动力学行为数值模拟

基于格子Boltzmann两相流大密度模型，研究气泡穿过多孔介质的动力学行为。研究发现：当孔隙率较大时，气泡只变形不破裂，能完整地通过多孔介质；而孔隙率较小时，气泡变形更加剧烈且发生破裂，穿过多孔介质所需的时间更长。另外，当障碍物表面接触角（θ）较小时，气泡均能完整地通过多孔介质，随着接触角的增大，气泡开始发生破裂，且θ越大，气泡破裂越严重，通过多孔介质的气泡剩余质量越小。数值结果还表明随着Eotvos数（Eo）增大，表面张力所占比重减小，气泡破裂越严重，穿过多孔介质的气泡剩余质量越小。最后，对比发现，多孔介质润湿性对气泡剩余质量影响最大，其次是孔隙率，而Eo数的影响最小。     

气泡碰撞过程中形变及破碎现象分析

以水气两相流为研究对象, 采用大涡模拟和流体体积法, 结合破碎准则和第三代涡识别方法对双气泡碰撞过程进行数值模拟。采用单一变量法, 研究气泡直径比、相对偏心距及气泡间相对距离对气泡破碎程度的影响。双气泡碰撞过程中, 两气泡的直径越接近时, 碰撞过程中的气泡破碎程度越弱, 偏心距的变化对气泡破碎程度影响不大。而相对距离小于1时, 随着其增大, 气泡碰撞过程中的破碎程度更加明显。当气泡相对距离大于1时, 气泡破碎程度趋于平缓。研究表明: 第三代涡识别方法能够很好地捕捉两相流湍流流场中漩涡位置, 并能敏锐地反映出湍流的变化。     

二维非定常Sine-Gordon方程辛算法及其孤子数值模拟

在矩形域[-a,a]×[-a,a]内对微分算子L=(ə²)/(əx²)+(ə²)/(əy²)用5点差分格式将二维非定常Sine Gordon方程离散化为一个2×799²阶非线性Hamilton系统.对该系统使用Euler中心格式,得到一个非线性方程组.对此方程组建立迭代解法并给出了这个迭代方法的收敛条件和收敛速度.Sine Gordon方程单孤子和双孤子的数值模拟试验显示该辛算法是有效的.     

周期力调制噪声驱动下单模激光系统的多重随机共振

在周期力调制噪声驱动下单模激光系统的光强方程中加入调幅波, 用线性化近似方法计算了系统的光强关联函数和输出信噪比, 并对信噪比进行数值计算和分析, 发现低频调制频率Ω、高频载波频率ω和周期力频率Ω_λ对系统的输出信噪比有很大的影响. 具体表现为信噪比R 随低频调制频率Ω 的变化过程中出现了多重随机共振和极强的单峰共振, 当Ω << ω 时, 系统出现的是多峰共振, 且随着Ω_λ 增加, 共振峰间的距离增大, 峰值位置不变; 当Ω → ω 时, 输出信噪比R迅速增大, 而Ω_λ 的影响被削弱甚至可以忽略, 多峰共振消失; 当Ω = ω 时, 系统出现了极强的单峰共振. 此外, 信噪比随周期力频率的变化呈现振幅减小的多重随机共振, 而随载流频率的变化出现单峰随机共振.     

非球形效应对强声场中次Bjerknes力的影响

考虑了非球形气泡在声场中的形状振动,推导了非球形气泡和球形气泡之间的次Bjerknes力方程,数值模拟了声场中非球形气泡和球形气泡之间的次Bjerknes力和两个球形气泡之间的次Bjerknes力,并对非球形气泡和球形气泡之间的次Bjerknes力的影响因素进行了分析讨论.研究结果表明:当驱动声压振幅大于非球形气泡的Black阈值且又能使得非球形气泡稳定振动时,在第一个声驱动周期内,非球形气泡和球形气泡之间的次Bjerknes力和两个球形气泡的次Bjerknes力方向差异较大,在大小上是两个球形气泡次Bjerkens力的数倍,且有着更长的作用距离.非球形气泡和球形气泡之间的次Bjerknes力取决于非球形气泡的形状模态、两个气泡初始半径的比值、驱动声压振幅、气泡间距和两个气泡的相对位置.     

两模腔场谱间的量子干涉

研究了高Q腔中单个二能级原子与两模二项式光场依赖强度耦合相互作用系统的腔场谱,给出了弱初始场条件下的数值结果,讨论了两模光场之间的量子干涉对腔场谱结构的影响. 发现当两模光场的频率差Δ>g（g为原子与腔场间的耦合常数）时,两模光场间的干涉效应对谱结构没有影响,系统的腔肠谱只是两模腔肠谱的简单叠加;当Δ≤g时两模腔场谱间的干涉比较明显. 在强初始场条件下,量子干涉效应可忽略. 关键词： 腔场谱 量子干涉 两模二项式光场     

外力诱导吸附高分子单链的拉伸分子动力学研究

采用拉伸分子动力学方法研究聚乙烯单链在真空中从一个强吸附表面被拉伸的动力学过程. 当聚乙烯单链从中间1/3处被拉伸时,通过理想弹簧可测量其平均力〈f〉,发现在拉力与拉伸距离的曲线上出现了两个平台,这与实验结果非常一致,同时研究了在这个过程中链的形状变化情况,发现形状变化经历了四个阶段. 当聚乙烯单链的一端与吸附表面成一定倾斜角拉伸时,发现通过理想弹簧测量到的平均吸附力〈f〉与倾斜角θ在θ>20°的范围里总保持线性关系,并且其斜率随着拉伸速度的增加而增 关键词： 拉伸分子动力学 聚乙烯链 吸附现象     

波浪破碎气体的卷入过程及相关统计量的估计

基于实验观测,导出了波浪破碎能量耗散ε_ed、气泡云卷入深度z_b、气体卷入速率Q(z)和湍流动能耗散率ε_T(z)的表达式,在此基础上建立了一种简单、实用的气泡粒径谱参数化模式N(a,z),揭示了波浪破碎气泡云卷入过程能量耗散、气泡破碎临界Hinze特征尺度和气泡粒径谱在不同海况下的变化. 研究表明：气泡云卷入过程能量耗 关键词： 波浪破碎能量耗散 气泡云卷入深度 气泡粒径谱     

欧拉圆盘不同能量耗散机理之间的关联

本文研究了欧拉圆盘运动过程中盘厚度以及盘面与水平面夹角α两因素对能量耗散的影响. 得出圆盘厚度与直径之比x对能量变化中各项因子的影响: x很小时, 质心在竖直方向上的动能变化和重力势能变化是系统能量耗散的主要因素; 当x>0.4142时, 圆盘绕与之平行的轴的转动动能变化成为主要因素, 并给出圆盘厚度可忽略的条件. 模拟了滚动摩擦、空气黏滞等不同能量耗散方式与x,α的关系, 导出各种耗散方式在圆盘运动的过程中的转变规律, 并指出x=0.1733, α>18°时能量耗散形式为纯滚动摩擦, 这修正了文献[26]结论.     

在有限温度下运用弱测量保护量子纠缠

根据单量子比特的映射, 提出一个在有限温度下运用弱测量保护两量子比特纠缠的方案. 在两个可通过局域幺正变化实现的等价的初始纠缠进入广义振幅阻尼信道前后, 对其分别进行弱测量, 并对四个弱测量的参数做全面的优化, 获取最大共生纠缠C_r 和弱测量参数m, n的解析表达式, 然后再进一步研究弱测量参数与信道参数的关系. 发现这种基于弱测量的纠缠保护方法在某些情况下可以有效地提高纠缠, 甚至可以避免纠缠的突然死亡. 当信道参数r一定时, 对不同的参数p,初始态|ψ>纠缠达到最大值时对应的弱测量参数m的取值一样, 共生纠缠关于p=0.5对称, 而初始态|φ>对应的参数m的取值不同; 在参数p一定、参数r不同、初始态|ψ>或|φ>的情况下, 当纠缠度取最大值时, 弱测量参数m的取值不变, 且随着r的增加, 纠缠度减少. 通过对信道参数的分析, 发现可以选择合适的信道参数和初始态来获得较大的纠缠.     

Effects of short-range attraction on Jamming transition

Enormous progresses to understand the jamming transition have been driven via simulating purely repulsive particles which were somehow idealized in the past two decades. While the attractive systems are both theoretical and practical compared with repulsive systems. By studying the statistics of rigid clusters, we find that the critical packing fraction φ_c varies linearly with attraction μ for different system sizes when the range of attraction is short. While for systems with long-range attractions, however, the slope of φ_c appears significantly different, which means that there are two distinct jamming scenarios. In this paper, we focus our main attention on short-range attractions scenario and define a new quantity named "short-range attraction susceptibility" χ_p, which describes the degree of response of the probability of finding jammed states pjto short-range attraction strength μ. Our central results are that χ_p diverges in the thermodynamic limit as χ_p ∝|φ-φ_c~∞|~(-γ_p), where φ_c~∞ is the packing fraction at the jamming transition for the infinite system in the absence of attraction. χ_p obeys scaling collapse with a scaling function in both two and three dimensions, illuminating that the jamming transition can be considered as a phase transition as proposed in previous work.     

次Bjerknes力作用下气泡的体积振动和散射声场

利用Lagrange方程得到了次Bjerknes力作用下气泡的体积振动方程,并探讨了次Bjerknes力作用下不同参数对气泡体积振动振幅和振动初相位的影响,研究了振动初相位差为π和0的气泡对在液体中形成的散射声场特征.结果表明:次Bjerknes作用力下,相邻气泡半径、气泡间距、多方指数均能影响气泡的体积振动振幅,气泡对的均衡半径、气泡间距和驱动频率则对气泡振动初相位产生明显影响;相距很近、相位相差为π的两个气泡的散射声压与气泡体积振动振幅、气泡间距、驱动频率和振动初相位有关,随声场距离成反比减小,与声场位置有关,其平均散射声功率是单个孤立气泡的1/6(kd_(12))~2半径相同、相距很近、相位相同的两个6气泡的散射声压与气泡振动初相位、体积振动振幅、气泡间距、驱动频率有关,随声场距离成反比减小,其平均散射声功率是单个孤立气泡的4倍.     

耦合双泡声空化特性的理论研究

当双泡中心间距足够小时,由于气泡间辐射压力波的存在,作用在气泡上的压力不等于外部驱动压力.通过考虑双泡之间的辐射压力波,利用改进的Keller-Miksis方程,分别计算了不同大小、不同间距、含不同惰性气体的双泡在声空化过程中半径的变化、次Bjerknes力的变化和双泡内温度的变化.计算结果表明,当双泡大小不同时,小气泡受到的抑制作用较强,温度变化也比较大.随着双泡间距离从100μm增大到1 cm时,气泡间的次Bjerknes力的数量级从10~(-4)N减小到10~(-8)N.含不同惰性气体的耦合双泡在回弹阶段表现出明显不同的振荡规律.     

窄线宽脉冲光纤激光的自相位调制预补偿研究

自相位调制(SPM)效应会展宽窄线宽脉冲光纤激光的光谱宽度,降低其相干性.通过相位调制对SPM引起的非线性相移进行预补偿,能够使脉冲激光在光纤中进行放大和传输后保持种子激光的光谱特性.基于三波耦合方程开展数值仿真,研究了在对SPM进行"欠补偿","完全补偿"和"过补偿"的情况下,SPM预补偿对受激布里渊散射阈值和激光光谱特性的影响.开展了SPM预补偿实验研究,将脉冲激光的光谱宽度从1.4 GHz压缩到120 MHz.研究内容可以为窄线宽脉冲光纤激光系统的设计搭建提供参考.     

In-situ TEM observation of the evolution of helium bubbles in Mo during He+ irradiation and post-irradiation annealing

The evolution of helium bubbles in purity Mo was investigated by in-situ transmission electron microscopy (TEM) during 30 keV He⁺ irradiation (at 673 K and 1173 K) and post-irradiation annealing (after 30 keV He⁺ irradiation with the fluence of 5.74×10¹⁶ He⁺/cm² at 673 K). Both He⁺ irradiation and subsequently annealing induced the initiation, aggregation, and growth of helium bubbles. Temperature had a significant effect on the initiation and evolution of helium bubbles. The higher the irradiation temperature was, the larger the bubble size at the same irradiation fluence would be. At 1173 K irradiation, helium bubbles nucleated and grew preferentially at grain boundaries and showed super large size, which would induce the formation of microcracks. At the same time, the geometry of helium bubbles changed from sphericity to polyhedron. The polyhedral bubbles preferred to grow in the shape bounded by {100} planes. After statistical analysis of the characteristic parameters of helium bubbles, the functions between the average size, number density of helium bubbles, swelling rate and irradiation damage were obtained. Meanwhile, an empirical formula for calculating the size of helium bubbles during the annealing was also provided.     

II. Recirculation flow induced by a bubble stream rising in a viscous liquid

The recirculation flow induced by the rising motion of a bubble stream in a viscous fluid within an open-top rectangular enclosure is studied. The three-dimensional volume averaged conservation equations are solved by a control-volume method using a hybrid finite differencing scheme to describe the liquid phase hydrodynamics. The momentum exhange between the bubbles and the liquid phase is modeled with a source term equals to the volumetric buoyancy force acting on the gas in the bubble stream. The volumetric buoyancy force accounts for in line interactions between bubbles through the average gas volume fraction in the gas liquid column which depends on the size and the rising velocity of bubbles. The fluid flow within an open-top rectangular enclosure is further investigated by particle image velocimetry for a bubble stream rising in a water-glycerol solution. The measured fluid velocities in a vertical plane are compared with the predictions of the numerical model over a wide range of fluid viscosity (43 mPa s-800 mPa s) and gas flow rates. Finally, the recirculation flows resulting from the interaction of two neighbouring vertical bubble streams are studied. Received: 23 July 1997 / Revised: 19 December 1997 / Accepted: 11 May 1998     

Ge-Sb-Se硫系玻璃的折射率和热光系数

Ge-Sb-Se硫系玻璃被认为是极佳的红外传输材料和有潜力的非线性光学材料.在光学设计中,玻璃的线性折射率(n)及其热光系数(ζ)是关键技术参数.以预测和调控Ge-Sb-Se玻璃的n和ζ为目的,考察了玻璃的n,ζ,密度(d)和体积膨胀系数(β)与化学参数dSe和拓扑网络结构参数r的内在联系.研究发现,玻璃的n随d的增加而增大;ζ随β的增大而近似线性减小;β随dSe的减小或r的增大而减小;当Ge含量固定时, d随dSe的减小或r的增大而增大,当Sb含量固定时, d在dSe=0时具有最小值.基于实测d和n,拟合获得了Ge, Sb和Se元素在2—12μm波段的摩尔折射度(R_i),分别为R_(Ge)=10.16—10.50 cm~3/mol,RSb=16.71—17.08 cm~3/mol和RSe=11.15—11.21 cm~3/mol,根据d和R_i计算得到的n与实测值的偏差小于1%.基于实测ζ和β,拟合得到了Ge, Sb和Se元素在2—12μm波段的摩尔折射度温度系数(φ_i),分别为φ_(Ge)=21.1—22.6 ppm/K,φ_(Sb)=7.2—8.4 ppm/K和φ_(Se)=90.2—94.2 ppm/K,根据β和φ_i计算得到的ζ与实测值的偏差小于6 ppm/K.     

Microstreaming velocity field and shear stress created by an oscillating encapsulated microbubble near a cell membrane

Sonoporation mediated by microbubbles is being extensively studied as a promising technology to facilitate gene/drug delivery to cells. However, the theoretical study regarding the mechanisms involved in sonoporation is still in its infancy. Microstreaming generated by pulsating microbubble near the cell membrane is regarded as one of the most important mechanisms in the sonoporation process. Here, based on an encapsulated microbubble dynamic model with considering nonlinear rheological effects of both shell elasticity and viscosity, the microstreaming velocity field and shear stress generated by an oscillating microbubble near the cell membrane are theoretically simulated. Some factors that might affect the behaviors of microstreaming are thoroughly investigated, including the distance between the bubble center and cell membrane (d), shell elasticity (χ), and shell viscosity (κ). The results show that (i) the presence of cell membrane can result in asymmetric microstreaming velocity field, while the constrained effect of the membrane wall decays with increasing the bubble-cell distance; (ii) the bubble resonance frequency increases with the increase in d and χ, and the decrease in κ, although it is more dominated by the variation of shell elasticity; and (iii) the maximal microstreaming shear stress on the cell membrane increases rapidly with reducing the d, χ, and κ. The results suggest that microbubbles with softer and less viscous shell materials might be preferred to achieve more efficient sonoporation outcomes, and it is better to have bubbles located in the immediate vicinity of the cell membrane.     

The effects of coupling and bubble size on the dynamical-systems behaviour of a small cluster of microbubbles

This study endeavours to apply a theoretical model for predicting the dynamics of a bubble cluster of various sizes, within which each bubble may assume different initial conditions from other bubbles in the cluster. The resulting system of coupled Keller-Miksis-Parlitz equations are solved numerically, and the effects of coupling and bubble size on bubble cluster dynamics are examined for a given set of ultrasound parameters. It has been found that the effects of coupling are significant, and a bubble cluster's bifurcation characteristics and route to chaos can be altered by inter-bubble interactions. This gives rise to the possibility of suppressing the chaotic oscillations of microbubbles by varying bubble cluster size. Small equilibrium radii bubbles have little influence on the dynamics of neighbouring bubbles in a cluster via coupling. Furthermore, a bubble system consisting of smaller-sized bubbles transitions from order to chaos at lower driving pressure amplitudes.     

Anisotropic elastic properties of microtubules

The present study aims at providing insight into the acceleration mechanism of a bubble chain rising in shear-thinning viscoelastic fluids. The experimental investigation by the Particle Image Velocimetry (PIV), birefringence visualisation and rheological simulation shows that two aspects are central to bubble interactions in such media: the stress creation by the passage of bubbles, and their relaxation due to the fluids memory forming an evanescent corridor of reduced viscosity. Interactions between bubbles were taken into account mainly through a linear superposition of the stress evolution behind each bubble. An analytical approach together with the rheological consideration was developed to compute the rise velocity of a bubble chain in function of the injection period and bubble volume. The model predictions compare satisfactorily with the experimental investigation.