search

SPECFEM / specfem2d

SPECFEM2D simulates forward and adjoint seismic wave propagation in two-dimensional acoustic, (an)elastic, poroelastic or coupled acoustic-(an)elastic-poroelastic media, with Convolution PML absorbing conditions.
https://specfem.org
GNU General Public License v3.0
203 stars 148 forks source link

what's wrong with the mesh? #1049

Open brioglade opened 5 years ago

brioglade commented 5 years ago

Hi,Someone could help me,please. Thanks.

ispec = 10943 ngnod = 9 node 1 x,y = 93.321346000000005 75.449723000000006
node 2 x,y = 93.047447000000005 74.612702999999996
node 3 x,y = 91.870892999999995 74.638728000000000
node 4 x,y = 92.377803999999998 75.698769999999996
node 5 x,y = 93.184396000000007 75.031212999999994
node 6 x,y = 92.459170000000000 74.625715999999997
node 7 x,y = 92.124347999999998 75.168749000000005
node 8 x,y = 92.849575000000002 75.574246000000002
node 9 x,y = 92.654371999999995 75.099981000000000
negative 2D Jacobian, element saved in DX_first_element_with_negative_jacobian.dx Error detected, aborting MPI... proc 3 application called MPI_Abort(MPI_COMM_WORLD, 30) - process 3

QuentinBrissaud commented 5 years ago

Hi,

These kind of issues usually happen when you either have a poor mesh quality with stretched element or, more likely, when you have created a mesh with a wrong node numbering.

If you are using gmsh, when dealing with coupled media, make sure to define your surfaces with segments in the same order (clockwise or counter clockwise).

Best regards

Quentin Brissaud

On Sat, Sep 21, 2019, 07:00 brioglade notifications@github.com wrote:

Hi,Someone could help me,please. Thanks. (for simplicity viscoelastic or poroelastic elements, if any, are counted as elastic in the above three estimates; in reality they have more degrees of freedom)

ispec = 10943 ngnod = 9 node 1 x,y = 93.321346000000005 75.449723000000006 node 2 x,y = 93.047447000000005 74.612702999999996 node 3 x,y = 91.870892999999995 74.638728000000000 node 4 x,y = 92.377803999999998 75.698769999999996 node 5 x,y = 93.184396000000007 75.031212999999994 node 6 x,y = 92.459170000000000 74.625715999999997 node 7 x,y = 92.124347999999998 75.168749000000005 node 8 x,y = 92.849575000000002 75.574246000000002 node 9 x,y = 92.654371999999995 75.099981000000000 negative 2D Jacobian, element saved in DX_first_element_with_negative_jacobian.dx Error detected, aborting MPI... proc 3 application called MPI_Abort(MPI_COMM_WORLD, 30) - process 3

— 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/1049?email_source=notifications&email_token=ABTH7TTRVU245I22CAC4CWLQKYSGFA5CNFSM4IY6KQG2YY3PNVWWK3TUL52HS4DFUVEXG43VMWVGG33NNVSW45C7NFSM4HM2NIAA, or mute the thread https://github.com/notifications/unsubscribe-auth/ABTH7TVDKUKZSFK4HWJZXG3QKYSGFANCNFSM4IY6KQGQ .