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

This is now an iterative process, which makes things better

Browse source
This commit is contained in:
Jorge Gonzalez 2026-03-09 18:07:06 +01:00
commit e7a47611d9
1 changed files with 6 additions and 4 deletions
Download patch file Download diff file Expand all files Collapse all files

10
src/modules/mesh/moduleMesh@boundaryParticle.f90
View file

@ -387,13 +387,15 @@ submodule(moduleMesh) boundaryParticle
end do end do
! Correction for this time step
alpha = 1.d0 - density_incident/density_rest alpha = 1.d0 - density_incident/density_rest
! Limit alpha between 0 and 1 ! Apply correction with a factor of 0.1 to avoid fast changes
alpha = min(alpha, 1.d0) self%alpha(edge%n) = self%alpha(edge%n) + 0.1d0 * alpha
alpha = max(alpha, 0.d0)
self%alpha(edge%n) = alpha ! Limit alpha between 0 and 1
self%alpha(edge%n) = min(self%alpha(edge%n), 1.d0)
self%alpha(edge%n) = max(self%alpha(edge%n), 0.d0)
deallocate(density_nodes) deallocate(density_nodes)