Feature/111 inflow outflow bcs

src/SELF_DGModel2D_t.f90

@@ -479,6 +479,33 @@
               this%hbc2d_NoNormalFlow(this%solution%boundary(i,j,iEl,1:this%nvar),nhat)
           enddo
 
+        elseif(bcid == SELF_BC_INFLOW) then
+
+          do i = 1,this%solution%interp%N+1 ! Loop over quadrature points
+            nhat = this%geometry%nhat%boundary(i,j,iEl,1,1:2)
+
+            this%solution%extBoundary(i,j,iEl,1:this%nvar) = &
+              this%hbc2d_Inflow(this%solution%boundary(i,j,iEl,1:this%nvar),nhat)
+          enddo
+
+        elseif(bcid == SELF_BC_OUTFLOW) then
+
+          do i = 1,this%solution%interp%N+1 ! Loop over quadrature points
+            nhat = this%geometry%nhat%boundary(i,j,iEl,1,1:2)
+
+            this%solution%extBoundary(i,j,iEl,1:this%nvar) = &
+              this%hbc2d_Outflow(this%solution%boundary(i,j,iEl,1:this%nvar),nhat)
+          enddo
+
+        elseif(bcid == SELF_BC_NOSLIP) then
+
+          do i = 1,this%solution%interp%N+1 ! Loop over quadrature points
+            nhat = this%geometry%nhat%boundary(i,j,iEl,1,1:2)
+
+            this%solution%extBoundary(i,j,iEl,1:this%nvar) = &
+              this%hbc2d_noslip(this%solution%boundary(i,j,iEl,1:this%nvar),nhat)
+          enddo
+
         endif
       endif
 
@@ -534,6 +561,33 @@
               this%pbc2d_NoNormalFlow(dsdx,nhat)
           enddo
 
+        elseif(bcid == SELF_BC_INFLOW) then
+
+          do i = 1,this%solution%interp%N+1 ! Loop over quadrature points
+            nhat = this%geometry%nhat%boundary(i,j,iEl,1,1:2)
+
+            this%solutiongradient%extBoundary(i,j,iEl,1:this%nvar,1:2) = &
+              this%pbc2d_Inflow(this%solution%boundary(i,j,iEl,1:this%nvar),nhat)
+          enddo
+
+        elseif(bcid == SELF_BC_OUTFLOW) then
+
+          do i = 1,this%solution%interp%N+1 ! Loop over quadrature points
+            nhat = this%geometry%nhat%boundary(i,j,iEl,1,1:2)
+
+            this%solutiongradient%extBoundary(i,j,iEl,1:this%nvar,1:2) = &
+              this%pbc2d_Outflow(this%solution%boundary(i,j,iEl,1:this%nvar),nhat)
+          enddo
+
+        elseif(bcid == SELF_BC_NOSLIP) then
+
+          do i = 1,this%solution%interp%N+1 ! Loop over quadrature points
+            nhat = this%geometry%nhat%boundary(i,j,iEl,1,1:2)
+
+            this%solutiongradient%extBoundary(i,j,iEl,1:this%nvar,1:2) = &
+              this%pbc2d_noslip(this%solution%boundary(i,j,iEl,1:this%nvar),nhat)
+          enddo
+
         endif
       endif
 

src/SELF_DGModel3D_t.f90

@@ -481,6 +481,39 @@
             enddo
           enddo
 
+        elseif(bcid == SELF_BC_INFLOW) then
+
+          do j = 1,this%solution%interp%N+1 ! Loop over quadrature points
+            do i = 1,this%solution%interp%N+1 ! Loop over quadrature points
+              nhat = this%geometry%nhat%boundary(i,j,k,iEl,1,1:3)
+
+              this%solution%extBoundary(i,j,k,iEl,1:this%nvar) = &
+                this%hbc3d_Inflow(this%solution%boundary(i,j,k,iEl,1:this%nvar),nhat)
+            enddo
+          enddo
+
+        elseif(bcid == SELF_BC_OUTFLOW) then
+
+          do j = 1,this%solution%interp%N+1 ! Loop over quadrature points
+            do i = 1,this%solution%interp%N+1 ! Loop over quadrature points
+              nhat = this%geometry%nhat%boundary(i,j,k,iEl,1,1:3)
+
+              this%solution%extBoundary(i,j,k,iEl,1:this%nvar) = &
+                this%hbc3d_Outflow(this%solution%boundary(i,j,k,iEl,1:this%nvar),nhat)
+            enddo
+          enddo
+
+        elseif(bcid == SELF_BC_NOSLIP) then
+
+          do j = 1,this%solution%interp%N+1 ! Loop over quadrature points
+            do i = 1,this%solution%interp%N+1 ! Loop over quadrature points
+              nhat = this%geometry%nhat%boundary(i,j,k,iEl,1,1:3)
+
+              this%solution%extBoundary(i,j,k,iEl,1:this%nvar) = &
+                this%hbc3d_NoSlip(this%solution%boundary(i,j,k,iEl,1:this%nvar),nhat)
+            enddo
+          enddo
+
         endif
       endif
 
