One-off geophysical detection of chlorinated DNAPL during remediation of an industrial site: a case study

1.   Introduction

Symmetries linked to Lie groups appear naturally in many control systems problems ^{[1,2,3,4,5,6,7,8,9,10,11,12,13]}. Numerical integrators evolving on Lie groups are commonly employed to improve its accuracy, as well as to avoid singularities by working with coordinate-free expressions in the Lie algebra associated to the Lie group, on which the system evolves. Typical tasks include trajectory tracking and estimation algorithms for the pose of mechanical systems.

Optimization problems on Lie groups ^[14] have applications in control engineering. Indeed, there are many examples of robotic systems that possess invariant quantities steaming from existing symmetries. Invariant quantities can be used as leverage to reduce the complexity of the equations, by projecting them into lower dimensional spaces. Techniques taking advantage of symmetries have been studied in ^{[15,16,17,18,19,20]}, among many others, mainly for applications in robotic and aerospace engineering and, in particular, for spacecraft attitude control and underwater vehicles ^[21]. Most of the applications provided in the literature focus on the single-agent situation and only a few works explore symmetry reduction in a multi-agent scenario (see for instance ^{[22,23,24,25]}). In this work, we employ symmetry reduction to study optimal control problems with broken symmetry for multi-agent systems on Lie groups, while agents avoid collisions and obstacles in the configuration space.

In ^[26], the authors studied symmetry reduction in optimal control problems with broken symmetry for single-agent systems. In this work, we advance on the results of ^[26] by considering a multi-agent scenario. Hence, a generalization of the previous variational principle and reduction by symmetries performance are needed. The results in this paper are the Lagrangian/variational counterpart of those in ^[22]; we also develop a discrete-time version of the results. At the continuous-time side, we obtain the Euler–Poincaré equations from a constrained variational principle while after discretizing the variational principle in time, we obtain the discrete-time Lie–Poisson equations.

The approach in this work to obtain the reduced optimality conditions is a generalization of the one followed in ^[26], and is as follows. First, we consider a collection of free agents evolving on a Lie group $G$ and an artificial potential $V^0$, used to prevent collisions with a fixed obstacle, which is not symmetry invariant. At the same time, we consider a representation of $G$ on a dual vector space $W^{*}$ and a graph $\mathcal{G}$ to describe interactions between neighboring agents. Though the artificial potential is not symmetry invariant, the interaction between neighbouring agents is, which is the standard situation when the interaction is done through a potential dependent on relative positions. Coupling the artificial potential a parameter depending on vectors in $W^{*}$, we obtain a symmetry invariant potential function under the action of $G$. Loosely speaking, the agents are coupled with vectors in $W^{*}$ that are acted by the representation. The associated action considered at this stage restores the full Lie group symmetry in the cost function from our optimal control problem and allows us to apply the semi-direct product reduction theory ^{[3,27,28,29,30,31,32]}, to obtain the corresponding Euler-Poincaré system on the semi-direct product Lie algebra $\mathfrak{g}\ltimes W^{*}$ at each node. This gives rise to a new system that finds no analogs in classical reduced-order models in optimal control of mechanical systems.

The paper is organized as follows, in Section 2, we introduce some preliminaries about geometric mechanics on Lie groups and Lie group actions. In Section 3.1, we present the problem together with a motivating example. In Section 4, we study the Euler–Poincaré reduction of optimal control problems for left-invariant multi-agent control systems on Lie groups with partial symmetry breaking cost functions. Furthermore, we consider the discrete-time framework and obtain the discrete-time Lie–Poisson equations in Section 5. In Section 6, an example is considered to illustrate the theory. Finally, some concluding remarks are given in Section 7.

2.   Lie group actions and representations

Consider a manifold $Q$ to be the configuration space of a mechanical system and let it be differentiable of dimension $d$ which can be described locally by coordinates $q = (q^1, \ldots, q^d)$. Denote by $TQ$ the tangent bundle of $Q$, with local coordinates described by positions and velocities for the system, $v_q = (q^1, \ldots, q^d, \dot{q}^{1}, \ldots, \dot{q}^{d})\in TQ$ with $\hbox{dim}(TQ) = 2d$. Let $T^{*}Q$ be its cotangent bundle, locally described by the positions and the momenta for the system, i.e., $(q, p)\in T^{*}Q$ with $\hbox{dim}(T^{*}Q) = 2d$. The tangent space at $q\in Q$ has a vector space structure, and it is denoted as $T_{q}Q$. The cotangent space at $q\in Q$ is just the dual space of $T_{q}Q$ and is denoted as $T_{q}^{*}Q$. The dynamics of a mechanical system is described by the equations of motion determined by a Lagrangian function $L:TQ\to \mathbb{R}$ given by $L(q, \dot{q}) = K(q, \dot{q})-V(q)$, where $K:TQ\to \mathbb{R}$ denotes the kinetic energy and $V:Q\to \mathbb{R}$ the potential energy of the system. The equations of motion are given by the Euler-Lagrange equations ${\frac{d}{dt}\bigg(\frac{\partial L}{\partial\dot{q}^i}\bigg) = \frac{\partial L}{\partial q^i}, \; \; i = 1, \dots, d}$, which determine a system of second-order differential equations. In case, the configuration space of the system is a Lie group, Euler-Lagrange equations can be reduced to a system of first-order ordinary differential equations.

Definition 2.1. Let $G$ be a Lie group and $Q$ a smooth manifold. A left-action of $G$ on $Q$ is a smooth map $\Phi:G\times Q\to Q$, such that $\Phi(\bar{e}, g) = g$ and $\Phi(h, \Phi(g, q)) = \Phi(hg, q)$ for all $g, h\in G$ and $q\in Q$, where $\bar{e}$ is the identity of the group $G$ and the map $\Phi_g:Q\to Q$ given by $\Phi_g(q) = \Phi(g, q)$ is a diffeomorphism for all $g\in G$.

Definition 2.2. A function $f:Q\to \mathbb{R}$ is called left invariant under $\Phi_g$ if $f\circ\Phi_g = f$ for any $g\in G$.

For a finite dimensional Lie group $G$, its Lie algebra $\mathfrak{g}$ is defined as the tangent space to $G$ at the identity, $\mathfrak{g}: = T_{\bar{e}}G$. Let $L_g:G\to G$ be the left translation of the element $g\in G$ given by $L_g(h) = gh$, where $h\in G$. $L_g$ is a diffeomorphism on $G$ and a left-action of $G$ on $G$ ^[33]. Its tangent map (i.e, the linearization or tangent lift) is denoted by $T_{h}L_{g}:T_{h}G\to T_{gh}G$. In a similar way, the cotangent map (cotangent lift), is denoted by $(T_{h}L_g)^{*}$, and is defined as the dual map of the tangent lift denoted by $T_{h}^{*}L_{g}:T^{*}_{gh}G\to T^{*}_{h}G$, and determined by the relation $\langle(T_hL_g)^{*}(\alpha_{gh}), y_h\rangle = \langle\alpha_{gh}, (T_hL_g)y_h\rangle$, $y_h\in T_{h}G$, $\alpha_{gh}\in T^{*}_{gh}G$. As it is known, the tangent and cotangent lift of a Lie group action are Lie group actions as well. Here, $\langle\cdot, \cdot\rangle:W^{*}\times W\to\mathbb{R}$ with $W$ a finite-dimensional vector space denotes the usual natural pairing between vectors and co-vectors, and defined by $\langle y, x\rangle: = y\cdot x$ for $y\in W^{*}$ and $x\in W$. For a matrix Lie algebra $\langle y, x\rangle = y^{T}x$ (see ^[33], Section 2.3). Using the natural pairing between vectors and co-vectors, for $g, h\in G$, $y\in\mathfrak{g}^{*}$ and $x\in\mathfrak{g}$, we write $\langle T^{*}_{g}L_{g^{-1}}(y), T_{{\bar{e}}}L_{g}(x)\rangle = \langle y, x\rangle$.

Denote by $\text{ad}^{*}:\mathfrak{g}\times\mathfrak{g}^{*}\to\mathfrak{g}^{*}$, $(\xi, \mu)\mapsto \text{ad}^{*}_{\xi}\mu$ the co-adjoint operator, defined by $\langle \text{ad}_{\xi}^{*}\mu, \eta\rangle = \langle\mu, \text{ad}_{\xi}\eta\rangle$ for all $\eta\in\mathfrak{g}$, where the $\text{ad}:\mathfrak{g}\times \mathfrak{g}\to \mathfrak{g}$ denotes the adjoint operator on $\mathfrak{g}$ given by the Lie-bracket, i.e., $\text{ad}_{\xi}\eta = [\xi, \eta]$, $\xi, \eta\in \mathfrak{g}$. We also define the adjoint action of $G$ on $\mathfrak{g}$, denoted by $\hbox{Ad}_{g}:\mathfrak{g}\to\mathfrak{g}$ and, in the case the Lie algebra is a matrix Lie group, it is given by $\hbox{Ad}_{g}\chi: = g\chi g^{-1}$, where $\chi\in\mathfrak{g}$, and the co-adjoint action of $G$ on $\mathfrak{g}^{*}$, denoted by $\hbox{Ad}_{g}^{*}:\mathfrak{g}^{*}\to\mathfrak{g}^{*}$, and given by $\langle\hbox{Ad}_{g}^{*}\alpha, \xi\rangle = \langle\alpha, \hbox{Ad}_{g}\xi\rangle$ with $\alpha\in \mathfrak{g}^{*}$.

