几乎所有有限的方法,利用提出的时间坚持Courant-Friedrichs-Lewy模型法计算一个报号码,把它比作一个Cmax决定稳定性,对于二维:

C = (Vxdt / dx) + (vdt / v) > = Cmax

在哪里C是报号码,V写在下面的方向的速度,dx和dy指定的方向和长度区间dt是在指定的时间步骤。如果你得到数值不稳定,增加数量的网格间距将低报,但然后你可能错过了你想要的动态捕获。如果你减少时间步,它要花很长时间动力学计算。

解决计算问题的一种方法是代码动态时间步:继续计算coraunt数字/ Cmax。如果C > Cmax,然后模型一步降低一半。再次检查,然后重置下一步,这样你就可以有更大的“动态”时代的步伐。我实现了这个多次在其他模型和它工作得很好。

另一种可能性可能是你高纵横比网:如果你的X和你的方向是10公里y1公里,是吗10:1比例非常高。在你的情况下,由于您所使用的模型包是三角形,你想要的角度尽可能急性三角形网格。配给钝角角度将意味着高方面,。

像llmari,我不熟悉水动力模型,所以我不做细节。我的专长是地幔对流和板弯曲模型。我希望这可以帮助。

几乎所有有限的方法,利用提出的时间坚持Courant-Friedrichs-Lewy模型法计算一个报号码,把它比作一个美元C_{马克斯}$决定稳定性,对于二维:

$ $ C = \ \压裂压裂{V_xdt} {dx} + {V_ydt} {dy} \通用电气C_{马克斯}$ $

在哪里美元加元是报号码,V_i美元写在下面的方向的速度,dx美元或dy美元指定的方向和长度区间dt美元是在指定的时间步骤。如果你得到数值不稳定,增加数量的网格间距将低报,但然后你可能错过了你想要的动态捕获。如果你减少时间步,它要花很长时间动力学计算。

解决计算问题的一种方法是代码动态时间步:继续计算美元加元和美元C_{马克斯}$。如果$ C > C_{马克斯}$,然后模型一步降低一半。再次检查,然后重置下一步,这样你就可以有更大的“动态”时代的步伐。我实现了这个多次在其他模型和它工作得很好。

另一种可能性可能是你高纵横比网:如果你的x美元和你的方向是10公里y美元1公里,是吗10:1美元比例非常高。在你的情况下,由于您所使用的模型包是三角形,你想要的角度尽可能急性三角形网格。配给钝角角度将意味着高方面,。

像llmari,我不熟悉水动力模型,所以我不做细节。我的专长是地幔对流和板弯曲模型。我希望这可以帮助。