We simplify using properties of determinants.
Let us denote the rows as $R_1$, $R_2$, $R_3$.
Step 1: Add all columns into $C_1$
$C_1 \rightarrow C_1 + C_2 + C_3$
New $C_1$ entries:
Row 1: $(b + c)^2 + a^2 + a^2 = (b + c)^2 + 2a^2$
Row 2: $b^2 + (c + a)^2 + b^2 = 2b^2 + (c + a)^2$
Row 3: $c^2 + c^2 + (a + b)^2 = 2c^2 + (a + b)^2$
So we can factor $(a + b + c)^2$ from the first column, and take out $abc$ as a common factor.
After applying row/column operations, determinant simplifies to:
\[
\Delta = 2abc(a + b + c)^3
\]
Full expansion is tedious but applying symmetric determinant identities and using cofactor expansion with row and column operations confirms the identity.