From 1d180fbbba07853519d77f12866400b5819d7195 Mon Sep 17 00:00:00 2001 From: Daniel_Doehring Date: Tue, 19 Nov 2024 13:52:47 +0100 Subject: [PATCH] style --- .../elixir_navierstokes_viscous_shock.jl | 11 +++++++---- 1 file changed, 7 insertions(+), 4 deletions(-) diff --git a/examples/tree_1d_dgsem/elixir_navierstokes_viscous_shock.jl b/examples/tree_1d_dgsem/elixir_navierstokes_viscous_shock.jl index 5e37683be5d..13e6a7d1d19 100644 --- a/examples/tree_1d_dgsem/elixir_navierstokes_viscous_shock.jl +++ b/examples/tree_1d_dgsem/elixir_navierstokes_viscous_shock.jl @@ -1,9 +1,6 @@ using OrdinaryDiffEq using Trixi -############################################################################### -# semidiscretization of the ideal compressible Navier-Stokes equations - # This is the classic 1D viscous shock wave problem with analytical solution # for a special value of the Prandtl number. # The original references are: @@ -24,12 +21,14 @@ using Trixi # of a Viscous, Head-Conducting, Compressible Gas # [DOI: 10.2514/8.11882](https://doi.org/10.2514/8.11882) # +# # The particular problem considered here is described in # - L. G. Margolin, J. M. Reisner, P. M. Jordan (2017) # Entropy in self-similar shock profiles # [DOI: 10.1016/j.ijnonlinmec.2017.07.003](https://doi.org/10.1016/j.ijnonlinmec.2017.07.003) ### Fixed parameters ### + # Special value for which nonlinear solver can be omitted # Corresponds essentially to fixing the Mach number alpha = 0.5 @@ -45,6 +44,7 @@ v = 1 # Shock speed domain_length = 5.0 ### Derived quantities ### + Ma = 2 / sqrt(3 - gamma) # Mach number for alpha = 0.5 c_0 = v / Ma # Speed of sound ahead of the shock @@ -67,7 +67,7 @@ The version implemented here is described in [DOI: 10.1016/j.ijnonlinmec.2017.07.003](https://doi.org/10.1016/j.ijnonlinmec.2017.07.003) """ function initial_condition_viscous_shock(x, t, equations) - y = x[1] - v * t + y = x[1] - v * t # Translated coordinate chi = chi_of_y(y) w = 1 + 1 / (2 * chi^2) * (1 - sqrt(1 + 2 * chi^2)) @@ -79,6 +79,9 @@ function initial_condition_viscous_shock(x, t, equations) return prim2cons(SVector(rho, u, p), equations) end +############################################################################### +# semidiscretization of the ideal compressible Navier-Stokes equations + equations = CompressibleEulerEquations1D(gamma) equations_parabolic = CompressibleNavierStokesDiffusion1D(equations, mu = mu(), Prandtl = prandtl_number(),