From a170e639b5c59a1fd3334cad36ef094c3455c285 Mon Sep 17 00:00:00 2001 From: dinatraykova Date: Sat, 2 Sep 2023 22:54:37 +0200 Subject: [PATCH] Fix oblique shock ID --- Examples/PerfectFluid/ObliqueShockID.hpp | 16 +++------------- .../ObliqueShockID_SimulationParameters.hpp | 2 +- 2 files changed, 4 insertions(+), 14 deletions(-) diff --git a/Examples/PerfectFluid/ObliqueShockID.hpp b/Examples/PerfectFluid/ObliqueShockID.hpp index c06346345..8de3f6cc1 100644 --- a/Examples/PerfectFluid/ObliqueShockID.hpp +++ b/Examples/PerfectFluid/ObliqueShockID.hpp @@ -66,16 +66,6 @@ class InitialFluidData data_t chi_regularised = simd_max(metric_vars.chi, 1e-6); - // Transition function - // template - // data_t transD(data_t d, double delta) - //{ - // return 0.5*(1.+tanh(d/delta)); - //} - - // data_t rho = m_params.emax - m_params.emin*transD(de,delta); - data_t rho = - m_params.emax - m_params.emin * 0.5 * (1. + tanh(de / delta)); x_p_L = pow(x + m_params.L, gamma); x_m_L = pow(x - m_params.L, gamma); y_p_L = pow(y + m_params.L, gamma); @@ -86,14 +76,14 @@ class InitialFluidData dvy = m_params.L - pow(x_m_L + y_p_L, one_o_gamma); dn = m_params.L - pow(x_p_L + y_p_L, one_o_gamma); + data_t rho = + m_params.emax - m_params.emin * 0.5 * (1. + tanh(de / delta)); + Tensor<1, data_t> vi; - // vi[0] = m_params.vx_in*transD(dvx,delta);; vi[0] = m_params.vx_in * 0.5 * (1. + tanh(dvx / delta)); - // vi[1] = m_params.vy_in*transD(dvy,delta);; vi[1] = m_params.vy_in * 0.5 * (1. + tanh(dvy / delta)); vi[2] = 0.; - // data_t nn = m_params.nn_in*transD(dn,delta) + 0.1; data_t nn = m_params.nn_in * 0.5 * (1. + tanh(dn / delta)) + 0.1; data_t eps = 0.; diff --git a/Examples/PerfectFluid/ObliqueShockID_SimulationParameters.hpp b/Examples/PerfectFluid/ObliqueShockID_SimulationParameters.hpp index 5dd55a713..7ee6163c0 100644 --- a/Examples/PerfectFluid/ObliqueShockID_SimulationParameters.hpp +++ b/Examples/PerfectFluid/ObliqueShockID_SimulationParameters.hpp @@ -37,7 +37,7 @@ class SimulationParameters : public SimulationParametersBase pp.load("fluid_emin", initial_params.emin, 2.97); pp.load("fluid_vx_in", initial_params.vx_in, 0.97); pp.load("fluid_vy_in", initial_params.vy_in, 0.97); - pp.load("fluid_nn_in", initial_params.nn_in, 0.2); + pp.load("fluid_nn_in", initial_params.nn_in, 0.5); pp.load("lambda", lambda, 1.); // eigenvalue for numerical flux pp.load("eos_w", eos_params.eos_w, 1. / 3.);