On the difference between the main effect in a one-factor and a two-factor regression

Adding to @RobertLong's answer, there is a slight conceptual mistake in the way $b_2$ is described in the question in the case where $X_1$ was centered. It is indeed true that $b_2$ becomes the average effect of $X_2$ over all possible values of $X_1$, in the sense that $overline{b_2+b_3X_1}=b_2$, but it should be emphasized that this is an average of simple effects. It may have nothing to do with the main effect of $X_2$ on the DV, which means that $b_2$ may be really far from $B_2$ even without interaction.

Here is an example where there is no interaction, and $b_2$ and $B_2$ have nothing in common: the vertical axis is the DV $Y$, the horizontal axis is for the covariate $X_2$, and the colors stand for levels of the covariate $X_1$. For any value of $X_1$, the simple effect $b_2+b_3X_1$ is around $-1$, while the main effect $B_2$ is clearly positive.

$b_2$ here corresponds to the conditional effect of $X_2$ when $X_1=0$. A common mistake is to understand $b_2$ as being the main effect of $X_2$, i.e. the average effect of $X_2$ over all possible values of $X_1$.

Indeed. I typically answer at least one question per week where this mistake is made. It it also worth pointing out for completeness that $b_1$ here corresponds to the conditional effect of $X_1$ when $X_2= 0 $ and not the main effect of $X_1$ which is easily seen by rearranging the formula

$$Y=(b_0+b_2X_2)+(b_1+b_3X_2)X_1$$

In practice, it seems that $b_2$ and $B_2$ are reasonably close to each other.

I think this is false in general for this model and will will only be true when the interaction term $b_3$ is very small.

Are there any "common knowledge" examples of situations where $B_2$ and $b_2$ are remarkably far from each other?

Yes, when the $b_3$ is meaningfully large then $B_2$ and $b_2$ will be meaningfully apart. I am thinking of how to show this algebraiclly and graphically but I don't have much time now, so I will resort to a simple simulation for now. First with no interaction:

> set.seed(25)
> N <- 100
> 
> dt <- data.frame(X1 = rnorm(N, 0, 1), X2 = rnorm(N, 5, 1))
> 
> X <- model.matrix(~ X1 + X2 + X1:X2, dt)
> 
> betas <- c(10, -2, 2, 0)
> 
> dt$Y <- X %*% betas + rnorm(N, 0, 1)
> 
> (m1 <- lm(Y ~ X1*X2, data = dt))$coefficients[3]
  X2 
2.06 
> (m2 <- lm(Y ~ X2, data = dt))$coefficients[2]
  X2 
1.96

as expected. And now with an interaction:

> set.seed(25)
> N <- 100
> 
> dt <- data.frame(X1 = rnorm(N, 0, 1), X2 = rnorm(N, 5, 1))
> 
> X <- model.matrix(~ X1 + X2 + X1:X2, dt)
> 
> betas <- c(10, -2, 2, 10)
> 
> dt$Y <- X %*% betas + rnorm(N, 0, 1)
> 
> (m1 <- lm(Y ~ X1*X2, data = dt))$coefficients[3]
  X2 
2.06 
> (m2 <- lm(Y ~ X2, data = dt))$coefficients[2]
  X2 
3.29 

Are there any known upper bounds to $|b_2-B_2|$

I don't think so. As you increase $|b_3|$ then $|b_2-B_2|$ should increase