A shape optimization formulation of weld pool determination

In this paper, we propose a shape optimization formulation for a problem modeling a process of welding. We show the existence of an optimal solution. The finite element method is used for the discretization of the problem. The discrete problem is solved by an identification technique using a parameterization of the weld pool by Bézier curves and Genetic algorithms.     

A global superconvergent $L^{\infty}$-error estimate of mixed finite element methods for semilinear elliptic optimal control problems

In this paper, we discuss the superconvergence of mixed finite element methods for a semilinear elliptic control problem with an integral constraint. The state and costate are approximated by the order $k=1$ Raviart-Thomas mixed finite element spaces and the control variable is approximated by piecewise constant functions. Approximation of the optimal control of the continuous optimal control problem will be constructed by a projection of the discrete adjoint state. It is proved that this approximation has convergence order $h^{2}$ in $L^{\infty}$-norm. Finally, a numerical example is given to demonstrate the theoretical results.     

Superconvergence and a posteriori error estimates of RT1 mixed methods for elliptic control problems with an integral constraint

In this paper, we investigate the superconvergence property and a posteriori error estimates of mixed finite element methods for a linear elliptic control problem with an integral constraint. The state and co-state are approximated by the order k = 1 Raviart-Thomas mixed finite element spaces and the control variable is approximated by piecewise constant functions. Approximations of the optimal control of the continuous optimal control problem will be constructed by a projection of the discrete adjoint state. It is proved that these approximations have convergence order h ². Moreover, we derive a posteriori error estimates both for the control variable and the state variables. Finally, a numerical example is given to demonstrate the theoretical results.     

Control the Fibre Orientation Distribution at the Outlet of Contraction

A two dimensional model of the orientation distribution of fibres in a paper machine headbox is studied. The goal is to control the fibre orientation distribution at the outlet of contraction by changing its shape. The mathematical formulation leads to an optimization problem with control in coefficients of a linear convection-diffusion equation as the state problem. Then, the problem is expressed as an optimal control problem governed by variational forms. By using an embedding method, the class of admissible shapes is replaced by a class of positive Radon measures. The optimization problem in measure space is then approximated by a linear programming problem. The optimal measure representing optimal shape is approximated by the solution of this linear programming problem. In this paper, we have shown that the embedding method (embedding the admissible set into a subset of measures), successfully can be applied to shape variation design to a one dimensional headbox. The usefulness of this idea is that the method is not iterative and it does not need any initial guess of the solution.      

Linear Formulation of a Distributed Boundary Control Problem

In this paper, we consider a distributed boundary control problem governed by an elliptic partial differential equation with state constraints and a minimax objective function. The continuous optimal control problem, discretized with the finite element method, is numerically approximated by a family of linear programming problems. Application to an optimal configuration problem is discussed.     

Coupling of FEM and BEM in Shape Optimization

In the present paper we consider the numerical solution of shape optimization problems which arise from shape functionals of integral type over a compact region of the unknown shape, especially L ²-tracking type functionals. The underlying state equation is assumed to satisfy a Poisson equation with Dirichlet boundary conditions. We proof that the shape Hessian is not strictly H ^1/2-coercive at the optimal domain which implies ill-posedness of the optimization problem under consideration. Since the adjoint state depends directly on the state, we propose a coupling of finite element methods (FEM) and boundary element methods (BEM) to realize an efficient first order shape optimization algorithm. FEM is applied in the compact region while the rest is treated by BEM. The coupling of FEM and BEM essentially retains all the structural and computational advantages of treating the free boundary by boundary integral equations.This research has been carried out when the second author stayed at the Department of Mathematics, Utrecht University, The Netherlands, supported by the EU-IHP project Nonlinear Approximation and Adaptivity: Breaking Complexity in Numerical Modelling and Data Representation     

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.     

Error estimates of mixed finite element methods for quadratic optimal control problems

In this paper, we investigate the error estimates for the solutions of optimal control problems by mixed finite element methods. The state and costate are approximated by Raviart-Thomas mixed finite element spaces of order k and the control is approximated by piecewise polynomials of order k. Under the special constraint set, we will show that the control variable can be smooth in the whole domain. We derive error estimates of optimal order both for the state variables and the control variable.     

