三维方腔介电液体电对流的数值模拟研究

本文对封闭方腔内介电液体电对流进行了三维数值模拟研究.方腔的6个边界为固壁;4个侧边界为电绝缘边界;上下界面为两个电极.直流电场作用在从底部电极注入的自由电荷上,从而对液体施加库伦体积力并驱动流体流动形成电对流.为了求解这一物理问题,发展了一种二阶精度的有限体积法来求解完整的控制方程,包括Navier-Stokes方程和一组简化的Maxwell方程.考虑到电荷密度方程的强对流占优特性,采用了全逆差递减格式来求解该方程,获得了准确有界的解.通过研究发现,该流动在有限振幅区内的分叉类型为亚临界,即系统存在一个线性和非线性临界值,分别对应流动的开始和终止.由于非线性临界值比线性值小,因此两个临界值之间有一个迟滞回线.与无限大域中的自由对流相比,侧壁施加的额外约束改变了流场结构,使这两个临界值均有所增大.此外,还讨论了电荷密度和速度场的空间分布特征,发现电荷密度分布中存在电荷空白区.最后对更小空间尺寸情况计算结果表明,流动的线性分叉类型为超临界.本文的结果拓展了已有的二维有限空间内电对流的研究,并为三维电对流的线性和弱非线性理论分析提供参考.     

离子选择性表面电对流的多块LB模拟

离子选择性表面(如纳米通道、离子交换膜等)复杂的动力学现象为微纳流控技术的发展提供了新思路. 向带有离子选择性表面的电解质溶液施加电压, 通过液体的电流密度会经历复杂的非线性变化过程; 当电压超过某一临界值时会引发对流现象, 这种流动被称为第二类电渗或离子选择性表面的电对流, 关于此类问题的数值求解引发了大量的研究. 本文提出一种基于多块网格加密的格子玻尔兹曼方法(lattice Boltzmann method, LBM)的数值模型, 用于模拟第二类电渗流动. 结合该算法, 给出了求解流动、电势和离子浓度的网格信息交换方程, 较好地解决了此类问题中大浓度梯度边界对计算分辨率的要求. 利用该数值模型模拟获得的电流?电压特性曲线先随着电压升高而迅速增大, 随后达到饱和状态, 与理论解吻合良好. 此外, 模拟结果还表明, 当流动发生后, 相对低电压下的流动倾向于形成大涡且流动呈指数趋势增强; 而较大电压则会先激发多个小涡, 并逐渐合并为大涡流动, 且大涡流动有更高的离子输运效率. 此外, 除了模拟离子选择性表面的电对流现象, 本文提出的数值格式还可拓展到其他电流体动力学问题的模拟.      

电流变液体黏性变化趋势的试验研究

以微观试验和流变性能试验为手段，分别研究零电场下和在电场作用下的电流变液体黏性变化规律．研究结果表明：零电场下电流变液体的黏性与Krieger-Dougherty公式具有很好的拟合效果，其中逾渗临界值强依赖于悬浮液体中固体颗粒的性质并随工作温度变化．在电场作用下，电流变悬浮液体的黏度随剪切速率的变化规律分为3个阶段：即呈线性的启动段、非线性的幂定律模型流动段和宾汉模型流动段．研究结果为电流变效应工程应用提供依据．     

进动圆筒内液体流动不稳定的非线性特征

在建立进动充液圆筒内液体偏差流动方程的基础上，结合液体惯性波和轴向二次流动线性解，通过对定常二次流动的线性稳定性分析，提出了函数空间表达的流动不稳定性非线性分岔分析方程. 对非惯性坐标系下液体流动的Navier-Stokes方程进行了数值求解，并对惯性波发生破裂(实验提供的3种主模态下得出的共振破裂现象)时的压力时间序列进行分析，得出了液体流动不稳定的基本非线性特征.     

边界拟合坐标系下的差分有限元破开算子法

在边界拟合坐标系下给出耦合了有限差分法和有限单元法的新型的破开算子法.利用该方法,Navier-Stokes方程被分解成对流方程和扩散方程,对流方程采用稳定性好的有限差分法求解,而扩散方程则采用有限单元法求解.由于计算是基于非均匀网格,采用边界加密的51×51网格就达到了前人在计算雷诺数为5000的方腔环流时采用的257×257均匀网格的效果,对于瑞利数Rt=107的竖腔自然对流的计算进一步表明,提出的方法是有效的.对于Re=105的高雷诺数方腔受驱环流,得到了非稳态、非周期的和带有随机特征的流场结构.     

求解对流扩散方程的一种高效的有限体积法

考虑无结构三角网格上求解对流扩散方程的有限体积法.引入一种梯度函数的计算方法,将现有方法中计算解变量在网格单元中心和网格单元边界的梯度的两个独立过程改造成一个过程来完成,发展了一种求解对流扩散方程的高效的有限体积法.数值实验结果表明,该方法完全达到了已有方法同样的精度,而在计算速度上有明显的提高.     

磁-电-弹性半空间在轴对称热载荷作用下的三维问题研究

考虑力-电-磁-热等多场耦合作用, 基于线性理论给出了磁-电-弹性半空间在表面轴对称温度载荷作用下的热-磁-电-弹性分析, 并得到了问题的解析解. 利用Hankel 积分变换法求解了磁-电-弹性材料中的热传导及控制方程, 讨论了在磁-电-弹性半空间在边界表面上作用局部热载荷时的混合边值问题, 利用积分变换和积分方程技术, 通过在边界表面上施加应力自由及磁-电开路条件, 推导得到了磁-电-弹性半空间中位移、电势及磁势的积分形式的表达式. 获得了磁-电-弹性半空间中温度场的解析表达式并且给出了应力, 电位移和磁通量的解析解. 数值计算结果表明温度载荷对磁-电-弹性场的分布有显著影响. 当温度载荷作用的圆域半径增大时, 最大正应力发生的位置会远离半无限大体的边界; 反之当温度载荷作用的圆域半径减小时, 最大应力发生的位置会靠近半无限大体的边界. 电场和磁场在温度载荷作用的圆域内在边界表面附近有明显的强化, 而磁-电-弹性场强化区域的强化程度跟温度载荷的大小和作用区域大小相关. 本研究的相关结果对智能材料和结构在热载荷作用下的设计和制造具有指导意义.      

一种新型介电弹性体仿生可调焦透镜的变焦分析

介电弹性体 (dielectric elastomer) 是电活性聚合物智能材料的一种,在外加电场作用下,可产生多种形式的响应.在驱动柔性透镜的变焦方面,相对于传统的机械操控变焦方法 显示出独特的优势.针对一款在电压激励下可高效调节焦距的介电弹性体仿人眼变焦透镜,该透镜由上下两层介电弹性薄膜和固定框架构成,并在封闭腔内充入盐水,上层薄膜涂覆环形柔性电极.在电压激励下,上层膜发生变形,由于盐水的体积保持恒定,引 起下层膜随之变形,使得透镜的焦距发生改变.采用 neo-Hookean 模型,利用变分原理导出了该透镜的控制方程、边界条件和连 续条件.利用打靶法求解了该非线性问题并高效地处理了非线性问题的界面连续条件. 理论分析结果与实验结果相吻合. 利用此模型开展了广泛的参数分析,研究表明,透镜的几何形状、初始焦距、介电弹性体薄膜的预拉伸率、涂覆的电极面积、材料的剪切模量等对透镜焦距的调节性能都有重要的影响.所建立的理论分析模型可为柔性仿生透镜的设计和参数优化提供有效的分析方法.      

直接边界元法三维自由面兴波时域模拟

