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

Fixed issue with volume in triangles

Browse source 
The right value in 2D Cartesian is used for calculating the volume.
This commit is contained in:
Jorge Gonzalez 2025-10-11 14:00:16 +02:00
commit 55e062a9ef
2 changed files with 3 additions and 3 deletions
Download patch file Download diff file Expand all files Collapse all files

4
src/modules/init/moduleInput.f90
View file

@ -968,9 +968,9 @@ MODULE moduleInput
!Read the 0D mesh !Read the 0D mesh
CALL mesh%readMesh(pathMeshParticle) CALL mesh%readMesh(pathMeshParticle)
!Get the volumne !Get the volume
CALL config%get(object // '.volume', volume, found) CALL config%get(object // '.volume', volume, found)
!Rescale the volumne !Rescale the volume
IF (found) THEN IF (found) THEN
mesh%cells(1)%obj%volume = mesh%cells(1)%obj%volume*volume / Vol_ref mesh%cells(1)%obj%volume = mesh%cells(1)%obj%volume*volume / Vol_ref
mesh%nodes(1)%obj%v = mesh%cells(1)%obj%volume mesh%nodes(1)%obj%v = mesh%cells(1)%obj%volume

2
src/modules/mesh/2DCart/moduleMesh2DCart.f90
View file

@ -577,7 +577,7 @@ MODULE moduleMesh2DCart
fPsi = self%fPsi(Xi, 4) fPsi = self%fPsi(Xi, 4)
!Compute total volume of the cell !Compute total volume of the cell
self%volume = detJ*4.D0/L_ref self%volume = detJ*4.D0
!Compute volume per node !Compute volume per node
self%n1%v = self%n1%v + fPsi(1)*self%volume self%n1%v = self%n1%v + fPsi(1)*self%volume
self%n2%v = self%n2%v + fPsi(2)*self%volume self%n2%v = self%n2%v + fPsi(2)*self%volume