butterfly mesh
 

hi... I'm a new starCD user, I'd to generate a tube with butterfly mesh (o-grid mesh)to avoid small cells near the tube center. Does anybody help me to perform that? tanks

Re: butterfly mesh
 

You can look into this thread: http://www.cfd-online.com/Forum/starcd.cgi?read=2639

Re: butterfly mesh
 

csys 1
cdel all
vdel all
ccom all
y



v 1
v 9 1
vfill 1 9
vgen 5 10 1 9 1 0 1 / 8
vset all
numb vert on
vsty 4 4
vplo
vgen 5 10 41 45 1 0 1 / 8
vset all
vplo
ctab 1 shell 2
c 1 2 12 11
cset,all
cdisp,on,vert $vstyle,4,4,1
cplot
cset all
cgen 8 1 cset
cset,all
cplot
cgen 4 10 cset
cset,all
cplot
cgen 5 10 25 28 1
cset,all
cplot
repl
zoom,off $replot
pldi,off,all $pldi,on,box $wind,0,0,13.3,10 $replot
csys 2
v 49 1 45
repl
vfill 45 49
repl
vfill 6 46 3
rp4 1 1 0
repl
v 85 1 45
vfill 45 85 3
repl
vfill 51 55
rp4 10 10
repl
vset all
vmer vset
cplot
vset all
cset all
vlis mxv
vgen 4 100 vset ,,, 0 90
cgen 4 100 cset
cset,all
cplot
zoom,off $replot
numb all off
repl
vset all
vmer vset
n
edge on
cplot
edge off
cplot
vset news edge
repl
vell 50
repl

Re: butterfly mesh
 

thanks a lot !!!!

Re: butterfly mesh
 

another script (using parameters):

!Creates a single butterfly pipe
!Set radius :
*set r1 1.5
!Set number of cells in radial direction :
!(This number should be even)
*set nr1 12
!Set length :
*set l 10
!Set no of z planes :
*set nz 25
*set nv1 1
*set nv2 nr1 + nv1
*set nr12 nr1 / 2
v nv1 0
v nv2 r1
vfil nv1 nv2
vset,,
vgen nv2 nv2 vset,,,0 r1 / nr1
ctab 2 shell 8
c nv1 nv1 + 1 nv1 + 1 + nv2 nv1 + nv2
cgen nr1 1 1
cset,,
cgen nr1 / 2 nv2 cset
*set nc1 nr1 / 2 - 1 * nr1 + 1
*set nc2 nr1 / 2 + nc1 - 1
cgen nr1 / 2 + 1 nv2 nc1 nc2
*set nvf2 nr1 / 2 * nv2 + nv2
csys 2
v nvf2 r1 45
vfil nv2 nvf2 nr1 / 2 - 1
vfil nr1 / 2 + 1 nv2 nr1 / 2 - 1,,,nr1 / 2 + 1 nv2
*set nvf3 nv2 * nr1 + nr12 + 1
vequ nvf2 nvf3
vfil nvf3 - nr12 nvf3
vfil nr12 * nv2 + 1 nvf3 - nr12 nr12 - 1,,,nr12 + 1 1
vset all
vgen 4 nv2 * nv2 vset,,,0 90
cgen 4 nv2 * nv2 all
cset all
vset news cset
vmer vset,,.001
cont
vset news edge
vsmoo 1 50
cset,,
vset news cset
vset inve
vdele vset
vcom all
y
*set dz l / nz
ctyp 1
vcex nz mxv - 1 cset,,,loca 0 0 dz
cset news flui
view 1 1 1
plty ehid
cplo