应用直接边界元法在时域中求解稳定航速运动的三维自由面兴波问题.基于格林定理,在所有边界面上划分网格,对边界积分方程进行数值离散,采用线性自由面边界条件,随时间步进更新自由面势.由于物体空间位置移动辐射条件不需要单独表述,迭代过程中自由面计算域保持不变.以割划水面NACA0024为例,计算模拟了自由面兴波稳定波形;提出了求解矩阵方程组奇异性的处理方法和解决割划问题的动网格技术.本文计算结果和有限体积法及有关试验结果对比表明,该方法是可靠的.     

突扩管流中微米颗粒分离效率的研究

助欧拉和拉格朗日方法数值模拟了突扩微尺度管道流中微米颗粒的分离情况. 在采用有限体积法求解电荷密度场、电场和流场的基础上,通过基于改进的Langevin方程研究了微管道中不同位置处的微米颗粒在水动力和介电电泳力综合作用下的运动轨迹. 研究发现:电渗流的驱动能力随着扩张比(ER)的增加而提高,然而其提高的趋势逐渐变小;当微米颗粒仅在水动力作用下时,随着ER的增加,颗粒之间的有效分离距离(ESL)随之线性增加,此时ESL与ER的比值约为5.9; 若是考虑介电电泳力对于微米颗粒运动的影响, ESL与ER的比值下降为4.79, 由此可以看出介电电泳力对突扩微管道流中的微米颗粒的分离效果有着一定的负面影响. 上述结论对于基于介电电泳技术设计的生物粒子分离芯片的优化设计有很大价值.      

多孔介质内黏弹性流体的热对流稳定性研究

基于修正的Darcy模型, 介绍了多孔介质内黏弹性流体热对流稳定性研究的现状和主要进展. 通过线性稳定性理论, 分析计算多孔介质几何形状(水平多孔介质层、多孔圆柱以及多孔方腔)、热边界条件(底部等温加热、底部等热流加热、底部对流换热以及顶部自由开口边界)、黏弹性流体的流动模型(Darcy-Jeffrey, Darcy-Brinkman-Oldroyd以及Darcy-Brinkman -Maxwell模型)、局部热非平衡效应以及旋转效应对黏弹性流体热对流失稳的临界Rayleigh数的影响. 利用弱非线性分析方法, 揭示失稳临界点附近热对流流动的分叉情况, 以及失稳临界点附近黏弹性流体换热Nusselt数的解析表达式. 采用数值模拟方法, 研究高Rayleigh数下黏弹性流体换热Nusselt数和流场的演化规律,分析各参数对黏弹性流体热对流失稳和对流换热速率的影响.主要结果:(1)流体的黏弹性能够促进振荡对流的发生;(2)旋转效应、流体与多孔介质间的传热能够抑制黏弹性流体的热对流失稳;(3)在临界Rayleigh数附近,静态对流分叉解是超临界稳定的, 而振荡对流分叉可能是超临界或者亚临界的,主要取决于流体的黏弹性参数、Prandtl数以及Darcy数;(4)随着Rayleigh数的增加,热对流的流场从单个涡胞逐渐演化为多个不规则单元涡胞, 最后发展为混沌状态.      

二维矩形腔体中混合流体行进波对流

Rayleigh-Benard模型是研究对流稳定性,时空结构和非线性特性的典型模型之一。本文的兴趣集中在二维矩形腔体中混合流体对流场的结构方面。利用SIMPLE算法数值求解流体力学方程组,模拟了充分发展的二维矩形腔体中混合流体对流。结果说明偏离传导失去稳定的系统经过亚临界分叉产生了振动对流。进一步,我们给出了分叉曲线及其沿分叉曲线的上部分支三个Rayleigh数对应的对流图案的垂直速度场,流线,温度场,浓度场和Shadowgraph强度的等值线图。所有场的结构分析表明浓度场及Shadowgraph强度的等值线图可以很好的特征行进波的运动特性。     

轮对系统的Hopf分岔研究

