Next Article in Journal
A Note on the Topological Group c0
Next Article in Special Issue
A Two Dimensional Discrete Mollification Operator and the Numerical Solution of an Inverse Source Problem
Previous Article in Journal
Selective Survey on Spaces of Closed Subgroups of Topological Groups
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On the Shape Differentiability of Objectives: A Lagrangian Approach and the Brinkman Problem

by
José Rodrigo González Granada
1,
Joachim Gwinner
2 and
Victor A. Kovtunenko
3,*,†
1
Department of Mathematics, Universidad Tecnológica de Pereira, 660003 Pereira, Colombia
2
Fakultät für Luft- und Raumfahrttechnik, Institut für Mathematik und Rechneranwendung, Universität der Bundeswehr München, Werner-Heisenberg-Weg 39, 85577 Neubiberg, Germany
3
Lavrentyev Institute of Hydrodynamics, Siberian Division of the Russian Academy of Sciences, Novosibirsk 630090, Russia
*
Author to whom correspondence should be addressed.
Current address: Institute for Mathematics and Scientific Computing, Karl-Franzens University of Graz, NAWI Graz, Heinrichstr, 36, 8010 Graz, Austria.
Axioms 2018, 7(4), 76; https://doi.org/10.3390/axioms7040076
Submission received: 26 September 2018 / Revised: 23 October 2018 / Accepted: 25 October 2018 / Published: 27 October 2018
(This article belongs to the Special Issue Applications of Differential Equations and Dynamical Systems)

Abstract

:
This paper establishes the shape derivative of geometry-dependent objective functions for use in constrained variational problems. Using a Lagrangian approach, our differentiablity result is based on the theorem of Delfour–Zolésio on directional derivatives with respect to a parameter of shape perturbation. As the key issue of the paper, we analyze the bijection under the kinematic transport of geometries that is needed for function spaces and feasible sets involved in variational problems. Our abstract theoretical result is applied to the Brinkman flow problem under incompressibility and mixed Dirichlet–Neumann boundary conditions, and provides an analytic formula of the shape derivative based on the velocity method.

1. Introduction

Within an abstract theoretical framework, this paper is devoted to the shape differentiablity of geometry-dependent objective functions as they are considered in constrained variational problems.
Constraints arise in a variety of applications. The constraint operator (see (3)) may become: the trace operator under contact conditions [1,2,3], the jump operator for cracks and anticracks [4,5,6], the gradient operator in plasticity [7], the divergence operator under incompressibility conditions [8,9,10], a state-constraint in mathematical programs with equilibrium constraints [11,12], and the like. The constraint problems are related to parameter identification problems (see the theory in References [13,14,15] and application to biological systems in Reference [16]), to inverse problems by the mean of observation data used in mathematical physics [17,18] and in acoustics [19,20,21], to overdetermined and free-boundary problems [22,23]. As an application, in the current paper we focus on the incompressible Brinkman flow problem under a divergence-free condition (see the related modeling of porous medium in References [24,25], well-posedness analysis in Reference [26], and fluid–porous coupling with numerics in References [27,28]).
For a general theory of (nonlinear) optimization in infinite dimensions we refer to References [7,29,30,31,32] for methods of shape optimization and, in particular, shape derivatives. Specifically in fracture mechanics, the shape derivative of energy functionals with respect to a crack extension is called the energy release rate and is of primary importance to engineers as a fracture criterion [33]. In References [4,7], a nonlinear theory of cracks subject to contact conditions and their shape derivatives was developed. Shape perturbation needs a bijective property for function spaces and feasible sets, which, however, fails for curvilinear cracks being in contact. To overcome this obstacle, References [5,6] suggest a Γ -type convergence, and References [34,35] propose the use of primal-dual Lagrangian reformulation of the crack problems. In incompressible fluid dynamics, to treat the divergence-free condition, divergence-preserving Piola transformation was employed in Reference [10], but lacking a mathematical foundation. The Lagrangian approach to shape optimization was developed further for abstract quadratic objective functions, and the direct proof of shape differentiability was given in Reference [36]. In Reference [9], we applied the primal-dual technique to the Stokes problem under a divergence-free constraint.
In the current work, we extend the Lagrangian approach to abstract geometry-dependent objective functions based on the theorem of Delfour–Zolésio. Thus, we provide the directional derivative with respect to a parameter of shape perturbation in a general setting. We stress that, while bijection fails for the primal cone (a feasible set of primal variables), we obtained the shape derivative in virtue of the bijection property of the dual cone (a feasible set of dual variables). In Section 2, we set a geometry-dependent constrained optimization problem in abstract form, and its shape derivative is derived in Section 3. We applied our theory to a Brinkman flow problem subject to the divergence-free condition in Section 4 and provide an analytic formula of the shape derivative based on the velocity method.

2. Geometry-Dependent CO Problem

