高精度HWENO格式与浸入边界法求解可压缩流问题

在物面边界处采用浸入边界法并构造三阶有限体积HWENO(Hermite Weighted Essentially Non-oscillatory)格式,可在较简单的笛卡尔网格上有效处理上述带复杂物面边界的可压缩流动问题。经典的定常和非定常数值算例验证了该方法的有效性。     

基函数法在黏性可压缩流动中的研究与应用

研究基函数法在二维和三维轴对称情形下可压缩黏性流动中的应用, 并利用此法数值地计算了二维圆柱的超音速黏性绕流和三维球头的轴对称超音速黏性绕流两个算例。取一阶三角函数为基函数, 构造出导数的中心格式和迎风格式。对于可压缩N-S 方程中的对流项, 采用通量分裂法及中心格式与迎风格式相结合的技术;对于可压缩N-S方程中的黏性项, 则采用中心格式进行处理。由此, 构造出了数值求解黏性可压缩流动的一阶三角函数类型的基函数格式。两个算例的数值计算结果表明, 基函数法不仅在处理无黏可压缩流动时, 是一种高精度、高分辨率的新型计算方法, 而且在处理黏性可压缩流动问题时也一样行之有效。     

纽结补中的不可压缩边界不可压缩曲面

利用扭转交叉数给出了几乎交错纽结补中不可压缩边界不可压缩曲面下的欧拉示性数估值，进而讨论了Ｃａｂｌｉｎｇ猜想。     

Ghost Fluid方法在可压缩多介质流动中的应用

在对可压多介质流动进行数值模拟时,界面两边流体状态方程不同会给计算带来很大麻烦,该文采用Ghost Fluid方法(GFM)成功地解决了这一问题.利用间断有限元Galerkin方法和Level Set方法,对2种不同比热比介质激波管问题进行了模拟,其数值结果和精确解非常吻合;对气体/水激波管、空气/氦气泡相互作用、空气/R22相互作用等问题进行了模拟,其结果与相关文献的结果一致,表明利用该方法进行可压缩多介质流动数值模拟是可行的.     

基于嵌入式迭代的高精度浸入边界法

针对传统浸入边界法施加边界条件精度不高的情况,将附加体积力的求解内嵌到压强泊松方程的迭代当中,提出了基于嵌入式迭代的浸入边界法.将所提出的算法应用到圆柱绕流数值模拟中,并与他人的实验值和数值解进行了对比,验证了所提算法具有较高的计算精度.分析了插值函数是导致浸入边界法高估圆柱阻力系数的原因,并对阻力系数进行了修正,修正值与采用贴体网格得到的结果非常吻合.比较了本文算法和传统算法的圆柱迎流侧驻点附近的速度矢量和圆柱周围的流线,证实了本文算法在施加边界条件方面的优越性.比较了两种算法的计算效率,迭代计算附加体积力所导致的额外耗时不显著.     

不可压缩流体周期性流动格子波尔兹曼的边界处理

提出了一种新的周期性边界处理方法,用于处理格子波尔兹曼(lattice boltzmann)方法中压力驱动不可压缩流体周期性流动的边界条件。当压力驱动的不可压缩流体充分发展时,可以在周期性截面上直接用压力梯度而不是采用当量体积力的常规周期性边界处理方法,来确定周期性流动截面上的粒子分布密度。应用所提出的压力驱动不可压缩流体周期性流动的边界处理方法,对二维压力驱动下充分发展的泊肃叶流(poiseuille flow)和周期性方柱扰流进行了数值模拟。结果表明:该周期边界条件处理方法能够很好地保持系统周期性,在其宏观的周期性边界截面上产生了连续的压力分布,真正地实现了直接压力驱动不可压缩流体周期性流动条件下的所有特征。     

LU-AUSM~ 混合格式在三维可压缩流动数值分析中的应用

将一维AUSM 格式扩展到三维任意曲线坐标 .为了获得高阶精度 ,采用了三阶MUSCL格式 .将AUSM 格式与LU SGS格式结合 ,对凸包通道跨音速无粘流动、平面叶栅跨音速和超音速无粘流动 ,以及喷管超音速粘性流动进行了数值模拟 .计算结果与文献计算结果和试验数据相符很好 .     

具有运动边界不可压缩黏性流动的CBS有限元解法

基于特征分裂(CBS)算法,提出了一种适用于具有运动边界的不可压缩黏性流动问题的有限元解法,给出了不可压缩黏性流动问题的控制方程,并推导了其无量纲形式.将传统CBS算法与动网格技术相结合,提出了一种动边界流动问题的求解流程,并对所采用的基于弹簧近似的动网格方法进行了说明.采用该方法对静水和均匀来流中的简谐振动圆柱绕流问题进行了求解,计算所得圆柱受力、圆柱表面压力和涡量分布以及流场速度分布等结果均与已有数值或实验结果相吻合,初步验证了所提方法的可行性.     

边界单元法在断裂力学中的应用

本文从弹性力基本方程出发,叙述了用边界单元法求解断裂力学问题的基本原理和数值计算方法.在微机上用常量单元计算了承拉双边裂纹板和中心裂纹板的应力强度因子K_1和裂纹张开位移,论证了边界单元法用于断裂问题的有效性.     

刚体-流体耦合运动的浸入边界法研究

