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

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 ${\langle ·,·\rangle }_U:U\times U^{*}\to \mathbb{R}$ (here we are referring to standard representations of dual spaces of Sobolev and Lebesgue spaces). Thus, given $f:U\to \mathbb{R}$ linear and continuous, we assume the existence of a unique $u^{*}\in U^{*}$ such that

$$\begin{array}{c}f\left(u\right)={\langle u,u^{*}\rangle }_U,\forall u\in U.\end{array}$$

(1)

The norm of f , denoted by ${\parallel f\parallel }_{U^{*}}$, is defined as

$$\begin{array}{c}{\parallel f\parallel }_{U^{*}}=\underset{u\in U}{sup}\left\{\right|{\langle u,u^{*}\rangle }_U{|:\parallel u\parallel }_U\le 1\}\equiv \parallel u^{*}{\parallel }_{U^{*}}.\end{array}$$

(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 $\mathsf{\Omega }\subset {\mathbb{R}}^n$ be an open, bounded, connected set with a regular (Lipschitzian) boundary denoted by $\partial \mathsf{\Omega }.$

Consider a functional $J:V\to \mathbb{R}$ where

$$J\left(u\right)=F(\nabla u_1,\cdots ,\nabla u_N)+G(u_1,\cdots ,u_N)-{\langle u_i,f_i\rangle }_{L^2},$$

and where

$$V=\{u=(u_1,\cdots ,u_N)\in W^{1,p}(\mathsf{\Omega };{\mathbb{R}}^N):u=u_0\mathrm{on}\partial \mathsf{\Omega }\},$$

$f\in L^2(\mathsf{\Omega };{\mathbb{R}}^N),$ and $1<p<+\infty .$

We assume there exists $\alpha \in \mathbb{R}$ such that

$$\alpha =\underset{u\in V}{inf}J\left(u\right).$$

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}(\mathsf{\Omega };{\mathbb{R}}^N)$, $Y_1=Y_1^{*}=L^2(\mathsf{\Omega };{\mathbb{R}}^{N\times n})$, $Y_2=Y_2^{*}=L^2(\mathsf{\Omega };{\mathbb{R}}^{N\times n})$, $Y_3=Y_3^{*}=L^2(\mathsf{\Omega };{\mathbb{R}}^N),$ at this point we define, $F_1:V\times V_0\to \mathbb{R}$, $G_1:V\to \mathbb{R}$, $G_2:V\to \mathbb{R}$, $G_3:V_0\to \mathbb{R}$ and $G_4:V\to \mathbb{R},$ by

$$\begin{array}{ccc}\hfill F_1(\nabla u,\nabla \varphi )& =& F(\nabla u_1+\nabla {\varphi }_1,\cdots ,\nabla u_N+\nabla {\varphi }_N)+\frac{K}{2}{\int }_{\mathsf{\Omega }}\nabla u_j·\nabla u_{j}dx\hfill \\ & & +\frac{K_2}{2}{\int }_{\mathsf{\Omega }}\nabla {\varphi }_j·\nabla {\varphi }_{j}dx\hfill \end{array}$$

(3)

and

$$G_1(u_1,\cdots ,u_n)=G(u_1,\cdots ,u_N)+\frac{K_1}{2}{\int }_{\mathsf{\Omega }}{u}_j{u}_{j}dx-{\langle u_i,f_i\rangle }_{L^2},$$

$$G_2(\nabla u_1,\cdots ,\nabla u_N)=\frac{K_1}{2}{\int }_{\mathsf{\Omega }}\nabla u_j·\nabla u_{j}dx,$$

$$G_3(\nabla {\varphi }_1,\cdots ,\nabla {\varphi }_N)=\frac{K_2}{2}{\int }_{\mathsf{\Omega }}\nabla {\varphi }_j·\nabla {\varphi }_{j}dx,$$

and

$$G_4(u_1,\cdots ,u_N)=\frac{K_1}{2}{\int }_{\mathsf{\Omega }}{u}_j{u}_{j}dx.$$

Define now $J_1:V\times V_0\to \mathbb{R}$,

$$J_1(u,\varphi )=F(\nabla u+\nabla \varphi )+G\left(u\right)-{\langle u_i,f_i\rangle }_{L^2}.$$

Observe that

$$\begin{array}{ccc}\hfill J_1(u,\varphi )& =& F_1(\nabla u,\nabla \varphi )+G_1\left(u\right)-G_2(\nabla u)-G_3(\nabla \varphi )-G_4\left(u\right)\hfill \\ & \le & F_1(\nabla u,\nabla \varphi )+G_1\left(u\right)-{\langle \nabla u,z_1^{*}\rangle }_{L^2}-{\langle \nabla \varphi ,z_2^{*}\rangle }_{L^2}-{\langle u,z_3^{*}\rangle }_{L^2}\hfill \\ & & +\underset{v_1\in Y_1}{sup}\{{\langle v_1,z_1^{*}\rangle }_{L^2}-G_2\left(v_1\right)\}\hfill \\ & & +\underset{v_2\in Y_2}{sup}\{{\langle v_2,z_2^{*}\rangle }_{L^2}-G_3\left(v_2\right)\}\hfill \\ & & +\underset{u\in V}{sup}\{{\langle u,z_3^{*}\rangle }_{L^2}-G_4\left(u\right)\}\hfill \\ & =& F_1(\nabla u,\nabla \varphi )+G_1\left(u\right)-{\langle \nabla u,z_1^{*}\rangle }_{L^2}-{\langle \nabla \varphi ,z_2^{*}\rangle }_{L^2}-{\langle u,z_3^{*}\rangle }_{L^2}\hfill \\ & & +G_2^{*}\left(z_1^{*}\right)+G_3^{*}\left(z_2^{*}\right)+G_4^{*}\left(z_3^{*}\right)\hfill \\ & =& J_1^{*}(u,\varphi ,z^{*}),\hfill \end{array}$$

(4)

$\forall u\in V,\varphi \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{array}{ccc}\hfill \underset{(u,\varphi )\in V\times V_0}{inf}J(u,\varphi )& =& \underset{(u,\varphi ,z^{*})\in V\times V_0\times Y^{*}}{inf}{J}_1^{*}(u,\varphi ,z^{*}).\hfill \end{array}$$

(5)

On the other hand

$$\underset{u\in V}{inf}J\left(u\right)\ge \underset{(u,\varphi )\in V\times V_0}{inf}{J}_1(u,\varphi )\ge \underset{u\in V}{inf}{Q}_J\left(u\right)=\underset{u\in V}{inf}J\left(u\right),$$

where $Q_J\left(u\right)$ refers to a standard quasi-convex regularization of J.

From these last two results we may obtain

$$\underset{u\in V}{inf}J\left(u\right)=\underset{(u,\varphi ,z^{*})\in V\times V_0\times Y^{*}}{inf}{J}_1^{*}(u,\varphi ,z^{*}).$$

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

$$\begin{array}{ccc}\hfill \underset{u\in V}{inf}{J}_1^{*}(u,\varphi ,z^{*})& =& \underset{u\in V}{inf}\{F_1(\nabla u,\nabla \varphi )+G_1\left(u\right)\hfill \\ & & -{\langle \nabla u,z_1^{*}\rangle }_{L^2}-{\langle \nabla \varphi ,z_2^{*}\rangle }_{L^2}-{\langle u,z_3^{*}\rangle }_{L^2}\hfill \\ & & +G_2^{*}\left(z_1^{*}\right)+G_3^{*}\left(z_2^{*}\right)+G_4^{*}\left(z_3^{*}\right)\}\hfill \\ & =& \underset{(v_1^{*},v_2^{*})\in C^{*}}{sup}\{-F_1^{*}(v_1^{*}+z_1^{*},\nabla \varphi )-G_1^{*}(v_2^{*}+z_3^{*})-{\langle \nabla \varphi ,z_2^{*}\rangle }_{L^2}\hfill \\ & & +G_2^{*}\left(z_1^{*}\right)+G_3^{*}\left(z_2^{*}\right)+G_4^{*}\left(z_3^{*}\right)\},\hfill \end{array}$$

(6)

where

$$C^{*}=\{v^{*}=(v_1^{*},v_2^{*})\in Y_1^{*}\times Y_3^{*}:-\mathrm{div}{\left(v_1^{*}\right)}_i+{\left(v_2^{*}\right)}_i=\mathbf{0},\forall i\in \{1,\cdots ,N\}\},$$

$$F_1^{*}(v_1^{*}+z_1^{*},\nabla \varphi )=\underset{v_1\in Y_1}{sup}\{{\langle v_1,z_1^{*}+v_1^{*}\rangle }_{L^2}-F_1(v_1,\nabla \varphi )\},$$

and

$$G_1^{*}(v_2^{*}+z_2^{*})=\underset{u\in V}{sup}\{{\langle u,v_2^{*}+z_2^{*}\rangle }_{L^2}-G_1\left(u\right)\}.$$

Thus, defining

$$J_2^{*}(\varphi ,z^{*},v^{*})=F_1^{*}(v_1^{*}+z_1^{*},\nabla \varphi )-G_1^{*}(v_2^{*}+z_3^{*})-{\langle \nabla \varphi ,z_2^{*}\rangle }_{L^2}+G_2^{*}\left(z_1^{*}\right)+G_3^{*}\left(z_2^{*}\right)+G_4^{*}\left(z_3^{*}\right),$$

