A NUMERICAL PROCEDURE FOR PREDICTING MULTIPLE SOLUTIONS OF A SPHERICAL TAYLOR–COUETTE FLOW

A new numerical procedure for predicting multiple solutions of Taylor vortices in a spherical gap is presented. The steady incompressible Navier–Stokes equations in primitive variables are solved by a finite- difference method using a matrix preconditioning technique. Routes leading to multiple flow states are designed heuristically by imposing symmetric properties. Both symmetric and asymmetric solutions can be predicted in a deterministic way. The current procedure gives very fast convergence rate to the desired flow modes. This procedure provides an alternative way of finding all possible stable steady axisymmetric flow modes.     

NESTED GRID METHODS FOR AN OCEAN MODEL: A COMPARATIVE STUDY

In this paper a comparison is carried out between three correction methods for multigrid local mesh refinement in oceanic applications: FIC, LDC and the direct method (DM) proposed by Spall and Holland. This study is based on a nested primitive equation model developed by Laugier on the basis of the code OPA (LODYC). The external barotropic problem is solved using any of the three local grid correction algorithms yielding an interactive nested grid model. The non-linear elliptic equation for the barotropic streamfunction tendency is solved on two nested grids, called the global and the zoom grid, that interact between themselves. The zoom grid is entirely embedded within the global domain with a horizontal grid step ratio of 3:1. The computation on the global grid supplies the boundary conditions for the zoom grid region and the fine grid fields are used to correct the global coarse solution. The three local correction methods are tested on two problems relevant to oceanic circulation phenomena proposed by Spall and Holland: a barotropic modon and an anticyclonic vortex. The results show that the nesting technique is a very efficient way to solve these problems in terms of a gain in precision compared with the required CPU time. The two-domain model with local mesh refinement allows one both to manage effectively the open boundary conditions for the local grid and to correct the global solution thanks to the zoom solution. In the case of the modon propagation the three local correction methods provide approximately the same results. For the baroclinic vortex it appears that the two iterative methods are more efficient than the direct one.     

MULTILEVEL SCHEMES FOR SOLVING UNSTEADY EQUATIONS

In this paper we propose and study various multilevel schemes for solving unsteady equations. Numerical comparisons between the choices of discretization and discussions are made on the Burgers equation in one and two dimensions. In particular we prove the advantage of using a hierarchical ADI solver.     

Third-order Nonlinear Optical Properties of Octa-substituted Metal-free Phthalocyanine Thin Films

Third-order nonlinear optical susceptibility, χ⁽³⁾ of symmetrically octa-substituted metal-free phthalocyanine thin films measured by the third-harmonic generation technique are reported. The metal-free phthalocyanine has been found to show a χ⁽³⁾ (−3ω; ω,ω, ω) value as large as 7.73×10⁻¹² esu at 1.80 μm. The figure of merit, χ⁽³⁾/α, was estimated to be 4.17×10¹⁷ esu cm at 1.05 μm and 6.97×10¹⁶ esu cm at 1.65 μm. Both linear and third-order optical properties of liquid-crystalline metal-free phthalocyanines are discussed     

Dust-Acoustic Nonlinear Waves in a Nanoparticle Fraction of Ultracold (2K) Multicomponent Dusty Plasma

The nonlinear dust-acoustic instability in the condensed submicron fraction of dust particles in the low-pressure glow discharge at ultra-low temperatures is experimentally and theoretically investigated. The main discharge parameters are estimated on the basisof the dust-acoustic wave analysis. In particular, the temperature and density of ions, as well as the Debye radius, are determined. It is shown that the ion temperature exceeds the temperature of the neutral gas. The drift characteristics of all plasma fractions are estimated. The reasons for the instability excitation are considered.     

任意高差的多档输电线动刚度理论建模

针对目前多档输电线动力特性研究中存在的不足,将任意一档导线视为子结构并基于现有的索振动理论获得了其动张力;利用相邻两档导线的动张力建立了悬垂绝缘子串绕其悬挂点摆动的动力学方程;综合每一档导线的振动方程以及悬垂绝缘子串的摆动方程建立了两档及三档输电线的动力学方程,求解得到由端部位移激励产生的动刚度;通过对动刚度进行分析得到了多档输电线的固有频率以及模态的理论表达形式;建立了小高差的两档及三档导线的有限元模型,理论分析结果与有限元结果一致,验证了理论公式的正确性。本文提出的动刚度理论对于研究多档输电线的动态响应具有重要的应用价值。     

大型偏心结构吊装欠约束问题研究

吊装施工过程中被吊模块的水平度是作业要求的重要指标,通常需要增加配重调平。传统有限元方法需要补充约束以消除单元刚体位移,且需要重复计算平衡方程来求解调平载荷,效率不高。将模块的运动分解为随动坐标系的整体运动以及相对该坐标系的弹性变形,可将欠约束问题化为多体系统的静平衡问题。基于虚功率原理推导了吊装平顺时刻的节点力平衡方程以及相应的切线刚度矩阵,并将配重表示为基础配重与载荷系数相乘的形式。通过对节点力平衡方程求导,得到一组以载荷系数为自变量的微分方程,通过求解微分方程并结合水平度判据,可快速搜寻满足水平度要求的载荷系数。数值算例表明,该方法在解决偏心模块吊装欠约束问题方面具有明显的优势,在确定配重载荷方面具有较快的速度和合理的精度。     

时间尺度上二阶拟线性阻尼动力方程的振动性分析

研究二阶拟线性时滞阻尼动力方程[a(t)|x^Δ(t)|^λ-1x^Δ(t)]^Δ+b(t)|x^Δ(t)|^λ-1x^Δ(t)+p(t)|x(δ(t))|^λ-1x(δ(t))=0的振动性,其中t ₀∈T,而T为任意时间尺度,考虑方程是非正则情形,即$\int^{∞}_{t_0}$[a^-1(s)e_-b/a(s,t₀)]^1/λ Δs < ∞。 通过引入广义Riccati变换,借助时间尺度上的有关理论,并结合不等式技巧,建立了该方程振动的一些新的充分条件,推广、改进并丰富了现有文献中的结果。     

两空间耦合下齿轮传动系统多稳态特性研究

通过将系统参数定义为参数变量, 构成参数空间,研究齿轮传动系统在参数空间和状态空间耦合下的非线性全局动力学特性,以及多参数、多初值和多稳态行为之间的关联特性.首先设计了一个两空间耦合下非线性系统多稳态行为的计算和辨识方法.其次,基于该方法并结合相图、Poincaré映射图、分岔图、最大Lyapunov指数、吸引域等,研究齿轮传动系统在不同参数平面上多稳态行为的存在区域和分布特性,以及多稳态行为在状态平面上的分布特性,揭示了参数平面和状态平面上系统可能隐藏的多稳态行为和分岔,并分析了多稳态行为的形成机理. 结果发现,两空间耦合下系统在参数平面上存在大量多稳态行为并呈"带状"分布, 状态平面上多稳态行为出现两种不同的侵蚀现象, 即内部侵蚀和边界侵蚀.分岔点或分岔曲线对初值的敏感性导致多稳态行为的出现.当齿侧间隙和误差波动在较小的范围内变化时,系统全局动力学特性受间隙和误差扰动的影响较小,受啮合频率的影响较大.两空间耦合下系统全局动力学特性变得丰富和复杂.