Spherically Symmetric Solutions to Compressible Hydrodynamic Flow of Liquid Crystals in N Dimensions

The paper is concerned with the system modeling the compressible hydrodynamic flow of liquid crystals with radially symmetric initial data and non-negative initial density in dimension N (N ≥ 2).The au...     

随机Gauss粗糙面上三维导体目标散射差场的随机泛函解析计算方法

提出一种解析的随机泛函方法(SFA),计算导体Gauss粗糙面上三维导体目标的复合电磁散射.推导粗糙面的随机Green函数,用一种新的四路径模型描述面体复合散射机理,用SFA求解双站差场雷达散射截面.以导体球目标为算例,与其他数值计算方法比较后验证了SFA的有效性与准确性,同时讨论了粗糙度、体目标尺寸以及距离粗糙面高度等参量变化对结果的影响,给出复杂形状体目标的双站差场雷达散射截面的空间角分布. 关键词： 随机泛函方法 粗糙面随机Green函数 差场雷达散射截面 面体复合散射     

基于Mie散射理论对砷化镓光子晶体安德森定域化研究

基于Mie散射理论和低浓度近似,对砷化镓作为散射体光子晶体中的安德森定域化参量进行了理论计算,并分析了影响定域化现象的各种因素。结果表明:在散射体体积分数为10%,相对折射率大于3.8时,远红外区50~65 μm范围内出现严格的定域化现象;随着散射体半径的增大,定域化区向长波方向移动,且定域化参量先增大后减小。     

同心圆环阵列模态波束形成器设计

针对贝塞尔函数零点造成的模态波束形成器性能损失的问题,提出了基于均匀同心圆环阵列的模态波束形成器设计方法。从阵元域出发,建立了阵元域和圆谐波域之间的转换关系,推导了圆谐波域阵列信号处理表达式,其中重点推导了延时求和、最小方差无失真响应和多约束波束形成器的设计方法。理论证明了在平面各向同性噪声场中的基于同心环阵的圆谐波域最小方差无失真响应波束形成器等价于相位模式波束形成器。在该结论的基础上,综合考虑相位模式波束形成器在一定频率范围内的稳健性和指向性,提出了一种同心圆环阵列的内环半径优化方法。使用16元双环阵列对本文提出的波束形成算法进行仿真和实验,结果表明:同心圆环阵列能够较好地解决贝塞尔函数零点问题;其中多约束波束形成器能够在多个关联的性能指标之间取得合理的折衷,实现优于相位模式和延时求和波束形成器的声源定位性能。对内环半径优化方法进行仿真和分析结果表明,该方法得到了相较于传统方法更优的阵列结构。     

带形域内裂编译地SH波散射问题的边界元解

采用边界元法研究含裂纹的带形域各向性弹性体，裂纹对ＳＨ波的散射问题，推导出带形域情况下不同边界条件的各种Ｇｒｅｅｎ函数，导出了以裂纹张开位移为未知数的边界积分方程，计算出表面散射场和总位移，算例表明，利用所提供的格林函数和边界元格式解答带形域的散射问题比较方便     

多层地基条带基础动力刚度矩阵的精细积分算法

提出应用精细积分算法计算多层地基的动力刚度问题. 精细积分是计算层状介质中波传播的高效而精确的数值方法. 利用傅里叶积分变换将层状地基的波动方程转换为频率-波数域内的两点边值问题的常微分方程组, 运用精细积分方法求解格林函数, 最后再将得到的频率-波数域内地基表面的动力刚度矩阵转换到频率-空间域内, 进而得到刚性条带基础频率域的动力柔度或刚度矩阵. 所建议的精细积分算法, 可以避免一般传递矩阵计算中的指数溢出问题, 对各种情况有广泛的适应性, 计算稳定, 在高频段可以保障收敛性, 并能达到较高的计算精度.      

瞬态热传导问题的无网格法快速求解算法研究

为了提高基于Galerkin弱积分形式的无网格方法求解瞬态热传导问题的计算效率,提出了两种方案:第一种方案在空间离散上采用基于任意凸多边形节点影响域的无网格形函数,并通过选取适当的节点影响半径因子,使背景网格内的积分点仅对该背景网格内的无网格节点有贡献,从而避免了节点搜索问题,减少了系统刚度矩阵的带宽,且当节点影响半径因子为1.01时,无网格方法的形函数近似具有插值特性;第二种方案在求解线性方程组时,引入质量矩阵集中技术,从而避免了系统方程组的求解.二维矩形区域、二维圆形区域的瞬态热传导数值算例结果表明:在保证计算精度的同时,采用任意多边形节点影响域的无网格方法比传统无网格方法的计算时间至少节省44.09%,采用质量矩阵集中技术的无网格方法比传统无网格方法的计算时间至少节省76.15%,且当节点影响半径因子为1.01时,其本质边界条件的施加和有限元方法一样简单;由于采用质量矩阵集中技术的无网格方法比采用任意多边形节点影响域的无网格方法精度较低,因此如仅从计算效率考虑,对精度要求不是很高(误差在5%以内),建议采用质量矩阵集中技术,如同时考虑计算精度和效率,建议采用多边形节点影响域的技术.     

比例边界坐标插值方法在谱元法中的应用——无穷域Euler方程的数值模拟

将比例边界坐标插值方法引入谱元法, 构成比例边界谱单元, 对无穷域Euler方程进行数值模拟.阐述了比例边界谱单元的基本使用方法以及基于比例边界谱元的Runge-Kutta间断Galerkin方法求解Euler方程的过程;计算了无穷域圆柱和NACA0012翼型绕流问题, 并与已有结果进行了比较, 显示了计算结果的正确性.用基于比例边界谱元的间断Galerkin方法求解无穷域Euler方程时, 最多只需将求解域划分为2个子域, 避免了一般谱方法将求解域划分为9个或者27个子域的麻烦. 比例边界谱单元为无穷域Euler方程的直接求解提供了一个可供参考的方法.     

介观双粗糙弹塑性流体动力润滑界面法向接触刚度模型

弹性流体动力润滑状态通常出现在机械高副零部件的点/线接触部位,如齿轮、轴承和蜗轮蜗杆等.宏观上点/线接触在介观层面表现为两粗糙表面的接触,在微观层面上则又表现为微凸体间的接触.由于在中/重载荷作用下,粗糙表面上的微凸体发生接触后会产生弹塑性/塑性变形,从而使得两粗糙表面的弹流润滑接触转变为弹塑性流体动力润滑接触.此外,界面的接触刚度决定了机械装备的整机刚度.为了精确获得弹性流体动力润滑状态下界面法向接触刚度及其主要影响因素,基于界面的法向接触刚度由固体接触刚度和润滑油膜刚度两部分构成的思想,根据固体弹塑性理论和流体动力学理论,分别对界面间微凸体侧接触及部分膜流体动力润滑进行分析,从微观入手揭示双粗糙表面弹塑性流体动力润滑接触机理,进而建立考虑微凸体侧接触弹塑性变形的流体动力润滑界面法向接触刚度模型.通过仿真分析,揭示了法向载荷、卷吸速度、表面粗糙度及润滑介质特性等因素对润滑界面法向接触刚度的影响规律.研究表明：在相同速度、粗糙度及润滑油黏度的工况下,固体接触刚度和油膜接触刚度均随着法向接触载荷的增加呈非线性增大;在相同载荷、速度及润滑油黏度的工况下,接触表面粗糙度越大,表面形貌对于润滑...