<div dir="ltr"><div dir="ltr">On Mon, Jul 20, 2020 at 9:36 AM Jed Brown <<a href="mailto:jed@jedbrown.org">jed@jedbrown.org</a>> wrote:<br></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">How would you like to discretize the diffusive terms?  The example has a type of gradient reconstruction so you can have cellwise gradients, but there are many techniques for discretizing diffusive terms in FV.  It's simpler if you use an orthogonal grid, but I doubt that you are.<br>
<br>
As for terminology, the diffusive part is usually stiff and thus must be treated implicitly.  In TS terminology, this would be part of the IFunction, not the RHSFunction.<br></blockquote><div><br></div><div>At a high level, I would say that this is doable, but complicated. You can see me trying to do something much easier (advection + visco-elasticity) in TS ex18,</div><div>where I want to discretize the elliptic part with FEM and the advective part with FVM. I assume that is why Jed wants to know how you want to handle the</div><div>elliptic terms, since this has a large impact on how you would implement.</div><div><br></div><div>  Thanks,</div><div><br></div><div>     Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
Thibault Bridel-Bertomeu <<a href="mailto:thibault.bridelbertomeu@gmail.com" target="_blank">thibault.bridelbertomeu@gmail.com</a>> writes:<br>
<br>
> Dear all,<br>
><br>
> I have been studying ex11.c from ts/tutorials to understand how to solve an<br>
> hyperbolic system of equations using PETSCFV. I first worked on the Euler<br>
> equations for inviscid fluids and based on what ex11.c presents, I was able<br>
> to add the right PETSc instructions in an already existing in-house code<br>
> with different gas models  to solve the problems in parallel (MPI) and with<br>
> the AMR capabilities offered by P4EST.<br>
><br>
> Now my goal is to move to Navier-Stokes equations. Theoretically the system<br>
> is not completely hyperbolic and can be seen as one with an hyperbolic part<br>
> (identical to the Euler equations) and a parabolic part coming from the RHS<br>
> diffusion terms.<br>
> I have been looking into the manual and also the sources of PETSc around<br>
> the DM, DMPlex, DS and FV classes but I could not find anything that speaks<br>
> to me as "adding a RHS to an hyperbolic system of equations" or "adding a<br>
> source term to an hyperbolic system of equations". What's more, that source<br>
> term depends on the derivatives of the context variables ...<br>
><br>
> I wanted to know if anyone maybe had a suggestion regarding this issue ?<br>
><br>
> Thank you very much in advance,<br>
><br>
> Thibault Bridel-Bertomeu<br>
> —<br>
> Eng, MSc, PhD<br>
> Research Engineer<br>
> CEA/CESTA<br>
> 33114 LE BARP<br>
> Tel.: (+33)557046924<br>
> Mob.: (+33)611025322<br>
> Mail: <a href="mailto:thibault.bridelbertomeu@gmail.com" target="_blank">thibault.bridelbertomeu@gmail.com</a><br>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div>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</div><div><br></div><div><a href="http://www.cse.buffalo.edu/~knepley/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div></div></div></div>