分数阶FitzHugh-Nagumo模型神经元的动力学特性及其同步

通过对分数阶FitzHugh-Nagumo模型神经元的研究,当外加电流强度作为分岔参数时,发现这种模型神经元从静息态到周期放电态所经历的Hopf分岔点不同于相应的整数阶模型神经元的分岔点;而且分数阶FitzHugh-Nagumo模型神经元呈现周期放电的外加电流强度的范围比相应的整数阶模型神经元的范围小,然而放电频率却比相应的整数阶模型神经元的放电频率高.同时还揭示在周期放电的情况下分数阶FitzHugh-Nagumo模型神经元之间的同步速率比相应的整数阶模型神经元之间的同步速率快.在数值模拟分数阶微分方程 关键词： 分数阶 Hopf分岔 FitzHugh-Nagumo模型 同步     

过阻尼分数阶Langevin方程及其随机共振

通过对广义Langevin方程阻尼核函数的适当选取,在过阻尼的情形下, 推导出分数阶Langevin方程.给合反常扩散理论和分数阶导数的记忆性, 讨论了分数阶Langevin方程的物理意义,进而得出分数阶Langevin方程产生随机共振的内在机理.数值模拟表明,在一定的阶数范围内,分数阶Langevin方程可以产生随机共振, 并且分数阶下的信噪比增益好于整数阶情形.     

时间非对称外力驱动分数阶布朗马达的定向输运

本文研究了周期对称势中时间非对称外力驱动的布朗粒子输运现象, 建立了分数阶布朗马达输运模型. 其中外力是零均值的, 而分数阶阶数则刻画了客观环境的非均匀性程度. 通过将模型离散化进行数值模拟, 讨论了分数阶阶数、系统参量和外部参量与定向流之间的依赖关系. 研究表明, 即使没有倾斜势场的作用, 时间非对称外力也可以诱导系统产生定向输运; 输运速度随分数阶阶数的增大而单调递增; 当阶数固定时, 系统的输运速度会随着势垒高度、噪声强度非单调变化, 表现出广义随机共振现象. 分析指出, 分数阶郎之万方程所刻画的输运现象是在整数阶模型基础上的一个推广, 进而为输运现象提供了一个可能更为真实的模型.     

分数阶布朗马达在闪烁棘齿势中的合作输运现象

引入分数阶微积分理论,建立耦合分数阶布朗马达在闪烁棘齿势中的合作输运模型, 利用分数阶差分法求得模型数值解并分析了模型参数对合作定向输运性质的影响. 发现在具有记忆性的分数阶棘齿系统中, 系统阶数与粒子间耦合强度不仅可影响粒子链输运速度, 还可使粒子链出现与整数阶方向相反的定向流; 在阶数固定下, 定向输运速度将随参数(噪声强度、耦合强度、棘齿势峰值高度)变化出现广义随机共振现象. 关键词： 分数阶布朗马达 闪烁棘齿势 合作定向输运 广义随机共振     

分数阶Brown马达及其定向输运现象

选取幂函数作为广义Langevin方程的阻尼核函数,采用闪烁棘轮势,建立了过阻尼分数阶Brown马达模型.结合分数阶微积分的记忆性,分析了粒子在过阻尼分数阶Brown马达作用下的运动特性.研究发现,较之整数阶情形,过阻尼分数阶Brown马达也会产生定向输运现象,并且在某些阶数下会产生整数阶情形所不具有的反向定向流.此外,还讨论了阶数和噪声强度对系统输运速度的影响,发现当阶数固定时,其平均输运速度会随噪声变化出现随机共振;当噪声强度固定时,其输运速度会随阶数变化而振荡,即出现多峰的广义随机共振现象.     

分数阶Oldroyd-B黏弹性Poiseuille流的Laplace数值反演分析

采用Laplace数值反演的Stehfest算法研究了分数阶Oldroyd-B粘弹性流体在两平板间非定常的Poiseuille流动问题.首先,通过数值解与近似解析解的比较验证了Stehfest算法的有效性.其次,运用Stehfest算法对平板Poiseuille流动进行了研究,揭示了分数阶黏弹性平板流的速度过冲和应力过冲现象,指出这些现象对分数导数的阶数存在明显的依赖性.同时,数值结果表明,整数阶本构方程仅仅是分数阶本构方程的特例,分数阶本构方程较整数阶本构方程具有更广泛的适用性。     

