+++ /dev/null
-# Viscous layers construction
-
-import salome
-salome.salome_init_without_session()
-import GEOM
-from salome.geom import geomBuilder
-geompy = geomBuilder.New()
-
-import SMESH, SALOMEDS
-from salome.smesh import smeshBuilder
-smesh = smeshBuilder.New()
-
-X = geompy.MakeVectorDXDYDZ( 1,0,0 )
-O = geompy.MakeVertex( 100,50,50 )
-plane = geompy.MakePlane( O, X, 200 ) # plane YZ
-
-box = geompy.MakeBoxDXDYDZ(200,100,100)
-
-shape = geompy.MakeHalfPartition( box, plane )
-
-faces = geompy.SubShapeAllSorted(shape, geompy.ShapeType["FACE"])
-face1 = faces[1]
-ignoreFaces = [ faces[0], faces[-1]]
-
-geompy.addToStudy( shape, "shape" )
-geompy.addToStudyInFather( shape, face1, "face1")
-
-# 3D Viscous layers
-
-mesh = smesh.Mesh(shape, "CFD")
-
-mesh.Segment().NumberOfSegments( 4 )
-
-mesh.Triangle()
-mesh.Quadrangle(face1)
-algo3D = mesh.Tetrahedron()
-
-thickness = 20
-numberOfLayers = 10
-stretchFactor = 1.5
-groupName = "Boundary layers"
-layersHyp = algo3D.ViscousLayers(thickness,numberOfLayers,stretchFactor,
- ignoreFaces, # optional
- groupName = groupName) # optional
-
-mesh.Compute()
-
-# retrieve boundary prisms created by mesh.Compute()
-boundaryGroup = mesh.GetGroupByName( layersHyp.GetGroupName() )[0]
-print( "Nb boundary prisms", boundaryGroup.Size() )
-
-mesh.MakeGroup("Tetras",SMESH.VOLUME,SMESH.FT_ElemGeomType,"=",SMESH.Geom_TETRA)
-mesh.MakeGroup("Pyras",SMESH.VOLUME,SMESH.FT_ElemGeomType,"=",SMESH.Geom_PYRAMID)
-mesh.MakeGroup("Prims",SMESH.VOLUME,SMESH.FT_ElemGeomType,"=",SMESH.Geom_PENTA)
-
-# 2D Viscous layers
-
-# 3 edges of the 4 edges of face1
-edgeIds = geompy.SubShapeAllIDs( face1, geompy.ShapeType["EDGE"])[:-1]
-
-mesh = smesh.Mesh(face1,"VicsousLayers2D")
-mesh.Segment().NumberOfSegments( 5 )
-
-# viscous layers will be created on 1 edge, as we set 3 edges to ignore
-vlHyp = mesh.Triangle().ViscousLayers2D( 2, 3, 1.5,
- edgeIds, isEdgesToIgnore=True, # optional
- groupName=groupName) # optional
-mesh.Compute()
-
-# retrieve boundary elements created by mesh.Compute()
-quadrangles = mesh.GetGroupByName( vlHyp.GetGroupName() )[0]
-print( "Nb boundary quadrangles", quadrangles.Size() )
-
-# viscous layers will be created on 3 edges, as we pass isEdgesToIgnore=False
-vlHyp.SetEdges( edgeIds, False )
-
-mesh.Compute()