@@ -542,6 +575,45 @@
             enddo
           enddo
 
+        elseif(bcid == SELF_BC_INFLOW) then
+
+          do j = 1,this%solutiongradient%interp%N+1 ! Loop over quadrature points
+            do i = 1,this%solutiongradient%interp%N+1 ! Loop over quadrature points
+              nhat = this%geometry%nhat%boundary(i,j,k,iEl,1,1:3)
+
+              dsdx = this%solutiongradient%boundary(i,j,k,iEl,1:this%nvar,1:3)
+
+              this%solutiongradient%extBoundary(i,j,k,iEl,1:this%nvar,1:3) = &
+                this%pbc3d_Inflow(dsdx,nhat)
+            enddo
+          enddo
+
+        elseif(bcid == SELF_BC_OUTFLOW) then
+
+          do j = 1,this%solutiongradient%interp%N+1 ! Loop over quadrature points
+            do i = 1,this%solutiongradient%interp%N+1 ! Loop over quadrature points
+              nhat = this%geometry%nhat%boundary(i,j,k,iEl,1,1:3)
+
+              dsdx = this%solutiongradient%boundary(i,j,k,iEl,1:this%nvar,1:3)
+
+              this%solutiongradient%extBoundary(i,j,k,iEl,1:this%nvar,1:3) = &
+                this%pbc3d_Outflow(dsdx,nhat)
+            enddo
+          enddo
+
+        elseif(bcid == SELF_BC_NOSLIP) then
+
+          do j = 1,this%solutiongradient%interp%N+1 ! Loop over quadrature points
+            do i = 1,this%solutiongradient%interp%N+1 ! Loop over quadrature points
+              nhat = this%geometry%nhat%boundary(i,j,k,iEl,1,1:3)
+
+              dsdx = this%solutiongradient%boundary(i,j,k,iEl,1:this%nvar,1:3)
+
+              this%solutiongradient%extBoundary(i,j,k,iEl,1:this%nvar,1:3) = &
+                this%pbc3d_NoSlip(dsdx,nhat)
+            enddo
+          enddo
+
         endif
       endif
 

src/SELF_Mesh.f90

@@ -96,9 +96,13 @@ subroutine SELF_FreeMesh(this)
   integer,parameter :: SELF_BC_PRESCRIBED = 100
   integer,parameter :: SELF_BC_RADIATION = 101
   integer,parameter :: SELF_BC_NONORMALFLOW = 102
+  integer,parameter :: SELF_BC_INFLOW = 103
+  integer,parameter :: SELF_BC_OUTFLOW = 104
+  integer,parameter :: SELF_BC_NOSLIP = 105
 
   ! Conditions on the solution gradients
   integer,parameter :: SELF_BC_PRESCRIBED_STRESS = 200
-  integer,parameter :: SELF_BC_NOSTRESS = 201
+  integer,parameter :: SELF_BC_NOSTRESS = 201 ! Defaults to setting the external gradient field as the negative of the internal gradient field so the average is zero.
+  integer,parameter :: SELF_BC_STRESS = 202 ! Defaults to prolonging the internal gradient field to the external gradient field
 
 endmodule SELF_Mesh

src/SELF_Mesh_3D_t.f90

@@ -90,7 +90,8 @@ module SELF_Mesh_3D_t
 !    2 - Global Side ID
 !    3 - Neighbor Element ID
 !    4 - 10*( neighbor local side )  + flip
-!    5 - Boundary Condition ID
+!    5 - Boundary Condition ID (hyperbolic part)
+!    6 - Boundary Condition ID (parabolic part) ! TODO
 !
 !
 ! ========================================================================= !

src/SELF_Model.f90

@@ -160,9 +160,15 @@ module SELF_Model
     procedure :: hbc2d_Prescribed => hbc2d_Prescribed_Model
     procedure :: hbc2d_Radiation => hbc2d_Generic_Model
     procedure :: hbc2d_NoNormalFlow => hbc2d_Generic_Model
+    procedure :: hbc2d_Inflow => hbc2d_Generic_Model
+    procedure :: hbc2d_Outflow => hbc2d_Generic_Model
+    procedure :: hbc2d_NoSlip => hbc2d_Generic_Model
     procedure :: hbc3d_Prescribed => hbc3d_Prescribed_Model
     procedure :: hbc3d_Radiation => hbc3d_Generic_Model
     procedure :: hbc3d_NoNormalFlow => hbc3d_Generic_Model
