Exact and Positive Controllability of Boundary Control Systems

Abstract: Using the semigroup approach to abstract boundary control problems we characterize the space of all exactly reachable states. Moreover, we study the situation when the controls of the system are required to be positive. The abstract results are applied to flows in networks with static as well as dynamic boundary conditions.

\section{Introduction}
This paper is a continuation of \cite{EKNS:08, EKKNS:10} where we introduced a semigroup approach to boundary control problems and applied it to the control of flows in networks. While in these previous works we concentrated on maximal \emph{approximate} controllability, we now focus on the \emph{exact}- and \emph{positive}  controllability spaces. In particular, this will generalize and refine results given in \cite{BBEAM:13, EKNS:08, EKKNS:10} where further references to the related literature can be found.

\smallbreak
As a simple motivation, we consider as in \cite{EKNS:08} a transport process along the edges of a finite network. This system is governed through the transmission conditions in the vertices of the network which represent the ``\emph{boundary space}'' for our problem. We then like to control the behavior of this system by acting upon a single node only. In this context it is reasonable to ask the following questions.

\begin{itemize}
\item Can we reach all possible states in final time?
\item If not, can we describe the maximal possible set of reachable states?
\item Is the choice of a particular control node important?
\item Which states can be reached if only positive controls are allowed?
\end{itemize}

In Section~\ref{examples} we will address all these questions. To this end we first recall in Section~\ref{TAF} our abstract framework from \cite{EKKNS:10} as well as some basic results concerning boundary control systems. In Section~\ref{sec:EC} we then characterize boundary admissible control operators and describe the corresponding exact reachability space. In Section~\ref{Sec:pos} we turn our attention to positive boundary control systems on Banach lattices. Finally, in Section~\ref{examples} we apply our results and explicitly compute the exact (positive)  reachability spaces for three different examples of a transport equation controlled at the boundary: in $\mathbb R^m$, in a network, and in a network with dynamic boundary conditions.

\section{The Abstract Framework}\label{TAF}
We start by recalling our general setting from \cite{EKKNS:10}.

\begin{af}\label{af-bcs}
We consider
\begin{enumerate}[(i)]
\item three Banach spaces $X$, $\partial X$ and $U$, called the
\emph{state}, \emph{boundary} and \emph{control space}, resp.;
\item a closed, densely defined \emph{system operator} $A_m:D(A_m)\subseteq X\to
X$;
\item a \emph{boundary operator} $Q\in\mathcal{L}([D(\Am)],\partial X)$;
\item a \emph{control operator} $B\in\mathcal{L}(U,\partial X)$.
\end{enumerate}
\end{af}

For these operators and spaces and a \emph{control function}
$u\in\rL^1_{\mathrm{loc}}(\mathbb R_+,U)$ we then consider the \emph{abstract Cauchy problem
with boundary control}\footnote{We denote by $\dot x(t)$ the derivative of $x$ with respect to the ``time'' variable $t$.}
\begin{alignat}{2}\label{ACPBC}
\begin{cases}
\dot x(t)=A_m x(t),&t\ge0,\\
Qx(t)=Bu(t), &t\ge0,\\
x(0)=x_0.
\end{cases}
\end{alignat}
A function $x(\cdot)=x(\cdot,x_0,u)\in\rC^1(\mathbb R_+,X)$ with $x(t)\in D(A_m)$
for all $t\ge0$ satisfying \eqref{ACPBC} is called a \emph{classical
solution}. Moreover, we denote the \emph{abstract boundary control
system} associated to \eqref{ACPBC} by $\Sigma_{\textrm{BC}}(\Am,B,Q)${}.

\smallbreak In order to investigate \eqref{ACPBC} we make the following standing assumptions which in particular ensure that the \emph{un}controlled abstract Cauchy problem, i.e., \eqref{ACPBC} with $B=0$, is well-posed.

\begin{ma}
\makeatletter
\hyper@anchor{\@currentHref}%
\makeatother
\label{ma-bcs}
\begin{enumerate}[(i)]
\item The restricted operator $A\subsetA_m$ with domain $D(A):=\ker Q$ generates
a strongly continuous semigroup $(T(t))_{t\ge0}$ on $X$;
\item the boundary
operator $Q:D(A_m)\to\partial X$ is surjective.
\end{enumerate}
\end{ma}

Under these assumptions the following properties have been shown in \cite[Lem.~1.2]{Gre:87}.

\begin{lemma}\label{lem-Gre} Let Assumptions~\ref{ma-bcs} be
satisfied. Then the following assertions are true for all $\lambda,\mu\in\rho(A)$.
\begin{enumerate}[(i)]
\item $D(A_m)=D(A)\oplus\ker(\lambda-A_m)$;
\item $Q|_{\ker(\lambda-A_m)}$ is invertible and the operator
\[Q_\lambda:=(Q|_{\ker(\lambda-A_m)})^{-1}:\partial X\to\ker(\lambda-A_m)\subseteq X\]
is bounded;
\item $\RmAQ_\lambda=
\RlAQ_\mu$.
\end{enumerate}
\end{lemma}

The following operators are essential to obtain explicit
representations of the solutions of the boundary control problem \eqref{ACPBC}.

\begin{definition}
For $\lambda\in\rho(A)$ we call the operator $Q_\lambda$, introduced in
Lemma~\ref{lem-Gre}.(ii), abstract \emph{Dirichlet operator} and define
\[B_\lambda:=Q_\lambda B \in \mathcal{L}\bigl(U,\ker(\lambda-A_m)\bigr)\subset\mathcal{L}(U,X).\]
\end{definition}

By \cite[Prop.~2.7]{EKKNS:10} the solutions of \eqref{ACPBC} can be represented by the following extrapolated version of the variation of parameters formula.

\begin{proposition}\label{prop-vpf-ex} Let $x_0\in X$, $u\in\rL^1_{\mathrm{loc}}(\mathbb R_+,U)$ and $\lambda\in\rho(A)$. If
$x(\cdot)=x(\cdot,x_0,u)$ is a classical solution \eqref{ACPBC}, then it
is given by the \emph{variation of parameters formula}
\begin{equation}\label{eq-vpf-1}
x(t)=T(t)x_0+(\lambda-A_{-1})\int_0^t T(t-s)B_\lambda u(s)\,ds,\quad
t\ge0.
\end{equation}
\end{proposition}

Our aim in the sequel is to investigate which states in $X$ can be \emph{exactly} reached from $x_0=0$ by solutions of \eqref{ACPBC}. To this end we have to impose an additional assumption which, by \eqref{eq-vpf-1}, ensures that solutions for $\mathrm{L}^p$-controls have values in $X$.

\begin{definition}
Let $1\le p\le+\infty$.
Then the control operator $B\in\mathcal{L}(U,\partial X)$ is called \emph{$p$-boundary admissible} if there exist $t>0$  and $\lambda\in\rho(A)$ such that
\begin{equation}\label{b-admiss}
\int_0^{t} T(t-s)B_\lambda u(s)\,ds\in D(A)\quad\text{for all }u\in{\rL^p\bigl([0,t],U\bigr)}.
\end{equation}
\end{definition}

\begin{remark}\label{rem:Bdd-op}
From Lemma~\ref{lem-Gre}.(iii) it follows that $(\lambda-A_{-1})Q_\lambda\in\mathcal{L}(\partial X,X_{-1})$,  hence also
\[
B_A:=(\lambda-A_{-1})B_\lambda\in\mathcal{L}(U,X_{-1})
\]
is independent of $\lambda\in\rho(A)$. Then $B\in\mathcal{L}(U,\partial X)$ is $p$-boundary admissible if and only if $B_A$ is $p$-admissible in the usual sense, cf. \cite[Def.~4.1]{Wei:89a}.
This implies that if \eqref{b-admiss} is satisfied for some $t>0$ then it is satisfied for every $t>0$. Moreover, we note that $B$ is $1$-boundary admissible if $\ker(\lambda-A_m)\subset{F_1^A}$, see \cite[Lem.~A.3]{EKKNS:10}. Finally, since ${\rL^p\bigl([0,t],U\bigr)}\subset\rL^1([0,t],U)$ it follows that $1$-boundary admissibility implies $p$-boundary admissibility for all $p>1$.
\end{remark}

\smallbreak
Now assume that $B\in\mathcal{L}(U,\partial X)$ is $p$-boundary admissible. Then for fixed $\lambda\in\rho(A)$ and $t>0$ the operators $\sB_{t}^{\textrm{BC}}:{\rL^p\bigl([0,t],U\bigr)}\to X$ given by
\begin{equation}\label{bt-bc}
\sB_{t}^{\textrm{BC}} u:= (\lambda- A)\int_0^tT(t-s)B_\lambda u(s)\,ds\\
=\int_0^tT_{-1}(t-s)B_A u(s)\,ds
\end{equation}
are called the \emph{controllability maps} of the system $\Sigma_{\textrm{BC}}(\Am,B,Q)$, where the second integral initially is taken in the extrapolation space $X_{-1}$. Note that by the closed graph theorem $\sB_{t}^{\textrm{BC}}\in\mathcal{L}({\rL^p\bigl([0,t],U\bigr)},X)$. Hence,  this definition is independent of the particular choice of $\lambda\in\rho(A)$ and gives the (unique) classical solution of \eqref{ACPBC} for given $u\in \rW^{2,1}([0,t],U)$ and  $x_0=0$. This motivates the following definition.

