A Discontinuous Galerkin Method for the Intrinsic Beam Model with Kelvin-Voigt Damping

We present an energy-stable discontinuous Galerkin (DG) spatial discretization of the intrinsic beam equations with Kelvin-Voigt damping. We reformulate the governing equations by splitting the total sectional resultants into elastic and viscous parts, thereby exposing a mixed system without mixed space-time derivatives. For homogeneous cantilever boundary data, the continuous system obeys an energy dissipation identity. With characteristic upwinding and alternating auxiliary traces the DG discretization is energy stable at the semidiscrete level for the stated homogeneous split boundary closure. The method is implemented in the Julia-based simulation framework Trixi.jl. A physically compatible manufactured solution with all state and viscous components active assesses convergence, while a rotating-beam test verifies the analytic constant-spin branch and the actuated energy balance.