Double aortic arch (DAA) malformation is one of the reasons for symptomatic vascular rings, the hemodynamics of which is still poorly understood. This study aims to investigate the blood flow characteristics in patient-specific double aortic arches using computational fluid dynamics (CFD). Seven cases of infantile patients with DAA were collected and their computed tomography images were used to reconstruct 3D computational models. A modified Carreau model was used to consider the non-Newtonian effect of blood and a three-element Windkessel model taking the effect of the age of patients into account was applied to reproduce physiological pressure waveforms. Numerical results show that blood flow distribution and energy loss of DAA depends on relative sizes of the two aortic arches and their angles with the ascending aorta. Ligation of either aortic arch increases the energy loss of blood in the DAA, leading to the increase in cardiac workload. Generally, the rising rate of energy loss before and after the surgery is almost linear with the area ratio between the aortic arch without ligation and the ascending aorta.