-
sms
学习笔记
1
摘自
水动力数值模拟
nmap
2005-06-06 15:47
/ccb/
恩，最近放下了
mike21,
改成琢磨
sms
了。
sms
感觉上手简单一些
刚开始作了一个实例，是一个概化的河床水动力学的。
1
、制作散点地形文件。
比如我有一个河段，分成
11
个断面（包括上下游
断面）
，每个断面有
15
个点，
把这些
10*15
个点的地形按照每个断面排
列，保存为
*.xyz
文件（本例中为
）
。
2
、导入网格。
在
sms
中导入
< br>，
会出现数据导入对话框，
就是
“file import wizard
step
1
of 2
，
缺省就好了，
”
下一步
“
，
这一步需要设
置一下，
”sms data
中改为
m
esh
，
确定
“mapping op
tions
中的
这个选项选中，这是确
定你的散点
数据能够生成三角形网格，
”
完成
“
。如果这一步生成的网格不满意，就需要进行
网格的修改优化。
如果刚开始学习怕麻烦，
可以
采用规则间距的散点来生成质量
较好的网格。
3
、确定上下游边界
恩，这个步骤我也是琢磨了好久，一直没有办法翻译
Nodestring
的意思（呵呵，
太笨了）
。
单击
”create nodestring
，在边界岸边单击，然后按住
shift
键在
对岸上双击。至
于先单击哪一岸，有一条规则：面向下游从右到左（
from
right
to
left
when
looking
downstream
）
，就是面向下
游时先点右岸然后点击左岸。对另一边界做
同样操作。
至于哪一
个边界先进性也有条规则：
对于分汊河道先进行整个河段的
那一
个边界。
（
which spans across the
whole river section
）如果是没有分汊，
估计也没有必要区分那一边界先进行。
4
、赋边界条件
首先选定要设置边界条件的边界，
单击工具栏的
“sele
ct Nodestring
这时你的边
界上会出现标记（我的
sms
是一个字母
G
< br>）
，单击在
G
与边界的交界处稍
向下一
些，字母会变为
B
，表示已被选
中边界，此时边界上会在两端出现箭头，如果按
照
3
步的规则，
箭头应该指向下游，
如果不是也不要紧
，
菜单
”Nodestrings“
有
个
选项可以反转箭头方向。选择了边界线后可以设置边界条件了，在菜单
”rma2
中
“assign BC
< br>此时会变成可以使用的高亮度状态，
点击即可以设置。
另
外边界同
样操作即可。
5
、赋值初始条件
6
、模型检验
7
、运行
rma2
8
、错误检查更正
9
、后处理
sms
学习笔记
2
摘自
水动力数值模拟
nmap
2005-06-07 22:54
/ccb/
续
sms
学习笔记
1
5
、赋值初始条件
menu
--
control
里
进
行
设
置
：
主
< br>要
在
“geeral
、
、
里面设置
里
是你的工作的名称，会在结果中出现，随便设置一
个吧，无所谓的。重要的是
里的
应该是指定初始
水位高程，根据你的初始条件自己看需不需要选中（
checked
）
，
如果选中，注意该值
不能低于
节点的
地形最高高程
（？未确定）
。
中
type
是
模
拟
恒
定
还
是
非
恒
定
的
边
界
条
件
，
convergence
parameters
是水位差值的收敛参数，两个计算时断的差值小
于该
值即表示收敛就可以停止迭代。
computer
Time
及下面的相关项
为设定计算时
间。
页面中
wet/dry
testing
是处理干湿点的选项，相关详细内容见后续<
/p>
笔记，也可在
manual
中寻找，在这
里打开这个项目就好了，表示模拟工程中需
要进行干湿转化。
g
lobal roughness assignment
设置糙率，
< br>可以设置恒定糙率或
者其他方案的糙率，本例中缺省
6
、模型检验
menu
中
rma2-----
model
check...
打开模型检验，可以检查输入参
数是否合理，
一般来说可以检查很傻瓜的错误。我们缺省就好了
7
、运行
rma2
menu
中
rma2---run
rma2
这个大家都会。呵呵
8
、错误检查更正
< br>如果在
gfgen
、
rma2<
/p>
运行完了后出现
“
-->
STOP
depth
convergence
exceeds
25.0 ”
这样的提示，那么恭喜你的参数设置错误了。恩，
目前无法提供相关方面
的解决信息，见后续笔记。
如果提示
” RMA2 Finished “
那么应该计算完毕，但是
结果真实性有待验证
。
恩，
本来是想今
天和大家一起完成一个概化模型的数据处理
---
计算
-----
后处理的
全过程的，让大家首先能够
不管结果是否合理，先把
sms
运行起来有点成就感，
但是现在事务缠身，
只好以后继续了。
后续笔记
的内容目前的想法是针对具体的
参数设置进行讨论，比如网格优化、非恒定边界条件的处
理、糙率设置等等。我
也是刚上手，写笔记也是为了更好的和大家交流，水平有限，错误
惨不忍睹，希
望大家指正，大家共同进步
2
<
/p>
重新生成的动画，确实很麻烦。汗！感谢
cjk3d
、
youknow
的关注
对于利用
sms
的后处理功能的问题，据我的了解，
sms
比之
mike
、
d3d
的优势
就在于其可移植的前后处理模块，但是呢，
这种
移植时基于程序级别
的，
就是要
你用自己的计算方法来修改
rma2
的计算程序。
（这个过程很麻烦）
3
9
、后处理
sms
的后处理功能很强大，它可以创建
scale/vector
（水位流速场）
、
flow
trace
（流动示踪）和粒子示踪（
particle/drogue
trace
）三种动画。其中
scale/vecto
（水位流速场）
简单一点，后面两种稍微麻烦一点。
首先将计算后生成的网格文件
geo
和解决方案文件
sol
打开（
geo
文
件先打开，
sol
文件是基于
geo<
/p>
中的网格节点的）
，就可以开始了。
我上面的动画主要是
scale/v
ector
和
particle/drogue
trace
，
flow
trace
和
particle/drogue t
race
类似，但是简单一些，就不做介绍，详细的请查阅
ma
nual
中的
tutorials
几个问题及补充：
问题：
1.
地形散点文件是怎么做的？因为我做的时候是用日
ras
里调过
来的河道地形文
件（
.sdf
格式）<
/p>
，但由于格式不一样，需要转换成
.xyz
，这里有个转换程序，但不
是我自己的，所以不方便，自己也捉摸了半天，没搞定。<
/p>
2.
干湿判断到底干什么用啊，我觉得
也可以不选用，但这样的话会产生模型运行
不痛快，程序里提到两种干湿判断方法，可双
选也可单选，第一种判断为，当单
元其中一个节点的高程低于水深时即判断该单元为湿，
第二种
marsh
porosity
是
渐进的，即只有当单元全部节点低于计算水深时单元才变湿；
到底这干湿在专业里是个怎么解释？呵呵糊涂的说
als
里的
涡流粘性转换系数
（
Eddy
viscosity
conversion
coefficient
）问题，程序里
有建议几种取值范围，一般我们用的就是
1000~50000
，可我真不知道这个系数
和结果存在
什么关系，只知道变换这个结果影响很大，这需要学连续性方程，
但
这就是我最头疼的：
）
4.
粗率值得选择也对结果有很大影响，问题这个系数的取值国内也有相对应的取
值范围，
但这个和结果的对应关系是什么，
我想也是连
续方程里的参数，
所以。
。
。
补充：
1.
边界定义好，在付值前，在
nodestring,
< br>里选择
renumber,
这里注意多次尝试在每
个边界作一次，同时查看，
gei
info
里的
mash
info
里的几项，选择最小的那次
编号。
3.
稳态流分析，下游水位的初始定义要大于
所有节点最高高程（可以在
get
info
< br>里查看网格信息）
，但也不能太大，否则会导致程序无法运行，可以通过设置
re
vsion
（也即
.bc
文件里的
rev
card
）
选项设置一系列循环值直至真是值（一般
1
m
< br>差不多了）
。
4.
模型运行时可以用工作路径里的
.bc
文件进行控
制，这样比较方便，但需要熟
悉
rma2
程序中各
card
的设置规则。
<
/p>
5.
迭代循环，一般建议次数
3-5
次就成，
6.
收
敛系数，我们一般的研究建议
0.005-0.0001
4
-
-
-
-
-
-
-
-
-
上一篇：GMAT数学术语大全(含数学概念解释)数学小白适用!
下一篇：CGAL引用文件路径