On Wed, Oct 08, 2008 at 12:22:29PM +0100, Heitor Pina wrote: > Hi > > > > I am a Dolfin newbie trying to solve a steady heat conduction problem > > div(K grad u) + f = 0 > > For constant scalar K the variational form (I am using the Python interface) > > K*dot(grad(v),grad(u))*dx > > works OK. > > The problem is that I do not know how to instruct Dolfin to deal with the > case > where K=K(x) is a (symmetric) tensor field (anisotropic heat conduction). > > > > Any help is appreciated.
The support for tensor-valued elements is currently very limited (and in particular not handled by FFC), so you'll have to use a work-around and treat the tensor field as a vector field. Take a look at the demo TensorWeightedPoisson.form in FFC. Tensor-valued elements will be supported by UFL, the new form language that will replace the current form language in FFC. -- Anders
signature.asc
Description: Digital signature
_______________________________________________ DOLFIN-dev mailing list [email protected] http://www.fenics.org/mailman/listinfo/dolfin-dev
