Suppose a smooth function $f : \mathbb{R}^n \to \mathbb{R}$ satisfies the Helmholtz equation, the PDE $\Delta f + k^2 f = 0$.
A while ago someone showed me a trick: Define a function $g:\mathbb{R}^{n+1} \to \mathbb{R}$ by: $$g(x_1, \ldots, x_n, x_{n+1}) = f(x_1, \ldots, x_n) \cosh (k x_{n+1})$$ It turns out that $\Delta g = 0$, so $g$ (thus, $f$) can be analyzed using tools from harmonic function theory.
My question: Does this trick have a name? Is it well known? Can it be found in books?
Here is a perspective on this method. Suppose we treat $k$ itself as a coordinate, and perform a Fourier transform with respect to it. This converts the Helmholtz equation $(\Delta+k^2)f=0$ to D'Alembert's equation $(\Delta-\partial_0^2)g=0$. (I have taken $x_0$ to be the conjugate of $k$, and $g$ to be the corresponding Fourier transform of $f$.) Upon defining $x_{n+1}\equiv i x_0$, this becomes $(\Delta+\partial_{n+1}^2)g=0$ which is equivalent to the form identified above.
Separating the variable $x_{n+1}$ by writing $g(\mathbf{x},x_{n+1})=f(\mathbf{x})\Phi(x_{n+1})$ then obtains $$\frac{\Delta f}{f}=-\frac{\Phi''}{\Phi}=-k^2$$ where the separation constant is chosen so as to regain the Helmholtz equation. This yields $\Phi(x_{n+1})=A\cosh(kx_{n+1})+B\sinh(kx_{n+1})$, generalizing the case of $\Phi(x_{n+1})=\cosh(kx_{n+1})$ given in the OP.