\begin{definition}
\begin{enumerate}[(a)]
\item
The \emph{exact reachability space in time $t\ge0$} of $\Sigma_{\textrm{BC}}(\Am,B,Q)${} is defined by\footnote{By $\operatorname{rg}(T)$ we denote the range $TX\subseteq Y$ of an operator $T:X\to Y$.}
\begin{equation}\label{eq-reach-sp-t}
e\sR_t^{\textrm{BC}}:=\operatorname{rg}(\sB_{t}^{\textrm{BC}}).
\end{equation}
Moreover, we define the \emph{exact reachability space} (in arbitrary time) by
\begin{equation}\label{eq-reach-sp}
e\sR^{\textrm{BC}}:=\bigcup_{t\ge0}\operatorname{rg}(\sB_{t}^{\textrm{BC}})
\end{equation}
and call $\Sigma_{\textrm{BC}}(\Am,B,Q)${} \emph{exactly controllable} (in arbitrary time) if $e\sR^{\textrm{BC}}=X$.
\item
The \emph{approximate reachability space in time $t\ge0$} of $\Sigma_{\textrm{BC}}(\Am,B,Q)${} is defined by
\begin{equation}\label{eq-app-reach-sp-t}
a\sR_t^{\textrm{BC}}:=\overline{e\sR_t^{\textrm{BC}}}.
\end{equation}
Moreover, we define the \emph{approximate reachability space} (in arbitrary time) by
\begin{equation}\label{eq-app-reach-sp}
a\sR^{\textrm{BC}}:=\overline{\;\bigcup_{t\ge0}a\sR_t^{\textrm{BC}}\,}
\end{equation}
and call $\Sigma_{\textrm{BC}}(\Am,B,Q)${} \emph{approximately controllable} if $a\sR^{\textrm{BC}}=X$.
\end{enumerate}
\end{definition}

From \cite[Thm.~2.12 \& Cor.~2.13]{EKKNS:10} we obtain the following properties and representations of the approximate reachability space.

\begin{proposition}\label{prop:aR}
Assume that $B\in\mathcal{L}(U,\partial X)$ is $p$-boundary admissible.
Then the following holds.
\begin{enumerate}[(i)]
\item $a\sR^{\textrm{BC}}$ is a closed linear subspace, invariant under $(T(t))_{t\ge0}$ and $R(\lambda,A)$ for $\lambda>\operatorname{{\omega_0}}(A)$.
\item $a\sR^{\textrm{BC}}=\overline{\operatorname{span}}\, \bigcup_{\lambda>\omega} \operatorname{rg}(B_\lambda)$ for some $\omega>\operatorname{{\omega_0}}(A)$.
\item $a\sR^{\textrm{BC}}\subseteq\overline{\operatorname{span}}\, \bigcup_{\lambda>\operatorname{{\omega_0}}(A)}\ker(\lambda-A_m)$.
\end{enumerate}
\end{proposition}

Part~(iii) shows that there is an upper bound for the reachability space
depending on the eigenvectors of $A_m$ only, independent of the
control operator $B$. This justifies the following notion.

