I'm confused by the following statement from the Evans text on Partial Differential Equations:
If $U \subseteq \mathbb{R}^{n}$ is open, and $u \in C^{1}(U)$ (i.e., $u$ is continuously differentiable on $U$), then given any $\phi \in C^{\infty}_{C}(U)$ ($\phi$ infinitely differentiable with compact support in $U$), we have via integration by parts:
$$\int \limits_{U} u \phi_{x_{i}}(x) \,dx = - \int \limits_{U} u_{x_{i}} \phi(x) \,dx$$
How does this follow via integration by parts? If we let $u = u(x)$ and $dv = \phi_{x_{i}}(x) \,dx$, then we get:
$$\int \limits_{U} u\,dv = uv - \int \limits_{U} v \,du$$
Applying our substitutions, we get:
$$ \int \limits_{U} u(x) \phi_{x_{i}}(x) \,dx = u(x) \phi(x) - \int \limits_{U} \phi(x) u_{x_{i}}(x) \,dx $$
What happens to the $u(x)\phi(x)$ in the formula? Why is it $0$? I understand that near the boundary of $U$, we get that $\phi(x)$ is $0$, but I don't know how to rigorously show anything.
First observe that the "partial integration" formula that you are using strictly speaking does only apply in the case $n=1$. In the case $n>1$, the "substitute" for the partial integration rule that is ordinarily invoked to prove your formula is the divergence theorem, which states $$ \int_{U}{\rm div}f\, dx=\int_{\partial U} \langle f(x), \nu(x) \rangle \, dS(x) \qquad (\ddagger) $$ for $f\in C^{1}\left(\overline{U};\mathbb{R}^{n}\right)$ and $U\subset\mathbb{R}^{n}$ open and bounded with $C^{1}$-boundary (the exact prerequisites may vary). Here, $\nu(x)$ is the outward pointing unit normal vector and $S$ denotes the surface measure. This is then applied to $$ f=\left(\begin{matrix}0\\ \vdots\\ u\phi\\ 0\\ \vdots\\ 0 \end{matrix}\right), $$ where the $u\phi$ is in the $i$th row. Observe that the right hand side of $(\ddagger)$ vanishes because of $\phi\in C_{c}^{\infty}\left(U\right)$. Using the product rule to compute ${\rm div}\,f$ (see below), this implies your claim.
In your case, this can be applied with $U = B_R(0)$, where $R$ is chosen so large that ${\mathrm supp} \, \phi \subset B_R(0)$. Using the compact support of $\phi$ it is easy to see that even $u\phi \in C^1(\Bbb{R}^n)$, so that everything works fine.
This approach has the drawback that one has to know all the machinery of the divergence theorem (open sets with $C^1$ boundaries, the divergence, the outward pointing unit normal field, integration on manifolds and the divergence theorem itself).
Hence, I propose to give more love to the following Lemma which only uses
Lemma: Let $\emptyset\neq V\subset\mathbb{R}^{d}$ be open and bounded. Let $\varphi\in C\left(\overline{V}\right)$ with $\varphi|_{\partial V}\equiv c$ for some $c\in\mathbb{R}$. Furthermore, let $i\in\left\{ 1,\dots,d\right\} $ and assume that $\partial_{i}\varphi\in C\left(V\right)\cap L^{1}\left(V\right)$. Then $$ \int_{V}\left(\partial_{i}\varphi\right)\left(x\right)\, dx=0. $$ Proof: Without loss of generality assume $i=d$ (the proof is essentially the same in all other cases). Because of $\partial_{i}\varphi\in L^{1}\left(V\right)$, we can apply Fubini's theorem to conclude $$ (\dagger) \qquad \int_{V}\left(\partial_{i}\varphi\right)\left(x\right)\, dx=\int_{\mathbb{R}^{d-1}}\int_{\mathbb{R}}\left(\partial_{d}\varphi\right)\left(x_{1},\dots,x_{d-1},x_{d}\right)\cdot\chi_{V}\left(x_{1},\dots,x_{d-1},x_{d}\right)\, dx_{d}\, d\left(x_{1},\dots,x_{d-1}\right), $$ where the inner integral converges absolutely for almost all $\left(x_{1},\dots,x_{d-1}\right)\in\mathbb{R}^{d-1}$.
Let us fix some $x'=\left(x_{1},\dots,x_{d-1}\right)\in\mathbb{R}^{d-1}$. Note that the set $$ V_{x'}:=\left\{ x_{d}\in\mathbb{R}\,\mid\,\left(x',x_{d}\right)\in V\right\} =\left(x_{d}\mapsto\left(x',x_{d}\right)\right)^{-1}\left(V\right) $$ is an open subset of $\mathbb{R}$. In the case $V_{x'}\neq\emptyset$, this implies $$ V_{x'}=\biguplus_{n\in\mathbb{N}}\left(a_{n},b_{n}\right) $$ with $-\infty<a_{n}\leq b_{n}<\infty$ for $n\in\mathbb{N}$. The intervals $(a_n, b_n)$ are the components of $V$, optionally "filled up" with $a_n = b_n$ for $n> n_0$ if the number $n_0$ of components is actually finite (this simplifies the notation).
Here, we used the boundedness of $V$. Otherwise, $a_{n}$ or $b_{n}$ could be infinite. Observe that the $a_{n},b_{n}$ do depend on $x'$. What we use here is that each set is the disjoint union of its components. In $\mathbb{R}$, every component of an open set is an open interval. Finally, $\mathbb{R}$ is second countable, so that the union is countable.
Finally, observe that in the case $a_{n}<b_{n}$, we can not have $a_{n}\in V_{x'}$ or $b_{n}\in V_{x'}$ (why?). Hence, $\left(x',a_{n}\right),\left(x',b_{n}\right)\in\overline{V}\setminus V=\partial V$, so that $\varphi\left(x',a_{n}\right)=c=\varphi\left(x',b_{n}\right)$, which implies $$ \int_{a_{n}}^{b_{n}}\frac{\partial\varphi}{\partial x_{d}}\left(x',x_{d}\right)\, dx_{d}=\varphi\left(x',b_{n}\right)-\varphi\left(x',a_{n}\right)=0. $$ Here, we used that $x_d \mapsto \varphi(x', x_d)$ is continuously differentiable on $(a_n ,b_n)$ and extends continuously to $[a_n,b_n]$. In the case $a_{n}=b_{n}$, the integral trivially vanishes.
In summary, we conclude \begin{eqnarray*} & & \int_{\mathbb{R}}\left(\partial_{d}\varphi\right)\left(x_{1},\dots,x_{d-1},x_{d}\right)\cdot\chi_{V}\left(x_{1},\dots,x_{d-1},x_{d}\right)\, dx_{d}\\ & = & \sum_{n\in\mathbb{N}}\int_{a_{n}}^{b_{n}}\frac{\partial\varphi}{\partial x_{d}}\left(x',x_{d}\right)\, dx_{d}\\ & = & 0. \end{eqnarray*} Plugging this into $\left(\dagger\right)$ above yields the claim.$\square$
Now choose some bounded open set $V\subset U$ with ${\rm supp}\left(\phi\right)\subset V$ and apply the above Lemma to $\varphi=u\cdot\phi$ (check that all assumptions are satisfied with $c=0$). Using the product rule, this implies $$ 0=\int_{V}\left(\partial_{i}\varphi\right)\left(x\right)\, dx=\int_{V}\phi\cdot\partial_{i}u+u\cdot\partial_{i}\phi\, dx, $$ which yields your claim.