> I fixed the error in my code by following step-39. There was a mistake in > attaching the right_hand_side to MeshWorker assembler. Now the RKDG scheme > works. The solution is however oscillatory since it is discontinuous. This > is true of the original step-12 code also. Some type of limiting has to be > implemented. However limiters are very ad-hoc fixes and dont fit well in > the FEM context, and I dont know yet how to implement them inside deal.ii
But there are other approaches, such as artificial diffusion, SUPG, least squares terms, ..., that all fit well into the finite element method. The literature has ample examples for stabilization methods. W. ------------------------------------------------------------------------- Wolfgang Bangerth email: [email protected] www: http://www.math.tamu.edu/~bangerth/ _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
