1,079 research outputs found
Sparse Stabilization and Control of Alignment Models
From a mathematical point of view self-organization can be described as
patterns to which certain dynamical systems modeling social dynamics tend
spontaneously to be attracted. In this paper we explore situations beyond
self-organization, in particular how to externally control such dynamical
systems in order to eventually enforce pattern formation also in those
situations where this wished phenomenon does not result from spontaneous
convergence. Our focus is on dynamical systems of Cucker-Smale type, modeling
consensus emergence, and we question the existence of stabilization and optimal
control strategies which require the minimal amount of external intervention
for nevertheless inducing consensus in a group of interacting agents. We
provide a variational criterion to explicitly design feedback controls that are
componentwise sparse, i.e. with at most one nonzero component at every instant
of time. Controls sharing this sparsity feature are very realistic and
convenient for practical issues. Moreover, the maximally sparse ones are
instantaneously optimal in terms of the decay rate of a suitably designed
Lyapunov functional, measuring the distance from consensus. As a consequence we
provide a mathematical justification to the general principle according to
which "sparse is better" in the sense that a policy maker, who is not allowed
to predict future developments, should always consider more favorable to
intervene with stronger action on the fewest possible instantaneous optimal
leaders rather than trying to control more agents with minor strength in order
to achieve group consensus. We then establish local and global sparse
controllability properties to consensus and, finally, we analyze the sparsity
of solutions of the finite time optimal control problem where the minimization
criterion is a combination of the distance from consensus and of the l1-norm of
the control.Comment: 33 pages, 5 figure
Approximation of high-dimensional parametric PDEs
Parametrized families of PDEs arise in various contexts such as inverse
problems, control and optimization, risk assessment, and uncertainty
quantification. In most of these applications, the number of parameters is
large or perhaps even infinite. Thus, the development of numerical methods for
these parametric problems is faced with the possible curse of dimensionality.
This article is directed at (i) identifying and understanding which properties
of parametric equations allow one to avoid this curse and (ii) developing and
analyzing effective numerical methodd which fully exploit these properties and,
in turn, are immune to the growth in dimensionality. The first part of this
article studies the smoothness and approximability of the solution map, that
is, the map where is the parameter value and is the
corresponding solution to the PDE. It is shown that for many relevant
parametric PDEs, the parametric smoothness of this map is typically holomorphic
and also highly anisotropic in that the relevant parameters are of widely
varying importance in describing the solution. These two properties are then
exploited to establish convergence rates of -term approximations to the
solution map for which each term is separable in the parametric and physical
variables. These results reveal that, at least on a theoretical level, the
solution map can be well approximated by discretizations of moderate
complexity, thereby showing how the curse of dimensionality is broken. This
theoretical analysis is carried out through concepts of approximation theory
such as best -term approximation, sparsity, and -widths. These notions
determine a priori the best possible performance of numerical methods and thus
serve as a benchmark for concrete algorithms. The second part of this article
turns to the development of numerical algorithms based on the theoretically
established sparse separable approximations. The numerical methods studied fall
into two general categories. The first uses polynomial expansions in terms of
the parameters to approximate the solution map. The second one searches for
suitable low dimensional spaces for simultaneously approximating all members of
the parametric family. The numerical implementation of these approaches is
carried out through adaptive and greedy algorithms. An a priori analysis of the
performance of these algorithms establishes how well they meet the theoretical
benchmarks
Numerical controllability of the wave equation through primal methods and Carleman estimates
This paper deals with the numerical computation of boundary null controls for
the 1D wave equation with a potential. The goal is to compute an approximation
of controls that drive the solution from a prescribed initial state to zero at
a large enough controllability time. We do not use in this work duality
arguments but explore instead a direct approach in the framework of global
Carleman estimates. More precisely, we consider the control that minimizes over
the class of admissible null controls a functional involving weighted integrals
of the state and of the control. The optimality conditions show that both the
optimal control and the associated state are expressed in terms of a new
variable, the solution of a fourth-order elliptic problem defined in the
space-time domain. We first prove that, for some specific weights determined by
the global Carleman inequalities for the wave equation, this problem is
well-posed. Then, in the framework of the finite element method, we introduce a
family of finite-dimensional approximate control problems and we prove a strong
convergence result. Numerical experiments confirm the analysis. We complete our
study with several comments
A multidomain spectral method for solving elliptic equations
We present a new solver for coupled nonlinear elliptic partial differential
equations (PDEs). The solver is based on pseudo-spectral collocation with
domain decomposition and can handle one- to three-dimensional problems. It has
three distinct features. First, the combined problem of solving the PDE,
satisfying the boundary conditions, and matching between different subdomains
is cast into one set of equations readily accessible to standard linear and
nonlinear solvers. Second, touching as well as overlapping subdomains are
supported; both rectangular blocks with Chebyshev basis functions as well as
spherical shells with an expansion in spherical harmonics are implemented.
Third, the code is very flexible: The domain decomposition as well as the
distribution of collocation points in each domain can be chosen at run time,
and the solver is easily adaptable to new PDEs. The code has been used to solve
the equations of the initial value problem of general relativity and should be
useful in many other problems. We compare the new method to finite difference
codes and find it superior in both runtime and accuracy, at least for the
smooth problems considered here.Comment: 31 pages, 8 figure
Spectral cluster bounds for orthonormal systems and oscillatory integral operators in Schatten spaces
We generalize the spectral cluster bounds of Sogge for the
Laplace-Beltrami operator on compact Riemannian manifolds to systems of
orthonormal functions. The optimality of these new bounds is also discussed.
These spectral cluster bounds follow from Schatten-type bounds on oscillatory
integral operators.Comment: 30 page
- …