A Posteriori Error Estimates of Mixed Methods for Parabolic Optimal Control Problems

In this article, we shall give a brief review on the fully discrete mixed finite element method for general optimal control problems governed by parabolic equations. The state and the co-state are approximated by the lowest order Raviart–Thomas mixed finite element spaces and the control is approximated by piecewise constant elements. Furthermore, we derive a posteriori error estimates for the finite element approximation solutions of optimal control problems. Some numerical examples are given to demonstrate our theoretical results.     

Shape Optimization for Navier–Stokes Equations with Algebraic Turbulence Model: Existence Analysis

We study a shape optimization problem for the paper machine headbox which distributes a mixture of water and wood fibers in the paper making process. The aim is to find a shape which a priori ensures the given velocity profile on the outlet part. The mathematical formulation leads to an optimal control problem in which the control variable is the shape of the domain representing the header, the state problem is represented by a generalized stationary Navier–Stokes system with nontrivial mixed boundary conditions. In this paper we prove the existence of solutions both to the generalized Navier–Stokes system and to the shape optimization problem.     

An approach for solving of a moving boundary problem

In this paper we shall study moving boundary problems, and we introduce an approach for solving a wide range of them by using calculus of variations and optimization. First, we transform the problem equivalently into an optimal control problem by defining an objective function and artificial control functions. By using measure theory, the new problem is modified into one consisting of the minimization of a linear functional over a set of Radon measures; then we obtain an optimal measure which is then approximated by a finite combination of atomic measures and the problem converted to an infinite-dimensional linear programming. We approximate the infinite linear programming to a finite-dimensional linear programming. Then by using the solution of the latter problem we obtain an approximate solution for moving boundary function on specific time. Furthermore, we show the path of moving boundary from initial state to final state.     

On an optimal shape design problem to control a thermoelastic deformation under a prescribed thermal treatment

A shape optimization problem concerned with thermal deformation of elastic bodies is considered. In this article, measure theory approach in function space is derived, resulting in an effective algorithm for the discretized optimization problem. First the problem is expressed as an optimal control problem governed by variational forms on a fixed domain. Then by using an embedding method, the class of admissible shapes is replaced by a class of positive Borel measures. The optimization problem in measure space is then approximated by a linear programming problem. The optimal measure representing optimal shape is approximated by the solution of this finite-dimensional linear programming problem. Numerical examples are also given.     

Superconvergence of triangular mixed finite elements for optimal control problems with an integral constraint

In this paper, we shall investigate the superconvergence property of quadratic elliptical optimal control problems by triangular mixed finite element methods. The state and co-state are approximated by the order k = 1 Raviart-Thomas mixed finite elements and the control is discretized by piecewise constant functions. We prove the superconvergence error estimate of h² in L²-norm between the approximated solution and the interpolation of the exact control variable. Moreover, by postprocessing technique, we find that the projection of the discrete adjoint state is superclose (in order h²) to the exact control variable.     

n-sided polygonal hybrid finite elements with unified fundamental solution kernels for topology optimization

In topology optimization, the optimized design can be obtained based on spatial discretization of design domain using natural polygonal finite elements to reduce the influence of mesh geometry on topology optimization solutions. However, the natural polygonal finite elements require separate interpolants for each type of elements and involve troublesome domain integrals. In this study, an alternative n-sided polygonal hybrid finite element possessing multiple-node connection is formulated in a unified form to compress the checkerboard patterns caused by numerical instability in topology optimization. Different from the natural polygonal finite elements, the present polygonal hybrid finite elements involve two sets of independent displacement fields. The intra-element displacement field defined inside the element is approximated by the linear combination of the fundamental solution of the problem to achieve the purpose of the local satisfaction of the governing equations of the problem, but not the specific boundary conditions and the inter-element continuity conditions. To overcome such drawback, the inter-element displacement field defined over the entire element boundary is independently approximated by means of the conventional shape function interpolation. As a result, only line integrals along the element boundary are involved in the computation, whose dimension is reduced by one compared to the domain integrals in the natural polygonal finite elements, and more importantly, allowing us to flexibly construct any polygons from Voronoi tessellations in discretizing complex design domains using same fundamental solution kernels. Numerical results obtained indicate that the present n-sided polygonal hybrid finite elements can produce more accurate displacement solutions and smaller mean compliance, compared to the standard finite elements and the natural polygonal finite elements.     

