几乎所有有限的方法,利用提出的时间坚持Courant-Friedrichs-Lewy模型法计算一个报号码,把它比作一个Cmax美元C_{马克斯}$决定稳定性,对于二维:
C = (Vxdt / dx) + (vdt / v) > = Cmax$ $ C = \ \压裂压裂{V_xdt} {dx} + {V_ydt} {dy} \通用电气C_{马克斯}$ $
在哪里C美元加元是报号码,VV_i美元写在下面的方向的速度,dx和dydx美元或dy美元指定的方向和长度区间dtdt美元是在指定的时间步骤。如果你得到数值不稳定,增加数量的网格间距将低报,但然后你可能错过了你想要的动态捕获。如果你减少时间步,它要花很长时间动力学计算。
解决计算问题的一种方法是代码动态时间步:继续计算coraunt数字/ Cmax美元加元和美元C_{马克斯}$。如果C > Cmax$ C > C_{马克斯}$,然后模型一步降低一半。再次检查,然后重置下一步,这样你就可以有更大的“动态”时代的步伐。我实现了这个多次在其他模型和它工作得很好。
另一种可能性可能是你高纵横比网:如果你的Xx美元和你的方向是10公里yy美元1公里,是吗10:110:1美元比例非常高。在你的情况下,由于您所使用的模型包是三角形,你想要的角度尽可能急性三角形网格。配给钝角角度将意味着高方面,。
像llmari,我不熟悉水动力模型,所以我不做细节。我的专长是地幔对流和板弯曲模型。我希望这可以帮助。