For $q\in Q$, the isotropy (or stabilizer or symmetry) group of $\Phi$ at $q$ is given by $G_q: = \{g\in G|\Phi_g(q) = q\}\subset G$. Since the map $\Phi(\cdot, q) = \Phi^q:G\to Q$ is a continuous, $G_q = (\Phi^q)^{-1}(q)$ is a closed subgroup, and thus, a Lie subgroup of $G$ (see ^[34] Sec. 9.3 for instance).

Example 1. Consider the special Euclidean group $SE(2)$ of rotations and translations on the plane. Elements on $SE(2)$ can be expressed by transformations of $\mathbb{R}^2$ of the form $z \mapsto Rz+r$, with $r\in \mathbb{R}^2$ and $R\in SO(2)$. This transformation can be represented by $g = (R, r)$, for ${R = \left(\begin{array}{cc} \cos \theta & -\sin \theta\\ \sin \theta & \cos \theta \\ \end{array} \right)}$ and $r = [x, y]^{T}$. The composition law is $(R, r)\cdot(S, s) = (RS, Rs+r)$ with identity element $(I, 0)$ and inverse $(R, r)^{-1} = (R^{-1}, -R^{-1}r)$. Under this composition rule, $SE(2)$ has the structure of the semidirect product Lie group $SO(2)\ltimes\mathbb{R}^2$. Here, as usual in the literature, we denote by $\ltimes$ the semidirect product of Lie groups.

The Lie algebra $\mathfrak{se}(2)$ of $SE(2)$ is determined by ${\mathfrak{se}(2) = \Big{\{}\left(\begin{array}{cc} A & b\\ 0 & 0 \\ \end{array} \right)\Big{|} A\in \mathfrak{so}(2)\simeq\mathbb{R}, \, b\in \mathbb{R}^2\Big{\}}}$. In the following, for simplicity, we write $A = -a \mathbb{J}$, $a\in\mathbb{R}$, where $\mathbb{J} = \left(\begin{array}{cc} 0 & 1\\ -1 & 0 \\ \end{array} \right).$ Therefore, we denote $\xi = (a, b)\in\mathfrak{se}(2)$. The adjoint action of $SE(2)$ on $\mathfrak{se}(2)$ is given by $\hbox{Ad}_{(R, r)}(a, b) = (a, a\mathbb{J}r+Rb)$ (see ^[35], pp. 153 for instance), so, $\hbox{Ad}_{(R, r)^{-1}}(a, b) = (a, R^{T}(b-a\mathbb{J}r))$.

Next, we provide the infinitesimal description of a Lie group action, that will be an important concept in the remainder of the paper.

Definition 2.3. Given a Lie group action $\Phi:G\times Q\to Q$, for $\xi\in\mathfrak{g}$, the map $\Phi_{\hbox{exp}(t\xi)}:Q\to Q$ is a flow on $Q$, where $exp$ is the exponential map of $G$. The corresponding vector field on $Q$, given by $\xi_Q(q): = \frac{d}{dt}\mid_{t = 0}\Phi_{\hbox{exp} (t\xi)}(q)$ is called the infinitesimal generator of the action corresponding to $\xi$.

Definition 2.4. Let us denote the set of vector fields on a Lie group $G$ by $\mathfrak{X}(G)$. A left invariant vector field is an element $X$ of $\mathfrak{X}(G)$ such that $T_{h}L_g(X(h)) = X(L_g(h)) = X(gh)$ $\forall\, g, h\in G$.

Especially, for $h = {\bar{e}}$, we have that a vector field $X$ is left-invariant if $X(g) = T_{{\bar{e}}}L_{g}\xi$ for $\xi = X({\bar{e}})\in\mathfrak{g}$. As a consequence, left invariant vector fields are identified with its corresponding element of $\mathfrak{g}$. Thus, for every $g\in G$ and $\xi \in \mathfrak{g}$, we define the left-invariant vector field as ${\xi}(g): = T_{{\bar{e}}}L_g(\xi)$.

Example 2. Consider the Euclidean Lie group $\mathbb{R}^{n}$ with the sum as group operation. For all $g\in\mathbb{R}^{n}$, the left translation $L_g$ is the usual translation on $\mathbb{R}^n$, that is, $L_g(h) = g+h$, $h\in\mathbb{R}^n$. So that, the tangent map to $L_g$ is the identity map on $\mathbb{R}^n$, that is, $T_{0}L_g = \hbox{id}_{T_{0}\mathbb{R}^n}$, where we are using that $T_{h}\mathbb{R}^n\simeq\mathbb{R}^n$ for all $h\in\mathbb{R}^n$, since $\mathbb{R}^n$ is a vector space. Therefore, left-invariant vector fields are constant vector fields, that is, $X = v_1\frac{\partial}{\partial x_1}+\ldots+v_n\frac{\partial}{\partial x_n}$ for $v = (v_1, \ldots, v_n)\in T_{0}\mathbb{R}^n$ and $x = (x_1, \ldots, x_n)\in\mathbb{R}^n$.

Consider a Lie group $G$, a vector space $W$ and the representation of $G$ on $W$, given by $\rho:G\times W\to W$, $(g, v)\mapsto\rho_g(v)$, which is a left action, and it is defined by the relation $\rho_{g_1}(\rho_{g_2}(v)) = \rho_{g_1g_2}(v)$, $g_1, g_2\in G$. Its dual is given by $\rho^*:G\times W^*\to W^*,$ $(g, \alpha)\mapsto\rho^*_{g^{-1}}(\alpha),$ satisfying $\langle\rho^*_{g^{-1}}(\alpha), v\rangle = \langle\alpha, \rho_{g^{-1}}(v)\rangle$, and it is also a left action of $G$ now on $W^{*}$.

The infinitesimal generator of the left action of $G$ on $W$ is $\rho': \mathfrak{g}\times W\to W,$ $(\xi, v)\mapsto\rho'(\xi, v) = \frac{d}{dt}|_{t = 0}\rho_{\exp{(t\xi)}}(v).$ For every $v\in W$, consider the linear transformation $\rho'_v: \mathfrak{g}\to W,$ $\xi\mapsto\rho'_v(\xi) = \rho'(\xi, v)$ and its dual $\rho'^*_v: W^*\to \mathfrak{g}^*,$ $\alpha\mapsto\rho'^*_v(\alpha).$ The last transformation defines the momentum map ${{\bf{J}}}_W:W\times W^*\to \mathfrak{g}^*,$ $(v, \alpha)\mapsto {{\bf{J}}}_W(v, \alpha): = \rho'^*_v(\alpha)$, such that for every $\xi\in \mathfrak{g}$, ${\langle {{\bf{J}}}_W(v, \alpha), \xi\rangle = \langle\rho'^*_v(\alpha), \xi\rangle = \langle \alpha, \rho'_v(\xi)\rangle = \langle\alpha, \rho'(\xi, v)\rangle}$.

For $\xi\in \mathfrak{g}$, consider the map $\rho'_\xi:W\to W$, $v\mapsto\rho'_\xi(v) = \rho'(\xi, v)$ and its dual $\rho'^*_\xi:W^*\to W^*$, $\alpha\mapsto\rho'^*_\xi(\alpha)$, such that $\langle\rho'^*_\xi(\alpha), v\rangle = \langle\alpha, \rho'_\xi(v)\rangle$. Then, this satisfies ${\langle{\bf{J}}_W(v, \alpha), \xi\rangle = \langle\alpha, \rho'(\xi, v)\rangle = \langle\alpha, \rho'_\xi(v)\rangle = \langle\rho'^*_\xi(\alpha), v\rangle}$. See ^[33] and ^[34] for more details on the momentum map.

