<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Fri, Sep 5, 2014 at 10:55 PM, Justin Chang <span dir="ltr"><<a href="mailto:jychang48@gmail.com" target="_blank">jychang48@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">Hi all,<br>
<br>
So I understand how the FEM code works in the DMPlex examples (ex12 and 62). Pardon me if this is a silly question.<br>
<br>
1) If I wanted to solve either the poisson or stokes using the discontinuous Galerkin method, is there a way to do this with the built-in DMPlex/FEM functions? Basically each cell/element has its own set of degrees of freedom, and jump/average operations would be needed to "connect" the dofs across element interfaces.<br>
<br>
2) Or how about using something like Raviart-Thomas spaces (we'll say lowest order for simplicity). Where the velocity dofs are not nodal quantities, instead they are denoted by edge fluxes (or face fluxes for tetrahedrals). Pressure would be piecewise constant.<br>
<br>
Intuitively these should be doable if I were to write my own DMPlex/PetscSection code, but I was wondering if the above two discretizations are achievable in the way ex12 and ex62 are.<br></blockquote><div><br></div><div>Lets do RT first since its easier. The primal space is</div><div><br></div><div>  P_K = Poly_{q--1}(K) + x Poly_{q-1}(K)</div><div><br></div><div>so at lowest order its just Poly_1. The dual space is moments of the normal component</div><div>of velocity on the edges. So you would write a dual space where the functionals integrated</div><div>the normal component. This is the tricky part:</div><div><br></div><div>  <a href="http://www.math.chalmers.se/~logg/pub/papers/KirbyLoggEtAl2010a.pdf">http://www.math.chalmers.se/~logg/pub/papers/KirbyLoggEtAl2010a.pdf</a><br></div><div><br></div><div>DG is just a generalization of this kind of thing where you need to a) have some geometric</div><div>quantities available to the pointwise functions (like h), and also some field quantities (like</div><div>the jump and average).</div><div><br></div><div>I understand exactly how I want to do the RT, BDM, BDMF, and NED elements, and those</div><div>will be in soon. I think DG is fairly messy and am not completely sure what I want here.</div><div><br></div><div>   Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">
Thanks,<br>
Justin<br>
</blockquote></div><br><br clear="all"><div><br></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
</div></div>