[ander]

Dear SU2 devs,

I'm trying to modify the source code to account for time varying far field boundary conditions for the Euler Equations. (In the future I would also like to extend them to vary spatially, but this is a first step).

I have modified the function CEulerSolver::BC_Far_Field, and the edited piece of the code looks like this:

      bool UsingTraject = config->UsingTrajectory();
      Vel2_Infty = 0.0; Vn_Infty = 0.0;
      for (iDim = 0; iDim < nDim; iDim++) {
          Vel_Infty[iDim] = UsingTraject ? config->GetVelocityFS_Traject(iDim) : 
          GetVelocity_Inf(iDim);
          Vel2_Infty     += Vel_Infty[iDim]*Vel_Infty[iDim];
          Vn_Infty       += Vel_Infty[iDim]*UnitNormal[iDim];
      }
      Density_Infty = UsingTraject ? config->GetDensityFS_Traject() : GetDensity_Inf();
      Pressure_Infty   = UsingTraject ? config->GetPressureFS_Traject() : 
          GetPressure_Inf();
      SoundSpeed_Infty = sqrt(Gamma*Pressure_Infty/Density_Infty);
      Entropy_Infty    = pow(Density_Infty, Gamma)/Pressure_Infty;

Here I have added some additionals values and functions to CConfig. The bool UsingTraject specifies if the time varying conditions are used. The functions Get_XXX_FS_Traject() reads in time varying free stream values from the config file, and these free stream values are updated once each timestep.

However my implementation doesn't seem to work. I have ran a supersonic case on a wedge-like object where I gradually changed the direction of the freestream velocity, but the incoming field doesn't seem to change direction.

I would like some guidance on possible reasons for why my implementation is failing. Is there is something wrong in my assumption that I can modify the BC_Far_Field function in this way?

Thanks in advance!