Example 3. Let $W = \mathfrak{g}$ and $\rho$ the adjoint representation of $G$ on $W$, i.e., $\rho_g = \hbox{Ad}_g$, for any $g\in G$. So, for $\alpha\in W^{*} = \mathfrak{g}^{*}$, $\rho^{*}$ is the coadjoint representation of $G$ on $W^{*}$, i.e., $\rho_{g}^{*} = \hbox{Ad}_{g}^{*}$. We also have that the infinitesimal generator for the adjoint representation is $\rho_{\xi}' = \hbox{ad}_{\xi}$ for any $\xi\in W$ (see ^[33], Def. $6.4$, pp. $225$), and it follows that $\langle{\bf{J}}_{W}(x, \alpha), \xi\rangle = \langle\alpha, \text{ad}_{\xi}x\rangle = \langle\alpha, - \text{ad}_{x}\xi\rangle = \langle- \text{ad}_{x}^{*}\alpha, \xi\rangle$, for $x\in W = \mathfrak{g}$, which gives ${\bf{J}}_{W}(x, \alpha) = - \text{ad}_{x}^{*}\alpha$.

Similarly, if now $W = \mathfrak{g}^{*}$ and the coadjoint representation of $G$ on $W$ is $\rho$, i.e., $\rho_g = \hbox{Ad}^{*}_{g^{-1}}$, for any $g\in G$, then $\rho_{\xi}' = - \text{ad}_{\xi}^{*}$ for any $\xi\in\mathfrak{g}$. So, if $x, \xi\in\mathfrak{g}$ and $\alpha\in \mathfrak{g}^{*}$, it follows that $\langle{\bf{J}}_{W}(\alpha, x), \xi\rangle = \langle- \text{ad}_{\xi}^{*}\alpha, x\rangle = \langle \alpha, \text{ad}_{x}\xi\rangle = \langle \text{ad}_{x}^{*}\alpha, \xi\rangle$, which gives ${\bf{J}}_{W}(\alpha, x) = \text{ad}_{x}^{*}\alpha$.

3.   Problem formulation

3.1.   Left-invariant multi-agent control systems

Denote, by $\mathcal{N}$, a set consisting of $s\geq 2$ free agents, and by $\mathcal{E}\subset\mathcal{N}\times\mathcal{N}$ the set describing the interaction between them. The neighboring relationships are illustrated by an undirected graph $\mathcal{G} = (\mathcal{N}, \mathcal{E})$, where $\mathcal{N}$ is the set of vertices and $\mathcal{E}$ the set of edges for $\mathcal{G}$. We further assume $\mathcal{G}$ is static and connected. For every agent $i\in \mathcal{N}$ the set $\mathcal{N}_i = \{j\in\mathcal{N}:(i, j)\in\mathcal{E}\}$ denotes the neighbors of that agent. The agent $i\in\mathcal{N}$ evolves on an $n$-dimensional Lie group $G$, and its configuration is denoted by $g_i\in G$. We denote by $G^s$ and by $T_{e}G^s = : \mathfrak{g}^s$ the cartesian products of $s$ copies of $G$ and $\mathfrak{g}$, respectively, where $e = (\overline{e}, \overline{e}, \dots, \overline{e})$ is the identity of $G^s$, and $\overline{e}$ the identity element of $G$.

For each agent $i\in\mathcal{N}$, there is an associated left-invariant control system described by the kinematic equations

where $g_i(t)\in C^1([0, T], G), \; T\in\mathbb{R}$ fixed, $u_i\in \mathfrak{g}$ is the control input and $g^i_0\in G$ is considered as the initial state condition. Letting $\dim \mathfrak{g} = n$, we can have as many basis for the Lie algebra as the number of agents. Thus, for each $i\in \mathcal{N}$, we have $\mathfrak{g} = \hbox{span}\{e_1, e_2, \ldots, e_n\}$ and the control inputs may be described by coordinates $u_i = [u_i^1, u_i^2, \dots, u_i^m]^T$, where $u_i(t)\in C^1([0, T], \mathfrak{g})$, with $m\leq n$. Hence, the control input for each agent is given by ${u_i(t) = e_0^i+\sum\limits_{k = 1}^{m}u^k_i(t)e_k}$, where $e_0^i\in \mathfrak{g}.$ Thus, the left-invariant control systems (3.1) for each agent $i\in\mathcal{N}$ can be written as

Note that the class of control systems described by (3.2) belongs to the family of underactuated control systems.

3.2.   Motivating example

Following with Example 1, consider the agents $i\in\mathcal{N}$ and $j\in\mathcal{N}_i$ represented as $g_k = (R_k, r_k)$, $k\in\{i, j\}$. Note that $g^{-1}_ig_j = (R_i^{T}R_j, -R_i^{T}(r_i-r_j))$, then, $\hbox{Ad}_{g_i^{-1}g_j}(1, 0) = (1, \mathbb{J}R_i^{T}(r_j-r_i))$. The inner product on $\mathfrak{se}(2)$ is given by $\langle\xi_1, \xi_2\rangle = \hbox{tr}(\xi_1^T\xi_2)$ for $\xi_1, \xi_2\in\mathfrak{se}(2)$ and, hence, the norm is given by $||\xi|| = \sqrt{\hbox{tr}(\xi^T\xi)},$ for any $\xi\in\mathfrak{se}(2)$. For $\xi = (a, b)\in\mathfrak{se}(2)$ we can write the norm of $\xi$ as $||(a, b)|| = \sqrt{2a^2+b^{T}b}$. Therefore, $||\hbox{Ad}_{g_i^{-1}g_j}(1, 0)||^{2} = 2+|\mathbb{J}R_i^{T}(r_j-r_i)|^{2} = 2+|r_j-r_i|^2$, where we have used that $R, \mathbb{J}\in SO(2)$ for the last equality. Hence, it follows that $|r_i-r_j| = \sqrt{||\hbox{Ad}_{g_i^{-1}g_j}(1, 0)||^{2}-2}$.

The previous computation shows that, if the interaction between agents is determined by a function depending on the distances between them, that is, $V_{ij}:G\times G\to\mathbb{R}$, is such that $V_{ij}(g_i, g_j) = V(|r_i-r_j|)$ for some $V:\mathbb{R}_{\geq 0}\to\mathbb{R}$; then, $V_{ij}$ is $SE(2)$-invariant, that is $V_{ij}(hg_i, hg_j) = V_{ij}(g_i, g_j)$. An alternative reasoning of this invariance property has been shown in ^[36].

Next, suppose that we wish to write the distance from an arbitrary point $r\in\mathbb{R}^2$ to a fixed point $x_0\in\mathbb{R}^2$ in terms of the adjoint action. Consider $\xi_0 = (1, \mathbb{J}x_0)\in\mathfrak{se}(2)$. Then, for any $(R, r)\in SE(2)$, we have $\hbox{Ad}_{(R, r)^{-1}}(1, \mathbb{J}x_0) = (1, R^{T}\mathbb{J}(x_0-r))$, and, therefore, $||\hbox{Ad}_{(R, r)^{-1}}(1, \mathbb{J}x_0)||^2 = 2+|x_0-r|^2$. Next, assume we have an obstacle avoidance function $V_i^{0}:SE(2)\to\mathbb{R}$ for each agent $i\in\mathcal{N}$ which can be written as $V_i^{0}(R_i, r_i) = V(|r_i|)$ with $V:\mathbb{R}_{\geq 0}\to\mathbb{R}$. Note that under this assumption, $V$ may be chosen arbitrarily. Then, $V_i^0$ is not $SE(2)$-invariant, but it is $SO(2)$ invariant, i.e., $V_i^{0}(R_{0}R_i, r_i) = V_i^{0}(R_i, r_i)$ for any $R_0\in SO(2)$. Note also that $SO(2)$ is the isotropy group for the coadjoint action, that is, $SO(2)\simeq\{g\in SE(2)\mid Ad_{g}(\xi_{0}) = \xi_{0}\}$, therefore, the obstacle avoidance potential functions $V_i^{0}$ are invariant under the left action of the isotropy group.

In this situation, one can redefine the potential function $V_i^{0}$ to make it $SE(2)$-invariant as follows. Consider $x_0 = 0$, so, $\xi_0 = (R_0, x_0) = (-\mathbb{J}, 0)\in\mathfrak{se}(2)$, and $\|\hbox{Ad}_{(R_i, r_i)^{-1}}\xi_0\| = 2+|r_i|^2$. Hence, ${V_i^{0}(R_i, r_i) = V(|r_i|) = V\left(\sqrt{||\hbox{Ad}_{(R_i, r_i)^{-1}}\xi_0||^2-2}\right)}$. This gives a motivation to define an extended obstacle avoidance function $V_{i, 0}^{ \text{ext}}:G\ltimes W^{*}\to\mathbb{R}$ with $G\ltimes W^{*} = SE(2)\ltimes\mathfrak{se}(2)$ as ${V_{i, 0}^{ \text{ext}}(g_i, \xi): = V\left(\sqrt{||\hbox{Ad}_{g_i^{-1}}\xi||^2-2}\right)}$.