\begin{definition}\label{RBCmax}
The \emph{maximal reachability space} of $\Sigma_{\textrm{BC}}(\Am,B,Q)${} is defined by
\begin{equation*
\sR^{\textrm{BC}}_{\textrm{max}}:=\overline{\operatorname{span}}\,\bigcup_{\lambda>\operatorname{{\omega_0}}(A)} \ker
(\lambda-A_m).
\end{equation*}
The system $\Sigma_{\textrm{BC}}(\Am,B,Q)${} is called \emph{maximally controllable} if $e\sR^{\textrm{BC}}=\sR^{\textrm{BC}}_{\textrm{max}}$.
\end{definition}

We stress that $\sR^{\textrm{BC}}_{\textrm{max}}\ne X$ may happen (some basic examples are provided in \cite[Sec.~5]{EKNS:08}), hence the relevant question about exact or approximate
controllability is indeed to compare $e\sR^{\textrm{BC}}$ or $a\sR^{\textrm{BC}}$ to the space
$\sR^{\textrm{BC}}_{\textrm{max}}$ and not to the whole space $X$, as it is usually
done in the classical situation.

After this short summary on boundary control systems  $\Sigma_{\textrm{BC}}(\Am,B,Q)${} taken mainly from \cite{EKKNS:10} in the context of approximate controllability, we now turn our attention to the case of \emph{exact} controllability.

\section{Exact controllability}
\label{sec:EC}

\smallbreak

We start this section by giving two characterizations of $p$-boundary admissibility for a control operator $B$ which frequently simplifies the explicit computation of the associated controllability map $\sB_t^{\textrm{BC}}$.
Here for $\lambda\in\mathbb C$ we introduce the function $\varepsilon_{\la}:\mathbb R\to\mathbb C$ by $\varepsilon_{\la}(s):=e^{\lambda s}$.
Moreover, for $f\in{\rL^p[0,t]}$ and $u\in U$ we define
\[
f\otimes u\in{\rL^p\bigl([0,t],U\bigr)}
\quad\text{by}\quad
(f\otimes u)(s):=f(s)\cdot u.
\]
Finally, we denote by $\eins_{[\alpha,\beta]}$ the characteristic function of the interval $[\alpha,\beta]\subset[0,t]$.

\begin{proposition}\label{prop:range-B1}
For a control operator $B\in\mathcal{L}(U,\partial X)$ the following are equivalent.
\begin{enumerate}[(a)]
\item $B$ is $p$-boundary admissible.
\item There exist $\lambda\in\rho(A)$, $t>0$ and $M\in\mathcal{L}\bigl({\rL^p\bigl([0,t],U\bigr)},X\bigr)$ such that for all $0\le\alpha\le\beta\le t$ and $v\in U$
\begin{equation}\label{eq:strange1}
\bigl(e^{\lambda\beta}\,T(t-\beta)-e^{\lambda\alpha}\,T(t-\alpha)\bigr)B_\lambda v=
M(\varepsilon_{\la}\cdot\eins_{[\alpha,\beta]}\otimes v).
\end{equation}
\item There exist $t>0$, $\lambda_0>\operatorname{{\omega_0}}(A)$ and $M\in\mathcal{L}\bigl({\rL^p\bigl([0,t],U\bigr)},X\bigr)$ such that for all $\lambda\ge\lambda_0$ and $v\in U$
\begin{equation}\label{eq:strange2}
\bigl(e^{\lambda t}-T(t)\bigr)B_\lambda v=
M(\varepsilon_{\la}\otimes v).
\end{equation}
\end{enumerate}
Moreover, in this case the controllability map is given by $\sB_t^{\textrm{BC}}=M$.
\end{proposition}

\begin{proof}  Let $u=\varepsilon_{\la}\cdot\eins_{[\alpha,\beta]}\otimes v$ for some $\lambda\in\rho(A)$, $0\le\alpha\le\beta\le t$ and $v\in U$. Then
\begin{align}
\int_0^tT(t-s)B_\lambda u(s)\,ds\notag
&=e^{\lambda t}\int_\alpha^\beta e^{-\lambda(t-s)}\, T(t-s)B_\lambda v\,ds\\\notag
&=e^{\lambda t}\int_{t-\beta}^{t-\alpha} e^{-\lambda s}\, T(s)B_\lambda v\,ds\\\label{eq:B-zul0}
&=R(\lambda,A)\cdot\bigl(e^{\lambda\beta}\,T(t-\beta)-e^{\lambda\alpha}\,T(t-\alpha)\bigr)B_\lambda v.
\end{align}

(a)$\Rightarrow$(b). Since by assumption $B$ is $p$-boundary admissible we have $\sB_{t}^{\textrm{BC}}\in\mathcal{L}({\rL^p\bigl([0,t],U\bigr)},X)$. Hence, \eqref{bt-bc} and \eqref{eq:B-zul0} imply \eqref{eq:strange1} for $M=\sB_t^{\textrm{BC}}$.

\smallbreak
(b)$\Rightarrow$(a). We start by proving \eqref{b-admiss}. The idea is to show this first for functions of the type $u=\varepsilon_{\la}\cdot\eins_{[\alpha,\beta]}\otimes v$. Then by linearity it also holds for linear combinations of such functions and a density argument implies \eqref{b-admiss} for arbitrary $u\in{\rL^p\bigl([0,t],U\bigr)}$. To this end let $u=\varepsilon_{\la}\cdot\eins_{[\alpha,\beta]}\otimes v$ for $[\alpha,\beta]\subset[0,t]$ and $v\in U$. Then \eqref{eq:strange1} and \eqref{eq:B-zul0} imply
\begin{align}
\int_0^tT(t-s)B_\lambda u(s)\,ds\notag
&=R(\lambda,A)\cdot M(\varepsilon_{\la}\cdot\eins_{[\alpha,\beta]}\otimes v)\\
&=R(\lambda,A)\cdot Mu\label{eq:B-zul}.
\end{align}
Note that the multiplication operator $\sM_\lambda\in\mathcal{L}\bigl({\rL^p\bigl([0,t],U\bigr)}\bigr)$ defined by $\sM_\lambda u:=\varepsilon_{\la}\cdot u$ is an isomorphism (with bounded inverse $\mathcal{M}_{-\lambda}$). Hence, it maps dense sets of ${\rL^p\bigl([0,t],U\bigr)}$ into dense sets. Since the step functions are dense in ${\rL^p\bigl([0,t],U\bigr)}$ (see \cite[p.14]{ABHN:01}), the linear combinations of functions of the type $\varepsilon_{\la}\cdot\eins_{[\alpha,\beta]}\otimes v$ for $[\alpha,\beta]\subset[0,t]$ and $v\in U$ form a dense subspace of ${\rL^p\bigl([0,t],U\bigr)}$. Thus, we conclude that \eqref{eq:B-zul} holds for all $u\in{\rL^p\bigl([0,t],U\bigr)}$. Clearly this implies that $B$ is $p$-boundary admissible and $\sB_t^{\textrm{BC}}=M$.

\smallbreak
Recall that $B_A=(\lambda-A_{-1})B_\lambda$ is independent of $\lambda\in\rho(A)$. Hence, the equivalence (a)$\Leftrightarrow$(c) follows by similar arguments replacing the total set $\{\varepsilon_{\la}\cdot\eins_{[\alpha,\beta]}\otimes v:0\le\alpha<\beta\le t, v\in U\}$ by the set $\{\varepsilon_{\la}\otimes v:\lambda\ge\lambda_0, v\in U\}$ which by the Stone--Weierstra{\ss} theorem is total as well in ${\rL^p\bigl([0,t],U\bigr)}$ for all $\lambda_0>\operatorname{{\omega_0}}(A)$.
\end{proof}

We note that by linearity it would suffice that Part~(b) of Proposition~\ref{prop:range-B1} is satisfied for $\alpha=0$ and all $0\le\beta\le t$ (or for all $0\le\alpha\le t$ and $\beta=t$).

\begin{corollary}\label{cor:Bn}
Let\footnote{We use the notation $\mathbb N_l:=\{l,l+1,l+2,\ldots\}$ for the set of natural numbers starting at $l\in\mathbb N$.} $n\in\mathbb N_1$ and assume that $B$ is $p$-boundary admissible.
Then for all $u\in{\rL^p\bigl([0,nt],U\bigr)}$
\begin{equation}\label{eq:sBnt}
\sB_{nt}^{\textrm{BC}} u=\sum_{k=0}^{n-1}T(t)^k M u_{k}
\end{equation}
where $u_k\in{\rL^p\bigl([0,t],U\bigr)}$ is defined by
\begin{equation}\label{eq:def-uk}
u_{k}(s)=u\bigl((n-k-1)t+s\bigr)
\end{equation}
and $M\in\mathcal{L}\bigl({\rL^p\bigl([0,t],U\bigr)},X\bigr)$ is the operator from Proposition~\ref{prop:range-B1}.
\end{corollary}

\begin{proof} Let $u\in{\rL^p\bigl([0,nt],U\bigr)}$. Then  by \eqref{bt-bc}
\begin{align*}
\sB_{nt}^{\textrm{BC}} u
&=(\lambda-A)\int_0^{nt}T(nt-s)B_\lambda\, u(s)\,ds\\
&=(\lambda-A)\sum_{k=1}^n T\bigl((n-k)t\bigr)\int_{(k-1)t}^{kt} T(kt-s)B_\lambda\, u(s)\,ds\\
&=\sum_{k=1}^n T\bigl((n-k)t\bigr)\cdot(\lambda-A)\int_{0}^t T(t-s)B_\lambda\, u_{n-k}(s)\,ds\\
&=\sum_{k=0}^{n-1} T(t)^k\, \sB_t^{\textrm{BC}} u_{k}.
\qedhere
\end{align*}
\end{proof}

In Section \ref{examples} we will see that \eqref{eq:strange1}, \eqref{eq:strange2}, and  \eqref{eq:sBnt} allow us to easily compute the controllability map in the situations studied in \cite[Sect.~4]{EKNS:08} and \cite[Sect.~3]{EKKNS:10} dealing with the control of flows in networks.

\begin{corollary}\label{cor:Reach}
If $B$ is $p$-boundary admissible,
 then the exact reachability space in time $nt$ for $n\in\mathbb N_1$ is given by
\begin{equation*}
e\sR_{nt}^{\textrm{BC}}=\l\{\sum_{k=0}^{n-1}T(t)^k M u_k : u_k\in{\rL^p\bigl([0,t],U\bigr)},\; 1\le k\le n-1 \r\},
\end{equation*}
where $M\in\mathcal{L}\bigl({\rL^p\bigl([0,t],U\bigr)},X\bigr)$ is the operator from Proposition~\ref{prop:range-B1}.
\end{corollary}

\section{Positive controllability}\label{Sec:pos}

In this section we are interested in positive control functions yielding positive states. To this end we will make the following

\begin{ada}
The spaces
$X$
and $U$ are Banach lattices.
\end{ada}
Moreover, by $Y^+:=\{y\in Y:Y\ge0\}$ we denote the positive cone in a Banach lattice $Y$.

\smallbreak
Note that in the sequel we do \emph{not} make any positivity assumptions on $(T(t))_{t\ge0}$, $B$ or $Q_\lambda$ if not stated otherwise.

\begin{definition}
\begin{enumerate}[(a)]
\item The \emph{exact positive reachability space in time $t\ge0$} of system $\Sigma_{\textrm{BC}}(\Am,B,Q)${} is defined by
\begin{equation}\label{eq-reach-sp-t+}
e^+\sR_t^{\textrm{BC}}:=\Bigl\{\sB_{t}^{\textrm{BC}} u:u\in{\rL^p\bigl([0,t],U^+\bigr)}\Bigr\}.
\end{equation}
Moreover, we define the \emph{exact positive reachability space} (in arbitrary time) by
\begin{equation}\label{eq-reach-sp+}
e^+\sR^{\textrm{BC}}:=\bigcup_{t\ge0}e^+\sR_t^{\textrm{BC}}
\end{equation}
and call $\Sigma_{\textrm{BC}}(\Am,B,Q)${}  \emph{exactly positive controllable} (in arbitrary time) if \\$e^+\sR^{\textrm{BC}}=X^+$.
\smallbreak
\item The \emph{approximate positive reachability space in time $t\ge0$}  of $\Sigma_{\textrm{BC}}(\Am,B,Q)${} is defined by
\begin{equation}\label{eq-a-reach-sp+}
a^+\sR_t^{\textrm{BC}}:=\overline{e^+\sR_t^{\textrm{BC}}}.
\end{equation}
Moreover, we define the \emph{approximate positive reachability space} (in arbitrary time) by
\begin{equation}\label{eq-app-reach-sp+}
a^+\sR^{\textrm{BC}}:=\overline{\;\bigcup_{t\ge0}a^+\sR_t^{\textrm{BC}}\,}
\end{equation}
and call $\Sigma_{\textrm{BC}}(\Am,B,Q)${}  \emph{approximately positive controllable} if $a^+\sR^{\textrm{BC}}=X^+$.
\end{enumerate}
\end{definition}

First we give necessary and sufficient conditions implying that starting from the initial state $x_0=0$ positive controls result in positive states.

\begin{proposition}\label{prop:reach-pos}
 Assume that $B\in\mathcal{L}(U,\partial X)$ is $p$-boundary admissible. Then
\begin{equation}\label{eq:pos-R-1}
e^+\sR_t^{\textrm{BC}}\subset X^+
\end{equation}
if and only if
\begin{equation}\label{eq:pos-R-1.5}
a^+\sR_t^{\textrm{BC}}\subset X^+
\end{equation}
if and only if there exists $\lambda\in\mathbb R\cap\rho(A)$ such that
\begin{equation}\label{eq:pos-R-2}
\bigl(e^{\lambda\beta}\,T(t-\beta)-e^{\lambda\alpha}\,T(t-\alpha)\bigr)B_\lambda\ge0
\qquad\text{for all $0\le\alpha\le\beta\le t$.}
\end{equation}
Moreover, if $(T(t))_{t\ge0}$ is positive, then the above assertions are satisfied if and only if
\begin{equation}\label{eq:pos-R-1x}
e^+\sR^{\textrm{BC}}\subset X^+
\end{equation}
if and only if
\begin{equation}\label{eq:pos-R-1.5x}
a^+\sR^{\textrm{BC}}\subset X^+
\end{equation}
if and only if there exists $\lambda>\operatorname{{\omega_0}}(A)$ and $t>0$ such that
\begin{equation}\label{eq:pos-R-3}
\bigl(e^{\lambda s}-T(s)\bigr)B_\lambda\ge0
\qquad\text{for all $0\le s\le t$}
\end{equation}
if and only if there exists  $\lambda_0>\operatorname{{\omega_0}}(A)$ such that
\begin{equation}\label{eq:pos-R-4}
B_\lambda\ge0
\qquad\text{for all $\lambda\ge\lambda_0$}.
\end{equation}
\end{proposition}

\begin{proof} The equivalence of \eqref{eq:pos-R-1} and  \eqref{eq:pos-R-1.5} follows from the closedness of $X^+$. To show the equivalence of \eqref{eq:pos-R-1} and \eqref{eq:pos-R-2} recall that by \cite[p.14]{ABHN:01} the step functions are dense in ${\rL^p\bigl([0,t],U\bigr)}$. Since the map $u\mapsto u^+$ on ${\rL^p\bigl([0,t],U\bigr)}$ is continuous, we conclude that  the positive step functions are dense in ${\rL^p\bigl([0,t],U^+\bigr)}$.
The claim then follows from (the proof of) Proposition~\ref{prop:range-B1} using the boundedness of the controllability map $\sB_t^{\textrm{BC}}$.

Now assume that $(T(t))_{t\ge0}$ is positive. Then the equivalences of \eqref{eq:pos-R-1}, \eqref{eq:pos-R-1.5}  with \eqref{eq:pos-R-1x}, \eqref{eq:pos-R-1.5x} follow from Corollary~\ref{cor:Reach} using the fact that the reachability spaces are growing in time. In particular, this implies that if \eqref{eq:pos-R-2} holds for some $t>0$ it holds for arbitrary $t>0$ and choosing $\beta=t$ and $\alpha=0$ we obtain \eqref{eq:pos-R-3} for arbitrary $t>0$.

To show the remaining assertions we fix some $\lambda>\operatorname{{\omega_0}}(A)$ and define on $\mathcal{X}:=X\times X$ the operator matrix
\[
\mathcal{A}:=
\begin{pmatrix}
A-\lambda&0\\0&0
\end{pmatrix},\quad
D(\mathcal{A}):=\Bigl\{\tbinom{x}{y}\in D(A_m)\times\partial X:Qx=By\Bigr\}.
\]
Then by \cite[Cor.~3.4]{Eng:99} the matrix $\mathcal{A}$ generates a $C_0$-semigroup $(\mathcal{T}(t))_{t\ge0}$ given by
\begin{equation}\label{eq:sTt}
\mathcal{T}(s)=\begin{pmatrix}
e^{-\lambda s}T(s)&\bigl(I-e^{-\lambda s}T(s)\bigr)B_{\lambda}\\0&I
\end{pmatrix},
\quad s\ge0.
\end{equation}
Moreover, by \cite[Lem.~3.1]{Eng:99} we have $(0,+\infty)\subset\rho(\mathcal{A})$ and
\begin{equation} \label{eq:RsA}
R(\mu,\mathcal{A})=
\begin{pmatrix}
R(\mu+\lambda,A)&\frac1\mu B_{\mu+\lambda}\\
0&\frac1\mu
\end{pmatrix}
\quad\text{for $\mu>0$}.
\end{equation}
Now, if \eqref{eq:pos-R-3} holds then $\mathcal{T}(s)\ge0$ for all $0\le s\le t$ which implies that $(\mathcal{T}(t))_{t\ge0}$ is positive which is equivalent to the fact that $\mathcal{A}$ is resolvent positive. However, by \eqref{eq:RsA} the latter is the case if and only if \eqref{eq:pos-R-4} is satisfied which shows the equivalence of \eqref{eq:pos-R-3} and \eqref{eq:pos-R-4}.
Finally, if \eqref{eq:pos-R-3} holds, then
\begin{align*}
\bigl(e^{\lambda\beta}\,T(t-\beta)-e^{\lambda\alpha}\,T(t-\alpha)\bigr)B_\lambda
&=e^{\lambda\beta}\,T(t-\beta)\cdot\bigl(I-e^{-\lambda(\beta-\alpha)}\,T(\beta-\alpha)\bigr)B_\lambda\\
&\ge0
\end{align*}
for all $0\le\alpha\le\beta\le t$. This proves \eqref{eq:pos-R-2} and completes the proof.
\end{proof}

In the sequel we use the notation $\operatorname{co} M$ and $\overline{\operatorname{co}}\, M$ to indicate the  convex hull and the closed convex hull of a set $M\subset X$, respectively.

\begin{proposition}\label{prop:aRp-char}
Assume that $B\in\mathcal{L}(U,\partial X)$ is $p$-boundary admissible and that $e^+\sR_t^{\textrm{BC}}\subset X^+$.
Then the following holds.
\begin{enumerate}[(i)]
\item $a^+\sR^{\textrm{BC}}$ is a closed convex cone, invariant under $(T(t))_{t\ge0}$ and $R(\lambda,A)$ for $\lambda>\operatorname{{\omega_0}}(A)$.
\item
$
a^+\sR^{\textrm{BC}}=\overline{\operatorname{co}}\,\l\{\bigl(e^{\lambda \beta}T(t-\beta)-e^{\lambda \alpha}T(t-\alpha)\bigr)B_\lambda v:0\le \alpha\le\beta\le t,\,v\in U^+\r\}
$
for all $\lambda>\operatorname{{\omega_0}}(A)$.
\item $a^+\sR^{\textrm{BC}}=\overline{\operatorname{co}}\,\{T(t)B_\lambda v:t\ge0,\,\lambda>w,\,v\in U^+\}$ for all $w>\operatorname{{\omega_0}}(A)$.
\item $a^+\sR^{\textrm{BC}}=\overline{\operatorname{co}}\,\{R(\lambda,A)^nB_\lambda v:n\in\mathbb N_0,\,\lambda>w,\,v\in U^+\}$ for some/all $w>\operatorname{{\omega_0}}(A)$.
\end{enumerate}
\end{proposition}

\begin{proof} (i). Clearly, $a^+\sR^{\textrm{BC}}$ is a closed convex cone.  Its invariance  under $(T(t))_{t\ge0}$ and $R(\lambda,A)$ for $\lambda>\operatorname{{\omega_0}}(A)$ follows from the representations in (iii) and (iv).

To show (ii) we note that by \eqref{bt-bc} and \eqref{eq:B-zul0} the inclusion ``$\supseteq$'' holds.
Now recall that  the positive step functions are dense in ${\rL^p\bigl([0,t],U^+\bigr)}$ and invariant under positive convey combinations. Hence, the boundedness of the controllability maps implies equality of the spaces in (ii).

\smallbreak

To obtain (iii) we note that by \eqref{bt-bc} and \eqref{eq:B-zul0} we have
\[
\bigl(e^{\lambda\beta}\,T(t-\beta)-e^{\lambda\alpha}\,T(t-\alpha)\bigr)B_\lambda v\ine^+\sR^{\textrm{BC}}
\]
for all $0\le\alpha\le\beta\le t$ and $v\in U^+$.
Multiplying this inclusion by $e^{-\lambda\beta}>0$ and putting $s:=t-\beta$ and $r:=t-\alpha$ implies
\[
\bigl(T(s)-e^{\lambda(s-r)}\,T(r)\bigr)B_\lambda v\ine^+\sR^{\textrm{BC}}
\]
for all $0\le s\le r$ and $v\in U^+$. Since $\lambda>\operatorname{{\omega_0}}(A)$ we obtain
\[
\lim_{r\to+\infty}e^{\lambda(s-r)}\,\|T(r)\|=0
\]
and hence
\[T(s)B_\lambda v\ina^+\sR^{\textrm{BC}}\]
for all $s\ge0$ and $v\in U^+$. This shows the inclusion ``$\supseteq$'' in (iii).
For the converse inclusion in (iii) it suffices to prove that
\begin{equation}\label{eq:inc2.(ii)}
e^+\sR_t^{\textrm{BC}}\subset\overline{\operatorname{co}}\,\bigl\{T(s)B_\mu y:s\ge0,\,\mu>w,\,y\in U^+\bigr\}
\end{equation}
for all $t>0$ and $w>\operatorname{{\omega_0}}(A)$. Since $B\in\mathcal{L}(U,\partial X)$ is $p$-boundary admissible the controllability map $\sB_t^{\textrm{BC}}$ is continuous. Moreover, the positive step functions are dense in ${\rL^p\bigl([0,t],U^+\bigr)}$ and $\overline{\operatorname{co}}\,\{T(s)B_\mu y:s\ge0,\,\mu>w,\,y\in U^+\}$ is a convex cone. Combining these facts and \eqref{eq:B-zul0}  it follows that \eqref{eq:inc2.(ii)} holds if
\begin{equation}\label{eq::inc2.(ii)-1}
\bigl(e^{\lambda\beta} T(t-\beta)-e^{\lambda\alpha} T(t-\alpha)\bigr)B_\lambda v\in \overline{\operatorname{co}}\,\bigl\{T(s)B_\mu y:s\ge0,\,\mu>w,\,y\in U^+\bigr\}
\end{equation}
for all  $0\le\alpha\le\beta\le t$, $k\in\mathbb N_0$ and $v\in X^+$.
Since $(T(t))_{t\ge0}$ is strongly continuous the following integral is the limit of Riemann sums, hence for $\nu>\max\{0,w\}$ we obtain using Lemma~\ref{lem-Gre}.(iii)
\begin{align*}
\overline{\operatorname{co}}\,\bigl\{T(s)B_\mu y:s\ge0,\,\mu>w,\,y\in U^+\bigr\}&\ni\nu\int_{t-\beta}^{t-\alpha}e^{\lambda(t-r)}T(r)B_\nu v\,dr\\
&=\bigl(e^{\lambda\beta}\,T(t-\beta)-e^{\lambda\alpha}\,T(t-\alpha)\bigr)\nu\RlAB_\nu v\\
&=\bigl(e^{\lambda\beta}\,T(t-\beta)-e^{\lambda\alpha}\,T(t-\alpha)\bigr)\nu\RnAB_\lambda v\\
&\to\bigl(e^{\lambda\beta}\,T(t-\beta)-e^{\lambda\alpha}\,T(t-\alpha)\bigr)B_\lambda v,
\end{align*}
as $\nu\to+\infty$.
This proves \eqref{eq::inc2.(ii)-1} and completes the proof of (iii).
\smallbreak
That the right-hand-sides of the equalities in (iii) and (iv) coincide follows from the integral representation of the resolvent (see \cite[Cor.~II.1.11]{EN:00}) and the Post--Widder inversion formula (see \cite[Cor.~III.5.5]{EN:00}). For the details we refer to the proof of \cite[Prop.~3.3]{BBEAM:13}.
\end{proof}

\begin{corollary}
Assume that $B\in\mathcal{L}(U,\partial X)$ is $p$-boundary admissible and that $a^+\sR^{\textrm{BC}}\subset X^+$. Then the following are equivalent.
\begin{enumerate}[(a)]
\item The system $\Sigma_{\textrm{BC}}(\Am,B,Q)${} is approximately positive controllable.
\item There exists $w>\operatorname{{\omega_0}}(A)$ such that  the following implication holds for all $\phi\in X'$
\[
\big<T(s)B_\lambda v,\phi\big>\ge0\text{ for all $v\in U^+$, $s\ge0$ and $\lambda>w$ $\Rightarrow$ $\phi\ge0$.}
\]
\item There exists $w>\operatorname{{\omega_0}}(A)$ such that the following implication holds for all $\phi\in X'$
\[
\bigl<R(\lambda,A)^nB_\lambda v,\phi\bigr>\ge0 \text{ for all $v\in U^+$, $n\in\mathbb N$ and $\lambda>w$ $\Rightarrow$ $\phi\ge0$.}
\]
\end{enumerate}
\end{corollary}

\begin{proof}
This follows from the proof of \cite[Thm.~3.4]{BBEAM:13} by replacing \cite[Prop.~3.3]{BBEAM:13}  with our Proposition~\ref{prop:aRp-char}.
\end{proof}

\begin{remark}
The previous two results generalize \cite[Prop.~3.3 and Thm.~3.4]{BBEAM:13}, respectively, where it is assumed that $(T(t))_{t\ge0}$, $B$ and $Q_{\lambda}$ for all $\lambda>\lambda_0$ are all positive and, in particular, the additional hypothesis
\begin{itemize}
\item[(H)]  \emph{There exists $\gamma>0$ and $\lambda_0\in\mathbb R$ such that $\|Qx\|\ge\gamma\lambda\|x\|$ for all $\lambda>\lambda_0$ and $x\in\ker(\lambda-A_m)$}
\end{itemize}
is made. We note that Hypothesis~(H) in reflexive state spaces $X$ implies that $A=A_m$, cf. \cite[Lem.~A.1]{ABE:16}. Hence, the results of \cite{BBEAM:13} are, e.g., not applicable to  state space like $X=\mathrm{L}^p([a,b],Y)$ for $p\in(1,+\infty)$ and reflexive $Y$.
\end{remark}

\smallskip
Combining Corollary~\ref{cor:Bn} and Proposition \ref{prop:reach-pos}
we finally obtain the following characterization of an exact positive reachability space.

\begin{corollary}\label{cor:Reach_pos}
Assume that $B$ is $p$-boundary admissible, $t>0$ and $n\in\mathbb N_1$.
Then the exact positive reachability space in time $nt$ is given by
\begin{equation*}
e^+\sR_{nt}^{\textrm{BC}}=\l\{\sum_{k=0}^{n-1}T(t)^k M u_k : u_k\in{\rL^p\bigl([0,t],U^+\bigr)},\; 1\le k\le n-1 \r\},
\end{equation*}
where $M\in\mathcal{L}\bigl({\rL^p\bigl([0,t],U^+\bigr)},X\bigr)$ is the operator from Proposition~\ref{prop:range-B1}. Moreover, the operator $M$ is positive  if and only if $a^+\sR_t^{\textrm{BC}}\subset X^+$.
\end{corollary}

\section{Examples}\label{examples}

In this section we will show how our abstract results can be applied to a transport equation with boundary control and to the vertex control of flows in networks.

\subsection{Exact Boundary Controllability for a Transport Equation}

In this subsection we study a transport equation in $\mathbb R^m$ given by\footnote{We denote by $x'(t,s)$ the derivative of $x(t,s)$ with respect to the ``space'' variable $s$.}
\begin{alignat}{2}\label{eq:TE}
\begin{cases}
\dot x(t,s)=x'(t,s),& s\in[0,1],\ t\ge0,\\
x(t,1)=\mathbb B x(t,0)+u(t)\cdot b,&t\ge0,\\
x(0,s)=0,& s\in[0,1].
\end{cases}
\end{alignat}
Here $x:\mathbb R_+\times[0,1]\to\mathbb C^m$, $\mathbb B\in\rM_{m}(\mathbb C)$, $u:\mathbb R_+\to\mathbb C$ and $b\in\mathbb C^m$. In order to fit this system in our general framework we choose

\begin{itemize}
\item the  state space $X:={\rL^p\bigl([0,1],\Cm\bigr)}$ where $1\le p<+\infty$,
\item the boundary space $\partial X:=\mathbb C^m$,
\item the control space $U:=\mathbb C$,
\item the control operator $B:=b\in\mathbb C^m\simeq\mathcal{L}(U,\partial X)=\mathcal{L}(\mathbb C,\mathbb C^m)$,
\item the system operator
\[A_m:=\operatorname{diag}\l(\frac{d}{ds}\r)_{m\times m}\quad\text{with domain}\quad D(A_m):=\rW^{1,p}\bigl([0,1],\Cm\bigr),\]
\item the boundary operator $Q:\rW^{1,p}\bigl([0,1],\Cm\bigr)\to\mathbb C^m$, $Qf:=f(1)-\mathbb B f(0)$,
\item the operator $A\subsetA_m$ with domain $D(A)=\ker Q$,
\item the state trajectory $x:\mathbb R_+\to{\rL^p\bigl([0,1],\Cm\bigr)}$, $x(t):=x(t,\cdot)$.
\end{itemize}
With these choices the controlled transport equation~\eqref{eq:TE} can be reformulated as an abstract Cauchy problem with boundary control of the form \eqref{ACPBC}. Clearly, the boundary operator $Q$ is surjective.

\smallbreak
By \cite[Cor.~18.4]{BKR:17}
we know that for $\lambda\in\mathbb C$ and $A=A_m|_{\ker(Q)}$ as above we have
\[
\lambda\in\rho(A)\iff
e^\lambda\in\rho(\mathbb B).
\]
Moreover, by  \cite[Prop.~18.7]{BKR:17}
 the operator $A$ generates a strongly continuous semigroup given by
\begin{equation}\label{eq:HG-TE}
\bigl(T(t)f\bigr)(s)=\mathbb B^k f(t+s-k)\quad \text{if }\;t+s\in[k,k+1) \text{ for } k\in\mathbb N_0,
\end{equation}
where $\mathbb B^0:=Id$.
This shows that the Assumptions~\ref{ma-bcs} are satisfied. To proceed we have to compute the associated Dirichlet operator.

\begin{lemma} For $\lambda\in\rho(A)$ the Dirichlet operator $Q_\lambda\in\mathcal{L}\bigl(\mathbb C^m,{\rL^p\bigl([0,1],\Cm\bigr)}\bigr)$ is given by
\begin{equation}\label{eq:Dirich-TE}
Q_\lambda=\varepsilon_{\la}\otimes R(e^{\lambda},\mathbb B).
\end{equation}
\end{lemma}

\begin{proof}
By Lemma~\ref{lem-Gre}.(ii) we know that $Q:{\ker(\lambda-A_m)}\to\partial X$ is invertible. Moreover, for $d\in\mathbb C^m=\partial X$ we have
\begin{align*}
Q\bigl(\varepsilon_{\la}\otimes R(e^{\lambda},\mathbb B)\,d\bigr)
=e^{\lambda}\cdot R(e^{\lambda},\mathbb B)\,d-\mathbb B\cdot R(e^{\lambda},\mathbb B)\,d=d
\end{align*}
which proves \eqref{eq:Dirich-TE}.
\end{proof}

Next we verify that in this context \eqref{eq:strange1} holds.

\begin{lemma}\label{lem:strange-TE}
Let $\lambda\in\rho(A)$. Then for all $0\le\alpha\le1$
\begin{equation}\label{eq:strange-TE}
\bigl(e^{\lambda\alpha}\cdot T(1-\alpha)B_\lambda \bigr)(s)=
\begin{cases}
\varepsilon_{\la}(1+s)\cdotR(e^{\lambda},\mathbb B)\, b&\text{if }0\le s<\alpha,\\
\varepsilon_{\la}(1+s)\cdotR(e^{\lambda},\mathbb B)\, b-\varepsilon_{\la}(s)\cdot b&\text{if }\alpha\le s\le1.
\end{cases}
\end{equation}
Hence, \eqref{eq:strange1} is satisfied for
\begin{align*}
M&=b\in\mathcal{L}\bigl({\rL^p[0,1]},{\rL^p\bigl([0,1],\Cm\bigr)}\bigr),\
(M u)(s)=u(s)\cdot b.
\end{align*}
\end{lemma}
\begin{proof} The claim follows from \eqref{eq:Dirich-TE} and \eqref{eq:HG-TE} by the following simple computation.
\begin{align*}
\bigl(e^{\lambda\alpha}\cdot T(1-\alpha)B_\lambda \bigr)(s)
&=e^{\lambda\alpha}\cdot\Bigl(T(1-\alpha)\bigl(\varepsilon_{\la}\otimesR(e^{\lambda},\mathbb B)\,b\bigr)\Bigr)(s)
\\
&=e^{\lambda\alpha}\cdot
\begin{cases}
\varepsilon_{\la}(1-\alpha+s)\cdotR(e^{\lambda},\mathbb B)\, b&\text{if }0\le s<\alpha,\\
\varepsilon_{\la}(s-\alpha)\cdot\BBR(e^{\lambda},\mathbb B)\, b&\text{if }\alpha\le s\le1,
\end{cases}
\\
&=
\begin{cases}
\varepsilon_{\la}(1+s)\cdotR(e^{\lambda},\mathbb B)\, b&\text{if }0\le s<\alpha,\\
\varepsilon_{\la}(1+s)\cdotR(e^{\lambda},\mathbb B)\, b-\varepsilon_{\la}(s)\cdot b&\text{if }\alpha\le s\le1.
\end{cases}
\qedhere
\end{align*}
\end{proof}
Thus $B$ is $p$-boundary admissible. Next we compute the appropriate reachability space.

\begin{corollary}\label{cor:Reach-TE}
If $t\ge m$ then the exact reachability space of the controlled transport equation~\eqref{eq:TE} is given by
\begin{equation*
e\sR_t^{\textrm{BC}}=e\sR^{\textrm{BC}}={\rL^p[0,1]}\otimes\operatorname{span}\l\{b,\mathbb B b,\ldots,\mathbb B^{m-1}b\r\}.
\end{equation*}
\end{corollary}

\begin{proof}
Note that by \eqref{eq:HG-TE} we have $T(1)f=\mathbb B f$. Hence, for $t=m$ the assertion follows immediately from Corollary~\ref{cor:Reach} and Lemma~\ref{lem:strange-TE}. Clearly, $e\sR_t^{\textrm{BC}}$ increases in time $t\ge0$. However, by the Cayley--Hamilton theorem $\operatorname{span}\{b,\mathbb B b,\ldots,\mathbb B^{l}b\}=\operatorname{span}\{b,\mathbb B b,\ldots,\mathbb B^{m-1}b\}$ for all $l\ge m-1$ and the claim follows.
\end{proof}

\begin{remark}\label{rem:C-H}
Let $l\le m$ be the degree of the minimal polynomial of $\mathbb B$. Then the previous proof shows that for all $t\ge l$ we even have
\begin{equation*
e\sR_t^{\textrm{BC}}=e\sR^{\textrm{BC}}={\rL^p[0,1]}\otimes\operatorname{span}\l\{b,\mathbb B b,\ldots,\mathbb B^{l-1}b\r\}.
\end{equation*}
\end{remark}

\begin{corollary}\label{cor:TE}
The following assertions are equivalent.
\begin{enumerate}[(a)]
\item Equation~\eqref{eq:TE} is exactly boundary controllable in time $t\ge m$, i.e.,
$e\sR_t^{\textrm{BC}}=X$.
\item Equation~\eqref{eq:TE} is maximally controllable in time $t\ge m$, i.e.,
$e\sR_t^{\textrm{BC}}=\sR^{\textrm{BC}}_{\textrm{max}}\,$.
\item $\operatorname{span}\l\{b,\mathbb B b,\ldots,\mathbb B^{m-1}b\r\}=\mathbb C^m$.
\end{enumerate}
\end{corollary}

\begin{proof} Note that $\ker(\lambda-A_m)=\varepsilon_{\la}\otimes\mathbb C^m$. Since by the Stone--Weierstra\ss{} theorem we have
\[\overline{\operatorname{span}}\,\bigcup_{\lambda>\operatorname{{\omega_0}}(A)} \{\varepsilon_{\la}\}={\rL^p[0,1]},\]
the maximal reachability space equals
\[
\sR^{\textrm{BC}}_{\textrm{max}}={\rL^p[0,1]}\otimes\mathbb C^m = X
\]
and the assertions follow immediately from Corollary~\ref{cor:Reach-TE}.
\end{proof}

\begin{remark}
The previous result characterizes the \emph{exact} maximal boundary controllability by a one-dimensional control in terms of a Kalman-type condition  which is well-known in control theory.
\end{remark}

Combining Remark~\ref{rem:C-H} and Corollary~\ref{cor:TE} we furthermore obtain the following

\begin{corollary} Let $l\in\mathbb N$ be the degree of the minimal polynomial of $\mathbb B$. If $l < m$, the transport equation~\eqref{eq:TE} is not
maximally controllable, i.e., $e\sR^{\textrm{BC}} \subsetneq\sR^{\textrm{BC}}_{\textrm{max}}$.
\end{corollary}

\smallbreak

Finally, we investigate positive controllability and consider
\begin{itemize}
\item the positive cone  $X^+:={\rL^p\bigl([0,1],\mathbb R_+^m\bigr)}$ in the state space $X$,
\item the positive cone  $U^+:=\mathbb R_+$ in the control space $U$,
\item a positive matrix $\mathbb B\in\mathrm{M}_m(\mathbb R_+)$,
\item a positive control operator $B:=b\in\mathbb R_+^m$.
\end{itemize}
Then by \eqref{eq:HG-TE}--\eqref{eq:Dirich-TE} the operators $T(t)\in\mathcal{L}(X)$ for $t\ge0$ and $B_\lambda\in\mathcal{L}(U,X)$ for $\lambda>\operatorname{{\omega_0}}(A)$ are positive.
Thus arguing as above using Proposition \ref{prop:reach-pos} and Corollary~\ref{cor:Reach_pos} we obtain the following.

\begin{corollary}\label{cor:pos-Reach-TE}
The exact positive reachability space of the controlled transport equation~\eqref{eq:TE} is given by
\begin{equation*}
e^+\sR^{\textrm{BC}}={\rL^p\bigl([0,1],\mathbb R_+\bigr)}\otimes\operatorname{co}\bigl\{\mathbb B^{k}b \;:\; k\in\mathbb N_0\bigr\}.
\end{equation*}
Hence, the problem is exactly positive controllable if and only if
\begin{equation*}
\operatorname{co}\bigl\{\mathbb B^{k}b \;:\; k\in\mathbb N_0\bigr\} = \mathbb R_+^m.
\end{equation*}
\end{corollary}

\subsection{Vertex control of flows in networks}

The previous example can be easily adapted to cover a transport problem on a network controlled in a single vertex. More precisely, consider a network consisting of $n$ vertices $\{v_1,\dots ,v_n\}$ and $m$ edges $\{e_1,\dots ,e_m\}$. As shown in \cite[Sec.~18.1]{BKR:17}, its structure can be described by
either the transposed weighted adjacency matrix
$\mathbb A\in\rM_{n}(\mathbb C)$ given by
\begin{equation*}
 \mathbb A_{ij}:=\begin{cases}
w _{jk} & \mbox{if $v_j\stackrel{e_k}{\longrightarrow}v_i$},  \\
0 & \text{otherwise,}
\end{cases}
\end{equation*}
or by the transposed weighted adjacency matrix of the line graph
$\mathbb B\in\rM_{m}(\mathbb C)$ where
\begin{equation*}
\mathbb B_{ij}:=\begin{cases}
w_{ki} & \mbox{if $\stackrel{e_j}{\longrightarrow}v_k \stackrel{e_i}{\longrightarrow}$}\ ,  \\
0 & \text{otherwise.}
\end{cases}
\end{equation*}
To proceed we also need the transposed weighted outgoing incidence matrix
$(\Fiwm)^\top=:\Psi\in\rM_{m\times n}(\mathbb C)$ defined by
\begin{equation*}
    \Psi_{ij} :=
    \begin{cases}
   w_{ij} & \mbox{if $v_j\stackrel{e_i}{\longrightarrow}$}\ , \\
    0 & \text{otherwise}
    \end{cases}
\end{equation*}
and the corresponding unweighted outgoing
incidence matrix denoted by $\Phi^{-}\in\rM_{n\times m}(\mathbb C)$.
For the weights we assume $0\le w_{ij} \le 1$, thus all these matrices are positive. Moreover, we assume that $\Psi$ is column stochastic (i.e., the weights on all the outgoing edges from a given vertex sum up to $1$).
For a detailed account of the various graph matrices we refer to \cite[Sec.~18.1]{BKR:17}. Here we only mention the following relations
\begin{equation}\label{eq:ABrel}
\Psi\mathbb A = \mathbb B\Psi,\quad
\Psi R(\lambda,\mathbb A) = R(\lambda,\mathbb B)\Psi,\quad
\text{and}\quad
\Phi^{-}\Psi = Id_{\mathbb C^n}
\end{equation}
which we will need in the sequel.

\smallbreak
We then consider a transport equation on the $m$ edges imposing $n$ boundary conditions in the vertices, controlled in a single vertex $v_i$, i.e.,
\begin{alignat}{2}\label{eq:TE-net}
\begin{cases}
\dot x(t,s)=x'(t,s),& s\in[0,1],\ t\ge0,\\
x(t,1)=\mathbb B x(t,0)+u(t)\cdot \Psi v,&t\ge0,\\
x(0,s)=0,& s\in[0,1],
\end{cases}
\end{alignat}
where $x:\mathbb R_+\times[0,1]\to\mathbb C^m$, $u:\mathbb R_+\to\mathbb C$,  and  the vertex $v=v_i$ is represented by the $i$-th canonical basis vector in $\mathbb C^n$. To rewrite this equation in an abstract form we take the same state space $X:={\rL^p\bigl([0,1],\Cm\bigr)}$, control space $U:=\mathbb C$  and boundary space $\partial X:=\mathbb C^m$  as above. Adapting the domain of $A_m$ as
\[D(A_m):=\l\{f\in\rW^{1,p}\bigl([0,1],\Cm\bigr): f(1) \in \operatorname{rg}\Psi \r\}\]
and choosing the control operator $B=b:=\Psi v\in\mathbb C^m$  we are in the situation considered in  \cite{EKNS:08} and \cite{BBEAM:13}, see also \cite[Sec.~18.4]{BKR:17}.

Then the \emph{approximate} controllability space  for the network flow problem computed in \cite[Cor.~4.3]{EKNS:08} by our Corollary~\ref{cor:Reach-TE} indeed coincides with the \emph{exact} controllability space.

\begin{corollary}
If $t\ge \min\{m,n\}=:l$ then the exact reachability space of the controlled transport in network problem \eqref{eq:TE-net} equals
\begin{align*}
e\sR_t^{\textrm{BC}}=e\sR^{\textrm{BC}}&={\rL^p[0,1]}\otimes\operatorname{span}\l\{\Psi v,\mathbb B \Psi v,\ldots,\mathbb B^{l-1}\Psi v\r\}\\
& = {\rL^p[0,1]}\otimes \Psi \operatorname{span}\l\{v,\mathbb A v,\ldots,\mathbb A^{l-1} v\r\} .
\end{align*}
\end{corollary}

Note that in big connected networks one usually has $n\le m$, hence the latter space is the more important one for applications.

\smallbreak

Positive control for this problem was already treated in \cite{BBEAM:13} and the approximate positive reachability spaces was computed. However, our approach even yields the \emph{exact} reachability space.

\begin{corollary}
The exact positive reachability space of the controlled transport in network problem \eqref{eq:TE-net} is given by
\begin{align*}
e^+\sR^{\textrm{BC}}&={\rL^p\bigl([0,1],\mathbb R_+\bigr)}\otimes\operatorname{co}\l\{\mathbb B^{k} \Psi v : k\in\mathbb N_0\r\}\\
& = {\rL^p\bigl([0,1],\mathbb R_+\bigr)}\otimes \Psi \operatorname{co}\l\{\mathbb A^k v : k\in\mathbb N_0\r\}.
\end{align*}
\end{corollary}

\subsection{Exact Boundary Controllability for Flows in Networks with Dynamical Boundary Conditions}

In this subsection we investigate exact controllability in the situation of \cite[Sect.~3]{EKKNS:10}. Without going much into details we only introduce the necessary facts to state the problem and to compute the exact reachability space $e\sR_t^{\textrm{BC}}$.

\smallskip
We start from the transport problem in the network introduced in the previous example, but now change the transmission process in the vertices allowing for dynamical boundary conditions.
To encode the structure of the underlying network and the imposed boundary conditions we use the incidence matrices introduced above as well as the weighted  incoming incidence matrix $\Phi_w^{+}$ given by
\[
\bigl(\Phi_w^{+}\bigr)_{ij}:=
  \begin{cases}
   w^+_{ij} & \mbox{if $\stackrel{e_j}{\longrightarrow}v_i$}\ , \\
    0 & \text{otherwise,}
    \end{cases}
\]
for some $0\le w^+_{ij} \le 1$. Defining
\begin{equation}\label{eq:ABrel2}
\mathbb A := \Phi_w^+\Psi  \quad\text{and}\quad \mathbb B:=\Psi\Phi_w^+
\end{equation}
we obtain the adjacency matrices as above (with different nonzero weights).
We mention that the relations \eqref{eq:ABrel} remain valid also in this case.
\smallbreak

We are then interested in the network transport problem with dynamical boundary conditions in $s=1$ considered already in \cite{Sik:05} and \cite[Sect.~3]{EKKNS:10}, i.e.,
\begin{alignat}{2}\label{eq:TE-net-db}
\begin{cases}
\dot x(t,s)=x'(t,s),& s\in[0,1],\ t\ge0,\\
\dot x(t,1)=\mathbb B x(t,0)+u(t)\cdot \Psi v,&t\ge0,\\
x(0,s)=0,& s\in[0,1],\\
\Phi^- x(1,0) = 0.
\end{cases}
\end{alignat}

To embed this example in our setting we introduce
\begin{itemize}
\item the state space $X:={\rL^p\bigl([0,1],\Cm\bigr)}\times\mathbb C^n$ where $1\le p<+\infty$,
\item the boundary space $\partial X:=\mathbb C^m$,
\item the control space $U:=\mathbb C$,
\item the control operator $B:=\Psi v\in\mathbb C^m\simeq\mathcal{L}(U,\partial X)=\mathcal{L}(\mathbb C,\mathbb C^m)$
where $v=v_i$ denotes the $i$-th canonical basis vector of $\mathbb C^n$ meaning that the control acts in the $i$-th vertex of the network,
\item the system operator\footnote{By $\delta_s$ we denote the point evaluation in $s\in[0,1]$, i.e., $\delta_s(f)=f(s)$.}
\begin{align*}
A_m:&=\begin{pmatrix}
\operatorname{diag}\bigl(\frac{d}{ds}\bigr)_{m\times m}&0\\\Phi_w^{+}\delta_0&0
\end{pmatrix}
\quad\text{with domain}\\
D(A_m):&=\l\{\tbinom fd\in\rW^{1,p}\bigl([0,1],\Cm\bigr)\times\mathbb C^n:f(1)\in\operatorname{rg}\Psi\r\},
\end{align*}
\item  the boundary operator $Q:D(A_m)\times\mathbb C^n\to\mathbb C^m$, $Q\binom fd:=\Phi^{-} f(1)-d$,
\item  the operator $A\subsetA_m$ with domain $D(A)=\ker Q$.
\end{itemize}

As is shown in \cite[Prop.~3.4]{EKKNS:10} these spaces and operators satisfy all assumptions of Section~\ref{TAF}.
To proceed we first need to compute the associated Dirichlet operator $Q_\lambda$ and an explicit representation of the semigroup operators $T(t)$ for $t\in[0,1]$.

\begin{lemma}
\makeatletter
\hyper@anchor{\@currentHref}%
\makeatother
\label{lem:Diri-SGR-dBC}
\begin{enumerate}[(i)]
\item For each $0\neq \lambda\in\rho(A)$, the Dirichlet operator $Q_{\lambda}\in \mathcal{L}(\mathbb C^n, {X})$ is given by
\begin{equation}\label{Ql}
   Q_{\lambda} =
   \binom{
   \lambda \varepsilon_{\lambda}\otimes\Psi R(\lambdae^{\lambda},\mathbb A)}
   {\mathbb A R(\lambdae^{\lambda},\mathbb A)}.
\end{equation}
\item The semigroup $(T(t))_{t\ge0}$ generated by $A$ is given by\footnote{We use the notations $\bigl[\binom fd\bigr]_1:=f$ and $\bigl[\binom fd\bigr]_2:=d$ for the canonical projections of $\binom fd\in X$.}
\begin{align}
\l[T(t)\tbinom fd\r]_1(s)&=
\begin{cases}\label{eq:[T(t)]1}
f(t+s)&\text{if }\kern22pt 0\le t<1-s,\\
\mathbb B\, V_{t+s-1} f +\Psi d&\text{if }\,1-s\le t\le1,
\end{cases}
\\
\l[T(t)\tbinom fd\r]_2&=\Phi_w^{+}\, V_t f +d\kern45pt\text{for }0\le t\le1,\label{eq:[T(t)]2}
\end{align}
where
\begin{equation}\label{def:V_s}
V_s f:=\int_0^s f(r)\,dr\quad \text{for }f\in{\rL^p\bigl([0,1],\Cm\bigr)}.
\end{equation}
\end{enumerate}
\end{lemma}

\begin{proof}
Assertion~(i) is proved in \cite[Prop.~3.8]{EKKNS:10}. Equation~\eqref{eq:[T(t)]2} is shown in the proof of \cite[Prop.~3.4.(iii)]{EKKNS:10}.  The statement \eqref{eq:[T(t)]1} for the first coordinate then follows from \cite[Lem.~6.1]{Sik:05}.
\end{proof}

Next  we apply Proposition~\ref{prop:range-B1} to the present situation.

\begin{lemma}\label{lem:strange-dBC}
Let $\lambda\in\rho(A)$. Then for all $0\le\alpha\le1$
\begin{align}\label{eq:strange-dBC}
\bigl[e^{\lambda\alpha}\cdot T(1-\alpha)B_\lambda \bigr]_1(s)&=
\begin{cases}
\lambda\varepsilon_{\la}(1+s)\cdot\PsiR(\lambdae^{\lambda},\mathbb A)\, v&\text{if }0\le s<\alpha,\\
\lambda\varepsilon_{\la}(1+s)\cdot\PsiR(\lambdae^{\lambda},\mathbb A)\, v-\varepsilon_{\la}(s)\cdot \Psi v&\text{if }\alpha\le s\le1.
\end{cases}
\\
\bigl[e^{\lambda\alpha}\cdot T(1-\alpha)B_\lambda \bigr]_2
&=e^{\lambda}\mathbb A R(\lambdae^{\lambda},\mathbb A)\, v
\end{align}
Hence the equality in \eqref{eq:strange1} is satisfied with
\begin{align*}
M&=\binom{\Psi v}{0}\in\mathcal{L}\Bigl({\rL^p[0,1]},{\rL^p\bigl([0,1],\Cm\bigr)}\times \mathbb C^n\Bigr),\
(M u)(s)=\binom{ u(s)\cdot\Psi v}{0}.
\end{align*}
\end{lemma}
\begin{proof}
Using the explicit representations of $Q_\lambda$ and $T(t)$ given in Lemma~\ref{lem:Diri-SGR-dBC} and the relations \eqref{eq:ABrel} we obtain
\begin{align*}
\bigl[e^{\lambda\alpha}\cdot T(1&-\alpha)B_\lambda \bigr]_1(s)=\\
&=e^{\lambda\alpha}\cdot
\begin{cases}
\lambda\varepsilon_{\la}(1-\alpha+s)\cdot\PsiR(\lambdae^{\lambda},\mathbb A)\, v&\text{if }0\le s<\alpha,\\
\lambda\mathbb B\, V_{s-\alpha}\,\varepsilon_{\la}\cdot\PsiR(\lambdae^{\lambda},\mathbb A)\, v+\Psi\mathbb A R(\lambdae^{\lambda},\mathbb A) v&\text{if }\alpha\le s\le1,
\end{cases}
\\
&=e^{\lambda\alpha}\cdot
\begin{cases}
\lambda\varepsilon_{\la}(1-\alpha+s)\cdot\PsiR(\lambdae^{\lambda},\mathbb A)\, v&\text{if }0\le s<\alpha,\\
\bigl(\varepsilon_{\la}(s-\alpha)-1\bigr)\cdot\Psi\,\AAAR(\lambdae^{\lambda},\mathbb A)\, v+\Psi\mathbb A R(\lambdae^{\lambda},\mathbb A) v&\text{if }\alpha\le s\le1,
\end{cases}
\\
&=
\begin{cases}
\lambda\varepsilon_{\la}(1+s)\cdot\PsiR(\lambdae^{\lambda},\mathbb A)\, v&\text{if }0\le s<\alpha,\\
\varepsilon_{\la}(s)\cdot\Psi\bigl(\lambda\elR(\lambdae^{\lambda},\mathbb A)-Id\bigr)\, v&\text{if }\alpha\le s\le1.
\end{cases}
\\
&=\begin{cases}
\lambda\varepsilon_{\la}(1+s)\cdot\PsiR(\lambdae^{\lambda},\mathbb A)\, v&\text{if }0\le s<\alpha,\\
\lambda\varepsilon_{\la}(1+s)\cdot\PsiR(\lambdae^{\lambda},\mathbb A)\, v-\varepsilon_{\la}(s)\cdot \Psi v&\text{if }\alpha\le s\le1.
\end{cases}
\end{align*}
Similarly, for the second coordinate we have
\begin{align*}
\bigl[e^{\lambda\alpha}\cdot T(1-\alpha)B_\lambda \bigr]_2
&=e^{\lambda\alpha}\Bigl(\lambda\Phi_w^{+}\,V_{1-\alpha}\,\varepsilon_{\la}\cdot\PsiR(\lambdae^{\lambda},\mathbb A)\,v+\mathbb A R(\lambdae^{\lambda},\mathbb A)\, v\Bigr)\\
&=e^{\lambda\alpha}\Bigl(\bigl(\varepsilon_{\la}(1-\alpha)-1\bigr)\cdot\AAAR(\lambdae^{\lambda},\mathbb A)\,v+\mathbb A R(\lambdae^{\lambda},\mathbb A)\, v\Bigr)\\
&=e^{\lambda}\mathbb A R(\lambdae^{\lambda},\mathbb A)\, v
\end{align*}
where we used \eqref{eq:ABrel2}.
\end{proof}

We note that by \cite[Prop.~3.5]{EKKNS:10} the states of the controlled flow at time $t\ge0$ are given by the first coordinate of the states in our ``extended'' state space $X={\rL^p\bigl([0,1],\Cm\bigr)}\times\mathbb C^n$. For this reason we also need to compute the first coordinate of $T(1)^k\binom{\Psi g}{0}$.

\begin{lemma}
We have
\begin{equation*}
\l(T(1)\tbinom fd \r)(s)=
\begin{pmatrix}
\Psi\,\Phi_w^{+}\, V_s&\Psi\\\kern10pt \Phi_w^{+}\,V_1&Id
\end{pmatrix}
\binom{f}{d}
=\binom{\mathbb B\,V_s f+\Psi d}{\Phi_w^{+}\, V_1 f +d},
\end{equation*}
where the operator $V_s\in\mathcal{L}\bigl({\rL^p\bigl([0,1],\Cm\bigr)},\rW^{1,p}\bigl([0,1],\Cm\bigr)\bigr)$ is defined in \eqref{def:V_s}.
Moreover, for $k\in\mathbb N_1$ we have
\begin{equation}\label{eq:T(k)-dBC}
\l[T(1)^k\tbinom{\Psi g}0\r]_1 (s)
=\Psi(\mathbb A V_s +\delta_1)^{k-1}\mathbb A\,V_sg
=(\mathbb B V_s+\delta_1)^{k-1}\mathbb B\Psi\,V_sg.
\end{equation}
\end{lemma}

\begin{proof}
The formula for $T(1)$ follows immediately from Lemma~\ref{lem:Diri-SGR-dBC}.(ii). Since $\Psi\mathbb A = \mathbb B\Psi$, it suffices to show the second equality in \eqref{eq:T(k)-dBC}. Obviously this equation holds for $k=1$. To verify it for $k>1$ we note that by \eqref{eq:ABrel} the matrix $\Psi$ is left invertible with left inverse $\Phi^{-}$. Hence, we obtain
\[
\l[T(1)\tbinom fd\r]_2=\Phi^{-}\delta_1 \l[T(1)\tbinom fd\r]_1.
\]
If $\binom fd\in\operatorname{rg} T(1)$ we can write $f=\Psi h$ and the previous equation implies
\begin{align*}
\l[T(1)\tbinom fd\r]_1 (s)
=\l[T(1)\tbinom {\Psi h}{\delta_1h}\r]_1(s)
=\mathbb B V_s\Psi h+\Psi\delta_1 h
=(\mathbb B V_s +\delta_1)f.
\end{align*}
Now assume that \eqref{eq:T(k)-dBC} holds for some $k\ge1$. Then for $\binom fd=T(1)^k\binom{\Psi g}{0}\in\operatorname{rg} T(1)$ we conclude
\begin{align*}
\l[T(1)^{k+1}\tbinom{\Psi g}0\r]_1(s)
&=\l[T(1)\cdot T(1)^{k}\tbinom{\Psi g}0\r](s)\\
&=(\mathbb B V_s +\delta_1)\cdot(\mathbb B V_s+\delta_1)^{k-1}\,\mathbb B\Psi\,V_sg\\
&=(\mathbb B V_s+\delta_1)^{k}\,\mathbb B\Psi\,V_s g.\qedhere
\end{align*}

\end{proof}

The previous two lemmas together with Corollary~\ref{cor:Bn} imply the following result.

\begin{corollary}\label{cor:Bn-dBC}
For $l\in\mathbb N_2$ and $u\in{\rL^p[0,l]}$ we have
\begin{align}
\bigl[\sB_l u\bigr]_1(s)
&=\Psi\biggl( u_0\otimes v+\sum_{k=1}^{l-1}\bigl(\mathbb A V_s+\delta_1\bigr)^{k-1}\, V_s(u_{k}\otimes\mathbb A v)\biggr) \nonumber\\
&=u_0\otimes \Psi v+\sum_{k=1}^{l-1}\bigl(\mathbb B V_s+\delta_1\bigr)^{k-1}\,V_s(u_{k}\otimes\mathbb B\Psi\, v)\label{eq:edBC}
\end{align}
where $u_k\in{\rL^p[0,1]}$ is defined as in \eqref{eq:def-uk}.
\end{corollary}

Using this explicit representation of the controllability map we now compute the exact reachability space for the control problem given in \eqref{eq:TE-net-db}.

\begin{corollary}\label{cor:Reach-dBC}
If $t\ge\min\{m,n\}=:l$ then the exact reachability space of the controlled flow with dynamic boundary conditions \eqref{eq:TE-net-db} is given by\footnote{Here we define $\mathrm{W}^{0,p}[0,1]:={\rL^p[0,1]}$.}
\begin{align*
\bigl[e\sR_t^{\textrm{BC}}\bigr]_1
&\subseteq
\l\{
\Psi\sum_{k=0}^{l}\l(u_{k}\otimes\mathbb A^k\, v\r):u_k\in\rW^{k,p}[0,1]\text{ for }0\le k\le l
\r\}\\
&=
\l\{
\sum_{k=0}^{l}\l(u_{k}\otimes\mathbb B^k\Psi\, v\r):u_k\in\rW^{k,p}[0,1]\text{ for }0\le k\le l
\r\}.
\end{align*}
\end{corollary}

\begin{proof}
The equality of the two sets on the right-hand-side follows immediately from \eqref{eq:ABrel}.
To show the inclusion in the second set we combine Corollaries~\ref{cor:Reach} and  \ref{cor:Bn-dBC}.
First observe, that for the operators $\mathbb B$, $V_s$,  and $\delta_1$ we have
\[
\mathbb B V_s f = V_s \mathbb B f, \quad \mathbb B \delta_1 f = \delta_1\mathbb B f,\quad
\delta_1 V_s f = V_1f
\]
for every $f\in{\rL^p\bigl([0,1],\Cm\bigr)}$ while
\[ \delta_1^k  f = \delta_1f =f(1)\quad\text{for } k\ge 1.\]
So, when expanding $(\mathbb B V_s+\delta_1)^{k-1}V_s$ we can rearrange the terms to obtain expressions of the form
\[\alpha_i \mathbb B^i V_{s_1}\cdots V_{s_{i+1}},\quad 0\le i\le k-1,\]
where $\alpha_i$ are scalar coefficients and $s_j\in\{s,1\}$, $1\le j\le i+1$.
Next, for arbitrary $u\in{\rL^p[0,1]}$
and $0\le k\le l$ we have
\[ V_{s_1}\cdots V_{s_k} u \in \rW^{k,p}[0,1], \quad s_j\in\{s,1\}, 1\le j\le k. \]
Combining these facts we obtain the desired result by considering \eqref{eq:edBC} for all $u\in{\rL^p[0,l]}$.
\end{proof}

From the previous Corollary we immediately obtain the following result which improves \cite[Thm.~3.10]{EKKNS:10} and shows that $ \bigl[a\sR_t^{\textrm{BC}}\bigr]_1$ is constant for $t\ge\min\{m,n\}=:l$.

\begin{corollary}\label{cor:Reach-approx-dBC}
If $t\ge\min\{m,n\}=:l$ then the approximate controllability space of the controlled flow with dynamic boundary conditions \eqref{eq:TE-net-db} is given by
\begin{align*}
\bigl[a\sR_t^{\textrm{BC}}\bigr]_1
&={\rL^p[0,1]}\otimes\operatorname{span}\l\{\Psi v,\mathbb B\,\Psi v,\ldots,\mathbb B^{l-1}\Psi v\r\}\\
&={\rL^p[0,1]}\otimes\Psi\operatorname{span}\l\{v,\mathbb A v,\ldots,\mathbb A^{l-1}v\r\}.
\end{align*}
\end{corollary}

\smallskip
In the same manner as before we also obtain the following result on positive controllability.

\begin{corollary}\label{cor:Reach-pos-dBC}
The approximate positive controllability space of the controlled flow with dynamic boundary conditions \eqref{eq:TE-net-db} is given by
\begin{align*}
\bigl[a^+\sR^{\textrm{BC}}\bigr]_1
&={\rL^p[0,1]}\otimes\overline{\operatorname{co}}\,\l\{\mathbb B^k\Psi v : k\in\mathbb N_0\r\}\\
&={\rL^p[0,1]}\otimes\Psi\,\overline{\operatorname{co}}\,\l\{\mathbb A^k v: k\in\mathbb N_0\r\}.
\end{align*}
\end{corollary}

\section*{Conclusion}
Using a new characterization of admissible boundary control operators (see Proposition~\ref{prop:range-B1}) we are able to describe explicitly the  \emph{exact} reachability space of the abstract boundary control system  $\Sigma_{\textrm{BC}}(\Am,B,Q)${}, cf. \eqref{ACPBC}. Moreover, this approach allows us also to determine the \emph{positive} reachability space obtained allowing only positive control functions. Our results generalize and improve the ones obtained in the former works \cite{BBEAM:13, EKNS:08, EKKNS:10} where only approximate controllability or positive controllability under quite restrictive assumptions are studied.

Functional Analysis (math.FA); Optimization and Control (math.OC)