紧急求助，河道水流模拟？？？

donglishi · 发表于 2004-4-10 21:41:12

[这个贴子最后由donglishi在 2004/04/11 07:41pm 第 1 次编辑]

[DISABLELBCODE]

donglishi · 发表于 2004-4-10 21:58:59

结果如图实在是太过均匀了，横向没有太大梯度

tomzijizhang · 发表于 2004-4-12 18:54:48

这个结果也不对。出口（X=500）处合适完全是因为强行应用连续方程的原因。在X=-250到X=250处，流体根本不动。你如果不用填色方式显示（建议用速度矢量图）会发现进口X=500的桔红色速度没有带动下游的流动。怀疑在此边界上差分方程没有表达正确。水总得从上游流到下游呀。总不能堆在进口出。即使堆积，一个正确的遵从物理规律的差分方程组也会让流体往没堆积的地方(下游)流，从而实现数值收敛。

donglishi · 发表于 2004-4-18 09:15:37

那么上游边界该如何处理呢，中心差分需用到边界处的水位、流速情况，假设上游边界流速已知，那水位该如何处理呢？？？

急需帮助，请指教！不胜感激！！

Tomzijizhang · 发表于 2004-4-20 21:16:28

通常边界不包括在计算域内，也就是说那个AX=B的矩阵没有边界元素。边界信息是通过计算域内的最靠近边界的单元/网格带近来的。怎么带内，通过那个单元与边界的接触面带来的。
通常的网格用通常的差分格式编完和，那只是表达了这些内部网格之间的物理关系。可边界约束不守这个物理关系，例如边界规定水深就不能用内部的差分方程表达。
接下来编程的工作是，对靠近边界的单元的那个与边界相邻的面进行特殊处理。例如A51X51+A52X52+A53X53=B5中的51单元是边界单元(差分中指单元)。因为51单元不出现在AX=B的矩阵A中但没说不能出现在B中，所以对最接近边界的内部单元方程变为A52X52+A53X53=B5-A51X51=新B5。这样第一个内部单元X52就有流动了。然后X53就也流起来了。
X可代表水位或流速。水位可能通过压力项导入。

donglishi · 发表于 2004-4-21 12:48:14

谢谢回复！！

		自动登录	找回密码
密码			注册

紧急求助，河道水流模拟？？？