Note that the extended obstacle avoidance function possesses an $SE(2)$-symmetry (i.e., $V_{i, 0}^{ \text{ext}}$ is invariant under a left action of $SE(2)$) since $V_{i, 0}^{ \text{ext}}\circ\tilde{\Phi}(h, (g, \xi)) = V_{i, 0}^{ \text{ext}}(L_h g, \hbox{Ad}_{h}\xi) = V_{i, 0}^{ \text{ext}}(g, \xi)$, for any $h\in SE(2)$, with left action $\tilde{\Phi}$ given by

3.3.   Problem setting

The problem under study consists on finding reduced necessary conditions for optimality in an optimal control problem for (3.1) (or equivalently (3.2)). These solution curves should minimize a cost function and prevent collisions among agents, while they should also avoid static obstacles in the workspace.

Problem (collision and obstacle avoidance): Find reduced optimality conditions on $g(t) = (g_1(t), \dots, g_s(t))\in G^s$ and the controls $u(t) = (u_1(t), \dots, u_s(t))\in \mathfrak{g}^s$ avoiding collision among the agents and obstacles (which will be defined shortly) in the workspace, and such that $(g(t), u(t))\in G^s\times \mathfrak{g}^s$ minimize the following cost function,

subject to the kinematics $\dot{g}_i(t) = T_{\overline{e}_i}L_{g_i(t)}(u_i(t))$, boundary conditions $g(0) = (g_1(0), \dots, g_n(0)) = (g_1^0, \dots, g_s^0)$ and $g(T) = (g_1(T), \dots, g_s(T)) = (g_1^T, \dots, g_s^T)$ with $T\in\mathbb{R}^{+}$ the final time, and under the following assumptions:

$({\rm{i}})$ There is a left representation $\rho$ of $G$ on a vector space $V$.

$({\rm{ii}})$ $C_i: G\times \mathfrak{g}\to\mathbb{R}$ are $G$-invariant functions for each $i\in\mathcal{N}$ (under a left action of $G$ on $G\times \mathfrak{g}$, which is given below) and are also differentiable almost everywhere.

$({\rm{iii}})$ $V_{ij}: G\times G\to\mathbb{R}$ (collision avoidance potential functions) satisfying $V_{ij} = V_{ji}$ are $G$-invariant functions under $\Phi$, defined by

i.e., $V_{ij}\circ\Phi_{g} = V_{ij}$, for any $g \in G$, that is, $V_{ij}(L_{g}(g_i), L_g(g_j)) = V_{ij}(g_i, g_j)$, for any $(g_i, g_j)\in \mathcal{E}$, $j\in\mathcal{N}_i$ and they are also differentiable almost everywhere.

$({\rm{iv}})$ $V_{i}^{0}: G\to\mathbb{R}$ (obstacle avoidance potential functions) are not $G$-invariant functions and they are also differentiables almost everywhere, for $i\in\mathcal{N}$.

$({\rm{v}})$ The obstacle avoidance functions $V_{i}^{0}$ depend on a parameter $\alpha_{0}\in W^{*}$. Hence, we can define the extended potential function as $V_{i, 0}^{ \text{ext}}: G\times W^{*}\to\mathbb{R}$, with $V_{i, 0}^{ \text{ext}}(\cdot, \alpha_{0}) = V_{i}^{0}$, by making the parameter evolve - due to the group action - with initial value $\alpha_0$.

$({\rm{vi}})$ The extended obstacle avoidance functions are $G$-invariant under $\tilde{\Phi}$, defined by

where $\rho_{g^{-1}}^{*}\in W^{*}$ is the adjoint of $\rho_{g^{-1}}\in W$, i.e., $V_{i, 0}^{ \text{ext}}\circ\tilde{\Phi}_{g} = V_{i, 0}^{ \text{ext}}$, for any $g \in G$, or $V_{i, 0}^{ \text{ext}}(L_{g}(h), \rho_{g^{-1}}^{*}(\alpha)) = V_{i, 0}^{ \text{ext}}(g, \alpha)$ where $\alpha\in W^{*}$.

$({\rm{vii}})$ The obstacle avoidance potential functions are invariant under the left action of the isotropy group

Note that $G\times \mathfrak{g}$ is a trivial vector bundle over $G$ and, in order to respect assumption (ⅱ), we define the left action of $G$ on $G\times \mathfrak{g}$ as follows,

We further assume that each $C_i: G\times \mathfrak{g}\to\mathbb{R}$ is $G$-invariant under (3.8), i.e., $C_i\circ\Psi_{g} = C_i$, for any $g\in G$. In addition, each agent $i\in\mathcal{N}$ occupies a disk of radius $\overline{r}$ on $G$. This radius is consider to be small enough in order that all agents can be packed on $G$ and, hence the potential functions are well defined and feasible for ${\bf{d}}(g_i(t), g_j(t)) > 2\overline{r}$ for all $t$, where ${\bf{d}}(\cdot, \cdot):G\times G\to\mathbb{R}$ denotes a distance function on $G$.

4.   Euler-Poincaré reduction for optimal control with broken symmetries

We next study reduced optimality conditions for extrema for the OCP. We address the problem as a constrained variational problem and obtain the Euler-Poincaré equations that extremal must satisfy in Theorem 4.1 and Proposition 4.2.

The optimal control problem (3.4) can be resolved as a constrained variational problem by utilizing the Lagrangian multipliers $\lambda_{g_i} = T^*_{g_i}L_{g_i^{-1}}(\lambda_i)\in T^*_{g _i}G$, where $\lambda_i\in C^1([0, T], \mathfrak{g}^*)$ into the cost functional.

Let $\{e^1, \dots, e^m, e^{m+1}, \dots, e^n\}$ be a basis of $\mathfrak{g}^{*}$ dual of the basis $\{e_{i}\}$ of $\mathfrak{g}$. Then ${\lambda_i = \sum\limits_{k = m+1}^n\lambda^i_k e^k}$, where $\lambda_k^i$ are the components of the vector $\lambda_i$ in the given basis of the Lie algebra $\mathfrak{g}$. Thus, we define the Lagrangian $L:G^s\times \mathfrak{g}^s\times(T^*_{g_i}G)^s\to\mathbb{R}$ by

where $V_{i}:G^{s}\rightarrow \mathbb{R}$,

By assumption (v), the obstacle avoidance potential functions $V^0_i:G\to\mathbb{R}$ depend on a parameter $\alpha_0\in W^{*}$, so we can extend {them} to $V_{i, 0}^{\text{ext}}:G\times W^{*}\to\mathbb{R}$ by making the parameter evolve under the Lie group action with $\alpha(0) = \alpha_0$, and, therefore, we can consider the extended Lagrangian function on $G^s\times \mathfrak{g}^s\times(T^*_{g_i}G)^s\times W^*$,

where $L_{\hbox{ext}}(g, u, \lambda, \alpha_0) = L(g, u, \lambda)$.

By assumptions (ⅰ) to (ⅳ), and by taking advantage of the $G$-invariance of $C_i$, $V_{i, 0}^{\text{ext}}$ and $V_{ij}$ (and so $V_i$), we can define the reduced extended Lagrangian $\ell:G^{s-1}\times \mathfrak{g}^s\times(\mathfrak{g}^*)^s\times W^*\to\mathbb{R}$ by

defining $g_1$ to be the identity $\bar{e}_1$ on $G$. This is equivalent to consider the Lagrangian $\tilde{L}_{\text{ext}}(h, u, \lambda, \alpha) = L_{\text{ext}}(\bar{e}_1, h, u, \lambda, \alpha)$ with $h\in G^{s-1}$. Then $\ell(g, u, \lambda, \alpha)$ is obtained by considering a left action in each term of the sum. In fact, if $L_{\text{ext}}^{i}$ is each term in the sum from $1$ to $s$ in the definition of $L_{\text{ext}}$, then

under the assumption that $g_{1} = \bar{e}_{1}$ and $\alpha_i = \rho_{g_i}^*(\alpha)$. Note here the slight abuse of notation regarding the positions $g$. In the definition of the reduced Lagrangian $g\in G^{s-1}$, while in that of the Lagrangian $L_{\text{ext}}$, $g\in G^s$.

Theorem 4.1. For $s\geq 2$, an extremal for the OCP (3.4) satisfies the following Euler-Poincaré equations

where $\mathit{\boldsymbol{J}}_W:T^{*}W\simeq W\times W^{*}\to \mathfrak{g}^{*}$ is the momentum map that corresponds to the left action of $G$ on $W$, and it is defined through the left representation $\rho$ of $G$ on $W$, and where $\Theta^{i}_{1} = 0$ if $i = 1$, otherwise $\Theta^{i}_{1} = 1$.

Proof. Consider the variational equation

which holds for variations of $g$, that vanish at the endpoints, and $u$. Also, consider the constrained variational equation

