As mentioned above, this is from the Bertsimas and Tsitsiklis, and the Phase I approach they are referring to is in Section 3.5. The standard form LP they use is
\begin{array}{ll}
\text{minimize} & c^T x \\
& A x = b \\
& x \geq 0
\end{array}
They assume that $b\geq 0$; if this is not the case, negate the corresponding rows to make it so. (And for simplicity, let's assume $b$ has at least one nonzero value.) The corresponding Phase I problem looks like this:
\begin{array}{ll}
\text{minimize} & \textstyle\sum_i y_i \\
& A x + y = b \\
& x \geq 0, y \geq 0
\end{array}
Now you see why $b\geq 0$ is important: $(x,y)=(0,b)$ constitutes a trivial feasible solution, so that's your starting point for the Phase I method, with an initial objective of $\sum_i b_i>0$. If the optimal value of this Phase I model is zero, then original model is feasible; otherwise, the original model is infeasible.
It is important to read the statement carefully. It is not claiming that an artificial variable will never re-enter the basis if you leave it in the tableau. In fact, it can. If you have the book, look at Example 3.8. In one of the steps, one of the nonbasic artificial variables has a negative reduced cost, making it an entirely valid candidate for selection, depending on the pivoting rule you use. (The example chose a different pivot, though.)
The statement, therefore, is not that it will never re-enter the basis. Rather, the statement is saying that it does not need to re-enter the basis. What we wish to show is that eliminating one of these variables from the tableau does not prevent the simplex method from terminating properly.
To see why this is the case, consider this "restarting" approach to solving the Phase I model:
- Initialize: basis $\mathcal{B}=(y_1,y_2,\dots,y_n)$, $(x,y)=(\vec{0},\vec{1})$.
- Begin the simplex algorithm simplex with the current basis $\mathcal{B}$ and current $(x,y)$.
- If the algorithm terminates before an artificial variable is eliminated:
- If the cost is zero, STOP. The problem is feasible, but extra steps need to be taken to drive the remaining artificial variables out of the basis. See Section 3.5 for details.
- If the cost is positive, STOP. the problem is infeasible.
- As soon as an artificial variable is removed from the basis, stop the simplex method, and obtain the current basis $\mathcal{B}$ and current point $(x,y)$. (The first time through, this will happen on the very first pivot. In subsequent iterations, it might take longer.)
- Remove the newly nonbasic artificial variable from the problem. If no basic artificial variables remain, STOP. Your problem is feasible and you have a basic feasible solution.
- Repeat steps 2-6.
The key to the success of this approach is noting that, after removing a nonbasic artificial variable in Step 5, you're still left with a valid pseudo-Phase I model and a feasible point for that model. The model still has an optimal value of $0$ if the model is feasible. It's also important to note that objective value reached in Step 4 is going to be the initial objective value at Step 2. So we do not lose progress with our restarts.
In fact, one way to look at it is this: we're not really restarting the simplex algorithm per se. Rather, we are just choosing to recompute the inverse basis matrix $B^{-1}$ from scratch every time an artificial variable exits the basis. Still, by expressing the approach in this way, it is hopefully easier to see why the elimination of nonbasic artificial variables does not prevent the algorithm from terminating.
Once you are satisfied with this, you don't need to actually restart at all. You will get exactly the same sequence of simplex pivots if you simply include in your pivot rule the clause: never choose a nonbasic artificial variable.