strange loop in the forward wavefield
Opened this issue · 8 comments
When I perform forward modeling of various homogeneous elastic media, after the source is excited for a period of time, a small circle with weak energy will appear again, independent of the initial wave.
By the way, my source is the moment tensor source.
Hi,
Please remind that the color level in the output images are relative to the max/min value at each time step. So after several steps, the main wave has been absorbed by PML, then a very week energy for example a very small reflection from PML becomes visible. So what we see on the image is it.
Hi,
Please remind that the color level in the output images are relative to the max/min value at each time step. So after several steps, the main wave has been absorbed by PML, then a very week energy for example a very small reflection from PML becomes visible. So what we see on the image is it.
Sorry, I still feel confused. In this picture, we can see that the faint energy reflected by the absorption boundary is exactly the opposite direction of the source propagation, but the small circle of weak energy that I am talking about is exactly the same as the direction of the source propagation. I don't know enough about this; can you give us an example to illustrate that?
Those are the ripples. Making the mesh size or time step size half or so will reduce it I think.
Those are the ripples. Making the mesh size or time step size half or so will reduce it I think.
I tried to cut the time step in half, but I didn't get the result we wanted!It doesn't seem to be any different from before. On the other hand, it's not weak, and in my other example it's very obvious.
I recreate the same situation as below:
You can see the same ripples after the main lobe.
Here is the recorded signal at the 3rd station from the bottom
Then if I limit the y axis range, you can see the source of the ripples, which is actually very small and thus can be regarded as a numerical artifact (as I explained at first).
If you don't want to see this on the image, you can check the parameter USE_CONSTANT_MAX_AMPLITUDE in Par_file.
If we check the other artifacts (which is happening just after the main lobe) at the 3rd station from the left of this image,
then the recorded signal is
This is caused by the high frequency component of source time function at the end of the main lobe is a bit too high for the current mesh design. If we make the mesh size half (use doubled number of nx etc.) this becomes a bit better.
Thank you very much for your detailed reply, it helped me a lot.
You're welcome!