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

MeshScale #50

Open
6 tasks
scopplestone opened this issue Nov 13, 2023 · 1 comment
Open
6 tasks

MeshScale #50

scopplestone opened this issue Nov 13, 2023 · 1 comment

Comments

@scopplestone
Copy link
Collaborator

scopplestone commented Nov 13, 2023

The small mesh dimensions require a smaller tolerance, e.g., jacobianTolerance=1e-20 (default jacobianTolerance=1e-16.).
Can this value be scaled automatically depending on the MeshScale=1e-6 or the bounding box of the complete Domain?

Also output the value of the determinant and the tolerance:

CALL abort(__STAMP__,'Element with null-negative detT found! Elem%ind: ',Elem%ind)

Todo

  • Define testcases/reggies for automatic testing the bug fix / new feature
    • 2D_Channel (see below)
    • mesh conversion from other formats?
  • Define jacobianTolerance and SpaceQuandt usage
  • Documenation
@BoltzplatzMirza
Copy link

BoltzplatzMirza commented Nov 13, 2023

ProjectName = 2D_Channel
Debugvisu = T
NVisu = 1
Mode = 1
meshScale = 1E-6

Corner =(/0.0,0.0,-2.5 ,, 1000.,0.0,-2.5 ,, 1000.,50. ,-2.5 ,, 0.0,50. ,-2.5 ,, 0.0,0.0,2.5 ,, 1000.,0.0,2.5 ,, 1000.,50. ,2.5 ,, 0.0,50. ,2.5 /)

nElems =(/400,50,1/)
elemtype =108

BCIndex =(/6 ,4 ,1 ,3 ,2 ,5/)
! =(/z-,y-,x+,y+,x-,z+/)
nZones = 1
nUserDefinedBoundaries=6
BoundaryName=BC_Xplus
BoundaryType=(/4,0,0,0/)
BoundaryName=BC_Xminus
BoundaryType=(/4,0,0,0/)
BoundaryName=BC_Yplus
BoundaryType=(/4,0,0,0/)
BoundaryName=BC_Yminus
BoundaryType=(/4,0,0,0/)
BoundaryName=BC_Zplus
BoundaryType=(/1,0,0,-1/)
BoundaryName=BC_Zminus
BoundaryType=(/1,0,0,1/)

VV=(/0. , 0. , 5E-6/) ! Displacement vector 1 (z-direction 5mü s.o.)

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