简支Mexwell模型粘弹性输流管道的稳定性分析

在弹性输流管道研究的基础上,采用递推格式的有限差分法,对简支Maxwell模型粘弹性输流管道（回转守恒系统）,探讨了其动力特性和稳定性问题,具体分析了材料的松弛时间对无量纲流速与前三阶模态的无量纲频率的实部及虚部之间的变化曲线的影响。发现发散临界流速随松弛时间的减小而降低,随后发生的耦合模态颤振临界流速随松弛时间的减小而增大;甚至在质量比较大时,随着松弛时间的减小,可推迟乃至不发生耦合模态颤振。当无量纲松弛时间达到10^3量级以上时,即可将其按弹性管道处理。甚至在H为10^2量级时,按弹性管道处理也不会带来太大的误差。     

具有水平流动的对流斑图成长和动力学特性

采用二维流体力学基本方程组对Prandtl数Pr=0.72具有水平流动的对流斑图成长和动力学特性进行了数值模拟.结果说明,对于给定的相对Rayleigh数Ra_r=5(Rayleigh数Ra=8 540)和Reynolds数Re=22.5,行波对流斑图的成长分为三个阶段,即对流发展阶段、指数成长阶段、周期变化阶段(过渡调整区、稳定周期变化区).行波对流的平均波数随着时间的发展或者对流斑图的成长而减小.随着相对Rayleigh数的增加,行波对流的指数成长阶段的时间变短,对流最大垂直流速的成长率变大.对于水平流动Re=5时,对流最大垂直流速的成长率γ_m与Ra_r的关系为γ_m=0.004 8Ra~(6.065 3)_r.在周期变化阶段,经过行波对流斑图和对流参数的过渡调整区后,对流进入斑图和对流参数的稳定周期变化区.对于给定的Ra_r=5时,行波对流的无量纲周期T_t随着Re变化的关系式为T_t=0.001 4Re~(2.363 5).     

竖直平板嵌入非Darcy多孔介质时,数值研究流过平板的不稳定MHD自然对流中的Dufour和Soret效应

就竖直平板嵌入非Darcy多孔介质中,导电流体流过平板时作不稳定的二维磁流体(MHD)双扩散对流,数值研究了Dufour和Soret效应对流动的影响.用Crank-Nicolson型的隐式有限差分法,按三对角矩阵处理,求解无量纲的非线性控制方程.详细地研究了问题中出现的各种参数对不稳定无量纲的速度、温度和浓度曲线的影响.进一步地,给出并分析了表面摩擦因数、Nus-selt数和Sherwood数随时间的变化.研究结果表明,不稳定速度、温度和浓度分布曲线,受Dufour和Soret的影响十分显著.随着Dufour数的增加或者Soret数的减小,表面摩擦因数和Sherwood数都在减小,而Nusselt数在增加.研究发现,当磁场参数增加时,边界层中的速度和温度在减小.     

水平来流对扰动成长和对流周期性的影响

对Pr=0.0272的纯流体在矩形腔体外加水平来流时,进行二维流体力学基本方程组的数值模拟.研究了该纯流体Rayleigh-Benard对流的一维行波斑图的成长及时空的演化.发现对流成长过程可以划分为3个阶段,即对流发展、对流指数成长和周期变化。在对流指数成长阶段对不同相对Rayleigh(瑞利)数Ra_r的最大垂直流速场随时间变化的情况进行分析,获得了最大垂直流速场指数成长阶段的线性成长率γ_m和相对Rayleigh数Ra_r的关系公式.研究了行波周期受水平来流Reynolds(雷诺)数的影响,揭示了行波对流周期性及其对水平来流Reynolds数的依赖性.     

一种新的混合流体对流竖向镜面对称对传波斑图

利用SIMPLE算法对混合流体对流的流体力学基本方程组进行了数值模拟,在混合流体分离比ψ=－0.6和矩形腔体长高比Γ=20的情况下,首次发现了一种新的竖向镜面对称对传波斑图,并初步探讨了它的动力学特性.竖向镜面对称对传波斑图的中心为驻波,随着时间的发展驻波的波长伸长.当波长增加到某个临界值时,一个滚动分裂成两个滚动,在这两个滚动之间产生一个具有180°相位差的新滚动.位于中心线上的滚动只有相位的突变及其波长的压缩或者伸长,没有对流滚动的移动,在它的两侧是向左右传播的对流滚动.驻波两次相位突变形成一个周期,驻波周期随着相对Rayleigh(瑞利)数Ra_r的增加而增加.这种对流结构存在于相对Rayleigh数Ra_r∈(3.6,4.3]的范围.当相对Rayleigh数Ra_r≤3.6时,系统出现具有缺陷的行波斑图;当Ra_r>4.3时系统过渡到行波斑图.说明竖向镜面对称对传波斑图是存在于具有缺陷的行波斑图和行波斑图之间的一种稳定的对流斑图.     

扁壳压曲的蠕变效应

