光滑粒子流体动力学方法及其在气动声学中的应用

来源 :天津大学 | 被引量 : 0次 | 上传用户:tatimess
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
在人体发音过程中,声道的形状会实时变化,而声道的形状决定着声音的音色以及内容,因此,当声道形状改变时声音也会发生相应的改变。此外,当声道中介质的流动速度较高时,介质的流动也可对声波的传播产生较大影响。因此,为了更加准确、真实地模拟人体发音过程中声波在声道中的传播,针对移动边界条件下流动介质中的声波传播问题的研究尤为重要。数值模拟是研究声学问题的重要手段。目前,应用在计算声学上的数值方法依据离散形式,可以分为网格法与无网格法。网格法主要包括有限差分法、有限元法、边界元法等。网格方法都必须按照一定规则对计算域进行网格划分并根据网格进行计算。本文所采用的光滑粒子动力学(Smoothed Particle Hydrodynamics,SPH)方法是一种无网格数值计算方法,它克服了网格方法中网格生成、网格畸变和网格移动引起的各种问题。与拉格朗日力学有机结合,该方法被认为是一种分析各类流体动力学问题的极具发展前景的数值计算方法。针对所关心的物理现象,即人体发音过程中声波在声道中的传播过程,本文首先建立了相应的数学模型。为了降低数值求解的复杂度对其进行了合理简化,得到了欧拉体系下的气动声学欧拉方程。然后结合物质导数的定义,推导了拉格朗日体系下声波传播的控制方程。分别利用SPH方法以及二阶精度的龙格-库塔(Runge-Kutta,RK-2)方法对控制方程中的空间导数和时间导数进行数值离散,建立了控制方程的数值离散格式。为了保证核函数计算的准确性,本文利用虚粒子方法处理因边界截断带来的边界信息的缺失问题,并简要介绍了在实际计算中虚粒子实施的过程。在计算过程中本文采用了效率较高的链表搜索算法进行邻域内邻近粒子的搜索,相比于通常采用的全局搜索算法大大提升了搜索效率。为了验证SPH方法在求解流动问题时的正确性及有效性,本文进行了溃坝问题和对称弯曲管问题的数值实验。通过与其它文章中的结果进行对比,验证了本文结果的正确性。为了应用于实际的气动声学问题,本文利用SPH方法成功求解了拉格朗日体系下的欧拉方程。通过对比静止流体中声波传播问题的SPH解和有限时域差分(FDTD)解,验证了SPH方法在声学计算中的准确性和可靠性。通过成功模拟一维和二维流动介质中的声波传播问题,阐明了利用SPH方法求解复杂气动声学问题的可行性。
其他文献
本文利用HPLC色素分析技术,于2013年5月~2014年6月调查了大亚湾大鹏澳海域微表层和次表层水体浮游植物光合色素的时空分布以及粒径特征,揭示了水产养殖对浮游植物群落结构和
极区夏季中层顶的尘埃等离子体有明显的分层现象,且回波强度与电子温度和观测雷达频率密切相关。本论文中,初步分析了极区夏季中层进行主动调制实验中产生的回波现象,利用电
本文主要探讨了Sturmian字中一类特别字中心字的一些性质情况,中心字是具有互素周期p,q,且其长度为p+q-2的字.我们在认真研读Luca与Carpi(2005)的理论基础之后,由其中对于字
本文引入了稳定的二阶全局极小值及tilt稳定的全局极小值概念.在f是正常的下半连续函数的假设下,利用变分分析的技巧和方法,证明了二阶全局极小值的稳定性、全局极小值的tilt
PIH1D1是一个在真菌与哺乳动物中进化保守的、隶属于R2TP复合物(Rvb1/Rvb2/Tah1/Pih)的支架蛋白,其在酵母中名为Pih1。伴侣样蛋白复合物R2TP被认为涉及诸如小核仁蛋白复合物
为了了解巢湖沉积物中溞属枝角类休眠卵的时空分布,2014年8月对巢湖沉积物进行了采样。根据巢湖的形态特征和营养状况,将其划分为三个采样湖区,即东湖区、中湖区和西湖区。经
本文通过基于密度泛函理论的第一性原理计算方法,研究了过渡金属原子吸附的单层过渡金属二硫属化物MX2(M=Mo,W;X=S,Se Te)以及过渡金属原子掺杂的单层反铁磁MnPX3(X=S,Se)中
纠缠在量子信息处理中起着重要作用,量子信息处理主要包括量子计算和量子通信。纠缠态的制备是量子纠缠研究的一个重要方面。某些特殊的量子纠缠态的制备方法已经相当成熟,例
内域介质波动数值模拟是大型结构-基础相互作用分析中的一个重要方面。本文基于有限元-有限差分法及ABAQUS/Explicit积分格式,研究离散网格中波动的传播规律,提出计算精度的
远离热力学平衡态的时空斑图广泛存在于物理、生物、生态和化学系统中。斑图动力学是从动力学角度对这类斑图的产生机理和规律进行研究的一门横向学科,是非线性科学领域的重