[cucs77]

Hello, I'm a new SU2 user. I use Shock-Fitting Methods, so I Need to add the boundary condition for the shock in Source code. Is It possible? If yes, which subroutine can I modify? Thanks.

[bigfootedrockmidget]

I don't understand what you are asking. The shocks will be present as a solution of the Euler equations, and we have shock capturing schemes available, like MUSCL schemes.
What is it that you want to achieve?

[cucs77]

In the shock-fitting technique the shock wave are considered as internal boundaries in the flow field. In particular in 3D, they are surfaces inside the field and their displacement is caused by the application of Jump relations. So, I have to impose a new boundary condition for the shock inside the source code. I don't know If Is possible.

[bigfootedrockmidget]

I had to look this up, it does not appear to be a widely used method anymore. This is Moretti's method, where you solve the Rankine-Hugoniot relations?
If it it is 'true' shock fitting it seems you need to track the location of the shock somehow, and possibly also modify the computational mesh. The best you can do in SU2 right now without a PhD effort is probably trying to find the location of the shock using the current mesh and move mesh points towards that location.

For using the RH relations, do you typically solve the RH relations for each cell or cell face and identify the shock in that way?

[cucs77]

It's an evolution of Moretti's technique. I developed a Shock-fitting code that works with a shock-capturing solver. It Is an hybrid approach. So, I don't need to use SU2 for the shock-fitting part. In particular, I have a mesh with a Wall and also with shock (upstream and downstream) surfaces. But I have to impose a new boundary condition. For example for an inviscid Wall I impose MARKER_EULER (...). Now, I need to impose a new bnd for example MARKER_SHOCK for the shock wave. I know that SU2 hasn't this condition in this moment. I'd like to know If there Is a possibility to modify the source code (the One on GitHub) tò add this new boundary condition. Finally, I'd likes to know in which subroutine I can find the code for the boundary condition.

[bigfootedrockmidget]

You can take for instance an existing boundary condition, like marker_symmetry, with the keyword SYMMETRY_PLANE, and have a look with vscode at the locations where they are called and what is being done with them. The implementation of boundary conditions is pretty centralized. so adding a new conditions should be straight-forward. You could also take an existing one like marker_internal, which is a passive marker, and see if you can use it for your purposes. You can couple it with external codes using the python interface, there are examples in Testcases/py_wrapper.

[cucs77]

Thank you very much.