本文讨论了混凝土扁壳压曲的蠕变效应.基于弹性薄壳的非线性理论,发现椭圆抛物面扁壳,其荷载-挠度曲线的上临界荷载将随时间而降低,而下临界荷载则随时间而上升.至于扁壳的局部失稳问题,其临界荷载仅取决于压曲发生瞬间材料的弹性模量.     

轴承-转子系统分岔参数的分维数识别法

广泛用于各行各业的转子系统的稳定性问题一直倍受关注.而在当今失稳多是由于一些非线性现象的出现所引起,这就对转子系统的设计提出了更高的要求:考虑非线性因素,避开会出现非线性现象的不稳定参数点或区域.若仅知未知系统的系列时间序列(有可能被噪声污染),如何识别系统运动性态的变化?为了探讨此问题,在本文中通过对一单盘Jeffcott转子的研究,得出了利用随参数变化的时间序列分维数趋势图,可以很好地识别轴承-转子动力系统发生分岔时的临界参数.     

Marangoni效应对等曲率弯曲界面凝固过程的影响

用渐近方法和相场数值模拟分别研究了球状晶体和柱状晶体凝固过程中,固液界面表面张力随温度的变化对界面运动的影响．结果表明Marangoni效应将增大临界晶核半径,减缓界面运动速度．在静止熔体中,枝晶尖端生长速度随Marangoni数线性下降．渐近展开和相场模拟得到的结果定性是一致的．     

扁球壳在热-机械荷载作用下的稳定性分析

基于扁壳几何非线性理论,应用虚功原理和变分法推导了均匀变温场中圆底扁薄球壳在均布外侧压力作用下的位移型几何非线性控制方程.考虑周边不可移简支边界条件,运用打靶法计算获得了不同几何参数的扁球壳轴对称弯曲变形的数值结果.定义了壳体临界几何参数.考察了壳体几何参数对平衡路径和临界荷载的影响.当壳体几何参数大于壳体临界几何参数时,上临界荷载随几何参数的增加单调增加,下临界荷载在很小范围内随几何参数的增加而增加,之后随几何参数的增加而减小.给定几何参数时,考察了不同均匀温度变化对壳体临界几何参数、临界荷载和平衡构型的影响.均匀升温使上临界荷载显著增加,使下临界荷载和临界几何参数显著减小.     

关于j_N方法的收敛性

jN方法是Asaoka提出的,故也称为Asaoka方法.由于使用它在计算均匀平板、各向同性散射条件下的临界参数时显示了很好的效果.因此引起了人们的重视,文献[2]—[6]从不同角度讨论了这个方法.本文将讨论此方法在计算临界参数和临界通量的合理性,并给出收敛速度.     

On the linear stability analysis of a Maxwell fluid with double-diffusive convection

The problem of double-diffusive convection and cross-diffusion in a Maxwell fluid in a horizontal layer in porous media is re-examined using the modified Darcy–Brinkman model. The effect of Dufour and Soret parameters on the critical Darcy–Rayleigh numbers is investigated. Analytical expressions of the critical Darcy–Rayleigh numbers for the onset of stationary and oscillatory convection are derived. Numerical simulations show that the presence of Dufour and Soret parameters has a significant effect on the critical Darcy–Rayleigh number for over-stability. In the limiting case some previously published results are recovered.     

An analytical study of linear and nonlinear double diffusive convection in a fluid saturated anisotropic porous layer with Soret effect

The double diffusive convection in a horizontal anisotropic porous layer saturated with a Boussinesq fluid, which is heated and salted from below in the presence of Soret coefficient is studied analytically using both linear and nonlinear stability analyses. The normal mode technique is used in the linear stability analysis while a weak nonlinear analysis based on a minimal representation of double Fourier series method is used in the nonlinear analysis. The generalized Darcy model including the time derivative term is employed for the momentum equation. The critical Rayleigh number, wavenumber for stationary and oscillatory modes and frequency of oscillations are obtained analytically using linear theory. The effect of anisotropy parameters, solute Rayleigh number, Soret parameter and Lewis number on the stationary, oscillatory, finite amplitude convection and heat and mass transfer are shown graphically.     

A nonlinear stability analysis of a rotating double-diffusive magnetized ferrofluid

A nonlinear stability result for a double-diffusive magnetized ferrofluid layer rotating about a vertical axis for stress-free boundaries is derived via generalized energy method. The mathematical emphasis is on how to control the nonlinear terms caused by magnetic body and inertia forces. The result is compared with the result obtained by linear instability theory. The critical magnetic thermal Rayleigh number given by energy theory is slightly less than those given by linear theory and thus indicates the existence of subcritical instability for ferrofluids. For non-ferrofluids, it is observed that the nonlinear critical stability thermal Rayleigh number coincides with that of linear critical stability thermal Rayleigh number. For lower values of magnetic parameters, this coincidence is immediately lost. The effect of magnetic parameter, M₃, solute gradient, S₁, and Taylor number, T_A1, on subcritical instability region have been analyzed. We also demonstrate coupling between the buoyancy and magnetic forces in the nonlinear stability analysis.     

