并行间断有限元算法求解Navier-Stokes方程

间断Galerkin有限元方法非常适合在非结构网格上高精度求解Navier-Stokes方程,然而其十分耗费计算资源.为了提高计算效率,提出了高效的MIMD并行算法.采用隐式时间离散GMRES+LU SGS格式,结合多重网格方法,当地时间步长加速算法收敛.为了保证各处理器间负载平衡,采用区域分解二级图方法划分网格,实现内存合理分配,数据只在相邻处理器间传递.数值模拟了RAE2822翼型和M6黏性绕流,加速比基本呈线性变化且接近理想值.结果表明了该算法能有效减少计算时间、合理分配内存,具有较高的加速比和并行效率,适合于MIMD粗粒度科学计算.     

浅海内孤立波动态传播过程中声波模态强度起伏规律

内孤立波是一种常见于浅海海域的非线性内波,具有振幅大、周期短和流速强等特点,它通过扰动水体中的温盐分布使声速剖面产生明显的距离依赖性,进而影响水下声传播特性.内波自生成后通常以1 m/s量级的速度传播,运动的内波使声传播路径上的声波模态能量在空间和时间上剧烈起伏.本文定义模态强度为模态系数模值(模态幅度)的平方,并用其衡量各阶模态所含声能量的大小.文中基于耦合简正波理论,推导了内波运动时声波模态强度起伏的表达式,将模态强度表征为振荡项和趋势项的线性叠加.以往的工作大多局限于单独从时域或频域研究内波运动时声波模态强度的时变规律,本文则结合短时傅里叶变换在时频平面上揭示了模态强度的起伏机理.理论推导和数值仿真均表明内孤立波使各阶声波模态之间发生能量交换,即模态耦合.内波的动态传播进一步引起模态干涉,这种干涉效应表现为模态强度中的振荡项并使模态强度随时间快速起伏.受模态剥离效应(不同阶模态之间衰减系数的差异)的影响,趋势项的幅度随时间不断变化,进而对模态干涉引起的振荡叠加了时变的偏置.模态强度的整体走势和振荡项中各频率分量振幅的时变特征均与模态衰减密切相关.同时,本文使用深度积分声强作为总接...     

通气空泡生成和溃灭特性试验研究

为了研究通气流量对超空泡生成和溃灭的影响,在西北工业大学高速水洞实验室进行了系列通气超空泡生成和渍灭特征的试验.结果表明:空泡生成过程中,空泡长度和直径的变化特性与空泡溃灭时的变化特性相反;通气流量对空泡生成和溃灭均有较大影响,数据显示通气量对空泡生成的影响大于对空泡溃灭的影响.此研究为进一步研究超空泡的快速稳定生成提供了实验基础.     

宽带密排平面阵声辐射分析

应用声有限元-边界元法计算了宽带基阵声场分布.建立了平面9元阵有限元-边界元模型,给出了阵元之间互阻抗数值计算方法,分析了阵元间互辐射对声性能影响的变化规律,计算了9元阵3个谐振频率点处指向性和波束宽度.利用3谐振换能器研制出一种宽带密排平面9元阵,并在消声水池中对其输入阻抗和指向性进行了实验测量.研究结果表明,按阵元...     

基于Shamma模型的车辆噪声听觉谱特性与目标分类

为了满足现代交通对汽车类型自动识别的需求,针对等响度车辆噪声提出基于Shamma模型的听觉谱特性识别特征。首先分析了Shamma模型中听觉早期阶段和听皮层多尺度分析阶段的转换原理,给出了中枢听觉模型中STRF模型及纹波密度-速率图的计算方法。利用三类车辆空转状态下噪声的幅频、听觉谱和密度-速率特性进行分类识别,发现密度-速率特征具有识别率高、计算量小的特点,更适宜用于作为目标自动识别特征。     

基于格子Boltzmann方法的二维气泡群熟化过程模拟

