From 04271cac347f5765f179b0fec2388cb29dc40d93 Mon Sep 17 00:00:00 2001
From: Florian Atteneder <florian.atteneder@uni-jena.de>
Date: Thu, 25 Apr 2024 18:12:08 +0200
Subject: [PATCH] implement bdry conditions for FV
---
src/EulerEq/rhs.jl | 10 ++++++++--
1 file changed, 8 insertions(+), 2 deletions(-)
diff --git a/src/EulerEq/rhs.jl b/src/EulerEq/rhs.jl
index 3d0681a5..0c315f62 100644
--- a/src/EulerEq/rhs.jl
+++ b/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
--
GitLab