In order to integrate partial differential equations from fluid dynamics at extremely high resolution, a numerical framework has been developed to obtain high precision solutions using an adaptive grid that increases the numerical resolution locally as needed. This framework is usable on massively parallel machines and shows close-to-optimal scaling on up to \(10^5\) CPUs.