Superconvergence for Optimal Control Problems Governed by Semi-linear Elliptic Equations

In this paper, we will investigate the superconvergence of the finite element approximation for quadratic optimal control problem governed by semi-linear elliptic equations. The state and co-state variables are approximated by the piecewise linear functions and the control variable is approximated by the piecewise constant functions. We derive the superconvergence properties for both the control variable and the state variables. Finally, some numerical examples are given to demonstrate the theoretical results. This work is supported by National Science Foundation of China, the National Basic Research Program under the Grant 2005CB321703, Scientific Research Fund of Hunan Provincial Education Department, and Hunan Provincial Innovation Foundation for Postgraduate (No. S2008yjscx04).     

A Priori Error Estimates for Optimal Control Problems Governed by Transient Advection-Diffusion Equations

In this paper, we investigate a characteristic finite element approximation of quadratic optimal control problems governed by linear advection-dominated diffusion equations, where the state and co-state variables are discretized by piecewise linear continuous functions and the control variable is approximated by piecewise constant functions. We derive some a priori error estimates for both the control and state approximations. It is proved that these approximations have convergence order , where h _U and h are the spatial mesh-sizes for the control and state discretization, respectively, and k is the time increment. Numerical experiments are presented, which verify the theoretical results. This research was supported by the National Basic Research Program of China (No. 2007CB814906) and the National Natural Science Foundation of China (No. 10771124).     

A Strongly Conservative Hybrid DG/Mixed FEM for the Coupling of Stokes and Darcy Flow

We consider the coupling of free and porous media flow governed by Stokes and Darcy equations with the Beavers–Joseph–Saffman interface condition. This model is discretized using a divergence-conforming finite element for the velocities in the whole domain. Hybrid discontinuous Galerkin techniques and mixed methods are used in the Stokes and Darcy subdomains, respectively. The discretization achieves mass conservation in the sense of \(H(\mathrm {div},\Omega )\), and we obtain optimal velocity convergence. Numerical results are presented to validate the theoretical findings.     

A Mixed Finite Element Scheme for Optimal Control Problems with Pointwise State Constraints

In this paper, we propose a mixed variational scheme for optimal control problems with point-wise state constraints, the main idea is to reformulate the optimal control problems to a constrained minimization problem involving only the state, which is characterized by a fourth order variational inequality. Then mixed form based on this fourth order variational inequality is formulated and a direct numerical algorithm is proposed without the optimality conditions of underlying optimal control problems. The a priori and a posteriori error estimates are proved for the mixed finite element scheme. Numerical experiments confirm the efficiency of the new strategy.     

A New Nonsymmetric Discontinuous Galerkin Method for Time Dependent Convection Diffusion Equations

We propose a discontinuous Galerkin finite element method for convection diffusion equations that involves a new methodology handling the diffusion term. Test function derivative numerical flux term is introduced in the scheme formulation to balance the solution derivative numerical flux term. The scheme has a nonsymmetric structure. For general nonlinear diffusion equations, nonlinear stability of the numerical solution is obtained. Optimal kth order error estimate under energy norm is proved for linear diffusion problems with piecewise P ^k polynomial approximations. Numerical examples under one-dimensional and two-dimensional settings are carried out. Optimal (k+1)th order of accuracy with P ^k polynomial approximations is obtained on uniform and nonuniform meshes. Compared to the Baumann-Oden method and the NIPG method, the optimal convergence is recovered for even order P ^k polynomial approximations.     

Regional Optimal Control Problem Governed by Distributed Bi-linear Hyperbolic Systems

This paper considers the regional bi-linear control problem of an important class of hyperbolic systems. The objective is to bring the state solutions at time T close to a desired observations w ^d only on a sub-region ω along the spatial domain Ω. We prove the existence of solution by minimizing sequence method. The adjoint system of this problem is introduced and used to characterize the optimal control. A numerical approach is developed and illustrated successfully by simulations.     

A Superconvergent HDG Method for Distributed Control of Convection Diffusion PDEs

We consider a distributed optimal control problem governed by an elliptic convection diffusion PDE, and propose a hybridizable discontinuous Galerkin method to approximate the solution. We use polynomials of degree \(k+1\) to approximate the state and dual state, and polynomials of degree \(k \ge 0\) to approximate their fluxes. Moreover, we use polynomials of degree k to approximate the numerical traces of the state and dual state on the faces, which are the only globally coupled unknowns. We prove optimal a priori error estimates for all variables when \( k \ge 0 \). Furthermore, from the point of view of the number of degrees of freedom of the globally coupled unknowns, this method achieves superconvergence for the state, dual state, and control when \(k\ge 1\). We illustrate our convergence results with numerical experiments.     

Two Mixed Finite Element Methods for Time-Fractional Diffusion Equations

Based on spatial conforming and nonconforming mixed finite element methods combined with classical L1 time stepping method, two fully-discrete approximate schemes with unconditional stability are first established for the time-fractional diffusion equation with Caputo derivative of order \(0<\alpha <1\). As to the conforming scheme, the spatial global superconvergence and temporal convergence order of \(O(h^2+\tau ^{2-\alpha })\) for both the original variable u in \(H^1\)-norm and the flux \(\vec {p}=\nabla u\) in \(L^2\)-norm are derived by virtue of properties of bilinear element and interpolation postprocessing operator, where h and \(\tau \) are the step sizes in space and time, respectively. At the same time, the optimal convergence rates in time and space for the nonconforming scheme are also investigated by some special characters of \(\textit{EQ}_1^{\textit{rot}}\) nonconforming element, which manifests that convergence orders of \(O(h+\tau ^{2-\alpha })\) and \(O(h^2+\tau ^{2-\alpha })\) for the original variable u in broken \(H^1\)-norm and \(L^2\)-norm, respectively, and approximation for the flux \(\vec {p}\) converging with order \(O(h+\tau ^{2-\alpha })\) in \(L^2\)-norm. Numerical examples are provided to demonstrate the theoretical analysis.     

A Mixed DG Method for Linearized Incompressible Magnetohydrodynamics

We introduce and analyze a discontinuous Galerkin method for the numerical discretization of a stationary incompressible magnetohydrodynamics model problem. The fluid unknowns are discretized with inf-sup stable discontinuous ? _k ³ ?? _k?1 elements whereas the magnetic part of the equations is approximated by discontinuous ? _k ³ ?? _k+1 elements. We carry out a complete a-priori error analysis of the method and prove that the energy norm error is convergent of order k in the mesh size. These results are verified in a series of numerical experiments.     

A Lagrangian Scheme for the Solution of Nonlinear Diffusion Equations Using Moving Simplex Meshes

A Lagrangian numerical scheme for solving nonlinear degenerate Fokker–Planck equations in space dimensions \(d\ge 2\) is presented. It applies to a large class of nonlinear diffusion equations, whose dynamics are driven by internal energies and given external potentials, e.g. the porous medium equation and the fast diffusion equation. The key ingredient in our approach is the gradient flow structure of the dynamics. For discretization of the Lagrangian map, we use a finite subspace of linear maps in space and a variational form of the implicit Euler method in time. Thanks to that time discretisation, the fully discrete solution inherits energy estimates from the original gradient flow, and these lead to weak compactness of the trajectories in the continuous limit. Consistency is analyzed in the planar situation, \(d=2\). A variety of numerical experiments for the porous medium equation indicates that the scheme is well-adapted to track the growth of the solution’s support.     

H_2/l_1混合优化问题的凸二次规划解法

采用上逼近算法求解 H2 / l1 混合优化问题。首先将其转化为有限维的凸二次规划问题 ,并利用L emke互补转轴算法求解 ;然后逐次进行逼近。计算示例表明所得结果是正确的     

A Semi-smooth Newton Method for Regularized State-constrained Optimal Control of the Navier-Stokes Equations

In this paper, we study semi-smooth Newton methods for the numerical solution of regularized pointwise state-constrained optimal control problems governed by the Navier-Stokes equations. After deriving an appropriate optimality system for the original problem, a class of Moreau-Yosida regularized problems is introduced and the convergence of their solutions to the original optimal one is proved. For each regularized problem a semi-smooth Newton method is applied and its local superlinear convergence verified. Finally, selected numerical results illustrate the behavior of the method and a comparison between the max-min and the Fischer-Burmeister as complementarity functionals is carried out.     

Multivariable Adaptive Control of the Rewinding Process of a Roll-to-roll System Governed by Hyperbolic Partial Differential Equations

In this paper, an active control scheme for the rewinding process of a roll-to-roll (R2R) system is investigated. The control objectives are to suppress the transverse vibration of the moving web, to track the desired velocity profile, and to keep the desired radius value of a rewind roller. The bearing coefficient in the rewind shaft is unknown and the rotating elements in the drive motor are various. The moving web is modeled as an axially moving beam system governed by hyperbolic partial differential equations (PDEs). The control scheme utilizes two control inputs: a control force exerted from a hydraulic actuator equipped with a damper, and a control torque applied to the rewind roller. Two adaptation laws are derived to estimate the unknown bearing coefficient and the bound of variations of the rotating elements. The Lyapunov method is employed to prove the robust stability of the rewind section, specifically the uniform and ultimate boundedness of all of the signals. The effectiveness of the proposed control schemes was verified by numerical simulations.     

On Parabolic Boundary Layers for Convection–Diffusion Equations in a Channel: Analysis and Numerical Applications

In this article we discuss singularly perturbed convection–diffusion equations in a channel in cases producing parabolic boundary layers. It has been shown that one can improve the numerical resolution of singularly perturbed problems involving boundary layers, by incorporating the structure of the boundary layers into the finite element spaces, when this structure is available; see e.g. [Cheng, W. and Temam, R. (2002). Comput. Fluid. V.31, 453–466; Jung, C. (2005). Numer. Meth. Partial Differ. Eq. V.21, 623–648]. This approach is developed in this article for a convection–diffusion equation. Using an analytical approach, we first derive an approximate (simplified) form of the parabolic boundary layers (elements) for our problem; we then develop new numerical schemes using these boundary layer elements. The results are performed for the perturbation parameter ε in the range 10⁻¹–10⁻¹⁵ whereas the discretization mesh is in the range of order 1/10–1/100 in the x-direction and of order 1/10–1/30 in the y-direction. Indications on various extensions of this work are briefly described at the end of the Introduction.Dedicated to David Gottlieb on his 60th birthday.     

A Discrete Scheme for Parametric Anisotropic Surface Diffusion

In this note we present, how anisotropic surface energies may be incorporated into the finite element method for parametric surface diffusion given by Bänsch et al. [2004. J. Comput. Phys. 203, 321–343]. We present the adapted variational formulation, and the resulting semi-implicit discretization. Finally several simulations with strong (convex) anisotropies are shown, where the corresponding Wulff shapes are approached as the steady state