Skip to content

Add MLRK Time scheme , expand pAdaptation, and Fix LES bug in multiphase #265

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.

Sign up for GitHub

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

Merged
hmarbona merged 28 commits into from
Jul 18, 2025
Merged

Add MLRK Time scheme , expand pAdaptation, and Fix LES bug in multiphase #265

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
28 commits
Select commit Hold shift + click to select a range
5631065
Add MultiLevel RK3
Jul 16, 2025
2ac9822
Introduce pAdaptation based on Variable range on RL type refinement
Jul 16, 2025
260dab4
Fixed bug in SpatialDiscretization and LES(Multiphase)
Jul 16, 2025
20d82a6
Debug MLRK
Jul 16, 2025
f6a065a
Debug MLRK
Jul 16, 2025
9fce2d9
Debug MLRK
Jul 16, 2025
20f1bc2
Debug pAdaptationRL
Jul 16, 2025
67289fb
Debug MLRK
Jul 16, 2025
b74f560
Debug MLRK
Jul 16, 2025
41bcdf0
Debug MLRK
Jul 16, 2025
5aeb872
Debug CI Makefile Parallel Debug GFortran
Jul 16, 2025
e9fecef
Debug MLRK
Jul 16, 2025
37f45bc
Merge branch 'master' of github.com:hmarbona/horses3d
Jul 16, 2025
5bd07c8
Debug MLRK
Jul 16, 2025
231e397
Debug MLRK
Jul 16, 2025
409c62e
Debug MLRK
Jul 16, 2025
b355f78
Debug MLRK
Jul 16, 2025
30f37f3
Debug MLRK
Jul 17, 2025
5cfab82
Debug MLRK
Jul 17, 2025
2e4eff9
Debug MLRK for NS due to AnisFAS time scheme
Jul 17, 2025
9e88e49
Debug MLRK for NS due to AnisFAS time scheme
Jul 17, 2025
303ab2a
Debug MLRK for NS due to AnisFAS time scheme
Jul 17, 2025
9c8bf76
Debug MLRK for NS
Jul 17, 2025
f022d28
Update Residual Multiphysics - LES Vreman
Jul 17, 2025
2a944e3
Update Residual Multiphas Vreman LES
Jul 17, 2025
3da8709
Debug LES Multiphase
Jul 17, 2025
75cf897
Debug - Update Residual TestCase Multiphase LES
Jul 17, 2025
c9f4f65
Debug - Update Residual TestCase Multiphase LES
Jul 17, 2025
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
8 changes: 5 additions & 3 deletions Solver/src/AcousticSolver/SpatialDiscretization.f90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -111,7 +111,7 @@ end subroutine Finalize_SpaceAndTimeMethods
!
!////////////////////////////////////////////////////////////////////////
!
SUBROUTINE ComputeTimeDerivative( mesh, particles, time, mode, HO_Elements, element_mask)
SUBROUTINE ComputeTimeDerivative( mesh, particles, time, mode, HO_Elements, element_mask, Level)
IMPLICIT NONE
!
! ---------
Expand All @@ -124,6 +124,7 @@ SUBROUTINE ComputeTimeDerivative( mesh, particles, time, mode, HO_Elements, elem
integer, intent(in) :: mode
logical, intent(in), optional :: HO_Elements
logical, intent(in), optional :: element_mask(:)
integer, intent(in), optional :: Level
!
! ---------------
! Local variables
Expand Down Expand Up @@ -167,7 +168,7 @@ END SUBROUTINE ComputeTimeDerivative
! This routine computes the time derivative element by element, without considering the Riemann Solvers
! This is useful for estimating the isolated truncation error
!
SUBROUTINE ComputeTimeDerivativeIsolated( mesh, particles, time, mode, HO_Elements, element_mask)
SUBROUTINE ComputeTimeDerivativeIsolated( mesh, particles, time, mode, HO_Elements, element_mask, Level)
IMPLICIT NONE
!
! ---------
Expand All @@ -179,7 +180,8 @@ SUBROUTINE ComputeTimeDerivativeIsolated( mesh, particles, time, mode, HO_Elemen
REAL(KIND=RP) :: time
integer, intent(in) :: mode
logical, intent(in), optional :: HO_Elements
logical, intent(in), optional :: element_mask(:)
logical, intent(in), optional :: element_mask(:)
integer, intent(in), optional :: Level
!
! ---------------
! Local variables
Expand Down
8 changes: 5 additions & 3 deletions Solver/src/CahnHilliardSolver/SpatialDiscretization.f90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -118,7 +118,7 @@ end subroutine Initialize_SpaceAndTimeMethods
!
!////////////////////////////////////////////////////////////////////////
!
subroutine ComputeTimeDerivative( mesh, particles, time, mode, HO_Elements, element_mask)
subroutine ComputeTimeDerivative( mesh, particles, time, mode, HO_Elements, element_mask, Level)
IMPLICIT NONE
!
! ---------
Expand All @@ -131,6 +131,7 @@ subroutine ComputeTimeDerivative( mesh, particles, time, mode, HO_Elements, elem
integer, intent(in) :: mode
logical, intent(in), optional :: HO_Elements
logical, intent(in), optional :: element_mask(:)
integer, intent(in), optional :: Level
!
! ---------------
! Local variables
Expand Down Expand Up @@ -372,7 +373,7 @@ subroutine ComputeTimeDerivative( mesh, particles, time, mode, HO_Elements, elem

end subroutine ComputeTimeDerivative

subroutine ComputeTimeDerivativeIsolated( mesh, particles, time, mode, HO_Elements, element_mask)
subroutine ComputeTimeDerivativeIsolated( mesh, particles, time, mode, HO_Elements, element_mask, Level)
IMPLICIT NONE
!
! ---------
Expand All @@ -384,7 +385,8 @@ subroutine ComputeTimeDerivativeIsolated( mesh, particles, time, mode, HO_Elemen
REAL(KIND=RP) :: time
integer, intent(in) :: mode
logical, intent(in), optional :: HO_Elements
logical, intent(in), optional :: element_mask(:)
logical, intent(in), optional :: element_mask(:)
integer, intent(in), optional :: Level

error stop 'ComputeTimeDerivativeIsolated not implemented for Cahn-Hilliard'
end subroutine ComputeTimeDerivativeIsolated
Expand Down
6 changes: 4 additions & 2 deletions Solver/src/IncompressibleNSSolver/SpatialDiscretization.f90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -173,7 +173,7 @@ end subroutine Finalize_SpaceAndTimeMethods
!
!////////////////////////////////////////////////////////////////////////
!
SUBROUTINE ComputeTimeDerivative( mesh, particles, time, mode, HO_Elements, element_mask)
SUBROUTINE ComputeTimeDerivative( mesh, particles, time, mode, HO_Elements, element_mask, Level)
IMPLICIT NONE
!
! ---------
Expand All @@ -186,6 +186,7 @@ SUBROUTINE ComputeTimeDerivative( mesh, particles, time, mode, HO_Elements, elem
integer, intent(in) :: mode
logical, intent(in), optional :: HO_Elements
logical, intent(in), optional :: element_mask(:)
integer, intent(in), optional :: Level
!
! ---------------
! Local variables
Expand Down Expand Up @@ -263,7 +264,7 @@ END SUBROUTINE ComputeTimeDerivative
! This routine computes the time derivative element by element, without considering the Riemann Solvers
! This is useful for estimating the isolated truncation error
!
SUBROUTINE ComputeTimeDerivativeIsolated( mesh, particles, time, mode, HO_Elements, element_mask)
SUBROUTINE ComputeTimeDerivativeIsolated( mesh, particles, time, mode, HO_Elements, element_mask, Level)
use EllipticDiscretizationClass
IMPLICIT NONE
!
Expand All @@ -277,6 +278,7 @@ SUBROUTINE ComputeTimeDerivativeIsolated( mesh, particles, time, mode, HO_Elemen
integer, intent(in) :: mode
logical, intent(in), optional :: HO_Elements
logical, intent(in), optional :: element_mask(:)
integer, intent(in), optional :: Level
!
! ---------------
! Local variables
Expand Down
Loading
Loading