For real parameter t I : = ( t 0 , t 1 ) with t 0 < t 1 , we introduce geometry by a kinematic flow
t Ω t R d , d N .
For every t I , let V ( Ω t ) and H ( Ω t ) be Hausdorff topological spaces with dual spaces V ( Ω t ) and H ( Ω t ) . Using the order relation for measurable functions in H ( Ω t ) , we specify a feasible set provided by the inequality constraint:
K ( Ω t ) : = { w V ( Ω t ) | G ( Ω t ) w 0 }
which is a convex closed cone determined by a linear continuous constraint operator
w G ( Ω t ) w : V ( Ω t ) H ( Ω t ) .
Extension of Function (3) to a nonlinear operator can be found in Reference [23]. For a geometry-dependent objective given by a continuous and generally nonlinear function
w F ( w ; Ω t ) : V ( Ω t ) R ¯ ,
we consider the constrained optimization (CO) problem: find u t K ( Ω t ) such that
F ( u t ; Ω t ) = inf w K ( Ω t ) F ( w ; Ω t ) = : f ( t ) ,
where optimal value (OV) function f : I R ¯ . The corresponding solution set implies:
F t : = { u K ( Ω t ) | F ( u ; Ω t ) = f ( t ) }
and may be empty, with exactly one (called a singleton) or more elements.
Proposition 1.
(i) 
Solution set F t (6) for CO problem (5) is nonempty under the assumption:
(F1) 
there exists a minimizing sequence u k K ( Ω t ) , such that F ( u k ; Ω t ) f ( t ) as k and an accumulation point u t K ( Ω t ) yielding the convergence
u k u t weakly in V ( Ω t ) as k
and the lower estimate:
lim   inf k F ( u k ; Ω t ) F ( u t ; Ω t ) .
(ii) 
If the next assumption holds:
(F2) 
objective function u t F ( · ; Ω t ) is Gâteaux-differentiable, that is, the following limit exists:
lim s 0 F ( u t + s w ; Ω t ) F ( u t ; Ω t ) s : = u F ( u t ; Ω t ) , w Ω t , w V ( Ω t )
with u F ( u t ; Ω t ) V ( Ω t ) and duality pairing · , · Ω t between V ( Ω t ) and V ( Ω t ) ,
then, a first-order optimality condition for CO problem (5) necessitates variational inequality (VI) stated in the form of two relations:
u F ( u t ; Ω t ) , u t Ω t = 0 , u F ( u t ; Ω t ) , w Ω t 0 for all w K ( Ω t ) .
Conversely, VI (10) is also sufficient for Problem (5) when
(F3) 
objective function w F ( · ; Ω t ) in (4) is convex.
(iii) 
Solution set F t = { u t } in (6) is a singleton if the condition (F3) is strict:
(F3’) 
the objective function w F ( · ; Ω t ) in (4) is strictly convex.
Proof. 
In a reflexive Banach space V ( Ω t ) , Condition (7) holds for coercive functions, and Condition (8) is satisfied for the weakly lower semicontinuous functions w F ( · ; Ω t ) in (4). The corresponding proof of assertions (i), (ii), and (iii) is standard, see e.g., (Reference [4] Theorem 1.11). □
For objective OV function f, we aimed at its shape derivative (SD):
f ( t ) : = lim s 0 + f ( t + s ) f ( t ) s ( one-sided ) ,
where, according to Problem (5), the perturbed OV is determined over perturbed geometry Ω t + s :
f ( t + s ) = inf v K ( Ω t + s ) F ( v ; Ω t + s ) = F ( u t + s ; Ω t + s ) .
The shape sensitivity of VI was investigated in Reference [32], Chapter 4. However, it can be advantageous to consider a Lagrangian formulation of the problem, see Remark 1.
Let us define the Lagrangian function as ( w , p ) L ( w , p ; Ω t ) : V ( Ω t ) × H ( Ω t ) R ¯ ,
L ( w , p ; Ω t ) : = F ( w ; Ω t ) ( p , G ( Ω t ) w ) Ω t
with duality pairing ( · , · ) Ω t between H ( Ω t ) and H ( Ω t ) . For a convex closed cone that is dual to the primal cone (2) (hence called the dual cone):
K ( Ω t ) : = { p H ( Ω t ) | ( p , G ( Ω t ) w ) Ω t 0 for all w K ( Ω t ) } ,
a saddle-point (SP) problem related to Problem (5) reads: find a pair ( u t , λ t ) V ( Ω t ) × K ( Ω t ) , such that
L ( u t , p ; Ω t ) L ( u t , λ t ; Ω t ) L ( w , λ t ; Ω t ) for all ( w , p ) V ( Ω t ) × K ( Ω t ) .
The Lagrangian OV function : I R ¯ is defined from the relations
( t ) = l t : = sup p K ( Ω t ) inf w V ( Ω t ) L ( w , p ; Ω t ) inf w V ( Ω t ) sup p K ( Ω t ) L ( w , p ; Ω t ) = : l t ,
and the corresponding solution sets are
V t = { u V ( Ω t ) | sup p K ( Ω t ) L ( u , p ; Ω t ) = l t } , K t = { λ K ( Ω t ) | inf w V ( Ω t ) L ( w , λ ; Ω t ) = l t } .
Using Notations (16) and (17), Inequalities (15) determine the set of SP for the Lagrangian (13):
S ( t ) : = { ( u t , λ t ) V t × K t , ( t ) = l t = L ( u t , λ t ; Ω t ) = l t } ,
which may be empty, a singleton, or contain more than one pair ( u t , λ t ) .
Proposition 2.
(i) 
Under the following assumption:
(G1) 
mapping V ( Ω t ) H ( Ω t ) , w G ( Ω t ) w in Problem (3) is surjective,
that is, for every ζ H ( Ω t ) there is at least one w V ( Ω t ) , such that G ( Ω t ) w = ζ , the dual cone (14) can be restated equivalently as the cone of non-negative elements in the dual space:
K ( Ω t ) = { p H ( Ω t ) | ( p , ζ ) Ω t 0 for all ζ H ( Ω t ) , ζ 0 } .
(ii) 
Under Assumptions (G1) and (F1), (F2), (F3), the set of SP in (18) is nonempty, and saddle-point ( u t , λ t ) S ( t ) satisfies the optimality system (OS) in the form of variational relations:
u F ( u t ; Ω t ) , w Ω t ( λ t , G ( Ω t ) w ) Ω t = 0 for all w V ( Ω t ) ,
( λ t , G ( Ω t ) u t ) Ω t = 0 , ( p , G ( Ω t ) u t ) Ω t 0 for all p K ( Ω t ) .
Primal component u t is unique when (F3’) holds. Uniqueness of dual component λ t takes place if
(G2) 
the Ladyzhenskaya–Babuška–Brezzi (LBB) condition holds: there exists c ̲ > 0 , such that
sup w V ( Ω t ) , w 0 ( p , G ( Ω t ) w ) Ω t w V ( Ω t ) c ̲ p H ( Ω t ) , p H ( Ω t ) .
In this case, sets V t = { u t } and K t = { λ t } in Problem (17) are singletons.
(iii) 
If a saddle-point in Problem (15) exists, then, under assumption (G1), the primal component u t K ( Ω t ) solves the CO problem (5), and the OV functions for objective (5) and for Lagrangian (16) coincide:
f ( t ) = ( t ) .
Proof. 
Assertion (i) is a consequence of the bipolar theorem. In Assertion (ii), using the Proposition 1 solution to the SP problem (15) can be derived in a standard way by determining u t from VI (10) and by setting Lagrange multiplier λ t from (20) (e.g., Reference [30] Chapter III, Proposition 3.1). Optimality Conditions (20) and (21) and the uniqueness assertion under LBB Condition (G2) are stated, for example, in Reference [3], Theorem 3.14. In Assertion (iii), Problem (5) follows straightforwardly by excluding λ t from OS (20) and (21), while Identity (23) is guaranteed by ( λ t , G ( Ω t ) u t ) Ω t = 0 in L ( u t , λ t ; Ω t ) . □
Based on Identity (23), the shape derivative of Objective (11) can be substituted by the SD of Lagrangian as follows:
f ( t ) = ( t ) : = lim s 0 + ( t + s ) ( t ) s ( one-sided ) .
According to Definition (16), the perturbed OV ( t + s ) = l t + s = sup μ K ( Ω t + s ) inf v V ( Ω t + s ) L ( v , μ ; Ω t + s ) in (24) is obtained for the perturbed Lagrangian ( v , μ ) L ( v , μ ; Ω t + s ) : V ( Ω t + s ) × H ( Ω t + s ) R ¯ ,
L ( v , μ ; Ω t + s ) = F ( v ; Ω t + s ) ( μ , G ( Ω t + s ) v ) Ω t + s ,
and the perturbed SP problem reads: find a pair ( u t + s , λ t + s ) V ( Ω t + s ) × K ( Ω t + s ) such that
L ( u t + s , μ ; Ω t + s ) L ( u t + s , λ t + s ; Ω t + s ) L ( v , λ t + s ; Ω t + s ) for all ( v , μ ) V ( Ω t + s ) × K ( Ω t + s ) .
However, the theorem of Delfour and Zolésio on differentiability of Lagrangians with respect to parameter t (Reference [31] Chapter 10, Theorem 5.1), which they call the Correa–Seeger theorem, is inapplicable to prove the limit in Problem (24). The reason is that the Lagrangian (25) is determined over geometry-dependent spaces, which, in turn, depend on s 0 + .

