[ff3d-users] How to solve when there are different regions in the domain

[Juzar Thingna]

Hello Everyone,
I'm trying to solve a set of coupled equations in a domain Omega. My domain
Omega consists of 3 regions and my functions u and v need to be continuous
across the regions. Also the constants of my problem vary across the
regions. I'm currently using the function "one" to define my constants in
different regions. But when I run my program it says that it has used 0 DOF
out of 10000 available. And thus it does not give an answer.

My domain syntax is the following

domain O =domain(S,inside(<1,0,0>) && inside(<0,0,0>) && inside (<0,1,0>));

when I remove the part && inside(<0,0,0>) && inside(<0,1,0>)and just keep
one region in my domain the program runs with some DOF(although
incorrectly).

I'm not quite sure why this is occurring? Any suggestion would be of great
help

Also I'm curious if across the regions the continuity of the function will
be maintained?

Thanks,
Juzar Thingna
Center for Computation Science and Engineering
National University of Singapore

[Stephane Del Pino]

Hello.

Are you sure that your structured mesh contains your computational domain?

Please send your command and pov files so I can help you.

Best regards,
Stéphane.

[Juzar Thingna]

Hi Stephane,
Here is my pov file

box {
<0,120,0>,<250,130,250>
pigment{
color rgb <0,0,0>}
}
box {
<3,130,3>,<123,250,123>
pigment{
color rgb <1,0,0>}
}
box {
<126,0,126>,<246,120,246>
pigment{
color rgb <0,1,0>}
}

And the computational domain is the foll

vector n=(100,100,100);
vector a=(0,0,0);
vector b=(250,250,250);
scene S=pov("graphene.pov");
mesh M=structured(n,a,b);

And the problem occurs when I specify

domain O =domain(S,inside(<1,0,0>) && inside(<0,0,0>)&& inside(<0,1,0>));

If I dont specify && inside(<0,0,0>)&& inside(<0,1,0>) it gives me some
output.
But I need to declare regions within my PDE domain. Any suggestions?

Regards,
Juzar Thingna
Center for Computation Science and Engineering
National University of Singapore

Post by Stephane Del Pino
Hello.
Are you sure that your structured mesh contains your computational domain?
Please send your command and pov files so I can help you.
Best regards,
Stéphane.
_______________________________________________
ff3d-users mailing list
http://lists.nongnu.org/mailman/listinfo/ff3d-users

[Stephane Del Pino]

Hello!

Post by Juzar Thingna
Hi Stephane,
Here is my pov file
box {
<0,120,0>,<250,130,250>
pigment{
color rgb <0,0,0>}
}
box {
<3,130,3>,<123,250,123>
pigment{
color rgb <1,0,0>}
}
box {
<126,0,126>,<246,120,246>
pigment{
color rgb <0,1,0>}
}
And the computational domain is the foll
vector n=(100,100,100);
vector a=(0,0,0);
vector b=(250,250,250);
scene S=pov("graphene.pov");
mesh M=structured(n,a,b);
And the problem occurs when I specify
domain O =domain(S,inside(<1,0,0>) && inside(<0,0,0>)&& inside(<0,1,0>));

I think you meant union of your boxes in previous mails.
So, it is not a 'and' (&&) but a 'or' (||) that you need for your domain
definition :

domain O =domain(S,inside(<1,0,0>) || inside(<0,0,0>) || inside(<0,1,0>));

Best regards,
Stéphane.