方腔中多孔介质自然对流分叉的数值模拟

本文用三次样条方法对充满饱和多孔介质的方腔的底部加热时出现的对流分叉现象进行了数值模拟,得到三个分叉点,观察到了分叉点邻域内流动发展状态,对稳定性与挑动幅度关系做了初步的探讨。     

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

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

底部热加载下两层多孔介质热流耦合现象研究

底部热加载条件下非线性流动传热问题的求解一直是多孔介质研究领域的难题。采用实验测试及数值模拟的方法,对底部热加载方式下两层多孔介质内热流耦合对流传热解的特性进行了研究。研究结果表明:在小瑞利数工况下,两层多孔介质内骨架内材料比例对温度场分布和非线性特性有重要影响;确定了非线性分叉、震荡解的特征值,并计算出了不同材料比例下的分叉及震荡所对应的临界瑞利数;通过实验验证了底部热加载时两层多孔介质内温度随时间呈非稳态震荡变化。结论可为实现热流分层、局部削弱或强化传热提供参考。     

基于能量释放率的Ⅱ型裂纹分叉研究

裂纹尖端存在奇异应力场,该类奇异应力场所具有的高度应力集中将导致裂纹开裂.本文应用典型的J积分理论来划分裂纹尖端的积分路径,基于能量释放率理论对Ⅱ型裂纹尖端的复杂分叉情况进行研究.通过所建立的断裂模型求解出Ⅱ型裂纹多种分叉形式的能量释放率的解析解;导出了各种分叉构型的能量驱动力;提出了各种分叉构型的K-型开裂准则;给出了裂纹分叉的临界开裂角;确定了裂纹的分叉韧性与断裂韧性之间的关系.通过本文研究方法得出的裂纹分叉形式及裂纹分叉临界开裂角与已有实验结果十分吻合.     

非线性问题和分叉问题及其数值方法

本文给出了一个一般性的分叉定义，说明了伪弧长算法在分叉计算中的应用，概述了静分叉点定位、用单纯形算法准确确定静分叉后各分叉解枝初始方向的算法，以及Ｈｏｐｆ分叉点定位和大范围连续追踪周期解轨道的数值方法。     

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

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

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

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

多孔介质中热对流的分叉机理研究

本文利用解析分析方法研究了数值模拟发现的多孔介质层中出现的对流分叉机理，指出控制方程中的Ｒａｙｌｅｉｇｈ数，是决定流动的特征参数。当Ｒａｙｌｅｉｇｈ数小于临界数值时，多孔介质内流动处于静止传热状态，并且这种状态是稳定的。如果Ｒａｙｌｅｉｇｈ数大于临界数值，非线性方程出现分叉解，文中指出，存在多个使平凡解失稳而分叉的临界Ｒａｙｌｅｉｇｈ数，当Ｒａｙｌｅｉｇｈ数由小到大经历这些临界数值时，其由平凡解发展起来的分叉解的流态，依次由单回流区转变为双回流区及三回流区。理论分析给出了分叉解和分叉解的振幅方程，阐明了分叉的机理，其结论和数值结果定性一致．     

非线性Mathieu方程1／2亚谐分叉解的实验研究

本文对一类Mathieu方程的1/2亚谐分叉特性进行了实验研究,得出了在整个参数平面上具有不同拓朴结构分叉图的实验曲线,研究了确定非线性系统衰减参数的方法。并对各种特定的物理系统,可能出现的不同拓朴结构的分叉图和所具有的不同参数区域进行了讨论。     

不可压缩流中机翼颤振稳定性研究

研究两个自由度的机翼在不可压缩流作用下颤振的分支问题.运用罗司-霍维茨判据确定系统的分叉点,应用中心流形理论将四维系统降为二维系统,用直接求周期解方法对分叉点的真假中心及稳定性问题进行了分析,并研究了系统的极限环颤振.结果表明,本文研究的分叉点不是中心,而是稳定或不稳定焦点.在两个分叉点处,系统发生了超临界和亚临界Hopf分叉,产生稳定或不稳定极限环.     

DC-DC开关功率变换器的非线性动力学行为研究

