dolfin team mailing list archive

[Anders Logg <logg@xxxxxxxxx>]

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

Attachment: signature.asc
Description: Digital signature