Nash Equilibria for the Multiobjective Control of Linear Partial Differential Equations

This article is concerned with the numerical solution of multiobjective control problems associated with linear partial differential equations. More precisely, for such problems, we look for the Nash equilibrium, which is the solution to a noncooperative game. First, we study the continuous case. Then, to compute the solution of the problem, we combine finite-difference methods for the time discretization, finite-element methods for the space discretization, and conjugate-gradient algorithms for the iterative solution of the discrete control problems. Finally, we apply the above methodology to the solution of several tests problems.     

On exact and approximate boundary controllabilities for the heat equation: A numerical approach

The present article is concerned with the numerical implementation of the Hilbert uniqueness method for solving exact and approximate boundary controllability problems for the heat equation. Using convex duality, we reduce the solution of the boundary control problems to the solution of identification problems for the initial data of an adjoint heat equation. To solve these identification problems, we use a combination of finite difference methods for the time discretization, finite element methods for the space discretization, and of conjugate gradient and operator splitting methods for the iterative solution of the discrete control problems. We apply then the above methodology to the solution of exact and approximate boundary controllability test problems in two space dimensions. The numerical results validate the methods discussed in this article and clearly show the computational advantage of using second-order accurate time discretization methods to approximate the control problems.     

Some iterative methods for the solution of a symmetric indefinite KKT system

This paper is concerned with the numerical solution of a Karush–Kuhn–Tucker system. Such symmetric indefinite system arises when we solve a nonlinear programming problem by an Interior-Point (IP) approach. In this framework, we discuss the effectiveness of two inner iterative solvers: the method of multipliers and the preconditioned conjugate gradient method. We discuss the implementation details of these algorithms in an IP scheme and we report the results of a numerical comparison on a set of large scale test-problems arising from the discretization of elliptic control problems. This research was supported by the Italian Ministry for Education, University and Research (MIUR), FIRB Project RBAU01JYPN.     

Proper efficiency and duality for a class of constrained multiobjective fractional optimal control problems containing arbitrary norms

We establish necessary and sufficient conditions for properly efficient solutions of a class of nonsmooth nonconvex optimal control problems with multiple fractional objective functions, linear dynamics, and nonlinear inequality constraints on both the state and control variables. Subsequently, we utilize these proper efficiency criteria to construct two multiobjective dual problems and prove appropriate duality theorems. Also, we specialize and discuss these results for a particular case of our principal problem which contains square roots of positivesemidefinite quadratic forms. As special cases of the main proper efficiency and duality results, this paper also contains similar results for control problems with multiple, fractional, and ordinary objective functions.     

Efficient solutions in V-KT-pseudoinvex multiobjective control problems: A characterization

In this paper, we introduce a new condition on functionals involved in a multiobjective control problem, for which we define the V-KT-pseudoinvex control problem. We prove that a V-KT-pseudoinvex control problem is characterized so that a Kuhn–Tucker point is an efficient solution. We generalize recently obtained optimality results of known mathematical programming problems and control problems. We illustrate these results with an example.     

Superconvergence analysis of fully discrete finite element methods for semilinear parabolic optimal control problems

We study the superconvergence property of fully discrete finite element approximation for quadratic optimal control problems governed by semilinear parabolic equations with control constraints. The time discretization is based on difference methods, whereas the space discretization is done using finite element methods. The state and the adjoint state are approximated by piecewise linear functions and the control is approximated by piecewise constant functions. First, we define a fully discrete finite element approximation scheme for the semilinear parabolic control problem. Second, we derive the superconvergence properties for the control, the state and the adjoint state. Finally, we do some numerical experiments for illustrating our theoretical results.     

Duality for Set-Valued Multiobjective Optimization Problems. Part 2: Optimal Control

The present paper is a continuation of a paper by Azimov (J. Optim. Theory Appl. 2007, accepted), where we derived duality relations for some general multiobjective optimization problems which include convex programming and optimal control problems. As a consequence, we established duality results for multiobjective convex programming problems. In the present paper (Part 2), based on Theorem 3.2 of Azimov (J. Optim. Theory Appl. 2007, accepted), we establish duality results for several classes of multiobjective optimal control problems.     