针对轮对系统中的非线性动力学问题, 本文基于Hopf分岔代数判据得到考虑陀螺效应的轮对系统Hopf分岔点解析表达式, 即轮对系统蛇形失稳的线性临界速度解析表达式. 基于分岔理论得到轮对系统的第一、第二Lyapunov系数表达式, 并结合打靶法分别得到不同纵向刚度下, 考虑陀螺效应与不考虑陀螺效应的轮对系统分岔图. 通过对比有无陀螺效应的轮对系统分岔图发现, 在同一纵向刚度下, 考虑陀螺效应的轮对系统线性临界速度和非线性临界速度均大于不考虑陀螺效应的轮对系统, 即陀螺效应可以提高轮对系统的运动稳定性. 基于Bautin分岔理论, 以纵向刚度和纵向速度作为参数, 分别得到考虑陀螺效应和不考虑陀螺效应的轮对系统, 从亚临界Hopf分岔到超临界Hopf分岔, 再从超临界Hopf分岔到亚临界Hopf分岔的迁移机理拓扑图. 通过对比有、无陀螺效应的轮对系统Bautin分岔拓扑图发现, 陀螺效应将改变轮对系统的退化Hopf分岔点, 但对于轮对系统Bautin分岔拓扑图的影响不大.      

Nonlinear dynamics of hardware-in-the-loop experiments on stick–slip phenomena

A single degree-of-freedom nonlinear mechanical model of the stick–slip phenomenon is studied when the Stribeck-type friction force is emulated by means of a digitally controlled actuator. The relative velocity of the slipping contact surfaces is considered as bifurcation parameter. The original physical system presents subcritical Hopf bifurcation with a wide bistable parameter region where stick–slip and steady-state slipping are both stable locally. Hardware-in-the-loop experiments are performed with a physical oscillatory system subjected to the emulated Stribeck forces. The effect of sampling time is studied with respect to the stability and nonlinear behavior of this experimental system. The existence of subcritical Neimark–Sacker bifurcations are proven in the digital system, the stability and bifurcation characteristics of the continuous and the digital systems are compared, and the counter-intuitive stabilizing effect of sampling time is shown both analytically and experimentally. The conclusions draw the attention to the limitations of hardware-in-the-loop experiments when the corresponding systems are strongly nonlinear.     

Linear and weakly nonlinear variable viscosity convection in spherical shells

A theoretical study of linear and weakly nonlinear thermal convection in a spherical shell is performed. The Boussinesq fluid is of infinite Prandtl number and its viscosity is temperature dependent. The linear stability eigenvalue problem is derived and solved by a shooting method assuming isothermal, stress-free boundaries, a self-gravitating fluid, and corresponding to two heating models. The first is heating from below, and the second is a model of combined heating from below and within, such that convection is described by a self-adjoint linear stability formulation. In addition, nonlinear, hemispherical, axisymmetric convection is computed by a finite volume technique for a shell with 0.5 aspect ratio. It is shown that 2-cell convection occurs as transcritical bifurcation for a viscosity constrast across the shell up to about 150. Motions with four cells are also possible. As expected, the subcritical range is found to increase with increasing viscosity contrast, even when the linear operator is self-adjoint.This research was supported by the AT&T Foundation.     

Convection due to the selective absorption of radiation in a porous medium

A model for convection due to the selective absorption of radiation in a fluid saturated porous medium is investigated. The model is based on a similar one introduced for a viscous fluid by Krishnamurti [x]. Employing this adapted model we show the growth rate for the linearised system is real. A linear instability analysis is performed. Global stability thresholds are also found using nonlinear energy theory. An excellent agreement is found between the linear instability and nonlinear stability Rayleigh numbers, so that the region of potential subcritical instabilities is very small, demonstrating that the linear theory accurately emulates the physics of the onset of convection. Received February 10, 2003 / Accepted February 10, 2003/ Published online May 9, 2003 / B. Straughan     

Nonlinear dynamics of extensible fluid-conveying pipes,supported at both ends