3. Shape Differentiability of Objectives for CO Problems

For the above reason, we further parametrized the geometry in Problem (1) as follows: For fixed t I let
( s , x ) ϕ s ( x ) , ( s , y ) ϕ s 1 ( y ) W 1 , ( I t ; W loc 1 , ( R d ) ) d
associate coordinate transformation (CT) y = ϕ s ( x ) and its inverse x = ϕ s 1 ( y ) :
( ϕ s 1 ϕ s ) ( x ) = x , ( ϕ s ϕ s 1 ) ( y ) = y
such that shape perturbation
Ω t + s = { y R d | y = ϕ s ( x ) , x Ω t }
builds the diffeomorphism
ϕ s : Ω t Ω t + s , x y , ϕ s 1 : Ω t + s Ω t , y x .
Following the velocity method (e.g., Reference [36]), a given kinematic velocity
Λ ( t , x ) L ( I ; W loc 1 , ( R d ) ) d
establishes Flow (27) by solving the nonautonomous ODE and the transport equation [35]:
d d s ϕ s = Λ ( t + s , ϕ s ) for s > 0 ,     ϕ s = x as s = 0 ;
s ϕ s 1 + Λ ( t + s , y ) · y ϕ s 1 = 0 for s > 0 , ϕ s 1 = y as s = 0 .
Assume that
(D1)
map ( v , μ ) ( v ϕ s , μ ϕ s ) is bijective between the function spaces
V ( Ω t + s ) V ( Ω t ) , H ( Ω t + s ) H ( Ω t ) ;
(D2)
map μ μ ϕ s is bijective between the dual cones
K ( Ω t + s ) K ( Ω t ) .
Assumption (D1) determines the transformed perturbed Lagrangian: I t × V ( Ω t ) × H ( Ω t ) R ¯ ,
( s , w , p ) L ˜ ( s , w , p ; Ω t ) : = F ˜ ( s , w ; Ω t ) ( p , G ˜ ( s ; Ω t ) w ) Ω t
with continuous function ( s , w ) F ˜ ( s , w ; Ω t ) : I t × V ( Ω t ) R ¯ and continuous operator ( s , w ) G ˜ ( s ; Ω t ) w : I t × V ( Ω t ) H ( Ω t ) , which are resulted from the application of CT (28):
F ˜ ( s , v ϕ s ; Ω t ) = F ( v ; Ω t + s ) , ( μ ϕ s , G ˜ ( s ; Ω t ) ( v ϕ s ) ) Ω t = ( μ , G ( Ω t + s ) v ) Ω t + s
in such a way that the identity holds:
L ˜ ( s , v ϕ s , μ ϕ s ; Ω t ) = L ( v , μ ; Ω t + s ) , ( v , μ ) V ( Ω t + s ) × H ( Ω t + s ) .
According to Problem (37), F ˜ ( 0 , w ; Ω t ) = F ( w ; Ω t ) and G ˜ ( 0 ; Ω t ) = G ( Ω t ) as s = 0 .
Remark 1.
We remark that Assumption (D2) on the dual cone stated in the (19) form is advantageous in comparison to the condition on Feasible Set (2):
(D2’) 
map v v ϕ s is bijective between the feasible sets
K ( Ω t + s ) K ( Ω t )
that is needed for OV perturbation for the objective in Problem (12).
In fact, Bijection (D2’) fails for integral and derivative-type operators in Feasible Set (2), which are generally not preserved G ( Ω t + s ) ϕ s G ( Ω t ) under velocity-induced geometry flow ϕ s (see the example in Remark 2).
Under Assumptions (D1) and (D2), we reset the OV function for the transformed perturbed Lagrangian as I t R ¯ , s ˜ ( · ; t ) :
˜ ( s ; t ) = l ˜ t + s : = sup p K ( Ω t ) inf w V ( Ω t ) L ˜ ( s , w , p ; Ω t ) inf w V ( Ω t ) sup p K ( Ω t ) L ˜ ( s , w , p ; Ω t ) = : l ˜ l + s
and the corresponding solution sets:
V ˜ t + s = { u V ( Ω t ) | sup p K ( Ω t ) L ˜ ( s , u , p ; Ω t ) = l ˜ t + s } , K ˜ t + s = { λ K ( Ω t ) | inf w V ( Ω t ) L ˜ ( s , w , λ ; Ω t ) = l ˜ t + s } .
If s = 0 , then l ˜ t = l t , l ˜ l = l l in (16), (40) and V ˜ t = V t , K ˜ t = K t for sets (17) and (41) according to (38). Since a solution to Perturbed SP Problem (26) exists, applying Identity (38), we get the solution
u ˜ t + s : = u t + s ϕ s V ( Ω t ) , λ ˜ t + s : = λ t + s ϕ s K ( Ω t )
satisfying the transformed perturbed SP inequalities:
L ˜ ( s , u ˜ t + s , p ; Ω t ) L ˜ ( s , u ˜ t + s , λ ˜ t + s ; Ω t ) L ˜ ( s , w , λ ˜ t + s ; Ω t ) for all ( w , p ) V ( Ω t ) × K ( Ω t ) ,
thus, it forms the SP set for transformed perturbed Lagrangian (36) similar to (18):
S ˜ ( t + s ) : = { ( u ˜ t + s , λ ˜ t + s ) V ˜ t + s × K ˜ t + s | l ˜ t + s = L ˜ ( s , u ˜ t + s , λ ˜ t + s ; Ω t ) = l ˜ t + s } .
Following Delfour–Zolésio, we assume that a small δ > 0 and a topology in V ( Ω t ) × H ( Ω t ) exist, such that the following hypotheses resulting from the specific representation (36) hold:
(H1)
for all s [ 0 , δ ) , the set of saddle-points S ˜ ( t + s ) in (44) is nonempty;
(H2)
for all u s [ 0 , δ ) V ˜ t + s × K ˜ t and τ [ 0 , δ ) , there exists a one-sided partial derivative of the transformed perturbed Lagrangian at perturbation parameter s = τ in the form:
s L ˜ ( τ , u , p ; Ω t ) = s F ˜ ( τ , u ; Ω t ) p , s G ˜ ( τ ; Ω t ) u Ω t , p H ( Ω t ) ,
obtained from the partial derivatives of the objective and constraint functions as follows
s F ˜ ( τ , u ; Ω t ) : = lim s 0 + F ˜ ( τ + s , u ; Ω t ) F ˜ ( τ , u ; Ω t ) s , s G ˜ ( τ ; Ω t ) : = lim s 0 + G ˜ ( τ + s ; Ω t ) G ˜ ( τ ; Ω t ) s ;
(H3)
as s 0 + , accumulation point u ˜ t V ˜ t and subsequence u ˜ t + s k V ˜ t + s k exist, such that
u ˜ t + s k u ˜ t strongly in V ( Ω t ) as s k 0 +
and the lower estimate holds:
lim   inf τ , s k 0 + s F ˜ ( τ , u ˜ t + s k ; Ω t ) s F ˜ ( 0 , u ˜ t ; Ω t ) ;
(H4)
as s 0 + , accumulation point λ ˜ t K ˜ t and subsequence λ ˜ t + s k K ˜ t + s k exist, such that
λ ˜ t + s k λ ˜ t strongly in H ( Ω t ) as s k 0 + .
The main result is stated in the next theorem.
Theorem 1.
(i) 
Under Hypotheses (H1)–(H4), accumulation point ( u ˜ t , λ ˜ t ) of Sequences (47) and (49) implies an SP of the partial derivative s L ˜ ( 0 , u , λ ; Ω t ) on V ˜ t × K ˜ t , such that
sup λ K ˜ t inf u V ˜ t s L ˜ ( 0 , u , λ ; Ω t ) = s L ˜ ( 0 , u ˜ t , λ ˜ t ; Ω t ) = inf u V ˜ t sup λ K ˜ t s L ˜ ( 0 , u , λ ; Ω t ) ;
the shape derivative exists and is represented by the partial derivative at the accumulation point:
s ˜ ( 0 ; t ) : = lim s 0 + ˜ ( s ; t ) ˜ ( 0 ; t ) s = s L ˜ ( 0 , u ˜ t , λ ˜ t ; Ω t ) ( one-sided ) .
(ii) 
Under Assumptions (D1) and (D2), the shape derivatives defined in (24) and (51) coincide:
f ( t ) = ( t ) = s ˜ ( 0 ; t ) .
Proof. 
(i) According to Hypothesis (H1), SP inequalities (43) hold for sequence s k 0 + from Hypotheses (H3) and (H4). In the next section, we used the accumulation point of Sequences (47) and (49) and insert it for a test function ( w , p ) = ( u ˜ t , λ ˜ t ) V ˜ t × K ˜ t V ( Ω t ) × K ( Ω t ) into Problem (43) for s = s k :
L ˜ ( s k , u ˜ t + s k , λ ˜ t ; Ω t ) L ˜ ( s k , u ˜ t + s k , λ ˜ t + s k ; Ω t ) L ˜ ( s k , u ˜ t , λ ˜ t + s k ; Ω t ) .
We also plugged ( w , p ) = ( u ˜ t + s k , λ ˜ t + s k ) V ˜ t + s k × K ˜ t + s k V ( Ω t ) × K ( Ω t ) into (43) as s = 0 :
L ˜ ( 0 , u ˜ t , λ ˜ t + s k ; Ω t ) L ˜ ( 0 , u ˜ t , λ ˜ t ; Ω t ) L ˜ ( 0 , u ˜ t + s k , λ ˜ t ; Ω t ) .
Subtracting L ˜ ( 0 , u ˜ t , λ ˜ t ; Ω t ) from left inequality (53), using right inequality (54) and the representation of L ˜ in (36), after division of the result with s k leads to
L ˜ ( s k , u ˜ t + s k , λ ˜ t + s k ; Ω t ) L ˜ ( 0 , u ˜ t , λ ˜ t ; Ω t ) s k L ˜ ( s k , u ˜ t + s k , λ ˜ t ; Ω t ) L ˜ ( 0 , u ˜ t + s k , λ ˜ t ; Ω t ) s k = F ˜ ( s k , u ˜ t + s k ; Ω t ) F ˜ ( 0 , u ˜ t + s k ; Ω t ) s k       λ ˜ t , G ˜ ( s k ; Ω t ) G ˜ ( 0 ; Ω t ) s k u ˜ t + s k Ω t = s F ˜ ( α ̲ k s k , u ˜ t + s k ; Ω t ) ( λ ˜ t , s G ˜ ( β ̲ k s k ; Ω t ) u ˜ t + s k ) Ω t ,
where mean value theorem guarantees existence of weights α ̲ k , β ̲ k ( 0 , 1 ) in (55) due to differentiability property (46). Passing to the limit as s k 0 + by the virtue of assumptions (47) and (48) proceeds (55) further with the lower estimate:
lim   inf s k 0 + L ˜ ( s k , u ˜ t + s k , λ ˜ t + s k ; Ω t ) L ˜ ( 0 , u ˜ t , λ ˜ t ; Ω t ) s k s F ˜ ( 0 , u ˜ t ; Ω t ) ( λ ˜ t , s G ˜ ( 0 ; Ω t ) u ˜ t ) Ω t = s L ˜ ( 0 , u ˜ t , λ ˜ t ; Ω t ) .
Similarly, subtracting L ˜ ( 0 , u ˜ t , λ ˜ t ; Ω t ) from right inequality (53) and using left inequality (54) provides the following relations with weights α ¯ k , β ¯ k ( 0 , 1 ) :
L ˜ ( s k , u ˜ t + s k , λ ˜ t + s k ; Ω t ) L ˜ ( 0 , u ˜ t , λ ˜ t ; Ω t ) s k L ˜ ( s k , u ˜ t , λ ˜ t + s k ; Ω t ) L ˜ ( 0 , u ˜ t , λ ˜ t + s k ; Ω t ) s k = F ˜ ( s k , u ˜ t ; Ω t ) F ˜ ( 0 , u ˜ t ; Ω t ) s k       λ ˜ t + s k , G ˜ ( s k ; Ω t ) G ˜ ( 0 ; Ω t ) s k u ˜ t Ω t = s F ˜ ( α ¯ k s k , u ˜ t ; Ω t ) ( λ ˜ t + s k , s G ˜ ( β ¯ k s k ; Ω t ) u ˜ t ) Ω t ,
hence the upper bound
lim   sup s k 0 + L ˜ ( s k , u ˜ t + s k , λ ˜ t + s k ; Ω t ) L ˜ ( 0 , u ˜ t , λ ˜ t ; Ω t ) s k s L ˜ ( 0 , u ˜ t , λ ˜ t ; Ω t ) .
Inequalities (56) and (58) together imply Equality (51) yielding the SD s ˜ of OV function ˜ given in Problem (40) for transformed Lagrangian L ˜ , thus proving Assertion (i).
(ii) When Bijection (34) and (35) holds, from Identity (38) we infer that
˜ ( s k ; t ) ˜ ( 0 ; t ) s k = L ˜ ( s k , u ˜ t + s k , λ ˜ t + s k ; Ω t ) L ˜ ( 0 , u ˜ t , λ ˜ t ; Ω t ) s k = L ( u t + s k , λ t + s k ; Ω t + s k ) L ( u t , λ t ; Ω t ) s k = ( t + s k ) ( t ) s k .
After passage (59) to the limit as s k 0 + , this follows the equivalence asserted in (52) between s ˜ and the SD = f from (24). □

4. Application to Brinkman Flow

Let domain Ω t in Problem (1) have Lipschitz continuous boundary Ω t obeying the unit normal vector n t = ( n 1 t , , n d t ) outward to Ω t and consisting of two nonempty, disjoint sets Γ t D and Γ t N . For stationary force f = ( f 1 , , f d ) H loc 1 ( R d ) d , we consider the Brinkman problem [24]: find a flow velocity u t = ( ( u t ) 1 , , ( u t ) d ) and a pressure λ t satisfying
μ Δ u t + α u t + λ t = f in Ω t ,
div u t = 0 in Ω t ,
u t = 0 on Γ t D , μ n t u t + λ t n t = 0 on Γ t N ,
under mixed Dirichlet–Neumann boundary conditions (62), where parameter μ > 0 denotes fluid viscosity, and α > 0 is the drag coefficient.
Incompressibility Condition (61) determines the operator
G ( Ω t ) = div : H Γ 1 ( Ω t ) d L 2 ( Ω t )
where L 2 ( Ω t ) = : H ( Ω t ) = H ( Ω t ) , and Sobolev space H Γ 1 ( Ω t ) d = : V ( Ω t ) accounts for the Dirichlet boundary condition in (62):
H Γ 1 ( Ω t ) d = { w = ( w 1 , , w d ) H 1 ( Ω t ) d | w = 0 a . e . Γ t D } .
The primal cone (2) is presented here by the equality-type constraint:
K ( Ω t ) = { w H Γ 1 ( Ω t ) d | div w = 0 a . e . Ω t } .
Corresponding duality pairings are
u , w Ω t = Ω t ( μ u · w + α u · w ) d x , u , w H Γ 1 ( Ω t ) d ,
( p , G ( Ω t ) w ) Ω t = Ω t p div w d x , p L 2 ( Ω t ) ,
and the latter one builds the dual cone
K ( Ω t ) = { p L 2 ( Ω t ) | Ω t p div w d x = 0 for all w K ( Ω t ) } .
The underlying objective function in (5) reads:
F ( w ; Ω t ) = Ω t μ 2 | w | 2 + α 2 | w | 2 f · w d x .
The quadratic function in (69) satisfies Assumptions (F1), (F2), and (F3’) of Proposition 1; hence, CO Problem (5) obeys unique solution u t K ( Ω t ) . The Lagrangian function in (13) takes the form:
L ( w , p ; Ω t ) = Ω t μ 2 | w | 2 + α 2 | w | 2 f · w p div w d x .
If meas ( Γ t N ) > 0 , then LBB condition (22) holds, (Reference [3], Theorem 7.2), namely:
sup w H Γ 1 ( Ω t ) d , w 0 Ω t p div w d x w H 1 ( Ω t ) d c ̲ p L 2 ( Ω t ) , c ̲ > 0 , p L 2 ( Ω t ) .
This means that the divergence operator in (63) is surjective; thus, Assumptions (G1), (G2) and Assertions (i)–(iii) of Proposition 2 hold true. In the case of equality constraint, the positive cone turns into the whole space K ( Ω t ) = L 2 ( Ω t ) according to (19). OS (20) and (21) implies the solution pair ( u t , λ t ) H Γ 1 ( Ω t ) d × L 2 ( Ω t ) satisfying the following variational equations:
Ω t μ u t · w + ( α u t f ) · w λ t div w d x = 0 for all w H Γ 1 ( Ω t ) d ,
Ω t p div u t d x = 0 for all p L 2 ( Ω t ) .
Note that (73) is equivalent to (61) by the fundamental lemma of the calculus of variations, and any smooth solution of (72) after integration by parts leads to Brinkman Equation (60) and the Neumann condition in (62).
Next we apply CT (27). It can be directly checked that bijective Properties (34) and (35), and hence Assumptions (D1) and (D2), hold true. Possible counterexamples are accounted below.
Remark 2.
Bijective property (39) in (D2’) fails for the primal cone defined in (65) because condition div y v = ( ϕ s T ϕ s ) · x ( v ϕ s ) = 0 is not equivalent to div x ( v ϕ s ) = 0 for v H Γ 1 ( Ω t + s ) d .
Remark 3.
If meas ( Γ t N ) = 0 , then the operator in (63) is not surjective since the divergence operator maps div : H 0 1 ( Ω t ) d L 0 2 ( Ω t ) , where
L 0 2 ( Ω t ) = { p L 2 ( Ω t ) | Ω t p d x = 0 } ,
and its topologically dual space L 0 2 ( Ω t ) excludes constants. In this case, Bijection (35) in (D2) fails, because Ω t + s μ d y = Ω t ( μ ϕ s ) det ( ϕ s ) d x = 0 is not equivalent to Ω t ( μ ϕ s ) d x = 0 for μ L 0 2 ( Ω t + s ) .
Based on CT y = ϕ s ( x ) with Jacobian matrix ϕ s and its determinant det ( ϕ s ) > 0 , from (70) we get the transformed perturbed objective and Lagrangian according to Relations (36)–(38):
L ˜ ( s , w , p ; Ω t ) = F ˜ ( s , w ; Ω t ) Ω t p ( ϕ s T ϕ s ) · w det ( ϕ s ) d x ,         F ˜ ( s , w ; Ω t ) = Ω t μ 2 | ( ϕ s T ϕ s ) w | 2 + α 2 | w | 2 ( f ϕ s ) · w det ( ϕ s ) d x .
By checking Hypotheses (H1)–(H4), we prove the result on shape differentiability in the following:
Theorem 2.
In the Brinkman problem, there exists the shape derivative expressed by the partial derivative:
f ( t ) = ( t ) = s ˜ ( 0 ; t ) = s L ˜ ( 0 , u t , λ t ; Ω t ) = Ω t ( μ 2 div Λ | u t | 2 μ u t · ( Λ u t )    + α 2 div Λ | u t | 2 ( ( div Λ ) f + ( Λ · ) f ) · u t + λ t Λ T · u t ) d x ,
where ( u t , λ t ) H Γ 1 ( Ω t ) d × L 2 ( Ω t ) is the unique SP solving OS (72) and (73), and kinematic velocity Λ ( t , x ) comes from (31).
Proof. 
(H1) For s [ 0 , δ ) with arbitrarily fixed δ > 0 , the set of SP (44) contains exactly one element:
S ˜ ( t + s ) = { ( u ˜ t + s , λ ˜ t + s ) : = ( u t + s ϕ s , λ t + s ϕ s ) H Γ 1 ( Ω t ) d × L 2 ( Ω t ) }
by transforming solution ( u t + s , λ t + s ) H Γ 1 ( Ω t + s ) d × L 2 ( Ω t + s ) to the perturbed at t + s system (72) and (73). This satisfies the transformed perturbed OS for all ( w , p ) H Γ 1 ( Ω t ) d × L 2 ( Ω t ) :
Ω t ( μ ( ϕ s T ϕ s ) u ˜ t + s · ( ϕ s T ϕ s ) w + ( α u ˜ t + s f ϕ s ) · w λ ˜ t + s ( ϕ s T ϕ s ) · w ) det ( ϕ s ) d x = 0 ,
Ω t p ( ϕ s T ϕ s ) · u ˜ t + s det ( ϕ s ) d x = 0 .
(H2) For the velocity Λ from (31), we expand as s 0 + at ( t , x ) (see Reference [37] Chapter 5):
ϕ s 1 ϕ s = I s Λ + o ( s ) , det ( ϕ s ) = 1 + s div Λ + o ( s ) , f ϕ s = f + s Λ · f + o ( s )
uniformly over t and Ω t . Plugging (79) into (74) we straightforwardly derived the partial derivative s L ˜ ( 0 , w , p ; Ω t ) and extended it to s L ˜ ( τ , w , p ; Ω t ) in (45) and (46) for τ [ 0 , δ ) by setting the time-shifted velocity Λ | t + τ : = Λ ( t + τ , x ) instead of Λ = Λ ( t , x ) . It has the following form:
s L ˜ ( τ , w , p ; Ω t ) = s F ˜ ( τ , w ; Ω t ) Ω t p ( div Λ | t + τ div w Λ | t + τ T · w ) d x , s F ˜ ( τ , w ; Ω t ) = Ω t ( μ 2 div Λ | t + τ | w | 2 μ w · ( Λ | t + τ w ) + α 2 div Λ | t + τ | w | 2 ( ( div Λ | t + τ ) f + ( Λ | t + τ · ) f ) · w ) d x .
(H3) & (H4) Testing (77) with w = u ˜ t + s , since ( ( ϕ s T ϕ s ) · u ˜ t + s ) det ( ϕ s ) = 0 due to (78), we have
Ω t μ | ( ϕ s T ϕ s ) u ˜ t + s | 2 + α | u ˜ t + s | 2 det ( ϕ s ) d x = Ω t ( f ϕ s ) · u ˜ t + s det ( ϕ s ) d x .
Using the asymptotic expansions (79) and applying Young’s inequality with a suitable weight to the right-hand side of (81), this follows the uniform in s [ 0 , δ ) estimate:
u ˜ t + s H 1 ( Ω t ) d 2 1 min ( μ , α ) Ω t μ | u ˜ t + s | 2 + α | u ˜ t + s | 2 d x = O ( 1 ) .
We divide (77) with the norm of w 0 and apply the Cauchy–Schwarz inequality such that
Ω t λ ˜ t + s ( ϕ s T ϕ s ) · w det ( ϕ s ) d x w H 1 ( Ω t ) d = Ω t μ ( ( ϕ s T ϕ s ) u ˜ t + s ) · ( ( ϕ s T ϕ s ) w ) + ( α u ˜ t + s f ϕ s ) · w det ( ϕ s ) d x w H 1 ( Ω t ) d = O u ˜ t + s H 1 ( Ω t ) d + f L 2 ( Ω t ) d .
Taking the supremum in (83) over admissible w, in the virtue of (82) and LBB condition (71), we get:
λ ˜ t + s L 2 ( Ω t ) = O ( 1 ) .
By the reflexivity of the underlying function spaces, from (82) and (84) it follows that there exists an accumulation point ( u ˜ t , λ ˜ t ) and subsequence ( u ˜ t + s k , λ ˜ t + s k ) , such that
( u ˜ t + s k , λ ˜ t + s k ) ( u ˜ t , λ ˜ t ) weakly in H Γ 1 ( Ω t ) d × L 2 ( Ω t ) as s k 0 + .
From asymptotic Relations (79) and (85), we get the limit of the system of linear equations (77) and (78) in the form of OS (72) and (73). Then, due to the uniqueness of its solution, the accumulation point implies
( u ˜ t , λ ˜ t ) = ( u t , λ t ) .
Using Representations (70), (72), (73), and algebra formula 1 2 ( a b ) 2 = ( a b ) b + 1 2 a 2 1 2 b 2 we rearrange the following terms:
Ω t μ 2 | ( u ˜ t + s k u t ) | 2 + α 2 | u ˜ t + s k u t | 2 d x = Ω t ( μ ( u ˜ t + s k u t ) · u t + α ( u ˜ t + s k u t ) · u t + f · ( u ˜ t + s k u t ) + λ ˜ t + s k div u ˜ t + s k ) d x + L ( u ˜ t + s k , λ ˜ t + s k ; Ω t ) L ( u t , λ ˜ t + s k ; Ω t ) .
In the virtue of asymptotic Formula (79), from (74) and (43), with test function w = u t , it follows
L ( u ˜ t + s k , λ ˜ t + s k ; Ω t ) L ( u t , λ ˜ t + s k ; Ω t ) L ˜ ( s k , u ˜ t + s k , λ ˜ t + s k ; Ω t ) L ˜ ( s k , u t , λ ˜ t + s k ; Ω t ) + c s k c s k
with some constant c > 0 , and div u ˜ t + s k = O ( s k ) according to (78), hence (87) yields that
lim   sup s k 0 + u ˜ t + s k u t H 1 ( Ω t ) d 2 0 .
Subtracting (72) from (77), due to (79) and (88) we have the asymptotic equality
Ω t ( λ ˜ t + s k λ t ) div w d x = Ω t μ ( u ˜ t + s k u t ) · w + α ( u ˜ t + s k u t ) · w d x + O ( s k ) .
Therefore, using (89) together with (88), the convergence in (85) is updated to a strong one:
( u ˜ t + s k , λ ˜ t + s k ) ( u ˜ t , λ ˜ t ) strongly in H Γ 1 ( Ω t ) d × L 2 ( Ω t ) as s k 0 + ,
implying (47) and (49) in Hypotheses (H3) and (H4).
Finally, due to the continuity of mapping ( τ , w ) s F ˜ ( τ , w ; Ω t ) in (80) and using the strong Convergence (90), we have:
lim τ , s k 0 + s F ˜ ( τ , u ˜ t + s k ; Ω t ) = s F ˜ ( 0 , u t ; Ω t ) ,
that proves (48). Applying Theorem 1, from Formula (80) at τ = 0 with ( w , p ) = ( u t , λ t ) and using div u t = 0 it follows Formula (75) of the shape derivative and finishes the proof. □
We remark a singularity of the mixed Dirichlet–Neumann boundary value Problem (72) and (73) at intersection Γ t D ¯ Γ t N ¯ , such that its solution ( u t , λ t ) is generally not in H 2 ( Ω t ) d × H 1 ( Ω t ) . Therefore:
Corollary 1.
Let the singular set be localized in a domain ω t ¯ Ω t ¯ such that the solution ( u t , λ t ) of (72) and (73) ( u t , λ t ) H 2 ( Ω t \ ω t ) d × H 1 ( Ω t \ ω t ) , force f const and velocity Λ const in ω t . In this case, using integration of (75) by parts, the following Hadamard representation over the boundary of Ω t \ ω t with the outward normal n t holds true:
s L ˜ ( 0 , u t , λ t ; Ω t ) = ( Ω t \ ω t ) Λ · n t ( μ 2 | u t | 2 + α 2 | u t | 2 f · u t ) u t ( μ n t u t λ t n t ) d S x .
Proof. 
In domain D with regular boundary D and unit normal vector n outward to D, for smooth functions w and p, we used the following formulas of integration by parts written component-wisely with the convention of summation over repeated indices i , j , k = 1 , , d :
D ( 1 2 Λ k , k w i , j 2 w i , k Λ k , j w i , j ) d x = D Λ k w i , j j w i , k d x + D ( Λ k n k 1 2 w i , j 2 Λ k w i , k w i , j n j ) d S x , D 1 2 Λ k , k w i 2 d x = D Λ k w i , k w i d x + D Λ k n k 1 2 w i 2 d S x , D ( Λ k , k f i + Λ k f i , k ) w i d x = D Λ k f i w i , k d x D Λ k n k f i w i d S x , D p Λ k , i w i , k d x = D Λ k ( p i w i , k + p w i , i k ) d x + D Λ k w i , k p n i d S x .
The summation of four equations in (93) (where the first equation is multiplied with μ and the second one with α ) provides the identity
D μ 2 div Λ | w | 2 μ w · ( Λ w ) + α 2 div Λ | w | 2 ( ( div Λ ) f + ( Λ · ) f ) · w + p Λ T · w d x = D ( Λ · w ) ( μ Δ w α w + f p ) p Λ · div w d x + D Λ · n ( μ 2 | u t | 2 + α 2 | w | 2 f · w ) w ( μ n w p n ) d S x .
The integral in (75) over ω t vanishes when f 0 and Λ 0 here. For the complement integral over Ω t \ ω t , we applied Formula (94) and Equations (60), (61) to derive (92). The proof is completed. □
Remark 4.
If the parameter α = 0 , then (60) turns into the Stokes flow equation. In this case, the shape differentiability result of Theorem 2 holds true for the Stokes problem and was proved by direct method in the earlier work [9].
Remark 5.
If parameter μ = 0 , then (60) turns into the equation describing Darcy flow [16]. In this case, function space V ( Ω t ) = { w H ( div ; Ω t ) | w = 0 a . e . Γ t D } , where
H ( div ; Ω t ) = { w = ( w 1 , , w d ) L 2 ( Ω t ) d | div w L 2 ( Ω t ) } .
It fails bijection property (34) in (D1), because div y v = ( ϕ s T ϕ s ) · x ( v ϕ s ) L 2 ( Ω t ) does not imply div x ( v ϕ s ) L 2 ( Ω t ) and vice versa (see also Remark 2).