we have got

$$\begin{array}{ccc}\hfill \underset{u\in V}{inf}J\left(u\right)& =& \underset{(u,\varphi )\in V\times V_0}{inf}{J}_1(u,\varphi )\hfill \\ & =& \underset{(u,\varphi ,z^{*})\in V\times V_0\times Y^{*}}{inf}{J}_1^{*}(u,\varphi ,z^{*})\hfill \\ & =& \underset{z^{*}\in Y^{*}}{inf}\left\{\underset{\varphi \in V_0}{inf}\left\{\underset{v^{*}\in C^{*}}{sup}{J}_2^{*}(\varphi ,z^{*},v^{*})\right\}\right\}.\hfill \end{array}$$

(7)

Finally, observe that

$$\begin{array}{ccc}& & \underset{u\in V}{inf}J\left(u\right)\hfill \\ & =& \underset{z^{*}\in Y^{*}}{inf}\left\{\underset{\varphi \in V_0}{inf}\left\{\underset{v^{*}\in C^{*}}{sup}{J}_2^{*}(\varphi ,z^{*},v^{*})\right\}\right\}\hfill \\ & \ge & \underset{v^{*}\in C^{*}}{sup}\left\{\underset{(z^{*},\varphi )\in Y^{*}\times V_0}{inf}{J}_2^{*}(\varphi ,z^{*},v^{*})\right\}.\hfill \end{array}$$

(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 $\mathsf{\Omega }=[0,1]\subset \mathbb{R}$ and consider a functional $J:V\to \mathbb{R}$ where

$$J\left(u\right)=\frac{1}{2}{\int }_{\mathsf{\Omega }}{({\left(u^{\prime }\right)}^2-1)}^{2}dx+\frac{1}{2}{\int }_{\mathsf{\Omega }}{u}^{2}dx-{\langle u,f\rangle }_{L^2},$$

and where

$$V=\{u\in W^{1,4}\left(\mathsf{\Omega }\right):u\left(0\right)=0\mathrm{and}u\left(1\right)=1/2\}$$

and $f\in L^2\left(\mathsf{\Omega }\right).$

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}\left(\mathsf{\Omega }\right)$, at this point we define, $F:V\to \mathbb{R}$ and $F_1:V\times V_0\to \mathbb{R}$ by

$$F\left(u\right)=\frac{1}{2}{\int }_{\mathsf{\Omega }}{({\left(u^{\prime }\right)}^2-1)}^{2}dx,$$

and

$$F_1(u,\varphi )=\frac{1}{2}{\int }_{\mathsf{\Omega }}{({(u^{\prime }+{\varphi }^{\prime })}^2-1)}^{2}dx.$$

Observe

$$F\left(u\right)\ge \underset{\varphi \in V_0}{inf}{F}_1(u,\varphi )\ge Q_F\left(u\right),\forall u\in V,$$

where $Q_F\left(u\right)$ refers to a quasi-convex regularization of $F.$

We define also

$$F_2:V\times V_0\to \mathbb{R},$$

$$F_3:V\times V_0\to \mathbb{R}$$

and

$$G:V\times V_0\to \mathbb{R}$$

by

$$F_2(u,\varphi )=\frac{1}{2}{\int }_{\mathsf{\Omega }}{({(u^{\prime }+{\varphi }^{\prime })}^2-1)}^{2}dx+\frac{1}{2}{\int }_{\mathsf{\Omega }}{u}^{2}dx-{\langle u,f\rangle }_{L^2},$$

$$\begin{array}{ccc}\hfill F_3(u,\varphi )& =& F_2(u,\varphi )+\frac{K}{2}{\int }_{\mathsf{\Omega }}{\left(u^{\prime }\right)}^{2}dx\hfill \\ & & +\frac{K_1}{2}{\int }_{\mathsf{\Omega }}{\left({\varphi }^{\prime }\right)}^{2}dx\hfill \end{array}$$

(9)

and

$$\begin{array}{ccc}\hfill G(u,\varphi )& =& \frac{K}{2}{\int }_{\mathsf{\Omega }}{\left(u^{\prime }\right)}^{2}dx\hfill \\ & & +\frac{K_1}{2}{\int }_{\mathsf{\Omega }}{\left({\varphi }^{\prime }\right)}^{2}dx\hfill \end{array}$$

(10)

Observe that if $K>0,K_1>0$ is large enough, both $F_3$ and G are convex.

Denoting $Y=Y^{*}=L^2\left(\mathsf{\Omega }\right)$ we also define the polar functional $G^{*}:Y^{*}\times Y^{*}\to \mathbb{R}$ by

$$G^{*}(v^{*},v_0^{*})=\underset{(u,\varphi )\in V\times V_0}{sup}\{{\langle u,v^{*}\rangle }_{L^2}+{\langle \varphi ,v_0^{*}\rangle }_{L^2}-G(u,\varphi )\}.$$

Observe that

$$\underset{u\in U}{inf}J\left(u\right)\ge \underset{((u,\varphi ),(v^{*},v_0^{*}))\in V\times V_0\times {\left[Y^{*}\right]}^2}{inf}\{G^{*}(v^{*},v_0^{*})-{\langle u,v^{*}\rangle }_{L^2}-{\langle \varphi ,v_0^{*}\rangle }_{L^2}+F_3(u,\varphi )\}.$$

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 {\left[Y^{*}\right]}^2\to \mathbb{R}$, where

$$J_1^{*}(u,\varphi ,v^{*},v_0^{*})=G^{*}(v^{*},v_0^{*})-{\langle u,v^{*}\rangle }_{L^2}-{\langle \varphi ,v_0^{*}\rangle }_{L^2}+F_3(u,\varphi ).$$

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<\epsilon \ll 1.$
• Choose $(u_1,{\varphi }_1)\in V\times V_0,$ such that $\parallel u_1{\parallel }_{1,\infty }\ll K/4$ and $\parallel {\varphi }_1{\parallel }_{1,\infty }\ll K/4.$
• Set $n=1.$
• Calculate $(v_n^{*},{\left(v_0^{*}\right)}_n)$ solution of the system of equations:

  $$\frac{\partial J_1^{*}(u_n,{\varphi }_n,v_n^{*},{\left(v_0^{*}\right)}_n)}{\partial v^{*}}=\mathbf{0}$$

  and

  $$\frac{\partial J_1^{*}(u_n,{\varphi }_n,v_n^{*},{\left(v_0^{*}\right)}_n)}{\partial v_0^{*}}=\mathbf{0},$$

  that is

  $$\frac{\partial G^{*}(v_n^{*},{\left(v_0^{*}\right)}_n)}{\partial v^{*}}-u_n=0$$

  and

  $$\frac{\partial G^{*}(v_n^{*},{\left(v_0^{*}\right)}_n)}{\partial v_0^{*}}-{\varphi }_n=0$$

  so that

  $$v_n^{*}=\frac{\partial G(u_n,{\varphi }_n)}{\partial u}$$

  and

  $${\left(v_0^{*}\right)}_n^{*}=\frac{\partial G(u_n,{\varphi }_n)}{\partial \varphi }$$
• Calculate $(u_{n+1},{\varphi }_{n+1})$ by solving the system of equations:

  $$\frac{\partial J_1^{*}(u_{n+1},{\varphi }_{n+1},v_n^{*},{\left(v_0^{*}\right)}_n)}{\partial u}=\mathbf{0}$$

  and

  $$\frac{\partial J_1^{*}(u_{n+1},{\varphi }_{n+1},v_n^{*},{\left(v_0^{*}\right)}_n)}{\partial \varphi }=\mathbf{0}$$

  that is

  $$-v_n^{*}+\frac{\partial F_3(u_{n+1},{\varphi }_{n+1})}{\partial u}=\mathbf{0}$$

  and

  $$-{\left(v_0^{*}\right)}_n+\frac{\partial F_3(u_{n+1},{\varphi }_{n+1})}{\partial \varphi }=\mathbf{0}$$
• If $max\{\parallel u_n-u_{n+1}{\parallel }_{\infty },\parallel {\varphi }_{n+1}-{\varphi }_n{\parallel }_{\infty }\}\le \epsilon$, then stop, else set $n:=n+1$ and go to item 4.

For the case in which $f\left(x\right)=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\left(x\right)=sin\left(\pi x\right)/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.

Let $\mathsf{\Omega }=[0,1]\subset \mathbb{R}$ and consider a functional $J:V\to \mathbb{R}$ where

$$J\left(u\right)=\frac{1}{2}{\int }_{\mathsf{\Omega }}min\{{(u^{\prime }-1)}^2,{(u^{\prime }+1)}^2\}dx+\frac{1}{2}{\int }_{\mathsf{\Omega }}{u}^{2}dx-{\langle u,f\rangle }_{L^2},$$

and where

$$V=\{u\in W^{1,2}\left(\mathsf{\Omega }\right):u\left(0\right)=0\mathrm{and}u\left(1\right)=1/2\}$$

and $f\in L^2\left(\mathsf{\Omega }\right).$

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 \mathbb{R}$ and $G:V\to \mathbb{R}$ by

