Skip to content
Snippets Groups Projects
Verified Commit 0f70e544 authored by Florian Atteneder's avatar Florian Atteneder
Browse files

rerun euler tests on latest version

parent 70f8fcff
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
Showing
with 14 additions and 105 deletions
eulereq_blast_wave/mono1_avmda/dirty
+ 1
− 1
View file @ 0f70e544
git commit hash: d20501b2654ae83d58570d05b14b04c3b58ccd8f
git commit hash: 9cbd54efc3eb0d92564b3bb207e56bf0a36d6b70
====================================================================================================
eulereq_blast_wave/mono2_avmda/dirty
+ 1
− 1
View file @ 0f70e544
git commit hash: d20501b2654ae83d58570d05b14b04c3b58ccd8f
git commit hash: 9cbd54efc3eb0d92564b3bb207e56bf0a36d6b70
====================================================================================================
eulereq_blast_wave/mono3_avmda/dirty
+ 1
− 1
View file @ 0f70e544
git commit hash: d20501b2654ae83d58570d05b14b04c3b58ccd8f
git commit hash: 9cbd54efc3eb0d92564b3bb207e56bf0a36d6b70
====================================================================================================
eulereq_blast_wave/mono4_avmda/dirty
+ 1
− 1
View file @ 0f70e544
git commit hash: d20501b2654ae83d58570d05b14b04c3b58ccd8f
git commit hash: 9cbd54efc3eb0d92564b3bb207e56bf0a36d6b70
====================================================================================================
eulereq_blast_wave/plot_dg_ts_combined.jl
+ 1
− 1
View file @ 0f70e544
......@@ -135,7 +135,7 @@ function make_plot(rootdir::String, dirs::Vector{String}, refsol_dir::String, ts
ax.yticks[] = 0.0:5.0:20.0
end
for ax in axs_μ
ylims!(ax, (-0.01, 0.11))
ylims!(ax, (-0.01, 0.07))
ax.yticks[] = 0.00:0.02:0.1
end
......
eulereq_blast_wave/ref_fv/dirty
+ 1
− 1
View file @ 0f70e544
git commit hash: b81c9c3af74bdb658dabd3202368e6f6fd441e93
git commit hash: 9cbd54efc3eb0d92564b3bb207e56bf0a36d6b70
====================================================================================================
eulereq_sod_shock_tube/mono1_avmda/dirty
+ 1
− 1
View file @ 0f70e544
git commit hash: d44f736f7a24257b4f4372fe141431269cc2cc47
git commit hash: 9cbd54efc3eb0d92564b3bb207e56bf0a36d6b70
====================================================================================================
eulereq_sod_shock_tube/mono2_avmda/dirty
+ 1
− 1
View file @ 0f70e544
git commit hash: d44f736f7a24257b4f4372fe141431269cc2cc47
git commit hash: 9cbd54efc3eb0d92564b3bb207e56bf0a36d6b70
====================================================================================================
eulereq_sod_shock_tube/mono3_avmda/dirty
+ 1
− 1
View file @ 0f70e544
git commit hash: d44f736f7a24257b4f4372fe141431269cc2cc47
git commit hash: 9cbd54efc3eb0d92564b3bb207e56bf0a36d6b70
====================================================================================================
eulereq_sod_shock_tube/mono4_avmda/dirty
+ 1
− 1
View file @ 0f70e544
git commit hash: d44f736f7a24257b4f4372fe141431269cc2cc47
git commit hash: 9cbd54efc3eb0d92564b3bb207e56bf0a36d6b70
====================================================================================================
eulereq_sod_shock_tube/plot_dg_ts_combined.jl
+ 2
− 2
View file @ 0f70e544
......@@ -133,8 +133,8 @@ function make_plot(rootdir::String, dirs::Vector{String}, refsol_dir::String, ts
ax.yticks[] = 0.0:0.2:1.0
end
for ax in axs_μ
ylims!(ax, (-0.001, 0.007))
ax.yticks[] = 0.000:0.002:0.006
ylims!(ax, (-0.0005, 0.0035))
ax.yticks[] = 0.000:0.001:0.003
end
return fig
......
eulereq_sod_shock_tube/ref_fv/dirty
+ 2
− 93
View file @ 0f70e544
git commit hash: d44f736f7a24257b4f4372fe141431269cc2cc47
git commit hash: 9cbd54efc3eb0d92564b3bb207e56bf0a36d6b70
====================================================================================================
diff --git a/src/EulerEq/EulerEq.jl b/src/EulerEq/EulerEq.jl
index 66c12fc..9800583 100644
--- a/src/EulerEq/EulerEq.jl
+++ b/src/EulerEq/EulerEq.jl
@@ -54,7 +54,7 @@ dg1d.@parameters EulerEq begin
"""
Boundary conditions. Available options:
- `periodic`
- - `from_id`
+ - `from_id`, assumes that any dynamically generated features in the bulk never reach the bdry
"""
bc = "periodic"
@check bc in [ "periodic", "from_id" ]
diff --git a/src/EulerEq/rhs.jl b/src/EulerEq/rhs.jl
index 3d0681a..6920c55 100644
--- a/src/EulerEq/rhs.jl
+++ b/src/EulerEq/rhs.jl
@@ -168,7 +168,7 @@ function step_fv_muscl!(P::Project, mesh::Mesh1d{FVElement}, state)
v = get_static_variables(mesh)
@unpack eos = equation
- periodic = isperiodic(mesh)
+ bc = P.prms.bc
prms_slope_limiter = (; method=P.prms.slope_limiter_method,
tvb_M=P.prms.slope_limiter_tvb_M)
K, = mesh.tree.dims
@@ -193,14 +193,12 @@ function step_fv_muscl!(P::Project, mesh::Mesh1d{FVElement}, state)
v₋, v₀, v₊ = v[k₋], v[k], v[k₊]
dv₊[k], dv₋[k] = v₊-v₀, v₀-v₋
end
- if periodic
+ if bc === :periodic
for k in (1,K)
k₋, k₊ = periodic_index(k-1,K), periodic_index(k+1,K)
v₋, v₀, v₊ = v[k₋], v[k], v[k₊]
dv₊[k], dv₋[k] = v₊-v₀, v₀-v₋
end
- else
- TODO()
end
# limit slope and compute left (L) and (R) interface values
@@ -259,7 +257,7 @@ function step_fv_muscl!(P::Project, mesh::Mesh1d{FVElement}, state)
end
# bdry conditions
- if periodic
+ if bc === :periodic
begin
k = 1
nx = -1
@@ -294,8 +292,6 @@ function step_fv_muscl!(P::Project, mesh::Mesh1d{FVElement}, state)
f₊, f₋ = (E₊+p₊)*q₊/rho₊, (E₋+p₋)*q₋/rho₋
nflx_E_R[k] = LLF(nx*f₊, nx*f₋, E₊, E₋, max_v)
end
- else
- TODO()
end
# evaluate rhs
@@ -306,6 +302,17 @@ function step_fv_muscl!(P::Project, mesh::Mesh1d{FVElement}, state)
rhs_E[k] = -(nflx_E_R[k] + nflx_E_L[k] )/dl
end
+ if bc === :from_id
+ # force rhs to be zero to preserve bdry values
+ # assuems that data never reaches bdry
+ rhs_rho[1] = 0
+ rhs_q[1] = 0
+ rhs_E[1] = 0
+ rhs_rho[K] = 0
+ rhs_q[K] = 0
+ rhs_E[K] = 0
+ end
+
return
end
diff --git a/src/EulerEq/setup.jl b/src/EulerEq/setup.jl
index cc051f0..9b1ff30 100644
--- a/src/EulerEq/setup.jl
+++ b/src/EulerEq/setup.jl
@@ -11,8 +11,9 @@ function Project(env::Environment, mesh::Mesh1d, prms)
id_smooth = prms["EulerEq"]["id_smooth"]
av_drag = prms["EulerEq"]["av_drag"]
av_recompute_substeps = prms["EulerEq"]["av_recompute_substeps"]
+ bc = Symbol(prms["EulerEq"]["bc"])
fixedprms = (; av_regularization, bernstein, hrsc, slope_limiter_method, slope_limiter_tvb_M,
- id_smooth, av_drag, av_recompute_substeps)
+ id_smooth, av_drag, av_recompute_substeps, bc)
Pρ = prms["EulerEq"]["av_navierstokes_prho"]
PT = prms["EulerEq"]["av_navierstokes_pt"]
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment