I am postgraduate student, and in my project I am working with a function in several variables. I would be surprised if it has not been studied extensively, but as I am not a mathematician, I would like to consult with you to find out how I can learn more about it.
Suppose $\vec{y}_1,\vec{y}_2,\cdots,\vec{y}_m \in \mathbb{R}^d$ are fixed, and $a_i,b_j > 0$ for all $i,j$ then the function is $f:\mathbb{R}^{n \cdot d} \to \mathbb{R}$ given by $$ f(\vec{x}_1,\vec{x}_2\cdots,\vec{x}_n) = \sum_{i=1}^n \sum_{j=1}^m \frac{a_j}{|\vec{x}_i - \vec{y}_j|} + \sum_{i\neq j} \frac{b_j}{|\vec{x}_i - \vec{x}_j |}. $$ Here $|\cdot|$ is the Euclidean norm. For a given collection $\{\vec{y}_1,\vec{y}_2,\cdots,\vec{y}_m\}$, I want to find out as much as possible about the inverse images $f^{-1}(c)$ for each real number $c$.
I am not very experienced, so my question is basically this: How would you search for information for this problem? 1. In which field does this question belong? 2. What would you call an equation of this form?

Do you know about Morse functions? One can say that $f^{-1}(c)$ is a smooth manifold for all but finitely many values of $c$. If you vary $c$ the $f^{-1}(c)$ will by homotopy equivalent as long as you don't hit any of the points where it is not a smooth manifold. At these values, the preimage will have a singularity and the topology of the preimage will change.
In a simple example with say $d=2$ and $n=1$ one can see what is going on, the general picture is similar but hard to visualize due to higher dimensions.
As suggested by william in comments, the book "Morse Theory" by John Milnor is considered a good source and introduction to the area.