Fix oblique shock ID

dinatraykova · Sep 2, 2023 · a170e63 · a170e63
1 parent 706134b
commit a170e63
Show file tree

Hide file tree

Showing 2 changed files with 4 additions and 14 deletions.
diff --git 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 <class data_t>
-        // 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
@@ -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.);