There was a question about this topic a while back, see here. It is known that there is no great way to design 3-dim complex numbers: any representation will have a number of drawbacks, such as impossibility to define an inverse, non-commuting products, or a solution that is trivial and of no value. I was first interested in this topic by looking at bicomplex numbers. Though being 4-dim and similar to quaternions, they do commute. But they don't always have an inverse.
The 3-dim complex numbers that I propose here commute, and many have an inverse. I don't know if they are associative and if they accept a matrix representation like complex numbers, quaternions, or bicomplex numbers. I am interested in an answer to these questions.
I define these 3-dim complex numbers as follows. Let $z_1 = a_1 + b_1 i + c_1 j$ and $z_2 = a_2 + b_2 i + c_2 j$. Then the product is defined as
$$z_1 * z_2 = (a_1a_2-b_1b_2+c_1c_2) + (b_1a_2 +a_1b_2) i + (c_1b_2+b_1c_2) j.$$
The first issue is that it is not compatible with the scalar product:
$$(a_1 + b_1 i + c_1 j) * a_2 = a_1 a_2 + b_1 a_2 i \neq (a_1 + b_1 i + c_1 j) \cdot a_2.$$
Note that $i^2=-1, ij = j$ and $j^2=1$. When the inverse exists, it is given by the following formula:
$$(a + bi + cj)^{-1} = \frac{a}{a^2+b^2+c^2} - \frac{bi}{a^2+b^2+c^2}+\frac{cj}{a^2+b^2+c^2}.$$
I used Mathematica to find the inverse, see here.
Update
Another interesting set of 3-dim complex numbers can be found here. The multiplication is associate and commutative, but again, not all of these numbers have an inverse. In fact, fewer than in my own set.