Discussion:
[ff3d-users] applying different BC's on the same domain
Juzar Thingna
2009-12-10 10:34:51 UTC
Permalink
Hello Everyone,
I'm new to FF3D and am trying to solve coupled partial differential
equations in 3-D. I'm quite unsure of how to implement different boundary
conditions on different parts of the domain. In the problem I intend to
solve the domain is a rectangular box and on one face of the box I need
Dirichlet boundary conditions but on the other faces I need Neumann boundary
conditions. Can anyone guide me as how this can be implemented using FF3D?
Also I would like to know how can we implement continuity across domains?

Eagerly awaiting a reply.

Regards,
Juzar
Research Scholar.
Center for Computational Science and Engineering.
National University of Singapore.
--
"At the end of all our searching we will arrive at the place we began and
know it for the first time."
Stephane Del Pino
2009-12-18 13:16:03 UTC
Permalink
Hello.

for instance to solve a Laplace problem in a unit square you can write

-------------------------------
vector n = (10,10,10);
vector a = (0,0,0);
vector b = (1,1,1);

mesh M = structured(n,a,b);
solve(u) in M
{
pde(u)
-div(grad(u))=1;
u = 0 on M xmin;
}
-------------------------------

Homogeneous Neumann conditions are automatically taken for the other
boundaries.

Best regards,
Stephane.
Juzar Thingna
2009-12-18 17:32:22 UTC
Permalink
Hi Stephane,
Thanks for your reply. If I have two domains and continuity of the function
is required across the domains. Does ff3d take it automatically? since the
mesh will be continuous across the domains. Please correct me if I'm wrong.

Cheers,
Juzar

On Fri, Dec 18, 2009 at 9:16 PM, Stephane Del Pino
Post by Stephane Del Pino
Hello.
for instance to solve a Laplace problem in a unit square you can write
-------------------------------
vector n = (10,10,10);
vector a = (0,0,0);
vector b = (1,1,1);
mesh M = structured(n,a,b);
solve(u) in M
{
pde(u)
-div(grad(u))=1;
u = 0 on M xmin;
}
-------------------------------
Homogeneous Neumann conditions are automatically taken for the other
boundaries.
Best regards,
Stephane.
_______________________________________________
ff3d-users mailing list
http://lists.nongnu.org/mailman/listinfo/ff3d-users
--
"At the end of all our searching we will arrive at the place we began and
know it for the first time."
Continue reading on narkive:
Loading...