Applied Mathematics
Vol.07 No.14(2016), Article ID:70028,24 pages
10.4236/am.2016.714134
On Henig Regularization of Material Design Problems for Quasi-Linear p-Biharmonic Equation
Peter Kogut1, Günter Leugering2, Ralph Schiel2
1Department of Differential Equations, Dnipropetrovsk National University, Dnipro, Ukraine
2Department Mathematik, Lehrstuhl II Universität Erlangen-Nürnberg Cauerstr, Erlangen, Germany

Copyright © 2016 by authors and Scientific Research Publishing Inc.
This work is licensed under the Creative Commons Attribution International License (CC BY).
http://creativecommons.org/licenses/by/4.0/


Received 25 June 2016; accepted 21 August 2016; published 24 August 2016
ABSTRACT
We study a Dirichlet optimal design problem for a quasi-linear monotone p-biharmonic equation with control and state constraints. We take the coefficient of the p-biharmonic operator as a design variable in
. In this article, we discuss the relaxation of such problem.
Keywords:
p-Biharmonic Problem, Optimal Design, Relaxation, Henig Dilating Cone, Existence Result

1. Introduction
The aim of this article is to analyze the following optimal design problem (OCP), which can be regarded as an optimal control problem, for quasi-linear partial differential equation (PDE) with mixed boundary conditions
(1)
subject to the quasi-linear equation
(2)
(3)
the pointwise state constraints
(4)
and the design (control) constraints
(5)
Here,
and
are the disjoint part of the boundary
(
),
stands for the control space,
, f, and
are given distributions. Problems of this type appear for p-power-like elastic isotropic flat plates of uniform thickness, where the design variable u is to be chosen such that the deflection of the plate matches a given profile. The model extends the classical weighted biharmonic equation, where the weight
involves the thickness a of the plate, see e.g. [1] - [3] , or u can be regarded as a rigidity parameter. The OCP (1)-(4) can be considered as a prototype of design problems for quasilinear state equations. For an interesting exposure to this subject we can refer to the monographs [4] - [6] .
A particular feature of OCP (1)-(4) is the restriction by the pointwise constraints (4) in
-space. In fact, the ordering cone of positive elements in 











As was shown in the recent publication [7] , the proposed approach is numerically viable for state-constrained optimal control problems with the state equation given by linear partial differential equations. In particular, using the finite element discretization of the Henig dilating cone of positive functions, it has been shown in [7] that the above approximation scheme, called conical regularization, where the regularization is done by replacing the ordering cone with a family of dilating cones, leads to a finite-dimensional optimization problem which can conveniently be treated by known numerical techniques. The non-emptiness of the feasible set for the state- constrained OCPs is an open question even for the simplest situation. Therefore, we consider a more flexible notion of solution to the boundary value problem (2)-(3). With that in mind we discuss a variant of the penalization approach, called the “variational inequality (VI) method”. Following this approach we weaken the requirements on admissible solutions to the original OCP and consider instead the family of penalized OCPs for appropriate variational inequalities
where the sets 

The outline of the paper is the following. In Section 2 we report some preliminaries and notation we need in the sequel. In Sections 3, we give a precise statement of the state constrained optimal control (or design) problem and describe the main assumptions on the initial data and control functions. In Section 4, we provide the results concerning solvability of the original problem with control and state constraints. We show that this problem admits at least one solution if and only if the corresponding set of feasible solutions is nonempty. In Section 5 we show that the pointwise state constraints can be replaced by the weakened conditions coming from Henig relaxation of ordering cones. As a result, we give a precise definition of the relaxed optimization problems and show that the solvability of the original OCP can be characterized by the associated relaxed problems. In particular, we prove that the optimal solution to the original problem can be attained in the limit by the optimal solution of the relaxed problem. We consider in Section 6 the “variational inequality method” as an approximation of the OCPs. Following this approach, we weaken the requirements on feasible solutions to the original OCP. In contrast to the Henig relaxation approach, the penalized optimal control problem for indicated variational inequality has a non-empty feasible set and this problem is always solvable. In conclusion, we show that some of the optimal solutions to the original problem can be attained in the limit by optimal solutions of the penalized problem. However, it is unknown whether the entire set of the optimal solutions can be attained in such way.
2. Definitions and Basic Properties
Let 













Let p be a real number such that






where
For any 
By ( [9] , Theorem 8.3), these linear operators can be extended continuously to the whole of space
as closed subspaces of 


