You are currently viewing a beta version of our website. If you spot anything unusual, kindly let us know.

Preprint

Article

Duality Principles and Numerical Procedures for a Large Class of Non-convex Models in the Calculus of Variations

Altmetrics

Downloads

3935

Views

3507

Comments

Fabio Botelho^*

This version is not peer-reviewed

Abstract

This article develops duality principles and numerical results for a large class of non-convex variational models. The main results are based on fundamental tools of convex analysis, duality theory and calculus of variations. More specifically the approach is established for a class of non-convex functionals similar as those found in some models in phase transition. Finally, in the last section we present a concerning numerical example and the respective software.

Keywords:

Subject: Computer Science and Mathematics - Applied Mathematics

1. Introduction

In this section we establish a dual formulation for a large class of models in non-convex optimization.

The main duality principle is applied to double well models similar as those found in the phase transition theory.

Such results are based on the works of J.J. Telega and W.R. Bielski [2,3,14,15] and on a D.C. optimization approach developed in Toland [16].

About the other references, details on the Sobolev spaces involved are found in [1]. Related results on convex analysis and duality theory are addressed in [5,6,7,9,13].

Finally, in this text we adopt the standard Einstein convention of summing up repeated indices, unless otherwise indicated.

In order to clarify the notation, here we introduce the definition of topological dual space.

Definition 1.1

(Topological dual spaces). Let U be a Banach space. We shall define its dual topological space, as the set of all linear continuous functionals defined on U. We suppose such a dual space of U, may be represented by another Banach space

U^{*}

, through a bilinear form

{〈 \cdot, \cdot 〉}_{U} : U \times U^{*} \to R

(here we are referring to standard representations of dual spaces of Sobolev and Lebesgue spaces). Thus, given

f : U \to R

linear and continuous, we assume the existence of a unique

u^{*} \in U^{*}

such that

\begin{matrix} f (u) = {〈 u, u^{*} 〉}_{U}, \forall u \in U . \end{matrix}

(1)

The norm of f , denoted by

{∥ f ∥}_{U^{*}}

, is defined as

\begin{matrix} {∥ f ∥}_{U^{*}} = sup_{u \in U} {| {〈 u, u^{*} 〉}_{U} {| : ∥ u ∥}_{U} \leq 1} \equiv ∥ u^{*} ∥_{U^{*}} . \end{matrix}

(2)

At this point we start to describe the primal and dual variational formulations.

2. A General Duality Principle Non-Convex Optimization

In this section we present a duality principle applicable to a model in phase transition.

This case corresponds to the vectorial one in the calculus of variations.

Let

Ω \subset R^{n}

be an open, bounded, connected set with a regular (Lipschitzian) boundary denoted by

\partial Ω .

Consider a functional

J : V \to R

where

J (u) = F (\nabla u_{1}, \dots, \nabla u_{N}) + G (u_{1}, \dots, u_{N}) - {〈 u_{i}, f_{i} 〉}_{L^{2}},

and where

V = {u = (u_{1}, \dots, u_{N}) \in W^{1, p} (Ω; R^{N}) : u = u_{0} on \partial Ω},

f \in L^{2} (Ω; R^{N}),

and

1 < p < + \infty .

We assume there exists

α \in R

such that

α = inf_{u \in V} J (u) .

Moreover, suppose F and G are Fréchet differentiable but not necessarily convex. A global optimum point may not be attained for J so that the problem of finding a global minimum for J may not be a solution.

Anyway, one question remains, how the minimizing sequences behave close the infimum of J.

We intend to use duality theory to approximately solve such a global optimization problem.

Denoting

V_{0} = W_{0}^{1, p} (Ω; R^{N})

Y_{1} = Y_{1}^{*} = L^{2} (Ω; R^{N \times n})

Y_{2} = Y_{2}^{*} = L^{2} (Ω; R^{N \times n})

Y_{3} = Y_{3}^{*} = L^{2} (Ω; R^{N}),

at this point we define,

F_{1} : V \times V_{0} \to R

G_{1} : V \to R

G_{2} : V \to R

G_{3} : V_{0} \to R

and

G_{4} : V \to R,

\begin{matrix} F_{1} (\nabla u, \nabla ϕ) & = & F (\nabla u_{1} + \nabla ϕ_{1}, \dots, \nabla u_{N} + \nabla ϕ_{N}) + \frac{K}{2} \int_{Ω} \nabla u_{j} \cdot \nabla u_{j} d x \\ + \frac{K_{2}}{2} \int_{Ω} \nabla ϕ_{j} \cdot \nabla ϕ_{j} d x \end{matrix}

(3)

and

G_{1} (u_{1}, \dots, u_{n}) = G (u_{1}, \dots, u_{N}) + \frac{K_{1}}{2} \int_{Ω} u_{j} u_{j} d x - {〈 u_{i}, f_{i} 〉}_{L^{2}},

G_{2} (\nabla u_{1}, \dots, \nabla u_{N}) = \frac{K_{1}}{2} \int_{Ω} \nabla u_{j} \cdot \nabla u_{j} d x,

G_{3} (\nabla ϕ_{1}, \dots, \nabla ϕ_{N}) = \frac{K_{2}}{2} \int_{Ω} \nabla ϕ_{j} \cdot \nabla ϕ_{j} d x,

and

G_{4} (u_{1}, \dots, u_{N}) = \frac{K_{1}}{2} \int_{Ω} u_{j} u_{j} d x .

Define now

J_{1} : V \times V_{0} \to R

J_{1} (u, ϕ) = F (\nabla u + \nabla ϕ) + G (u) - {〈 u_{i}, f_{i} 〉}_{L^{2}} .

Observe that

\begin{matrix} J_{1} (u, ϕ) & = & F_{1} (\nabla u, \nabla ϕ) + G_{1} (u) - G_{2} (\nabla u) - G_{3} (\nabla ϕ) - G_{4} (u) \\ \leq & F_{1} (\nabla u, \nabla ϕ) + G_{1} (u) - {〈 \nabla u, z_{1}^{*} 〉}_{L^{2}} - {〈 \nabla ϕ, z_{2}^{*} 〉}_{L^{2}} - {〈 u, z_{3}^{*} 〉}_{L^{2}} \\ + sup_{v_{1} \in Y_{1}} {{〈 v_{1}, z_{1}^{*} 〉}_{L^{2}} - G_{2} (v_{1})} \\ + sup_{v_{2} \in Y_{2}} {{〈 v_{2}, z_{2}^{*} 〉}_{L^{2}} - G_{3} (v_{2})} \\ + sup_{u \in V} {{〈 u, z_{3}^{*} 〉}_{L^{2}} - G_{4} (u)} \\ = & F_{1} (\nabla u, \nabla ϕ) + G_{1} (u) - {〈 \nabla u, z_{1}^{*} 〉}_{L^{2}} - {〈 \nabla ϕ, z_{2}^{*} 〉}_{L^{2}} - {〈 u, z_{3}^{*} 〉}_{L^{2}} \\ + G_{2}^{*} (z_{1}^{*}) + G_{3}^{*} (z_{2}^{*}) + G_{4}^{*} (z_{3}^{*}) \\ = & J_{1}^{*} (u, ϕ, z^{*}), \end{matrix}

(4)

\forall u \in V, ϕ \in V_{0}, z^{*} = (z_{1}^{*}, z_{2}^{*}, z_{3}^{*}) \in Y^{*} = Y_{1}^{*} \times Y_{2}^{*} \times Y_{3}^{*}

Here we assume

K, K_{1}, K_{2}

are large enough so that

F_{1}

and

G_{1}

are convex.

Hence, from the general results in [16], we may infer that

\begin{matrix} inf_{(u, ϕ) \in V \times V_{0}} J (u, ϕ) & = & inf_{(u, ϕ, z^{*}) \in V \times V_{0} \times Y^{*}} J_{1}^{*} (u, ϕ, z^{*}) . \end{matrix}

(5)

On the other hand

inf_{u \in V} J (u) \geq inf_{(u, ϕ) \in V \times V_{0}} J_{1} (u, ϕ) \geq inf_{u \in V} Q_{J} (u) = inf_{u \in V} J (u),

where

Q_{J} (u)

refers to a standard quasi-convex regularization of J.

From these last two results we may obtain

inf_{u \in V} J (u) = inf_{(u, ϕ, z^{*}) \in V \times V_{0} \times Y^{*}} J_{1}^{*} (u, ϕ, z^{*}) .

Moreover, from standards results on convex analysis, we may have

\begin{matrix} inf_{u \in V} J_{1}^{*} (u, ϕ, z^{*}) & = & inf_{u \in V} {F_{1} (\nabla u, \nabla ϕ) + G_{1} (u) \\ - {〈 \nabla u, z_{1}^{*} 〉}_{L^{2}} - {〈 \nabla ϕ, z_{2}^{*} 〉}_{L^{2}} - {〈 u, z_{3}^{*} 〉}_{L^{2}} \\ + G_{2}^{*} (z_{1}^{*}) + G_{3}^{*} (z_{2}^{*}) + G_{4}^{*} (z_{3}^{*})} \\ = & sup_{(v_{1}^{*}, v_{2}^{*}) \in C^{*}} {- F_{1}^{*} (v_{1}^{*} + z_{1}^{*}, \nabla ϕ) - G_{1}^{*} (v_{2}^{*} + z_{3}^{*}) - {〈 \nabla ϕ, z_{2}^{*} 〉}_{L^{2}} \\ + G_{2}^{*} (z_{1}^{*}) + G_{3}^{*} (z_{2}^{*}) + G_{4}^{*} (z_{3}^{*})}, \end{matrix}

(6)

where

C^{*} = {v^{*} = (v_{1}^{*}, v_{2}^{*}) \in Y_{1}^{*} \times Y_{3}^{*} : - div {(v_{1}^{*})}_{i} + {(v_{2}^{*})}_{i} = 0, \forall i \in {1, \dots, N}},

F_{1}^{*} (v_{1}^{*} + z_{1}^{*}, \nabla ϕ) = sup_{v_{1} \in Y_{1}} {{〈 v_{1}, z_{1}^{*} + v_{1}^{*} 〉}_{L^{2}} - F_{1} (v_{1}, \nabla ϕ)},

