✅作者简介:热爱科研的Matlab仿真开发者,修心和技术同步精进,代码获取、论文复现及科研仿真合作可私信。
🍎个人主页:Matlab科研工作室
🍊个人信条:格物致知。
更多Matlab完整代码及仿真定制内容点击👇
本文描述了利用浸没边界框架对流固耦合系统进行数值模拟。该框架由[3]和[6]提出,并在此基础上进行改进。
流体方程求解:
采用经典的分步法对流体方程进行求解。该方法基于交错网格上的二阶精度中心差分格式[4]。时间积分采用半隐式方案,对流项采用显式三阶龙格-库塔格式,扩散项采用隐式Crank-Nicholson格式进行积分。由半隐式动量方程积分产生的庞大五对角线性方程组,采用近似分解法进行求解。而为了保证质量守恒,用于强制满足连续性方程的泊松方程则采用直接求解器进行求解,以确保达到严格的容差。有关流体相数值方法的更多细节,可参见文献[7]。
浸没边界处理:
遵循文献[9]中的方法,通过在浸没边界上放置拉格朗日标记点,计算满足无滑移边界条件所需的力项。该力项以体积力场的形式表示,然后传递到欧拉节点。拉格朗日标记点上的信息通过移动最小二乘法进行插值[10],该方法被认为可以抑制流体动力载荷在移动界面处的虚假振荡,同时保持空间上的二阶精度。该求解器已在先前的工作[3, 5]中针对各种静止和移动边界问题进行了广泛的验证。
结构方程求解:
刚体运动方程通过状态空间公式转化为两个一阶常微分方程,然后用四阶显式龙格-库塔格式进行积分。由于涉及的相密度比不会影响松散耦合方法的稳定性和精度,因此流体和结构子系统以顺序方式进行求解。Borazjani、Ge和Sotiropoulos[1]的工作中提供了这一假设的数值证据。
边界条件:
在域入口处使用狄利克雷速度边界条件,在域出口处使用辐射外流边界条件。上下边界采用自由滑移条件。相邻拉格朗日标记点之间的相对间距设置为0.5∆x,其中∆x是局部欧拉网格间距,这是在界面条件精度和计算成本之间的数值折衷[6]。
边界层处理:
由于本文采用的浸没边界处理方法,流场在边界上呈现平滑过渡层,其厚度最多为两个欧拉单元,这在de Tullio和Pascazio[3]的数值实验中得到了验证。因此,通过在沿着从边界指向外侧的法线方向创建的探测点处插值场变量来计算粘性和压力载荷[11]。探测点的长度被选择为局部平均单元尺寸。
总结:
本文介绍了利用浸没边界框架对流固耦合系统进行数值模拟的方法。该方法采用经典的分步法和半隐式时间积分方案,并对边界条件和边界层进行了详细的处理。该求解器已在先前的工作中进行了广泛的验证,可以用于解决各种流固耦合问题。
⛳️ 运行结果
🔗 参考文献
🎈 部分理论引用网络文献,若有侵权联系博主删除
👇 关注我领取海量matlab电子书和数学建模资料
🎁 私信完整代码和数据获取及论文数模仿真定制
🌈 各类智能优化算法改进及应用
生产调度、经济调度、装配线调度、充电优化、车间调度、发车优化、水库调度、三维装箱、物流选址、货位优化、公交排班优化、充电桩布局优化、车间布局优化、集装箱船配载优化、水泵组合优化、解医疗资源分配优化、设施布局优化、可视域基站和无人机选址优化、背包问题、 风电场布局、时隙分配优化、 最佳分布式发电单元分配、多阶段管道维修、 工厂-中心-需求点三级选址问题、 应急生活物质配送中心选址、 基站选址、 道路灯柱布置、 枢纽节点部署、 输电线路台风监测装置、 集装箱调度、 机组优化、 投资优化组合、云服务器组合优化、 天线线性阵列分布优化、CVRP问题、VRPPD问题、多中心VRP问题、多层网络的VRP问题、多中心多车型的VRP问题、 动态VRP问题、双层车辆路径规划(2E-VRP)、充电车辆路径规划(EVRP)、油电混合车辆路径规划、混合流水车间问题、 订单拆分调度问题、 公交车的调度排班优化问题、航班摆渡车辆调度问题、选址路径规划问题、港口调度、港口岸桥调度、停机位分配、机场航班调度、泄漏源定位
🌈 机器学习和深度学习时序、回归、分类、聚类和降维
2.1 bp时序、回归预测和分类
2.2 ENS声神经网络时序、回归预测和分类
2.3 SVM/CNN-SVM/LSSVM/RVM支持向量机系列时序、回归预测和分类
2.4 CNN|TCN|GCN卷积神经网络系列时序、回归预测和分类
2.5 ELM/KELM/RELM/DELM极限学习机系列时序、回归预测和分类
2.6 GRU/Bi-GRU/CNN-GRU/CNN-BiGRU门控神经网络时序、回归预测和分类
2.7 ELMAN递归神经网络时序、回归\预测和分类
2.8 LSTM/BiLSTM/CNN-LSTM/CNN-BiLSTM/长短记忆神经网络系列时序、回归预测和分类
2.9 RBF径向基神经网络时序、回归预测和分类