|
发表于 2011-11-11 22:07:31
|
显示全部楼层
回复 1# shirazbj 的帖子
从边定义起,使用box(示例2D channel手动定义)
boundaryCondition.setVelocityConditionOnBlockBoundaries (
lattice, Box2D(0, nx-1, 0, 0) );
boundaryCondition.setVelocityConditionOnBlockBoundaries(
lattice, Box2D(0, nx-1, ny-1, ny-1) );
boundaryCondition.setVelocityConditionOnBlockBoundaries (
lattice, Box2D(0, 0, 1, ny-2) );
boundaryCondition.setVelocityConditionOnBlockBoundaries (
lattice, Box2D(nx-1, nx-1, 1, ny-2) );
对于内部角点和外部角点如果不落在整个大的box上(对于Channel流动所有的点都在外部的box,不用重新定义边界点),需要使用(2D)
addExternalVelocityCornerNP
来定义内部点和外部点,终于一点对应的边的法向
出口的几种定义方法:
boundary:: outflow
boundary::neumann
可以直接定义在boundary上,例
boundaryCondition.setVelocityConditionOnBlockBoundaries (
lattice, Box2D(nx-1, nx-1, 1, ny-2),boundary:: outflow );
[ 本帖最后由 xhbj 于 2011-11-11 14:09 编辑 ] |
|