and

G_{1}^{*} (v_{2}^{*} + z_{2}^{*}) = sup_{u \in V} {{〈 u, v_{2}^{*} + z_{2}^{*} 〉}_{L^{2}} - G_{1} (u)} .

Thus, defining

J_{2}^{*} (ϕ, z^{*}, v^{*}) = F_{1}^{*} (v_{1}^{*} + z_{1}^{*}, \nabla ϕ) - G_{1}^{*} (v_{2}^{*} + z_{3}^{*}) - {〈 \nabla ϕ, z_{2}^{*} 〉}_{L^{2}} + G_{2}^{*} (z_{1}^{*}) + G_{3}^{*} (z_{2}^{*}) + G_{4}^{*} (z_{3}^{*}),

we have got

\begin{matrix} inf_{u \in V} J (u) & = & inf_{(u, ϕ) \in V \times V_{0}} J_{1} (u, ϕ) \\ = & inf_{(u, ϕ, z^{*}) \in V \times V_{0} \times Y^{*}} J_{1}^{*} (u, ϕ, z^{*}) \\ = & inf_{z^{*} \in Y^{*}} \{inf_{ϕ \in V_{0}} \{sup_{v^{*} \in C^{*}} J_{2}^{*} (ϕ, z^{*}, v^{*})\}\} . \end{matrix}

(7)

Finally, observe that

\begin{matrix} inf_{u \in V} J (u) \\ = & inf_{z^{*} \in Y^{*}} \{inf_{ϕ \in V_{0}} \{sup_{v^{*} \in C^{*}} J_{2}^{*} (ϕ, z^{*}, v^{*})\}\} \\ \geq & sup_{v^{*} \in C^{*}} \{inf_{(z^{*}, ϕ) \in Y^{*} \times V_{0}} J_{2}^{*} (ϕ, z^{*}, v^{*})\} . \end{matrix}

(8)

This last variational formulation corresponds to a concave relaxed formulation in

v^{*}

concerning the original primal formulation.

3. Another Duality Principle for a Simpler Related Model in Phase Transition with a Respective Numerical Example

In this section we present another duality principle for a related model in phase transition.

Let

Ω = [0, 1] \subset R

and consider a functional

J : V \to R

where

