There is a mathematical formula, but it requires knowing the pressure and velocity distribution around the surface of the object:
$$D=\int_{S}pcos(\theta)dA +\int_S\tau_wsin(\theta)dA$$
To get pressure and velocity, you either need experiment data, a computational simulation, or an analytical solution to the navier stokes equations for your particular case of interest.