are compact.
Let
as the closure of 





of 

Throughout this paper, we use the notation


is a uniformly convex Banach space [10] . Moreover, the norm 






is uniquely solvable in 


This allows us to conclude the following. If 


and y is a solution of (8), then




for a suitable positive constant 
Thus, by the Closed Graph Theorem, we can conclude that 

By 

is finite.
We recall that a sequence 





It is well-known also the following compactness result for BV-spaces (Helly’s selection theorem, see [15] ).
Theorem 1. If 

strongly converging in 







3. Setting of the Optimal Control Problem
Let



where 
Let 


1) the function 


2) the function 

In addition, the following conditions of subcritical growth, monotonicity, and non-negativity are fulfilled:



for some
is the critical exponent for the Sobolev imbedding





Let






by choosing an appropriate function 
is the operator of fourth order called the generalized p-biharmonic operator, and the class of admissible controls 

It is clear that 

More precisely, we are concerned with the following optimal control problem

Before we will discuss the question of existence of admissible pairs to the problem (19), we note that the function 


Moreover, taking into account the growth condition (12) and the compactness of the Sobolev imbedding 


Definition 3.1. We say that an element 


Remark 3.1. Since the set 




Passing to the limit as 

Hence,

and we arrive at the standard definition of weak solution to the boundary value problem (15)-(16). However, in order to avoid some mathematical difficulties, we will mainly use the Minty inequality in our further analysis. It is
worth to note that having applied Green’s formula twice to operator 

we arrive at the identity
Hence, if y as an element of 
In particular, taking 

As a result, conditions (11), (18), and inequalities (14) and (9) lead us to the following a priori estimate

The existence of a unique weak solution to the boundary value problem (15)-(16) in the sense of Definition 3.1 follows from an abstract theorem on monotone operators.
Theorem 2 ( [17] ) Let V be a reflexive separable Banach space. Let V* be the dual space, and let 


Here, the above mentioned properties of the strict monotonicity, hemicontinuity, and coercivity of the operator A have respectively the following meaning:




In our case, we can define the operator 


In view of the properties (12)-(14) and compactness of the Sobolev imbedding 




for which 



Thus,
and, hence, in view of Remark 3.1, the Minty relation (21) holds true.
Taking this fact into account, we adopt the following notion.
Definition 3.2. We say that 





where 


We denote by 

Remark 3.2. Before we proceed further, we need to make sure that minimization problem (19) is meaningful, i.e. there exists at least one pair 





It is reasonably now to make use of the following Hypothesis.
(H1) There exists at least one pair 

4. Existence of Optimal Solutions
In this section we focus on the solvability of optimal control problem (15)-(19). Hereinafter, we suppose that the
space 


pology on the set 


We begin with a couple of auxiliary results.
Lemma 1. Let 



Proof. Since 








As a consequence, we have the following property.
Corollary 1. Let 

in 


Our next step concerns the study of topological properties of the feasible set 
The following result is crucial for our further analysis.
Theorem 3. Let 




Proof. By Theorem 1 and compactness properties of the space




Then by Lemma 1, we have
It remains to show that the limit pair 


In view of (34) and Lemma 1, we have
Moreover, due to the compactness of the Sobolev imbedding 

where Hölder’s inequality yields
We, thus, can pass to the limit in relation (35) as 

and the cone 

strongly in 
This fact together with 


In conclusion of this section, we give the existence result for optimal pairs to problem (19).
Theorem 4. Assume that, for given distributions



Proof. Since the set 


implies the existence of a constant 
Hence, in view of the definition of the class of admissible controls 
quence 









Thus,
Hence, 
5. Henig Relaxation of State-Constrainted OCP (19)
The main goal of this section is to provide a regularization of the pointwise state constraints by replacing the ordering cone 


We begin with some formal descriptions and abstract results. Let Z be a real normed space, and let 
Definition 5.1. A nonempty convex subset B of a nontrivial ordering cone 







In what follows, we always assume that the ordering cone 



In order to introduce a representation for a base of


and
we define the dual cone and the quasi-interior of the dual cone of

Due to Lemma 1.28 in [25] , we can give the following result.
Lemma 2. Let 








Remark 5.3. As follows from Lemma 2, the set

is a closed base of ordering cone
Now, we are prepared to introduce the definition of a so-called Henig dilating cone (see Zhuang, [24] ) which is based on the existence of a closed base of ordering cone
Definition 5.2. Let Z be a normed space, and let 

