komatits
commented
5 years ago Thanks! I will have a look next week. Let me cc Alexis Bottero on this.
Dimitri.
On 12/7/18 1:47 AM, JffStevens wrote:
Axisymmetric calculations are unstable in versions of SPECFEM2D released in 2018. The problem did not exist in the December 2017 version. The instability occurs on the axisymmetric axis where it intersects with the top of the absorbing boundary elements. The calculation appears to run normally for some period of time, then goes unstable. This can be observed in the EXAMPLEs: paper_axisymmetry_example and axisymmetric_case_AXISYM_option In the second case, the problem is only apparent in the last few cycles. I am attaching a figure from the last cycle of paper_axisymmetric example illustrating the problem. forward_image0012000 https://user-images.githubusercontent.com/45675229/49620592-70e03880-f976-11e8-9dd9-162b5cd10d13.jpg
— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub https://github.com/geodynamics/specfem2d/issues/1008, or mute the thread https://github.com/notifications/unsubscribe-auth/AFjDKfTgh1BNw1y5U6yKquAXN4ZhWvW2ks5u2bqEgaJpZM4ZHlcL.
-- Dimitri Komatitsch, CNRS Research Director (DR CNRS) Laboratory of Mechanics and Acoustics, Marseille, France http://komatitsch.free.fr
bottero
Axisymmetric calculations are unstable in versions of SPECFEM2D released in 2018. The problem did not exist in the December 2017 version. The instability occurs on the axisymmetric axis where it intersects with the top of the absorbing boundary elements. The calculation appears to run normally for some period of time, then goes unstable. This can be observed in the EXAMPLEs: paper_axisymmetry_example and axisymmetric_case_AXISYM_option In the second case, the problem is only apparent in the last few cycles. I am attaching a figure from the last cycle of paper_axisymmetric example illustrating the problem.