that holds for variations of $u_{i}$ and $\alpha_{i}$ with $\delta u_{i} = \dot{\eta}_{i}+ \text{ad}_{u_{i}}\eta_{i}$ and $\delta \alpha_{i} = \rho'^*_{\eta_{i}}(\alpha_{i})$, where $\eta_{i}$ is a path of $\mathfrak{g}^s$ that vanishes at the endpoints, i.e. $\eta_{i}(0) = \eta_{i} (T) = 0.$

The two variational equations are equivalent since the cost functions $C_i$ and the extended potential functions $V_{i, 0}^{\text{ext}}$ are $G$-invariant, i.e. $C_i\circ\Psi_g, \; \; V_{i, 0}^{\text{ext}}\circ\tilde{\Phi} = V_{i, 0}^{\text{ext}}$ and $\langle\lambda_{g_i}, T_{\bar{e}_i}L_{g_i}u_i\rangle = \langle T^{*}_{g_i}L_{g_i^{-1}}\lambda_{g_i}, u_i\rangle = \langle\lambda_i, u_i\rangle$. The variations $\delta g_i$ of $g_i$ induce and are induced by variations $\delta u_i = \dot{\eta}_i+ \text{ad}_{u_i}\eta_i$ with $\eta_i(0) = \eta_i(T) = 0$, where $\eta_i = T_{g_i}L_{g_i}^{-1}(\delta g_i)$ and $\delta u_i = T_{g_i}L_{g_i^{-1}}(\dot{g}_i)$. Variations of $\alpha_i$ are given by $\delta\alpha_i = \rho'^*_{\eta_i}(\alpha_i)$. Thus, we have

Using the variations of $u_i$, which are given by $\delta u_i = \dot{\eta}_i+ \text{ad}_{u_i}\eta_i$, applying integration by parts and by the definition of the co-adjoint action the first two terms, yield:

From the variations $\delta\alpha_i = \rho'^*_{\eta_i}(\alpha_i)$, the third term gives

Taking into account that $T(L_{g_i}\circ L_{g_i^{-1}}) = TL_{g_i}\circ TL_{g_i^{-1}}$ is equivalent to the identity map on $TG_i$ and $\eta_i = T_{g_i}L_{g_i^{-1}}(\delta g_i)$, the fourth term can be written as

Therefore, after performing a change of variables between indexes $i$ and $k$ in the fourth term, the above variational equation (4.4) yields

for $i = 1$. Otherwise,

Finally, by taking the time derivative of $\alpha_i = \rho_{g_i}(\alpha_0)$, we have $\dot{\alpha}_i = \rho_{u_i}^{\prime*}(\alpha_i)$, together with $\alpha(0) = \rho^{*}_{g_{0}^{i}}(\alpha_{0}^{i}).$

Note that the above Euler-Poincaré equations (4.2) cannot, in fact, describe the motion properly because there are more unknowns than equations. In particular, observe that equations (4.2) together with (3.1) (or equivalently (3.2)), give rise to only two equations for the three unknown variables $u_i$, $\lambda_i$ and $g_i$. However, we provide an additional structure to the Lie algebra, $\mathfrak{g}$, that allows one to decouple equations (4.2) into two equations. The next Proposition describes this process.

Proposition 4.2. If the structure of the Lie algebra permits a decomposition $\mathfrak{g} = \mathfrak{r}\oplus\mathfrak{s}$ where $\mathfrak{r} = {\rm{span}}\{e_1, \dots, e_m\}$ and $\mathfrak{s} = {\rm{span}}\{e_{m+1}, \dots, e_n\}$ such that

then the Euler-Poincaré equations of motion (4.2) are given by the following equations:

where $u_{i}\in \mathfrak{r}$ and the restrictions $\big|_{\mathfrak{r}}$ and $\big|_{\mathfrak{s}}$ give the projection onto $\mathfrak{r}^*$ and $\mathfrak{s}^*$, respectively, with respect to the associated splitting of the dual space $\mathfrak{g}^* = \mathfrak{r}^* \oplus \mathfrak{s}^*$.

Remark 4. Note that this is the case for semisimple Lie algebras, they admit a Cartan decomposition, i.e., if $\mathfrak{g}$ is semisimple, then $\mathfrak{g} = \mathfrak{r}\oplus\mathfrak{s}$ such that $[\mathfrak{r}, \mathfrak{r}]\subseteq\mathfrak{s}, \quad [\mathfrak{s}, \mathfrak{r}]\subseteq\mathfrak{r}, \quad [\mathfrak{s}, \mathfrak{s}]\subseteq\mathfrak{s}, \nonumber$ where the set $\mathfrak{r} = \{x\in\mathfrak{g}\mid\theta(x) = -x\}$ is the $-1$ eigenspace of the Cartan involution $\theta$, whereas the set $\mathfrak{s} = \{x\in\mathfrak{g}\mid\theta(x) = x\}$ is the $+1$ eigenspace of the Cartan involution $\theta$. Additionally, the Killing form is proven to be positive definite on $\mathfrak{r}$ and negative definite on $\mathfrak{s}$ (see, e.g., ^[15]). Thus, suitable candidates that satisfy the assumption of Proposition 4.2 are connected semisimple Lie groups. On the other hand, a Cartan decomposition determines a Cartan involution $\theta$ (see, e.g., ^[37]). In particular, the proposed decomposition for the Lie algebra is not restrictive in the sense that the usual manifolds/work-spaces used in applications as $SO(n)$ and $SE(n)$ allow such a decomposition.

Proof. Given $\mathfrak{g} = \mathfrak{r}\oplus\mathfrak{s}$ we get $\mathfrak{g}^* = \mathfrak{r}^*\oplus\mathfrak{s}^*$, where $\mathfrak{r}^* = {\rm{span}}\{e^1, \dots, e^m\}$ and $\mathfrak{s}^* = {\rm{span}}\{e^{m+1}, \dots, e^s\}$. Thus, from (4.6) we have that $\text{ad}^*_{\mathfrak{s}}\mathfrak{s}^*\subseteq\mathfrak{s}^*, \; \text{ad}^*_{\mathfrak{s}}\mathfrak{r}^*\subseteq\mathfrak{r}^*, \; \text{ad}^*_{\mathfrak{r}}\mathfrak{s}^*\subseteq\mathfrak{r}^*, \; \text{ad}^*_{\mathfrak{r}}\mathfrak{r}^*\subseteq\mathfrak{s}^*$ and given that $u_{i}\in \mathfrak{r}$, $\frac{\partial C_i}{\partial u_i}\in\mathfrak{r}^*$ and $\lambda_i\in\mathfrak{s}^*$ by definition we conclude that $\text{ad}^*_{u_i}\frac{\partial C_i}{\partial u_i}\in\mathfrak{s}^*$ and $\text{ad}^*_{u_i}\lambda_i\in\mathfrak{r}^*.$ Also, ${\bf{J}}_W\Big(\frac{\partial V_{i, 0}^{\text{ext}}}{\partial\alpha_i}, \alpha_i\Big) \in \mathfrak{g}^*$ and $\sum_{k = 1}^{s}T^*_{\overline{e}_i}L_{g_i}\bigg(\frac{\partial V_{k}}{\partial g_i}\bigg)\in \mathfrak{g}^*$ hence, they have a decomposition into $\mathfrak{r}^*$ and $\mathfrak{s}^*$. Thus, the equations (4.2) split into the following equations

Remark 5. For the initial value problem guaranteeing a solution for the previous system of equations, we look for a solution to the equations with the initial condition $u_i(0) = T_{g(0)}L_{g^{-1}(0)}(\dot{g}_i(0))$ and the kinematic equation $\dot{g}_i(t) = T_{\overline{e}}L_{g_i(t)}(u_i(t))$ with $g(0) = (g_1(0), \dots, g_s(0))$.

Remark 6. Assuming the reduced Lagrangian is hyperregular, we use the Legendre transformation and we can define the reduced Hamiltonian $h: G^{s-1}\times (\mathfrak{g}^*)^s\times (\mathfrak{g}^*)^s\times W^*\to \mathbb{R}$ given by

where $\mu = \frac{\partial \ell}{\partial u} = \left(\frac{\partial C}{\partial u}+\lambda\right)\in(\mathfrak{g}^*)^s$. Under this assumption, the Euler-Poincare equations (4.2), (4.3) can be written as the Lie-Poisson equations (see, e.g., ^[22])

5.   Discrete-time reduced necessary conditions

In this section, we study the discrete-time reduction by symmetries for necessary conditions in the collision and obstacle avoidance optimal control problem. The goal is to construct a variational integrator based on the discretization of the augmented cost functional. Such integrator inherits discrete-time symmetries from its continuous counterpart and generates a well-defined (local) flow for reduced necessary conditions characterizing (local) extremal in the optimal control problem.

