[petsc-users] CFD code with immersed boundary method

Matthew Knepley knepley at gmail.com
Wed Nov 7 08:55:02 CST 2012

On Wed, Nov 7, 2012 at 9:31 AM, Arun Kamath <arun.kamath.m at gmail.com> wrote:

> Hi!
> I am working with a CFD code  where I need to solve the Poisson Pressure
> equation. We solve a two phase flow problem which results in a jump in the
> coefficients and a non symmetric system.
> The boundary conditions are fulfilled using immersed boundary method ghost
> cells at the solid boundaries of the domain.
> I read about the DAs to deal with the boundary conditions, but they are
> applicable to only "IJK" matrices as per my understanding.
> So when I am using the "Aij" (CSR) matrix, how would you suggest I deal
> with the ghost cells of the physical domain while forming the matrix and
> the vector?

Are you saying that you have an unstructured grid with ghost cells, and you
would like PETSc to produce
a Vec of the correct size and Mat with the correct preallocation for that
grid and discretization, just as DA
does for structured grids?



> Thanks a lot.
> --
> Arun Kamath
> PhD Stipendiat
> NTNU, Trondheim

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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20121107/15336948/attachment.html>

More information about the petsc-users mailing list