Let's say that your surface is defined by $\vec{r}(\Theta) = [r_0 + s f(\Theta) ] \hat{r}$, where $\Theta$ represents the spherical angles. I'm going to conjecture that the normal vector in radial coordinates is proportional to $\vec{n} = (-s\nabla f,1)$, where the gradient is a 2-vector in the angle space, and the last component is the radial direction.
We can verify this by making sure the normal vector is orthogonal to the surface. Note that a tangent vector to the surface is $(1,s\nabla f)$. This tangent points along the gradient direction. A normal vector orthogonal to this one points along an isocontour, which by construction is $(\hat{\Theta}^\perp,0)$, where $\hat{\Theta}^\perp$ is 90 degree rotation of the unit angular vectors in the tangent plane to the sphere passing through the point. The dot product with both of these is zero; the first is obvious, while for the second, $\nabla f \cdot \hat{\Theta}^\perp = 0$.
TL;DR: In radial coordinates: $(-s\nabla f,1)$, then just normalize and convert to cartesian if needed.
Elaboration: Based on your comment, it seems that your $f(\vec{x})$ is a scalar function of a 3-vector. Then $\nabla f$ is a 3-vector. The notation I use lumps the spherical coordinates $(\phi,\theta)$ into a single abstract 2-vector-like quantity $\Theta$, so that when you say $f(\Theta)$, that means $f$ is a function of only the spherical coordinate angles (and not a function of radius). Then $\nabla f(\Theta)=\nabla_\Theta f(\Theta)$ is a gradient in this 2-dimensional angle space. Your $\nabla f$ needs to be projected onto the surface of a sphere first, since you only ever sample $f(\vec{x})$ for $\lVert \vec{x} \rVert = 1$.
So, to calculate your normals, given a point $\vec{P}(\vec{x}) = [R + s f(\vec{x})]\vec{x}$ with $\vec{x}$ such that $\lVert \vec{x} \rVert = 1$,
- Let $\vec{g}(\vec{x}) = \nabla f(\vec{x})$. Here, $\vec{g}$ is the true gradient of $f$ in 3-space, and you can calculate that in cartesian coordinates, so you get $\vec{g} = (g_x,g_y,g_z)$.
- Project out the radial component of $\vec{g}$ to get $\vec{h}(\vec{x})$. To do this, $\vec{h} = \vec{g} - \frac{\vec{g}\cdot\vec{x}}{\vec{x}\cdot\vec{x}}\vec{x}$. Note that the denominator should be 1. Here, $\vec{h}$ represents the component of $\vec{g}$ that should be tangential to a sphere centered at the origin and passing through $\vec{P}$.
- An outward normal vector to the surface is $\vec{n} = \vec{x} - s\cdot \vec{h}$. Normalize it to get a unit normal vector.
All these computations can be done in Cartesian coordinates, but notice that I never had to resort to referring to the Cartesian components of any vector; the basic operations are vector arithmetic and dot products.
If you have a surface embedded in a Euclidean space (for the sake of example, the radius 1 sphere of $R^3$ which we call $S^2$), there are a few things to disambiguate.
NB: In math, the "sphere" is the peel, the "ball" is the inside of the orange. We're talking only about the surface: the standard equation for an embedding of $S^2$ in $R^3$ is $S^2 = \{x^2 + y^2 + z^2 = 1 \space | \space (x,y,z) \in R^3\}$, not $B^3 = \{x^2 + y^2 + z^2 \leq 1 \space | \space (x,y,z) \in R^3\}$.
So let's start simple. We're in $R^3$, there is no sphere. Only a a scalar field over $R^3$. A scalar field can be defined as a function $f : R^n \to R$, here $n = 3$. Visually, it's like giving a color to each point of $R^3$, with:
- darker, redder colors for points of $R^3$ mapped to a value close to $+\infty$
- white for points of $R^3$ mapped to zero
- darker, bluer colors for points of $R^3$ mapped to a value close to $-\infty$
Let's take a function which "generates" the radius $1$ sphere, meaning $f(x,y,z) = x^2 + y^2 + z^2 - 1$. With this function, points inside $S^2$ are blue (but rather clear blue as the minimum is given by $f(0,0,0) = -1$, light blue), and the pure white points are $S^2$ itself. Use a grapher like Geogebra to see that this is the case, by testing $f(x,y,z) = 0$, which should be the level set corresponding to the sphere. Swap $0$ with another constant $c$, and you should get other spheres with other radii.
Your level sets (sets of points that map to the same scalar value) for this function $f$ are "each sphere centered in $0$", which is equivalent here to "the set of points for a fixed $c$ for $f(x,y,z) = c$", and also equivalent "each set of points of the same color".
In this context, your gradient (no matter the constant $c$, or equivalently the constant $C > 0$ in $f_C(x,y,z) = x^2 + y^2 + z^2 - C$) will be $\nabla f = (2x, 2y, 2z)$, by deriving in each independent direction.
Now ask yourself, what does this correspond to in this picture ? For each point $p = (x,y,z)$, $\nabla f(p)$ gives a vector that goes from the origin and passes through $p$, with double the length of the origin to $p$. These "rays" are precisely the vectors which are normal to your sphere (if you consider them starting from $p$, on each sphere, rather than from the origin).
[If you were to use another scalar function, generating another manifold (or rather, level sets each which is a manifold, and which is non-intersecting with the others), and calculated the gradient, you would find the same result, where the gradient vector at $p$ is normal to the manifold corresponding to the level set of the value $f(p)$. The function $f(x,y,z) = x^2 + y^2 - z^2 - c$ for the cone ($c = 0$) and hyperbolae ($c \neq 0$) is a highly instructive example, as well (especially if pseudo-euclidean metrics and spacetime are an objective of yours).]
In this context, you can see that the gradient describes "the direction of greatest increase". The direction a point "p" must evolve to get to go from "dark blue to dark red" the fastest possible. Since for every point $p$ you can define a gradient vector $\nabla f(p)$ for that point, the gradient operator actually turns a scalar field $f$, into a vector field $\nabla f$. Important point: this means that $\nabla f : R^n \to R^n$: it takes a point from our base euclidean space, and returns a vector of same dimension.
But what if we were to consider a scalar field, and its consequent gradient field, on the sphere $S^2$ itself ? Ie, a map $g : S^2 \to R$ ? Why not generalize scalar functions to all input manifolds ?
In this context, the rest of $R^3$ is completely ignored: if you look into the difference between "extrinsic geometry" and "intrinsic geometry", you'll see that we could technically represent our sphere as a distorted 2D map (like the famous Mercator projection, or others like stereographic, etc.). Geometry still works, only you need to keep in mind that some things change (sometimes it'll be angles, sometimes lengths, etc.: it actually depends on the type of projection you use, and how it removes curvature from the manifold to make it flat-like).
If we want to use the sphere, which is curved, embedded in $R^3$, and define vectors on the sphere, we need to define the "tangent bundle". Let's unpack the concept.
Basically, if we define a scalar field on the sphere, it's easy: each point on the sphere is given a color corresponding to its value, from dark blue (very negative) to dark red (very positive). Think of temperatures on a weather map, or on a globe: really simple.
For a vector field, it's bit more ambiguous. A vector starting from my sphere "goes out from the sphere into $R^3$", because it is a straight arrow, not a curved arrow. But the "winds" on your sphere model stay on your model, and don't go out into the surrounding space, right ? So what gives ?
Well, we'll decide that for you manifold M of dimension $n$ (here, $n = 2$), at each point, we'll have a version of $R^n$ that is tangent to the manifold at that point, noted $T_p(M)$. The set of all $T_p(M), p \in M$, is called the tangent bundle of M, and is noted $T(M)$. In the case of $S^2$, you can imagine each $T_p(S^2)$ as a plane tangent to the point $p$. The tangent bundle $T(M)$ is like a "hairy sphere", but instead of each "hair" being a vector, each hair is replaced by a tangent plane.
[What's very interesting, and hard to visualize, is that you can generally turn T(M) into a manifold of higher dimension with nice properties. The wikipedia article on tangent bundles gives the only visualization that I can provide: https://en.wikipedia.org/wiki/Tangent_bundle
This "turning the circle $S^1$ and its tangent lines into a cylinder" is very useful. Why ? Because a smooth vector field over $S^1$ is precisely an image (output) of the circle vertically on the cylinder (like a rubber band around the cylinder). The circle on the cylinder itself is precisely graphed by the scalar function $f(\theta) = 0, \theta \in [0, 2\pi[$. And you've been using such a scheme for years: the curves of continuous functions of $R \to R$, that you know very, very well, are precisely vector fields on $R$, and the space $R^2$ in which you graph/represent these curves is the tangent bundle $T(R)$. The function $f(x) = 0$ is graphed precisely to the line $R$ itself.
Do note that vectors on a 1D space and scalars on a 1D space are the same thing, so scalar field and vector field are indistinguishable in this example, as the output space is $R^1$ in both. There is a distinction to be made as soon as you start from a 2D space.
Sadly, $T(S^2)$ is 4-dimensional (and parallelizability isn't verified for $S^2$ because of the "hairy ball theorem", which is very problematic), so we have to stick with our hairy-sphere-but-hairs-are-tangent-planes image.]
Say we define the function $g(\theta, \rho) = \theta$ where $\theta$ is latitude on the sphere, ie, $\theta = - \pi$ at the north pole, and $\theta = + \pi$ at the south pole. Your north pole is the bluest point on your sphere, the south pole is the reddest point, and the circles of fixed latitude are your level sets. The equator is the white circle. Then, $\nabla g$ is the vector field over the sphere that has two singularities, one at the south pole and one at the north pole, and where the arrows all follow the meridians from the north pole to the south pole.
In this context, your gradient field is made up of vectors on the tangent bundle - ie, at each point $p$, the vector $\nabla g(p)$ exists in $T_p(S^2)$.
In a nutshell, this is the distinction one needs to be very clear about: what is the nature of the manifold you use as your input space (especially in the case of an embedding) ?
As an added bonus: do you know that all function spaces are vector spaces ? If so, then you understand that since vector fields from an $n$-manifold $M^n$ to $R^n$ are functions, they form a function space $(M^n \to R^n)$ that behaves like a vector space. As in, you can add and scale elements of $(M^n \to R^n)$, which are vector fields, like you would vectors. You can also define a multiplication between the vector fields (the "commutator") on this function space, turning it into what is called a Lie algebra. Lie algebras plays a fundamental role in physics and the theory of differential equations over manifolds.
Other added bonus: look up the notion of Hodge dual/Hodge star: in the context of an $(n-1)$-space $M$ embedded in an $n$-space $E$, the normal at $p \in M$ is dual to the tangent plane at $p$.
Best Answer
The surface gradient $\nabla_s f$ has 3 components (for a 2-dimensional hypersurface in $\mathbb{R}^3$). Let's call them $$ \nabla_s f = ( \underline{D}_1 f, \underline{D}_2 f, \underline{D}_3 f ). $$ Then the surface Laplacian is simply: $$ \Delta_s f = \sum_{i=1}^3 \underline{D}_i \underline{D}_i f. $$