DC-DC开关功率变换器是一种典型的分段光滑动力学系统, 在一定的工作和参数条件下, 系统会出现各种分岔如倍周期分岔、Hopf分岔、边界碰撞分岔和混沌运动. 系统评述了DC-DC开关功率变换器的非线性动力学行为的研究进展;介绍了离散非线性映射、分段线性模型、平均值模型等3种建模方法;分析了这种电路系统中的分岔特点及通向混沌的途径与机制;结合我们的研究工作, 讨论了对这种电路系统进行混沌控制的必要性及相关策略;最后, 从应用的角度提出了未来的若干研究方向.      

宏观交通流模型的余维2分岔分析

交通流特性是混合交通流建模的一个重要因素. 交通流模型中的分岔现象是导致复杂交通现象的因素之一. 交通流的分岔, 涉及复杂的动力学特征且研究较少. 因此, 提出了一个最优速度模型来研究驾驶员记忆对驾驶行为的影响. 基于带有记忆的最优速度连续交通流模型, 利用非线性动力学, 分析和预测了复杂交通现象. 推导了鞍结 (LP) 分岔存在条件, 并通过数值计算得到了余维1 Hopf (H) 分岔、LP分岔和同宿轨 (HC) 分岔以及余维2广义Hopf (GH) 分岔、尖点 (CP) 分岔和Bogdanov-Takens (BT) 分岔等多种分岔结构. 根据双参数分岔区域的特点, 研究了记忆参数对单参数分岔结构的影响, 分析了不同分岔结构对交通流的影响, 并用相平面描述了平衡点附近轨迹的变化特征. 选择Hopf分岔和鞍结分岔作为密度演化的起点, 描述了均匀流、稳定和不稳定的拥挤流以及走走停停现象. 结果表明, 驾驶员记忆对交通流的稳定性有重要影响; 动力学行为很好地解释了交通拥堵现象; 考虑余维2分岔的影响, 能更好地理解交通拥堵产生的根源, 并为制定有效抑制拥堵的方法提供一定的理论依据.      

Bifurcations and chaos in a discrete predator–prey model with Crowley–Martin functional response

In this paper, a discrete-time predator–prey model with Crowley–Martin functional response is investigated based on the center manifold theorem and bifurcation theory. It is shown that the system undergoes flip bifurcation and Neimark–Sacker bifurcation. An explicit approximate expression of the invariant curve, caused by Neimark–Sacker bifurcation, is given. The fractal dimension of a strange attractor and Feigenbaum’s constant of the model are calculated. Moreover, numerical simulations using AUTO and MATLAB are presented to support theoretical results, such as a cascade of period doubling with period-2, 4, 6, 8, 16, 32 orbits, period-10, 20, 19, 38 orbits, invariant curves, codimension-2 bifurcation and chaotic attractor. Chaos in the sense of Marotto is also proved by both analytical and numerical methods. Analyses are displayed to illustrate the effect of magnitude of interference among predators on dynamic behaviors of this model. Further the chaotic orbit is controlled to be a fixed point by using feedback control method.     

Bifurcation and buckling analysis of a unilaterally confined self-rotating cantilever beam

A nonlinear dynamic model of a simple non-holonomic system comprising a self-rotating cantilever beam subjected to a unilateral locked or unlocked constraint is established by employing the general Hamilton's Variational Principle. The critical values, at which the trivial equilibrium loses its stability or the unilateral constraint is activated or a saddle-node bifurcation occurs, and the equilibria are investigated by approximately analytical and numerical methods. The results indicate that both the buckled equilibria and the bifurcation mode of the beam are different depending on whether the distance of the clearance of unilateral constraint equals zero or not and whether the unilateral constraint is locked or not. The unidirectional snap-through phenomenon (i.e. catastrophe phenomenon) is destined to occur in the system no matter whether the constraint is lockable or not. The saddle-node bifurcation can occur only on the condition that the unilateral constraint is lockable and its clearance is non-zero. The results obtained by two methods are consistent. The project supported by the National Natural Science Foundation of China (10272002) and the Doctoral Program from the Ministry of Education of China (20020001032) The English text was polished by Yunming Chen.     

Bifurcation boundary analysis as a nonlinear damage detection feature: does it work?

Many systems in engineering and science are inherently nonlinear and require damage detection. For such systems, nonlinear damage detection methods may be useful. A bifurcation boundary analysis method as a new nonlinear damage detection tool was previously introduced in the literature to track bifurcation boundary changes due to damages over a small region of an aeroelastic panel model. Results of this method based upon a finite difference solution showed higher sensitivities to the small amount of damage than methods based upon linear models. In this paper, four methods including Finite Difference, Finite Element (FEM), Rayleigh-Ritz and Galerkin Approach are used to further investigate the sensitivity of the bifurcation boundary for damage detection. Results of the FEM and Rayleigh-Ritz method agree with each other and also show that the sensitivity of the bifurcation boundary to damage is much less than what previously reported when using a finite difference solution method.     