where 
It is clear that 


Proposition 5. Let Z be a normed space, and let 


the following statements hold true.
1) 

2)
3)
4)
5) the implication

holds true with
In the context of constraint qualifications problem, the following result plays an important role.
Proposition 6. Let Z be a normed space, and let 



holds true.
Proof. Let 



holds true. Let’s assume for a moment that

Then we obtain
which completes the proof. In order to show (40), let 
Then
yields
As a result, (40) is satisfied.
Remark 5.2. The following property, coming from Proposition 6, turns out rather useful: in order to prove

The following result shows that Henig dilating cones 
Proposition 7. Let 



ence such that


with respect to the norm topology of Z as k tends to infinity, that is
where
Proof. Let 






Taking into account the inclusion (41) and the fact that
we get

To show that the sequence 


However, the inclusion (43) is equivalent to

Let 





This implies
Combining (42), (43), and (44), we arrive at the relation
Thus, 
Taking these results into account, we associate with OCP (19) the following family of Henig relaxed pro- blems

subject to the constraints

or in a more compact form each of these problems can be stated as follows

where

the base B takes the form (36), and the feasible set 






Here, 
Since, by Proposition 6, the inclusion 


We are now in a position to show that using the relaxation approach we can reduce the main suppositions of Theorem 4. In particular, we can characterize Hypothesis (

Theorem 8. Let 

for given distributions


the Henig relaxed problem (47) has a nonempty set of feasible solutions 




then the sequence s
Proof. Since the implication 


Let 





it follows by compactness arguments (see the proof of Theorem 4) that there exist a subsequence of
(still denoted by the same index) and a pair 
Closely following the proof of Theorem 3, it can be shown that the limit pair 




It remains to establish the inclusions

By contraposition, let us assume that





by Proposition 7 and definition of the Kuratowski limit, it is easy to conclude the existence of an index 

However, in view of the strong convergence property (51), there is an index 

Combining (53) and (54), we finally obtain
This, however, is a contradiction to
Thus,


As an obvious consequence of this Theorem and Theorem 4, we have the following noteworthy property of the Henig relaxed problems (47).
Corollary 2. Let




The next result is crucial in this section. We show that some optimal solutions for the original OCP (19) can be attained by solving the corresponding Henig relaxed problems (45)-(46). However, we do not claim that the entire set of the solutions to OCP (19) can be restored in such way.
Theorem 9. Let








Then there is a subsequence 




Proof. In view of a priory estimate (24), the uniform boundedness of optimal controls with respect to BV-norm (55) implies the fulfilment of condition (50)2. Hence, the compactness property (56) and the inclusion 





On the other hand, by Proposition 5 (see property (4)), we have 



As a result, taking into account the relations (58) and (59), and the lower semicontinuity property of the cost functional I with respect to the t-convergence, we finally get
Thus,
and we arrive at the desired property (57)2. The proof is complete.
Remark 5.3. It is worth to note that condition (55) can be omitted if the original OCP (19) is regular, that is when Hypothesis (




Since, by Proposition 6, the inclusion 


it follows that
As a result, (60) leads to the estimate
As was mentioned at the beginning of this section, the main benefit of the relaxed optimal control problems (45)-(46) comes from the fact that the Henig dilating cone 




6. Variational Inequality Approach to Regularization of OCP (19)
As follows from Theorem 4, the existence of optimal solutions to the problem (19) can be obtained by using compactness arguments and the Hypothesis (

and the corresponding state constraint
with some accuracy. Here, the operator 



for each


Here, 
Note that the reverse statement is not true in general. In fact, we discuss a variant of the penalization approach, called the “variational inequality (VI) method”. This idea was first studied in [27] . Thus, if a pair 

equation

measure in an associated cost functional. As a result, we arrive at the following penalized OCP:

subject to the constraints

or in a more compact form this problem can be stated as follows

where 


In this section we show that penalized OCP (65) is solvable for each 


Lemma 3. Under assumptions (11)-(14), for every fixed 



Proof. Let 





where 



Thus, the operator 



As an obvious consequence of Lemma 3, we have the following noteworthy property of penalized OCP (63) - (64).
Corollary 3 For each 

To proceed further, we introduce the following notion.
Definition 6.1. An operator 





implies the relation

for all
Definition 6.2. We say that an operator 

any sequence 



imply the relation
Our next intention is to prove the following crucial result.
Theorem 10. The operator
Proof. Let 




and divide our proof onto several steps.
Step 1. We show that, for each

Indeed, since 





As for the first term in (70), we note that 





by the Lebesgue Dominated Theorem. Since the sequence 

it follows from (72) that 


Step 2. Let us show that

By analogy with the previous step, we note that 





This means that
But we also have that the sequence 


for each



by definition of the weak convergence in

In view of the subcritical growth condition (12), we have the following estimate
where 





Thus, to conclude this step, we have to show that


So, taking into account (76) and the fact that 


for all positive
Since the function 

Step 3. This is the final step of our proof. As follows from (69), for every element 


Let 



Since

Passing to the limit in (79) as

where
and
Hence, for each

Since the convergence 




As a result, we deduce from (81) and (82) that
that is, the inequality (68) is valid.
Remark 6.1. In fact (see [19] , Remark 3.13), we have the following implication:
Hence, in view of Theorem 10, we can claim that the operator

We are now in a position to show that the penalized optimal control problem in the coefficient of variational inequality (63)-(64) is solvable for each value
Lemma 4 If the assumptions (11)-(14) are valid, then the OCP (63)-(64) admits at least one solution 




Proof. Let 

immediately imply that the sequence 



On the other hand, from (83) it follows that
So, comparing these two chains of relations, we arrive at the existence of a constant 



Since
and the set 


as

Having put here
Hence,
by the quasi-monotonicity property of the operator A. Combining this inequality with (84), we come to the re- lation
Thus, 
Let us show that 




Substituting 

Since the quasi-monotone operator possesses the 

Thus, 
The next step of our analysis is to consider a sequence of optimal pairs 

Theorem 11. Let 


Then the sequence 



i.e. 
Proof. Let 

Hence, the following estimate for the optimal states takes place

Let us show that the sequence of corresponding optimal controls 



From this, we immediately conclude that
and estimate (88), we may assume that there exists a pair 




Let us show that the pair 



Thus, 


It remains to prove that 


Therefore, passing to the limit in this inequality as 
This contradiction immediately leads us to the conclusion: The 
Remark 6.2. As follows from the proof of Theorem 11, whatever the sequence of optimal solutions 

Remark 6.3. It is easy to see that in the case if the feasible set to the original OCP is nonempty, it suffices to guarantee the fulfilment of assumption (85). Indeed, let 



and we arrive at the inequality (85).
Acknowledgements
Research is funded by DFG-Excellence Cluster Engineering for Advanced Materials.
Cite this paper
Peter Kogut,Günter Leugering,Ralph Schiel, (2016) On Henig Regularization of Material Design Problems for Quasi-Linear p-Biharmonic Equation. Applied Mathematics,07,1547-1570. doi: 10.4236/am.2016.714134
References
- 1. Sprekels, J. and Tiba, D. (2005) Optimal Design of Mechanical Structures. In: Imanuvilov, O., Leugering, G., Triggiani, R., et al., Eds., Control Theory of Partial Differential Equations, Chapman& Hall/CRC, Boca Raton, 259-271.
- 2. Sprekels, J. and Tiba, D. (2009) Optimization Problems for Thin Elastic Structures. Number 158 in International Series of Numerical Mathematics, Birkhauser Verlag, Basel.
- 3. Neitaanmaki, P., Sprekles, J. and Tiba, D. (2006) Optimization of Elliptic Systems, Theory and Applications. Springer Monographs in Mathematics, Springer-Verlag, Berlin.
- 4. Gazzola, F., Grunau, H.-C. and Sweers, G. (2010) Polyharmonic Boundary Value Problems: Positivity Preserving and Nonlinear Higher Order Elliptic Equations in Bounded Domains. Springer-Verlag, Berlin.
http://dx.doi.org/10.1007/978-3-642-12245-3 - 5. Litvinov, V.G. (1987) Optimization in Elliptic Boundary Value Problems with Applications to Machanics. MIR, Moscow.
- 6. Lurie, K.A. (1993) Applied Optimal Control Theory of Distributed Systems. Planum Press.
- 7. Khan, A. and Sama, M. (2013) A New Conical Regularization for Some Optimization and Optimal Control Problems: Convergence Analysis and Finite Element Discretization. Numerical Functional Analysis and Optimization, 34, 861-895.
http://dx.doi.org/10.1080/01630563.2013.806546 - 8. Adams, R. (1975) Sobolv Spaces. Academic Press, Cambridge.
- 9. Lions, J.-L. and Magenes, E. (1968) Problemes aux limites non homogenes et applications. Vol. 1, Dunod, Paris.
- 10. Bucur, D. and Gazzola, F. (2011) The First Biharmonic Steklov Eigenvalue: Positivity Preserving and Shape Optimization. Milan Journal of Mathematics, 79, 247-258.
http://dx.doi.org/10.1007/s00032-011-0143-x - 11. Gilbarg, D. and Trudinger, N.S. (2001) Elliptic Partial Differential Equations of Second Order. Springer-Verlag, Berlin.
- 12. Colasuonno, F. and Pucci, P. (2001) Multiplicity of Solutions for p(x)-Polyharmonic Elliptic Kirchoff Equations. Nonlinear Analysis, Theory, Methods, and Applications, 74, 5962-5974.
http://dx.doi.org/10.1016/j.na.2011.05.073 - 13. Lubyshev, V.F. (2011) Multiple Solutions of an Even-Order Nonlinear Problems with Convex-Concave Nonlinearity. Nonlinear Analysis, Theory, Methods, and Applications, 74, 1345-1354.
http://dx.doi.org/10.1016/j.na.2010.10.007 - 14. Giusti, E. (1984) Minimal Surfaces and Functions of Bounded Variation. Birkhauser Verlag, Basel.
http://dx.doi.org/10.1007/978-1-4684-9486-0 - 15. Attouch, H., Buttazzo, G. and Michaille, G. (2006) Variational Analysis in Sobolev and BV Spaces: Application to PDE and Optimization. SIAM, Philadelphia.
- 16. Roubicek, T. (2013) Nonlinear Partial Differential Equations with Applications. Birkhauser Verlag, Basel.
http://dx.doi.org/10.1007/978-3-0348-0513-1 - 17. Lions, J.-L. (1969) Some Methods of Solving Non-Linear Boundary Value Problems. Dunod-Gauthier-Villars, Paris.
- 18. Casas, E. (1992) Optimal Control in the Coefficients of Elliptic Equations with State Constraints. Applied Mathematics & Optimization, 26, 21-37.
http://dx.doi.org/10.1007/BF01218394 - 19. Kogut, P.I. and Leugering, G. (2011) Optimal Control Problems for Partial Differential Equations on Reticulated Domains: Approximation and Asymptotic Analysis, Birkhauser Verlag, Basel.
http://dx.doi.org/10.1007/978-0-8176-8149-4 - 20. Roubicek, T. (1997) Relaxation in Optimization Theory and Variational Calculus. De Gruyter, Berlin, New York.
http://dx.doi.org/10.1515/9783110811919 - 21. Borwein, J.M. and Zhuang, D.M. (1993) Super Efficiency in Vector Optimization. Transactions of the American Mathematical Society, 338, 105-122.
http://dx.doi.org/10.1090/S0002-9947-1993-1098432-5 - 22. Henig, M.I. (1982) Proper Efficiency with Respect to Cones. Journal of Optimization Theory and Applications, 36, 387-407.
http://dx.doi.org/10.1007/BF00934353 - 23. Makarov, E.K. and Rackkovski, N.N. (1999) Unified Representation of Proper Efficiency by Means of Dilating Cones. Journal of Optimization Theory and Applications, 101, 141-165.
http://dx.doi.org/10.1023/A:1021775112119 - 24. Zhuang, D.M. (1994) Density Result for Proper Efficiencies. SIAM Journal on Control and Optimization, 32, 51-58.
http://dx.doi.org/10.1137/S0363012989171518 - 25. Jahn, J. (2004) Vector Optimization: Theory, Applications, and Extensions. Springer-Verlag, Berlin.
http://dx.doi.org/10.1007/978-3-540-24828-6 - 26. Schiel, R. (2014) Vector Optimization and Control with PDEs and Pointwise State Constraints. Phd Thesis, Friedrich-Alexander-University Erlangen-Nurnberg, Nurnberg.
- 27. Mel’nik, V.S. (1986) Method of Monotone Operators in the Theory of Constrained Optimal System. Rep. Ukrain. Acad. Sci., 7, 64-67.
- 28. Barbu, V. (1993) Analysis and Control of Infinite Dimensional Systems. Academic Press, Cambridge.
- 29. Ekeland, I. and Temam, R. (1976) Convex Analysis and Variational Problems. Elsevier, Amsterdam.












































































































