Optimal spatial field control for controlled release

Most distributed parameter control problems involve manipulation within the spatial domain. Such problems arise in a variety of applications including epidemiology, tissue engineering, and cancer treatment. This paper proposes an approach to solve a state-constrained spatial field control problem that is motivated by a biomedical application. In particular, the considered manipulation over a spatial field is described by partial differential equations (PDEs) with spatial frequency constraints. The proposed optimization algorithm for tracking a reference spatial field combines three-dimensional Fourier series, which are truncated to satisfy the spatial frequency constraints, with exploitation of structural characteristics of the PDEs. The computational efficiency and performance of the optimization algorithm are demonstrated in a numerical example. In the example, the spatial tracking error is shown to be almost entirely due to the limitation on the spatial frequency of the manipulated field. The numerical results suggest that the proposed optimal control approach has promise for controlling the release of macromolecules in tissue engineering applications.