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.     

An Approach to Shape Optimization with State Constraints

We introduce a new approach to solve shape optimization problems with state constraints. The problem is reformulated on a fixed reference domain using conformal pull-back. The shape dependence is then hidden in the conformal parameter. The problem is discretized using FEM and solved as an NLP. Finally the optimal shape can be reconstructed from the optimal conformal parameter. (© 2012 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Gradient evaluation in DAE optimal control problems by sensitivity equations and adjoint equations

The problem of calculating gradients for discretized optimal control problems is addressed. These gradients are needed by a sequential quadratic programming (SQP) method, which is used to solve the discretized optimal control problem. Two approaches for gradient calculation are discussed and compared numerically. The first approach uses sensitivity equations, the second uses adjoint equations. It turns out that it depends on the problem dimensions which approach is more preferable compared to the other. (© 2005 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Optimal shape design for systems governed by variational inequalities,part 2: Existence theory for the evolution case

Some general existence results for optimal shape design problems for systems governed by parabolic variational inequalities are established by the mapping method and variational convergence theory. Then, an existence theorem is given for the optimal shape for an electrochemical machining problem, in which the cost functional is not lower semicontinuous, by extending the general results to this case. Furthermore, this problem is approximated by a set of optimal shape design problems which have more smooth cost functionals and are easier to handle computationally.The authors wish to express their sincere thanks to the reviewers for supplying additional references and for their valuable comments, which made the paper more readable.     

Mathematical analysis of the optimal location of wastewater outfalls

In this work we deal with the design of wastewater treatmentsystems, mainly the optimal placement of underwater outfalls.This problem can be formulated as a state constrained optimalcontrol problem where the control is the position of the outfalls,the cost function is the sum of the distances to the wastewaterfarms and the state equations are those modelling dissolvedoxygen and biochemical oxygen demand concentrations. We discretizethe problem by means of a characteristic Galerkin method andwe propose an interior point algorithm for the numerical resolutionof the discretized optimization problem.     

Optimal shape design for systems governed by variational inequalities,part 1: Existence theory for the elliptic case

Some general existence results for optimal shape design problems for systems governed by elliptic variational inequalities are established by the mapping method and variational convergence theory. Then, an existence theorem is given for the optimal shape for an electrochemical machining problem, in which the cost functional is not lower semicontinuous, by extending the general results to this case. Furthermore, this problem is approximated by a set of optimal shape design problems which have more smooth cost functionals and are easier to handle computationally.The authors with to express their sincere thanks to the reviewers for supplying additional references and for their valuable comments, which made the paper more readable.     

Implementation of Some Methods of Shape Design for Variational Inequalities

In this paper we present some results concerning the optimal shape design problem governed by the fourth-order variational inequalities. The problem can be considered as a model example for the design of the shapes for elastic-plastic problem. The computations are done by finite element method, and the performance criterion is minimized by the material derivative method. We also discuss the error estimates in the appropriate norm and present some numerical results. An example is used to clearly illustrate the essential elements of shape design problems.     

The shape optimization of the arterial graft design by level set methods

The goal of the arterial graft design problem is to find an optimal graft built on an occluded artery, which can be mathematically modeled by a fluid based shape optimization problem. The smoothness of the graft is one of the important aspects in the arterial graft design problem since it affects the flow of the blood significantly. As an attractive design tool for this problem, level set methods are quite efficient for obtaining better shape of the graft. In this paper, a cubic spline level set method and a radial basis function level set method are designed to solve the arterial graft design problem. In both approaches, the shape of the arterial graft is implicitly tracked by the zero-level contour of a level set function and a high level of smoothness of the graft is achieved. Numerical results show the efficiency of the algorithms in the arterial graft design.     

Optimality criteria coupled adaptive mesh method for optimal shape design of Stokes flow

An adaptive mesh method combined with the optimality criteria algorithm is applied to optimal shape design problems of fluid dynamics. The shape sensitivity analysis of the cost functional is derived. The optimization problem is solved by a simple but robust optimality criteria algorithm, and an automatic local adaptive mesh refinement method is proposed. The mesh adaptation, with an indicator based on the material distribution information, is itself shown as a shape or topology optimization problem. Taking advantages of this algorithm, the optimal shape design problem concerning fluid flow can be solved with higher resolution of the interface and a minimum of additional expense. Details on the optimization procedure are provided. Numerical results for two benchmark topology optimization problems are provided and compared with those obtained by other methods. Copyright © 2016 John Wiley & Sons, Ltd.     

Finite Difference Method for an Optimal Control Problem for a Nonlinear Time-dependent Schrödinger Equation

The finite difference method is applied to an optimal control problem for a system governed by a nonlinear Schrödinger equation with a complex coe?cient. The optimal control problem is discretized by the finite difference method, the error estimate for the finite difference scheme is established and the convergence of difference approximations of the optimal control according to the functional is proved.     

Optimal shape design of systems governed by variational inequalities

In this paper, we present some results concerning the optimal shape design problem governed by the variational inequalities of the fourth order. This problem can be considered as a model example for the design of the shape for elastic-plastic problem. The performance criterion is minimized by the material derivative method.     

Finite element analysis for general elastic multi-structures

A finite element method is introduced to solve the general elastic multi-structure problem, in which the displacements on bodies, the longitudinal displacements on plates and the longitudinal displacements on beams are discretized using conforming linear elements, the rotational angles on beams are discretized using conforming elements of second order, the transverse displacements on plates and beams are discretized by the Morley elements and the Hermite elements of third order, respectively. The generalized Korn's inequality is established on related nonconforming element spaces, which implies the unique solvability of the finite element method. Finally, the optimal error estimate in the energy norm is derived for the method.     

Error Estimates for the Approximation of a Class of Optimal Control Systems Governed by Linear PDEs

This paper deals with a class of optimal control problems in which the system is governed by a linear partial differential equation and the control is distributed and with constraints. The problem is posed in the framework of the theory of optimal control of systems. A numerical method is proposed to approximate the optimal control. In this method, the state space as well as the convex set of admissible controls are discretized. An abstract error estimate for the optimal control problem is obtained that depends on both the approximation of the state equation and the space of controls. This theoretical result is illustrated by some numerical examples from the literature.     

On the convergence of a numerical method for optimal control problems

Under suitable restrictions, convergence to the solutions of a class of optimal control problems is proved for a method by Rosen in which the differential equations, constraints, and cost functional are discretized, and the resulting mathematical programming problem is solved approximately by a penalty-function approach.     

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.     

Numerical simulation of fluid-structure interaction problems on hybrid meshes with algebraic multigrid methods

Fluid-structure interaction problems arise in many fields of application such as flows around elastic structures and blood flow in arteries. The method presented in this paper for solving such a problem is based on a reduction to an equation at the interface, involving the so-called Steklov-Poincaré operators. This interface equation is solved by a Newton iteration, for which directional derivatives involving shape derivatives with respect to the interface perturbation have to be evaluated appropriately. One step of the Newton iteration requires the solution of several decoupled linear sub-problems in the structure and the fluid domains. These sub-problems are spatially discretized by a finite element method on hybrid meshes. For the time discretization, implicit first-order methods are used for both sub-problems. The discretized equations are solved by algebraic multigrid methods.     

Sub-exponential graph coloring algorithm for stencil-based Jacobian computations

Partial differential equations can be discretized using a regular Cartesian grid and a stencil-based method to approximate the partial derivatives. The computational effort for determining the associated Jacobian matrix can be reduced. This reduction can be modeled as a (grid) coloring problem. Currently, this problem is solved by using a heuristic approach for general graphs or by developing a formula for every single stencil. We introduce a sub-exponential algorithm using the Lipton–Tarjan separator in a divide-and-conquer approach to compute an optimal coloring. The practical relevance of the algorithm is evaluated when compared with an exponential algorithm and a greedy heuristic.     

Suboptimal control of a 2D molten carbonate fuel cell PDAE model

Numerical results for suboptimal boundary control of an integro partial differential algebraic equation system of dimension 28 are presented. The application is a molten carbonate fuel cell power plant. The technically and economically important fast tracking of the new stationary cell voltage during a load change is optimized. The nonstandard optimal control problem s.t. degenerated PDE is discretized by the method of lines yielding a very large DAE constrained standard optimal control problem. An index analysis is performed to identify consistent initial conditions.     

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.      

An Alternative to Tikhonov Regularization for Linear Sampling Methods

The problem of determining the shape of an obstacle from far-field measurements is considered. It is well known that linear sampling methods have been widely used for shape reconstructions obtained via the singular system of an ill conditioned discretized far-field operator. For our reconstructions we assume that the far-field data are noisy and we employ a novel regularization method that does not require determination of a regularization parameter.