[petsc-users] Guidelines for solving the euler equations with an implicit matrix free approach

[Matthew Knepley]

On Fri, Feb 11, 2011 at 10:37 AM, Kontsantinos Kontzialis <
ckontzialis at lycos.com> wrote:

> Dear Petsc team,
>
>  I'm new in Petsc and I'm trying to solve the euler equations of fluid
> dynamics
> using an implicit matrix free approach with a spatial discontinues galerkin
> discretization.
> I need some directions about how can I solve the following system:
>
> (M/dt+dR/du)*DU=R
>
> where R denotes the residual of the system and dR/du the residual jacobian.
> Please help.
>

Petsc provides linear algebra and nonlinear solvers. This is fine once you
have discretized.
It sounds like you will use DG:

  a) on a structured or unstructured grid?

The PETSc DA supports structured grids in any dimension. After this, you
want to use the
TS object to present your system. There are many examples in the TS, e.g.
ex10 for
radiation-diffusion or ex14 for hydrostatic ice flow.

Once you have your problem producing the correct residual and Jacobian, we
can talk
about solvers.

   Matt


> Kostas
>
-- 
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/20110212/19a43c53/attachment.htm>