Abstract
Idiopathic pulmonary fibrosis (IPF) is a chronic lung disease characterized by excessive scarring and fibrosis due to the abnormal accumulation of extracellular matrix components, primarily collagen. This study aims to design and solve an optimal control problem to regulate M2 macrophage activity in IPF, thereby preventing fibrosis formation by controlling the anti-M1 agent. The research models the diffusion of M2 macrophages in inflamed tissue using a novel dynamical system with partial differential equation (PDE) constraints. The control problem is formulated to minimize fibrosis by regulating an anti-M1 agent. The study employs a two-step process of discretization followed by optimization, utilizing the Galerkin spectral method to transform the M2 diffusion PDE into an algebraic system of ordinary differential equations (ODEs). The optimal control problem is then solved using Pontryagin/s minimum principle, canonical Hamiltonian equations, and extended Riccati differential equations. The numerical simulations indicate that without control, M2 macrophage levels increase and stabilize, contributing to fibrosis. In contrast, the optimal control strategy effectively reduces M2 macrophages, preventing fibrosis formation within 120 days. The results highlight the potential of the proposed optimal control approach in modulating tissue repair processes and mitigating the progression of IPF. This study underscores the significance of targeting M2 macrophages and employing mathematical methods to develop innovative therapies for lung fibrosis.