Big coupled systems are nothing new. Sounds like you may have a differential-algebraic equation (DAE)? There are solvers for this too.
Using a coupled solver, e.g. a DAE solver, lets you focus on formulating the continuous-time equations, and let the solver handle discretizing into time steps (often at high order accuracy, often with adaptive sizes).