Hi all, This might not be relevant, but I am working on boundary conditions that not only depend on spatial positions but also a couple of variables. For example,
mesh=... u=CellVariable(mesh=mesh,...) v=CellVariable(mesh=mesh,...) bcs=FixedValue(faces=mesh.getFacesLeft(),value=u.getFaceGrad()*v.getFaceValue()) Do you think this will cause convergence difficulties? Thanks, Zhiwen On Fri, Nov 7, 2008 at 3:03 PM, Jonathan Guyer <[EMAIL PROTECTED]> wrote: > > > On Nov 7, 2008, at 1:21 PM, Daniel Lewis wrote: > > This information would be nice for developing spatially varying boundary >> conditions... Just a thought. >> > > I think you can already do this. > > >>> X, Y = mesh.getFaceCenters() > >>> bcs = FixedValue(faces=mesh.getFacesLeft(), value=Y) > > Or did you have something else in mind? > >
