Let's say I have the Diophantine equation $$ x^2+3n^2 = y^2+3z^2. \tag{$\star$} $$ where $n$ is a known integer, and we're trying to determine solutions in integers $x,y,z \ge 1$. Rewrite ($\star$) as $$ x^2 - 3z^2 = y^2 - 3n^2 = k, \tag{$\dagger$} $$ where $k$ is some unknown integer. For simplicity's sake, let's consider the case $n=1$, so we have $$ x^2 - 3z^2 = y^2 - 3(1)^2 = k, \tag{$\Diamond$} $$ Evidently, $(y,1)$ must be the fundamental solution of the equation $U^2-3V^2=k$. Now either $\lvert z \rvert=1$ implying $\lvert x \rvert=\lvert y \rvert$, or $\lvert z \rvert>1$ implying $\lvert x \rvert > \lvert y \rvert$. Let's assume the more interesting latter case.
My question is: can Pell equation machinery — in particular, some form of descent/ascent — be used to solve the original equation, or at least give formulas for $x$ and $z$ in terms of $y$ and $1$ (or some other permutation of the elements)?
I realise we don't know how far “up the ladder” the solution $(x,z)$ is, but I thought there might be some known and effective results in this direction.
EDIT: Also of interest to me would be an answer in the special case where $y$ is an unknown integer with a known parameterization, e.g. $y=ab+cd$ for positive [but possibly unknown] integers $a,b,c,d$.
The complete solution to,
$$x_1^2+mx_2^2 = x_3^2+mx_4^2\tag{1}$$
is given by,
$$(a c + b d m)^2 + m(a d - b c)^2 = (a c - b d m)^2 + m(a d + b c)^2\tag{2}$$
known as the Brahmagupta-Fibonacci identity.
We wish to solve $ad-bc = \pm n$ in the integers for some constant $n$. We can do this as a Pell-like eqn by using the substitution $a,b,c = q^2,\, p+q,\, p-q$ to get,
$$ad-bc = p^2-(d+1)q^2 = \pm n\tag{3}$$
for some arbitrary $d$. For some constant $\pm n$, there is an indefinite number of $d$ that will do, each one of which will lead to a different Pell-like eqn with an infinite number of integer solutions. (Note that this Pell-like eqn doesn't depend on $m$ at all which, in your example, was $m=3$.) The $x_i$ can then be expressed in terms of $p,q$.
More simply, we can just solve $ad-bc = n$ for $d$, let $a=1$, then an explicit formula is,
$$(b^2 c m + b m n + c)^2 + m n^2 = (b^2 c k + b m n - c)^2 + m(2b c + n)^2\tag{4}$$
where one term ($x_2$) remains a constant n regardless of the three free variables $b,c,m$.