  Massive stars in binary systems (as WR140, WR147 or $\eta$ Carinae) have long
been regarded as potential sources of high-energy $\gamma$-rays. The emission
is thought to arise in the region where the stellar winds collide and produce
relativistic particles which subsequently might be able to emit $\gamma$-rays.
Detailed numerical hydrodynamic simulations have already offered insight in the
complex dynamics of the wind collision region (WCR), while independent
analytical studies, albeit with simplified descriptions of the WCR, have shed
light on the spectra of charged particles. In this paper, we describe a
combination of these two approaches. We present a 3D-hydrodynamical model for
colliding stellar winds and compute spectral energy distributions of
relativistic particles for the resulting structure of the WCR. The hydrodynamic
part of our model incorporates the line-driven acceleration of the winds,
gravity, orbital motion and the radiative cooling of the shocked plasma. In our
treatment of charged particles we consider diffusive shock acceleration in the
WCR and the subsequent cooling via inverse Compton losses (including
Klein-Nishina effects), bremsstrahlung, collisions and other energy loss
mechanisms.
