I am simulating in COMSOL a system of 3 coupled PDEs (parabolic & elliptic) along with 10 stiff ODEs. In order to have the system working, I am downsizing the time step size too much to achieve convergence. I am lucky enough that I am doing 1D simulations which don't take too much time, but this drives me to think about more complex 3D cases in which the distretized system will result in hugh DOF to be solved, with the very low step size this takes too much time.
So, is there computationally efficient method to solve such system (e.g. downsizing step size for ODEs but having larger ones for PDEs. Just imagining)