一种适合迭代求解的反馈力浸入边界法

对Goldstein提出的反馈力浸入边界法进行了新的思考，改进了其对力源项的计算，拓展了该浸入边界法的使用范围。传统的反馈力浸入边界法在进行力源项的计算时，含有对速度误差的时间积分项，只能用于含时间项的Navier-Stokes （N-S）方程的求解，且在显式时间推进时有严格的时间步长限制。本文改进的方法则直接通过迭代过程中的速度误差求和来计算力源项，避免了时间相关的参数，使其不仅能适合非定常隐式时间推进，还能与不含时间项的定常N-S方程求解方法结合。为了验证改进方法的可靠性，对二维静止圆柱绕流、静止流体中的振荡圆柱、运动椭圆翼以及三维静止圆球的流场进行了计算，计算结果均与文献结果符合较好，表明本文改进的方法是有效的。得出的结论为：可以直接基于迭代次数进行反馈力源项的计算，改进的反馈力浸入边界法不仅可与非定常N-S方程结合，进行隐式求解，还可以与定常N-S方程结合用于定常流动的模拟，可将改进的方法运用到更多的流动问题当中。     

基于浸入边界法的低雷诺数流固耦合数值模拟(英文)

提出一种基于SIMPLE算法的非定常流固耦合计算方法。流体Navier-Stokes方程空间采用非结构化网格有限体积法离散,时间项采用了欧拉隐式方法。利用浸入式边界方法模拟静止或者运动固体区域,流固界面作用力通过流体体积(VOF)方法进行处理。从而可以用固定网格求解任意复杂区域中的流固耦合作用。本文模拟了低雷诺数静止及振荡圆柱绕流,所得结果与文献中贴体网格计算结果吻合,从而验证了本文方法的合理性和正确性。     

基于浸入式边界方法的二维结冰机翼的数值模拟

应用浸入式边界方法模拟干净/结冰NACA0012翼形绕流问题.采用SIMPLEC方法离散求解二维不可压雷诺平均Navier-Stokes方程, 湍流模型采用一方程Spalart-Allmaras模型.计算结果与以往的实验结果吻合得很好, 证明了浸入式边界方法和所发展的程序的可靠性.      

浸入式贴体网格边界方法

提出一种基于重叠网格的浸入式边界方法,用以模拟复杂的三维黏性跨声速流动。该方法采用一套固定的笛卡儿正交网格,用以主流的求解;采用一套可以移动的贴体等距面网格,用以拟合或者离散物面的作用力,通过空间插值,实现两套网格重叠部分的信息传递。分析了浸入式贴体网格边界方法的优势,介绍了贡献单元的寻找策略和物理通量的插值方法。流场求解采用Spalart-Allmaras带湍流模型的Navier-Stokes方程组,其中对流项采用流通矢量分裂和5阶WENO(weighted essentially non-oscillatory)-Z格式离散,黏性项采用6阶中心差分格式离散,时间项采用龙格库塔显式格式离散。数值验证算例表明:该方法具备高于4阶的空间求解精度,并适用于刚体动网格非定常流场模拟,且无需更新网格形状。等距面贴体网格生成过程简单,避免了繁琐的人工调节过程,与笛卡儿网格结合,可提供足够的壁面附近网格密度,同时有效减少了网格总量需求。     

浸入边界法模拟小圆柱对主圆柱涡脱落的抑制

使用浸入边界法研究了小圆柱对主圆柱涡脱落的抑制.方法中使用非贴体笛儿尔网格,易于处理包含复杂边界的流动问题.采用离散附加力直接加入边界条件方法对虚单元进行重构,使边界条件在浸入边界上精确满足.使用隐式分步法解二维非定常不可压Navier-Stokes方程,通过速度和压力解耦提高计算效率.数值模拟单圆柱绕流及不同位置小圆柱和主圆柱的流动干扰,通过分析流场涡结构和升、阻力系数,得到小圆柱对主圆柱涡脱落的延迟和抑制作用.计算结果与已有实验结论和数值结果对比,计算误差不超过5%,说明浸入边界法可以简单有效地处理圆柱涡脱落抑制这类流动干扰问题.     

基于浸入式边界方法的运动物体声辐射数值模拟

