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