弱的参数周期扰动对一非线性系统安全域的影响与分形侵蚀安全域的控制

讨论了弱参数周期扰动对于非线性系统安全域的影响,在一定频率下的参数周期扰动将加速安全域的分形侵蚀,而在另一些频率的扰动下,将抑制安全域的分形侵蚀,并且存在着增进安全域的最优频率.提出了用弱参数周期扰动控制受到分形侵蚀的安全域的方法,并用Melnikov方法进行了分析.最后讨论了这种控制安全域的方法在实际环境中当具有外加噪声时的鲁棒性. 关键词： 分形吸引域边界 参数的周期扰动 控制     

耗散粒子动力学处理复杂固体壁面的一种有效方法

耗散粒子动力学(dissipative particle dynamics,DPD)作为一种介观尺度拉格朗日型粒子方法,已经成功地应用于微纳米流动和生化科技的研究中.复杂固体壁面的处理和壁面边界条件的实施一直是DPD方法发展及应用的一个障碍.提出了处理复杂固体壁面的一种新的方法.复杂固体区域通过冻结随机分布并且达到平衡状态的DPD粒子代表;所冻结的DPD粒子位于临近流动区域的一个截距内;在靠近固体壁面的流动区域中设置流动反弹层,当流动DPD粒子进入此流动层后反弹回流动区域.应用这种固体壁面处理方法对简单流动区域的Poiseuille流动和复杂多孔介质内的流动进行了分析.研究表明,这种新的固体壁面处理方法能够有效模拟复杂固体区域,准确实施壁面边界条件.     

消声器声学性能预测的子结构快速多极子边界元法

鉴于快速多极子边界元法的应用主要局限于单区域声学问题计算,发展基于子结构技术的快速多极子边界元法以计算多区域声场问题,介绍基本原理、具体实施过程以及优缺点.以带有插进口管的膨胀腔消声器为例,应用子结构快速多极子边界元法和传统边界元法计算其传递损失,通过与实验测量结果的比较,验证方法的有效性和计算精度.研究表明,快速多极子边界元法与传统边界元法相比,节点数越多,其在节省计算时间,减少计算量等方面的优势越明显.     

HT-7超导托卡马克边界等离子体参量及其涨落的实验研究

利用径向可移动朗缪尔三探针和马赫探针对HT-7超导托卡马克边界等离子体参量及其涨落进行了空时空分辨测量.给出了欧姆放电及其与低杂波电流驱动共同作用下边界等离子体电位φ_p、电子温度T_e和电子密度n_e及其涨落的径向分布.实验表明,在限制器附近,存在一由E_r×B确定的极向旋转速度剪切层.在剪切层内,T_e和n_e分布较陡,且φ_p,T_e和n_e的相对涨落水平下降明显.这说明剪切层对边界区的等离子体涨落具有抑制作用.低杂波驱动使径向电场梯度变陡,从而使剪切程度加深但对剪切层宽度无影响.此外,测量表明等离子体环向速度马赫数M_φ存在径向梯度.环向流的这种径向梯度可能是形成径向电场所需的平均极向流的一种重要驱动源 关键词： 托卡马克 边界涨落 低杂波驱动     

颗粒与顺列管束磨损的数值模拟

采用直接数值模拟方法模拟了在管道中直径80μm的颗粒对顺列10×10管束的碰撞和磨损.圆管与流场流体之间的相互作用采用内嵌边界方法进行计算并得到流场结构.颗粒的运动采用了拉格朗日跟踪方法计算并与流场之间进行了双向耦合.计算分析了沿流向以及垂直于流向各排管束受到颗粒的碰撞产生的磨损量.计算结果表明沿流向第三排管束以后的各排管束应采取防磨措旌.     

有限元/边界元耦合方法在海洋生物目标强度特性研究的应用*

目标强度特性是海洋生物声学识别与资源量评估的重要依据,其中,基于近似几何体和声阻抗特性的理论模型法是研究海洋生物目标强度的重要手段。由于对几何形态近似处理以及数值求解方法的限制,传统理论模型对声波频率、入射方位以及目标声阻抗、形态尺寸等均有各自不同的适用范围,单一模型难以满足不同种类或同一种类但不同尺寸海洋生物的目标强度求解。本文尝试将逐渐见诸应用的有限元/边界元耦合方法用于海洋生物目标强度特性研究,分别以球形生物、纺锤形鱼类尾明角灯鱼(Ceratoscopelus warmingii)和细长形浮游动物南极大磷虾(Euphausia superba)为例进行仿真计算,并与相适应的经典理论模型进行对比分析。结果表明,对于球形生物,有限元/边界元耦合方法与解析模型的目标强度频响曲线完全吻合;对于纺锤形鱼类,有限元/边界元耦合方法可有效弥补基于模态级数解的形变圆柱体模型在中低频和两端入射时的准确性问题;对于细长形浮游动物,有限元/边界元耦合方法与畸变波玻恩近似模型高度吻合。综上,有限元/边界元耦合方法对多种海洋生物目标强度求解均具有较好的适用性,未来有待进一步结合实验测定进行验证。     

双波段能谷光子晶体的拓扑边界态传输特性

谷赝自旋是构筑光学拓扑绝缘体的重要手段,其性质决定了光波拓扑传输特性。双波段能谷结构为同时调控多波段范围内的光波提供了重要平台。构筑了具有双波段能谷结构的光子晶体超胞,研究了zigzag型和armchair型界面在两个波段内拓扑边界态的传输特性及鲁棒性。结果表明：对于zigzag型界面,对称型界面与反对称型界面的光波传输特性差异显著,对称型界面允许平面波传输而反对称型界面大幅抑制平面波传输;而对于armchair型界面,对称型界面与反对称型界面的光波传输未表现出差异性。对于允许平面波传输的界面,无论是zigzag型还是armchair型界面均对杂质、缺陷、尖锐拐角等具有良好的鲁棒性。双波段拓扑边界态传输特性的系统性研究对于拓展光学拓扑绝缘体的应用空间具有重要价值。     

基于拓扑描述函数的连续体结构拓扑优化方法

提出了一种利用拓扑描述函数(TDF)作为拓扑设计变量求解连续体结构拓扑优化问题 的新方法. 优化问题的目标函数是结构的整体柔顺性，约束条件为对于可利用材料的体积限 制. 这种方法不仅可以消除拓扑优化中经常出现的棋盘格式等数值不稳定现象，而且能够有 效地抑制传统算法处理此类优化问题时所引发的边界扩散效应. 与其它的基于水平集描述函 数的拓扑优化方法相比，所提出的算法不仅无需求解控制水平集函数演化的双曲守恒方 程，而且合理地考虑了目标函数的拓扑导数信息，因而使得算法的计算效率有了显著的提高.     

椎体射频消融的温度场分析

为了探究椎体组织射频消融中的损伤范围,提出利用恒温边界作为热源激励的仿真模型来分析温度场,并且针对温度场构建中血液灌流难以确定的问题,提出一种根据仿真和实验差异用最小二乘法估计血液灌流项的方法,然后设计了动物在体实验并建立相应的有限元模型,针对消融温度为100℃时估计出的血液灌流项,设计了消融温度为95℃的实验进行验证。结果表明,引入血液灌流项后,仿真和实验之间的最大误差由7℃降至5℃,均方根误差由5.3℃降至3.9℃,可以对构建椎体组织射频消融的温度场提供参考。     

不完全边值问题的极小二乘解及其稳定性

定义了一类在部分边界上函数值为0函数的Sobolev空间，并用以讨论部分边界缺乏边值的二阶散度型椭圆型微分方程与其最小范数极小二乘解的稳定性．