Is there a good way to find the fan and polytope of the blowup of $\mathbb{P}^3$ along the union of two invariant intersecting lines?
Everything I find in the literature is for blowups along smooth invariant centers.
Thanks!
Is there a good way to find the fan and polytope of the blowup of $\mathbb{P}^3$ along the union of two invariant intersecting lines? Everything I find in the literature is for blowups along smooth invariant centers. Thanks! 


The question is local near the intersection of the two lines, so the more basic question is to work this out for $\mathbb{A}^3$. So we want to blow up $k[x,y,z]$ at $\langle z, xy \rangle$. There are two maximal charts: $$\mathrm{Spec} \ k[x,y,z, (xy)/z] \ \mbox{and} \ \mathrm{Spec} \ k[x,y,z,z/(xy)].$$ Note that the first chart is isomorphic to $k[a,b,c,d]/(adbc)$, so we expect it to correspond to a cone with four rays in the fan, while the second chart is isomorphic to $k[u,v,w]$, so we expect it to give an ordinary simplicial cone. Let's work out exactly how these fit together. Our first semigroup is the integer points in the cone $\mathrm{Span}_{\mathbb{R}_{+}} ((1,0,0), (0,1,0), (0,0,1), (1,1,1))$. The dual cone is $C_1 := \{ (u,v,w) : u,v,w \geq 0, \ u+v \geq w \}$. So we expect that cone to appear in the fan. Similarly, the second chart corresponds to the cone $C_2 := \{ (u,v,w) : u,v \geq 0, \ u+v \leq w \}$. We have omitted the inequality $w \geq 0$ because it is redundant. So $\mathbb{A}^3$ blown up at two crossing lines in the toric variety for the fan which consists of the two cones $C_1$ and $C_2$ above (and their faces). In order to do $\mathbb{P}^3$, you'll want to insert the picture I just described in place of the cone which corresponds to the torus invariant $\mathbb{A}^3$ around the intersection. You'll also have to modify the two other torus charts which contain one of the lines (each). Since that's a blow up at a smooth center, I'll leave it to you. 


To find the polytope associated to a toric variety directly you have to realize the variety as the closure of a map from the torus. In this case at least, it is not too hard to get such a description. Let the homogeneous coordinates of $\mathbb{P}^3$ be $[w:x:y:z]$ and the two lines be $C_1 := \lbrace w = x = 0 \rbrace$ and $C_2 := \lbrace w = y = 0 \rbrace$. Then the blowup $B$ of $\mathbb{P}^3$ along $C_1 \cup C_2$ is the closure of
PS: There is a detail to be filled: blowups along singular subvarieties are not in general normal, so a priori $B$ might not be a normal toric variety (i.e. the polytope is associated not to $B$ but the normalization of $B$). But as David shows in his answer (and probably proved for a general torus invariant subspaces in the article he mentions in the comment), $B$ is indeed normal. 


Since you already know how to blow up along either ${\mathbb P}^1$ individually, we can concentrate on what's happening nearby the intersection. Which means we can work affinely. Then the polyhedron for ${\mathbb A}^3$ is the octant $({\mathbb R}_{\geq 0})^3$, and your two lines correspond to two of the three edges. To blow them up, take a carpenter's plane to those two edges, and shave them off exactly the same amount. The result will have an edge connecting two vertices, one of degree 3, one of degree 4 (the isolated singularity on the blowup). In coordinates, the polyhedron is {$(x,y,z) : x,y,z \geq 0, x+y \geq 1, x+z \geq 1$}. If you want the blowup of ${\mathbb P}^3$ not ${\mathbb A}^3$, also impose $x+y+z \leq N$ for some large $N$ (three will do). 


If you are interested not in the toric picture, but in the structure of the variety itself, the following approach is very useful. Note that the ideal of the pair of lines has the following resolution on $P^3$: $$ 0 \to O(3) \to O(2) \oplus O(1) \to I \to 0. $$ Hence the graded algebra $\oplus I^k$ is the quotient of the symmetric algebra $\oplus S^k(O(2) \oplus O(1))$. This means that the blowup is a subscheme in $Proj(\oplus S^k(O(2) \oplus O(1))) = P_{P^3}(O(2) \oplus O(1))$. Moreover, the equation of the blowup in this projectivization of the vector bundle is given by the embedding $O(3) \to O(2) \oplus O(1)$ in the resolution above. Namely, the map $O(3) \to O(2) \oplus O(1)$ can be thought of as a global section of the line bundle $L\otimes \pi^*O(3)$, where $L$ is the Grothendieck relative $O(1)$ bundle on the projectivization and $\pi:P_{P^3}(O(2) \oplus O(1)) \to P^3$ is the projection. So, we conclude that the blowup is the zero locus of a section of a line bundle $L\otimes \pi^*O(3)$ on $P_{P^3}(O(2) \oplus O(1))$. 

