implement bdry conditions for FV

src/EulerEq/rhs.jl

@@ -200,7 +200,7 @@ function step_fv_muscl!(P::Project, mesh::Mesh1d{FVElement}, state)
         dv₊[k], dv₋[k] = v₊-v₀, v₀-v₋
       end
     else
-      TODO()
+      # TODO Atm we assume that evolution does not reach bdry
     end
 
     # limit slope and compute left (L) and (R) interface values
@@ -295,7 +295,7 @@ function step_fv_muscl!(P::Project, mesh::Mesh1d{FVElement}, state)
       nflx_E_R[k]   = LLF(nx*f₊, nx*f₋, E₊, E₋, max_v)
     end
   else
-    TODO()
+    # TODO Atm we assume that evolution does not reach bdry
   end
 
   # evaluate rhs
@@ -305,6 +305,12 @@ function step_fv_muscl!(P::Project, mesh::Mesh1d{FVElement}, state)
     rhs_q[k]   = -(nflx_q_R[k]   + nflx_q_L[k]  )/dl
     rhs_E[k]   = -(nflx_E_R[k]   + nflx_E_L[k]  )/dl
   end
+  if !periodic
+    # TODO Atm we assume that evolution does not reach bdry
+    rhs_rho[1] = rhs_rho[end] = 0
+    rhs_q[1]   = rhs_q[end]   = 0
+    rhs_E[1]   = rhs_E[end]   = 0
+  end
 
   return
 end