[aspect-devel] Anisotropic Viscosity?

Timo Heister timo.heister at gmail.com
Sat May 2 18:04:52 PDT 2015

This should be done in a backwards compatible way. One option would be
to allow a material model to fill in a tensor that is set to the
identity as a default in addition to the current scalar value. After
the material model is executed, the tensor will be scaled by the
scalar viscosity.

On Sat, May 2, 2015 at 9:00 PM, Jonathan Perry-Houts
<jperryh2 at uoregon.edu> wrote:
> Could that be done outside of the material models for backwards
> compatibility? Like, if anisotropy is turned off then the material is
> only responsible for returning a scalar, and otherwise then it needs to
> return a tensor? Or is that too convoluted, and better to keep it general?
> On 05/01/2015 06:41 PM, Wolfgang Bangerth wrote:
>> On 05/01/2015 03:16 PM, Thieulot, C. (Cedric) wrote:
>>> for some of the work he would like to carry out in Utrecht in a
>>> near-ish future,
>>> a simple anistropic model would be a great asset.
>>> Let’s talk at the hackathon.
>> It should actually only require two small changes in the assembly --
>> essentially two-lines. The bigger issue is that you'd have to change all
>> of the material models to multiply the viscosity they compute by the
>> identity tensor. I'd say all in all less than an hour's work.
>> Best
>>  W.
> _______________________________________________
> Aspect-devel mailing list
> Aspect-devel at geodynamics.org
> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel

More information about the Aspect-devel mailing list