Type 1 diabetes (T1D) is a chronic condition affecting pancreatic insulin production, thus requiring lifelong exogenous insulin therapy to maintain blood glucose within a safe range. While continuous glucose monitoring (CGM) and other technological advancements have improved its management, clinicians still struggle to achieve good glycemic control in pediatric patients due to their evolving physiology and metabolic variability, which necessitate frequent therapy adjustments. To address this challenge, this work present a therapy optimization algorithm leveraging digital twins (DTs) to biweekly adapt the state-of-art meal insulin bolus (MIB) formula proposed by Noaro et al., ensuring its personalization based on patients’ changing metabolic dynamics estimated by DTs. The algorithm has been testes in a 24-week in silico trial involving five pediatric patients from the Tidepool Big Data Donation project, comprising 724 daily CGM traces. Biweekly optimization was performed using a derivative-free algorithm to minimize the glycemia risk index (GRI), an overall risk metric for the patient, by adjusting MIBs via personalized modulation parameters. The proposed approach demonstrated an increasing trend of the time spent in the euglycemic range (+0.63% every two weeks), and a decreasing trend of the time spent in hyperglycemia (-0.18% every two weeks) and of GRI (-0.64 every two weeks). These findings highlight the potential of the proposed approach to improve glucose control over time, thus to support clinicians in providing tailored, frequent datadriven therapy adjustments for pediatric patients, addressing the unique demands of this population.