dolfin team mailing list archive

[Anders Logg <logg@xxxxxxxxx>]

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