Let me translate this question into the adèlic language, because this question is more naturally posed there. Let $\pi$ be a cuspidal automorphic representation of $\mathrm{GL}_3(\mathbb{A_Q})$ and let $\pi'$ be a cuspidal automorphic representation of $\mathrm{GL}_2(\mathbb{A_Q})$. For vectors $\varphi \in \pi$ and $\varphi' \in \pi'$, we consider the global Eulerian integral
\[\int\limits_{\mathrm{GL}_2(\mathbb{Q}) \backslash \mathrm{GL}_2(\mathbb{A_Q})} \varphi\begin{pmatrix} h & \\ & 1 \end{pmatrix} \varphi'(h) |\det h|^{s - \frac{1}{2}} \, dh.\]
This is precisely the integral $\mathcal{L}(s,F \times u)$ if $\varphi$ is the adèlic lift of $F$ and $\varphi'$ is the adèlic lift of $u$.

The vector $\varphi$ has the Fourier expansion
\[\varphi(g) = \sum_{\gamma \in N_2(\mathbb{Q}) \backslash \mathrm{GL}_2(\mathbb{Q})} W_{\varphi} \left(\begin{pmatrix} \gamma & \\ & 1 \end{pmatrix} g\right),\]
where

\[W_{\varphi}(g) = \int\limits_{N_3(\mathbb{Q}) \backslash N_3(\mathbb{A_Q})} \varphi(ng) \psi^{-1}(n) \, dn.\]
By unfolding, we find that the global Eulerian integral is equal to
\[\int\limits_{N_2(\mathbb{Q}) \backslash \mathrm{GL}_2(\mathbb{A_Q})} W_{\varphi}\begin{pmatrix} h & \\ & 1 \end{pmatrix} \varphi'(h) |\det h|^{s - \frac{1}{2}} \, dh,\]
which in turn is equal to
\[\int\limits_{N_2(\mathbb{A_Q}) \backslash \mathrm{GL}_2(\mathbb{A_Q})} W_{\varphi}\begin{pmatrix} h & \\ & 1 \end{pmatrix} W'_{\varphi'}(h) |\det h|^{s - \frac{1}{2}} \, dh,\]
where

\[W'_{\varphi'}(h) = \int\limits_{N_2(\mathbb{Q}) \backslash N_2(\mathbb{A_Q})} \varphi(ng) \psi(n) \, dn.\]

The point of this is that if we choose $\varphi$ and $\varphi'$ to be pure tensors $\bigotimes_{p \leq \infty} \varphi_p$ and $\bigotimes_{p \leq \infty} \varphi_p'$, then this global Eulerian integral factorises into the product of local Eulerian integrals of the form

\[\int\limits_{N_2(\mathbb{Q}_p) \backslash \mathrm{GL}_2(\mathbb{Q}_p)} W_{\varphi_p}\begin{pmatrix} h_p & \\ & 1 \end{pmatrix} W'_{\varphi_p'}(h_p) |\det h_p|^{s - \frac{1}{2}} \, dh_p.\]

Now if $\pi_p$ and $\pi_p'$ are both unramified and $\varphi_p$ and $\varphi_p'$ are both chosen to be spherical vectors (and normalised appropriately), then this local Eulerian integral is equal to $L(s, \pi_p \otimes \pi_p')$, the local component of the Rankin-Selberg $L$-function $L(s,\pi \otimes \pi')$ associated to $\pi$ and $\pi'$. So if both $\pi$ and $\pi'$ are unramified everywhere, then the global Eulerian integral is equal to
\[\prod_{p \leq \infty}L(s, \pi_p \otimes \pi_p') = \Lambda(s, \pi \otimes \pi'),\]
the completed Rankin-Selberg $L$-function (i.e. the finite part $L(s,\pi \otimes \pi) = \prod_p L(s, \pi_p \otimes \pi_p')$ multiplied by the infinite part, which is a product of gamma factors). In particular, this is the case if $\varphi$ is the adèlic lift of an even Hecke-Maaß cusp form $F$ of level $1$ on $\mathrm{GL}_3$ and $\varphi'$ is the adèlic lift of an even Hecke-Maaß cusp form $u$ of level $1$ on $\mathrm{GL}_2$.

What if either $\pi_p$ or $\pi_p'$ is ramified? Well, if $p$ is a prime (i.e. $p \neq \infty$), $\pi_p$ is ramified, but $\pi_p'$ is unramified, then Jacquet, Piatetski-Shapiro, and Shalika prove the existence of a vector $\varphi_p$ for which the local Eulerian integral is still equal to $L(s, \pi_p \otimes \pi_p')$ whenever $\varphi_p'$ is equal to the (appropriately normalised) spherical vector. They call $\varphi_p$ the *essential vector*, though a better name would be *local newform*, since this is the generalisation to $\mathrm{GL}_n$ of the $\mathrm{GL}_2$-notion of a newform. In particular, it is invariant under a congruence subgroup.

What about if $\pi_p'$ is ramified? In this case, if we take $\varphi_p$ and $\varphi_p'$ both to be local newforms, then the local Eulerian integral vanishes! So these are not the right choice of vectors. In fact, very little is known about "natural" choices of vectors for which the local Eulerian integral is equal to the local $L$-function.

In particular, let's take $p = \infty$ and $\pi_{\infty}'$ to be the ramified principal series representation $\mathrm{sgn} |\cdot|^{it} \boxplus \mathrm{sgn} |\cdot|^{-it}$ of $\mathrm{GL}_2(\mathbb{R})$; this is the component at infinity of an odd Hecke-Maaß cusp form on $\mathrm{GL}_2$ with Laplacian eigenvalue $1/4 + t^2$, while $\pi_{\infty}$ will be taken to be a spherical representation of $\mathrm{GL}_3(\mathbb{R})$, so that it is the local component at infinity of an even Hecke-Maaß cusp form on $\mathrm{GL}_3$. Then your question boils down to the following:

Do there exist vectors $\varphi_{\infty} \in \pi_{\infty}$ and $\varphi_{\infty}' \in \pi_{\infty}'$ such that the local Eulerian integral

\[\int\limits_{N_2(\mathbb{R}) \backslash \mathrm{GL}_2(\mathbb{R})} W_{\varphi_{\infty}}\begin{pmatrix} h_{\infty} & \\ & 1 \end{pmatrix} W'_{\varphi_{\infty}'}(h_{\infty}) |\det h_{\infty}|^{s - \frac{1}{2}} \, dh_{\infty}\]
is equal to $L(s, \pi_{\infty} \otimes \pi_{\infty}')$?

Jacquet has shown that such vectors do exist (and that the vectors are $K$-finite, where $K$ is the maximal compact subgroup of the relevant reductive group), but the proof is completely nonconstructive. Moreover, there is no reason to expect uniqueness.

Very recently, Hirano, Ishii, and Miyazaki have announced that they have discovered an explicit description of a pair of vectors satisfying this. The description is unfortunately rather complicated (it involves the decomposition of $\pi_{\infty}$ and $\pi_{\infty'}$ into $K$-types, and choosing explicit vectors in particular $K$-types). Moreover, they don't explicitly state the relation between their vectors and the "natural" choice of vectors; here the "natural" vector is the spherical vector for $\pi_{\infty}$ and the unique vector (up to scalar) lying in the minimal $K$-type of $\pi_{\infty'} = \mathrm{sgn} |\cdot|^{it} \boxplus \mathrm{sgn} |\cdot|^{-it}$, which is the determinant representation of $K = \mathrm{O}(2)$ (which is one-dimensional). Finally, their result must be taken with a pinch of salt, since the proofs are yet to appear.