Optimal control theory and advanced optimality conditions for a diffuse interface model of tumor growth

Matthias Ebenbeck,Patrik Knopf

DOI: https://doi.org/10.1051/cocv/2019059

2020-10-14

Abstract:In this paper, we study a distributed optimal control problem for a diffuse interface model for tumor growth. The model consists of a Cahn-Hilliard type equation for the phase field variable coupled to a reaction diffusion equation for the nutrient and a Brinkman type equation for the velocity. The system is equipped with homogeneous Neumann boundary conditions for the tumor variable, the chemical potential and the nutrient as well as a "no-friction" boundary condition for the velocity. The control acts as a medication by cytotoxic drugs and enters the phase field equation. The cost functional is of standard tracking type and is designed to track the phase field variable during the evolution and at some fixed final time. We prove that the model satisfies the basics for calculus of variations and we establish first-order and second-order conditions for local optimality. Moreover, we present a globality condition for critical controls and we show that the optimal control is unique on small time intervals.

Optimization and Control,Analysis of PDEs

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to optimize the diffuse - interface model of tumor growth through optimal control theory and advanced optimality conditions. Specifically, the authors studied a distributed optimal control problem based on a diffuse - interface model that describes tumor growth. This model consists of three main equations:

1. **Cahn - Hilliard - type equation**: It is used to describe the phase - field variable, that is, the difference in volume fraction between tumor tissue and healthy tissue.
2. **Reaction - diffusion equation**: It is used to describe the change in nutrient concentration.
3. **Brinkman - type equation**: It is used to describe the velocity field.

These partial differential equations (PDEs) are equipped with homogeneous Neumann boundary conditions (for the phase - field variable, chemical potential, and nutrients) and "friction - free" boundary conditions (for velocity). The control variable represents the administration of cytotoxic drugs, which acts in the phase - field equation. The goal is to minimize a standard tracking - type cost function that aims to track the change of the phase - field variable during time evolution and its state at a certain fixed final time point.

The main contributions of the paper include:
- Proving that the model satisfies the basic requirements of the calculus of variations.
- Establishing the conditions for first - order and second - order local optimality.
- Proposing a global condition for determining whether a critical control is globally optimal.
- Proving that the optimal control is unique within a small time interval.

Through these studies, the authors hope to find a strategy such that the administration of cytotoxic drugs can most effectively achieve the expected treatment goals while minimizing harm to patients. This not only helps to understand the dynamic process of tumor growth but also provides a theoretical basis for formulating personalized treatment plans.