Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Question about "edge noise" #6

Open
emprice opened this issue Apr 12, 2020 · 1 comment
Open

Question about "edge noise" #6

emprice opened this issue Apr 12, 2020 · 1 comment

Comments

@emprice
Copy link

emprice commented Apr 12, 2020

Hi Prof. Warburton!

I am a graduate student at Harvard University working towards my PhD in astrophysics. I have really enjoyed your nodal DG book, and I've been trying to implement the codes in it, using both the book and this repo. I've run into a bit of a snag, however, and I'm not sure what might cause this issue. I thought you might have some insight to share.

I've already checked for obvious bugs in my code, printing out matrices and such and comparing them to the output from yours (using Octave instead of MATLAB for now). So far, everything seems consistent. When I run my code on a simplified version of my target problem, I get a lot of noise at the cell edges, such as the plots I've chosen below. I'm showing both the 2-dimensional view and a slice at y = 0; the variable in both cases is the x-velocity in the Euler equations. You can see that the amplitude of the noise is quite small, but that is a bit misleading -- this particular frame is captured after an absurdly small timestep. In reality, I expect to take timesteps on the order of 5-6 orders of magnitude greater than this, and the noise only grows over time.

bokeh_plot_2
bokeh_plot

This simulation used just 2 cells, and the noise clearly follows the edge of the cell. There should be no shocks in this system, and the solution is expected to be quite smooth. I've pushed filtering to the maximum I can, and it helps a bit, but the noise always comes back through after just a couple of timesteps. I've been debugging for about a week and am really at a loss. Is this a feature inherent to DG-FEM? Is there a workaround? Or a debugging step I could try?

Thanks in advance for your time,
Ellen Price

@tcew
Copy link
Owner

tcew commented Feb 1, 2021

Hi - unfortunately I only just saw your issue post. I am missing some information:

  1. what is your PDE ?
  2. what are you boundary conditions ?
  3. what form of DG are you using ?
  4. what is your time integrator ?

TW

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants