9,526 research outputs found
A variational framework for flow optimization using semi-norm constraints
When considering a general system of equations describing the space-time
evolution (flow) of one or several variables, the problem of the optimization
over a finite period of time of a measure of the state variable at the final
time is a problem of great interest in many fields. Methods already exist in
order to solve this kind of optimization problem, but sometimes fail when the
constraint bounding the state vector at the initial time is not a norm, meaning
that some part of the state vector remains unbounded and might cause the
optimization procedure to diverge. In order to regularize this problem, we
propose a general method which extends the existing optimization framework in a
self-consistent manner. We first derive this framework extension, and then
apply it to a problem of interest. Our demonstration problem considers the
transient stability properties of a one-dimensional (in space) averaged
turbulent model with a space- and time-dependent model "turbulent viscosity".
We believe this work has a lot of potential applications in the fluid
dynamics domain for problems in which we want to control the influence of
separate components of the state vector in the optimization process.Comment: 30 page
Quantitative Stability of Linear Infinite Inequality Systems under Block Perturbations with Applications to Convex Systems
The original motivation for this paper was to provide an efficient
quantitative analysis of convex infinite (or semi-infinite) inequality systems
whose decision variables run over general infinite-dimensional (resp.
finite-dimensional) Banach spaces and that are indexed by an arbitrary fixed
set . Parameter perturbations on the right-hand side of the inequalities are
required to be merely bounded, and thus the natural parameter space is
. Our basic strategy consists of linearizing the parameterized
convex system via splitting convex inequalities into linear ones by using the
Fenchel-Legendre conjugate. This approach yields that arbitrary bounded
right-hand side perturbations of the convex system turn on constant-by-blocks
perturbations in the linearized system. Based on advanced variational analysis,
we derive a precise formula for computing the exact Lipschitzian bound of the
feasible solution map of block-perturbed linear systems, which involves only
the system's data, and then show that this exact bound agrees with the
coderivative norm of the aforementioned mapping. In this way we extend to the
convex setting the results of [3] developed for arbitrary perturbations with no
block structure in the linear framework under the boundedness assumption on the
system's coefficients. The latter boundedness assumption is removed in this
paper when the decision space is reflexive. The last section provides the aimed
application to the convex case
On the stability of the optimal value and the optimal set in optimization problems
The paper develops a stability theory for the optimal value and the optimal set mapping of optimization problems posed in a Banach space. The problems considered in this paper have an arbitrary number of inequality constraints involving lower semicontinuous (not necessarily convex) functions and one closed abstract constraint set. The considered perturbations lead to problems of the same type as the nominal one (with the same space of variables and the same number of constraints), where the abstract constraint set can also be perturbed. The spaces of functions involved in the problems (objective and constraints) are equipped with the metric of the uniform convergence on the bounded sets, meanwhile in the space of closed sets we consider, coherently, the Attouch-Wets topology. The paper examines, in a unified way, the lower and upper semicontinuity of the optimal value function, and the closedness, lower and upper semicontinuity (in the sense of Berge) of the optimal set mapping. This paper can be seen as a second part of the stability theory presented in [17], where we studied the stability of the feasible set mapping (completed here with the analysis of the Lipschitz-like property)
Semi-infinite optimization: Structure and stability of the feasible set
The problem of the minimization of a functionf: ℝn→ℝ under finitely many equality constraints and perhaps infinitely many inequality constraints gives rise to a structural analysis of the feasible setM[H, G]={x∈ℝn¦H(x)=0,G(x, y)≥0,y∈Y} with compactY⊂ℝr. An extension of the well-known Mangasarian-Fromovitz constraint qualification (EMFCQ) is introduced. The main result for compactM[H, G] is the equivalence of the topological stability of the feasible setM[H, G] and the validity of EMFCQ. As a byproduct, we obtain under EMFCQ that the feasible set admits local linearizations and also thatM[H, G] depends continuously on the pair (H, G). Moreover, EMFCQ is shown to be satisfied generically
Sum-of-squares of polynomials approach to nonlinear stability of fluid flows: an example of application
With the goal of providing the first example of application of a recently proposed method, thus demonstrating its ability to give results in principle, global stability of a version of the rotating Couette flow is examined. The flow depends on the Reynolds number and a parameter characterising the magnitude of the Coriolis force. By converting the original Navier-Stokes equations to a finite-dimensional uncertain dynamical system using a partial Galerkin expansion, high-degree polynomial Lyapunov functionals were found by sum-of-squares-of-polynomials optimization. It is demonstrated that the proposed method allows obtaining the exact global stability limit for this flow in a range of values of the parameter characterising the Coriolis force. Outside this range a lower bound for the global stability limit was obtained, which is still better than the energy stability limit. In the course of the study several results meaningful in the context of the method used were also obtained. Overall, the results obtained demonstrate the applicability of the recently proposed approach to global stability of the fluid flows. To the best of our knowledge, it is the first case in which global stability of a fluid flow has been proved by a generic method for the value of a Reynolds number greater than that which could be achieved with the energy stability approach
Low Complexity Regularization of Linear Inverse Problems
Inverse problems and regularization theory is a central theme in contemporary
signal processing, where the goal is to reconstruct an unknown signal from
partial indirect, and possibly noisy, measurements of it. A now standard method
for recovering the unknown signal is to solve a convex optimization problem
that enforces some prior knowledge about its structure. This has proved
efficient in many problems routinely encountered in imaging sciences,
statistics and machine learning. This chapter delivers a review of recent
advances in the field where the regularization prior promotes solutions
conforming to some notion of simplicity/low-complexity. These priors encompass
as popular examples sparsity and group sparsity (to capture the compressibility
of natural signals and images), total variation and analysis sparsity (to
promote piecewise regularity), and low-rank (as natural extension of sparsity
to matrix-valued data). Our aim is to provide a unified treatment of all these
regularizations under a single umbrella, namely the theory of partial
smoothness. This framework is very general and accommodates all low-complexity
regularizers just mentioned, as well as many others. Partial smoothness turns
out to be the canonical way to encode low-dimensional models that can be linear
spaces or more general smooth manifolds. This review is intended to serve as a
one stop shop toward the understanding of the theoretical properties of the
so-regularized solutions. It covers a large spectrum including: (i) recovery
guarantees and stability to noise, both in terms of -stability and
model (manifold) identification; (ii) sensitivity analysis to perturbations of
the parameters involved (in particular the observations), with applications to
unbiased risk estimation ; (iii) convergence properties of the forward-backward
proximal splitting scheme, that is particularly well suited to solve the
corresponding large-scale regularized optimization problem
- …