openfoam viv 网格变形

[msm]

我在算桥梁涡激振动的时候，最大库郎数=0.9，振幅小的可以运行，振幅大的由于部分边界层网格变形过大导致时间步长无穷小，然后停止运行，有没有大佬指点一下，我下边应该怎么做

[Mina Lee]

你应该放几个动网格的运动图片上来，并且描述一下你的动网格方法。

[msm]

dynamicMeshDict设置
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object dynamicMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dynamicFvMesh dynamicMotionSolverFvMesh;

motionSolverLibs ("libfvMotionSolvers.so");

solver displacementLaplacian;

displacementLaplacianCoeffs {
diffusivity inverseDistance (DECK);
accelerationRelaxation 0.4;
}

[msm]

[msm]

@Mina-Lee 您好，我把动网格设置贴上来了，您可以看一下，网格运行的时候一直挺正常，振幅超过5
cm突然就变形很大，只有这一小块变形大，其他部分正常

[ice_flow]

1.根据现有的结果计算一下局部库朗数，看看是否是突变的速度导致的网格变形。（尝试调整一下时间项的离散方案）
2.用比较稀疏的网格尝试计算，看看能否承受变形（去掉边界层）
3.松弛因子accelerationRelaxation 0.4;可以再调整，网格变形方式也可以尝试调整。
大概是这些方向，如果都不行。那么可以考虑使用重叠网格方案，对于刚体运动是一定可行的。

[Mina Lee]

@msm 不好意思，根据你所提供的网格信息，无法判断是哪里的问题。
你提供的图片没看出来是哪个位置的。

1、可以参考fluent一样，把边界层提取出来，给边界层也指定位移。
2、也可以使用二次逆距离扩散系数试一下，你现在用的是一次逆距离扩散系数。
3、你需要关注一下，是不是因为单个时间步长的位移过大引起的网格问题，如果不是的话，可以调整动网格方法。

[msm]

@ice_flow 谢谢您，我去试试

[msm]

@Mina-Lee 好的，谢谢您的建议