+    procedure :: hbc3d_Inflow => hbc3d_Generic_Model
+    procedure :: hbc3d_Outflow => hbc3d_Generic_Model
+    procedure :: hbc3d_NoSlip => hbc3d_Generic_Model
 
     ! Boundary condition functions (parabolic)
     procedure :: pbc1d_Prescribed => pbc1d_Prescribed_Model
@@ -171,9 +177,15 @@ module SELF_Model
     procedure :: pbc2d_Prescribed => pbc2d_Prescribed_Model
     procedure :: pbc2d_Radiation => pbc2d_Generic_Model
     procedure :: pbc2d_NoNormalFlow => pbc2d_Generic_Model
+    procedure :: pbc2d_Inflow => pbc2d_Generic_Model
+    procedure :: pbc2d_Outflow => pbc2d_Generic_Model
+    procedure :: pbc2d_NoSlip => pbc2d_Generic_Model
     procedure :: pbc3d_Prescribed => pbc3d_Prescribed_Model
     procedure :: pbc3d_Radiation => pbc3d_Generic_Model
     procedure :: pbc3d_NoNormalFlow => pbc3d_Generic_Model
+    procedure :: pbc3d_Inflow => pbc3d_Generic_Model
+    procedure :: pbc3d_Outflow => pbc3d_Generic_Model
+    procedure :: pbc3d_NoSlip => pbc3d_Generic_Model
 
     procedure :: ReportEntropy => ReportEntropy_Model
     procedure :: ReportMetrics => ReportMetrics_Model

src/cpu/SELF_CIGLES3D.f90

@@ -0,0 +1,36 @@
+! //////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////// !
+!
+! Maintainers : support@fluidnumerics.com
+! Official Repository : https://github.com/FluidNumerics/self/
+!
+! Copyright © 2024 Fluid Numerics LLC
+!
+! Redistribution and use in source and binary forms, with or without modification, are permitted provided that the following conditions are met:
+!
+! 1. Redistributions of source code must retain the above copyright notice, this list of conditions and the following disclaimer.
+!
+! 2. Redistributions in binary form must reproduce the above copyright notice, this list of conditions and the following disclaimer in
+!    the documentation and/or other materials provided with the distribution.
+!
+! 3. Neither the name of the copyright holder nor the names of its contributors may be used to endorse or promote products derived from
+!    this software without specific prior written permission.
+!
+! THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS “AS IS” AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
+! LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT
+! HOLDER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT
+! LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY
+! THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF
+! THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
+!
+! //////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////// !
+
+module self_CIGLES3D
+
+  use self_CIGLES3D_t
+
+  implicit none
+
+  type,extends(CIGLES3D_t) :: CIGLES3D
+  endtype CIGLES3D
+
+endmodule self_CIGLES3D

src/gpu/SELF_CIGLES3D.f90

