On Fri, Jun 5, 2009 at 4:36 AM, Burak Atakan <[email protected]>wrote:
>
> For testing i just thought about: u=[x,y] for each cell
>
> However we fail.
>
> It is no problem if a fixed value for u is given (by the way, if a constant
> is
> used, in which direction does the flow go in a 2D system?).
>
> Also a list like [0.,-.1] is no problem.
>
> However, if I try to use the positions or the transposed positions from
> pos= mesh.getCellCenters() or
The problem is right here. Instead of using the cell centers, you should be
using the face centers. Try using something like this:
pos = mesh.getFaceCenters()
eq1 = 0 == ExponentialConvectionTerm(coeff = pos)\
+ ImplicitDiffusionTerm(coeff=a)
Mike Atkins