J (u) = \frac{1}{2} \int_{Ω} {({(u^{'})}^{2} - 1)}^{2} d x + \frac{1}{2} \int_{Ω} u^{2} d x - {〈 u, f 〉}_{L^{2}},

and where

V = {u \in W^{1, 4} (Ω) : u (0) = 0 and u (1) = 1 / 2}

and

f \in L^{2} (Ω) .

A global optimum point is not attained for J so that the problem of finding a global minimum for J has no solution.

Anyway, one question remains, how the minimizing sequences behave close the infimum of J.

We intend to use duality theory to approximately solve such a global optimization problem.

Denoting

V_{0} = W_{0}^{1, 4} (Ω)

, at this point we define,

F : V \to R

and

F_{1} : V \times V_{0} \to R

F (u) = \frac{1}{2} \int_{Ω} {({(u^{'})}^{2} - 1)}^{2} d x,

and

F_{1} (u, ϕ) = \frac{1}{2} \int_{Ω} {({(u^{'} + ϕ^{'})}^{2} - 1)}^{2} d x .

Observe

F (u) \geq inf_{ϕ \in V_{0}} F_{1} (u, ϕ) \geq Q_{F} (u), \forall u \in V,

where

Q_{F} (u)

refers to a quasi-convex regularization of

F .

We define also

F_{2} : V \times V_{0} \to R,

F_{3} : V \times V_{0} \to R

and

G : V \times V_{0} \to R

F_{2} (u, ϕ) = \frac{1}{2} \int_{Ω} {({(u^{'} + ϕ^{'})}^{2} - 1)}^{2} d x + \frac{1}{2} \int_{Ω} u^{2} d x - {〈 u, f 〉}_{L^{2}},

\begin{matrix} F_{3} (u, ϕ) & = & F_{2} (u, ϕ) + \frac{K}{2} \int_{Ω} {(u^{'})}^{2} d x \\ + \frac{K_{1}}{2} \int_{Ω} {(ϕ^{'})}^{2} d x \end{matrix}

(9)

and

\begin{matrix} G (u, ϕ) & = & \frac{K}{2} \int_{Ω} {(u^{'})}^{2} d x \\ + \frac{K_{1}}{2} \int_{Ω} {(ϕ^{'})}^{2} d x \end{matrix}

(10)

Observe that if

K > 0, K_{1} > 0

is large enough, both

F_{3}

and G are convex.

Denoting

Y = Y^{*} = L^{2} (Ω)

we also define the polar functional

G^{*} : Y^{*} \times Y^{*} \to R

G^{*} (v^{*}, v_{0}^{*}) = sup_{(u, ϕ) \in V \times V_{0}} {{〈 u, v^{*} 〉}_{L^{2}} + {〈 ϕ, v_{0}^{*} 〉}_{L^{2}} - G (u, ϕ)} .

Observe that

inf_{u \in U} J (u) \geq inf_{((u, ϕ), (v^{*}, v_{0}^{*})) \in V \times V_{0} \times {[Y^{*}]}^{2}} {G^{*} (v^{*}, v_{0}^{*}) - {〈 u, v^{*} 〉}_{L^{2}} - {〈 ϕ, v_{0}^{*} 〉}_{L^{2}} + F_{3} (u, ϕ)} .

With such results in mind, we define a relaxed primal dual variational formulation for the primal problem, represented by

J_{1}^{*} : V \times V_{0} \times {[Y^{*}]}^{2} \to R

, where

J_{1}^{*} (u, ϕ, v^{*}, v_{0}^{*}) = G^{*} (v^{*}, v_{0}^{*}) - {〈 u, v^{*} 〉}_{L^{2}} - {〈 ϕ, v_{0}^{*} 〉}_{L^{2}} + F_{3} (u, ϕ) .

Having defined such a functional, we may obtain numerical results by solving a sequence of convex auxiliary sub-problems, through the following algorithm.

Set $K \approx 150$ and $K_{1} = K / 20$ and $0 < ε ≪ 1 .$
Choose $(u_{1}, ϕ_{1}) \in V \times V_{0},$ such that $∥ u_{1} ∥_{1, \infty} ≪ K / 4$ and $∥ ϕ_{1} ∥_{1, \infty} ≪ K / 4 .$
Set $n = 1 .$
Calculate $(v_{n}^{*}, {(v_{0}^{*})}_{n})$ solution of the system of equations:

$\frac{\partial J_{1}^{*} (u_{n}, ϕ_{n}, v_{n}^{*}, {(v_{0}^{*})}_{n})}{\partial v^{*}} = 0$

and

$\frac{\partial J_{1}^{*} (u_{n}, ϕ_{n}, v_{n}^{*}, {(v_{0}^{*})}_{n})}{\partial v_{0}^{*}} = 0,$

that is

$\frac{\partial G^{*} (v_{n}^{*}, {(v_{0}^{*})}_{n})}{\partial v^{*}} - u_{n} = 0$

and

$\frac{\partial G^{*} (v_{n}^{*}, {(v_{0}^{*})}_{n})}{\partial v_{0}^{*}} - ϕ_{n} = 0$

so that

$v_{n}^{*} = \frac{\partial G (u_{n}, ϕ_{n})}{\partial u}$

and

${(v_{0}^{*})}_{n}^{*} = \frac{\partial G (u_{n}, ϕ_{n})}{\partial ϕ}$
Calculate $(u_{n + 1}, ϕ_{n + 1})$ by solving the system of equations:

$\frac{\partial J_{1}^{*} (u_{n + 1}, ϕ_{n + 1}, v_{n}^{*}, {(v_{0}^{*})}_{n})}{\partial u} = 0$

and

$\frac{\partial J_{1}^{*} (u_{n + 1}, ϕ_{n + 1}, v_{n}^{*}, {(v_{0}^{*})}_{n})}{\partial ϕ} = 0$

that is

$- v_{n}^{*} + \frac{\partial F_{3} (u_{n + 1}, ϕ_{n + 1})}{\partial u} = 0$

and

$- {(v_{0}^{*})}_{n} + \frac{\partial F_{3} (u_{n + 1}, ϕ_{n + 1})}{\partial ϕ} = 0$
If $max {∥ u_{n} - u_{n + 1} ∥_{\infty}, ∥ ϕ_{n + 1} - ϕ_{n} ∥_{\infty}} \leq ε$ , then stop, else set $n : = n + 1$ and go to item 4.

For the case in which

f (x) = 0

, we have obtained numerical results for

K = 1500

and

K_{1} = K / 20

. For such a concerning solution

u_{0}

obtained, please see Figure 1. For the case in which

f (x) = sin (π x) / 2

, we have obtained numerical results for

K = 100

and

K_{1} = K / 20

. For such a concerning solution

u_{0}

obtained, please see Figure 2.

Remark 3.1.

Observe that the solutions obtained are approximate critical points. They are not, in a classical sense, the global solutions for the related optimization problems. Indeed, such solutions reflect the average behavior of weak cluster points for concerning minimizing sequences.

4. A Convex Dual Variational Formulation for a Third Similar Model

In this section we present another duality principle for a third related model in phase transition.

Figure 2. solution

u_{0} (x)

for the case

f (x) = sin (π x) / 2

Figure 2. solution

u_{0} (x)

for the case

f (x) = sin (π x) / 2

Let

Ω = [0, 1] \subset R

and consider a functional

J : V \to R

where

J (u) = \frac{1}{2} \int_{Ω} min {{(u^{'} - 1)}^{2}, {(u^{'} + 1)}^{2}} d x + \frac{1}{2} \int_{Ω} u^{2} d x - {〈 u, f 〉}_{L^{2}},

and where

V = {u \in W^{1, 2} (Ω) : u (0) = 0 and u (1) = 1 / 2}

and

f \in L^{2} (Ω) .

A global optimum point is not attained for J so that the problem of finding a global minimum for J has no solution.

Anyway, one question remains, how the minimizing sequences behave close to the infimum of J.

We intend to use the duality theory to solve such a global optimization problem in an appropriate sense to be specified.

At this point we define,

F : V \to R

and

G : V \to R

\begin{matrix} F (u) & = & \frac{1}{2} \int_{Ω} min {{(u^{'} - 1)}^{2}, {(u^{'} + 1)}^{2}} d x \\ = & \frac{1}{2} \int_{Ω} {(u^{'})}^{2} d x - \int_{Ω} | u^{'} | d x + 1 / 2 \\ \equiv & F_{1} (u^{'}), \end{matrix}

(11)

and

G (u) = \frac{1}{2} \int_{Ω} u^{2} d x - {〈 u, f 〉}_{L^{2}} .

Denoting

Y = Y^{*} = L^{2} (Ω)

we also define the polar functional

F_{1}^{*} : Y^{*} \to R

and

G^{*} : Y^{*} \to R

\begin{matrix} F_{1}^{*} (v^{*}) & = & sup_{v \in Y} {{〈 v, v^{*} 〉}_{L^{2}} - F_{1} (v)} \\ = & \frac{1}{2} \int_{Ω} {(v^{*})}^{2} d x + \int_{Ω} | v^{*} | d x, \end{matrix}

(12)

and

\begin{matrix} G^{*} ({(v^{*})}^{'}) & = & sup_{u \in V} {- {〈 u^{'}, v^{*} 〉}_{L^{2}} - G (u)} \\ = & \frac{1}{2} \int_{Ω} {({(v^{*})}^{'} + f)}^{2} d x - \frac{1}{2} v^{*} (1) . \end{matrix}

(13)

Observe this is the scalar case of the calculus of variations, so that from the standard results on convex analysis, we have

inf_{u \in V} J (u) = max_{v^{*} \in Y^{*}} {- F_{1}^{*} (v^{*}) - G^{*} (- {(v^{*})}^{'})} .

Indeed, from the direct method of the calculus of variations, the maximum for the dual formulation is attained at some

{\hat{v}}^{*} \in Y^{*}

Moreover, the corresponding solution

u_{0} \in V

is obtained from the equation

u_{0} = \frac{\partial G ({({\hat{v}}^{*})}^{'})}{\partial {(v^{*})}^{'}} = {({\hat{v}}^{*})}^{'} + f .

Finally, the Euler-Lagrange equations for the dual problem stands for

\{\begin{matrix} {(v^{*})}^{″} + f^{'} - v^{*} - sign (v^{*}) = 0, & in Ω, \\ {(v^{*})}^{'} (0) = 0, {(v^{*})}^{'} (1) = 1 / 2, \end{matrix}

(14)

where

sign (v^{*} (x)) = 1

v^{*} (x) > 0,

sign (v^{*} (x)) = - 1,

v^{*} (x) < 0

and

- 1 \leq sign (v^{*} (x)) \leq 1,

v^{*} (x) = 0 .

We have computed the solutions

v^{*}

and corresponding solutions

u_{0} \in V

for the cases in which

f (x) = 0

and

f (x) = sin (π x) / 2 .

For the solution

u_{0} (x)

for the case in which

f (x) = 0

, please see Figure 3.

For the solution

u_{0} (x)

for the case in which

f (x) = sin (π x) / 2

, please see Figure 4.

Remark 4.1.

Observe that such solutions

u_{0}

obtained are not the global solutions for the related primal optimization problems. Indeed, such solutions reflect the average behavior of weak cluster points for concerning minimizing sequences.

4.1. The Algorithm through Which We Have Obtained the Numerical Results

In this subsection we present the software in MATLAB through which we have obtained the last numerical results.

This algorithm is for solving the concerning Euler-Lagrange equations for the dual problem, that is, for solving the equation

\{\begin{matrix} {(v^{*})}^{″} + f^{'} - v^{*} - sign (v^{*}) = 0, & in Ω, \\ {(v^{*})}^{'} (0) = 0, {(v^{*})}^{'} (1) = 1 / 2 . \end{matrix}

(15)

Here the concerning software in MATLAB. We emphasize to have used the smooth approximation

| v^{*} | \approx \sqrt{{(v^{*})}^{2} + e_{1}},

where a small value for

e_{1}

is specified in the next lines.

*************************************

clear all
$m_{8} = 800;$ (number of nodes)
$d = 1 / m_{8};$
$e_{1} = 0.00001;$
$f o r i = 1 : m_{8}$

$y o (i, 1) = 0.01;$

$y_{1} (i, 1) = sin (π * i / m_{8}) / 2;$

$e n d;$
$f o r i = 1 : m_{8} - 1$

$d y_{1} (i, 1) = (y_{1} (i + 1, 1) - y_{1} (i, 1)) / d;$

$e n d;$
$f o r k = 1 : 3000$ (we have fixed the number of iterations)

$i = 1;$

$h_{3} = 1 / \sqrt{v o {(i, 1)}^{2} + e_{1}};$

$m_{12} = 1 + d^{2} * h_{3} + d^{2};$

$m_{50} (i) = 1 / m_{12};$

$z (i) = m_{50} (i) * (d y_{1} (i, 1) * d^{2});$
$f o r i = 2 : m_{8} - 1$

$h_{3} = 1 / \sqrt{v o {(i, 1)}^{2} + e_{1}};$

$m_{12} = 2 + h_{3} * d^{2} + d^{2} - m 50 (i - 1);$

$m 50 (i) = 1 / m_{12};$

$z (i) = m_{50} (i) * (z (i - 1) + d y_{1} (i, 1) * d^{2});$

$e n d;$
$v (m_{8}, 1) = (d / 2 + z (m_{8} - 1)) / (1 - m_{50} (m_{8} - 1));$
$f o r i = 1 : m_{8} - 1$

$v (m_{8} - i, 1) = m_{50} (m_{8} - i) * v (m_{8} - i + 1) + z (m_{8} - i);$

$e n d;$
$v (m_{8} / 2, 1)$
$v o = v;$

$e n d;$
$f o r i = 1 : m_{8} - 1$

$u (i, 1) = (v (i + 1, 1) - v (i, 1)) / d + y_{1} (i, 1);$

$e n d;$
$f o r i = 1 : m_{8} - 1$

$x (i) = i * d;$

$e n d;$

$p l o t (x, u (:, 1))$

********************************

5. An Improvement of the Convexity Conditions for a Non-Convex Related Model through an Approximate Primal Formulation

In this section we develop an approximate primal dual formulation suitable for a large class of variational models.

Here, the applications are for the Kirchhoff-Love plate model, which may be found in Ciarlet, [10].

At this point we start to describe the primal variational formulation.

Let

Ω \subset R^{2}

be an open, bounded, connected set which represents the middle surface of a plate of thickness h. The boundary of

Ω

, which is assumed to be regular (Lipschitzian), is denoted by

\partial Ω

. The vectorial basis related to the cartesian system

{x_{1}, x_{2}, x_{3}}

is denoted by

(a_{α}, a_{3})

, where

α = 1, 2

(in general Greek indices stand for 1 or 2), and where

a_{3}

is the vector normal to

Ω

, whereas

a_{1}

and

a_{2}

are orthogonal vectors parallel to

Ω .

Also,

n

is the outward normal to the plate surface.

The displacements will be denoted by

\hat{u} = {{\hat{u}}_{α}, {\hat{u}}_{3}} = {\hat{u}}_{α} a_{α} + {\hat{u}}_{3} a_{3} .

The Kirchhoff-Love relations are

\begin{matrix} {\hat{u}}_{α} (x_{1}, x_{2}, x_{3}) = u_{α} (x_{1}, x_{2}) - x_{3} w {(x_{1}, x_{2})}_{, α} \\ and {\hat{u}}_{3} (x_{1}, x_{2}, x_{3}) = w (x_{1}, x_{2}) . \end{matrix}

(16)

Here

- h / 2 \leq x_{3} \leq h / 2

so that we have

u = (u_{α}, w) \in U

where

\begin{matrix} U & = & \{u = (u_{α}, w) \in W^{1, 2} (Ω; R^{2}) \times W^{2, 2} (Ω), \\ u_{α} = w = \frac{\partial w}{\partial n} = 0 on \partial Ω\} \\ = & W_{0}^{1, 2} (Ω; R^{2}) \times W_{0}^{2, 2} (Ω) . \end{matrix}

It is worth emphasizing that the boundary conditions here specified refer to a clamped plate.

We also define the operator

Λ : U \to Y \times Y

, where

Y = Y^{*} = L^{2} (Ω; R^{2 \times 2})

, by

Λ (u) = {γ (u), κ (u)},

γ_{α β} (u) = \frac{u_{α, β} + u_{β, α}}{2} + \frac{w_{, α} w_{, β}}{2},

κ_{α β} (u) = - w_{, α β} .

The constitutive relations are given by

N_{α β} (u) = H_{α β λ μ} γ_{λ μ} (u),

(17)

M_{α β} (u) = h_{α β λ μ} κ_{λ μ} (u),

(18)

where:

{H_{α β λ μ}}

and

\{h_{α β λ μ} = \frac{h^{2}}{12} H_{α β λ μ}\}

, are symmetric positive definite fourth order tensors. From now on, we denote

{{\bar{H}}_{α β λ μ}} = {H_{α β λ μ}}^{- 1}

and

{{\bar{h}}_{α β λ μ}} = {h_{α β λ μ}}^{- 1}

Furthermore

{N_{α β}}

denote the membrane force tensor and

{M_{α β}}

the moment one. The plate stored energy, represented by

(G \circ Λ) : U \to R

is expressed by

(G \circ Λ) (u) = \frac{1}{2} \int_{Ω} N_{α β} (u) γ_{α β} (u) d x + \frac{1}{2} \int_{Ω} M_{α β} (u) κ_{α β} (u) d x

(19)

and the external work, represented by

F : U \to R

, is given by

F (u) = {〈 w, P 〉}_{L^{2}} + {〈 u_{α}, P_{α} 〉}_{L^{2}},

(20)

where

P, P_{1}, P_{2} \in L^{2} (Ω)

are external loads in the directions

a_{3}

a_{1}

and

a_{2}

respectively. The potential energy, denoted by

J : U \to R

is expressed by:

J (u) = (G \circ Λ) (u) - F (u)

Define now

J_{3} : \tilde{U} \to R

J_{3} (u) = J (u) + 10 \int_{Ω} \frac{e^{K b w}}{K^{3 / 2}} d x .

In such a case for

b = P / | P |

Ω

K = 100

and

\tilde{U} = {{u \in U : ∥ w ∥}_{\infty} \leq 0.01 and P w \geq 0 a . e . in Ω},

we get

\begin{matrix} \frac{\partial J_{3} (w)}{\partial w} & = & \frac{\partial J (w)}{\partial w} + 10 b \frac{e^{K b w}}{K^{1 / 2}} \\ \approx & \frac{\partial J (w)}{\partial w} + O (1.5), \end{matrix}

(21)

and

\begin{matrix} \frac{\partial^{2} J_{3} (w)}{\partial w^{2}} & = & \frac{\partial^{2} J (w)}{\partial w^{2}} + 10 b^{2} e^{K b w} K^{1 / 2} \\ \approx & \frac{\partial^{2} J (w)}{\partial w^{2}} + O (50) . \end{matrix}

(22)

This new functional

J_{3}

has a relevant improvement in the convexity conditions concerning the previous functional J.

Indeed, we have obtained a gain in positiveness for the second variation

\frac{\partial^{2} J (u)}{\partial w^{2}},

which has increased of order

O (50) .

Moreover the difference between the approximate and exact equation

\frac{\partial J (u)}{\partial w} = 0

is of order

O (1.5)

which corresponds to a small perturbation in the original equation for a load of

P = 1500 N / m^{2},

for example. Summarizing, the exact equation may be approximate solved in an appropriate sense.

6. An Approximate Convex Variational Formulation for Another Related Model

In this section, we obtain an approximate convex variational formulation for a related model, more specifically, for a Ginzburg-Landau type equation.

Let

Ω \subset R^{3}

be an open, bounded and connected set with a regular (Lipschitzian) boundary denoted by

\partial Ω .

Consider a functional

J : V \to R

where

\begin{matrix} J (u) & = & \frac{γ}{2} \int_{Ω} \nabla u \cdot \nabla u d x + \frac{α}{2} \int_{Ω} {(u^{2} - β)}^{2} d x \\ - {〈 u, f 〉}_{L^{2}}, \end{matrix}

(23)

where

γ > 0

α > 0

β > 0

V = W_{0}^{1, 2} (Ω)

and

f \in L^{2} (Ω) .

We define

A^{+} = {u \in V : u f \geq 0, a . e . in Ω},

V_{2} = {{u \in V : ∥ u ∥}_{\infty} \leq 1},

and

V_{1} = V_{2} \cap A^{+} .

At this point we define

v = u / 10

so that

\begin{matrix} J (u) & = & J_{1} (v) \\ = & \frac{10^{2} γ}{2} \int_{Ω} \nabla v \cdot \nabla v d x + \frac{α}{2} \int_{Ω} ({(10 v)}^{2} - β) d x \\ - {〈 10 v, f 〉}_{L^{2}} . \end{matrix}

(24)

Moreover we define

\begin{matrix} J_{2} (v) & = & \frac{1}{10} J_{1} (v) \\ = & \frac{10 γ}{2} \int_{Ω} \nabla v \cdot \nabla v d x + \frac{α}{20} \int_{Ω} ({(10 v)}^{2} - β) d x \\ - {〈 v, f 〉}_{L^{2}}, \end{matrix}

(25)

and

J_{3} : U_{3} \to R

where

J_{3} (v) = J_{2} (v) + K_{1} \int_{Ω} \frac{a^{(K b v / 300)}}{ln (a) K^{3 / 2}} 300 d x,

where

K_{1} = 50 / 12

a = 2.71

K = 4755

b = f / | f |

Ω

and

U_{2} = {v \in V : f v \geq 0, a . e . in Ω},

U_{4} = {{v \in V : ∥ v ∥}_{\infty} \leq 1 / 10},

and

U_{3} = U_{2} \cap U_{4} .

Thus, with such numerical values, we may obtain

\begin{matrix} \frac{\partial J_{3} (v)}{\partial v} & = & \frac{\partial J_{2} (v)}{\partial v} + 0.0145019 K_{1} b 2 . 71^{(317 b v / 2) 0} \\ \approx & \frac{\partial J_{2} (v)}{\partial v} + O (0.25), \end{matrix}

(26)

and

\begin{matrix} \frac{\partial^{2} J_{3} (v)}{\partial v^{2}} & = & \frac{\partial^{2} J_{2} (v)}{\partial v^{2}} + 0.229154 K_{1} b^{2} 2 . 71^{(317 b v / 20)} \\ \approx & \frac{\partial^{2} J_{2} (v)}{\partial v^{2}} + O (2.5) . \end{matrix}

(27)

Remark 6.1.

This new functional

J_{1}

has a relevant improvement in the convexity conditions concerning the previous functional J.

Indeed, we have obtained a gain in positiveness for the second variation

\frac{\partial^{2} J_{2} (v)}{\partial v^{2}},

which has increased of order

O (2.5) .

Moreover the difference between the approximate and exact equation

\frac{\partial J_{2} (v)}{\partial v} = 0

is of order

O (0.25)

which for appropriate parameters

γ > 0, α > 0

and

β > 0

, corresponds to a small perturbation in the original equation. Summarizing, the exact equation may be approximately solved in an appropriate sense.

Finally, for this last example, we highlight it is relatively easy to improve even more both such an approximation quality and the convexity conditions concerning the original variational model.

With such statements and results in mind, we may prove the following theorem.

Theorem 6.2.

Suppose

γ > 0,

α > 0

and

β > 0

are such that

\frac{\partial^{2} J_{3} (v)}{\partial v^{2}} > 0,

U_{3}

Assume also,

v_{0} \in U_{3}

is such that

δ J_{3} (v_{0}) = 0 .

Under such hypotheses,

J_{3}

is convex on

U_{3}

so that

J_{3} (v_{0}) = min_{v \in U_{3}} J_{3} (v) .

Moreover,

δ J (u_{0}) = 0 + O (0.25),

where

u_{0} = 10 v_{0} \in V_{1}

Proof.

From the hypotheses

\frac{\partial^{2} J_{3} (v)}{\partial v^{2}} > 0

U_{3}

, so that

J_{3}

is convex on the convex set

U_{3} .

Consequently, since

δ J_{3} (v_{0}) = 0

, we obtain

J_{3} (v_{0}) = min_{v \in U_{3}} J_{3} (v) .

Finally, from the approximation indicated in the last remark and

u_{0} \in V_{1}

we get

δ J (u_{0}) = 0 + O (0.25) .

The proof is complete. □

References

R.A. Adams and J.F. Fournier, Sobolev Spaces, 2nd edn. (Elsevier, New York, 2003).
W.R. Bielski, A. Galka, J.J. Telega, The Complementary Energy Principle and Duality for Geometrically Nonlinear Elastic Shells. I. Simple case of moderate rotations around a tangent to the middle surface. Bulletin of the Polish Academy of Sciences, Technical Sciences, Vol. 38, No. 7-9, 1988.
W.R. Bielski and J.J. Telega, A Contribution to Contact Problems for a Class of Solids and Structures, Arch. Mech., 37, 4-5, pp. 303-320, Warszawa 1985.
J.F. Annet, Superconductivity, Superfluids and Condensates, 2nd edn. ( Oxford Master Series in Condensed Matter Physics, Oxford University Press, Reprint, 2010).
F.S. Botelho, Functional Analysis, Calculus of Variations and Numerical Methods in Physics and Engineering, CRC Taylor and Francis, Florida, 2020.
F.S. Botelho, Variational Convex Analysis, Ph.D. thesis, Virginia Tech, Blacksburg, VA -USA, (2009).
F. Botelho, Topics on Functional Analysis, Calculus of Variations and Duality, Academic Publications, Sofia, (2011).
F. Botelho, Existence of solution for the Ginzburg-Landau system, a related optimal control problem and its computation by the generalized method of lines, Applied Mathematics and Computation, 218, 11976-11989, (2012). [CrossRef]
F. Botelho, Functional Analysis and Applied Optimization in Banach Spaces, Springer Switzerland, 2014. [CrossRef]
P.Ciarlet, Mathematical Elasticity, Vol. II – Theory of Plates, North Holland Elsevier (1997).
J.C. Strikwerda, Finite Difference Schemes and Partial Differential Equations, SIAM, second edition (Philadelphia, 2004).
L.D. Landau and E.M. Lifschits, Course of Theoretical Physics, Vol. 5- Statistical Physics, part 1. (Butterworth-Heinemann, Elsevier, reprint 2008).
R.T. Rockafellar, Convex Analysis, Princeton Univ. Press, (1970).
J.J. Telega, On the complementary energy principle in non-linear elasticity. Part I: Von Karman plates and three dimensional solids, C.R. Acad. Sci. Paris, Serie II, 308, 1193-1198; Part II: Linear elastic solid and non-convex boundary condition. Minimax approach, ibid, pp. 1313-1317 (1989).
A.Galka and J.J.Telega Duality and the complementary energy principle for a class of geometrically non-linear structures. Part I. Five parameter shell model; Part II. Anomalous dual variational priciples for compressed elastic beams, Arch. Mech. 47 (1995) 677-698, 699-724.
J.F. Toland, A duality principle for non-convex optimisation and the calculus of variations, Arch. Rat. Mech. Anal., 71, No. 1 (1979), 41-61. [CrossRef]

Figure 1. solution

u_{0} (x)

for the case

f (x) = 0

Figure 1. solution

u_{0} (x)

for the case

f (x) = 0

Figure 3. solution

u_{0} (x)

for the case

f (x) = 0

Figure 3. solution

u_{0} (x)

for the case

f (x) = 0

Figure 4. solution

u_{0} (x)

for the case

f (x) = sin (π x) / 2

Figure 4. solution

u_{0} (x)

for the case

f (x) = sin (π x) / 2

Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

© 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.

Submitted:

16 February 2023

Posted:

20 February 2023

Read the latest preprint version here

Alerts

Abstract

Keywords:

Subject: Computer Science and Mathematics - Applied Mathematics

1. Introduction

In this section we establish a dual formulation for a large class of models in non-convex optimization.

The main duality principle is applied to double well models similar as those found in the phase transition theory.

Such results are based on the works of J.J. Telega and W.R. Bielski [2,3,14,15] and on a D.C. optimization approach developed in Toland [16].

About the other references, details on the Sobolev spaces involved are found in [1]. Related results on convex analysis and duality theory are addressed in [5,6,7,9,13].

Finally, in this text we adopt the standard Einstein convention of summing up repeated indices, unless otherwise indicated.

In order to clarify the notation, here we introduce the definition of topological dual space.

Definition 1.1

U^{*}

, through a bilinear form

{〈 \cdot, \cdot 〉}_{U} : U \times U^{*} \to R

(here we are referring to standard representations of dual spaces of Sobolev and Lebesgue spaces). Thus, given

f : U \to R

linear and continuous, we assume the existence of a unique

u^{*} \in U^{*}

such that

\begin{matrix} f (u) = {〈 u, u^{*} 〉}_{U}, \forall u \in U . \end{matrix}

(1)

The norm of f , denoted by

{∥ f ∥}_{U^{*}}

, is defined as

\begin{matrix} {∥ f ∥}_{U^{*}} = sup_{u \in U} {| {〈 u, u^{*} 〉}_{U} {| : ∥ u ∥}_{U} \leq 1} \equiv ∥ u^{*} ∥_{U^{*}} . \end{matrix}

(2)

At this point we start to describe the primal and dual variational formulations.

2. A General Duality Principle Non-Convex Optimization

In this section we present a duality principle applicable to a model in phase transition.

This case corresponds to the vectorial one in the calculus of variations.

Let

Ω \subset R^{n}

be an open, bounded, connected set with a regular (Lipschitzian) boundary denoted by

\partial Ω .

Consider a functional

J : V \to R

where

J (u) = F (\nabla u_{1}, \dots, \nabla u_{N}) + G (u_{1}, \dots, u_{N}) - {〈 u_{i}, f_{i} 〉}_{L^{2}},

and where

V = {u = (u_{1}, \dots, u_{N}) \in W^{1, p} (Ω; R^{N}) : u = u_{0} on \partial Ω},

f \in L^{2} (Ω; R^{N}),

and

1 < p < + \infty .

We assume there exists

α \in R

such that

α = inf_{u \in V} J (u) .

Anyway, one question remains, how the minimizing sequences behave close the infimum of J.

We intend to use duality theory to approximately solve such a global optimization problem.

Denoting

V_{0} = W_{0}^{1, p} (Ω; R^{N})

Y_{1} = Y_{1}^{*} = L^{2} (Ω; R^{N \times n})

Y_{2} = Y_{2}^{*} = L^{2} (Ω; R^{N \times n})

Y_{3} = Y_{3}^{*} = L^{2} (Ω; R^{N}),

at this point we define,

F_{1} : V \times V_{0} \to R

G_{1} : V \to R

G_{2} : V \to R

G_{3} : V_{0} \to R

and

G_{4} : V \to R,

\begin{matrix} F_{1} (\nabla u, \nabla ϕ) & = & F (\nabla u_{1} + \nabla ϕ_{1}, \dots, \nabla u_{N} + \nabla ϕ_{N}) + \frac{K}{2} \int_{Ω} \nabla u_{j} \cdot \nabla u_{j} d x \\ + \frac{K_{2}}{2} \int_{Ω} \nabla ϕ_{j} \cdot \nabla ϕ_{j} d x \end{matrix}

(3)

and

G_{1} (u_{1}, \dots, u_{n}) = G (u_{1}, \dots, u_{N}) + \frac{K_{1}}{2} \int_{Ω} u_{j} u_{j} d x - {〈 u_{i}, f_{i} 〉}_{L^{2}},

G_{2} (\nabla u_{1}, \dots, \nabla u_{N}) = \frac{K_{1}}{2} \int_{Ω} \nabla u_{j} \cdot \nabla u_{j} d x,

G_{3} (\nabla ϕ_{1}, \dots, \nabla ϕ_{N}) = \frac{K_{2}}{2} \int_{Ω} \nabla ϕ_{j} \cdot \nabla ϕ_{j} d x,

and

G_{4} (u_{1}, \dots, u_{N}) = \frac{K_{1}}{2} \int_{Ω} u_{j} u_{j} d x .

Define now

J_{1} : V \times V_{0} \to R

J_{1} (u, ϕ) = F (\nabla u + \nabla ϕ) + G (u) - {〈 u_{i}, f_{i} 〉}_{L^{2}} .

Observe that

\begin{matrix} J_{1} (u, ϕ) & = & F_{1} (\nabla u, \nabla ϕ) + G_{1} (u) - G_{2} (\nabla u) - G_{3} (\nabla ϕ) - G_{4} (u) \\ \leq & F_{1} (\nabla u, \nabla ϕ) + G_{1} (u) - {〈 \nabla u, z_{1}^{*} 〉}_{L^{2}} - {〈 \nabla ϕ, z_{2}^{*} 〉}_{L^{2}} - {〈 u, z_{3}^{*} 〉}_{L^{2}} \\ + sup_{v_{1} \in Y_{1}} {{〈 v_{1}, z_{1}^{*} 〉}_{L^{2}} - G_{2} (v_{1})} \\ + sup_{v_{2} \in Y_{2}} {{〈 v_{2}, z_{2}^{*} 〉}_{L^{2}} - G_{3} (v_{2})} \\ + sup_{u \in V} {{〈 u, z_{3}^{*} 〉}_{L^{2}} - G_{4} (u)} \\ = & F_{1} (\nabla u, \nabla ϕ) + G_{1} (u) - {〈 \nabla u, z_{1}^{*} 〉}_{L^{2}} - {〈 \nabla ϕ, z_{2}^{*} 〉}_{L^{2}} - {〈 u, z_{3}^{*} 〉}_{L^{2}} \\ + G_{2}^{*} (z_{1}^{*}) + G_{3}^{*} (z_{2}^{*}) + G_{4}^{*} (z_{3}^{*}) \\ = & J_{1}^{*} (u, ϕ, z^{*}), \end{matrix}

(4)

\forall u \in V, ϕ \in V_{0}, z^{*} = (z_{1}^{*}, z_{2}^{*}, z_{3}^{*}) \in Y^{*} = Y_{1}^{*} \times Y_{2}^{*} \times Y_{3}^{*}

Here we assume

K, K_{1}, K_{2}

are large enough so that

F_{1}

and

G_{1}

are convex.

Hence, from the general results in [16], we may infer that

\begin{matrix} inf_{(u, ϕ) \in V \times V_{0}} J (u, ϕ) & = & inf_{(u, ϕ, z^{*}) \in V \times V_{0} \times Y^{*}} J_{1}^{*} (u, ϕ, z^{*}) . \end{matrix}

(5)

On the other hand

inf_{u \in V} J (u) \geq inf_{(u, ϕ) \in V \times V_{0}} J_{1} (u, ϕ) \geq inf_{u \in V} Q_{J} (u) = inf_{u \in V} J (u),

where

Q_{J} (u)

refers to a standard quasi-convex regularization of J.

From these last two results we may obtain

inf_{u \in V} J (u) = inf_{(u, ϕ, z^{*}) \in V \times V_{0} \times Y^{*}} J_{1}^{*} (u, ϕ, z^{*}) .

Moreover, from standards results on convex analysis, we may have

\begin{matrix} inf_{u \in V} J_{1}^{*} (u, ϕ, z^{*}) & = & inf_{u \in V} {F_{1} (\nabla u, \nabla ϕ) + G_{1} (u) \\ - {〈 \nabla u, z_{1}^{*} 〉}_{L^{2}} - {〈 \nabla ϕ, z_{2}^{*} 〉}_{L^{2}} - {〈 u, z_{3}^{*} 〉}_{L^{2}} \\ + G_{2}^{*} (z_{1}^{*}) + G_{3}^{*} (z_{2}^{*}) + G_{4}^{*} (z_{3}^{*})} \\ = & sup_{(v_{1}^{*}, v_{2}^{*}) \in C^{*}} {- F_{1}^{*} (v_{1}^{*} + z_{1}^{*}, \nabla ϕ) - G_{1}^{*} (v_{2}^{*} + z_{3}^{*}) - {〈 \nabla ϕ, z_{2}^{*} 〉}_{L^{2}} \\ + G_{2}^{*} (z_{1}^{*}) + G_{3}^{*} (z_{2}^{*}) + G_{4}^{*} (z_{3}^{*})}, \end{matrix}

(6)

where

C^{*} = {v^{*} = (v_{1}^{*}, v_{2}^{*}) \in Y_{1}^{*} \times Y_{3}^{*} : - div {(v_{1}^{*})}_{i} + {(v_{2}^{*})}_{i} = 0, \forall i \in {1, \dots, N}},

F_{1}^{*} (v_{1}^{*} + z_{1}^{*}, \nabla ϕ) = sup_{v_{1} \in Y_{1}} {{〈 v_{1}, z_{1}^{*} + v_{1}^{*} 〉}_{L^{2}} - F_{1} (v_{1}, \nabla ϕ)},

and

G_{1}^{*} (v_{2}^{*} + z_{2}^{*}) = sup_{u \in V} {{〈 u, v_{2}^{*} + z_{2}^{*} 〉}_{L^{2}} - G_{1} (u)} .

Thus, defining

J_{2}^{*} (ϕ, z^{*}, v^{*}) = F_{1}^{*} (v_{1}^{*} + z_{1}^{*}, \nabla ϕ) - G_{1}^{*} (v_{2}^{*} + z_{3}^{*}) - {〈 \nabla ϕ, z_{2}^{*} 〉}_{L^{2}} + G_{2}^{*} (z_{1}^{*}) + G_{3}^{*} (z_{2}^{*}) + G_{4}^{*} (z_{3}^{*}),

we have got

\begin{matrix} inf_{u \in V} J (u) & = & inf_{(u, ϕ) \in V \times V_{0}} J_{1} (u, ϕ) \\ = & inf_{(u, ϕ, z^{*}) \in V \times V_{0} \times Y^{*}} J_{1}^{*} (u, ϕ, z^{*}) \\ = & inf_{z^{*} \in Y^{*}} \{inf_{ϕ \in V_{0}} \{sup_{v^{*} \in C^{*}} J_{2}^{*} (ϕ, z^{*}, v^{*})\}\} . \end{matrix}

(7)

Finally, observe that

\begin{matrix} inf_{u \in V} J (u) \\ = & inf_{z^{*} \in Y^{*}} \{inf_{ϕ \in V_{0}} \{sup_{v^{*} \in C^{*}} J_{2}^{*} (ϕ, z^{*}, v^{*})\}\} \\ \geq & sup_{v^{*} \in C^{*}} \{inf_{(z^{*}, ϕ) \in Y^{*} \times V_{0}} J_{2}^{*} (ϕ, z^{*}, v^{*})\} . \end{matrix}

(8)

This last variational formulation corresponds to a concave relaxed formulation in

v^{*}

concerning the original primal formulation.

3. Another Duality Principle for a Simpler Related Model in Phase Transition with a Respective Numerical Example

In this section we present another duality principle for a related model in phase transition.

Let

Ω = [0, 1] \subset R

and consider a functional

J : V \to R

where

J (u) = \frac{1}{2} \int_{Ω} {({(u^{'})}^{2} - 1)}^{2} d x + \frac{1}{2} \int_{Ω} u^{2} d x - {〈 u, f 〉}_{L^{2}},

and where

V = {u \in W^{1, 4} (Ω) : u (0) = 0 and u (1) = 1 / 2}

and

f \in L^{2} (Ω) .

A global optimum point is not attained for J so that the problem of finding a global minimum for J has no solution.

Anyway, one question remains, how the minimizing sequences behave close the infimum of J.

We intend to use duality theory to approximately solve such a global optimization problem.

Denoting

V_{0} = W_{0}^{1, 4} (Ω)

, at this point we define,

F : V \to R

and

F_{1} : V \times V_{0} \to R

F (u) = \frac{1}{2} \int_{Ω} {({(u^{'})}^{2} - 1)}^{2} d x,

and

F_{1} (u, ϕ) = \frac{1}{2} \int_{Ω} {({(u^{'} + ϕ^{'})}^{2} - 1)}^{2} d x .

Observe

F (u) \geq inf_{ϕ \in V_{0}} F_{1} (u, ϕ) \geq Q_{F} (u), \forall u \in V,

where

Q_{F} (u)

refers to a quasi-convex regularization of

F .

We define also

F_{2} : V \times V_{0} \to R,

F_{3} : V \times V_{0} \to R

and

G : V \times V_{0} \to R

F_{2} (u, ϕ) = \frac{1}{2} \int_{Ω} {({(u^{'} + ϕ^{'})}^{2} - 1)}^{2} d x + \frac{1}{2} \int_{Ω} u^{2} d x - {〈 u, f 〉}_{L^{2}},

\begin{matrix} F_{3} (u, ϕ) & = & F_{2} (u, ϕ) + \frac{K}{2} \int_{Ω} {(u^{'})}^{2} d x \\ + \frac{K_{1}}{2} \int_{Ω} {(ϕ^{'})}^{2} d x \end{matrix}

(9)

and

\begin{matrix} G (u, ϕ) & = & \frac{K}{2} \int_{Ω} {(u^{'})}^{2} d x \\ + \frac{K_{1}}{2} \int_{Ω} {(ϕ^{'})}^{2} d x \end{matrix}

(10)

Observe that if

K > 0, K_{1} > 0

is large enough, both

F_{3}

and G are convex.

Denoting

Y = Y^{*} = L^{2} (Ω)

we also define the polar functional

G^{*} : Y^{*} \times Y^{*} \to R

G^{*} (v^{*}, v_{0}^{*}) = sup_{(u, ϕ) \in V \times V_{0}} {{〈 u, v^{*} 〉}_{L^{2}} + {〈 ϕ, v_{0}^{*} 〉}_{L^{2}} - G (u, ϕ)} .

Observe that

inf_{u \in U} J (u) \geq inf_{((u, ϕ), (v^{*}, v_{0}^{*})) \in V \times V_{0} \times {[Y^{*}]}^{2}} {G^{*} (v^{*}, v_{0}^{*}) - {〈 u, v^{*} 〉}_{L^{2}} - {〈 ϕ, v_{0}^{*} 〉}_{L^{2}} + F_{3} (u, ϕ)} .

With such results in mind, we define a relaxed primal dual variational formulation for the primal problem, represented by

J_{1}^{*} : V \times V_{0} \times {[Y^{*}]}^{2} \to R

, where

J_{1}^{*} (u, ϕ, v^{*}, v_{0}^{*}) = G^{*} (v^{*}, v_{0}^{*}) - {〈 u, v^{*} 〉}_{L^{2}} - {〈 ϕ, v_{0}^{*} 〉}_{L^{2}} + F_{3} (u, ϕ) .

Having defined such a functional, we may obtain numerical results by solving a sequence of convex auxiliary sub-problems, through the following algorithm.

Set $K \approx 150$ and $K_{1} = K / 20$ and $0 < ε ≪ 1 .$
Choose $(u_{1}, ϕ_{1}) \in V \times V_{0},$ such that $∥ u_{1} ∥_{1, \infty} ≪ K / 4$ and $∥ ϕ_{1} ∥_{1, \infty} ≪ K / 4 .$
Set $n = 1 .$
Calculate $(v_{n}^{*}, {(v_{0}^{*})}_{n})$ solution of the system of equations:

$\frac{\partial J_{1}^{*} (u_{n}, ϕ_{n}, v_{n}^{*}, {(v_{0}^{*})}_{n})}{\partial v^{*}} = 0$

and

$\frac{\partial J_{1}^{*} (u_{n}, ϕ_{n}, v_{n}^{*}, {(v_{0}^{*})}_{n})}{\partial v_{0}^{*}} = 0,$

that is

$\frac{\partial G^{*} (v_{n}^{*}, {(v_{0}^{*})}_{n})}{\partial v^{*}} - u_{n} = 0$

and

$\frac{\partial G^{*} (v_{n}^{*}, {(v_{0}^{*})}_{n})}{\partial v_{0}^{*}} - ϕ_{n} = 0$

so that

$v_{n}^{*} = \frac{\partial G (u_{n}, ϕ_{n})}{\partial u}$

and

${(v_{0}^{*})}_{n}^{*} = \frac{\partial G (u_{n}, ϕ_{n})}{\partial ϕ}$
Calculate $(u_{n + 1}, ϕ_{n + 1})$ by solving the system of equations:

$\frac{\partial J_{1}^{*} (u_{n + 1}, ϕ_{n + 1}, v_{n}^{*}, {(v_{0}^{*})}_{n})}{\partial u} = 0$

and

$\frac{\partial J_{1}^{*} (u_{n + 1}, ϕ_{n + 1}, v_{n}^{*}, {(v_{0}^{*})}_{n})}{\partial ϕ} = 0$

that is

$- v_{n}^{*} + \frac{\partial F_{3} (u_{n + 1}, ϕ_{n + 1})}{\partial u} = 0$

and

$- {(v_{0}^{*})}_{n} + \frac{\partial F_{3} (u_{n + 1}, ϕ_{n + 1})}{\partial ϕ} = 0$
If $max {∥ u_{n} - u_{n + 1} ∥_{\infty}, ∥ ϕ_{n + 1} - ϕ_{n} ∥_{\infty}} \leq ε$ , then stop, else set $n : = n + 1$ and go to item 4.

For the case in which

f (x) = 0

, we have obtained numerical results for

K = 1500

and

K_{1} = K / 20

. For such a concerning solution

u_{0}

obtained, please see Figure 1. For the case in which

f (x) = sin (π x) / 2

, we have obtained numerical results for

K = 100

and

K_{1} = K / 20

. For such a concerning solution

u_{0}

obtained, please see Figure 2.

Remark 3.1.

4. A Convex Dual Variational Formulation for a Third Similar Model

In this section we present another duality principle for a third related model in phase transition.

Figure 2. solution

u_{0} (x)

for the case

f (x) = sin (π x) / 2

Figure 2. solution

u_{0} (x)

for the case

f (x) = sin (π x) / 2

Let

Ω = [0, 1] \subset R

and consider a functional

J : V \to R

where

J (u) = \frac{1}{2} \int_{Ω} min {{(u^{'} - 1)}^{2}, {(u^{'} + 1)}^{2}} d x + \frac{1}{2} \int_{Ω} u^{2} d x - {〈 u, f 〉}_{L^{2}},

and where

V = {u \in W^{1, 2} (Ω) : u (0) = 0 and u (1) = 1 / 2}

and

f \in L^{2} (Ω) .

A global optimum point is not attained for J so that the problem of finding a global minimum for J has no solution.

Anyway, one question remains, how the minimizing sequences behave close to the infimum of J.

We intend to use the duality theory to solve such a global optimization problem in an appropriate sense to be specified.

At this point we define,

F : V \to R

and

G : V \to R

\begin{matrix} F (u) & = & \frac{1}{2} \int_{Ω} min {{(u^{'} - 1)}^{2}, {(u^{'} + 1)}^{2}} d x \\ = & \frac{1}{2} \int_{Ω} {(u^{'})}^{2} d x - \int_{Ω} | u^{'} | d x + 1 / 2 \\ \equiv & F_{1} (u^{'}), \end{matrix}

(11)

and

G (u) = \frac{1}{2} \int_{Ω} u^{2} d x - {〈 u, f 〉}_{L^{2}} .

Denoting

Y = Y^{*} = L^{2} (Ω)

we also define the polar functional

F_{1}^{*} : Y^{*} \to R

and

G^{*} : Y^{*} \to R

\begin{matrix} F_{1}^{*} (v^{*}) & = & sup_{v \in Y} {{〈 v, v^{*} 〉}_{L^{2}} - F_{1} (v)} \\ = & \frac{1}{2} \int_{Ω} {(v^{*})}^{2} d x + \int_{Ω} | v^{*} | d x, \end{matrix}

(12)

and

\begin{matrix} G^{*} ({(v^{*})}^{'}) & = & sup_{u \in V} {- {〈 u^{'}, v^{*} 〉}_{L^{2}} - G (u)} \\ = & \frac{1}{2} \int_{Ω} {({(v^{*})}^{'} + f)}^{2} d x - \frac{1}{2} v^{*} (1) . \end{matrix}

(13)

Observe this is the scalar case of the calculus of variations, so that from the standard results on convex analysis, we have

inf_{u \in V} J (u) = max_{v^{*} \in Y^{*}} {- F_{1}^{*} (v^{*}) - G^{*} (- {(v^{*})}^{'})} .

Indeed, from the direct method of the calculus of variations, the maximum for the dual formulation is attained at some

{\hat{v}}^{*} \in Y^{*}

Moreover, the corresponding solution

u_{0} \in V

is obtained from the equation

u_{0} = \frac{\partial G ({({\hat{v}}^{*})}^{'})}{\partial {(v^{*})}^{'}} = {({\hat{v}}^{*})}^{'} + f .

Finally, the Euler-Lagrange equations for the dual problem stands for

\{\begin{matrix} {(v^{*})}^{″} + f^{'} - v^{*} - sign (v^{*}) = 0, & in Ω, \\ {(v^{*})}^{'} (0) = 0, {(v^{*})}^{'} (1) = 1 / 2, \end{matrix}

(14)

where

sign (v^{*} (x)) = 1

v^{*} (x) > 0,

sign (v^{*} (x)) = - 1,

v^{*} (x) < 0

and

- 1 \leq sign (v^{*} (x)) \leq 1,

v^{*} (x) = 0 .

We have computed the solutions

v^{*}

and corresponding solutions

u_{0} \in V

for the cases in which

f (x) = 0

and

f (x) = sin (π x) / 2 .

For the solution

u_{0} (x)

for the case in which

f (x) = 0

, please see Figure 3.

For the solution

u_{0} (x)

for the case in which

f (x) = sin (π x) / 2

, please see Figure 4.

Remark 4.1.

Observe that such solutions

u_{0}

obtained are not the global solutions for the related primal optimization problems. Indeed, such solutions reflect the average behavior of weak cluster points for concerning minimizing sequences.

4.1. The Algorithm through Which We Have Obtained the Numerical Results

In this subsection we present the software in MATLAB through which we have obtained the last numerical results.

This algorithm is for solving the concerning Euler-Lagrange equations for the dual problem, that is, for solving the equation

\{\begin{matrix} {(v^{*})}^{″} + f^{'} - v^{*} - sign (v^{*}) = 0, & in Ω, \\ {(v^{*})}^{'} (0) = 0, {(v^{*})}^{'} (1) = 1 / 2 . \end{matrix}

(15)

Here the concerning software in MATLAB. We emphasize to have used the smooth approximation

| v^{*} | \approx \sqrt{{(v^{*})}^{2} + e_{1}},

where a small value for

e_{1}

is specified in the next lines.

*************************************

clear all
$m_{8} = 800;$ (number of nodes)
$d = 1 / m_{8};$
$e_{1} = 0.00001;$
$f o r i = 1 : m_{8}$

$y o (i, 1) = 0.01;$

$y_{1} (i, 1) = sin (π * i / m_{8}) / 2;$

$e n d;$
$f o r i = 1 : m_{8} - 1$

$d y_{1} (i, 1) = (y_{1} (i + 1, 1) - y_{1} (i, 1)) / d;$

$e n d;$
$f o r k = 1 : 3000$ (we have fixed the number of iterations)

$i = 1;$

$h_{3} = 1 / \sqrt{v o {(i, 1)}^{2} + e_{1}};$

$m_{12} = 1 + d^{2} * h_{3} + d^{2};$

$m_{50} (i) = 1 / m_{12};$

$z (i) = m_{50} (i) * (d y_{1} (i, 1) * d^{2});$
$f o r i = 2 : m_{8} - 1$

$h_{3} = 1 / \sqrt{v o {(i, 1)}^{2} + e_{1}};$

$m_{12} = 2 + h_{3} * d^{2} + d^{2} - m 50 (i - 1);$

$m 50 (i) = 1 / m_{12};$

$z (i) = m_{50} (i) * (z (i - 1) + d y_{1} (i, 1) * d^{2});$

$e n d;$
$v (m_{8}, 1) = (d / 2 + z (m_{8} - 1)) / (1 - m_{50} (m_{8} - 1));$
$f o r i = 1 : m_{8} - 1$

$v (m_{8} - i, 1) = m_{50} (m_{8} - i) * v (m_{8} - i + 1) + z (m_{8} - i);$

$e n d;$
$v (m_{8} / 2, 1)$
$v o = v;$

$e n d;$
$f o r i = 1 : m_{8} - 1$

$u (i, 1) = (v (i + 1, 1) - v (i, 1)) / d + y_{1} (i, 1);$

$e n d;$
$f o r i = 1 : m_{8} - 1$

$x (i) = i * d;$

$e n d;$

$p l o t (x, u (:, 1))$

********************************

5. An Improvement of the Convexity Conditions for a Non-Convex Related Model through an Approximate Primal Formulation

In this section we develop an approximate primal dual formulation suitable for a large class of variational models.

Here, the applications are for the Kirchhoff-Love plate model, which may be found in Ciarlet, [10].

At this point we start to describe the primal variational formulation.

Let

Ω \subset R^{2}

be an open, bounded, connected set which represents the middle surface of a plate of thickness h. The boundary of

Ω

, which is assumed to be regular (Lipschitzian), is denoted by

\partial Ω

. The vectorial basis related to the cartesian system

{x_{1}, x_{2}, x_{3}}

is denoted by

(a_{α}, a_{3})

, where

α = 1, 2

(in general Greek indices stand for 1 or 2), and where

a_{3}

is the vector normal to

Ω

, whereas

a_{1}

and

a_{2}

are orthogonal vectors parallel to

Ω .

Also,

n

is the outward normal to the plate surface.

The displacements will be denoted by

\hat{u} = {{\hat{u}}_{α}, {\hat{u}}_{3}} = {\hat{u}}_{α} a_{α} + {\hat{u}}_{3} a_{3} .

The Kirchhoff-Love relations are

\begin{matrix} {\hat{u}}_{α} (x_{1}, x_{2}, x_{3}) = u_{α} (x_{1}, x_{2}) - x_{3} w {(x_{1}, x_{2})}_{, α} \\ and {\hat{u}}_{3} (x_{1}, x_{2}, x_{3}) = w (x_{1}, x_{2}) . \end{matrix}

(16)

Here

- h / 2 \leq x_{3} \leq h / 2

so that we have

u = (u_{α}, w) \in U

where

\begin{matrix} U & = & \{u = (u_{α}, w) \in W^{1, 2} (Ω; R^{2}) \times W^{2, 2} (Ω), \\ u_{α} = w = \frac{\partial w}{\partial n} = 0 on \partial Ω\} \\ = & W_{0}^{1, 2} (Ω; R^{2}) \times W_{0}^{2, 2} (Ω) . \end{matrix}

It is worth emphasizing that the boundary conditions here specified refer to a clamped plate.

We also define the operator

Λ : U \to Y \times Y

, where

Y = Y^{*} = L^{2} (Ω; R^{2 \times 2})

, by

Λ (u) = {γ (u), κ (u)},

γ_{α β} (u) = \frac{u_{α, β} + u_{β, α}}{2} + \frac{w_{, α} w_{, β}}{2},

κ_{α β} (u) = - w_{, α β} .

The constitutive relations are given by

N_{α β} (u) = H_{α β λ μ} γ_{λ μ} (u),

(17)

M_{α β} (u) = h_{α β λ μ} κ_{λ μ} (u),

(18)

where:

{H_{α β λ μ}}

and

\{h_{α β λ μ} = \frac{h^{2}}{12} H_{α β λ μ}\}

, are symmetric positive definite fourth order tensors. From now on, we denote

{{\bar{H}}_{α β λ μ}} = {H_{α β λ μ}}^{- 1}

and

{{\bar{h}}_{α β λ μ}} = {h_{α β λ μ}}^{- 1}

Furthermore

{N_{α β}}

denote the membrane force tensor and

{M_{α β}}

the moment one. The plate stored energy, represented by

(G \circ Λ) : U \to R

is expressed by

(G \circ Λ) (u) = \frac{1}{2} \int_{Ω} N_{α β} (u) γ_{α β} (u) d x + \frac{1}{2} \int_{Ω} M_{α β} (u) κ_{α β} (u) d x

(19)

and the external work, represented by

F : U \to R

, is given by

F (u) = {〈 w, P 〉}_{L^{2}} + {〈 u_{α}, P_{α} 〉}_{L^{2}},

(20)

where

P, P_{1}, P_{2} \in L^{2} (Ω)

are external loads in the directions

a_{3}

a_{1}

and

a_{2}

respectively. The potential energy, denoted by

J : U \to R

is expressed by:

J (u) = (G \circ Λ) (u) - F (u)

Define now

J_{3} : \tilde{U} \to R

J_{3} (u) = J (u) + 10 \int_{Ω} \frac{e^{K b w}}{K^{3 / 2}} d x .

In such a case for

b = P / | P |

Ω

K = 100

and

\tilde{U} = {{u \in U : ∥ w ∥}_{\infty} \leq 0.01 and P w \geq 0 a . e . in Ω},

we get

\begin{matrix} \frac{\partial J_{3} (w)}{\partial w} & = & \frac{\partial J (w)}{\partial w} + 10 b \frac{e^{K b w}}{K^{1 / 2}} \\ \approx & \frac{\partial J (w)}{\partial w} + O (1.5), \end{matrix}

(21)

and

\begin{matrix} \frac{\partial^{2} J_{3} (w)}{\partial w^{2}} & = & \frac{\partial^{2} J (w)}{\partial w^{2}} + 10 b^{2} e^{K b w} K^{1 / 2} \\ \approx & \frac{\partial^{2} J (w)}{\partial w^{2}} + O (50) . \end{matrix}

(22)

This new functional

J_{3}

has a relevant improvement in the convexity conditions concerning the previous functional J.

Indeed, we have obtained a gain in positiveness for the second variation

\frac{\partial^{2} J (u)}{\partial w^{2}},

which has increased of order

O (50) .

Moreover the difference between the approximate and exact equation

\frac{\partial J (u)}{\partial w} = 0

is of order

O (1.5)

which corresponds to a small perturbation in the original equation for a load of

P = 1500 N / m^{2},

for example. Summarizing, the exact equation may be approximate solved in an appropriate sense.

6. An Approximate Convex Variational Formulation for Another Related Model

In this section, we obtain an approximate convex variational formulation for a related model, more specifically, for a Ginzburg-Landau type equation.

Let

Ω \subset R^{3}

be an open, bounded and connected set with a regular (Lipschitzian) boundary denoted by

\partial Ω .

Consider a functional

J : V \to R

where

\begin{matrix} J (u) & = & \frac{γ}{2} \int_{Ω} \nabla u \cdot \nabla u d x + \frac{α}{2} \int_{Ω} {(u^{2} - β)}^{2} d x \\ - {〈 u, f 〉}_{L^{2}}, \end{matrix}

(23)

where

γ > 0

α > 0

β > 0

V = W_{0}^{1, 2} (Ω)

and

f \in L^{2} (Ω) .

We define

A^{+} = {u \in V : u f \geq 0, a . e . in Ω},

V_{2} = {{u \in V : ∥ u ∥}_{\infty} \leq 1},

and

V_{1} = V_{2} \cap A^{+} .

At this point we define

v = u / 10

so that

\begin{matrix} J (u) & = & J_{1} (v) \\ = & \frac{10^{2} γ}{2} \int_{Ω} \nabla v \cdot \nabla v d x + \frac{α}{2} \int_{Ω} ({(10 v)}^{2} - β) d x \\ - {〈 10 v, f 〉}_{L^{2}} . \end{matrix}

(24)

Moreover we define

\begin{matrix} J_{2} (v) & = & \frac{1}{10} J_{1} (v) \\ = & \frac{10 γ}{2} \int_{Ω} \nabla v \cdot \nabla v d x + \frac{α}{20} \int_{Ω} ({(10 v)}^{2} - β) d x \\ - {〈 v, f 〉}_{L^{2}}, \end{matrix}

(25)

and

J_{3} : U_{3} \to R

where

J_{3} (v) = J_{2} (v) + K_{1} \int_{Ω} \frac{a^{(K b v / 300)}}{ln (a) K^{3 / 2}} 300 d x,

where

K_{1} = 50 / 12

a = 2.71

K = 4755

b = f / | f |

Ω

and

U_{2} = {v \in V : f v \geq 0, a . e . in Ω},

U_{4} = {{v \in V : ∥ v ∥}_{\infty} \leq 1 / 10},

and

U_{3} = U_{2} \cap U_{4} .

Thus, with such numerical values, we may obtain

\begin{matrix} \frac{\partial J_{3} (v)}{\partial v} & = & \frac{\partial J_{2} (v)}{\partial v} + 0.0145019 K_{1} b 2 . 71^{(317 b v / 2) 0} \\ \approx & \frac{\partial J_{2} (v)}{\partial v} + O (0.25), \end{matrix}

(26)

and

\begin{matrix} \frac{\partial^{2} J_{3} (v)}{\partial v^{2}} & = & \frac{\partial^{2} J_{2} (v)}{\partial v^{2}} + 0.229154 K_{1} b^{2} 2 . 71^{(317 b v / 20)} \\ \approx & \frac{\partial^{2} J_{2} (v)}{\partial v^{2}} + O (2.5) . \end{matrix}

(27)

Remark 6.1.

This new functional

J_{1}

has a relevant improvement in the convexity conditions concerning the previous functional J.

Indeed, we have obtained a gain in positiveness for the second variation

\frac{\partial^{2} J_{2} (v)}{\partial v^{2}},

which has increased of order

O (2.5) .

Moreover the difference between the approximate and exact equation

\frac{\partial J_{2} (v)}{\partial v} = 0

is of order

O (0.25)

which for appropriate parameters

γ > 0, α > 0

and

β > 0

, corresponds to a small perturbation in the original equation. Summarizing, the exact equation may be approximately solved in an appropriate sense.

Finally, for this last example, we highlight it is relatively easy to improve even more both such an approximation quality and the convexity conditions concerning the original variational model.

With such statements and results in mind, we may prove the following theorem.

Theorem 6.2.

Suppose

γ > 0,

α > 0

and

β > 0

are such that

\frac{\partial^{2} J_{3} (v)}{\partial v^{2}} > 0,

U_{3}

Assume also,

v_{0} \in U_{3}

is such that

δ J_{3} (v_{0}) = 0 .

Under such hypotheses,

J_{3}

is convex on

U_{3}

so that

J_{3} (v_{0}) = min_{v \in U_{3}} J_{3} (v) .

Moreover,

δ J (u_{0}) = 0 + O (0.25),

where

u_{0} = 10 v_{0} \in V_{1}

Proof.

From the hypotheses

\frac{\partial^{2} J_{3} (v)}{\partial v^{2}} > 0

U_{3}

, so that

J_{3}

is convex on the convex set

U_{3} .

Consequently, since

δ J_{3} (v_{0}) = 0

, we obtain

J_{3} (v_{0}) = min_{v \in U_{3}} J_{3} (v) .

Finally, from the approximation indicated in the last remark and

u_{0} \in V_{1}

we get

δ J (u_{0}) = 0 + O (0.25) .

The proof is complete. □

References

R.A. Adams and J.F. Fournier, Sobolev Spaces, 2nd edn. (Elsevier, New York, 2003).
W.R. Bielski, A. Galka, J.J. Telega, The Complementary Energy Principle and Duality for Geometrically Nonlinear Elastic Shells. I. Simple case of moderate rotations around a tangent to the middle surface. Bulletin of the Polish Academy of Sciences, Technical Sciences, Vol. 38, No. 7-9, 1988.
W.R. Bielski and J.J. Telega, A Contribution to Contact Problems for a Class of Solids and Structures, Arch. Mech., 37, 4-5, pp. 303-320, Warszawa 1985.
J.F. Annet, Superconductivity, Superfluids and Condensates, 2nd edn. ( Oxford Master Series in Condensed Matter Physics, Oxford University Press, Reprint, 2010).
F.S. Botelho, Functional Analysis, Calculus of Variations and Numerical Methods in Physics and Engineering, CRC Taylor and Francis, Florida, 2020.
F.S. Botelho, Variational Convex Analysis, Ph.D. thesis, Virginia Tech, Blacksburg, VA -USA, (2009).
F. Botelho, Topics on Functional Analysis, Calculus of Variations and Duality, Academic Publications, Sofia, (2011).
F. Botelho, Existence of solution for the Ginzburg-Landau system, a related optimal control problem and its computation by the generalized method of lines, Applied Mathematics and Computation, 218, 11976-11989, (2012). [CrossRef]
F. Botelho, Functional Analysis and Applied Optimization in Banach Spaces, Springer Switzerland, 2014. [CrossRef]
P.Ciarlet, Mathematical Elasticity, Vol. II – Theory of Plates, North Holland Elsevier (1997).
J.C. Strikwerda, Finite Difference Schemes and Partial Differential Equations, SIAM, second edition (Philadelphia, 2004).
L.D. Landau and E.M. Lifschits, Course of Theoretical Physics, Vol. 5- Statistical Physics, part 1. (Butterworth-Heinemann, Elsevier, reprint 2008).
R.T. Rockafellar, Convex Analysis, Princeton Univ. Press, (1970).
J.J. Telega, On the complementary energy principle in non-linear elasticity. Part I: Von Karman plates and three dimensional solids, C.R. Acad. Sci. Paris, Serie II, 308, 1193-1198; Part II: Linear elastic solid and non-convex boundary condition. Minimax approach, ibid, pp. 1313-1317 (1989).
A.Galka and J.J.Telega Duality and the complementary energy principle for a class of geometrically non-linear structures. Part I. Five parameter shell model; Part II. Anomalous dual variational priciples for compressed elastic beams, Arch. Mech. 47 (1995) 677-698, 699-724.
J.F. Toland, A duality principle for non-convex optimisation and the calculus of variations, Arch. Rat. Mech. Anal., 71, No. 1 (1979), 41-61. [CrossRef]

Figure 1. solution

u_{0} (x)

for the case

f (x) = 0

Figure 1. solution

u_{0} (x)

for the case

f (x) = 0

Figure 3. solution

u_{0} (x)

for the case

f (x) = 0

Figure 3. solution

u_{0} (x)

for the case

f (x) = 0

Figure 4. solution

u_{0} (x)

for the case

f (x) = sin (π x) / 2

Figure 4. solution

u_{0} (x)

for the case

f (x) = sin (π x) / 2

MDPI Initiatives

Important Links

Choose an area of interest and we will send you notifications of new preprints at your preferred frequency.

Disclaimer

Duality Principles and Numerical Procedures for a Large Class of Non-convex Models in the Calculus of Variations

Abstract

1. Introduction

2. A General Duality Principle Non-Convex Optimization

3. Another Duality Principle for a Simpler Related Model in Phase Transition with a Respective Numerical Example

4. A Convex Dual Variational Formulation for a Third Similar Model

4.1. The Algorithm through Which We Have Obtained the Numerical Results

5. An Improvement of the Convexity Conditions for a Non-Convex Related Model through an Approximate Primal Formulation

6. An Approximate Convex Variational Formulation for Another Related Model

References

Abstract

1. Introduction

2. A General Duality Principle Non-Convex Optimization

3. Another Duality Principle for a Simpler Related Model in Phase Transition with a Respective Numerical Example

4. A Convex Dual Variational Formulation for a Third Similar Model

4.1. The Algorithm through Which We Have Obtained the Numerical Results

5. An Improvement of the Convexity Conditions for a Non-Convex Related Model through an Approximate Primal Formulation

6. An Approximate Convex Variational Formulation for Another Related Model

References

MDPI Initiatives

Important Links

Subscribe