Galerkin FEM for elliptic PDEs
Nelze vybrat více než 25 témat Téma musí začínat písmenem nebo číslem, může obsahovat pomlčky („-“) a může být dlouhé až 35 znaků.

sphere2.sh 1.0KB

1234567891011121314151617181920212223242526272829303132
  1. #!/usr/bin/env bash
  2. #gmsh -3 -format vtk -o sphere.vtk sphere.geo
  3. #gmsh -2 -format stl -o sphereBox.stl sphereBox.geo
  4. #paraview --state=boundarySphere.pvsm
  5. #IMPORTANT!
  6. # Select ResampleWithDataset1
  7. # -->File -> SaveData -> MergedSphere.vtu
  8. # Set up problem
  9. # Gravity problem
  10. ../VTKEdgeGsphere MergedSphere.vtu 2.25 gravity sphereGrav.vtu
  11. ../FEM4EllipticPDE_bhmag sphereGrav.vtu sphereGrav.vtu sphereOutGrav.vtu
  12. #paraview --state=sliceSphere.pvsm
  13. # --data=sphereOut.vtu
  14. # even finer meshing
  15. # radius = 4, u = \pm .9 analytic = \pm 1.33
  16. # Finer meshing
  17. # radius = 5, u = \pm 1.65 analytic = \pm 1.67
  18. # radius = 4, u = \pm 1.09 analytic = \pm 1.33
  19. # radius = 5, u = \pm 2.5 analytic = \pm 1.67
  20. # radius = 4, u = \pm 1.66 analytic = \pm 1.33
  21. # radius = 3, u= \pm .926 -.985 analytic = \pm 1
  22. # radius = 2, u = \pm.42 analytic = \pm .667
  23. # for radius=1, u = \pm .111 analytic = \pm .333
  24. # for radius=.5, u = \pm .0477 -.0462 analytic = \pm .167