In this paper, the post-divergence behaviour of extensible fluid-conveying pipes supported at both ends is studied using the weakly nonlinear equations of motion of Semler, Li and Païdoussis. The two coupled nonlinear partial differential equations are discretized via Galerkin's method and the resulting set of ordinary differential equations is solved either by Houbolt's finite difference method or via AUTO. Typically, the pipe is stable at its original static equilibrium position up to the flow velocity where it loses stability by static divergence via a supercritical pitchfork bifurcation. The amplitude of the resultant buckling increases with increasing flow, but no secondary instability occurs beyond the pitchfork bifurcation. The effects of the system parameters on pipe behaviour as well as the possibility of a subcritical pitchfork bifurcation have also been studied.     

Convection due to the selective absorption of radiation in a porous medium

A model for convection due to the selective absorption of radiation in a fluid saturated porous medium is investigated. The model is based on a similar one introduced for a viscous fluid by Krishnamurti [x]. Employing this adapted model we show the growth rate for the linearised system is real. A linear instability analysis is performed. Global stability thresholds are also found using nonlinear energy theory. An excellent agreement is found between the linear instability and nonlinear stability Rayleigh numbers, so that the region of potential subcritical instabilities is very small, demonstrating that the linear theory accurately emulates the physics of the onset of convectionReceived: 10 February 2003, Accepted: 11 March 2003, Published online: 12 September 2003     

Weakly nonlinear stability analysis of triple diffusive convection in a Maxwell fluid saturated porous layer

The weakly nonlinear stability of the triple diffusive convection in a Maxwell fluid saturated porous layer is investigated. In some cases, disconnected oscillatory neutral curves are found to exist, indicating that three critical thermal Darcy-Rayleigh numbers are required to specify the linear instability criteria. However, another distinguishing feature predicted from that of Newtonian fluids is the impossibility of quasi-periodic bifurcation from the rest state. Besides, the co-dimensional two bifurcation points are located in the Darcy-Prandtl number and the stress relaxation parameter plane. It is observed that the value of the stress relaxation parameter defining the crossover between stationary and oscillatory bifurcations decreases when the Darcy-Prandtl number increases. A cubic Landau equation is derived based on the weakly nonlinear stability analysis. It is found that the bifurcating oscillatory solution is either supercritical or subcritical, depending on the choice of the physical parameters. Heat and mass transfers are estimated in terms of time and area-averaged Nusselt numbers.     

Nonlinear dynamics of functionally graded pipes conveying hot fluid

For improved stability of fluid-conveying pipes operating under the thermal environment, functionally graded materials (FGMs) are recommended in a few recent studies. Besides this advantage, the nonlinear dynamics of fluid-conveying FG pipes is an important concern for their engineering applications. The present study is carried out in this direction, where the nonlinear dynamics of a vertical FG pipe conveying hot fluid is studied thoroughly. The FG pipe is considered with pinned ends while the internal hot fluid flows with the steady or pulsatile flow velocity. Based on the Euler–Bernoulli beam theory and the plug-flow model, the nonlinear governing equation of motion of the fluid-conveying FG pipe is derived in the form of the nonlinear integro-partial-differential equation that is subsequently reduced as the nonlinear temporal differential equation using Galerkin method. The solutions in the time or frequency domain are obtained by implementing the adaptive Runge–Kutta method or harmonic balance method. First, the divergence characteristics of the FG pipe are investigated and it is found that buckling of the FG pipe arises mainly because of temperature of the internal fluid. Next, the dynamic characteristics of the FG pipe corresponding to its pre- and post-buckled equilibrium states are studied. In the pre-buckled equilibrium state, higher-order parametric resonances are observed in addition to the principal primary and secondary parametric resonances, and thus the usual shape of the parametric instability region deviates. However, in the post-buckled equilibrium state of the FG pipe, its chaotic oscillations may arise through the intermittent transition route, cyclic-fold bifurcation, period-doubling bifurcation and subcritical bifurcation. The overall study reveals complex dynamics of the FG pipe with respect to some system parameters like temperature of fluid, material properties of FGM and fluid flow velocity.