COMSOL Multiphysics files corresponding to the numerical models presented in the paper “Modelling magma flow within dyke-fed sill geometries: A coupled thermal and fluid dynamics approach” (https://doi.org/10.1016/j.jvolgeores.2025.108480) are provided. A detailed description of the model setup, including the selected physics, geometry, material properties, and boundary conditions, can be found in the paper. The models are distributed in COMSOL Multiphysics file format (.mph) and require COMSOL Multiphysics version 5.4 or later to be opened. For compactness of file size, computed solutions are not included; the models must be run to generate results. Only models with different geometries are included, excluding those that share the same geometry but differ in injection velocity. A table of the included models is provided below, and the definitions of the variables can be found in the referenced paper. Williams, K.M., A. Geyer, C. Annen, and J. Kavanagh. “Modelling Magma Flow within Dyke-Fed Sill Geometries: A Coupled Thermal and Fluid Dynamics Approach.” Journal of Volcanology and Geothermal Research 469 (January 2026): 108480. https://doi.org/10.1016/j.jvolgeores.2025.108480. Model Sill thickness (a) (m) Dyke(s) thickness (b) (m) Dyke1 thickness (b1) (m) Dyke2 thickness (b2) (m) Dyke spacing (s) (m) Injection velocity (uinlet) (m/s) 1FD_10_5_M 10 5 0.01 1FD_100_5_M 100 5 0.01 1FD_100_50_M 100 50 0.01 2FD_10_5_10_M 10 5 10 0.01 2FD_100_50_100_M 100 50 100 0.01 2FD_100_50_50_M 100 50 50 0.01 2FD_100_50/10_50_M 100 50 10 50 0.01 3FD_10_5_10_F 10 5 10 0.1 3FD_10_5_10_FS 10 5 10 0.1, 0.001 4FD_10_5_10_F 10 5 10 0.1 1FD_end_10_5_M 10 5 - 0.01 1FD_end100_10_5_M 10 5 - 0.01