← Back to team overview

dolfin team mailing list archive

How to apply different initial conditions to a mesh domain?

 

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?

Regards,

Alessandro



-- 
Alessandro Frigeri, PhD
University of Perugia, Italy



Follow ups