$$\begin{array}{ccc}\hfill F\left(u\right)& =& \frac{1}{2}{\int }_{\mathsf{\Omega }}min\{{(u^{\prime }-1)}^2,{(u^{\prime }+1)}^2\}dx\hfill \\ & =& \frac{1}{2}{\int }_{\mathsf{\Omega }}{\left(u^{\prime }\right)}^{2}dx-{\int }_{\mathsf{\Omega }}\left|u^{\prime }\right|dx+1/2\hfill \\ & \equiv & F_1\left(u^{\prime }\right),\hfill \end{array}$$

(11)

and

$$G\left(u\right)=\frac{1}{2}{\int }_{\mathsf{\Omega }}{u}^{2}dx-{\langle u,f\rangle }_{L^2}.$$

Denoting $Y=Y^{*}=L^2\left(\mathsf{\Omega }\right)$ we also define the polar functional $F_1^{*}:Y^{*}\to \mathbb{R}$ and $G^{*}:Y^{*}\to \mathbb{R}$ by

$$\begin{array}{ccc}\hfill F_1^{*}\left(v^{*}\right)& =& \underset{v\in Y}{sup}\{{\langle v,v^{*}\rangle }_{L^2}-F_1\left(v\right)\}\hfill \\ & =& \frac{1}{2}{\int }_{\mathsf{\Omega }}{\left(v^{*}\right)}^{2}dx+{\int }_{\mathsf{\Omega }}\left|v^{*}\right|dx,\hfill \end{array}$$

(12)

and

$$\begin{array}{ccc}\hfill G^{*}\left({\left(v^{*}\right)}^{\prime }\right)& =& \underset{u\in V}{sup}\{-{\langle u^{\prime },v^{*}\rangle }_{L^2}-G\left(u\right)\}\hfill \\ & =& \frac{1}{2}{\int }_{\mathsf{\Omega }}{({\left(v^{*}\right)}^{\prime }+f)}^{2}dx-\frac{1}{2}{v}^{*}\left(1\right).\hfill \end{array}$$

(13)

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

$$\underset{u\in V}{inf}J\left(u\right)=\underset{v^{*}\in Y^{*}}{max}\{-F_1^{*}\left(v^{*}\right)-G^{*}(-{\left(v^{*}\right)}^{\prime })\}.$$

Indeed, from the direct method of the calculus of variations, the maximum for the dual formulation is attained at some ${\widehat{v}}^{*}\in Y^{*}$.

Moreover, the corresponding solution $u_0\in V$ is obtained from the equation

$$u_0=\frac{\partial G\left({\left({\widehat{v}}^{*}\right)}^{\prime }\right)}{\partial {\left(v^{*}\right)}^{\prime }}={\left({\widehat{v}}^{*}\right)}^{\prime }+f.$$

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