为了在声学计算中简化网格的处理,从而发展出适用于运动复杂边界问题模拟的计算格式,基于浸入式边界方法,提出了一种适用于无黏流动、无穿透边界条件的处理方法,并对一些简单的运动物体声辐射问题进行了数值模拟.数值结果与解析解得到了令人满意的对比,初步验证了所提出的研究方案的可行性.     

基于深度神经网络的含运动边界非定常流场预测方法研究

针对飞行器设计过程中对流固耦合系统快速预测的需求,探索基于数据驱动的非定常流场建模策略,缩短流场演化求解耗时,从而加快流固耦合系统模拟速度。流固耦合系统中流场演化求解部分等价于含运动边界的非定常流场演化。本文提出了一种基于神经网络的流场预测模型,来学习并预测运动边界附近的非定常流场演化过程。此神经网络可以基于当前流场状态及边界运动信息预测出下一个时刻的流场状态。通过不同振动频率及振幅下的运动圆柱绕流问题,测试了本文提出的神经网络模型的预测精度及泛化能力。神经网络的预测结果和计算流体力学仿真结果中流场结构吻合度较高,且通过对预测的流场数据中边界上的压力积分得到的气动力也具有较高的精度。测试结果展示了该神经网络模型具有良好的预测性能,因此该方法可以用于快速、准确获得运动边界周围非定常流场状态。     

一种改进的浸入运动边界算法

浸入运动边界是数值处理流固耦合作用的主要方法之一，其原始权重函数根据经验确定且形式固定，难以准确地定量描述复杂的流固现象。本文提出了一种新的权重函数对格子Boltzmann方法框架下的浸入运动边界进行修正，通过假定零固含率处的权重因子多阶导数为0，引入可调参数b增强权重函数对固含率的关联，实现了对固体边界渗透性能的改进。对静止圆柱绕流、Taylor-Couette流和振动圆柱绕流3个算例开展数值实验，结果表明：随着雷诺数的增加，原始权重函数预测的升阻力系数均高于文献结果，严重低估了边界处流体组分的作用；当新权重函数中的参数b为3时，改进的浸入运动边界效果最佳；新权重函数增强了固体边界的渗透率，提高了流场的预测精度，使平均误差降低了38.5%；最后，运用新权重函数成功地复现了振动圆柱绕流的锁频现象，表明改进的浸入运动边界在处理复杂的流固耦合问题方面具有广阔的应用前景。     

基于机器学习的非定常流场网格自适应

现有针对非定常流场数值模拟的网格自适应方法,通常每隔一段时间步就进行一次网格调整,增加了计算复杂度和精度损失。针对这一问题,本文基于间断伽辽金（discontinuity Galerkin, DG）有限元法提出了结合BPNN(backpropagation neural network)和MMPDE(moving mesh partial differential equation)的非定常流场网格自适应方法。该方法首先采用DG有限元法对Navier-Stokes方程进行非定常计算,得到统计意义上的网格间断量;然后以初始网格和间断量训练BPNN回归模型,用于预测任意位置节点的间断量;接着使用MMPDE变分法移动网格节点,使其符合统计意义的间断量分布;最终通过Laplacian网格平滑法保证网格单元质量。圆柱绕流非定常流场算例的验证结果表明,该方法能够在不改变网格拓扑结构和不增加节点数的情况下完成一次性网格自适应,显著提高了非定常流场数值计算的精度和效率。     

基于IBVCM的二维非定常动边界数值模拟研究

主要研究了影响非定常动边界问题的一些基本参数并分析了其与推进力的关系。数值模拟采用格子波尔兹曼方法求解流场信息,并用浸入边界速度修正方法进行边界处速度的修正。研究对象为二维情况下摆动的鱼,影响推进力的参数有鱼体的形状、频率、波长、雷诺数。通过对各个参数的一系列模拟,得出其对推进力的不同影响关系。本文总共选取五种形状的鱼体,其运动模式为鱼体中线以行波的方式做横向摆动。最后,从涡的角度分析鱼的尾尖在不同位置时鱼体表面涡的分布情况,得到旋涡变化与鱼身运动的相互关系。通过研究发现,不同形状的鱼推进力不同,在高频、增大波长、增大雷诺数时,推进力增大。     

基于指数插值的浸没边界法在可压缩流模拟中的应用研究

