I am trying to determine the order of my numerical method for resolving a fluid-structure interaction problem using the immersed boundary method. I am using Crank-Nicolson to resolve the fluid equations and the immersed body is updated in time using an explicit operator splitting + projection method. The CFL condition for the fluid equations is
$$ CFL = |u_\max(t)| \frac{\Delta t}{\Delta x} \,. $$
I ran one set of simulations where I keep the timestep size fixed and refine the spatial resolution of the fluid and body meshes. This way, I'm able to determine the spatial order of accuracy of my scheme.
Now, when I try to keep my spatial resolution fixed and vary the timestep, I sometimes violate the CFL condition. Not wanting to resort to very minuscule timestep sizes, I did a convergence study where I held the CFL number fixed and varied the spatial and temporal resolutions proportionality. This tells me the overall order of my entire numerical method, however, does not distinguish between the spatial and temporal convergence rates. Is it true that in this case, if the temporal accuracy is less than the spatial accuracy, then the results will tell me the temporal convergence rate?