为避免复杂贴体网格的生成和动、静界面滑移网格技术插值带来的误差,提出采用浸入边界法模拟流体刚体的耦合运动.借助求解不可压缩N-S方程组的分步投影方法的思想,来求解基于浸入边界法的耦合系统方程.其中刚体边界离散点的作用力密度通过强制满足刚体边界的无滑移条件（位移和速度一致）导出,结合刚体定轴旋转的动量矩定理,采用Newmark时间积分方法,实现对刚体运动特性的预测.同时,通过δ光滑函数,将刚体边界离散点的作用力密度转换到流体的欧拉网格点上,实现对流场速度的修正.在空间离散上,对流项采用Quick迎风格式,扩散项采用中心差分格式,采用2阶显式Adams-Bashforth法离散时间项.以雷诺数为300的圆柱绕流为基准数值算例,验证数值计算结果的可靠性.      

基于浸入边界法的二维流场数值模拟

采用浸入边界方法进行了圆柱绕流的数值模拟.浸入边界方法是通过在Navier-Stokes方程中加入一项体积力,以满足物面边界的无滑移条件.整个流场在笛卡尔网格上求解,物面边界用均匀分布的拉格朗日点表示.拉格朗日点上的力通过离散了的Navier-Stokes方程求解后再分布到周围的网格点上.在有限差分法的基础上利用浸入边界方法分别计算了Re=40和Re=100时圆柱绕流,并将升力系数,阻力系数,斯特鲁哈尔数等结果与文献中的数值和试验结果进行了比较.     

边界元素法和矩量法在终端开路同轴线中的应用

本文应用边界元素法和矩最法研究无损测量中终端开路同轴线作为传感器时的边缘电容及场分布，把边界元素法用在了轴对称的电磁场问题及多介质区域。文中给出了两种方法求解该问题的主要步骤和公式，并对计算细节作了说明，尤其是对奇点的处理。对边缘电容Ｃ（ε）和介电常数ε之间所建立的线性模型进行了全面讨论，给出了齐全的数据和结果，并同有关文献的结果及实验值进行了比较，验证了方法的可靠性及精确性，为实验工作提供了理论     

沉浸边界法在湍流燃烧中的应用

沉浸边界法可以在直角坐标系中处理各种形状的边界问题．将此用于三维的圆柱火焰绕流问题的研究,其中,湍流模型采用大涡模拟（LES）,化学反应速率采用EBU漩涡破碎模型．计算结果与实验结果比较相符,基本反映了湍流火焰的绕流特征,说明该方法用于讨论湍流火焰的绕流问题是有效的．     

随机边界远法及其应用

本文以边界无法为基础，导出结构荷载效应特征函数算式，提出一种新的可靠度计算方法－随机边界元法，从算例表明，应用西文方法分析和求解复杂工程结构的可靠度时，具有程序实施简便，计算度高等优点。     

边界单元法的改进及其在求解三维问题中的应用

本文提出虚边界单元法和点载荷法,有效地提高了边界上和边界附近域解的精度,并将这两种方法推广到求解三维问题,使问题的求解得到了简化。     

LU—AUSM^＋混合格式在三维可压缩流动数值分析中的应用

将一维ＡＵＳＭ＾＋格式扩展到三维任意曲线坐标。为了获得高阶精度,采用了三阶ＭＵＳＣＩ格式。将ＡＵＳＭ＾＋格式与ＬＵ－ＳＧＳ格式结合,对凸包通道跨音速无粘流动、平面叶栅跨音速和超音速无粘流动,以及喷管超音速粘性流动进行了数值模拟。计算结果与文献计算结果和试验数据相符很好。     

基于浸入边界法的双层柔性蒙皮减阻性能研究

为了降低水下航行器受到的阻力,本文根据柔性表面减阻技术以及双层隔振技术提出了一种双层柔性蒙皮结构.利用基于浸入边界法的开源软件包IB2d对所提出的双层柔性蒙皮结构进行了数值模拟,分析了各参数对减阻性能的影响.分析结果表明:双层柔性蒙皮在特定间距范围内具有减阻效果,并且在一定的范围内增大刚度比会减小双层柔性蒙皮近壁面流速...     

边界元法及其在岩石力学中的应用

边界元法是晚近十数年来迅速发展起来的一种数值计算方法,特别适用于解决岩石力学方面的问题.在平面问题Kelvin解的基础上,本文讨论了应力不连续法和直接边界积分法,并对应于应力不连续法给出了位移时不连续法的解。最后给出了几个在岩石力学中应用实例。     

边界单元法在地基动力分析中的应用

本文系采用边界单元法对半无限地基进行三维动力分析.分析表明,在半无限体的三维动力分析中,该法远比有限单元法优越得多.文中有关基础与地基间相互作用的动力分析问题的讨论,仅限于三维刚性基础和均质线弹性(或粘弹性)土体.本文将动力问题用边界元法表示成时间域中的公式,并以逐步积分法求得相应的动力反应.最后,将圆形和矩形刚性基础的计算结果与理论值相比,两者十分吻合.     

边界拟合坐标变换在有限单元法中的应用

用有限分析解实现边界拟合坐标变换，构造四节点超参元来代替四节点等参元． 拉普拉斯方程和泊松方程的算例表明，方法简便且有较高精度。