由于浸没边界法在处理复杂外形和动态问题时具有极强的灵活性和便利性,因而受到广泛的关注。在当前浸没边界法的研究和应用中,主要工作集中在关于中/低雷诺数和不可压流动的模拟。本文引入 Choi 等提出的指数插值法对目标点速度进行插值,并根据多项式插值和 Crocce-Busemann 关系式等确定目标点压力和温度,并借鉴直角网格法中的虚拟网格法计算物面内的流动变量,建立了适用于可压缩流动的浸没边界法。通过亚、跨声速段翼型绕流和超声速圆柱绕流算例对方法进行验证,结果显示在 k 值取值合适时模拟结果与实验数据较为吻合。通过对相关算例网格收敛性分析,证明了本文发展的浸没边界法具有较高的精度。文中各算例的模拟结果表明,此方法能准确捕捉亚、跨、超声速绕流问题中的各类流动结构,适用于可压缩流动的模拟。     

非紧致气动噪声半空间传播基本解的边界积分方法

为了获得低马赫数流动诱发的非紧致气动噪声在半空间内传播的基本解,结合复等效源方法和边界元方法建立了半空间精确格林函数的边界积分方程,当半空间边界为刚度型阻抗边界时可避免奇异性积分.同时基于等效源方法提出一种半空间二维非紧致圆柱声散射模型,推导了静止介质中声散射基本解的理论表达式.对静止介质中的二维圆柱声散射,数值解在研...     

含孔边裂纹板的弯曲断裂计算

 采用复变函数理论和边界配置方法 ,对含孔边裂纹板的弯曲断裂进行了分析计算。首先假设挠度的复变函数式 ,进而可以求板的内力。它们能满足一系列的基本方程和支配条件 ,仅板的边界条件需要考虑 ,并且可用边界配置法和最小二乘法近似满足。对孔双边裂纹问题进行了应力强度因子计算。数值算例表明 ,本文方法精度较高 ,计算量小 ,是一种有效的半解析、半数值计算方法。     

含运动边界二维多体干扰流场的数值模拟

本文采用二阶精度Godonov扩展方法求解非定常：Euler方程，以基于格心的自适应非结构网格的显式有限体积法为基础，空间和时间都是二阶精度，使用HLLC近似黎曼解的方法计算网格单元边界处的守恒量通量，对作俯仰振动的NACA0012翼型绕流及激波管中运动、静止圆柱的流场进行了数值模拟，并对数值模拟的结果进行了分析，为进一步对含多体分离流场进行数值模拟打下了基础。     

自适应笛卡尔网格超声速黏性流动数值模拟

复杂外形/流场的高质量网格的生成往往需要占用大量人力资源,而自适应笛卡尔网格方法能够自动化生成高质量网格,具有很好的工程实用价值和应用前景。基于笛卡尔网格方法,采用叉树数据结构进行数据的存储和访问,分别从几何特征和流场解特征出发进行网格的自适应加密和粗化,发展了一种二维情况下自动、高效的自适应笛卡尔网格生成方法。从浸入边界方法出发,结合虚拟镜像对称方法和曲率修正技术进行黏性物面边界条件的处理,同时建立了多值点问题的处理技术,发展了一种在笛卡尔网格下可有效模拟黏性物面边界条件的方法。针对自适应笛卡尔网格非均匀的特点,发展了悬挂网格的处理方法,并构建了适用于自适应笛卡尔网格的黏性数值求解器。通过典型算例的考核,验证了所发展的自适应笛卡尔网格生成技术和构建的数值求解器具有较高的精度和可靠性。     

考虑地面影响的半空间气动噪声数值计算方法

含非紧致边界的半空间气动噪声问题普遍存在于航空领域,在数值计算过程中不仅要考虑非紧致边界的散射效应,还需要考虑地面对声波的反射。基于Lighthill声模拟理论和广义格林函数方法,引入边界声压积分项和可压缩修正项可以考虑非紧致边界以及流体可压缩性对声波散射的影响。然而在半空间声场的计算过程中,包含沿无限大地面的积分项,为了克服沿地面积分所带来的巨大工作量,根据无穷大平板对声波的反射特性,推导得出频域下基于半空间格林函数的声学积分方程。同时开展了半空间中应力源噪声与圆柱绕流噪声的数值预测与分析,验证了该方法的准确性与高效性。