In this paper, we investigate a Local Discontinuous Galerkin (LDG) approximation for systems with balanced Orlicz-structure. We propose a new numerical flux, which yields optimal convergence rates for linear ansatz functions. In particular, our approach yields a unified treatment for problems with $(p,\delta)$-structure for arbitrary $p \in (1,\infty)$ and $\delta\ge 0$.