二维明渠计算边界条件的给定？

liruohua

小弟在用荷兰的一个商用河道计算软件Delft-3D计算一个山区性河道，60km长，平均300m宽，坡降0.0005，地形很复杂，滩地很多，河道弯弯曲曲。用其中的一个二维模型计算。在计算时，
（一）上游给流量边界，下游给水位，算出的结果是上有没有水流，也就是水流光了，河床露出来了，没有水！
（二）上下游给水位边界，可以算，结果较合理，但流量不好控制，很难模拟真实情况。
（三）上游给水位，下游给流量，可以算，结果较合理，但这种边界条件用得很少，因为项目要专家评审，不知能不能通过。
请各位老兄帮忙解决！
SOS！！

g0water

从物理概念上说：上游给流量边界，下游给水位，这种给法是正确的
反过来概念上不合理，在数值方法和理论上可行。
可以现象一下试验吗，怎么来控制边界条件的？
可以检查以下上面给得流量方向是否正确？

tomzijizhang

一般要上游流量下游水位。但上游是水库溃坝时，就上下都是水位，但流量由特殊断面控制。
几个可能性需要注意
DELFT3D也是解明渠方程组的。二维是叫圣维南方程。三维也叫吗，不知道。
问题是，象别的同类程序一样，有数值计算上的不可理解的解出现。可能性是
（1）要让程序自动判定流动型态，急，缓，水跃的位置，过渡等等。
（2）在局部建筑物如桥梁上下游的短距离内加设断面来让程序适应建筑物。这是对程序的善意欺骗。
（3）出现坡度激变的地方要多设断面
（4）断面参数不能离水位太远，如水位只有半米，断面资料最好集中在0到1米范围。有人图省事，在地图上读不去现场测。结果地图是一米等高线分辨率的，程序插值困难和不精确。
（5）汇流断面特别是瓶颈，常常在那里不过关。建议先一条河计算通过然后加支流，一条一条的加。加一条算一条。有人一次架设整个河道网，常常不知在那里出问题了。我审查过瑞典最大的水电咨询公司的MIKE11计算结果。汇流断面不过关，出现数值振荡，就在下游设限流断面。这属于做弊。
（6）注意使用较小的时间步长。但不能太太太小，否则含时间步长的项太小，数值计算错误。
明渠计算所用的模型各家的类似，没有一家的过关。有自由水面的，流动形态又变化，不好算。要多试。祝好运。可与TOMZIJIYAHOO＠YAHOO．SE联系，可提供我的一个含有溃坝，溃堤以及洪水的算例和调查报告。也希望交流切磋。

liruohua

很高兴能在流体中文网上认识你，我在流体中文网上发的的那篇关于边界给定的问题，我试过好多次，我认为上游给流量，下游给水位计算结果不合理，原因主要是地形问题，河道太长，弯道太多，滩地极多，且有深潭，一句话，地形太复杂，太乱！我现在在试用上下游给水位边界的结果作为热启动的初始条件，然后上游给流量，下游给水位进行计算！
     你认为这样合理吗？
    我很想看看你的溃坝，溃堤以及洪水的算例和调查报告，多谢！
     你的邮箱发送不过去！

tomzijizhang

对不起，TOMZIJIYAHOO＠YAHOO．SE错了。该是TOMZIJIZHANG＠YAHOO．SE。我叫张子冀。报告有约2M大，是英文的但夹了瑞典文的证词（访问当地居民关于洪水的情况，是个实际的洪水灾害例子）。希望将来能拜读你的报告。
地形复杂只能是从上游一段一段地加。看是那个弯道，滩地出问题。如先算河道只有1/4长，看能否通过。此时下游条件只能假定--是给定水位或程序自动求的正常水深。然后加1/4。慢慢来比一次成型要好。
滩地不过关是你的断面变化太突然，程序检查出水跃出现等。该加断面--一个滩地要三个断面合适加上下游各一。

liruohua

非常感谢你的建议，我现在已经在试了，我的邮箱是：liruohua2000@163.com
我叫李若华，希望以后多交流。