On Fri, Feb 29, 2008 at 06:40:18PM +0100, Alessandro Frigeri wrote:
> Hello,
>
> I'd like to apply different initial conditions (temperature) to a single
> xml mesh.
>
> At the moment I load the mesh at first:
>
> Mesh mesh("mymesh.xml");
>
> then I assign the initial condition 200 K to the whole mesh:
>
> Function u0(mesh,200);
>
> The idea is to assign a different temperatures to a 'subdomain' of the
> main mesh - ideally by giving another mesh as a topological constrain.
>
> Is there any procedure you would suggest?
Take a look at the Stokes demos in src/demo/pde/stokes which
demonstrate how to set BCs based on sub domains.
But the easiest thing is probably to define a Function for the
boundary condition where you return different values depending on x in
eval().
--
Anders
_______________________________________________
DOLFIN-dev mailing list
[email protected]
http://www.fenics.org/mailman/listinfo/dolfin-dev
