On Thu, Aug 13, 2009 at 2:02 AM, Stephen Wornom
stephen.wornom at sophia.inria.fr wrote:
Matthew Knepley wrote:
On Wed, Aug 12, 2009 at 10:39 AM, Stephen Wornom
stephen.wornom at sophia.inria.fr mailto:stephen.wornom at sophia.inria.fr
wrote:
Matthew Knepley wrote:
1) You
Stephen,
There are two ways that I know of to deal with pressure checker
boarding: staggered grids or some form of Rhie-Chow interpolation.
IMO, these are simple only for uniform, Cartesian grids. For grids
that are curvilinear, unstructured, non-uniform, and/or
non-orthogonal, things get real
1) You should really handle this by creating the constant vector on the
pressure
space and using MatNullSpaceCreate()
2) You can also easily handle this by fixing the pressure at one point
Matt
On Wed, Aug 12, 2009 at 10:19 AM, William A. Perkins
william.perkins at pnl.gov wrote:
Matthew Knepley wrote:
1) You should really handle this by creating the constant vector on
the pressure
space and using MatNullSpaceCreate()
2) You can also easily handle this by fixing the pressure at one point
At what indices or location does one fix the pressure? What value is it
On Wed, Aug 12, 2009 at 10:39 AM, Stephen Wornom
stephen.wornom at sophia.inria.fr wrote:
Matthew Knepley wrote:
1) You should really handle this by creating the constant vector on the
pressure
space and using MatNullSpaceCreate()
2) You can also easily handle this by fixing the