JGonzalez/fpakc
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

Final commit for volume nodes

Browse source 
1D Rad is fixed.

Still there might be an issue with the triangles in 2D Cyl.
This commit is contained in:
Jorge Gonzalez 2022-12-24 11:41:12 +01:00
commit a8bea6102f
2 changed files with 14 additions and 3 deletions
Download patch file Download diff file Expand all files Collapse all files

9
src/modules/mesh/1DRad/moduleMesh1DRad.f90
View file

@ -260,10 +260,17 @@ MODULE moduleMesh1DRad
Xii = 0.D0 Xii = 0.D0
fPsi = self%fPsi(Xii) fPsi = self%fPsi(Xii)
detJ = self%detJac(Xii) detJ = self%detJac(Xii)
!Computes total volume of the cell
r = DOT_PRODUCT(fPsi, self%r) r = DOT_PRODUCT(fPsi, self%r)
l = 2.D0*detJ l = 2.D0*detJ
self%volume = r*l self%volume = r*l
self%arNodes = fPsi*r*l !Computes volume per node
xi = (/-5.D-1, 0.D0, 0.D0/)
r = DOT_PRODUCT(self%fPsi(xi),self%r)
self%arNodes(1) = fPsi(1)*r*l
xi = (/ 5.D-1, 0.D0, 0.D0/)
r = DOT_PRODUCT(self%fPsi(xi),self%r)
self%arNodes(2) = fPsi(2)*r*l
END SUBROUTINE areaRad END SUBROUTINE areaRad

4
src/modules/mesh/2DCyl/moduleMesh2DCyl.f90
View file

@ -318,8 +318,10 @@ MODULE moduleMesh2DCyl
xi = 0.D0 xi = 0.D0
detJ = self%detJac(xi)*PI8 !4*2*pi detJ = self%detJac(xi)*PI8 !4*2*pi
fPsi = self%fPsi(xi) fPsi = self%fPsi(xi)
!Computes total volume of the cell
r = DOT_PRODUCT(fPsi,self%r) r = DOT_PRODUCT(fPsi,self%r)
self%volume = r*detJ self%volume = r*detJ
!Computes volume per node
xi = (/-5.D-1, -5.D-1, 0.D0/) xi = (/-5.D-1, -5.D-1, 0.D0/)
r = DOT_PRODUCT(self%fPsi(xi),self%r) r = DOT_PRODUCT(self%fPsi(xi),self%r)
self%arNodes(1) = fPsi(1)*r*detJ self%arNodes(1) = fPsi(1)*r*detJ
@ -712,8 +714,10 @@ MODULE moduleMesh2DCyl
xi = (/1.D0/3.D0, 1.D0/3.D0, 0.D0 /) xi = (/1.D0/3.D0, 1.D0/3.D0, 0.D0 /)
detJ = self%detJac(xi)*PI !2PI*1/2 detJ = self%detJac(xi)*PI !2PI*1/2
fPsi = self%fPsi(xi) fPsi = self%fPsi(xi)
!Computes total volume of the cell
r = DOT_PRODUCT(fPsi,self%r) r = DOT_PRODUCT(fPsi,self%r)
self%volume = r*detJ self%volume = r*detJ
!Computes volume per node
self%arNodes = fPsi*r*detJ self%arNodes = fPsi*r*detJ
END SUBROUTINE areaTria END SUBROUTINE areaTria