二维浅水方程边界条件的困惑

lijixuan · 发表于 2005-3-14 11:03:51

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？注册

x

小弟正在做用有限体积法在三角网格上求解二维浅水方程的程序，程序已经写好了，但是边界条件出现了问题，想向各位大侠请教。边界分为三种情况：给定单宽流量、给定水位和给定水位流量关系。比较理想的方式是：上游给定单宽流量，下游给定水位。上游给定单宽流量时是把Qb=Hb*Ub和相容条件Ub+2*Sqrt(g*Hb)=UL+2*Sqrt(g*HL)联立求解，求得边界处的水深Hb和流速Ub.这是一个非线形方程，一般采用迭代法求解，或者也可变形为求解一元三次方程，直接用求根公式，避免迭代。但我的程序经常求得负根，也就是说边界处的流速和水深都是负的。迭代法和求根公式都没有问题，请问这是怎么回事?怎么解决？

cjk3d · 发表于 2005-3-22 09:57:41

你这个肯定是单向流，上游的流量和下游的水位边界都是强制给定的，那么换句话说，你的河床糙率如果给得不合适，两头的边界条件是不能封闭的。

jimzhang · 发表于 2005-3-23 11:25:11

三角形网格不好做。

cnturbulence · 发表于 2005-3-23 12:23:25

建议计算边界上作处理，考虑滩头的露出、淹没等，这样就可以避免出现负的情况。

wlh · 发表于 2005-4-12 17:40:25

检查网格的质量，我也遇到过这种情况

bdsu · 发表于 2005-6-8 09:02:48

对于网格的质量用什么来衡量、检查啊？用sms中有这功能，但其他的就不了解了，请指点下？；P

nhrisydeng · 发表于 2006-10-27 11:41:07

能给你的程序传份给我，可以吗？我的EMail:sydeng@nhri.cn

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

二维浅水方程边界条件的困惑

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。