5. Conclusions

The current study extends the class of shape differential models to the incompressible Brinkman flow problem under mixed Dirichlet–Neumann boundary conditions. This result can be continued to flow problems under pure no-slip boundary conditions for area-preserving maps within L 0 2 ( Ω t ) -spaces (see Remark 3) that form special linear group SL ( d ) in R d . Future research in this field might use divergence-preserving maps (see Remark 2 and Remark 5) like Piola transformation, which challenges a rigorous mathematical tool. The other challenging direction concerns examining the geometric singularity of the mixed Dirichlet–Neumann boundary value problem (72) and (73) at intersection point Γ t D ¯ Γ t N ¯ .

Author Contributions

Supervision, J.G.; Project Administration, V.A.K.; Funding Acquisition, J.R.G.G.

Funding

This research was funded by the Austrian Science Fund (FWF) grant number P26147-N26: “Object identification problems: numerical analysis” (PION).

Acknowledgments

V.A.K. thanks the Austrian Academy of Sciences (OeAW) for their partial support, the Colombian Institute for Education and Technical Studies Abroad (ICETEX), and Universidad Tecnológica de Pereira (UTP) for supporting his visit of UTP and the Conference ICAMI 2017 at San Andrés Island, Colombia; and the Chair of the Conference Organizing Committee O.O. Vasilieva for their help.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
COconstrained optimization
CTcoordinate transformation
LBBLadyzhenskaya–Babuška–Brezzi
ODEordinary differential equation
OSoptimality system
OVoptimal value
SDshape derivative
SPsaddle-point
VIvariational inequality

