变分问题存在唯一性和对应边值问题的导出

梯度和格林公式

$\phi (x,y,z)$ 是数量函数，那么有$$\nabla \cdot \nabla \phi =(\frac{\partial }{\partial x},\frac{\partial }{\partial y},\frac{\partial }{\partial z})\cdot (\frac{\partial \phi }{\partial x},\frac{\partial \phi }{\partial y},\frac{\partial \phi }{\partial z})=\frac{{\partial }^2\phi }{\partial x^2}+\frac{{\partial }^2\phi }{\partial y^2}+\frac{{\partial }^2\phi }{\partial z^2}=\Delta \phi$$

格林第一公式

$${\int }_{\Omega }(\nabla u\cdot \nabla v+v\Delta u)dx={\int }_{\Gamma }v\frac{\partial u}{\partial n}ds$$

书上第四章的格林公式

$\Omega$ 是二维区域
$${\iint }_{\Omega }(u\Delta v-v\Delta u)dxdy={\int }_{\partial \Omega }(u\frac{\partial v}{\partial \mathbf{n}}-v\frac{\partial u}{\partial \mathbf{n}})dl$$其中 $\mathbf{n}$ 是 $\partial \Omega$ 的单位外法向量

题目一

设 $\Omega$ 是 ${\mathbb{R}}^2$ 中的一个有界开区域，其光滑边界用 $\Gamma$ 表示。$\Omega$ 上定义的泛函为：
$$J(v)=\frac{1}{2}||v|{|}_{H_1(\Omega )}^2+{\int }_{\Gamma }(\frac{1}{2}\sigma v^2-gv)ds$$其中 $\sigma \in C(\Gamma )，g\in C(\Gamma )$，且 $0<{\sigma }_0<\sigma <{\sigma }_1<+\infty$。变分问题的描述为：求 $u\in H^1(\Omega )$，使得 $$J(u)=mi{n}_{v\in H^1(\Omega )}J(v)$$请详细地证明这个变分问题的解存在并且唯一

$||v|{|}_{H_1(\Omega )}^2$ 表示函数 $v$ 在 $H_1(\Omega )$ 空间中的范数的平方。

$H_1(\Omega )$ 是一个Sobolev空间，通常用来表示具有某些平滑性的函数空间。对于一个定义在域 $\Omega$ 上的函数 $v$，其在 $H_1(\Omega )$ 空间中的范数定义为：

$||v|{|}_{H_1(\Omega )}={\left({\int }_{\Omega }|v{|}^2+|\nabla v{|}^{2}dx\right)}^{1/2}$

其中，$|v{|}^2$ 表示函数 $v$ 的平方，$|\nabla v{|}^2$ 表示 $v$ 的梯度的平方（即各个偏导数的平方和），${\int }_{\Omega }$ 表示对整个域 $\Omega$ 的积分。

因此，$||v|{|}_{H_1(\Omega )}^2$ 就是上述范数的平方，它考虑了函数 $v$ 本身的大小和它的变化率（通过梯度）。

证明：
(1) 证明泛函 $J$ 在 $H^1(\Omega )$ 上有下确界

为了证明泛函 $J$ 在 $H^1(\Omega )$ 上有下确界，我们需要证明对于任意 $v\in H^1(\Omega )$，$J(v)$ 的值有一个下界。根据泛函 $J$ 的定义：

$$J(v)=\frac{1}{2}||v|{|}_{H_1(\Omega )}^2+{\int }_{\Gamma }(\frac{1}{2}\sigma v^2-gv)ds$$

首先，我们考虑泛函的第一部分，即 $||v|{|}_{H_1(\Omega )}^2$，这是 $H^1(\Omega )$ 中元素的平方范数，显然是非负的。所以，这一部分对泛函 $J(v)$ 的值起着一个非负的贡献。

所以$$J(v)=\frac{1}{2}||v|{|}_{H_1(\Omega )}^2+{\int }_{\Gamma }(\frac{1}{2}\sigma v^2-gv)ds$$$$\ge {\int }_{\Gamma }(\frac{1}{2}\sigma v^2-gv)ds={\int }_{\Gamma }\frac{1}{2}(\sigma v^2-2\sqrt{\sigma }v\cdot g\frac{1}{\sqrt{\sigma }}+\frac{g^2}{\sigma }）-\frac{g^2}{2\sigma }ds$$$$={\int }_{\Gamma }\frac{1}{2}(\sqrt{\sigma }v-g\frac{1}{\sqrt{\sigma }}{)}^2-\frac{g^2}{2\sigma }ds$$$$\ge -{\int }_{\Gamma }-\frac{1}{2\sigma }{g}^{2}ds(因为\frac{1}{2}(\sqrt{\sigma }v-g\frac{1}{\sqrt{\sigma }}{)}^2是非负项)$$

所以 $J$ 在 $H_1(\Omega )$ 有下确界，不妨设 $m={\inf }_{v\in H^1(\Omega )}J(v)$，那么由下确界的定义，$\forall k\in {\mathbb{N}}_{+},\exists u_k\in H^1(\Omega )$ 使得$$m\le J(u_k)=\frac{1}{2}||u_k|{|}_{H_1(\Omega )}^2+{\int }_{\Gamma }(\frac{1}{2}\sigma u_k^2-g{u}_k)ds<m+\frac{1}{k}$$

(2) 证明 { u k } \{u_k\} {uk​} 为 $H^1(\Omega )$ 中的基本列，也即变分问题解的存在性

首先，根据引理 $3.4$ 的等式，有 $$J(u_k)+J(u_l)-2J(\frac{u_k+u_l}{2})=||\frac{\nabla (u_k-u_l)}{2}|{|}_{L_2(\Omega )}^2\ge 0$$

$\forall k,l\in {\mathbb{N}}_{+}$ 有
$$J(u_k)+J(u_l)-2J(\frac{u_k+u_l}{2})\le m+\frac{1}{k}+m+\frac{1}{l}-2m=\frac{1}{k}+\frac{1}{l}$$这里是因为第一步证明出了下确界后，根据下确界的定义有 $J(u_k)<m+\frac{1}{k}$ 同理对 $u_l$ 也是，而又因为 $J$ 的下确界是 $m$ 所以$$2J(\frac{u_k+u_l}{2})\ge 2m$$$$-2J(\frac{u_k+u_l}{2})\le -2m$$ 所以得到上面的不等式

得到了 $J(u_k)+J(u_l)-2J(\frac{u_k+u_l}{2})$ 的不等式关系后，再把其利用平行四边形等式进行展开，

并将 $dx$ 和 $ds$ 进行合并化简，拼凑 $||u_k-u_l|{|}_{L_2(\Omega )}^2$

这里贴一下公式 $$||u|{|}_{H_1(\Omega )}^2={\int }_{\Omega }{u}^{2}dx+{\int }_{\Omega }|\nabla u{|}^{2}dx$$$x$ 是向量，这里把泛函 $J(v)$ 的表达式形式用这种形式表示$$J(v)=\frac{1}{2}{\int }_{\Omega }{v}^2+|\nabla v{|}^{2}dx+{\int }_{\Gamma }\frac{1}{2}\sigma v^2-gvds$$把 $J(u_k),J(u_l),-2J(\frac{u_k+u_l}{2})$ 分别带进去，计算量比较大，化简后凑完全平方 $(u_k-u_l{)}^2$ 最终得到

$$J(u_k)+J(u_l)-2J(\frac{u_k+u_l}{2})=\frac{1}{4}{\int }_{\Omega }(u_k-u_l{)}^2+|\nabla (u_k-u_l){|}^{2}dx$$$$+\frac{1}{4}\sigma {\int }_{\Gamma }(u_k-u_l{)}^{2}ds=\frac{1}{4}||u_k-u_l|{|}_{H^1(\Omega )}^2+\frac{1}{4}\sigma {\int }_{\Gamma }(u_k-u_l{)}^{2}ds$$因为上面推导得$$J(u_k)+J(u_l)-2J(\frac{u_k+u_l}{2})\le \frac{1}{k}+\frac{1}{l}$$所以也即$$\frac{1}{4}||u_k-u_l|{|}_{H^1(\Omega )}^2+\frac{1}{4}\sigma {\int }_{\Gamma }(u_k-u_l{)}^{2}ds\le \frac{1}{k}+\frac{1}{l}$$又因为$\frac{1}{4}\sigma {\int }_{\Gamma }(u_k-u_l{)}^{2}ds$是非负项，可得 $$0\le ||u_k-u_l|{|}_{H^1(\Omega )}^2\le 4(\frac{1}{k}+\frac{1}{l})\to 0(k,l\to \infty )$$所以 { u k } \{u_k\} {uk​} 是 $H^1(\Omega )$ 的基本列

由 $H^1(\Omega )$ 的完备性知，$\exists u\in H^1(\Omega ),s.t.||u_k-u|{|}_{H^1(\Omega )}\to 0(k\to \infty )$，

从而 $||u_k-u|{|}_{L_2(\Omega )}\to 0(k\to \infty )$ 这是因为根据 $H^1(\Omega )$ 和 $L_2(\Omega )$ 各自范数定义的数学表达式，$H^1(\Omega )$ 范数是两个平方项之和取 $1/2$ 次幂，而 $L^2(\Omega )$ 范数是一个平方项取 $1/2$ 次幂且该项在 $H^1(\Omega )$ 范数中也出现

从而 $||u_k-u|{|}_{L_2(\Omega )}\to 0(k\to \infty )$

对这个等式$$m\le J(u_k)=\frac{1}{2}||u_k|{|}_{H_1(\Omega )}^2+{\int }_{\Gamma }(\frac{1}{2}\sigma u_k^2-g{u}_k)ds<m+\frac{1}{k}$$取极限 $k\to \infty$

得到 $J(u)=m={\inf }_{v\in H^1(\Omega )}J(v)$

(3) 证明变分问题解的唯一性

假设 $u,w$ 都是变分问题的解，即 $J(u)=J(w)=m={\inf }_{v\in H^1(\Omega )}J(v)$

那么按照 (2) 中凑出的表达式 $$0\le ||\frac{\nabla (u-w)}{2}|{|}_{L_2(\Omega )}^2=J(u)+J(w)-2J(\frac{u+w}{2})$$$$=m+m-2J(\frac{u+w}{2})$$这里是因为 $u,w$ 都是变分问题的解，即 $J(u)=J(w)=m$

而 $J(\frac{u+w}{2})\ge m$，因为 $m$ 是下确界，所以 $-2J(\frac{u+w}{2})\le -2m$

所以 $$0\le ||\frac{\nabla (u-w)}{2}|{|}_{L_2(\Omega )}^2=J(u)+J(w)-2J(\frac{u+w}{2})$$$$=m+m-2J(\frac{u+w}{2})\le m+m-2m=0$$

所以 $u=w$

然后导出与这个变分问题等价的边值问题

若 $u$ 是该变分问题的解， $v\in H_0^1(\Omega )$，令$$j(ϵ)=J(u+ϵv)=\frac{1}{2}{\int }_{\Omega }(u+ϵv{)}^2+|\nabla (u+ϵv){|}^{2}dx+{\int }_{\Gamma }\frac{1}{2}\sigma (u+ϵv{)}^2-g(u+ϵv)ds$$对 $ϵ$ 求导

附注 $|\nabla (u+ϵv){|}^2$ 对 $ϵ$ 求导的过程

$\nabla (u+ϵv)=\nabla u+\nabla ϵv=\nabla u+ϵ\nabla v$

$|\vec{\mathbf{x}}{|}^2=\vec{\mathbf{x}}\cdot \vec{\mathbf{x}}$，向量模的平方等于自身和自身做点积
$$|\nabla (u+ϵv){|}^2=(\nabla u+ϵ\nabla v)\cdot (\nabla u+ϵ\nabla v)$$

展开后得到：

$$|\nabla (u+ϵv){|}^2=\nabla u\cdot \nabla u+2ϵ\nabla u\cdot \nabla v+{ϵ}^2\nabla v\cdot \nabla v$$

对 $ϵ$ 求一阶导数，我们得到：

$$\frac{\partial }{\partial ϵ}|\nabla (u+ϵv){|}^2=2\nabla u\cdot \nabla v+2ϵ\nabla v\cdot \nabla v$$

当 $ϵ=0$ 时，我们有：

$${\frac{\partial }{\partial ϵ}|\nabla (u+ϵv){|}^2|}_{ϵ=0}=2\nabla u\cdot \nabla v$$

对 $ϵ$ 求导得$$j^{\prime }(ϵ)=\frac{1}{2}{\int }_{\Omega }2(u+ϵv)v+2\nabla u\cdot \nabla v+2ϵ|\nabla v{|}^{2}dx+{\int }_{\Gamma }\sigma (u+ϵv)v-gvds$$那么$$j^{\prime }(0)={\int }_{\Omega }uv+\nabla u\cdot \nabla vdx+{\int }_{\Gamma }(\sigma u-g)vds$$$$={\int }_{\Omega }(u-\Delta u)vdx+{\int }_{\Gamma }(\sigma u-g+\frac{\partial u}{\partial n})vds$$

这一步的推导过程：
这个等式是通过应用格林公式和部分积分来推导的。格林公式将一个域内的散度项转化为边界上的通量项。我们可以按照以下步骤来推导这个等式：

首先，考虑函数 $u$ 和 $v$ 在定义在域 $\Omega$ 上，并且 $\Gamma$ 是 $\Omega$ 的边界。

1. 应用格林公式:

   根据格林第一公式，$${\int }_{\Omega }(\nabla u\cdot \nabla v+v\Delta u)dx={\int }_{\Gamma }v\frac{\partial u}{\partial n}ds$$我们有：
   $${\int }_{\Omega }\nabla u\cdot \nabla vdx=-{\int }_{\Omega }v\Delta udx+{\int }_{\Gamma }v\frac{\partial u}{\partial n}ds$$
   其中，$\Delta u$ 是 $u$ 的拉普拉斯算子，$\frac{\partial u}{\partial n}$ 是 $u$ 沿边界法线方向的导数。

2. 整合项:

   将上述公式的右侧加到等式的左侧，我们得到：
   $${\int }_{\Omega }uvdx+{\int }_{\Omega }\nabla u\cdot \nabla vdx+{\int }_{\Gamma }(\sigma u-g)vds={\int }_{\Omega }uvdx-{\int }_{\Omega }v\Delta udx+{\int }_{\Gamma }v\frac{\partial u}{\partial n}ds+{\int }_{\Gamma }(\sigma u-g)vds$$

   这里，我们添加了 ${\int }_{\Omega }uvdx$ 和 ${\int }_{\Gamma }(\sigma u-g)vds$ 这两项，它们分别代表了在域 $\Omega$ 和边界 $\Gamma$ 上的积分。

3. 简化:

   现在，我们可以将等式的右侧重新组织，得到：
   $${\int }_{\Omega }uv+\nabla u\cdot \nabla vdx+{\int }_{\Gamma }(\sigma u-g)vds={\int }_{\Omega }(u-\Delta u)vdx+{\int }_{\Gamma }(\sigma u-g+\frac{\partial u}{\partial n})vds$$
   这就是最终的等式。

这个推导利用了格林公式和部分积分的基本性质，将微分算子转化为边界上的积分，从而简化了问题的处理。

由于 $u$ 是变分问题的解，根据必要性条件可知 $j^{\prime }(0)=0$ 即
$$j^{\prime }(0)={\int }_{\Omega }uv+\nabla u\cdot \nabla vdx+{\int }_{\Gamma }(\sigma u-g)vds$$$$={\int }_{\Omega }(u-\Delta u)vdx+{\int }_{\Gamma }(\sigma u-g+\frac{\partial u}{\partial n})vds=0$$

在这个表达式中，$j^{\prime }(0)$ 是某个泛函 $j$ 在 $0$ 处的导数。根据给出的条件，对于所有在 $\Omega$ 上光滑且在边界上消失的函数 $v$（即 $v\in C_0^{\infty }(\Omega )$），都有：

$${\int }_{\Omega }(u-\Delta u)vdx+{\int }_{\Gamma }(\sigma u-g+\frac{\partial u}{\partial n})vds=0$$

由于 $v$ 在边界 $\Gamma$ 上为零（即 $v{|}_{\Gamma }=0$），边界上的积分项消失，我们得到：

$${\int }_{\Omega }(u-\Delta u)vdx=0$$

这个等式对于所有 $v\in C_0^{\infty }(\Omega )$ 都成立。而任意给定的 $v$ 不恒为 $0$，所以$$u-\Delta u=0\text{in}\Omega$$根据变分原理，如果一个线性泛函对于一个充分大的函数空间中的所有函数都为零，那么被积函数必须在该空间的对偶空间中几乎处处为零。在这种情况下，因为 $C_0^{\infty }(\Omega )$ 是一个在 $\Omega$ 上的密集子空间，我们可以推断出：

$$u-\Delta u=0\text{in}\Omega$$

这意味着在 $\Omega$ 上有 ${\int }_{\Omega }(u-\Delta u)vdx=0$ 成立。这个结果表明 $u$ 是一个满足泊松方程 $u-\Delta u=0$ 的函数。

从上一个推导步骤中，我们已经知道对于所有 $v\in C_0^{\infty }(\Omega )$，有：

$${\int }_{\Omega }(u-\Delta u)vdx=0$$

并且我们推断出 $u-\Delta u=0$ 在 $\Omega$ 内成立。现在，我们回到原始的等式：

$${\int }_{\Omega }(u-\Delta u)vdx+{\int }_{\Gamma }(\sigma u-g+\frac{\partial u}{\partial n})vds=0$$

由于我们已经知道第一个积分项为零，所以只剩下边界项：

$${\int }_{\Gamma }(\sigma u-g+\frac{\partial u}{\partial n})vds=0$$

这个等式对所有 $v\in C_0^{\infty }(\Omega )$ 成立，但是 $C_0^{\infty }(\Omega )$ 中的函数在边界上是零。为了让这个等式对于边界上非零的函数也成立，我们需要考虑一个更广泛的函数空间，比如在边界上也光滑的函数空间。

然而，即使在这种情况下，为了让上述积分为零，被积函数（即边界项中的 $\sigma u-g+\frac{\partial u}{\partial n}$）必须在边界 $\Gamma$ 上几乎处处为零。因此，我们可以推导出：

$$\sigma u-g+\frac{\partial u}{\partial n}=0\text{on}\Gamma$$

这表明 $u$ 在边界 $\Gamma$ 上满足这个边界条件。结合内部的泊松方程 $u-\Delta u=0$，我们得到了一个完整的边值问题。

所以变分问题的等价边值问题为：
$$\{\begin{array}{ll}u-\Delta u=0& \text{in }\Omega ,\\ \sigma u+\frac{\partial u}{\partial n}{|}_{\Gamma }=g& \end{array}$$

题目二

设 $\Omega$ 是 ${\mathbb{R}}^2$ 中的有界开区域，其光滑边界为 $\Gamma$，$\Omega$ 上定义的泛函为：$$J(v)={\iint }_{\Omega }\frac{1}{2}[(\frac{\partial v}{\partial x}{)}^2+(\frac{\partial v}{\partial y}{)}^2]dxdy+{\int }_{\Gamma }(\frac{1}{2}\sigma v^2-gv)ds$$变分问题的描述为：求 $u\in V,s.t.$$$J(u)=\underset{v\in V}{\min }J(v)$$其中 $V=C^2(\Omega )\cap C^1(\overline{\Omega })$ 请导出与此变分问题等价的边值问题

解：

令 $j(ϵ)=J(u+ϵv)={\iint }_{\Omega }\frac{1}{2}\left[{\left(\frac{\partial (u+ϵv)}{\partial x}\right)}^2+{\left(\frac{\partial (u+ϵv)}{\partial y}\right)}^2\right]dxdy+{\int }_{\Gamma }\left(\frac{1}{2}\sigma (u+ϵv{)}^2-g(u+ϵv)\right)ds$

对 $ϵ$ 求导，使用链式法则和乘积法则：

$\frac{\partial }{\partial ϵ}{\left(\frac{\partial (u+ϵv)}{\partial x}\right)}^2=2\left(\frac{\partial (u+ϵv)}{\partial x}\right)\left(\frac{\partial v}{\partial x}\right)$

同样，

$\frac{\partial }{\partial ϵ}{\left(\frac{\partial (u+ϵv)}{\partial y}\right)}^2=2\left(\frac{\partial (u+ϵv)}{\partial y}\right)\left(\frac{\partial v}{\partial y}\right)$

因此，$\Omega$ 上双重积分的导数为：

${\iint }_{\Omega }\left[\left(\frac{\partial (u+ϵv)}{\partial x}\right)\left(\frac{\partial v}{\partial x}\right)+\left(\frac{\partial (u+ϵv)}{\partial y}\right)\left(\frac{\partial v}{\partial y}\right)\right]dxdy$

对于 $\Gamma$ 上的单重积分部分，我们有：

$\frac{1}{2}\sigma (u+ϵv{)}^2-g(u+ϵv)$

对 $ϵ$ 求导，我们得到：

$\sigma (u+ϵv)v-gv$

因此，$\Gamma$ 上单重积分的导数为：

${\int }_{\Gamma }\left(\sigma (u+ϵv)v-gv\right)ds$

将这两部分合并，我们得到 $j^{\prime }(ϵ)$ 的表达式：

$j^{\prime }(ϵ)={\iint }_{\Omega }\left[\left(\frac{\partial (u+ϵv)}{\partial x}\right)\left(\frac{\partial v}{\partial x}\right)+\left(\frac{\partial (u+ϵv)}{\partial y}\right)\left(\frac{\partial v}{\partial y}\right)\right]dxdy+{\int }_{\Gamma }\left(\sigma (u+ϵv)v-gv\right)ds$

那么 $j^{\prime }(0)={\iint }_{\Omega }\left[\left(\frac{\partial u}{\partial x}\right)\left(\frac{\partial v}{\partial x}\right)+\left(\frac{\partial u}{\partial y}\right)\left(\frac{\partial v}{\partial y}\right)\right]dxdy+{\int }_{\Gamma }\left(\sigma uv-gv\right)ds$

$$v_x{u}_x+v_y{u}_y=\nabla u\cdot \nabla v=(u_x,u_y)\cdot (v_x,v_y)$$

即$j^{\prime }(0)={\iint }_{\Omega }\nabla u\cdot \nabla vdxdy+{\int }_{\Gamma }\left(\sigma uv-gv\right)ds$

利用格林第一公式可以得到 $$j^{\prime }(0)=-{\iint }_{\Omega }v\Delta udxdy+{\int }_{\Gamma }\frac{\partial u}{\partial n}vds+{\int }_{\Gamma }\left(\sigma uv-gv\right)ds$$

如果我们取 $v\in C_0^{\infty }(\Omega )$，即 $v$ 在 $\Omega$ 内部光滑且在边界 $\Gamma$ 上为零，那么边界项都将消失，因为 $v{|}_{\Gamma }=0$。

因此，我们有
$$-{\iint }_{\Omega }v\Delta udxdy=0$$
对于所有的 $v\in C_0^{\infty }(\Omega )$。

由于 $v$ 是任意的，根据基本引理（du Bois-Reymond Lemma），我们可以推导出
$$-\Delta u=0,\text{ in }\Omega$$
这就是所谓的拉普拉斯方程，表示 $u$ 在 $\Omega$ 内部是调和函数。

当我们得到 $-\Delta u=0,\text{ in }\Omega$ 之后，如果我们取 $v\in C_0^{\infty }(\partial \Omega )$，即 $v$ 在边界 $\partial \Omega$ 上光滑且在 $\Omega$ 的内部为零，我们可以进一步分析边界条件。

由于 $v$ 在 $\Omega$ 的内部为零，内部的积分项 ${\iint }_{\Omega }v\Delta udxdy$ 自然为零。因此，变分公式简化为：
$$j^{\prime }(0)={\int }_{\Gamma }\frac{\partial u}{\partial n}vds+{\int }_{\Gamma }\left(\sigma uv-gv\right)ds=0$$

由于 $v$ 是在边界 $\partial \Omega$ 上任意的，我们可以推断出：
$$\frac{\partial u}{\partial n}+\sigma u-g=0,\text{ on }\partial \Omega$$

这就是所谓的诺依曼边界条件或第三类边界条件，它描述了函数 $u$ 在边界上的导数（即法向导数）与 $u$ 本身和给定函数 $g$ 的关系。