博文

Salome 构建 Sajben Diffuser的几何型线

已有 3536 次阅读 2016-7-31 22:43 |个人分类:salome|系统分类:科研笔记

Sajben Diffuser的型线控制方程如下图所示。

取自Benjamin的硕士论文：

“Simulation and validation of compressible ﬂow in nozzle geometries and validation of OpenFOAM for this application”。

他也是引用了Bogar的文章
Bogar, T. J., Sajben, M., & Kroutil, J. C. 1983. Characteristic Frequencies of Transonic Diffuser Flow Oscillations. AIAA Journal, 21(9), 1232–1240.

下面说一下怎么构建。

用TUI命令流，为什么不用GUI呢？

因为上面有公式！而且还有一点麻烦。GUI实现起来更麻烦。索性，用TUI吧。

一些介绍参考：在Salome中使用python命令流创建curve

下面贴一下.py代码

注意：curve有三种：polyline，Bezier和Interpolation，

第1种生成多段直线，第2种是贝塞尔曲线，第3种生成一条通过个点的曲线。

# #################################################################################
# Using SALOME
# To creat the Sajben diffuser geometry, which can be found in Benjamin's master thesis,
# i.e. "Simulation and validation of compressible flow in nozzle geometries and validation of OpenFOAM for this application",
# which is located in Sec.4.2.1
# Then it will be meshed in Salome, and is send to OpenFoam for a CFD simulation
# #################################################################################

import salome
salome.salome_init()
import GEOM
from salome.geom import geomBuilder
geompy = geomBuilder.New(salome.myStudy)
gg = salome.ImportComponentGUI("GEOM")

# ###################
# set variables
# ###################

a_c= 1.4114
a_d=1.5000
l_c= -2.598
l_d= 7.216
C1_c= 0.81
C1_d= 2.25
C2_c=1.0
C2_d=0.0
C3_c= 0.5
C3_d= 0.5
C4_c= 0.6
C4_d= 0.6
pie = 3.1415926

xi_c='C1_c*(t/l_c)*pow(1.0+C2_c*(t/l_c),C3_c)/pow(1.0-(t/l_c),C4_c)'
xi_d='C1_d*(t/l_d)*pow(1.0+C2_d*(t/l_d),C3_d)/pow(1.0-(t/l_d),C4_d)'
h_c='a_c*cosh('+xi_c+')/(a_c-1.0+cosh('+xi_c+'))'
h_d='a_d*cosh('+xi_d+')/(a_d-1.0+cosh('+xi_d+'))'

x_min=4.55
x_max=9.1

# ###################
# create several polylines using parametric definition of the basic points
# ###################

# 1st: The inlet bc of Sajben diffuser. It is a line.
t_min=0.0
t_max=1.4114
t_num=1
x_constant='-4.55'
param_polyline_1 = geompy.MakeCurveParametric(x_constant, "t", "0.", t_min, t_max, t_num, GEOM.Polyline, theNewMethod=True)

# 2nd: The constant section on the top wall with a downstream covergent section. It is a line.
t_min=-4.55
t_max=-2.597
t_num=1
h_constant='1.4114'
param_polyline_2 = geompy.MakeCurveParametric("t", h_constant, "0.", t_min, t_max, t_num, GEOM.Polyline, theNewMethod=True)

# 3Rd: The covergent section. It is a curve.
t_min=-2.597
t_max=0.0
t_num=300
param_polyline_c = geompy.MakeCurveParametric("t", h_c, "0.", t_min, t_max, t_num, GEOM.Interpolation, theNewMethod=True)

# 4th: The divergent section. It is a curve.
t_min=0.0
t_max=7.215
t_num=300
param_polyline_d = geompy.MakeCurveParametric("t", h_d, "0.", t_min, t_max, t_num, GEOM.Interpolation, theNewMethod=True)

# 5th: The constant section on the top wall downstream of the divergennt section. It is a line.
t_min=7.215
t_max=9.1
t_num=1
h_constant='1.5'
param_polyline_3 = geompy.MakeCurveParametric("t", h_constant, "0.", t_min, t_max, t_num, GEOM.Polyline, theNewMethod=True)

# 6th: The outlet bc of Sajben diffuser. It is a line.
t_min=0.0
t_max=1.5
t_num=1
x_constant='9.1'
param_polyline_4 = geompy.MakeCurveParametric(x_constant, "t", "0.", t_min, t_max, t_num, GEOM.Polyline, theNewMethod=True)

# 7th: The bottom wall of Sajben diffuser. It is a line.
t_min=-4.55
t_max=9.1
t_num=1
h_constant='0.'
param_polyline_5 = geompy.MakeCurveParametric("t", h_constant, "0.", t_min, t_max, t_num, GEOM.Polyline, theNewMethod=True)

# ###################
# add objects in the study
# ###################

id_param_polyline_1 = geompy.addToStudy(param_polyline_1, "Polyline Parametric Inlet ")
id_param_polyline_2 = geompy.addToStudy(param_polyline_2, "Polyline Parametric TopConstant 1")
id_param_polyline_c = geompy.addToStudy(param_polyline_c, "Polyline Parametric Convergent")
id_param_polyline_d = geompy.addToStudy(param_polyline_d, "Polyline Parametric Divergent")
id_param_polyline_3 = geompy.addToStudy(param_polyline_3, "Polyline Parametric TopConstant 2")
id_param_polyline_4 = geompy.addToStudy(param_polyline_4, "Polyline Parametric Outlet")
id_param_polyline_5 = geompy.addToStudy(param_polyline_5, "Polyline Parametric Bottom")

# ###################
# display the objects
# ###################
gg.createAndDisplayGO(id_param_polyline_1)
gg.createAndDisplayGO(id_param_polyline_2)
gg.createAndDisplayGO(id_param_polyline_c)
gg.createAndDisplayGO(id_param_polyline_d)
gg.createAndDisplayGO(id_param_polyline_3)
gg.createAndDisplayGO(id_param_polyline_4)
gg.createAndDisplayGO(id_param_polyline_5)

具体怎么用呢？

打开Salome，

界面是这样的：