@@ -0,0 +1,151 @@
+! //////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////// !
+!
+! Maintainers : support@fluidnumerics.com
+! Official Repository : https://github.com/FluidNumerics/self/
+!
+! Copyright © 2024 Fluid Numerics LLC
+!
+! Redistribution and use in source and binary forms, with or without modification, are permitted provided that the following conditions are met:
+!
+! 1. Redistributions of source code must retain the above copyright notice, this list of conditions and the following disclaimer.
+!
+! 2. Redistributions in binary form must reproduce the above copyright notice, this list of conditions and the following disclaimer in
+!    the documentation and/or other materials provided with the distribution.
+!
+! 3. Neither the name of the copyright holder nor the names of its contributors may be used to endorse or promote products derived from
+!    this software without specific prior written permission.
+!
+! THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS “AS IS” AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
+! LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT
+! HOLDER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT
+! LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY
+! THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF
+! THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
+!
+! //////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////// !
+
+module self_CIGLES3D
+
+  use self_CIGLES3D_t
+
+  implicit none
+
+  type,extends(CIGLES3D_t) :: CIGLES3D
+  contains
+    ! procedure :: setboundarycondition => setboundarycondition_CIGLES3D
+    ! procedure :: boundaryflux => boundaryflux_CIGLES3D
+    ! procedure :: fluxmethod => fluxmethod_CIGLES3D
+
+  endtype CIGLES3D
+
+  ! interface
+  !   subroutine setboundarycondition_CIGLES3D_gpu(extboundary,boundary,sideinfo,nhat,N,nel) &
+  !     bind(c,name="setboundarycondition_CIGLES3D_gpu")
+  !     use iso_c_binding
+  !     type(c_ptr),value :: extboundary,boundary,sideinfo,nhat
+  !     integer(c_int),value :: N,nel
+  !   endsubroutine setboundarycondition_CIGLES3D_gpu
+  ! endinterface
+
+  ! interface
+  !   subroutine fluxmethod_CIGLES3D_gpu(solution,flux,rho0,c,N,nel,nvar) &
+  !     bind(c,name="fluxmethod_CIGLES3D_gpu")
+  !     use iso_c_binding
+  !     use SELF_Constants
+  !     type(c_ptr),value :: solution,flux
+  !     real(c_prec),value :: rho0,c
+  !     integer(c_int),value :: N,nel,nvar
+  !   endsubroutine fluxmethod_CIGLES3D_gpu
+  ! endinterface
+
+  ! interface
+  !   subroutine boundaryflux_CIGLES3D_gpu(fb,fextb,nhat,nscale,flux,rho0,c,N,nel) &
+  !     bind(c,name="boundaryflux_CIGLES3D_gpu")
+  !     use iso_c_binding
+  !     use SELF_Constants
+  !     type(c_ptr),value :: fb,fextb,flux,nhat,nscale
+  !     real(c_prec),value :: rho0,c
+  !     integer(c_int),value :: N,nel
+  !   endsubroutine boundaryflux_CIGLES3D_gpu
+  ! endinterface
+
+contains
+
+  ! subroutine boundaryflux_CIGLES3D(this)
+  !   implicit none
+  !   class(CIGLES3D),intent(inout) :: this
+
+  !   call boundaryflux_CIGLES3D_gpu(this%solution%boundary_gpu, &
+  !                                       this%solution%extBoundary_gpu, &
+  !                                       this%geometry%nhat%boundary_gpu, &
+  !                                       this%geometry%nscale%boundary_gpu, &
+  !                                       this%flux%boundarynormal_gpu, &
+  !                                       this%rho0,this%c,this%solution%interp%N, &
+  !                                       this%solution%nelem)
+
+  ! endsubroutine boundaryflux_CIGLES3D
+
+  ! subroutine fluxmethod_CIGLES3D(this)
+  !   implicit none
+  !   class(CIGLES3D),intent(inout) :: this
+
+  !   call fluxmethod_CIGLES3D_gpu(this%solution%interior_gpu, &
+  !                                     this%flux%interior_gpu, &
+  !                                     this%rho0,this%c,this%solution%interp%N,this%solution%nelem, &
+  !                                     this%solution%nvar)
+
+  ! endsubroutine fluxmethod_CIGLES3D
+
+  ! subroutine setboundarycondition_CIGLES3D(this)
+  !   !! Boundary conditions are set to periodic boundary conditions
+  !   implicit none
+  !   class(CIGLES3D),intent(inout) :: this
+  !   ! local
+  !   integer :: i,iEl,j,k,e2,bcid
+  !   real(prec) :: x(1:3)
+
+  !   if(this%prescribed_bcs_enabled) then
+  !     call gpuCheck(hipMemcpy(c_loc(this%solution%extboundary), &
+  !                             this%solution%extboundary_gpu,sizeof(this%solution%extboundary), &
+  !                             hipMemcpyDeviceToHost))
+
+  !     ! Prescribed boundaries are still done on the CPU
+  !     do iEl = 1,this%solution%nElem ! Loop over all elements
+  !       do k = 1,6 ! Loop over all sides
+
+  !         bcid = this%mesh%sideInfo(5,j,iEl) ! Boundary Condition ID
+  !         e2 = this%mesh%sideInfo(3,j,iEl) ! Neighboring Element ID
+
+  !         if(e2 == 0) then
+  !           if(bcid == SELF_BC_PRESCRIBED) then
+
+  !             do j = 1,this%solution%interp%N+1 ! Loop over quadrature points
+  !               do i = 1,this%solution%interp%N+1 ! Loop over quadrature points
+  !                 x = this%geometry%x%boundary(i,j,k,iEl,1,1:3)
+
+  !                 this%solution%extBoundary(i,j,k,iEl,1:this%nvar) = &
+  !                   this%hbc3D_Prescribed(x,this%t)
+  !               enddo
+  !             enddo
+
+  !           endif
+  !         endif
+
+  !       enddo
+  !     enddo
+
+  !     call gpuCheck(hipMemcpy(this%solution%extBoundary_gpu, &
+  !                             c_loc(this%solution%extBoundary), &
+  !                             sizeof(this%solution%extBoundary), &
+  !                             hipMemcpyHostToDevice))
+  !   endif
+  !   call setboundarycondition_CIGLES3D_gpu(this%solution%extboundary_gpu, &
+  !                                               this%solution%boundary_gpu, &
+  !                                               this%mesh%sideInfo_gpu, &
+  !                                               this%geometry%nhat%boundary_gpu, &
+  !                                               this%solution%interp%N, &
+  !                                               this%solution%nelem)
+
+  ! endsubroutine setboundarycondition_CIGLES3D
+
+endmodule self_CIGLES3D