5.1.   Trajectory discretization

Given the set $\mathcal{T} = \{t_k\in\mathbb{R}^{+}, \, t_{k} = kh\mid k = 0, \ldots, N\}$, $Nh = T$, with $T$ fixed (recall that $T\in\mathbb{R}^{+}$ is the end point of the cost functional - see for instance equation (3.4)), a discrete trajectory for the agent $i$ is determined by a set of $N+1$ points equally spaced in time, $g^{0:N}_i = \{g^0_i, \ldots, g^{N}_i\}$, where $g^k_i\simeq g_i(kh)\in G$, and $h = T/N$ is the time step. The path between two adjacent points $g_i^k$ and $g_i^{k+1}$ must be given by a curve lying on the Lie group $G$. To construct such a curve we make use of a retraction map $\mathcal{R}:\mathfrak{g}\to G$.

Definition 5.1. A retraction map $\mathcal{R}:\mathfrak{g}\to G$ is an analytic local diffeomorphism assigning a neighborhood $\mathcal{O}\subset\mathfrak{g}$ of $0\in {\mathfrak g}$ to a neighborhood of the identity $\overline{e}\in G$.

The retraction map (see Figure 1) is used to express small discrete changes in the group configuration through unique Lie algebra elements given by $u^k = \mathcal{R}^{-1}((g^{k})^{-1}g^{k+1})/h$, where $u^k\in\mathfrak{g}$ (see ^[38,39] for further details). That is, if $u^{k}$ were regarded as an average velocity between $g^{k}$ and $g^{k+1}$, then $\mathcal{R}$ is an approximation to the integral flow of the dynamics. The difference $\xi^{k, k+1}: = (g^{k})^{-1}\, g^{k+1}\in G$, which is an element of a nonlinear space, can now be represented by a vector space element $u^{k}$. For the derivation of the discrete equations of motion, the right trivialized tangent retraction map will be used. It is the function $d\mathcal{R}:\mathfrak{g}\times\mathfrak{g}\rightarrow \mathfrak{g}$ given by

where $\eta\in\mathfrak{g}$ and $R:G\to G$ the right translation on $G$ and $T\mathcal{R}(\xi)\cdot\eta$ is the directional derivative of $\mathcal{R}$ at $\xi$ in the direction of $\eta$ (see ^[38,39] for the derivation of such a map). Here we use the following notation, $d\mathcal{R}_\xi: = d\mathcal{R}(\xi):\mathfrak{g}\rightarrow \mathfrak{g}.$ The function $d\mathcal{R}$ is linear, but only on one argument.

Remark 7. The natural choice of a retraction map is the exponential map at the identity $\overline{e}$ of the group $G,$ $\mbox{exp}_{\overline{e}}:\mathfrak{g}\rightarrow G$. Bear in mind that, for a finite-dimensional Lie group, $\mbox{exp}_{\overline{e}}$ is locally a diffeomorphism and gives rise to a natural chart ^[40]. Then, there exists a neighborhood $U$ of $\overline{e}\in G$ such that $\mbox{exp}_{\overline{e}}^{-1}:U\rightarrow \mbox{exp}_{\overline{e}}^{-1}(U)$ is a local $\mathcal{C}^{\infty}-$diffeomorphism. For an element $g\in G$ a chart is given by $\Psi_{g} = \mbox{exp}_{\overline{e}}^{-1}\circ L_{g^{-1}}.$

Generally, it is not an easy task to work with the exponential map since the differential of the exponential map involves power series expansions with iterated Lie-brackets. Consequently, it will be more convenient to use a different retraction map. More concretely, the Cayley map, which is usually used in numerical integration with matrix Lie-groups configurations (see ^[38,39] for further details), will provide to us a proper framework in the application shown in the next Section.

5.2.   Discretization of the Lagrangian function

Next, we consider a discrete cost function to construct variational integrators in the same way as in discrete mechanics ^[41]. In other words, consider the continuous-time Lagrangian ${\bf{L}}:G^s\times\mathfrak{g}^s\to\mathbb{R}$ defined by the cost functional (3.4), that is,

and for a given $h > 0$, we define the discrete Lagrangian $L_d:G^s\times \mathfrak{g}^s\to\mathbb{R}$ as an approximation of the cost functional (3.4) along each discrete segment between $g^{k}$ and $g^{k+1} = g^{k}\mathcal{R}(hu^{k})$, that is,

where $\kappa$ and $\zeta$ are functions of $(g^k, u^{k})\in G^s\times \mathfrak{g}^s$ which approximate the configuration $g(t)$ and the control input $u(t)$, respectively, in that interval. In the following, for simplicity, $\kappa$ will be the projection onto $G^{s}$ and $\zeta$ will be the projection onto $\mathfrak{g}^{s}$. Thus, we consider

We remark that different choices of $\kappa$ and $\zeta$ could result in higher-order numerical methods, such as, using the middle point rule with $\kappa(g^k, u^{k}) = g^{k}\mathcal{R}(\frac{h}{2}u^{k})$.

5.3.   Discrete-time optimal control problem and reduction of discrete-time necessary conditions for optimality

Next, we are going to define the optimal control problem for discrete-time systems and derive a variational integrator for $L_{d}\colon G^s\times \mathfrak{g}^s\to \mathbb{R}$, in a similar fashion as the variational equation presented in Theorem 4.1.

Problem: Consider the discrete-time optimal control problem for collision and obstacle avoidance of left-invariant multi agent control systems, which is given by finding the discrete configurations $\{g^k\}_{k = 0}^{N} = \{(g_1^k, \ldots, g_s^k)\}_{k = 0}^{N}$ and discrete control inputs $\{u^k\}_{k = 0}^{N} = \{(u_1^k, \ldots, u_s^k)\}_{k = 0}^{N}$ minimizing the discrete cost functional

subject to $g_i^{k+1} = g_i^k\mathcal{R}(hu_i^k)$ (i.e., a first order approximation of equation (3.1)) with given boundary conditions $g^0$ and $g^N$, where $h > 0$ denotes the time step, $\mathcal{R}:\mathfrak{g}\to G$ is a retraction map, $u_i(0)$ and $u_i(T)$ are given, and each cost function $C_{i}:G\times\mathfrak{g}\to\mathbb{R}$, potential functions $V_{i}^{0}$ and $V_{i}$ satisfy properties (ⅰ) - (ⅶ).

The discrete-time optimal control problem (5.3) can be considered as a discrete constrained variational problem by introducing the Lagrange multipliers $\mu_{i}^k\in \mathfrak{g}$ into the cost functional. Consider the augmented discrete Lagrangian $\mathcal{L}_{d}:G^s\times G^s\times \mathfrak{g}^s\times (\mathfrak{g}^{*})^s\to\mathbb{R}$ given by

where $\xi^{k, k+1}_i = (g_i^k)^{-1}g_i^{k+1}\in G$, for each $i\in\mathcal{N}$. Note that the last term in the augmented Lagrangian represents a first-order discretization of the kinematic constraint paired with a Lagrange multiplier in analogy with the variational equation presented in Section 4.

Now, extending the potential $V_{i}^{0}$, we obtain an extended Lagrangian $\mathcal{L}_{ext, d}:G^s\times G^s\times\mathfrak{g}^s\times(\mathfrak{g}^{*})^s\times (W^{*})^{s}\to\mathbb{R}$ given by

which is invariant under the left action of $G$ on $G^s\times G^s\times \mathfrak{g}^s \times (\mathfrak{g}^{*})^s \times (W^{*})^s$ given by $\tilde{\Phi}_{g}(h_{1}, h_{2}, u, \mu, \alpha) = (gh_{1}, gh_{2}, u, \mu, \rho_{g^{-1}}^{*}(\alpha))$ by assumption (ⅵ). In particular, under assumptions (ⅳ)-(ⅵ), the extended discrete Lagrangian $\mathcal{L}_{ext, d}(\cdot, \cdot, \cdot, \cdot, \alpha_0) = :\mathcal{L}_{ext, d, \alpha_0}$ is $G_{\alpha_0}$-invariant under $\tilde{\Phi}_g$.

The following result (Theorem 5.2) derives a variational integration for reduced optimality conditions for the discrete-time optimal control (5.3) in analogy with the results presented in Section 4. To derive the numerical algorithm, first we need the following result describing variations for elements on the Lie algebra and its relation with variations on the Lie group by using the retraction map, in addition to a property used in the proof of Theorem 5.2. Though it is a well-known result in the literature, we include it here to make the paper self-contained.

Lemma 5.1 (adapted from ^[38,39]). The following properties hold

(i)

where $\eta^{k} = T_{g^k}L_{(g^k)^{-1}}(\delta g^k)\in\mathfrak{g}$ and $\xi^{k, k+1} = (g^k)^{-1}g^{k+1}$.