Ostwald熟化(ripening)是指局部热力学平衡状态下,颗粒/液滴/气泡系统为了减小界面能而自发地进行颗粒群尺度分布调整的过程,具有重要研究价值.针对目前数值模拟研究不充分的现状,本文采用格子Boltzmann方法,对相变速率主控的二维蒸气泡系统演化开展了数值模拟研究.模拟结果与本文推导的二维气泡群演化标度律符合较好,证实了格子Boltzmann方法对复杂相变-物质输运过程捕捉的准确性.研究同时表明,蒸气泡系统演化过程中物质输运为液相压力不平衡所驱动,并且在小气泡“溃灭”过程中水动力学作用会影响气泡群半径分布函数的局部细节;气-液状态方程参数对熟化过程的影响效果分析显示,气液两相比内能差是驱动相变的核心要素,此差异越大相变速率越快,该结论进一步诠释了化学势驱动熟化过程的物理图像.     

浅海负跃层中利用互相关输出峰值迁移曲线的声源深度判别

针对浅海负跃层波导中的声源深度判别问题,通过研究不同深度声源信号到达接收阵的简正波结构特征,提出了一种利用互相关输出峰值迁移曲线的声源深度判别方法.该方法利用浅海负跃层中主导简正波类型由声源深度决定这一先验信息,根据峰值迁移曲线中互相关时延在0时刻的位置进行声源深度判别.首先,利用垂直线列阵波束输出与单阵元接收信号做互相关处理,得到互相关时延—俯仰角输出二维图.然后,在图中提取不同时延上角度维峰值,根据时延小于等于0、大于0的情况对所提取的峰值点进行分段线性拟合,得到互相关输出峰值迁移曲线.最后,以曲线位于互相关时延为0时的位置进行声源深度判断.理论分析表明,曲线受波导中声速随深度变化的影响较小,互相关时延0时刻位置的判别阈值由低阶简正波决定,区分声源深度的区间由跃层的位置和厚度决定,且跃层上下声速差越大越有利于声源深度判别.仿真实验和海试数据分析结果表明,该方法可有效进行声源深度判别,并且无需模态分离、声源相对接收阵移动、精确声速剖面测量等前提条件.     

纳米通道内液态微流动密度分布特性数值模拟研究

微通道内流动因表面积/体积比值极大, 造成许多微尺度效应, 进而使微通道内出现完全不同于宏观流动的流体密度分布特性. 本文以纳米通道内液态Poiseuille流为对象, 采用非平衡分子动力学模拟方法研究了流体原子间相互作用强度ε_LL, 流体原子间平衡距离σ_LL以及壁面原子与流体原子间平衡距离σ_LS对通道内流体密度分布的影响规律. 数值模拟中, 统计系综取微正则系综, 势能函数选用LJ/126模型, 壁面设为Rigid-atom壁面, 温度校正使用速度定标法, 牛顿运动方程的求解则采用Verlet算法. 模拟结果表明, 随ε_LL的减弱, 近壁面区密度分布的振荡幅度则逐渐增大; 而σ_LL 则同时影响流体原子的存在形态和密度分布, 较大的σ_LL 会造成流体原子在整个通道内呈现面心立方结构的类似固体排列, 较小的σ_LL会使得流体原子呈现不断变化的 "团簇" 结构; 随σ_LS的变大, 近壁面区流体密度振荡幅度增大, 且流体密度分布起点离壁面越远. 另外, 本文还从近壁面区流体原子的 "俘获-逃逸" 行为角度, 初步解释了原子间相互作用强度对密度分布的影响规律. 关键词： 纳米通道 微流动 密度分布 分子动力学     

空间飞行器轨道估计与误差分析

卫星探测是当今空间飞行器观测与轨道估计的重要方式.针对空间飞行器轨道估计及卫星无源探测的误差分析等问题,在双星观测条件下,使用样条插值及拟合的方法对包含随机误差和系统误差的初始数据进行光滑处理,然后基于立体几何和逐点交汇的思想建立了运动方程参数模型,采用最小二乘法估计参数,得到了空间飞行器的估计轨道,并使用LMF法进行了系统误差分析.针对单星观测的情况,一方面,认为飞行器轨道上的点共面且均在过地心的平面内,另一方面,通过建立飞行器轨道模版数据库与观测相对比,了解识别是何种飞行器.添加了这些约束和先验信息后,可以得到大致的轨道估计.