On Fri, Feb 11, 2011 at 10:37 AM, Kontsantinos Kontzialis <span dir="ltr"><<a href="mailto:ckontzialis@lycos.com">ckontzialis@lycos.com</a>></span> wrote:<br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
Dear Petsc team,<br>
<br>
I'm new in Petsc and I'm trying to solve the euler equations of fluid dynamics<br>
using an implicit matrix free approach with a spatial discontinues galerkin discretization.<br>
I need some directions about how can I solve the following system:<br>
<br>
(M/dt+dR/du)*DU=R<br>
<br>
where R denotes the residual of the system and dR/du the residual jacobian. Please help.<br></blockquote><div><br></div><div>Petsc provides linear algebra and nonlinear solvers. This is fine once you have discretized.</div>
<div>It sounds like you will use DG:</div><div><br></div><div> a) on a structured or unstructured grid?</div><div><br></div><div>The PETSc DA supports structured grids in any dimension. After this, you want to use the</div>
<div>TS object to present your system. There are many examples in the TS, e.g. ex10 for</div><div>radiation-diffusion or ex14 for hydrostatic ice flow.</div><div><br></div><div>Once you have your problem producing the correct residual and Jacobian, we can talk</div>
<div>about solvers.</div><div><br></div><div> Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
Kostas<br></blockquote></div>-- <br>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener<br>