On Thu, Apr 18, 2013 at 11:06 AM, Kirk, Benjamin (JSC-EG311) < [email protected]> wrote:
> > Well you did say inviscid, correct? Then I would think the > discretizations would be equivalent - modulo your choice of tau. Which is > what exactly? > > I am using the same tau definition as in your dissertation. tau_mat = diag(tau_c, tau_m, tau_m, tau_m, tau_e) tau_c = tau_m = tau_e = tau tau = ( (2/dt)^2 + ( (2\|u\| + a) / h_u )^2 )^-0.5 The discontinuity capturing operator, delta, is also the same as in your work. That brings me to another question. I am looking at using higher order elements and using this tau and delta leads to spurious oscillations in problems with shocks. Do you have a recommendation on a better tau and delta definition for higher order elements? I was considering deriving these matrices based on the residual-free bubbles approach, but could also experiment with definitions. Manav ------------------------------------------------------------------------------ Precog is a next-generation analytics platform capable of advanced analytics on semi-structured data. The platform includes APIs for building apps and a phenomenal toolset for data science. Developers can use our toolset for easy data analysis & visualization. Get a free account! http://www2.precog.com/precogplatform/slashdotnewsletter _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