Neumann Control of Unstable Parabolic Systems: Numerical Approach

The present article is concerned with the Neumann control of systems modeled by scalar or vector parabolic equations of reaction-advection-diffusion type with a particular emphasis on systems which are unstable if uncontrolled. To solve these problems, we use a combination of finite-difference methods for the time discretization, finite-element methods for the space discretization, and conjugate gradient algorithms for the iterative solution of the discrete control problems. We apply then the above methodology to the solution of test problems in two dimensions, including problems related to nonlinear models.     

Automatic differentiation of explicit Runge-Kutta methods for optimal control

This paper considers the numerical solution of optimal control problems based on ODEs. We assume that an explicit Runge-Kutta method is applied to integrate the state equation in the context of a recursive discretization approach. To compute the gradient of the cost function, one may employ Automatic Differentiation (AD). This paper presents the integration schemes that are automatically generated when differentiating the discretization of the state equation using AD. We show that they can be seen as discretization methods for the sensitivity and adjoint differential equation of the underlying control problem. Furthermore, we prove that the convergence rate of the scheme automatically derived for the sensitivity equation coincides with the convergence rate of the integration scheme for the state equation. Under mild additional assumptions on the coefficients of the integration scheme for the state equation, we show a similar result for the scheme automatically derived for the adjoint equation. Numerical results illustrate the presented theoretical results.     

Mixed Discontinuous Galerkin Time-Stepping Method for Linear Parabolic Optimal Control Problems

In this paper, we discuss the mixed discontinuous Galerkin (DG) finite element approximation to linear parabolic optimal control problems. For the state variables and the co-state variables, the discontinuous finite element method is used for the time discretization and the Raviart-Thomas mixed finite element method is used for the space discretization. We do not discretize the space of admissible control but implicitly utilize the relation between co-state and control for the discretization of the control. We derive a priori error estimates for the lowest order mixed DG finite element approximation. Moveover, for the element of arbitrary order in space and time, we derive a posteriori $L^2(0, T ;L^2(Ω))$ error estimates for the scalar functions, assuming that only the underlying mesh is static. Finally, we present an example to confirm the theoretical result on a priori error estimates.     

η-Approximation Method for Non-convex Multiobjective Variational Problems

In this paper, we use the η-approximation method for a class of non-convex multiobjective variational problems with invex functionals. In this approach, for the considered multiobjective variational problem, the associated η-approximated multiobjective variational problem is constructed at the given feasible solution. The equivalence between (weakly) e?cient solutions in the original multiobjective variational problem and its associated η-approximated multiobjective variational problem is established under invexity hypotheses.     

SQP-methods for solving optimal control problems with control and state constraints: adjoint variables, sensitivity analysis and real-time control

Parametric nonlinear optimal control problems subject to control and state constraints are studied. Two discretization methods are discussed that transcribe optimal control problems into nonlinear programming problems for which SQP-methods provide efficient solution methods. It is shown that SQP-methods can be used also for a check of second-order sufficient conditions and for a postoptimal calculation of adjoint variables. In addition, SQP-methods lead to a robust computation of sensitivity differentials of optimal solutions with respect to perturbation parameters. Numerical sensitivity analysis is the basis for real-time control approximations of perturbed solutions which are obtained by evaluating a first-order Taylor expansion with respect to the parameter. The proposed numerical methods are illustrated by the optimal control of a low-thrust satellite transfer to geosynchronous orbit and a complex control problem from aquanautics. The examples illustrate the robustness, accuracy and efficiency of the proposed numerical algorithms.     

Boundary integral approximation of a heat-diffusion problem in time-harmonic regime

In this paper we propose and analyse numerical methods for the approximation of the solution of Helmholtz transmission problems in the half plane. The problems we deal with arise from the study of some models in photothermal science. The solutions to the problem are represented as single layer potentials and an equivalent system of boundary integral equations is derived. We then give abstract necessary and sufficient conditions for convergence of Petrov–Galerkin discretizations of the boundary integral system and show for three different cases that these conditions are satisfied. We extend the results to other situations not related to thermal science and to non-smooth interfaces. Finally, we propose a simple full discretization of a Petrov–Galerkin scheme with periodic spline spaces and show some numerical experiments.     