$$\left\{\begin{array}{cc}{\left(v^{*}\right)}^{″}+f^{\prime }-v^{*}-\mathrm{sign}\left(v^{*}\right)=0,\hfill & \mathrm{in}\mathsf{\Omega },\hfill \\ {\left(v^{*}\right)}^{\prime }\left(0\right)=0,{\left(v^{*}\right)}^{\prime }\left(1\right)=1/2,\hfill \end{array}\right.$$

(14)

where $\mathrm{sign}\left(v^{*}\left(x\right)\right)=1$ if $v^{*}\left(x\right)>0,$$\mathrm{sign}\left(v^{*}\left(x\right)\right)=-1,$ if $v^{*}\left(x\right)<0$ and

$$-1\le \mathrm{sign}\left(v^{*}\left(x\right)\right)\le 1,$$

if $v^{*}\left(x\right)=0.$

We have computed the solutions $v^{*}$ and corresponding solutions $u_0\in V$ for the cases in which $f\left(x\right)=0$ and $f\left(x\right)=sin\left(\pi x\right)/2.$

For the solution $u_0\left(x\right)$ for the case in which $f\left(x\right)=0$, please see Figure 3.

For the solution $u_0\left(x\right)$ for the case in which $f\left(x\right)=sin\left(\pi x\right)/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

$$\left\{\begin{array}{cc}{\left(v^{*}\right)}^{\prime \prime }+f^{\prime }-v^{*}-\mathrm{sign}\left(v^{*}\right)=0,\hfill & \mathrm{in}\mathsf{\Omega },\hfill \\ {\left(v^{*}\right)}^{\prime }\left(0\right)=0,{\left(v^{*}\right)}^{\prime }\left(1\right)=1/2.\hfill \end{array}\right.$$

(15)

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

$$|v^{*}|\approx \sqrt{{\left(v^{*}\right)}^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;$
• $fori=1:m_8$
  $yo(i,1)=0.01;$
  $y_1(i,1)=sin(\pi \ast i/m_8)/2;$
  $end;$
• $fori=1:m_8-1$
  $d{y}_1(i,1)=(y_1(i+1,1)-y_1(i,1))/d;$
  $end;$
• $fork=1:3000$ (we have fixed the number of iterations)
  $i=1;$
  $h_3=1/\sqrt{vo{(i,1)}^2+e_1};$
  $m_{12}=1+d^2\ast h_3+d^2;$
  $m_{50}\left(i\right)=1/m_{12};$
  $z\left(i\right)=m_{50}\left(i\right)\ast (d{y}_1(i,1)\ast d^2);$
• $fori=2:m_8-1$
  $h_3=1/\sqrt{vo{(i,1)}^2+e_1};$
  $m_{12}=2+h_3\ast d^2+d^2-m50(i-1);$
  $m50\left(i\right)=1/m_{12};$
  $z\left(i\right)=m_{50}\left(i\right)\ast (z(i-1)+d{y}_1(i,1)\ast d^2);$
  $end;$
• $v(m_8,1)=(d/2+z(m_8-1))/(1-m_{50}(m_8-1));$
• $fori=1:m_8-1$
  $v(m_8-i,1)=m_{50}(m_8-i)\ast v(m_8-i+1)+z(m_8-i);$
  $end;$
• $v(m_8/2,1)$
• $vo=v;$
  $end;$
• $fori=1:m_8-1$
  $u(i,1)=(v(i+1,1)-v(i,1))/d+y_1(i,1);$
  $end;$
• $fori=1:m_8-1$
  $x\left(i\right)=i*d;$
  $end;$
  $plot(x,u(:,1\left)\right)$

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

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 $\mathsf{\Omega }\subset {\mathbb{R}}^2$ be an open, bounded, connected set which represents the middle surface of a plate of thickness h. The boundary of $\mathsf{\Omega }$, which is assumed to be regular (Lipschitzian), is denoted by $\partial \mathsf{\Omega }$. The vectorial basis related to the cartesian system $\{x_1,x_2,x_3\}$ is denoted by $({\mathbf{a}}_{\alpha },{\mathbf{a}}_3)$, where $\alpha =1,2$ (in general Greek indices stand for 1 or 2), and where ${\mathbf{a}}_3$ is the vector normal to $\mathsf{\Omega }$, whereas ${\mathbf{a}}_1$ and ${\mathbf{a}}_2$ are orthogonal vectors parallel to $\mathsf{\Omega }.$ Also, $\mathbf{n}$ is the outward normal to the plate surface.

The displacements will be denoted by

$$\widehat{\mathbf{u}}=\{{\widehat{u}}_{\alpha },{\widehat{u}}_3\}={\widehat{u}}_{\alpha }{\mathbf{a}}_{\alpha }+{\widehat{u}}_3{\mathbf{a}}_3.$$

The Kirchhoff-Love relations are

$$\begin{array}{ccc}& & {\widehat{u}}_{\alpha }(x_1,x_2,x_3)=u_{\alpha }(x_1,x_2)-x_{3}w{(x_1,x_2)}_{,\alpha }\hfill \\ & & \mathrm{and}{\widehat{u}}_3(x_1,x_2,x_3)=w(x_1,x_2).\hfill \end{array}$$

(16)

Here $-h/2\le x_3\le h/2$ so that we have $u=(u_{\alpha },w)\in U$ where

$$\begin{array}{ccc}\hfill U& =& \left\{u=(u_{\alpha },w)\in W^{1,2}(\mathsf{\Omega };{\mathbb{R}}^2)\times W^{2,2}\left(\mathsf{\Omega }\right),\right.\hfill \\ & & u_{\alpha }=w=\frac{\partial w}{\partial \mathbf{n}}=0\left.\mathrm{on}\partial \mathsf{\Omega }\right\}\hfill \\ & =& W_0^{1,2}(\mathsf{\Omega };{\mathbb{R}}^2)\times W_0^{2,2}\left(\mathsf{\Omega }\right).\hfill \end{array}$$

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

We also define the operator $\mathsf{\Lambda }:U\to Y\times Y$, where $Y=Y^{*}=L^2(\mathsf{\Omega };{\mathbb{R}}^{2\times 2})$, by

$$\mathsf{\Lambda }\left(u\right)=\left\{\gamma \right(u),\kappa (u\left)\right\},$$

$${\gamma }_{\alpha \beta }\left(u\right)=\frac{u_{\alpha ,\beta }+u_{\beta ,\alpha }}{2}+\frac{w_{,\alpha }{w}_{,\beta }}{2},$$

$${\kappa }_{\alpha \beta }\left(u\right)=-w_{,\alpha \beta }.$$

The constitutive relations are given by

$$N_{\alpha \beta }\left(u\right)=H_{\alpha \beta \lambda \mu }{\gamma }_{\lambda \mu }\left(u\right),$$

(17)

$$M_{\alpha \beta }\left(u\right)=h_{\alpha \beta \lambda \mu }{\kappa }_{\lambda \mu }\left(u\right),$$

(18)

where: $\left\{H_{\alpha \beta \lambda \mu }\right\}$ and $\left\{h_{\alpha \beta \lambda \mu }=\frac{h^2}{12}{H}_{\alpha \beta \lambda \mu }\right\}$, are symmetric positive definite fourth order tensors. From now on, we denote $\left\{{\overline{H}}_{\alpha \beta \lambda \mu }\right\}={\left\{H_{\alpha \beta \lambda \mu }\right\}}^{-1}$ and $\left\{{\overline{h}}_{\alpha \beta \lambda \mu }\right\}={\left\{h_{\alpha \beta \lambda \mu }\right\}}^{-1}$.

Furthermore $\left\{N_{\alpha \beta }\right\}$ denote the membrane force tensor and $\left\{M_{\alpha \beta }\right\}$ the moment one. The plate stored energy, represented by $(G\circ \mathsf{\Lambda }):U\to \mathbb{R}$ is expressed by

$$(G\circ \mathsf{\Lambda })\left(u\right)=\frac{1}{2}{\int }_{\mathsf{\Omega }}{N}_{\alpha \beta }\left(u\right){\gamma }_{\alpha \beta }\left(u\right)dx+\frac{1}{2}{\int }_{\mathsf{\Omega }}{M}_{\alpha \beta }\left(u\right){\kappa }_{\alpha \beta }\left(u\right)dx$$

(19)

and the external work, represented by $F:U\to \mathbb{R}$, is given by

$$F\left(u\right)={\langle w,P\rangle }_{L^2}+{\langle u_{\alpha },P_{\alpha }\rangle }_{L^2},$$

(20)

where $P,P_1,P_2\in L^2\left(\mathsf{\Omega }\right)$ are external loads in the directions ${\mathbf{a}}_3$, ${\mathbf{a}}_1$ and ${\mathbf{a}}_2$ respectively. The potential energy, denoted by $J:U\to \mathbb{R}$ is expressed by:

$$J\left(u\right)=(G\circ \mathsf{\Lambda })\left(u\right)-F\left(u\right)$$

Define now $J_3:\tilde{U}\to \mathbb{R}$ by

$$J_3\left(u\right)=J\left(u\right)+J_5\left(w\right).$$

where

$$J_5\left(w\right)=10{\int }_{\mathsf{\Omega }}\frac{a^{Kbw}}{ln\left(a\right)K^{3/2}}dx+10{\int }_{\mathsf{\Omega }}\frac{a^{-K(bw-1/100)}}{ln\left(a\right)K^{3/2}}dx.$$

In such a case for $a=2.71$, $K=185$, $b=P/\left|P\right|$ in $\mathsf{\Omega }$ and

$$\tilde{U}={\{u\in U:\parallel w\parallel }_{\infty }\le 0.01\mathrm{and}Pw\ge 0\mathrm{a}.\mathrm{e}.\mathrm{in}\mathsf{\Omega }\},$$

we get

$$\begin{array}{ccc}\hfill \frac{\partial J_3\left(u\right)}{\partial w}& =& \frac{\partial J\left(u\right)}{\partial w}+\frac{\partial J_5\left(u\right)}{\partial w}\hfill \\ & \approx & \frac{\partial J\left(u\right)}{\partial w}+\mathcal{O}(\pm 3.0),\hfill \end{array}$$

(21)

and

$$\begin{array}{ccc}\hfill \frac{{\partial }^2{J}_3\left(u\right)}{\partial w^2}& =& \frac{{\partial }^{2}J\left(u\right)}{\partial w^2}+\frac{{\partial }^2{J}_5\left(u\right)}{\partial w^2}\hfill \\ & \approx & \frac{{\partial }^{2}J\left(u\right)}{\partial w^2}+\mathcal{O}\left(850\right).\hfill \end{array}$$

(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\left(u\right)}{\partial w^2},$ which has increased of order $\mathcal{O}(700-1000).$

Moreover the difference between the approximate and exact equation

$$\frac{\partial J\left(u\right)}{\partial w}=\mathbf{0}$$

is of order $\mathcal{O}(\pm 3.0)$ which corresponds to a small perturbation in the original equation for a load of $P=1500N/m^2,$ for example. Summarizing, the exact equation may be approximately 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 $\mathsf{\Omega }\subset {\mathbb{R}}^3$ be an open, bounded and connected set with a regular (Lipschitzian) boundary denoted by $\partial \mathsf{\Omega }.$

Consider a functional $J:V\to \mathbb{R}$ where

$$\begin{array}{ccc}\hfill J\left(u\right)& =& \frac{\gamma }{2}{\int }_{\mathsf{\Omega }}\nabla u·\nabla udx+\frac{\alpha }{2}{\int }_{\mathsf{\Omega }}{(u^2-\beta )}^{2}dx\hfill \\ & & -{\langle u,f\rangle }_{L^2},\hfill \end{array}$$

(23)

where $\gamma >0$, $\alpha >0$, $\beta >0$, $V=W_0^{1,2}\left(\mathsf{\Omega }\right)$ and $f\in L^2\left(\mathsf{\Omega }\right).$

We define

$$A^{+}=\{u\in V:uf\ge 0,\mathrm{a}.\mathrm{e}.\mathrm{in}\mathsf{\Omega }\},$$

$$V_2={\{u\in V:\parallel u\parallel }_{\infty }\le 1\},$$

and

$$V_1=V_2\cap A^{+}.$$

At this point we define $v=u/10$ so that

$$\begin{array}{ccc}\hfill J\left(u\right)& =& J_1\left(v\right)\hfill \\ & =& \frac{{10}^2\gamma }{2}{\int }_{\mathsf{\Omega }}\nabla v·\nabla vdx+\frac{\alpha }{2}{\int }_{\mathsf{\Omega }}({\left(10v\right)}^2-\beta )dx\hfill \\ & & -{\langle 10v,f\rangle }_{L^2}.\hfill \end{array}$$

(24)

Moreover we define

$$\begin{array}{ccc}\hfill J_2\left(v\right)& =& \frac{1}{10}{J}_1\left(v\right)\hfill \\ & =& \frac{10\gamma }{2}{\int }_{\mathsf{\Omega }}\nabla v·\nabla vdx+\frac{\alpha }{20}{\int }_{\mathsf{\Omega }}({\left(10v\right)}^2-\beta )dx\hfill \\ & & -{\langle v,f\rangle }_{L^2},\hfill \end{array}$$

(25)

and $J_3:U_3\to \mathbb{R}$ where

$$J_3\left(v\right)=J_2\left(v\right)+J_5\left(v\right)$$

where

$$J_5\left(v\right)=K_1\left({\int }_{\mathsf{\Omega }}\frac{a^{K^3\left(5bw\right)}}{ln\left(a\right)K^4}dx+{\int }_{\mathsf{\Omega }}\frac{a^{-K^3(5bw-0.5)}}{ln\left(a\right)K^4}dx\right).$$

Here $K_1=1/360$, $a=2.71$, $K=2$, $b=f/\left|f\right|$ in $\mathsf{\Omega }$ and

$$U_2=\{v\in V:fv\ge 0,\mathrm{a}.\mathrm{e}.\mathrm{in}\mathsf{\Omega }\},$$

$$U_4={\{v\in V:\parallel v\parallel }_{\infty }\le 1/10\},$$

and

$$U_3=U_2\cap U_4.$$

Thus, with such numerical values, we may obtain

$$\begin{array}{ccc}\hfill \frac{\partial J_3\left(v\right)}{\partial v}& =& \frac{\partial J_2\left(v\right)}{\partial v}+\frac{\partial J_5\left(v\right)}{\partial v}\hfill \\ & \approx & \frac{\partial J_2\left(v\right)}{\partial v}+\mathcal{O}(\pm 0.3),\hfill \end{array}$$

(26)

and

$$\begin{array}{ccc}\hfill \frac{{\partial }^2{J}_3\left(v\right)}{\partial v^2}& =& \frac{{\partial }^2{J}_2\left(v\right)}{\partial v^2}+\frac{{\partial }^2{J}_5\left(v\right)}{\partial v^2}\hfill \\ & \approx & \frac{{\partial }^2{J}_2\left(v\right)}{\partial v^2}+\mathcal{O}\left(7.0\right).\hfill \end{array}$$

(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\left(v\right)}{\partial v^2},$ which has increased of order $\mathcal{O}(5-14).$

Moreover the difference between the approximate and exact equation

$$\frac{\partial J_2\left(v\right)}{\partial v}=\mathbf{0}$$

is of order $\mathcal{O}(\pm 0.3)$ which for appropriate parameters $\gamma >0,\alpha >0$ and $\beta >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 $\gamma >0,$$\alpha >0$ and $\beta >0$ are such that

$$\frac{{\partial }^2{J}_3\left(v\right)}{\partial v^2}>\mathbf{0},$$

in $U_3$

Assume also, $v_0\in U_3$ is such that

$$\delta J_3\left(v_0\right)=\mathbf{0}.$$

Under such hypotheses, $J_3$ is convex on $U_3$ so that

$$J_3\left(v_0\right)=\underset{v\in U_3}{min}{J}_3\left(v\right).$$

Moreover,

$$\delta J\left(u_0\right)=\mathbf{0}+\mathcal{O}(\pm 0.3),$$

where $u_0=10v_0\in V_1$

Proof. 

From the hypotheses

$$\frac{{\partial }^2{J}_3\left(v\right)}{\partial v^2}>\mathbf{0}$$

in $U_3$, so that $J_3$ is convex on the convex set $U_3.$

Consequently, since $\delta J_3\left(v_0\right)=\mathbf{0}$, we obtain

$$J_3\left(v_0\right)=\underset{v\in U_3}{min}{J}_3\left(v\right).$$

Finally, from the approximation indicated in the last remark and $u_0\in V_1$ we get

$$\delta J\left(u_0\right)=\mathbf{0}+\mathcal{O}(\pm 0.3).$$

The proof is complete.

□

7. An Exact Convex Dual Variational Formulation for a Non-Convex Primal One

In this section we develop a convex dual variational formulation suitable to compute a critical point for the corresponding primal one.

Let $\mathsf{\Omega }\subset {\mathbb{R}}^2$ be an open, bounded, connected set with a regular (Lipschitzian) boundary denoted by $\partial \mathsf{\Omega }.$

Consider a functional $J:V\to \mathbb{R}$ where

$$J\left(u\right)=F(u_x,u_y)-{\langle u,f\rangle }_{L^2},$$

$V=W_0^{1,2}\left(\mathsf{\Omega }\right)$ and $f\in L^2\left(\mathsf{\Omega }\right).$

Here we denote $Y=Y^{*}=L^2\left(\mathsf{\Omega }\right)$ and $Y_1=Y_1^{*}=L^2\left(\mathsf{\Omega }\right)\times L^2\left(\mathsf{\Omega }\right).$

Defining

$$V_1={\{u\in V:\parallel u\parallel }_{1,\infty }\le K_1\}$$

for some appropriate $K_1>0$, suppose also F is twice Fréchet differentiable and

$$det\left\{\frac{{\partial }^{2}F(u_x,u_y)}{\partial v_1\partial v_2}\right\}\ne 0,$$

$\forall u\in V_1.$

Define now $F_1:V\to \mathbb{R}$ and $F_2:V\to \mathbb{R}$ by

$$F_1(u_x,u_y)=F(u_x,u_y)+\frac{\epsilon }{2}{\int }_{\mathsf{\Omega }}{u}_x^{2}dx+\frac{\epsilon }{2}{\int }_{\mathsf{\Omega }}{u}_y^{2}dx,$$

and

$$F_2(u_x,u_y)=\frac{\epsilon }{2}{\int }_{\mathsf{\Omega }}{u}_x^{2}dx+\frac{\epsilon }{2}{\int }_{\mathsf{\Omega }}{u}_y^{2}dx,$$

where here we denote $dx=d{x}_{1}d{x}_2.$

Moreover, we define the respective Legendre transform functionals $F_1^{*}$ and $F_2^{*}$ as

$$F_1^{*}\left(v^{*}\right)={\langle v_1,v_1^{*}\rangle }_{L^2}+{\langle v_2,v_2^{*}\rangle }_{L^2}-F_1(v_1,v_2),$$

where $v_1,v_2\in Y$ are such that

$$v_1^{*}=\frac{\partial F_1(v_1,v_2)}{\partial v_1},$$

$$v_2^{*}=\frac{\partial F_1(v_1,v_2)}{\partial v_2},$$

and

$$F_2^{*}\left(v^{*}\right)={\langle v_1,v_1^{*}+f_1\rangle }_{L^2}+{\langle v_2,v_2^{*}\rangle }_{L^2}-F_2(v_1,v_2),$$

where $v_1,v_2\in Y$ are such that

$$v_1^{*}+f_1=\frac{\partial F_2(v_1,v_2)}{\partial v_1},$$

$$v_2^{*}=\frac{\partial F_2(v_1,v_2)}{\partial v_2}.$$

Here $f_1$ is any function such that

$${\left(f_1\right)}_x=f,\mathrm{in}\mathsf{\Omega }.$$

Furthermore, we define

$$\begin{array}{ccc}\hfill J^{*}\left(v^{*}\right)& =& -F_1^{*}\left(v^{*}\right)+F_2^{*}\left(v^{*}\right)\hfill \\ & =& -F_1^{*}\left(v^{*}\right)+\frac{1}{2\epsilon }{\int }_{\mathsf{\Omega }}{(v_1^{*}+f_1)}^{2}dx+\frac{1}{2\epsilon }{\int }_{\mathsf{\Omega }}{\left(v_2^{*}\right)}^{2}dx.\hfill \end{array}$$

(28)

Observe that through the target conditions

$$v_1^{*}+f_1=\epsilon u_x,$$

$$v_2^{*}=\epsilon u_y,$$

we may obtain the compatibility condition

$${(v_1^{*}+f_1)}_y-{\left(v_2^{*}\right)}_x=0.$$

Define now

$$A^{*}=\{v^{*}=(v_1^{*},v_2^{*})\in B_r(0,0)\subset Y_1^{*}:{(v_1^{*}+f_1)}_y-{\left(v_2^{*}\right)}_x=0,\mathrm{in}\mathsf{\Omega }\},$$

for some appropriate $r>0$ such that $J^{*}$ is convex in $B_r(0,0).$

Consider the problem of minimizing $J^{*}$ subject to $v^{*}\in A^{*}.$

Assuming $r>0$ is large enough so that the restriction in r is not active, at this point we define the associated Lagrangian

$$J_1^{*}(v^{*},\phi )=J^{*}\left(v^{*}\right)+{\langle \phi ,{(v_1^{*}+f)}_y-{\left(v_2^{*}\right)}_x\rangle }_{L^2},$$

where $\phi$ is an appropriate Lagrange multiplier.

Therefore

$$\begin{array}{ccc}\hfill J_1^{*}\left(v^{*}\right)& =& -F_1^{*}\left(v^{*}\right)+\frac{1}{2\epsilon }{\int }_{\mathsf{\Omega }}{(v_1^{*}+f_1)}^{2}dx+\frac{1}{2\epsilon }{\int }_{\mathsf{\Omega }}{\left(v_2^{*}\right)}^{2}dx\hfill \\ & & +{\langle \phi ,{(v_1^{*}+f)}_y-{\left(v_2^{*}\right)}_x\rangle }_{L^2}.\hfill \end{array}$$

(29)

The optimal point in question will be a solution of the corresponding Euler-Lagrange equations for $J_1^{*}.$

From the variation of $J_1^{*}$ in $v_1^{*}$ we obtain

$$-\frac{\partial F_1^{*}\left(v^{*}\right)}{\partial v_1^{*}}+\frac{v_1^{*}+f}{\epsilon }-\frac{\partial \phi }{\partial y}=0.$$

(30)

From the variation of $J_1^{*}$ in $v_2^{*}$ we obtain

$$-\frac{\partial F_1^{*}\left(v^{*}\right)}{\partial v_2^{*}}+\frac{v_2^{*}}{\epsilon }+\frac{\partial \phi }{\partial x}=0.$$

(31)

From the variation of $J_1^{*}$ in $\phi$ we have

$${(v_1^{*}+f)}_y-{\left(v_2^{*}\right)}_x=0.$$

From this last equation, we may obtain $u\in V$ such that

$$v_1^{*}+f=\epsilon u_x,$$

and

$$v_2^{*}=\epsilon u_y.$$

From this and the previous extremal equations indicated we have

$$-\frac{\partial F_1^{*}\left(v^{*}\right)}{\partial v_1^{*}}+u_x-\frac{\partial \phi }{\partial y}=0,$$

and

$$-\frac{\partial F_1^{*}\left(v^{*}\right)}{\partial v_2^{*}}+u_y+\frac{\partial \phi }{\partial x}=0.$$

so that

$$v_1^{*}+f=\frac{\partial F_1(u_x-{\phi }_y,u_y+{\phi }_x)}{\partial v_1},$$

and

$$v_2^{*}=\frac{\partial F_1(u_x-{\phi }_y,u_y+{\phi }_x)}{\partial v_2}.$$

From this and Equation (30) and (31) we have

$$\begin{array}{ccc}& & -\epsilon {\left(\frac{\partial F_1^{*}\left(v^{*}\right)}{\partial v_1^{*}}\right)}_x-\epsilon {\left(\frac{\partial F_1^{*}\left(v^{*}\right)}{\partial v_2^{*}}\right)}_y\hfill \\ & & +{(v_1^{*}+f_1)}_x+{\left(v_2^{*}\right)}_y\hfill \\ & =& -\epsilon u_{xx}-\epsilon u_{yy}+{\left(v_1^{*}\right)}_x+{\left(v_2^{*}\right)}_y+f=0.\hfill \end{array}$$

(32)

Replacing the expressions of $v_1^{*}$ and $v_2^{*}$ into this last equation, we have

$$-\epsilon u_{xx}-\epsilon u_{yy}+{\left(\frac{\partial F_1(u_x-{\phi }_y,u_y+{\phi }_x)}{\partial v_1}\right)}_x+{\left(\frac{\partial F_1(u_x-{\phi }_y,u_y+{\phi }_x)}{\partial v_2}\right)}_y+f=0,$$

so that

$${\left(\frac{\partial F(u_x-{\phi }_y,u_y+{\phi }_x)}{\partial v_1}\right)}_x+{\left(\frac{\partial F(u_x-{\phi }_y,u_y+{\phi }_x)}{\partial v_2}\right)}_y+f=0,\mathrm{in}\mathsf{\Omega }.$$

(33)

Observe that if

$${\nabla }^2\phi =0$$

then there exists $\widehat{u}$ such that u and $\phi$ are also such that

$$u_x-{\phi }_y={\widehat{u}}_x$$

and

$$u_y+{\phi }_x={\widehat{u}}_y.$$

The boundary conditions for $\phi$ must be such that $\widehat{u}\in W_0^{1,2}.$

From this and Equation (33) we obtain

$$\delta J\left(\widehat{u}\right)=\mathbf{0}.$$

Summarizing, we may obtain a solution $\widehat{u}\in W_0^{1,2}$ of equation $\delta J\left(\widehat{u}\right)=\mathbf{0}$ by minimizing $J^{*}$ on $A^{*}$.

Finally, observe that clearly $J^{*}$ is convex in an appropriate large ball $B_r(0,0)$ for some appropriate $r>0$

8. Another Primal Dual Formulation for a Related Model

Let $\mathsf{\Omega }\subset {\mathbb{R}}^3$ be an open, bounded and connected set with a regular boundary denoted by $\partial \mathsf{\Omega }.$

Consider the functional $J:V\to \mathbb{R}$ where

$$\begin{array}{ccc}\hfill J\left(u\right)& =& \frac{\gamma }{2}{\int }_{\mathsf{\Omega }}\nabla u·\nabla udx+\frac{\alpha }{2}{\int }_{\mathsf{\Omega }}{(u^2-\beta )}^{2}dx\hfill \\ & & -{\langle u,f\rangle }_{L^2},\hfill \end{array}$$

(34)

$\alpha >0,\beta >0,\gamma >0$, $V=W_0^{1,2}\left(\mathsf{\Omega }\right)$ and $f\in L^2\left(\mathsf{\Omega }\right).$

Denoting $Y=Y^{*}=L^2\left(\mathsf{\Omega }\right)$, define now $J_1^{*}:V\times Y^{*}\to \mathbb{R}$ by

$$\begin{array}{ccc}\hfill J_1^{*}(u,v_0^{*})& =& -\frac{\gamma }{2}{\int }_{\mathsf{\Omega }}\nabla u·\nabla udx-{\langle u^2,v_0^{*}\rangle }_{L^2}\hfill \\ & & +\frac{K_1}{2}{\int }_{\mathsf{\Omega }}{(-\gamma {\nabla }^{2}u+2v_0^{*}u-f)}^{2}dx+{\langle u,f\rangle }_{L^2}\hfill \\ & & +\frac{1}{2\alpha }{\int }_{\mathsf{\Omega }}{\left(v_0^{*}\right)}^{2}dx+\beta {\int }_{\mathsf{\Omega }}{v}_0^{*}dx,\hfill \end{array}$$

(35)

Define also

$$A^{+}=\{u\in V:uf\ge 0,\mathrm{a}.\mathrm{e}.\mathrm{in}\mathsf{\Omega }\},$$

$$V_2={\{u\in V:\parallel u\parallel }_{\infty }\le K_3\},$$

and

$$V_1=V_2\cap A^{+}$$

for some appropriate $K_3>0$ to be specified.

Moreover define

$$B^{*}=\{v_0^{*}\in Y^{*}:\parallel v_0^{*}{\parallel }_{\infty }\le K\}$$

for some appropriate $K>0$ to be specified.

Observe that, denoting

$$\phi =-\gamma {\nabla }^{2}u+2v_0^{*}u-f$$

we have

$$\frac{{\partial }^2{J}_1^{*}(u,v_0^{*})}{\partial {\left(v_0^{*}\right)}^2}=\frac{1}{\alpha }+4K_1{u}^2$$

$$\frac{{\partial }^2{J}_1^{*}(u,v_0^{*})}{\partial u^2}=\gamma {\nabla }^2-2v_0^{*}+K_1{(-\gamma {\nabla }^2+2v_0^{*})}^2$$

and

$$\frac{{\partial }^2{J}_1^{*}(u,v_0^{*})}{\partial u\partial v_0^{*}}=K_1(2\phi +2(-\gamma {\nabla }^{2}u+2v_0^{*}u))-2u$$

so that

$$\begin{array}{ccc}& & det\left\{{\delta }^2{J}_1^{*}(u,v_0^{*})\right\}\hfill \\ & =& \frac{{\partial }^2{J}_1^{*}(u,v_0^{*})}{\partial {\left(v_0^{*}\right)}^2}\frac{{\partial }^2{J}_1^{*}(u,v_0^{*})}{\partial u^2}-{\left(\frac{{\partial }^2{J}_1^{*}(u,v_0^{*})}{\partial u\partial v_0^{*}}\right)}^2\hfill \\ & =& \frac{K_1{(-\gamma {\nabla }^2+2v_0^{*})}^2}{\alpha }-\frac{\gamma {\nabla }^2+2v_0^{*}+4\alpha u^2}{\alpha }\hfill \\ & & -4K_1^2{\phi }^2-8K_1\phi (-\gamma {\nabla }^2+2v_0^{*})u+8K_1\phi u\hfill \\ & & +4K_1(-\gamma {\nabla }^{2}u+2v_0^{*}u)u.\hfill \end{array}$$

(36)

Observe now that a critical point $\phi =0$ and $(-\gamma {\nabla }^{2}u+2v_0^{*}u)u=fu\ge 0$ in $\mathsf{\Omega }$.

Therefore, for an appropriate large $K_1>0$, also at a critical point, we have

$$\begin{array}{ccc}& & det\left\{{\delta }^2{J}_1^{*}(u,v_0^{*})\right\}\hfill \\ & =& 4K_{1}fu-\frac{{\delta }^{2}J\left(u\right)}{\alpha }+K_1\frac{{(-\gamma {\nabla }^2+2v_0^{*})}^2}{\alpha }>\mathbf{0}.\hfill \end{array}$$

(37)

Remark 8.1.

From this last equation we may observe that $J_1^{*}$ has a large region of convexity about any critical point $(u_0,{\widehat{v}}_0^{*})$, that is, there exists a large $r>0$ such that $J_1^{*}$ is convex on $B_r(u_0,{\widehat{v}}_0^{*}).$

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

Theorem 8.2.

Assume $K_1\gg max\{1,K,K_3\}$ and suppose $(u_0,{\widehat{v}}_0^{*})\in V_1\times B^{*}$ is such that

$$\delta J_1^{*}(u_0,{\widehat{v}}_0^{*})=\mathbf{0}.$$

Under such hypotheses, there exists $r>0$ such that $J_1^{*}$ is convex in $E^{*}=B_r(u_0,{\widehat{v}}_0^{*})\cap (V_1\times B^{*})$,

$$\delta J\left(u_0\right)=\mathbf{0},$$

and

$$-J\left(u_0\right)=J_1(u_0,{\widehat{v}}_0^{*})=\underset{(u,v_0^{*})\in E^{*}}{inf}{J}_1^{*}(u,v_0^{*}).$$

9. A Third Primal Dual Formulation for a Related Model

Let $\mathsf{\Omega }\subset {\mathbb{R}}^3$ be an open, bounded and connected set with a regular boundary denoted by $\partial \mathsf{\Omega }.$

Consider the functional $J:V\to \mathbb{R}$ where

$$\begin{array}{ccc}\hfill J\left(u\right)& =& \frac{\gamma }{2}{\int }_{\mathsf{\Omega }}\nabla u·\nabla udx+\frac{\alpha }{2}{\int }_{\mathsf{\Omega }}{(u^2-\beta )}^{2}dx\hfill \\ & & -{\langle u,f\rangle }_{L^2},\hfill \end{array}$$

(38)

$\alpha >0,\beta >0,\gamma >0$, $V=W_0^{1,2}\left(\mathsf{\Omega }\right)$ and $f\in L^2\left(\mathsf{\Omega }\right).$

Denoting $Y=Y^{*}=L^2\left(\mathsf{\Omega }\right)$, define now $J_1^{*}:V\times Y^{*}\times Y^{*}\to \mathbb{R}$ by

$$\begin{array}{ccc}\hfill J_1^{*}(u,v_1^{*},v_0^{*})& =& \frac{\gamma }{2}{\int }_{\mathsf{\Omega }}\nabla u·\nabla udx+\frac{1}{2}{\int }_{\mathsf{\Omega }}K{u}^{2}dx\hfill \\ & & -{\langle u,v_1^{*}\rangle }_{L^2}+\frac{1}{2}{\int }_{\mathsf{\Omega }}\frac{{\left(v_1^{*}\right)}^2}{(-2v_0^{*}+K)}dx\hfill \\ & & +\frac{K_1}{2}{\int }_{\mathsf{\Omega }}{(v_1^{*}-(-\gamma {\nabla }^{2}u)-Ku)}^{2}dx\hfill \\ & & \frac{1}{2}{\int }_{\mathsf{\Omega }}{K}_2\left(u-\frac{v_1^{*}}{(-2v_0^{*}+K)}\right)dx+dx+{\langle u,f\rangle }_{L^2}\hfill \\ & & -\frac{1}{2\alpha }{\int }_{\mathsf{\Omega }}{\left(v_0^{*}\right)}^{2}dx-\beta {\int }_{\mathsf{\Omega }}{v}_0^{*}dx.\hfill \end{array}$$

(39)

Define also

$$A^{+}=\{u\in V:uf\ge 0,\mathrm{a}.\mathrm{e}.\mathrm{in}\mathsf{\Omega }\},$$

$$V_2={\{u\in V:\parallel u\parallel }_{\infty }\le K_3\},$$

and

$$V_1=V_2\cap A^{+}$$

for some appropriate $K_3>0$ to be specified.

Moreover define

$$B^{*}=\{v_0^{*}\in Y^{*}:\parallel v_0^{*}{\parallel }_{\infty }\le K_4\}$$

and

$$D^{*}=\{v_1^{*}\in Y^{*}:\parallel v_1^{*}\parallel \le K_5\},$$

for some appropriate real constants $K_4,K_5>0$ to be specified.

Observe that we may obtain

$$\frac{{\partial }^2{J}_1^{*}(u,v_1^{*},v_0^{*})}{\partial {\left(v_1^{*}\right)}^2}=K_1+\frac{1}{-2v_0^{*}+K}+\frac{K_2}{{(-2v_0^{*}+K)}^2}$$

$$\frac{{\partial }^2{J}_1^{*}(u,v_1^{*},v_0^{*})}{\partial u^2}=-\gamma {\nabla }^2+K+K_1{(-\gamma {\nabla }^2+K)}^2+K_2$$

and

$$\frac{{\partial }^2{J}_1^{*}(u,v_1^{*},v_0^{*})}{\partial u\partial v_1^{*}}=-1-K_1(-\gamma {\nabla }^2+K)-\frac{K_2}{-2v_0^{*}+K}$$

so that

$$\begin{array}{ccc}& & det\left\{\frac{{\partial }^2{J}_1^{*}(u,v_1^{*},v_0^{*})}{\partial u\partial v_1^{*}}\right\}\hfill \\ & =& \frac{{\partial }^2{J}_1^{*}(u,v_1^{*},v_0^{*})}{\partial {\left(v_1^{*}\right)}^2}\frac{{\partial }^2{J}_1^{*}(u,v_1^{*},v_0^{*})}{\partial u^2}-{\left(\frac{{\partial }^2{J}_1^{*}(u,v_1^{*},v_0^{*})}{\partial u\partial v_1^{*}}\right)}^2\hfill \\ & =& \frac{6K_1(-\gamma {\nabla }^2){(-\gamma {\nabla }^2+2v_0^{*})}^2+(-\gamma {\nabla }^2+2v_0^{*})(-\gamma {\nabla }^2(-2v_0^{*})K_1-v_0-4\gamma {\nabla }^2)}{2{(-\gamma {\nabla }^2+v_0^{*})}^2}\hfill \\ & \equiv & H\left(v_0^{*}\right).\hfill \end{array}$$

(40)

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

Theorem 9.1.

Assume $K=2\gamma {\nabla }^2,$$K_2=10(-\gamma {\nabla }^2)$ and $K_1\gg max\{1,|K|,K_2,K_3,K_4,K_5\}$ and suppose $(u_0,{\widehat{v}}_1^{*},{\widehat{v}}_0^{*})\in V_1\times D^{*}\times B^{*}$ is such that

$$\delta J_1^{*}(u_0,{\widehat{v}}_1^{*},{\widehat{v}}_0^{*})=\mathbf{0}$$

and

$$H\left({\widehat{v}}_0^{*}\right)>\mathbf{0}.$$

Assume also

$$\frac{{\partial }^2{J}_1^{*}(u_0,{\widehat{v}}_1^{*},v_0^{*})}{\partial {\left(v_0^{*}\right)}^2}<\mathbf{0},$$

$\forall v_0^{*}\in B^{*}.$

Under such hypotheses, we have that

$$\delta J\left(u_0\right)=\mathbf{0},$$

and

$$\begin{array}{ccc}\hfill J\left(u_0\right)& =& J_1^{*}(u_0,{\widehat{v}}_1^{*},{\widehat{v}}_0^{*})\hfill \\ & =& \underset{(u,v_1^{*})\in V_1\times D^{*}}{inf}\left\{\underset{v_0^{*}\in B^{*}}{sup}{J}_1^{*}(u,v_1^{*},v_0^{*})\right\}\hfill \\ & =& \underset{v_0^{*}\in B^{*}}{sup}\left\{\underset{(u,v_1^{*})\in V_1\times D^{*}}{inf}{J}_1^{*}(u,v_1^{*},v_0^{*})\right\}.\hfill \end{array}$$

(41)

Proof. 

The proof that

$$\delta J\left(u_0\right)=\mathbf{0}$$

and

$$J\left(u_0\right)=J_1^{*}(u_0,{\widehat{v}}_1^{*},{\widehat{v}}_0^{*})$$

may be easily made similarly as in the previous sections.

Moreover, from the hypotheses,

$$J_1^{*}(u_0,{\widehat{v}}_1^{*},{\widehat{v}}_0^{*})=\underset{(u,v_1^{*})\in V_1\times D^{*}}{inf}{J}_1^{*}(u,v_1^{*},{\widehat{v}}_0^{*})$$

and

$$J_1^{*}(u_0,{\widehat{v}}_1^{*},{\widehat{v}}_0^{*})=\underset{v_0^{*}\in B^{*}}{sup}{J}_1^{*}(u_0,{\widehat{v}}_1^{*},v_0^{*}).$$

From this, a standard saddle point theorem and the remaining hypotheses, we may infer that

$$\begin{array}{ccc}\hfill J\left(u_0\right)& =& J_1^{*}(u_0,{\widehat{v}}_1^{*},{\widehat{v}}_0^{*})\hfill \\ & =& \underset{(u,v_1^{*})\in V_1\times D^{*}}{inf}\left\{\underset{v_0^{*}\in B^{*}}{sup}{J}_1^{*}(u,v_1^{*},v_0^{*})\right\}\hfill \\ & =& \underset{v_0^{*}\in B^{*}}{sup}\left\{\underset{(u,v_1^{*})\in V_1\times D^{*}}{inf}{J}_1^{*}(u,v_1^{*},v_0^{*})\right\}.\hfill \end{array}$$

(42)

The proof is complete. □

10. A Fourth Primal Dual Formulation for a Related Model

Let $\mathsf{\Omega }\subset {\mathbb{R}}^3$ be an open, bounded and connected set with a regular boundary denoted by $\partial \mathsf{\Omega }.$

Consider the functional $J:V\to \mathbb{R}$ where

$$\begin{array}{ccc}\hfill J\left(u\right)& =& \frac{\gamma }{2}{\int }_{\mathsf{\Omega }}\nabla u·\nabla udx+\frac{\alpha }{2}{\int }_{\mathsf{\Omega }}{(u^2-\beta )}^{2}dx\hfill \\ & & -{\langle u,f\rangle }_{L^2},\hfill \end{array}$$

(43)

$\alpha >0,\beta >0,\gamma >0$, $V=W_0^{1,2}\left(\mathsf{\Omega }\right)$ and $f\in L^2\left(\mathsf{\Omega }\right).$

Denoting $Y=Y^{*}=L^2\left(\mathsf{\Omega }\right)$, define now $J_1^{*}:V\times Y^{*}\times Y^{*}\to \mathbb{R}$ by

$$\begin{array}{ccc}\hfill J_1^{*}(u,v_0^{*},v_1^{*})& =& \frac{\gamma }{2}{\int }_{\mathsf{\Omega }}\nabla u·\nabla udx+\frac{1}{2}{\int }_{\mathsf{\Omega }}K{u}^{2}dx\hfill \\ & & -{\langle u,v_1^{*}\rangle }_{L^2}+\frac{1}{2}{\int }_{\mathsf{\Omega }}\frac{{\left(v_1^{*}\right)}^2}{(-2v_0^{*}+K)}dx\hfill \\ & & +\frac{1}{2(\alpha +\epsilon )}{\int }_{\mathsf{\Omega }}{(v_0^{*}-\alpha (u^2-\beta ))}^{2}dx+{\langle u,f\rangle }_{L^2}\hfill \\ & & -\frac{1}{2\alpha }{\int }_{\mathsf{\Omega }}{\left(v_0^{*}\right)}^{2}dx-\beta {\int }_{\mathsf{\Omega }}{v}_0^{*}dx,\hfill \end{array}$$

(44)

where $\epsilon >0$ is a small real constant.

Define also

$$A^{+}=\{u\in V:uf\ge 0,\mathrm{a}.\mathrm{e}.\mathrm{in}\mathsf{\Omega }\},$$

$$V_2={\{u\in V:\parallel u\parallel }_{\infty }\le K_3\},$$

and

$$V_1=V_2\cap A^{+}$$

for some appropriate $K_3>0$ to be specified.

Moreover define

$$B^{*}=\{v_0^{*}\in Y^{*}:\parallel v_0^{*}{\parallel }_{\infty }\le K_4\}$$

and

$$D^{*}=\{v_1^{*}\in Y^{*}:\parallel v_1^{*}\parallel \le K_5\},$$

for some appropriate real constants $K_4,K_5>0$ to be specified.

Observe that, defining

$$\phi =v_0^{*}-\alpha (u^2-\beta )$$

we may obtain

$$\frac{{\partial }^2{J}_1^{*}(u,v_0^{*},v_1^{*})}{\partial u^2}=-\gamma {\nabla }^2+K+\frac{\alpha }{\alpha +\epsilon }4u^2-2\phi \frac{\alpha }{\alpha +\epsilon }$$

$$\frac{{\partial }^2{J}_1^{*}(u,v_0^{*},v_1^{*})}{\partial {\left(v_1^{*}\right)}^2}=\frac{1}{-2v_0^{*}+K}$$

and

$$\frac{{\partial }^2{J}_1^{*}(u,v_0^{*},v_1^{*})}{\partial u\partial v_1^{*}}=-1$$

so that

$$\begin{array}{ccc}& & det\left\{\frac{{\partial }^2{J}_1^{*}(u,v_0^{*},v_1^{*})}{\partial u\partial v_0^{*}}\right\}\hfill \\ & =& \frac{{\partial }^2{J}_1^{*}(u,v_1^{*},v_0^{*})}{\partial {\left(v_1^{*}\right)}^2}\frac{{\partial }^2{J}_1^{*}(u,v_1^{*},v_0^{*})}{\partial u^2}-{\left(\frac{{\partial }^2{J}_1^{*}(u,v_1^{*},v_0^{*})}{\partial u\partial v_1^{*}}\right)}^2\hfill \\ & =& \frac{-\gamma {\nabla }^2+2v_0^{*}+4\frac{{\alpha }^2}{\alpha +\epsilon }{u}^2-2\frac{\alpha }{\alpha +\epsilon }\phi }{-2v_0^{*}+K}\hfill \\ & \equiv & H(u,v_0^{*}).\hfill \end{array}$$

(45)

However, at a critical point, we have $\phi =\mathbf{0}$ so that, we define

$$C_{v_0^{*}}^{*}=\{u\in V:\phi \le \mathbf{0}\}.$$

From such results, assuming $K\gg max\{K_3,K_4,K_5\},$ define now

$$E_{v_0^{*}}=\{u\in V:H(u,v_0^{*})>\mathbf{0}\}.$$

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

Theorem 10.1.

Suppose $(u_0,{\widehat{v}}_1^{*},{\widehat{v}}_0^{*})\in E^{*}=(V_1\cap C_{{\widehat{v}}_0^{*}}\cap E_{{\widehat{v}}_0^{*}})\times D^{*}\times B^{*}$ is such that

$$\delta J_1^{*}(u_0,{\widehat{v}}_0^{*},{\widehat{v}}_1^{*})=\mathbf{0}.$$

Under such hypotheses, we have that

$$\delta J\left(u_0\right)=\mathbf{0}$$

and

$$\begin{array}{ccc}\hfill J\left(u_0\right)& =& \underset{u\in V_1}{inf}J\left(u\right)\hfill \\ & =& J_1^{*}(u_0,{\widehat{v}}_1^{*},{\widehat{v}}_0^{*})\hfill \\ & =& \underset{(u,v_1^{*})\in V_1\times D^{*}}{inf}\left\{\underset{v_0^{*}\in B^{*}}{sup}{J}_1^{*}(u,v_1^{*},v_0^{*})\right\}\hfill \\ & =& \underset{v_0^{*}\in B^{*}}{sup}\left\{\underset{(u,v_1^{*})\in V_1\times D^{*}}{inf}{J}_1^{*}(u,v_1^{*},v_0^{*})\right\}.\hfill \end{array}$$

(46)

Proof. 

The proof that

$$\delta J\left(u_0\right)=\mathbf{0}$$

and

$$J\left(u_0\right)=J_1^{*}(u_0,{\widehat{v}}_1^{*},{\widehat{v}}_0^{*})$$

may be easily made similarly as in the previous sections.

Moreover, from the hypotheses, we have

$$J_1^{*}(u_0,{\widehat{v}}_1^{*},{\widehat{v}}_0^{*})=\underset{(u,v_1^{*})\in V_1\times D^{*}}{inf}{J}_1^{*}(u,v_1^{*},{\widehat{v}}_0^{*})$$

and

$$J_1^{*}(u_0,{\widehat{v}}_1^{*},{\widehat{v}}_0^{*})=\underset{v_0^{*}\in B^{*}}{sup}{J}_1^{*}(u_0,{\widehat{v}}_1^{*},v_0^{*}).$$

From this, from a standard saddle point theorem and the remaining hypotheses, we may infer that

$$\begin{array}{ccc}\hfill J\left(u_0\right)& =& J_1^{*}(u_0,{\widehat{v}}_1^{*},{\widehat{v}}_0^{*})\hfill \\ & =& \underset{(u,v_1^{*})\in V_1\times D^{*}}{inf}\left\{\underset{v_0^{*}\in B^{*}}{sup}{J}_1^{*}(u,v_1^{*},v_0^{*})\right\}\hfill \\ & =& \underset{v_0^{*}\in B^{*}}{sup}\left\{\underset{(u,v_1^{*})\in V_1\times D^{*}}{inf}{J}_1^{*}(u,v_1^{*},v_0^{*})\right\}.\hfill \end{array}$$

(47)

Moreover, observe that

$$\begin{array}{ccc}\hfill J_1^{*}(u_0,{\widehat{v}}_1^{*},{\widehat{v}}_0^{*})& =& \underset{(u,v_1^{*})\in V_1\times D^{*}}{inf}{J}_1^{*}(u,v_1^{*},{\widehat{v}}_0^{*})\hfill \\ & \le & \frac{\gamma }{2}{\int }_{\mathsf{\Omega }}\nabla u·\nabla udx+\frac{K}{2}{\int }_{\mathsf{\Omega }}{u}^{2}dx\hfill \\ & & +\langle u^2,{\widehat{v}}_0^{*}\rangle -\frac{K}{2}{\int }_{\mathsf{\Omega }}{u}^{2}dx\hfill \\ & & -\frac{1}{2\alpha }{\int }_{\mathsf{\Omega }}{\left({\widehat{v}}_0^{*}\right)}^{2}dx-\beta {\int }_{\mathsf{\Omega }}{\widehat{v}}_0^{*}dx\hfill \\ & & +\frac{1}{2(\alpha +\epsilon )}{\int }_{\mathsf{\Omega }}{({\widehat{v}}_0^{*}-\alpha (u^2-\beta ))}^{2}dx-{\langle u,f\rangle }_{L^2}\hfill \\ & \le & \underset{v_0^{*}\in Y^{*}}{sup}\left\{\frac{\gamma }{2}{\int }_{\mathsf{\Omega }}\nabla u·\nabla udx+\langle u^2,v_0^{*}\rangle \right.\hfill \\ & & -\frac{1}{2\alpha }{\int }_{\mathsf{\Omega }}{\left(v_0^{*}\right)}^{2}dx-\beta {\int }_{\mathsf{\Omega }}{v}_0^{*}dx\hfill \\ & & \left.+\frac{1}{2(\alpha +\epsilon )}{\int }_{\mathsf{\Omega }}{(v_0^{*}-\alpha (u^2-\beta ))}^{2}dx-{\langle u,f\rangle }_{L^2}\right\}\hfill \\ & =& \frac{\gamma }{2}{\int }_{\mathsf{\Omega }}\nabla u·\nabla udx+\frac{\alpha }{2}{\int }_{\mathsf{\Omega }}{(u^2-\beta )}^{2}dx\hfill \\ & & -{\langle u,f\rangle }_{L^2},\forall u\in V_1.\hfill \end{array}$$

(48)

Summarizing, we have got

$$J\left(u_0\right)=J_1^{*}(u_0,{\widehat{v}}_1^{*},{\widehat{v}}_0^{*})\le \underset{u\in V_1}{inf}J\left(u\right).$$

From such results, we may infer that

$$\begin{array}{ccc}\hfill J\left(u_0\right)& =& \underset{u\in V_1}{inf}J\left(u\right)\hfill \\ & =& J_1^{*}(u_0,{\widehat{v}}_1^{*},{\widehat{v}}_0^{*})\hfill \\ & =& \underset{(u,v_1^{*})\in V_1\times D^{*}}{inf}\left\{\underset{v_0^{*}\in B^{*}}{sup}{J}_1^{*}(u,v_1^{*},v_0^{*})\right\}\hfill \\ & =& \underset{v_0^{*}\in B^{*}}{sup}\left\{\underset{(u,v_1^{*})\in V_1\times D^{*}}{inf}{J}_1^{*}(u,v_1^{*},v_0^{*})\right\}.\hfill \end{array}$$

(49)

The proof is complete. □