(ii)

where $\mu^k\in(\mathfrak{g}^{*})$ and $d\mathcal{R}^{-1}$ is the inverse right trivialized tangent of the retraction map $\mathcal{R}$ defined in (5.1).

Theorem 5.2. Under assumptions (i)-(vii), an extremal for the discrete-time optimal control problem (5.3) satisfies the following equations

for $k = 1, \ldots, N-1$; where $\Theta^{i}_{1} = 0$ if $i = 1$, otherwise $\Theta^{i}_{1} = 1$.

Proof. Since the cost functions and the potential functions satisfy assumptions (ⅰ) - (ⅶ), as in the continuous-time case, it is possible to induce the reduced augmented discrete Lagrangian $\ell_{ext, d}:G^{s-1}\times G^s \times\mathfrak{g}^s\times(\mathfrak{g}^{*})^s\times (W^{*})^s\to\mathbb{R}$ as

where $\bar\alpha^{k}_i = \rho_{g^{k}_i}^{*}(\alpha_i^0)$ for a fixed $\alpha_i^0\in W^{*}$ satisfying $\alpha_i^0 = \rho^*_{g_i^0}(\alpha_i^0)$ and, with a slight abuse of notation, $C_{i}(u_{i}^k) = C_{i}(\overline{e}, u_{i}^k)$ and $V_{i}^{0, ext}(\alpha^{k}_i) = V_{i}^{0, ext}(\overline{e}, \alpha^{k}_i)$. Notice that, also here, $g_{1}^{k}$ is set to be the identity element, so that we have $g^{k}\in G^{s-1}$.

As in the proof for Theorem 4.1, the technical part is to show that an extremal of the reduced variational equation

