On 5 October 2016 at 18:49, Matthew Knepley <knep...@gmail.com> wrote:
> On Wed, Oct 5, 2016 at 11:19 AM, E. Tadeu <e.ta...@gmail.com> wrote: > >> Matt, >> >> Do you know if there is any example of solving Navier Stokes using a >> staggered approach by using a different DM object such as DMPlex? >> > > SNES ex62 can do P2/P1 Stokes, which is similar. Is that what you want to > see? > > For real structured grid, staggered mesh stuff like MAC, I would just do > this on a single DMDA, but think of it as being staggered, and expand my > stencil as necessary. > Following that up, for a DMDA example using a staggered grid, take a look at snes/ex30.c http://www.mcs.anl.gov/petsc/petsc-current/src/snes/examples/tutorials/ex30.c.html Thanks, Dave > > Thanks, > > Matt > > >> >> Thanks, >> Edson >> >> >> On Tue, Oct 4, 2016 at 11:12 PM, Matthew Knepley <knep...@gmail.com> >> wrote: >> >>> On Tue, Oct 4, 2016 at 9:02 PM, Somdeb Bandopadhyay <sb020...@gmail.com> >>> wrote: >>> >>>> Dear all, >>>> I want to write a solver for incompressible navier stokes >>>> using python and I want to use PETsc (particularly dmda & ksp) for this. >>>> May I know if this type of work is feasible/already done? >>>> >>> >>> How do you plan to discretize your system? DMDA supports only >>> collocation discretizations, so some sort of penalty for pressure would >>> have to be employed. >>> >>> Thanks, >>> >>> Matt >>> >>> >>>> I intend to run my solver in a cluster and so am slightly >>>> concerned about the performance if I use python with petsc. >>>> My deepest apologies if this mail of mine caused you any >>>> inconvenience. >>>> >>>> Somdeb >>>> >>> >>> >>> >>> -- >>> What most experimenters take for granted before they begin their >>> experiments is infinitely more interesting than any results to which their >>> experiments lead. >>> -- Norbert Wiener >>> >> >> > > > -- > What most experimenters take for granted before they begin their > experiments is infinitely more interesting than any results to which their > experiments lead. > -- Norbert Wiener >