Prove the curvature of a level set equals divergence of the normalized gradient

From $\| \gamma'(s) \| = 1$, it follows that $$0 = \frac{d}{ds} \langle \gamma'(s),\gamma'(s) \rangle = 2\langle \gamma''(s),\gamma'(s)\rangle.$$ So $\gamma'$ and $\gamma''$ are perpendicular. Assume that $\gamma''(s) \neq 0$. Thus $N(s) := \gamma''(s)/\|\gamma''(s)\|$ locally defines a unit normal field to $\gamma$. Since $\gamma$ parametrizes a level set of $\phi$ and $\nabla \phi$ is orthogonal to its level sets it follows that $\nabla \phi / \|\nabla \phi\|(\gamma(s)) = \pm N(s)$. Consequently $$\kappa(s) = \|\gamma''(s)\| = |\langle \gamma''(s),N(s)\rangle| = |\frac{d}{ds} \langle \gamma'(s),N(s)\rangle - \langle \gamma'(s),\frac{d}{ds}N(s)\rangle| = |\langle \gamma'(s),\frac{d}{ds} N(s)\rangle|.$$ Since $\left\|N(s)\right\| = 1$ it follows as for $\gamma'$ that $\langle N,\frac{dN}{ds} \rangle = 0$, thus $\frac{dN}{ds}$ is a multiple of $\gamma'$ and $\left\|\frac{dN}{ds}\right\| = |\langle \gamma'(s),\frac{dN}{ds}\rangle|$. Thus $$\kappa(s) = \left\|\frac{dN}{ds}(s)\right\| = \left\|\frac{d}{ds} \frac{\nabla \phi}{\|\nabla \phi\|}\right\|.$$ Thus it remains to see that $|\nabla \cdot N| = \left\|\frac{dN}{ds}\right\|$: For this one needs the following fact: Let $\{b_i\}$ be an orthonormal basis of $\mathbb R^n$. Then $$\nabla \cdot X = \sum_{i = 1}^n \langle DX(b_i),b_i\rangle$$ for any vectorfield $X : \mathbb R^n \to \mathbb R^n$ with differential $DX$. Unfortunately (via google) I did not find a reference for this Edit: Using Riemannian geometry this follows since divergence is defined as the trace of the Levi Civita connection $\nabla$, which on $\mathbb R^n$ is given by $\nabla_XY = DY(X)$ for vector fields $X,Y$. From this we see $$\left\|\frac{dN}{ds}(s)\right\| = \langle \gamma'(s),\frac{dN}{ds}(s)\rangle = \langle \gamma'(s),DN(\gamma'(s))\rangle + \langle N(s),DN(N(s))\rangle = \nabla \cdot N(s).$$ By $DN$ we mean the differential of $N$ (note that strictly speaking this does not exists, as $N$ is only defined along $\gamma$. But we can extend $N$ to a vectorfield defined on an open subset using $N = \pm \nabla \phi/||\nabla \phi||$ as above). The second equality holds since $\langle N,DN(X))\rangle = 1/2\partial_X\langle N,N \rangle = 0$ for any vector field $X$ ($\partial_X$ denotes the directional derivative in direction $X$).

Maybe these calculations are more complicated than they could be, but I don't see a simplification right away.