In this paper we study domain walls of \(\mathcal {N}=1\)  \(4d\)  massive SQCD that preserve half of the supersymmetry, these are called BPS domain walls [1]}, [2]}, [3]}. These walls interpolate between two vacuum configuration at the two ends of the Universe, at \(x=\pm \infty \) . When the low energy physics of the model is described by a Wess-Zumino model with superpotential \(W\) , the tension of the walls is fixed by supersymmetry and it is given by
\(T=2{W(v_i)-W(v_j)},\) 
