if $p=q+\varepsilon v$, where $v=\dfrac{\nabla f(q)}{\|\nabla f(q)\|}$.
$f(q)=0$. $\nabla f(x) $ is not $0$ on the domain.
How to prove there exists $\varepsilon_1>0$ such that if $\varepsilon_1>\varepsilon>0$, then $f(p)>0$?
I know gradient direction is the fastest increasing direction, so by the condition given, it seems obvious.
How to write it rigorously?
How to write this may depend on what has already been established. Suppose we can show that $$ \left.\frac d {d\varepsilon} f(p) \right|_{\varepsilon=0} >0. \tag 1 $$ Then the problem reduces to this: If $g(0)=0$ and $g'(0)>0$, how do we show that there is some $\varepsilon_1>0$ such that for all $\varepsilon$ between $0$ and $\varepsilon_1$ we have $g(\varepsilon)>0$?
That $g'(0)>0$ means that $$ D=\lim_{\varepsilon\to0}\frac{g(\varepsilon)}{\varepsilon} > 0. $$ So consider the interval $\left(\dfrac D 2,\dfrac{3D} 2\right)$. Certainly there exists $\varepsilon_1>0$ such that for all $\varepsilon$ within a distance of $\varepsilon_1$ from $0$, the fraction $g(\varepsilon)/\varepsilon$ is in the interval $\left(\dfrac D 2,\dfrac{3D} 2\right)$. If $\varepsilon$ is positive, that certainly means $g(\varepsilon)>0$ and so we're done.
The question is then how to prove $(1)$. And there I would try looking at the sort of chain rule that says $$ \frac d {d\varepsilon} f(p) = \nabla f(p)\cdot \frac{dp}{d\varepsilon} \text{ (a dot-product of two vectors)}. $$