A priori error estimates for optimal control problems with pointwise constraints on the gradient of the state

We analyze a finite element approximation of an elliptic optimal control problem with pointwise bounds on the gradient of the state variable. We derive convergence rates if the control space is discretized implicitly by the state equation. In contrast to prior work we obtain these results directly from classical results for the W ^1,∞-error of the finite element projection, without using adjoint information. If the control space is discretized directly, we first prove a regularity result for the optimal control to control the approximation error, based on which we then obtain analogous convergence rates.     

Finite element approximation of elliptic control problems with constraints on the gradient

We consider an elliptic optimal control problem with control constraints and pointwise bounds on the gradient of the state. We present a tailored finite element approximation to this optimal control problem, where the cost functional is approximated by a sequence of functionals which are obtained by discretizing the state equation with the help of the lowest order Raviart–Thomas mixed finite element. Pointwise bounds on the gradient variable are enforced in the elements of the triangulation. Controls are not discretized. Error bounds for control and state are obtained in two and three space dimensions. A numerical example confirms our analytical findings.     

半线性抛物最优控制问题全离散插值系数有限元方法的收敛性分析

本文考虑全离散插值系数有限元方法求解半线性抛物最优控制问题,其中控制变量用分片常数函数逼近,状态变量和对偶状态变量用分片线性函数逼近.对于方程中的半线性项,先用插值系数技巧处理,再用牛顿迭代法求解.通过引入一些辅助变量和投影算子,并利用有限元空间的逼近性质,得到半线性抛物最优控制问题插值系数有限元方法的收敛性结果;数值算例结果验证了理论结果的正确性.     

Solving some optimal control problems using the barrier penalty function method

In this paper we present a new approach to solve a two-level optimization problem arising from an approximation by means of the finite element method of optimal control problems governed by unilateral boundary-value problems. The problem considered is to find a minimum of a functional with respect to the control variablesu. The minimized functional depends on control variables and state variablesx. The latter are the optimal solution of an auxiliary quadratic programming problem, whose parameters depend onu.Our main idea is to replace this QP problem by its dual and then apply the barrier penalty method to this dual QP problem or to the primal one if it is in an appropriate form. As a result we obtain a problem approximating the original one. Its good property is the differentiable dependence of state variables with respect to the control variables. Furthermore, we propose a method for finding an approximate solution of a penalized lower-level problem if the optimal solution of the original QP problem is known. We apply the result obtained to some optimal shape design problems governed by the Dirichlet-Signorini boundary-value problem.This research was supported by the Academy of Finland and the Systems Research Institute of the Polish Academy of Sciences.     

Shape sensitivity analysis of an elastic contact problem: Convergence of the Nitsche based finite element approximation

In a recent work, we introduced a finite element approximation for the shape optimization of an elastic structure in sliding contact with a rigid foundation where the contact condition (Signorini’s condition) is approximated by Nitsche’s method and the shape gradient is obtained via the adjoint state method. The motivation of this work is to propose an a priori convergence analysis of the numerical approximation of the variables of the shape gradient (displacement and adjoint state) and to show some numerical results in agreement with the theoretical ones. The main difficulty comes from the non-differentiability of the contact condition in the classical sense which requires the notion of conical differentiability.     

Study of a Numerical Approach for a Transient Flow Problem in Porous Media

In this work, we deal with the numerical study of the new approximation method proposed in [7] for a transient flow problem in porous media. The stationary problem, obtained from a time discretization of this transient problem, is considered as an optimal shape design formulation. We prove the existence of the solution of the discrete optimal shape problem obtained from finite element discretization. We study the convergence and give numerical results showing the efficiency of the proposed approach.