[aerojoe]

Hello,

I have been modifying the source code of SU2 for some time. I have created necessary classes and functions but I'm having some trouble calculating momentum source terms. It will be an implicit and steady state solution.

Just like the Body Force and Gravitational Force, I would like to apply a constant force in Newton to a certain zone. The problem is I can only have information about the force in Newton unit at each node. I can't figure out how to transform this force into momentum source term in flow_sources.cpp file.

For body forces, an equation as such below is used,

residual[iDim+1] = -Volume * U_i[0] * Body_Force_Vector[iDim] / Force_Ref;

In which the Force_Ref is calculated in CEuler_Solver.cpp by using reference density, reference pressure and consequently reference velocity. And Body_Force_Vector usually involves three components with 0,0,9.81. Can I use the same approach in my case?

Any help is appreciated,
Thank you in advance.