Sorry, ich habe den Vektorpfeil bei \(\vec f\) übersehen...
In diesem Fall wirkt der Operator auf jede einzelne Komponete von \(\vec f\):
$$(\vec r\,\vec\nabla)\vec f=\left(x\frac{\partial}{\partial x}+y\frac{\partial}{\partial y}+z\frac{\partial}{\partial z}\right)\begin{pmatrix}f_x\\f_y\\f_z\end{pmatrix}=\begin{pmatrix}x\frac{\partial f_x}{\partial x}+y\frac{\partial f_x}{\partial y}+z\frac{\partial f_x}{\partial z}\\[1ex]x\frac{\partial f_y}{\partial x}+y\frac{\partial f_y}{\partial y}+z\frac{\partial f_y}{\partial z}\\[1ex]x\frac{\partial f_z}{\partial x}+y\frac{\partial f_z}{\partial y}+z\frac{\partial f_z}{\partial z}\end{pmatrix}$$
Normalerweise lautet ja die Taylor-Reihe für Funktionen \(f(\vec r)\):$$f(\vec r_0+\Delta \vec r)=e^{\vec\Delta r\cdot\vec\nabla}f(\vec r_0)=\sum\limits_{n=0}^\infty\frac{(\Delta\vec r\cdot\vec\nabla)^n}{n!}\,f(\vec r_0)$$
Wenn du diese Entwicklung auf Vektorfunktionen \(\vec f(\vec r)\) überträgst, musst du jede Komponente der Funktion einzeln nähern. Das kann man vektoriell zusammenfassen:$$\vec f(\vec r_0+\Delta \vec r)=e^{\vec\Delta r\cdot\vec\nabla}\vec f(\vec r_0)=\sum\limits_{n=0}^\infty\frac{(\Delta\vec r\cdot\vec\nabla)^n}{n!}\,\vec f(\vec r_0)$$
Dabei taucht dann für \(n=1\) "dein" Ausdruck von oben auf.