Global bifurcations and chaos in a Van der Pol-Duffing-Mathieu system with three-well potential oscillator

Semi-analytical and semi-numerical method is used to investigate the global bifurcations and chaos in the nonlinear system of a Van der Pol-Duffing-Mathieu oscillator. Semi-analytical and semi-numerical method means that the autonomous system, called Van der Pol-Duffing system, is analytically studied to draw all global bifurcations diagrams in parameter space. These diagrams are called basic bifurcation diagrams. Then fixing parameter in every space and taking parametrically excited amplitude as a bifurcation parameter, we can observe the evolution from a basic bifurcation diagram to chaotic pattern by numerical methods. The project supported by the National Natural Science Foundation of China     

准周期响应对称破缺分岔点的一种快速计算方法

稳态响应如周期及准周期解的分岔计算, 是非线性动力学研究的难点问题之一. 与计算方法及分析理论相对完善的周期响应相比, 准周期响应的求解只是在近些年才得到较大进展, 而且其分岔分析更加棘手, 仍需要更有效的理论和方法. 目前, 稳态响应尤其是准周期响应的分岔计算, 一般需采用数值方法, 通过调节参数反复试算得到. 为此, 本文基于增量谐波平衡IHB法提出一种快速方法, 可以高效地确定准周期响应的对称破缺分岔点. 方法的理论基础是在准周期解的广义谐波级数表达基础上, 当响应发生对称破缺分岔时, 其偶次(含零次)谐波系数将逐渐由0变为小量. 基于此性质, 将零次谐波系数预先设定为小量, 同时将分岔控制参数视为可变的迭代变量, 进而通过IHB法构造迭代格式. 作为算例, 研究不可约频率作用下的双频激励Duffing系统以及Duffing-van der Pol耦合系统. 结果表明, 只要迭代格式收敛, 随着预设小量减小, 控制参数将逐渐接近分岔近似值; 同时, 通过提高谐波截断数可显著提高近似分岔值的计算精度. 所提方法无需反复试算, 只要迭代过程收敛、便可实现分岔点直接快速计算.      

Global analysis of secondary bifurcation of an elastic bar

In a three dimensional framework of finite deformation configurations, this paper investigates the secondary bifurcation of a uniform, isotropic and linearly elastic bar under compression in a large range of parameters. The governing differential equations and finite dimensional equations of this problem are discussed. It is found that, for a bar with two ends hinged, usually many secondary bifurcation points appear on the primary branches which correspond to the maximum bending stiffness. Results are shown on parameter charts. Secondary modes and branches are also calculated with numerical methods. The project supported in part by the National Natural Science Foundation of China     

Stability and bifurcation for a flexible beam under a large linear motion with a combination parametric resonance

Stability and bifurcation behaviors for a model of a flexible beam undergoing a large linear motion with a combination parametric resonance are studied by means of a combination of analytical and numerical methods. Three types of critical points for the bifurcation equations near the combination resonance in the presence of internal resonance are considered, which are characterized by a double zero and two negative eigenvalues, a double zero and a pair of purely imaginary eigenvalues, and two pairs of purely imaginary eigenvalues in nonresonant case, respectively. The stability regions of the initial equilibrium solution and the critical bifurcation curves are obtained in terms of the system parameters. Especially, for the third case, the explicit expressions of the critical bifurcation curves leading to incipient and secondary bifurcations are obtained with the aid of normal form theory. Bifurcations leading to Hopf bifurcations and 2-D tori and their stability conditions are also investigated. Some new dynamical behaviors are presented for this system. A time integration scheme is used to find the numerical solutions for these bifurcation cases, and numerical results agree with the analytic ones.     

Unfolding of multiparameter equivariant bifurcation problems with two groups of state variables under left-right equivalent group

Introduction Intheunfoldingtheoryofbifurcationproblems,Refs.[1-6]providedvariousversionsof theversalunfoldingtheorem.Itispointedoutthattheequivalentrelationadoptedintheabove mentionedreferencesiscontactequivalencederivedfromthesingularitytheoryofsmoothmap…