References

  1. Gwinner, J.; Ovcharova, N. From solvability and approximation of variational inequalities to solution of nondifferentiable optimization problems in contact mechanics. Optimization 2015, 64, 1683–1702. [Google Scholar] [CrossRef]
  2. Gwinner, J.; Stephan, E.P. Advanced Boundary Element Methods: Treatment of Boundary Value, Transmission and Contact Problems; Springer: Berlin, Germany, 2018; ISBN 978-3-319-92000-9. [Google Scholar]
  3. Kikuchi, N.; Oden, J.T. Contact Problems in Elasticity: A Study of Variational Inequalities and Finite Element Methods; SIAM: Philadelphia, PA, USA, 1988; ISBN 978-0-8987-1468-5. [Google Scholar]
  4. Khludnev, A.M.; Kovtunenko, V.A. Analysis of Cracks in Solids; WIT-Press: Southampton, UK, 2000; ISBN 978-1-8531-2625-3. [Google Scholar]
  5. Lazarev, N.P.; Rudoy, E.M. Shape sensitivity analysis of Timoshenko’s plate with a crack under the nonpenetration condition. Z. Angew. Math. Mech. 2014, 94, 730–739. [Google Scholar] [CrossRef]
  6. Shcherbakov, V.V. Shape derivative of the energy functional for the bending of elastic plates with thin defects. J. Phys. Conf. Ser. 2017, 894, 012084. [Google Scholar] [CrossRef] [Green Version]
  7. Khludnev, A.M.; Sokolowski, J. Modelling and Control in Solid Mechanics; Birkhäuser: Basel, Switzerland, 1997; ISBN 978-3-0348-8984-1. [Google Scholar]
  8. Chesnokov, A.A.; Kovtunenko, P.V. Weak discontinuities in solutions of long-wave equations for viscous flow. Stud. Appl. Math. 2014, 132, 50–64. [Google Scholar] [CrossRef]
  9. Kovtunenko, V.A.; Ohtsuka, K. Shape differentiability of Lagrangians and application to Stokes problem. SIAM J. Control Optim. 2018, 56, 3668–3684. [Google Scholar] [CrossRef]
  10. Moubachir, M.; Zolésio, J.-P. Moving Shape Analysis and Control; Chapman & Hall/CRC: Boca Raton, FL, USA, 2006; ISBN 978-1-5848-8611-2. [Google Scholar]
  11. De los Reyes, J.C. On the optimal control of some nonsmooth distributed parameter systems arising in mechanics. GAMM-Mitteilungen 2018, 40, 268–286. [Google Scholar] [CrossRef] [Green Version]
  12. Hintermüller, M.; Surowiec, T. First-order optimality conditions for elliptic mathematical programs with equilibrium constraints via variational analysis. SIAM J. Optim. 2011, 21, 1561–1593. [Google Scholar] [CrossRef]
  13. Gwinner, J. On two-coefficient identification in elliptic variational inequalities. Optimization 2018, 67, 1–4. [Google Scholar] [CrossRef]
  14. Gwinner, J.; Jadamba, B.; Khan, A.A.; Sama, M. Identification in variational and quasi-variational inequalities. J. Convex Anal. 2018, 25, 545–569. [Google Scholar]
  15. Radwan, A.; Vasilieva, O.; Enkhbat, R.; Griewank, A.; Guddat, J. Parametric approach to optimal control. Optim. Lett. 2012, 6, 1303–1316. [Google Scholar] [CrossRef]
  16. Batzel, J.J.; Bachar, M.; Kappel, F. Mathematical Modeling and Validation in Physiology: Applications to the Cardiovascular and Respiratory Systems; Springer: Berlin/Heidelberg, Germany, 2012; ISBN 978-3-6423-2882-4. [Google Scholar]
  17. Hasanov Hasanoğlu, A.; Romanov, V.G. Introduction to Inverse Problems for Differential Equations; Springer: Berlin, Germany, 2017; ISBN 978-3-319-62796-0. [Google Scholar]
  18. Lavrentiev, M.M.; Avdeev, A.V.; Lavrentiev, M.M., Jr.; Priimenko, V.I. Inverse Problems of Mathematical Physics; De Gruyter: Berlin, Germany, 2012; ISBN 978-9-0676-4396-2. [Google Scholar]
  19. Alekseev, G.V. Invisibility Problem in Acoustics Optics and Heat Transfer; Dalnauka: Vladivostok, Russia, 2016; ISBN 978-5-8044-1045-3. (In Russian) [Google Scholar]
  20. Cakoni, F.; Kovtunenko, V.A. Topological optimality condition for the identification of the center of an inhomogeneity. Inverse Probl. 2018, 34, 035009. [Google Scholar] [CrossRef] [Green Version]
  21. Kovtunenko, V.A.; Kunisch, K. High precision identification of an object: Optimality-conditions-based concept of imaging. SIAM J. Control Optim. 2014, 52, 773–796. [Google Scholar] [CrossRef]
  22. Haslinger, J.; Ito, K.; Kozubek, T.; Kunisch, K.; Peichl, G. On the shape derivative for problems of Bernoulli type. Interfaces Free Bound. 2009, 11, 317–330. [Google Scholar] [CrossRef]
  23. Ito, K.; Kunisch, K. Lagrange Multiplier Approach to Variational Problems and Applications; SIAM: Philadelphia, PA, USA, 2008; ISBN 978-0-8987-1649-8. [Google Scholar]
  24. Hill, A.A.; Rajagopal, K.R.; Vergori, L. On the stability and uniqueness of the flow of a fluid through a porous medium. Z. Angew. Math. Phys. 2016, 67, 49. [Google Scholar] [CrossRef]
  25. Nield, D.; Bejan, A. Convection in Porous Media; Springer: New York, NY, USA, 2006; ISBN 978-0387-29096-6. [Google Scholar]
  26. Kohr, M.; Lanza de Cristoforis, M.; Mikhailov, S.E.; Wendland, W.L. Integral potential method for transmission problem with Lipschitz interface in R 3 for the Stokes and Darcy–Forchheimer–Brinkman PDE systems. Z. Angew. Math. Phys. 2016, 67, 116. [Google Scholar] [CrossRef]
  27. Angot, P. Well-posed Stokes/Brinkman and Stokes/Darcy coupling revisited with new jump interface conditions. ESAIM Math. Model. Numer. Anal. 2018. [Google Scholar] [CrossRef]
  28. Bernardi, C.; Hecht, F.; Pironneau, O. Coupling Darcy and Stokes equations for porous media with cracks. ESAIM Math. Model. Numer. Anal. 2005, 39, 7–35. [Google Scholar] [CrossRef] [Green Version]
  29. Abello Muñoz, C.A.; Cárdenas Alzate, P.P.; González Granada, J.R. On the coercive functions and minimizers. Adv. Stud. Theor. Phys. 2017, 11, 709–715. [Google Scholar] [CrossRef]
  30. Ekeland, I.; Temam, R. Convex Analysis and Variational Problems; North-Holland: Amsterdam, The Netherlands, 1976; ISBN 978-0-4441-0898-2. [Google Scholar]
  31. Delfour, M.C.; Zolésio, J.-P. Shape and Geometries: Metrics, Analysis, Differential Calculus, and Optimization; SIAM: Philadelphia, PA, USA, 2011; ISBN 978-0-898719-36-9. [Google Scholar]
  32. Sokolowski, J.; Zolésio, J.-P. Introduction to Shape Optimization. Shape Sensitivity Analysis; Springer: Berlin/ Heidelberg, Germany, 1992; ISBN 978-3-6426-3471-0. [Google Scholar]
  33. Kovtunenko, V.A.; Leugering, G. A shape-topological control problem for nonlinear crack - defect interaction: the anti-plane variational model. SIAM J. Control Optim. 2016, 54, 1329–1351. [Google Scholar] [CrossRef]
  34. Kovtunenko, V.A. Primal-dual methods of shape sensitivity analysis for curvilinear cracks with non-penetration. IMA J. Appl. Math. 2006, 71, 635–657. [Google Scholar] [CrossRef]
  35. Kovtunenko, V.A.; Kunisch, K. Problem of crack perturbation based on level sets and velocities. Z. Angew. Math. Mech. 2007, 87, 809–830. [Google Scholar] [CrossRef]
  36. Hintermüller, M.; Kovtunenko, V.A. From shape variation to topology changes in constrained minimization: A velocity method-based concept. Optim. Methods Softw. 2011, 26, 513–532. [Google Scholar] [CrossRef]
  37. Walker, S.W. The Shapes of Things: A Practical Guide to Differential Geometry and the Shape Derivative; SIAM: Philadelphia, PA, USA, 2015; ISBN 978-1-6119-7395-2. [Google Scholar]

Share and Cite

MDPI and ACS Style

González Granada, J.R.; Gwinner, J.; Kovtunenko, V.A. On the Shape Differentiability of Objectives: A Lagrangian Approach and the Brinkman Problem. Axioms 2018, 7, 76. https://doi.org/10.3390/axioms7040076

AMA Style

González Granada JR, Gwinner J, Kovtunenko VA. On the Shape Differentiability of Objectives: A Lagrangian Approach and the Brinkman Problem. Axioms. 2018; 7(4):76. https://doi.org/10.3390/axioms7040076

Chicago/Turabian Style

González Granada, José Rodrigo, Joachim Gwinner, and Victor A. Kovtunenko. 2018. "On the Shape Differentiability of Objectives: A Lagrangian Approach and the Brinkman Problem" Axioms 7, no. 4: 76. https://doi.org/10.3390/axioms7040076

APA Style

González Granada, J. R., Gwinner, J., & Kovtunenko, V. A. (2018). On the Shape Differentiability of Objectives: A Lagrangian Approach and the Brinkman Problem. Axioms, 7(4), 76. https://doi.org/10.3390/axioms7040076

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop