Sphere in Cube
 

I am creating three dimensional cube with reactangular coordiate system and exactly at the center I want to creat sphere with spherical coordiate system. I can create them as and where I want but finally it given connectivity error, it shows some of the unconnected regions. Does anybody know about it. Thank you in advance.

Sameer

Re: Sphere in Cube
 

cdel all
vdel all
ctdel 4 20
ctab 2 flui 4
loca 10 sphere

! number of smoothing loops
*set smoo 5

csys 1
ctype 1
vc3d -1 1 20 -1 1 20 -1 1 20
cset all
cset subs gran -0.75 0.75 -0.75 0.75 -0.75 0.75
ctype 2
cmod cset
live surf create
cset news type 2
cdel cset

cset news type 4
vcex 5 ,, cset ,,, norm -0.05

cset all
live surf create
cset news type 1
vset news cset
cset add vset all
cset subs type 5
cdel cset

cset news type 5
vset news cset
vmod vset 0.5

*defi noexec
vset none
vell 20
vset news cset
vmod vset 0.5
*end
*loop 0 smoo

cset news flui
vset news surf
vell 200

cset subs gran -2 0.05
axis y
view 1 0.5 1
cplo

Re: Sphere in Cube
 

Dear Joern. Thanks a lot for your help. I do appreciate the time and consideration.

Re: Sphere in Cube
 

Dear Joern,

I experience little problem running this script. Can you please check the script again. I am very thankful for your kind help. I udestood most of the script. And finally it gives 'Connectivity' error. Thank you for your time and consideration. Regards,

Sameer

Re: Sphere in Cube
 

I moved the line "loca 10 sphere" to the top of the code. So coordinate system 10 was not active when it should be for the "vmod" command.

Just add a "csys 10" before the first "vmod vset 0.5"

Re: Sphere in Cube
 

Dear Joern,

Thanks a lot for your help. Can you please check the following command lines. I am creating solid sphere exactly at the center of the cube which is filled of the water. I have all the commands. i can successfully create the cube and sphere inside it. But it gives connectivity erroe and shows two groups of cells. Please run the following script

ctab 1 fluid 2

csys 1

ctype 1

vc3d -10 10 40 -10 10 40 -10 10 40

cset all

cplot

cset news gran -0.50 0.50 -0.50 0.50 -0.50 0.50

cplot

ctab,2,solid,3,,2

ctname,2,solid

cmod cset

cplot

crefine 5 5 5 cset 0 0 merg couple 1

cmod cset

cplot

cset news type 2

cplot

live surf create

cset news type 3

vcex 10,,cset,,,norm 0.05

cset news type 2

cplot

cset news type 3

cdel cset

cset news type 2

cplot

live surf create

cset news type 4

vset none

vell 15

loca 10 sphere

vset news cset

vmod vset 1

cset news type 2

vell 200

cplot

Can you please let me know how to remove this error of connectivity. And that will be all for my problem.

Thanks

Sameer.

Re: Sphere in Cube
 

Coupling between fluid and solid cells is not allowed.

Re: Sphere in Cube
 

Then what is the best idea to proceed further. How do I remove that connectivity error. I went through manual and applied all the possible commands like Vmerge and all but still it shows error in that. I am clueless at this point of time. Thanks

Sameer

Re: Sphere in Cube
 

Couples are not allowed between fluid and solid regions. So don't do the "crefine" for only the solid cells.

Run the macro without the crefine at all to see what you get and use some more cells for the initial block.

Re: Sphere in Cube
 

Cell are sometimes left-handed after a vcextrude command. Just open the check-tool and use the option "right-handed". Then click on "fix left-handed".

Here is the complete macro:


cdel all
vdel all
ctdel 4 20
ctab 2 flui 4

! number of smoothing loops
*set smoo 5

csys 1
ctype 1
vc3d -1 1 20 -1 1 20 -1 1 20
cset all
cset subs gran -0.75 0.75 -0.75 0.75 -0.75 0.75
ctype 2
cmod cset
live surf create
cset news type 2
cdel cset

cset news type 4
vcex 5 ,, cset ,,, norm -0.05

cset all
live surf create
cset news type 1
vset news cset
cset add vset all
cset subs type 5
cdel cset

cset news type 5
vset news cset
loca 10 sphere
csys 10
vmod vset 0.5

*defi noexec
vset none
vell 20
vset news cset
vmod vset 0.5
*end
*loop 0 smoo

cset news flui
vset news surf
vell 200

chec,all,0,righ,news,noli
cflip,cset
check,all

cset all
cset subs gran -2 0.05
axis y
view 1 0.5 1
cplo