Lithium-ion batteries show a complex thermo-electrochemical performance and aging behavior. This paper presents a modeling and simulation framework that is able to describe both multi-scale heat and mass transport and complex electrochemical reaction mechanisms. The transport model is based on a 1D + 1D + 1D (pseudo-3D or P3D) multi-scale approach for intra-particle lithium diffusion, electrode-pair mass and charge transport, and cell-level heat transport, coupled via boundary conditions and homogenization approaches. The electrochemistry model is based on the use of the open-source chemical kinetics code CAN-TERA, allowing flexible multi-phase electrochemistry to describe both main and side reactions such as SEI formation. A model of gas-phase pressure buildup inside the cell upon aging is added. We parameterize the model to reflect the performance and aging behavior of a lithium iron phosphate (LiFePO 4 , LFP)/graphite (LiC 6 ) 26650 battery cell. Performance (0.1-10 C discharge/charge at 25, 40 and 60 • C) and calendaric aging experimental data (500 days at 30 • C and 45 • C and different SOC) from literature can be successfully reproduced. The predicted internal cell states (concentrations, potential, temperature, pressure, internal resistances) are shown and discussed. The model is able to capture the nonlinear feedback between performance, aging, and temperature. Mathematical modeling and numerical simulation have become standard techniques in lithium-ion battery research and developmentfrom the atomistic scale up to the system scale. [1][2][3][4][5] Historically, most lithium-ion cell models are based on the work of John Newman and coworkers who developed a one-dimensional model of electrochemistry and mass and charge transport in porous electrodes on the ∼100 μm scale, 6 which was later extended by transport in the active materials particles on the ∼1 μm scale, giving rise to 1D + 1D or "pseudo-2D" (P2D) models. 7,8 This type of model is widely used today. Extensions include solid electrolyte interphase (SEI) formation, 9 aging mechanisms, 10,11 impedance simulations, 12 and multi-phase chemistry in lithium-air 13,14 and lithium-sulfur 15 cells. Temperature has a strong influence on the performance and lifetime of a lithium-ion battery. A straightforward approach has been to include heat sources and heat conductivity to the standard P2D type models. [16][17][18] However, temperature gradients typically occur on a higher scale, that is, the mm and cm scale of a single cell, as compared to electrode scale described by typical P2D models. Therefore, model extensions to the cell scale are necessary. Consequently, scale-coupling methods have been developed that describe both, electrochemistry on the electrode-pair scale, and heat transport on the cell scale more efficiently. Scale coupling usually uses independent computational domains on the various scales coupled through adequate boundary conditions. As a result, models with various dimensionalities have been presented, for example, 3D + 1D + 1D (cell scale...