  We perform 2D and 3D kinetic simulations of reconnection-mediated turbulent
flares in a magnetized electron-positron plasma, with weak and strong radiative
cooling. Such flares can be generated around neutron stars and accreting black
holes. We focus on the magnetically-dominated regime where tension of the
background magnetic field lines exceeds the plasma rest mass by a factor
$\sigma_0 > 1$. In the simulations, turbulence is excited on a macroscopic
scale $l_0$, and we observe that it develops by forming thin, dynamic current
sheets on various scales. The deposited macroscopic energy dissipates by
energizing thermal and nonthermal particles. The particle energy distribution
is shaped by impulsive acceleration in reconnecting current sheets, gradual
stochastic acceleration, and radiative losses. We parameterize radiative
cooling by the ratio $A$ of light-crossing time $l_0/c$ to a cooling timescale,
and study the effect of increasing $A$ on the flare. When radiative losses are
sufficiently weak, $A<\sigma_0^{-1}$, the produced emission is dominated by
stochastically accelerated particles, and the radiative power depends
logarithmically on $A$. The resulting flare radiation spectrum is broad and
anisotropic. In the strong-cooling regime, $A>\sigma_0^{-1}$, stochastic
acceleration is suppressed while impulsive acceleration in the current sheets
continues to operate. As $A$ increases further, the emission becomes dominated
by thermal particles. Our simulations offer a new tool to study particle
acceleration by turbulence, especially at high energies, where cooling competes
with acceleration. We find that the particle distribution is influenced by
strong intermittency of dissipation, and stochastic acceleration cannot be
described by a universal diffusion coefficient.