satisfies equations (5.6)-(5.8) for all variations of $\mathcal{R}^{-1}(\xi^{k, k+1})$ (induced by variations of $g^k$ vanishing at the endpoints), $u^k$ and $\bar{\alpha}^k$ of the form $\rho_{\eta^{k}}^{'*}(\bar\alpha^{k})$, where $\eta^k\in\mathfrak{g}^s$ vanishes at the endpoints. Then, similarly as in the proof for Theorem 4.1, it follows that an extremal for the optimal control problem (5.3) satisfies the variational equation

for all variations of $g^k$ (vanishing at the endpoints), $\mathcal{R}^{-1}(\xi^{k, k+1})$ (induced by variations of $g^k$) and $u^k$.

Note that

where we used Lemma 5.1 to obtain the last term. Since variations $\delta u_{i}^{k}$ are arbitrary, we obtain $\mu_{i}^{k} = \frac{\partial C_{i}}{\partial u_{i}^{k}}$. As for the second term, we have that

As we have show along the proof for Theorem 4.1, we have that

and the indexes might be interchanged. The last term to obtain equation (5.7) may be dealt with, using integration by parts in discrete-time, which is just rearranging the indexes, together with the second statement in Lemma 5.1 and the fact that $\eta_{i}^{0} = \eta_{i}^{N} = 0$. Therefore,

and the equations (5.6)-(5.8) follow.

Remark 8. Equations (5.6)-(5.8) are as a discrete approximation of the Lie-Poisson equations for the Hamiltonian version of the optimal control problem considered in ^[22]. The equation ${\mu_i^k = \left(\frac{\partial C_i}{\partial u^k_i}\right)}$ represents the discrete time version of the reduced Legendre transformation and the equation $g^{k+1}_i = g^k_i\mathcal{R}(hu^k_i)$ is the analogous of the reconstruction equation in the discrete time counterpart. These three equations are used to compute $u_i^{k}, \mu_i^{k}$ and $g^{k+1}_i$ given $u_i^{k-1}$, $\mu_i^{k-1}$, $g^{k-1}_i$ and $g^{k}_i$ from $k = 1$ to $k = N-1$.

To compute the discrete-time reduced necessary condition for the optimal control problem (5.3) we must enforce boundary conditions given by the continuous-time quantities. More precisely, we must set

relating the momenta at the initial and final times, and use to transform boundary values between the continuous and discrete representation. They follow from the principle that any variation with free boundary points of the action (5.10) along a solution of equations (5.6)-(5.9) equals the change in momentum $\langle \frac{\partial C_i}{\partial u_i}(u_i(T)), g_{i}^N \delta g_{i}^{N} \rangle - \langle \frac{\partial C_i}{\partial u_i}(u_i(0)), g_{i}^0 \delta g_{i}^{0} \rangle$ (see ^[39] for a discussion in the single agent case). Therefore, we must enforce

where the variation is taken along a family of solutions of the discrete equations (5.6)-(5.9) whose boundary values are not fixed. By taking the variations of the discrete action, we eventually get a vanishing term corresponding to the fact that this family of sequences satisfies the discrete equations together with boundary terms multiplying $\eta_{i}^{0}$ and $\eta_{i}^{N}$, and simplifying, to the above equations.

Remark 9. If we choose the midpoint rule to discretize the potential $V_{i}$, then we would obtain the following boundary conditions

The boundary condition $g_s(T)$ for agent $s$ is enforced by the relation

Recalling that $\mathcal{R}(0) = \overline{e}$, this last expression just means that $g^N_s = g_s(T)$. Moreover, by computing recursively the equation $g^{k+1}_s = g^{k}_s\mathcal{R}(hu_s^k)$ for $k = 1, \ldots, N-1$, using that $g_s^{0} = g_s(0)$ and $\mathcal{R}(0) = \overline{e}$, it is possible to translate the final configuration $g_s^N$ in terms of $u^k_s$, such that there is no need to optimize over any of the configurations $g^k_s$. In that sense, (5.7) for $i = s$, $l = 0$ together with

form a set of $(nN)$-equations (since $\dim\mathfrak{g} = n$) where $nN$ unknowns are for $u^{0:N-1}_s$, denoting the entire sequence of controls $\{u_{s}^{0}, \dots, u_{s}^{N-1}\}$ for each agent $s$.

The numerical algorithm to compute the reduced optimality conditions is summarized in Algorithm 1.

Note also that the exact form of equations (5.6)-(5.8) depends on the choice of $\mathcal{R}$. This choice will also effect the computational efficiency of the optimization framework in the case, the above equalities are imposed as constraints. For instance, in Section 6, we will employ the Cayley transform on the Lie group $SE(2)$ as a choice of $\mathcal{R}$ to write in a compact form the numerical integrator ^[38], ^[17], but another natural choice would be to employ the exponential map, as we explained in Section 5.1.

6.   Case study

In this case study, we apply the proposed reduction by symmetry strategy to an optimal control for autonomous surface vehicles (ASVs). The configuration space whose elements determine the motion of each ASV is $SE(2)\cong SO(2)\times \mathbb{R}^2$. An element $g_i\in SE(2)$ is given by $g_i = \begin{pmatrix} \cos\theta_i & -\sin\theta_i&x_i\\ \sin\theta_i & \cos\theta_i&y_i\\ 0 & 0 & 1 \end{pmatrix}$, where $(x_i, y_i)\in \mathbb{R}^2$ represents the center of mass of a planar rigid body describing the ASV and $\theta_i$ represents the angular orientation of the ASV. The control inputs, for each ASV, are given by $u_i = (u_i^1, u_i^2)$, where $u_i^1$ denotes the speed of the center of mass for the ASV and $u_i^2$ denotes the angular velocity of the ASV.

The kinematic equations for the multi-agent system are:

Using the notation of Example 1, the Lie algebra $\mathfrak{se}(2)$ is identified with $\mathbb{R}^2$ through the isomorphism $\begin{pmatrix} -a\mathbb{J}&b\\ 0 & 0 \end{pmatrix}\mapsto (a, b).$ The elements of the basis of the Lie algebra $\mathfrak{se}(2)$ are ${e_1 = \begin{pmatrix} 0 & -1 & 0\\ 1 & 0 & 0\\ 0 & 0 & 0 \end{pmatrix}, e_2 = \begin{pmatrix} 0 & 0 & 1\\ 0 & 0 & 0\\ 0 & 0 & 0 \end{pmatrix}, e_3 = \begin{pmatrix} 0 & 0 & 0\\ 0 & 0 & 1\\ 0 & 0 & 0 \end{pmatrix}}$,

which satisfy $[e_1, e_2] = e_3, \; [e_2, e_3] = 0, \; [e_3, e_1] = e_2$. Thus, the kinematic equations (6.1) take the form $\dot{g_i} = g_iu_i = g_i(u_i^1e_1+u_i^2e_2)$ and give rise to a left-invariant control system on $SE(2)^s\times \mathfrak{se}(2)^s$. The inner product on $\mathfrak{se}(2)$ is given by $\langle\langle\xi_1, \xi_2\rangle\rangle: = tr(\xi_1^T\xi_2)$ for $\xi_1, \xi_2\in \mathfrak{se}(2)$ and hence, the norm is given by $||\xi|| = \sqrt{tr(\xi^T\xi)},$ for any $\xi\in \mathfrak{se}(2)$. The dual Lie algebra $\mathfrak{se}(2)^*$ of $SE(2)$ is defined through the dual pairing, $\langle\alpha, \xi\rangle = tr(\alpha \xi)$, where $\alpha\in \mathfrak{se}(2)^*$ and $\xi\in \mathfrak{se}(2)$, hence, the elements of the basis of $\mathfrak{se}(2)^*$ are ${e^1 = \begin{pmatrix} 0 & \frac{1}{2} & 0\\ -\frac{1}{2} & 0 & 0\\ 0 & 0 & 0 \end{pmatrix}, e^2 = \begin{pmatrix} 0 & 0 & 0\\ 0 & 0 & 0\\ 1 & 0 & 0 \end{pmatrix}, e^3 = \begin{pmatrix} 0 & 0 & 0\\ 0 & 0 & 0\\ 0 & 1 & 0 \end{pmatrix}}$.

Consider the cost function $C_i(g_i, u_i) = \frac{1}{2}\langle u_i, u_i\rangle$ and the artificial potential function $V_{ij}:SE(2)\times SE(2)\to \mathbb{R}$ given by ${V_{ij}(g_i, g_j) = \frac{\sigma_{ij}}{2((x_i-x_j)^2+(y_i-y_j)^2-4\overline{r}^2)}}$, where $\sigma_{ij}\in \mathbb{R}_{\geq 0}$ and $\overline{r}$ is the radius of the disk each agent occupies as defined at the end of Section Ⅲ. Consider a spherical obstacle with unit radius and without loss of generality let it be centered at the origin. Hence, consider the obstacle avoidance potential function $V_i^0:SE(2)\to \mathbb{R}$, ${V_i^0(g_i) = \frac{\sigma_{i0}}{2(x^2+y^2-(\overline{r}+1)^2)}}$, where $\sigma_{i0}\in\mathbb{R}_{ > 0}$.

Note that the obstacle avoidance potential functions are not $SE(2)$-invariant but $SO(2)$-invariant, so they break the symmetry. Using the norm of $\mathfrak{se}(2)$ and for $\overline{r} = 1$, $V_{ij}$ and $V_i^0$ are equivalently given by ${V_{ij}(g_i, g_j) = \frac{\sigma_{ij}}{2(|| \text{Ad}_{g_i^{-1}g_j}e_1||^2-6)}}$ and ${V_i^0(g_i) = \frac{\sigma_{i0}}{2(|| \text{Ad}_{g^{-1}_i}e_1||^2-6)}}$.

Let $W = \mathfrak{se}(2)^*$, so we define the extended potential functions $V_{i, 0}^{\text{ext}}:SE(2)\times \mathfrak{se}(2)\to \mathbb{R}$ by $V_{i, 0}^{\text{ext}}(g_i, \alpha) = \frac{\sigma_{i0}}{2(|| \text{Ad}_{g^{-1}_i}\alpha||^2-6)}$, which are $SE(2)$-invariant under the action of $\tilde{\Phi}$ given by (3.6), i.e. $V_{i, 0}^{\text{ext}}\circ\tilde{\Phi} = V_{i, 0}^{\text{ext}},$ for any $g\in SE(2)$. Since, $W = \mathfrak{se}(2)^*$ we have ${{\bf{J}}_W\Big(\frac{\partial V_{i, 0}^{\text{ext}}}{\partial\alpha_i}, \alpha_i\Big) = \text{ad}^*_{\alpha_i}\Big(\frac{\partial V_{i, 0}^{\text{ext}}}{\partial\alpha_i}\Big)}$, and equations (4.2) and (4.3) yield

together with $\dot{\alpha}_i = - \text{ad}_{u_i}\alpha_i$ and $\alpha_i = \text{Ad}_{g_i^{-1}}\alpha_0$.

Note also that $u_i = u_i^1e_1+u_i^2e_2$, $\alpha_i = \alpha_i^1e_1+\alpha_i^2e_2+\alpha_i^3e_3$ and $\lambda_i = \lambda_3^ie^3$ thus,

where $x_{ij} = x_i-x_j, \; y_{ij} = y_i-y_j$ and ${ \text{Ad}_{g_i^{-1}}\alpha_0 = \begin{pmatrix} 0 & -1&x_i\sin\theta_i-y_i\cos\theta_i\\ 1 & 0&x_i\cos\theta_i+y_i\sin\theta_i\\ 0 & 0 & 0 \end{pmatrix}}$.

Therefore, by applying Proposition 4.2 for $\mathfrak{r} = \text{span} \{e_1, e_2\}$ and $\mathfrak{s} = \text{span} \{e_3\}$, Euler-Lagrange equations for the OCP (3.4) are

with

$\begin{array}{ll} \dot{\alpha}_i^1 = 0, & \alpha_i^1(0) = 1, \\ \dot{\alpha}_i^2 = u_i^1\alpha_i^3, & \alpha_i^2(0) = x_i^0\sin\theta_i^0-y_i^0\cos\theta_i^0, \\ \dot{\alpha}_i^3 = -u_i^1\alpha_i^2+u_i^2\alpha_i^1, & \alpha_i^3(0) = x_i^0\cos\theta_i^0+y_i^0\sin\theta_i^0. \end{array}$

For the discrete-time setting, one would choose

where ${g^{k}_i = \begin{bmatrix} \cos\theta^{k}_i & -\sin\theta^{k}_i & \phantom{-}x^{k}_i\\ \sin\theta^{k}_i & \phantom{-}\cos\theta^{k}_i & \phantom{-}y^{k}_i\\ 0 & \phantom{-}0 & \phantom{-}1 \end{bmatrix}\in SE(2)\nonumber}$ and ${u^{k}_i = \sum\limits_{j = 1}^{3}u^{k}_{ij}e_{j}\in\mathfrak{se}(2)}$. Also, in the discrete-time setting, the extended potential function $V_{d, \text{ext}}: SE(2)\times \mathfrak{se}(2)\to\mathbb{R}$ can be constructed in exactly the same way as in the above example, and is given by

where ${\alpha_{i}^{k} = \sum\limits_{j = 1}^{3} \alpha_{ij}^{k}e^{j}}$. We do not give all the details again and leave it up to reader to verify that the assumptions (ⅰ) - (ⅶ) from (3.3) are satisfied. The discrete-time equations are

where

For numerical purposes, we first choose a suitable retraction map, like the Cayley map or the exponential map, and then compute the quantities $(d\mathcal{R}_{hu^{k}_i}^{-1})^{*}\mu^{k}_i$ and $(d\mathcal{R}_{-hu^{k-1}_i}^{-1})^{*}\mu^{k-1}_i$. As an example, if we choose the Cayley map $\text{cay}: \mathfrak{se}(2)\to\mathrm{SE}(2)$ as the retraction map, (see ^[38] and ^[17] for instance) then we have

where

and ${\mu_{i}^{k} = \sum\limits_{j = 1}^{3} \mu_{ij}^{k}e^{j}}$. Note that for ${v = \sum\limits_{i = 1}^{3}v^{i}e_{i}\in\mathfrak{g}}$, the matrix representation for $d\text{cay}_{v}^{-1}$ is given by

7.   Conclusions

We studied the reduction by symmetry for optimality conditions of extremal in an OCP for collision and obstacle avoidance of left-invariant multi-agent control system on Lie groups, by exploiting the physical symmetries of the agents and obstacles. Reduced optimality conditions are obtained using techniques from variational calculus and Lagrangian mechanics on Lie groups, in the continuous-time and discrete-time settings. We applied the results to an OCP for multiple unmanned surface vehicles. The method proposed in this work allows the construction of position and velocity estimators, by discretizing the variational equation given in Theorem 4.1 - instead of discretizing the equations of motion - and by deriving variational integrators - see Theorem 5.2. The reduction of sufficient conditions for optimality will be also studied by using the notion of conjugate points, as in ^[42], in future work, as well as the reduction by symmetry of the variational obstacle avoidance problems ^[43] on semidirect products of Lie groups endowed with a bi-invariant metric on a Riemannian manifold.

Acknowledgments

L. Colombo is very grateful to A. Bloch, R. Gupta and T. Ohsawa for many useful comments and stimulating discussions during the last years on the topic of this paper, which is inspired by our common previous work.

The authors acknowledge financial support from the Spanish Ministry of Science and Innovation, under grants PID2019-106715GB-C21, MTM2016-76702-P.

Conflict of interest

We the authors declare that there are no conflicts of interest.