紧聚焦角向偏振分数阶涡旋光诱导磁化场特性

基于理查德-沃尔夫矢量衍射理论和逆法拉第效应, 首次推导出紧聚焦角向偏振分数阶涡旋光诱导磁化场表达式, 分析出角向偏振分数阶涡旋光诱导磁化场可近似等效为有限个相邻整数阶涡旋角向偏振光诱导磁化场与其交叉诱导磁化场的加权叠加. 数值模拟了紧聚焦条件下不同分数阶涡旋角向偏振光诱导磁化场的分布特性. 模拟结果表明角向偏振分数阶涡旋光诱导磁化场呈非对称分布. 当取分数阶涡旋拓扑荷α∈ [0.5, 1.5]时, 随着涡旋拓扑荷的增加, 磁化场横向平面分布出现分裂现象, 同时还引起磁斑垂直于光轴方向的自移效应. 当α = 0.5 或 1.5 时, 磁斑中心最大偏移量达0.24λ; 当分数阶涡旋拓扑荷α∈ (2, 3]时, 磁化场横向平面分布分裂出 2 个强度热点, 强度分布环径逐渐变大. 当分数阶涡旋拓扑荷趋于α = 3整数时, 磁化横向平面分布也趋于圆对称性分布. 特别有趣的是, 当分数阶涡旋拓扑荷取半整数, 尤其大于 3 时, 磁化场强度热点数与其包围的暗点数均与涡旋阶数存在正相关关系, 其中热点数为α - 0.5, 暗点数为α - 1.5. 这些研究结果将可应用于全光磁记录和磁性粒子的捕获与操控等.     

空时非对称分数阶类Langevin棘齿

研究了空时非对称分数阶类Langevin分子马达棘齿模型,其中势函数是空间对称破缺的周期势,时间非对称类Langevin噪声由Logistic映射生成,而分数阶则刻画了分子马达工作环境的非理想程度.通过将模型转化为离散映射,即研究其整时间点情形,数值模拟了噪声的时间非对称性、势函数的空间非对称性以及分数阶对模型定向输运行为的影响.数值模拟结果表明:噪声的时间非对称性是定向流产生的根源,而势函数的空间非对称性能够与其进行竞争与协作,并在适当的参数条件下导致定向流的逆转;分数阶仅影响定向流的大小而不改变其方向.与经典的整数阶分子马达模型或时间非对称分数阶分子马达棘齿模型相比,该模型可以更为真实地描述分子马达的噪声整流工作机理.     

分数阶双头分子马达的欠扩散输运现象

研究具有幂律记忆性的细胞液中双头分子马达的定向输运现象,选取幂函数作为广义Langevin方程的阻尼核函数,建立了分数阶过阻尼耦合Brown马达模型,讨论了阶数及耦合系数对双头分子马达定向输运速度的影响. 仿真结果表明,分数阶过阻尼双头分子马达也会产生定向输运现象,并且在某些阶数下会产生整数阶情形所不具有的反向定向流. 当噪声强度固定时,输运速度随着阶数以及耦合系数的变化均会出现广义随机共振现象. 特别地,研究发现双头分子马达在记忆闪烁棘轮势中具有某些单头分子马达所不具备的运动特性,定向流的大小和方向由噪声与双头间作用力相互耦合控制. 关键词： 分数阶双头分子马达 欠扩散 定向输运 广义随机共振     

间歇湍流的分数阶动力学

间歇湍流意味着湍流涡旋并不充满空间,其维数介于2和3之间.湍流扩散为超扩散,且概率密度分布具有长尾特征.本文将流体力学的Navier-Stokes(NS)方程中的黏性项用分数阶的拉普拉斯算子表达.分析表明,分数阶拉普拉斯的阶数α和间歇湍流的维数D相联系.对于均匀各向同性的Kolmogorov湍流α=2,即用整数阶NS方程描述.而对于间歇性湍流,一定用分数阶的NS方程来描述.对于Kolmogorov湍流,扩散方差正比于t3,即Richardson扩散.而对于间歇性湍流,扩散方差要比Richardson扩散更强.     

