在LS-DYNA中对流固耦合问题的求解提供了两种方法,ALE流固耦合和基于ICFD的流固耦合。本文将对两种方法的各自特点进行分析。
ALE方法的流固耦合
适用范围
ALE求解器主要用于层流。ALE求解器也不是完全的NS求解器,因此不考虑边界层的影响,比如阻力。流体粘性效应仅通过材料模型获得,也就是*mat_null的MU。ale求解器是为了求解高压高速度梯度的短时问题而开发的。该求解器并不适用于求解低压力梯度问题比如声学,也不适用于长时间问题(秒的尺度甚至更长)。时间的限制是显式时间积分造成的,时间步长受限于基于单元尺度和材料声速。在ALE的case中,lsdyna重启动时间步长会进一步被ALE-LAGRANGE耦合罚刚度限制。
能量耗散
ALE求解器的输运固有的存在一定程度的耗散,比如爆炸产生的压力幅值预测偏小,尤其是输运使用一节精度时。使用二阶精度通常可以减少非物理能量耗散,但是会增加额外的计算消耗。细化网格也可以帮助减少能量耗散,但也会再次增加计算消耗。
网格类型限制
lsdyna二次开发初始ALE网格应该使用合适长宽比的六面体单元。退化单元比如四面体、金字塔单元应避免使用,因为他们会导致降低精度,并且在输运过程会造成数值不稳定性。注意在求解过程中,*ale_reference_system可能会影响单元形状。如果单元形状变得不合理,*ale_reference_system需要进行调整以维持单元合理的形状lsdyna软件介绍
基于ALE方法的流固耦合
基于ICFD的流固耦合
ICFD求解器自R7版本开始出现,但是建议使用最新版本,因为最新版本求解性能更好。ICFD求解器必须使用双精度版本。ICFD求解器使用完全隐式时间步长和求解方法,对网格类型没有特殊要求。ICFD求解器是基于有限元法FEM的不可压缩CFD求解器。ls-dyna中基于FEM方法的ICFD求解器是CFD求解器市场一个补充,lsdyna材料参数其他CFD求解器多是基于有限差分法FDM和有限体积法FVM。
基于ICFD的柔性挡板流固耦合