Inequalities for a kinematic dynamo

A lower bound is determined for a functional occurring in dynamo theory. The significance of this bound for the critical Rayleigh number and heat source for a magnetohydrodynamic dynamo is then deduced from work of Proctor [4].     

Oscillatory convection and the Cattaneo law of heat conduction

The problem of thermal convection is investigated for a layer of fluid when the heat flux law of Cattaneo is adopted. The boundary conditions are those appropriate to two fixed surfaces. It is shown that for small Cattaneo number the critical Rayleigh number initially increases from its classical value of 1707.765 until a critical value of the Cattaneo number is reached. For Cattaneo numbers greater than this critical value a notable Hopf bifurcation is observed with convection occurring at lower Rayleigh numbers and by oscillatory rather than stationary convection. The aspect ratio of the convection cells likewise changes.     

Effect of dust particles on a layer of micropolar ferromagnetic fluid heated from below saturating a porous medium

The problem of the effect of dust particles on the thermal convection in micropolar ferromagnetic fluid saturating a porous medium subject to a transverse uniform magnetic field has been investigated theoretically. Linear stability analysis and normal mode analysis methods are used to find an exact solution for a flat micropolar ferromagnetic fluid layer contained between two free boundaries. In case of stationary convection, the effect of various parameters like medium permeability, dust particles, non-buoyancy magnetization, coupling parameter, spin-diffusion parameter and micropolar heat conduction parameter are analyzed. For sufficiently large values of magnetic parameter M₁, the critical magnetic thermal Rayleigh number for the onset of instability is determined numerically and results are depicted graphically. It is also observed that the critical magnetic thermal Rayleigh number is reduced solely because the heat capacity of clean fluid is supplemented by that of the dust particles. The principle of exchange of stabilities is found to hold true for the micropolar ferromagnetic fluid saturating a porous medium heated from below in the absence of micropolar viscous effect, microinertia and dust particles.     

Energy stability in the Bénard problem for a fluid of second grade

The critical Rayleigh number for stability of convection in a fluid of second grade is obtained by means of the energy method. If the thermodynamic restrictions derived by Dunn and Fosdick [1] are adopted the nonlinear energy boundary is shown to be the same as that of linear theory, and the nonexistence of subcritical instabilities is deduced.     

Electrohydrodynamic instability of a rotating layer of a viscoelastic fluid heated from below

The stability of a rotating layer of viscoelastic dielectric liquid (Walters liquid B) heated from below is considered. Linear stability theory is used to derive an eigenvalue system of ten orders and exact eigenvalue equation for a neutral instability is obtained. Under somewhat artificial boundary conditions, this equation can be solved exactly to yield the required eigenvalue relationship from which various critical values are determined in detail. Critical Rayleigh heat numbers and wavenumber for the onset of instability are presented graphically as function of the Taylor number for various values of electric Rayleigh number and the elastic parameters.     

Stability analysis of Rayleigh�CB��nard convection in a porous medium

In this paper we study the problem of Rayleigh?CBénard convection in a porous medium. Assuming that the viscosity depends on both the temperature and pressure and that it is analytic in these variables we show that the Rayleigh?CBénard equations for flow in a porous media satisfy the idea of exchange of stabilities. We also show that the static conduction solution is linearly stable if and only if the Rayleigh number is less than or equal to a critical Rayleigh number. Finally, we show that a measure of the thermal energy of the fluid decays exponentially which in turn implies that the L² norm of the perturbed temperature and velocity also decay exponentially.     

Effect of rotation on the onset of double diffusive convection in a Darcy porous medium saturated with a couple stress fluid

The effect of rotation on the onset of double diffusive convection in a horizontal couple stress fluid-saturated porous layer, which is heated and salted from below, is studied analytically using both linear and weak nonlinear stability analyses. The extended Darcy model, which includes the time derivative and Coriolis terms, has been employed in the momentum equation. The onset criterion for stationary, oscillatory and finite amplitude convection is derived analytically. The effect of Taylor number, couple stress parameter, solute Rayleigh number, Lewis number, Darcy–Prandtl number, and normalized porosity on the stationary, oscillatory, and finite amplitude convection is shown graphically. It is found that the rotation, couple stress parameter and solute Rayleigh number have stabilizing effect on the stationary, oscillatory, and finite amplitude convection. The Lewis number has a stabilizing effect in the case of stationary and finite amplitude modes, with a destabilizing effect in the case of oscillatory convection. The Darcy–Prandtl number and normalized porosity advances the onset of oscillatory convection. A weak nonlinear theory based on the truncated representation of Fourier series method is used to find the finite amplitude Rayleigh number and heat and mass transfer. The transient behavior of the Nusselt number and Sherwood number is investigated by solving the finite amplitude equations using Runge–Kutta method.