Continuum description of anomalous diffusion on a comb structure

Anomalous diffusion on a comb structure consisting of a one-dimensional backbone and lateral branches (teeth) of random length is considered. A well-defined classification of the trajectories of random walks reduces the original problem to an analysis of classical diffusion on the backbone, where, however, the time of this process is a random quantity. Its distribution is dictated by the properties of the random walks of the diffusing particles on the teeth. The feasibility of applying mean-field theory in such a model is demonstrated, and the equation for the Green’s function with a partial derivative of fractional order is obtained. The characteristic features of the propagation of particles on a comb structure are analyzed. We obtain a model of an effective homogeneous medium in which diffusion is described by an equation with a fractional derivative with respect to time and an initial condition that is an integral of fractional order. Zh. éksp. Teor. Fiz. 114, 1284–1312 (October 1998)     

一类分数阶混沌系统的研究

推广了一类分数阶混沌系统并证明了这类分数阶混沌系统的拓扑等价性,指出分数阶系统产生混沌吸引子的必要条件是系统平衡点的稳定性不变.通过数值模拟表明,此分数阶系统与整数阶系统一样仍然产生一对互不相交的2-涡卷混沌吸引子. 关键词： 分数阶 混沌 稳定性     

Oxidation of titanium between 25°C and 400°C

The oxidation of Ti between 25°C and 400°C has been studied in an ultra high vacuum system with ellipsometry, Auger spectroscopy and surface potential difference. The surface potential difference is primarily sensitive to the adsorbed layer of oxygen, the Auger spectrometer is sensitive to the oxygen and Ti concentration at the oxide-gas interface and the ellipsometer is sensitive to the film thickness. Consequently, it was possible to follow the various processes separately during oxidation and oxide dissolution. Oxygen adsorbs in a percusor state as a molecule which dissociates and adsorbs in the atomic state. Reaction of atomic oxygen with Ti at the oxide-gas interface causes the stoichiometry to vary with the temperature, oxygen pressure and film thickness. At low temperatures and oxygen pressure the initial film growth follows a logarithmic time law. The predominant transport species is cationic. After initial oxidation, at higher temperatures and oxygen pressure, the predominant transport species is anionic. At constant temperature and oxygen pressure the oxide grows until the rate of growth is equal to the rate of dissolution into the metal. At this point the film thickness remains constant with time. Doping the Ti surface with Au results in a decrease in the oxide dissolution process due to blocking of interstitial diffusion paths by Au atoms.     

色散渐减光纤组成的环形镜对高阶孤子的理想压缩

利用数值模拟方法证明,采用色散渐减光纤组成的非线性光纤环形镜不仅可压缩高阶光孤子,而且能有效地消除压缩后脉冲的底座、提高光脉冲的输出能量.研究表明,对于一确定的色散渐减光纤,只要选取不同环形镜长度,即可对不同阶数的高阶孤子进行理想压缩.且孤子阶数越高,所需最佳环形镜长度越短、压缩后光脉冲的峰值强度越大、能量透射率越低.研究同时发现,环形镜的功率分束比存在一最佳值,在此值下所得压缩后的光脉冲不仅具有较大的峰值及能量透射率,且不含底座.     

一种新的分数阶系统稳定理论及在back-stepping方法同步分数阶混沌系统中的应用

对阶次小于1的分数阶系统的稳定问题,提出了一种新的判据.基于该理论,将back-stepping方法拓展到分数阶系统中并设计控制器实现了分数阶Newton-Leipnik混沌系统的同步.数值仿真验证了稳定性理论的准确性及控制器设计方法的有效性. 关键词： 分数阶 同步 back-stepping方法 Newton-Leipnik     

On the dynamics,existence of chaos,control and synchronization of a novel complex chaotic system