Discretization methods for the solution of semi-infinite programming problems

In the first part of this paper, we prove the convergence of a class of discretization methods for the solution of nonlinear semi-infinite programming problems, which includes known methods for linear problems as special cases. In the second part, we modify and study this type of algorithms for linear problems and suggest a specific method which requires the solution of a quadratic programming problem at each iteration. With this algorithm, satisfactory results can also be obtained for a number of singular problems. We demonstrate the performance of the algorithm by several numerical examples of multivariate Chebyshev approximation problems.     

Generalized Homotopy Approach to Multiobjective Optimization

This paper proposes a new generalized homotopy algorithm for the solution of multiobjective optimization problems with equality constraints. We consider the set of Pareto candidates as a differentiable manifold and construct a local chart which is fitted to the local geometry of this Pareto manifold. New Pareto candidates are generated by evaluating the local chart numerically. The method is capable of solving multiobjective optimization problems with an arbitrary number k of objectives, makes it possible to generate all types of Pareto optimal solutions, and is able to produce a homogeneous discretization of the Pareto set. The paper gives a necessary and sufficient condition for the set of Pareto candidates to form a (k-1)-dimensional differentiable manifold, provides the numerical details of the proposed algorithm, and applies the method to two multiobjective sample problems.     

Optimal Control of Interface Problems with hp-finite Elements

We investigate the optimal control of elliptic partial differential equations with jumping coefficients. As discretization, we use interface concentrated finite elements on subdomains with smooth data. In order to apply convergence results, we prove higher regularity of the optimal solution using the concept of quasi-monotone coefficients and a domain that is injective modulo polynomials of degree 1 at each vertex. Numerical results are presented for a semi-linear control problem with a non-local radiation operator, which models the production process of silicon carbide single crystals.     

Asymptotic analysis in convex composite multiobjective optimization problems

In this paper, we present a unified approach for studying convex composite multiobjective optimization problems via asymptotic analysis. We characterize the nonemptiness and compactness of the weak Pareto optimal solution sets for a convex composite multiobjective optimization problem. Then, we employ the obtained results to propose a class of proximal-type methods for solving the convex composite multiobjective optimization problem, and carry out their convergence analysis under some mild conditions.     

Symmetric coupling of finite and boundary elements for exterior magnetic field problems

We consider a coupled finite element (fe)–boundary element (be) approach for three‐dimensional magnetic field problems. The formulation is based on a vector potential in a bounded domain (fe) and a scalar potential in an unbounded domain (be). We describe a coupled variational problem yielding a unique solution where the constraints in the trial spaces are replaced by appropriate side conditions. Then we discuss a Galerkin discretization of the coupled problem and prove a quasi‐optimal error estimate. Finally we discuss an efficient preconditioned iterative solution strategy for the resulting linear system. Copyright © 2002 John Wiley & Sons, Ltd.     

Optimization Techniques for Solving Elliptic Control Problems with Control and State Constraints: Part 1. Boundary Control

We study optimal control problems for semilinear elliptic equations subject to control and state inequality constraints. In a first part we consider boundary control problems with either Dirichlet or Neumann conditions. By introducing suitable discretization schemes, the control problem is transcribed into a nonlinear programming problem. It is shown that a recently developed interior point method is able to solve these problems even for high discretizations. Several numerical examples with Dirichlet and Neumann boundary conditions are provided that illustrate the performance of the algorithm for different types of controls including bang-bang and singular controls. The necessary conditions of optimality are checked numerically in the presence of active control and state constraints.     

Generalized asymptotic functions in nonconvex multiobjective optimization problems

In this paper, we use generalized asymptotic functions and second-order asymptotic cones to develop a general existence result for the nonemptiness of the proper efficient solution set and a sufficient condition for the domination property in nonconvex multiobjective optimization problems. A new necessary condition for a point to be efficient or weakly efficient solution is given without any convexity assumption. We also provide a finer outer estimate for the asymptotic cone of the weakly efficient solution set in the quasiconvex case. Finally, we apply our results to the linear fractional multiobjective optimization problem.