In this paper, the authors have studied the dynamics of a novel complex chaotic system with fractional order derivative and found the existence of chaos. The novel complex system is simulated for integer as well as fractional orders which shows some unusual phenomena. The main contribution of this effort is an implementation of the Largest Lyapunov Exponent (LLE) criteria based on Wolf’s algorithm. The conditions for chaos control based on the fractional Routh–Hurwitz stability conditions and feedback control are given. Also synchronization between a fractional order novel chaotic system and a controlled fractional order novel system using the modified adaptive projective synchronization method for different scaling matrices has been obtained. Numerical simulation results are carried out using the Adams–Bashforth–Moulton method.     

Matrix approach to discrete fractional calculus II: Partial fractional differential equations

A new method that enables easy and convenient discretization of partial differential equations with derivatives of arbitrary real order (so-called fractional derivatives) and delays is presented and illustrated on numerical solution of various types of fractional diffusion equation. The suggested method is the development of Podlubny’s matrix approach [I. Podlubny, Matrix approach to discrete fractional calculus, Fractional Calculus and Applied Analysis 3 (4) (2000) 359–386]. Four examples of numerical solution of fractional diffusion equation with various combinations of time-/space-fractional derivatives (integer/integer, fractional/integer, integer/fractional, and fractional/fractional) with respect to time and to the spatial variable are provided in order to illustrate how simple and general is the suggested approach. The fifth example illustrates that the method can be equally simply used for fractional differential equations with delays. A set of MATLAB routines for the implementation of the method as well as sample code used to solve the examples have been developed.     

Influence of sintering temperature and pH on the phase transformation,particle size and anti-reflective properties of RHA nano silica powders

Nano silica powders were synthesized from rice husk ash, the most silica-rich raw material, using alkaline extraction followed by acid precipitation. The phase transformation during sintering, the influence of sintering temperature and pH on the particle size and anti-reflective properties of nano silica were investigated. The results showed that the amorphous SiO₂ sintered at 600°C were transformed to a cristobalite structures completely during the sintering process at 800°C and 1100°C. With the increasing sintering temperature and pH, the particle size distributions (d₅₀) were increased respectively in the range of 62–84, 192–240, and 283–329?nm at sintering temperatures of 600°C, 800°C, and 1100°C. When the sintering temperatures were increased at 1100°C, 98.15% and 96.84% of transmittances were obtained respectively at the highest and lowest points of the anti-reflection band and could be used for anti-reflective applications.     

分数阶系统的一种稳定性判定定理及在分数阶统一混沌系统同步中的应用

根据Lyapunov稳定定理及其逆定理和分数阶系统稳定定理,提出了如果整数阶系统稳定,其对应的阶次小于1的分数阶形式的系统也稳定的分数阶系统稳定的判定定理,并给出了详细的证明过程.并将该理论运用于分数阶混沌系统的同步,实现了未知参数分数阶统一混沌系统的自适应同步,仿真结果证实了该理论的正确性. 关键词： 分数阶系统 混沌 Lyapunov稳定定理 Lyapunov稳定逆定理     

Analysing the sintering of heterogeneous powder structures by in situ microtomography

This work analyses the microstructure changes of various copper-based powder systems during sintering from 3D images provided by in situ synchrotron microtomography. The investigated systems include a copper powder with a wide particle size distribution of 0–63 µm poured into a quartz capillary, a pre-sintered copper compact with artificially created large pores and a mixture of copper and alumina particles. The experiments were carried out at the European Synchrotron in Grenoble, France. Powders were sintered up to 1060°C under reducing atmosphere in a furnace located between the X-ray source and the detector. During each experiment, 3D images were taken at various times of the thermal cycle. We have obtained images with a resolution of 1.5 µm and the time of acquisition of every image was ～1 min. Quantitative analysis of these images allowed the changes of various important parameters to be followed. Such parameters characterise the sintering process at the particle length scale: interparticle coordination, pore size distribution and particle centre-to-centre distance. Moreover, by tracking the displacement of each particle centre and comparing it to the displacement predicted by classical mean field assumption, we have been able to assess the magnitude of particle rearrangement occurring during sintering. From these data, the sintering behaviour of heterogeneous powder systems is discussed with particular emphasis of collective particle phenomena.