本篇文章, 我们继续对文献[16]中的问题进行深入研究, 对Dirichlet-Robin混合边值齐次化问题解的$ H^{1} $和$ L^{2} $收敛率进行分析.具体地, 令$ \Omega $是$ {{\Bbb R}} ^{n} $中有界的$ C^{1, 1} $区域, $ n\geq 3 $.假设$ \partial \Omega = \Gamma_{1}\bigcup\Gamma_{2} $, 这里的$ \Gamma_{1} $和$ \Gamma_{2} $是$ \partial \Omega $中的两个互不相交的子集.不妨设$ u_{\varepsilon}\in H^{1}(\Omega) $是以下问题的弱解

(1.1) $ \begin{equation} \left\{\begin{array}{ll} L_{\varepsilon}u_{\varepsilon} = -{\rm div}(A(x/\varepsilon)\triangledown u_{\varepsilon}) = -\frac{\partial}{\partial x_{i}}\left(a_{ij}(x/\varepsilon)\frac{\partial u_{\varepsilon} }{\partial x_{j}}\right) = F, \; \; &x\in\Omega, \\[3pt] u_{\varepsilon} = f, \; \; &x\in\Gamma_{1}, \\[3pt] \frac{\partial u_{\varepsilon}}{\partial \nu_{\varepsilon}}+k(x/\varepsilon)u_{\varepsilon} = g, \; \; &x\in\Gamma_{2}.\\[3pt] \end{array} \right. \end{equation} $

全篇文章使用了求和约定.这里假设矩阵$ A(y) = (a_{ij}(y)) $, 其中$ 1\leq i, j\leq n $, 函数$ k(y) $满足以下条件.

$ \bullet $周期性条件:对于任意的$ y\in {{\Bbb R}} ^{n} $, $ Y = [0, 1)^{n}\simeq {{\Bbb R}} ^{n}/{\mathbb Z}^{n} $, 有

(1.2) $ \begin{equation} A(y+Y) = A(y), \; \; k(y+Y) = k(y). \end{equation} $

$ \bullet $强制性条件:存在$ \lambda>0 $, 对于所有的$ y\in {{\Bbb R}} ^{n} $, 有

(1.3) $ \begin{equation} \lambda \mid\xi\mid ^{2}\leq a_{ij}(y)\xi_{i}\xi_{j}\leq \frac{1}{\lambda} \mid\xi\mid ^{2}. \end{equation} $

$ \bullet $光滑性条件:

(1.4) $ \begin{equation} F\in L^{2}(\Omega), \; f\in H^{1}(\partial\Omega), \; g\in H^{1}(\partial\Omega), \; k\in C^{\infty}({{\Bbb R}} ^{n}). \end{equation} $

$ \bullet $凸性条件:

(1.5) $ \begin{equation} \Omega\; \mbox{是紧的、一致凸光滑区域.} \end{equation} $

$ \bullet $相容性条件:

(1.6) $ \begin{equation} \int_{\partial\Omega}u_{\varepsilon}(x){\rm d}\sigma(x) = \int_{\Omega}F(x){\rm d}x+\int_{\Gamma_{2}}g(x){\rm d}\sigma(x) = 0. \end{equation} $

伴随着问题(1.1)的齐次化问题是

(1.7) $ \begin{equation} \left\{\begin{array}{ll} L_{0}u_{0} = -q_{ij}\frac{\partial^{2}u_{0}}{\partial x_{i}\partial x_{j}} = F, \; \; &x\in\Omega, \\[3pt] u_{0} = f, \; \; &x\in\Gamma_{1}, \\[3pt] \frac{\partial u_{0}}{\partial \nu_{0}}+k_{0}u_{0} = g, \; \; &x\in\Gamma_{2}, \\[3pt] \end{array} \right. \end{equation} $

其中$ k_{0} = \int_{Y}k(y){\rm d}y $.

令$ H_{0}^{1}(\Omega;\Gamma_{1}) $表示$ \Omega $中的光滑函数在$ H^{1}(\Omega) $中的闭包, 它们在$ \Gamma_{1} $上为0.则$ u_{\varepsilon}\in H^{1}(\Omega) $被称为问题(1.1)的弱解, 即在$ \Gamma_{1} $上$ u_{\varepsilon} = f $, 并且对任意的$ \phi\in H^{1}_{0}(\Omega;\Gamma_{1}) $, 都有

$ \begin{eqnarray*} \int_{\Omega}A(x/\varepsilon)\triangledown u_{\varepsilon}\cdot \triangledown \phi {\rm d}x+\int_{\Gamma_{2}}k(x/\varepsilon) u_{\varepsilon}\phi {\rm d}\sigma(x) = \int_{\Omega}F\phi {\rm d}x+\int_{\Gamma_{2}}g\phi {\rm d}\sigma(x). \end{eqnarray*} $

由Lax-Milgram定理, 结合条件(1.2)–(1.6), 可以得到问题(1.1)弱解的存在性和唯一性.众所周知, 当$ \varepsilon\rightarrow 0 $时, $ u_{\varepsilon}\rightharpoonup u_{0} $在$ H^{1}(\Omega) $中弱收敛.具体证明可参阅文献[4].

最近, 关于齐次化中的各种问题, 学者们研究得到了许多收敛率结果. 2011年, Gérard和Masmoudi^[3]针对边界分层的Neumann问题, 获得了解的$ L^{2} $收敛结果. 2012年, Kenig、Lin和Shen^[8]针对椭圆系统的Dirichlet或Neumann齐次化问题建立了$ L^{2} $及$ H^{\frac{1}{2}} $收敛估计. 2013年, Aleksanyan、Shahgholian和Sjölin^[1-2]针对光滑一致凸区域上固定算子、振荡Dirichlet边值问题, 研究得到了解的逐点收敛和$ L^{p} $收敛率. 2014年, Kenig、Lin和Shen^[9]对振荡Neumann边值齐次化问题, 利用振荡积分的估计, 得到了解的$ W^{k, p} $收敛率. 2015年, 本文作者针对固定算子、振荡Neumann边值的齐次化问题, 利用调和分析的方法, 证明得到了解的逐点收敛和$ W^{1, p} $收敛率^[15]. 2015年, Gu^[5]针对线性Stokes系统, 证明得到了解的$ H^{1} $和$ L^{2} $收敛. 2016年, Shen^[12]对于Dirichlet-Neumann混合边值的齐次化问题, 证明得到了解的$ H^{1} $和$ L^{2} $收敛率. 2017年, Gu和Shen^[6]对于线性Stokes问题, 通过Green函数的渐近性估计, 得到了解的$ L^{p} $与$ L^{\infty} $收敛率. 2018年, Niu和Xu^[7]针对高阶抛物方程的齐次化问题, 获得了经典的$ L^{2} $收敛估计.

2019年, 本文作者针对Dirichlet-Robin混合边值的齐次化问题, 研究得到了解的$ H^{1} $及$ L^{2} $收敛率^[16].那篇文章, Robin边界条件并没有出现振荡因子, 即$ k $是常数.但本工作却需要对边界上出现的振荡因子进行重新估计, 因此出现了本质上的困难.

本工作的主要难点在于两方面.一方面, 边界上出现的振荡因子$ k(x/\varepsilon) $, 它的出现必然导致$ \Gamma_{2} $上伴随的振荡积分更加难以估计.另一方面, 边界交叉项也会使得一阶逼近函数的估计变得更加困难.本文的研究受到了Geng和Zhuge^[4]的研究作为启发, 同时借鉴了Suslina^[13-14]的对偶方法和光滑算子作为研究手段.本文的创新之处在于将对偶方法和光滑算子延拓到混合边值问题的情形.

以下是本文的主要结果.

定理1.1  令$ \Omega $是$ {{\Bbb R}} ^{n} $中紧的、一致凸区域.假设$ u_{\varepsilon}\in H^{1}(\Omega) $、$ u_{0}\in H^{2}(\Omega) $分别是振荡Robin边值问题(1.1)和(1.7)的两个弱解.则在假设条件(1.2)–(1.6)下, 存在常数C, 对于任意的$ 0<\delta<3/2(n-1) $, 都有

$ \begin{eqnarray*} \parallel u_{\varepsilon}-u_{0}-\varepsilon\chi T_{\varepsilon}(\triangledown \widetilde{u}_{0})\parallel_{H^{1}(\Omega)}\leq C\varepsilon^{\frac{3}{2(n-1)}-\delta}\parallel u_{0}\parallel_{H^{2}(\Omega)}, \end{eqnarray*} $

其中$ T_{\varepsilon} $是光滑算子, $ \chi $是校正函数, $ \widetilde{u}_{0} $是$ u_{0} $的延拓.

定理1.2  在定理1.1的假设条件下, 存在常数C, 对于任意的$ 0<\delta<3/2(n-1) $, 满足

$ \begin{eqnarray*} \parallel u_{\varepsilon}-u_{0}\parallel_{L^{2}(\Omega)}\leq C\varepsilon^{\frac{3}{2(n-1)}-\delta}\parallel u_{0}\parallel_{H^{2}(\Omega)}. \end{eqnarray*} $

本文组织如下:第2部分, 包含了一些基本定义和有用的命题.第3部分, 利用对偶方法和光滑算子, 得到了混合边值问题解的$ H^{1} $和$ L^{2} $收敛率.最后一部分, 总结了文章的主要结果, 并给出了一些相关评论.

我们首先介绍一些定义.

令$ B_{r}(x) $表示球心是$ x $点, 半径是$ r $的开球.

$ \Omega_{\varepsilon} = \{x\in {{\Bbb R}} ^{n}: dist(x, \partial\Omega)\leq \varepsilon\}. $

因为$ \Omega $是光滑区域, 所以存在有界延拓算子$ E:H^{2}(\Omega)\rightarrow H^{2}({{\Bbb R}} ^{n}) $, 满足$ \widetilde{u}_{0} $是$ u_{0} $的延拓, 且有

$ \parallel \widetilde{u}_{0}\parallel_{H^{2}({{\Bbb R}} ^{n})}\leq C\parallel u_{0}\parallel_{H^{2}(\Omega)}. $

令$ H_{0}^{1}(\Omega;\Gamma_{1}) $表示$ \Omega $中的光滑函数在$ H^{1}(\Omega) $中的闭包, 并且在$ \Gamma_{1} $上为0.这里还假设$ \varphi\in H^{1}_{0}(\Omega;\Gamma_{1}) $, 且满足

$ \parallel \varphi\parallel_{H^{1}({{\Bbb R}} ^{n})}\leq C\parallel \varphi\parallel_{H^{1}(\Omega)}. $

全篇文章, $ C $表示正常数, 其在不同公式中可能取值不同.

伴随着问题(1.1)中的振荡算子$ L_{\varepsilon} $, 齐次化算子为

(2.1) $ \begin{equation} L_{0} = -q_{ij}\frac{\partial^{2}}{\partial x_{i}\partial x_{j}}, \; \; \; \; x\in\Omega. \end{equation} $

常数矩阵$ Q = (q_{ij}) $可表示为

$ \begin{eqnarray*} \begin{array}{rl} q_{ij} = & \int_{Y}\left[a_{ij}(y)+a_{ik}(y)\frac{\partial \chi^{j}(y)}{\partial y_{k}}\right]{\rm d}y, \end{array} \end{eqnarray*} $

其中$ y\in Y = [0, 1)^{n}\simeq {{\Bbb R}} ^{n}/{\mathbb Z}^{n} $.校正函数$ \chi(y) = (\chi^{j}(y)) $是以下细胞问题的解

$ \begin{eqnarray*} \left\{\begin{array}{ll} \frac{\partial}{\partial y_{i}}\left[a_{ik}(y)\frac{\partial \chi^{j}(y)}{\partial y_{k}}+a_{ij}(y)\right] = 0, \; \; \; \; y\in Y, \\[3pt] \chi^{j}(y+Y) = \chi^{j}(y), \\[3pt] \int_{Y}\chi^{j}(y){\rm d}y = 0. \end{array} \right. \end{eqnarray*} $

命题2.1  令$ F_{ij}(y)\in L^{2}(Y) $表示周期函数, 其中$ Y = [0, 1)^{n} $.假设$ \int_{Y}F_{ij}(y){\rm d}y = 0 $和$ \frac{\partial}{\partial y_{i}}(F_{ij}(y)) = 0 $.则存在$ \Phi_{kij}\in H^{1}(Y) $, 满足$ F_{ij} = \frac{\partial \Phi_{kij} }{\partial y_{k}} $并有$ \Phi_{kij} = -\Phi_{ikj} $.

上述命题是众所周知的, 读者可以参阅文献[8].如果假定

$ F_{ij}(y) = q_{ij}-a_{ij}(y)-a_{ik}(y)\frac{\partial\chi^{j}(y)}{\partial y_{k}}. $

注意到, 周期函数$ F_{ij}(y) $满足$ \int_{Y} F_{ij}(y){\rm d}y = 0 $和$ \frac{\partial}{\partial y_{i}}\left(F_{ij}\right) = 0 $.从上述命题可知, 存在函数$ \Phi_{kij}(y) $, 满足$ \Phi_{kij} = -\Phi_{ikj} $和$ F_{ij} = \frac{\partial \Phi_{kij}}{{\partial y_{k}}}. $这个命题描述了齐次化过程中矩阵$ Q $和$ A $的特殊关系, 它对于分析解的渐近性十分重要.

接下来, 我们介绍光滑算子的一些性质.作为本文的一个重要研究动机, 试图将光滑算子的用法推广到振荡Robin混合边值齐次化问题的情形.

固定$ \phi\in C^{\infty}_{0}(B_{1}(0)) $, 使其满足$ \phi \geq 0 $和$ \int_{{{\Bbb R}} ^{n}}\phi {\rm d}x = 1 $.在$ L^{2} $范数下定义$ T_{\varepsilon} $如下

$ T_{\varepsilon}(u)(x) = u\ast\phi_{\varepsilon} = \int_{{{\Bbb R}} ^{n}}u(x-y)\phi_{\varepsilon}(y){\rm d}y, $

其中$ \phi_{\varepsilon}(x) = \varepsilon^{-n}\phi(x/\varepsilon) $.通常称$ T_{\varepsilon} $为$ \varepsilon $尺度的光滑算子.

命题2.2  令$ u_{0}\in H^{2}({{\Bbb R}} ^{n}) $, 周期函数$ f\in L^{2}({{\Bbb R}} ^{n}) $.则有

$ \parallel \triangledown u_{0}-T_{\varepsilon}(\triangledown u_{0})\parallel_{L^{2}({{\Bbb R}} ^{n})}\leq C\varepsilon\parallel \triangledown^{2} u_{0}\parallel_{L^{2}({{\Bbb R}} ^{n})}, $

$ \parallel f(x/\varepsilon)T_{\varepsilon}(u_{0})\parallel_{L^{2}({{\Bbb R}} ^{n})}\leq C\parallel f\parallel_{L^{2}({{\Bbb R}} ^{n})}\parallel u_{0}\parallel_{L^{2}({{\Bbb R}} ^{n})} $

和

$ \parallel T_{\varepsilon}(\triangledown^{2} u_{0})\parallel_{L^{2}({{\Bbb R}} ^{n})}\leq C\parallel \triangledown^{2} u_{0}\parallel_{L^{2}({{\Bbb R}} ^{n})}. $

证  利用Parseval定理和Hölder不等式, 可得以上结果, 可参见文献[11].

命题2.3  如果$ u_{0}\in H^{2}({{\Bbb R}} ^{n}) $, 并且$ f\in L^{2}({{\Bbb R}} ^{n}) $是周期函数.则有

$ \parallel T_{\varepsilon}(\triangledown u_{0}) \parallel_{L^{2}(\Omega_{\varepsilon})} \leq C\varepsilon^{1/2}\parallel u_{0}\parallel_{H^{2}({{\Bbb R}} ^{n})} $

和

$ \parallel f(x/\varepsilon)T_{\varepsilon}(u_{0})\parallel_{L^{2}(\Omega_{\varepsilon})}\leq C\varepsilon^{1/2}\parallel f\parallel_{L^{2}({{\Bbb R}} ^{n})}\parallel u_{0}\parallel_{H^{1}({{\Bbb R}} ^{n})}. $

证  本估计可利用Fubini定理来证明.读者可参阅文献[10]或[12].

下面一个命题, 是奇异振荡积分中的经典结论.

命题2.4  令$ \Omega $表示$ {{\Bbb R}} ^{n} $中紧的、一致凸光滑区域.假设$ h(x/\varepsilon) $是光滑周期函数, $ v_{\varepsilon} $是以下问题的解

$ \left\{\begin{array}{ll} -{\rm div}Q(\triangledown v_{\varepsilon}) = 0, \; \; \; \; &x\in\Omega, \\[3pt] \frac{\partial v_{\varepsilon}}{\partial \nu_{0}} = h(x/\varepsilon), \; \; \; \; \; \; \; &x\in\partial\Omega. \end{array} \right. $

则由假设

$ \int_{\partial\Omega}v_{\varepsilon}{\rm d}\sigma(x) = \int_{\partial\Omega}h(x/\varepsilon){\rm d}\sigma(x) = 0, $

有

$ \parallel {v_\varepsilon }{\parallel _{{L^p}(\Omega )}} \le C{\varepsilon ^{1/p}}\;\;\;{\rm{和}}\;\;\;\parallel {v_\varepsilon }{\parallel _{{L^p}(\Omega )}} \le C{\varepsilon ^{\frac{1}{p} - \delta }},$

其中$ 1\leq p<\infty $并且$ 0<\delta<1/p $ .

证  这些估计是文献[1]中的主要结果.

本部分的目标是证明混合边值问题解的$ H^{1} $和$ L^{2} $收敛率.

令一阶逼近函数形式为

(3.1) $ \begin{equation} M_{\varepsilon} = u_{\varepsilon}-u_{0}-\varepsilon\chi T_{\varepsilon}(\triangledown \widetilde{u}_{0}). \end{equation} $

则对于任意的$ \varphi\in H^{1}_{0}(\Omega;\Gamma_{1}) $, 都有

$ \int_{\Omega}A(x/\varepsilon)\triangledown u_{\varepsilon}\cdot \triangledown \varphi {\rm d}x+ \int_{\partial\Omega}k(x/\varepsilon)u_{\varepsilon} \varphi {\rm d}\sigma(x) = \int_{\Omega}Q\triangledown u_{0}\cdot \triangledown \varphi {\rm d}x+ \int_{\partial\Omega}k_{0}u_{0} \varphi {\rm d}\sigma(x). $

显然, 对于任意的$ \varphi\in H_{0}^{1}(\Omega;\Gamma_{1}) $, 都有

$ \int_{\partial\Omega} k(x/\varepsilon)u_{\varepsilon}\varphi {\rm d}\sigma(x) = \int_{\Gamma_{2}} k(x/\varepsilon)u_{\varepsilon}\varphi {\rm d}\sigma(x). $

于是就有

$ \begin{eqnarray*} &&\int_{\Omega}A(x/\varepsilon)\triangledown M_{\varepsilon}\cdot \triangledown \varphi {\rm d}x+\int_{\partial\Omega} k(x/\varepsilon)M_{\varepsilon}\varphi {\rm d}\sigma(x)\\ & = &\int_{\Omega}[Q\triangledown u_{0} -A(x/\varepsilon)\triangledown u_{0}-\varepsilon A(x/\varepsilon)\triangledown(\chi T_{\varepsilon}(\triangledown \widetilde{u}_{0}))]\cdot \triangledown \varphi {\rm d}x\\ &&-\int_{\partial\Omega} \varepsilon k(x/\varepsilon) \chi T_{\varepsilon} (\triangledown \widetilde{u}_{0})\varphi {\rm d}\sigma(x)+\int_{\partial\Omega}(k_{0}-k(x/\varepsilon))u_{0}\varphi {\rm d}\sigma. \end{eqnarray*} $

直接计算, 可得

$ \begin{eqnarray*} && Q\triangledown u_{0} -A(x/\varepsilon)\triangledown u_{0}-\varepsilon A(x/\varepsilon)\triangledown (\chi T_{\varepsilon}(\triangledown \widetilde{u}_{0}))\\ & = &[Q(\triangledown u_{0}-T_{\varepsilon}(\triangledown \widetilde{u}_{0}))-A(x/\varepsilon)(\triangledown u_{0}-T_{\varepsilon}(\triangledown \widetilde{u}_{0}))-\varepsilon A(x/\varepsilon) \chi T_{\varepsilon}(\triangledown^{2} \widetilde{u}_{0})]\\ &&+[Q-A(x/\varepsilon)-A(x/\varepsilon)\triangledown \chi]T_{\varepsilon}(\triangledown \widetilde{u}_{0}). \end{eqnarray*} $

于是, 可得

(3.2) $ \begin{eqnarray} &&\int_{\Omega}A(x/\varepsilon)\triangledown M_{\varepsilon}\cdot \triangledown \varphi {\rm d}x+\int_{\partial\Omega} k(x/\varepsilon)M_{\varepsilon}\varphi {\rm d}\sigma(x) {}\\ & = &\int_{\Omega}[Q(\triangledown u_{0}-T_{\varepsilon}(\triangledown \widetilde{u}_{0}))-A(x/\varepsilon)(\triangledown u_{0}-T_{\varepsilon}(\triangledown \widetilde{u}_{0})) -\varepsilon A(x/\varepsilon) \chi T_{\varepsilon}(\triangledown^{2} \widetilde{u}_{0})]\cdot \triangledown \varphi {\rm d}x {}\\ &&+\int_{\Omega}[Q-A(x/\varepsilon)-A(x/\varepsilon)\triangledown \chi]T_{\varepsilon}(\triangledown \widetilde{u}_{0})\cdot \triangledown \varphi {\rm d}x {}\\ &&-\int_{\partial\Omega} \varepsilon k(x/\varepsilon) \chi T_{\varepsilon}(\triangledown \widetilde{u}_{0})\varphi {\rm d}\sigma(x) +\int_{\partial\Omega} [k_{0}-k(x/\varepsilon)]u_{0}\cdot\varphi {\rm d}\sigma(x) {}\\ &\doteq &I_{1}+I_{2}+I_{3}+I_{4}. \end{eqnarray} $

首先, 来控制$ I_{1} $项, 观察可得

(3.3) $ \begin{eqnarray} \mid I_{1}\mid & \leq& C\int_{\Omega}|\triangledown u_{0}-T_{\varepsilon}(\triangledown \widetilde{u}_{0})|\cdot |\triangledown \varphi| {\rm d}x+C\varepsilon\int_{\Omega}| T_{\varepsilon}(\triangledown^{2} \widetilde{u}_{0})|\cdot |\triangledown \varphi| {\rm d}x {}\\ &\leq& C \parallel \triangledown \widetilde{u}_{0}-T_{\varepsilon}(\triangledown \widetilde{u}_{0})\parallel_{L^{2}({{\Bbb R}} ^{n})}\parallel \triangledown \varphi\parallel_{L^{2}(\Omega)}+C\varepsilon\parallel T_{\varepsilon}(\triangledown^{2} \widetilde{u}_{0}) \parallel_{L^{2}({{\Bbb R}} ^{n})}\parallel \triangledown \varphi\parallel_{L^{2}(\Omega)} {}\\ &\leq& C\varepsilon \parallel u_{0}\parallel _{H^{2}(\Omega)}\parallel \varphi\parallel_{H^{1}(\Omega)}, \end{eqnarray} $

这里, 用到了命题2.2.

接下来, 处理$ I_{2} $项.注意到

$ F_{ij} = q_{ij}-a_{ij}(y)-a_{ik}(y)\frac{\partial\chi^{j}(y)}{\partial y_{k}}. $

$ F_{ij} $是周期函数, 并且满足命题2.1.则存在周期函数

$ \Phi = (\Phi_{kij})\in H^{1}(Y), $

满足

$ \Phi_{kij} = -\Phi_{ikj}\ \mbox{和}{\quad}   q_{ij}-a_{ij}(y)-a_{ik}(y)\frac{\partial\chi^{j}(y)}{\partial y_{k}} = \frac{\partial \Phi_{kij} }{\partial y_{k}}. $

此时, 假设$ \eta_{\varepsilon}\in C_{0}^{\infty}({{\Bbb R}} ^{n}) $是截断函数, 并且满足当$ x\in \Omega_{\varepsilon} $时$ \eta_{\varepsilon}(x)\equiv1 $; 当$ x\not\in \Omega_{2\varepsilon} $时$ \eta_{\varepsilon}(x)\equiv0 $.因此, 利用散度定理和命题2.1, 可得

$ \begin{eqnarray*} I_{2}& = &\int_{\Omega}(1-\eta_{\varepsilon})[Q-A(x/\varepsilon)-A(x/\varepsilon)\triangledown \chi] T_{\varepsilon}(\triangledown \widetilde{u}_{0})\cdot \triangledown \varphi {\rm d}x\\ &&+\int_{\Omega}\eta_{\varepsilon}[Q-A(x/\varepsilon)-A(x/\varepsilon)\triangledown \chi] T_{\varepsilon}(\triangledown \widetilde{u}_{0})\cdot \triangledown \varphi {\rm d}x \\ & = &\int_{\Omega}\frac{\partial}{\partial x_{k}}(\varepsilon \Phi_{kij})(1-\eta_{\varepsilon}) T_{\varepsilon}\left(\frac{\partial \widetilde{u}_{0}}{\partial x_{j}}\right)\cdot \frac{\partial \varphi}{\partial x_{i}}{\rm d}x\\ &&+\int_{\Omega_{2\varepsilon}} \eta_{\varepsilon}[Q-A(x/\varepsilon)-A(x/\varepsilon)\triangledown \chi] T_{\varepsilon}(\triangledown \widetilde{u}_{0}) \cdot \triangledown \varphi {\rm d}x\\ & = &-\int_{\Omega}\varepsilon \Phi_{kij}(1-\eta_{\varepsilon}) T_{\varepsilon}\left(\frac{\partial^{2} \widetilde{u}_{0}} {\partial x_{k}\partial x_{j}}\right)\frac{\partial \varphi}{\partial x_{i}}{\rm d}x\\ &&-\int_{\Omega}\varepsilon \Phi_{kij}\frac{\partial}{\partial x_{k}}(1-\eta_{\varepsilon}) T_{\varepsilon}\left(\frac{\partial \widetilde{u}_{0}}{\partial x_{j}}\right)\frac{\partial \varphi}{\partial x_{i}}{\rm d}x \\ &&-\int_{\Omega}\varepsilon \Phi_{kij}(1-\eta_{\varepsilon}) T_{\varepsilon}\left(\frac{\partial \widetilde{u}_{0}}{ \partial x_{j}}\right)\frac{\partial^{2} \varphi}{\partial x_{k}\partial x_{i}}{\rm d}x\\ &&+\int_{\Omega_{2\varepsilon}} \eta_{\varepsilon}[Q-A(x/\varepsilon)-A(x/\varepsilon)\triangledown \chi] T_{\varepsilon}(\triangledown \widetilde{u}_{0})\cdot \triangledown \varphi {\rm d}x \\ & = &-\int_{\Omega}\varepsilon \Phi_{kij}(1-\eta_{\varepsilon}) T_{\varepsilon}\left(\frac{\partial^{2} \widetilde{u}_{0}}{\partial x_{k}\partial x_{j}}\right)\frac{\partial \varphi}{\partial x_{i}}{\rm d}x \\ &&-\int_{\Omega}\varepsilon \Phi_{kij}\frac{\partial}{\partial x_{k}}(1-\eta_{\varepsilon}) T_{\varepsilon}\left(\frac{\partial \widetilde{u}_{0}}{\partial x_{j}}\right)\frac{\partial \varphi}{\partial x_{i}}{\rm d}x\\ &&+\int_{\Omega_{2\varepsilon}} \eta_{\varepsilon}[Q-A(x/\varepsilon)-A(x/\varepsilon)\triangledown \chi] T_{\varepsilon} (\triangledown \widetilde{u}_{0})\cdot \triangledown \varphi {\rm d}x \\ &\doteq& I_{21}+I_{22}+I_{23}, \end{eqnarray*} $

其中, 利用了函数$ \Phi_{kij} $的反对称性.

为了估计$ I_{21} $, 利用命题2.2, 可得

$ \begin{eqnarray*} \mid I_{21}\mid &\leq & C\varepsilon\int_{\Omega}\mid T_{\varepsilon}(\triangledown^{2} \widetilde{u}_{0})\mid\cdot\mid\triangledown\varphi\mid {\rm d}x\\ &\leq &C\varepsilon\parallel T_{\varepsilon}(\triangledown^{2} \widetilde{u}_{0}) \parallel_{L^{2}({{\Bbb R}} ^{n})}\parallel \triangledown \varphi\parallel_{L^{2}(\Omega)}\\ &\leq& C\varepsilon \parallel u_{0}\parallel_{H^{2}(\Omega)}\parallel \triangledown \varphi\parallel_{L^{2}(\Omega)}. \end{eqnarray*} $

要估计$ I_{22} $, 利用命题2.3, 就有

$ \begin{eqnarray*} \mid I_{22}\mid &\leq& C\int_{\Omega_{2\varepsilon}}\mid \Phi(x/\varepsilon) T_{\varepsilon}(\triangledown \widetilde{u}_{0})\mid\cdot\mid\triangledown\varphi\mid {\rm d}x\\ &\leq& C\varepsilon^{1/2}\parallel T_{\varepsilon}(\triangledown^{2} \widetilde{u}_{0}) \parallel_{L^{2}({{\Bbb R}} ^{n})}\parallel \triangledown \varphi\parallel_{L^{2}(\Omega)}\\ &\leq& C\varepsilon^{1/2} \parallel u_{0}\parallel_{H^{2}(\Omega)}\parallel \triangledown \varphi\parallel_{L^{2}(\Omega)}. \end{eqnarray*} $

接下来, 处理$ I_{23} $.类似地

$ \begin{eqnarray*} \mid I_{23}\mid &\leq& C\int_{\Omega_{2\varepsilon}}\mid [Q-A(x/\varepsilon) -A(x/\varepsilon)\triangledown \chi]T_{\varepsilon}(\triangledown \widetilde{u}_{0})\mid\cdot\mid\triangledown\varphi\mid {\rm d}x\\ &\leq& C\varepsilon^{1/2}\parallel T_{\varepsilon}(\triangledown \widetilde{u}_{0}) \parallel_{L^{2}({{\Bbb R}} ^{n})}\parallel \triangledown \varphi\parallel_{L^{2}(\Omega)}\\ &\leq& C\varepsilon^{1/2} \parallel u_{0}\parallel_{H^{2}(\Omega)}\parallel \triangledown \varphi\parallel_{L^{2}(\Omega)}. \end{eqnarray*} $

因此, 我们就有

(3.4) $ \begin{equation} \mid I_{2}\mid \leq C\varepsilon^{1/2} \parallel u_{0}\parallel_{H^{2}(\Omega)}\parallel \triangledown \varphi\parallel_{L^{2}(\Omega)}. \end{equation} $

下面, 处理$ I_{3} $.利用函数$ k $的光滑性, 再使用迹定理, 可得

(3.5) $ \begin{eqnarray} \mid I_{3}\mid & \leq &\varepsilon\int_{\partial\Omega} \mid T_{\varepsilon}(\bigtriangledown \widetilde{u}_{0}) \mid\cdot\mid \varphi \mid {\rm d}\sigma(x) {}\\ & \leq& C\varepsilon\left( \parallel T_{\varepsilon}(\bigtriangledown \widetilde{u}_{0}) \parallel_{L^{2}({{\Bbb R}} ^{n})}+\parallel T_{\varepsilon}(\bigtriangledown^{2} \widetilde{u}_{0}) \parallel_{L^{2}({{\Bbb R}} ^{n})}\right)\parallel \varphi\parallel_{H^{1}(\Omega)} {}\\ & \leq &C\varepsilon \parallel u_{0}\parallel_{H^{2}(\Omega)}\parallel \varphi\parallel_{H^{1}(\Omega)}. \end{eqnarray} $

还剩余$ I_{4} $的估计, 这是本工作的主要困难.正如前文所述, 这里将利用对偶的方法解决此问题.鉴于此, 需要分析以下振荡积分

$ I_{4} = \int_{\partial\Omega} [k_{0}-k(x/\varepsilon)]u_{0}\cdot\varphi {\rm d}\sigma(x). $

这里假设$ m_{\varepsilon} $是以下Neumann问题的解

$ \left\{\begin{array}{ll} -{\rm div}Q(\triangledown m_{\varepsilon}) = 0, \; \; \; \; \; \; \; \; &x\in\Omega, \\[3pt] \frac{\partial m_{\varepsilon}}{\partial \nu_{\varepsilon}} = k_{0}-k(x/\varepsilon), \; \; \; \; \; & x\in\partial\Omega, \end{array} \right. $

$ w_{\varepsilon} $是以下Dirichlet问题的解

$ \left\{\begin{array}{ll} -{\rm div}Q(\triangledown w_{\varepsilon}) = 0, \; \; \; \; \; \; \; \; &x\in\Omega, \\ w_{\varepsilon} = u_{0}\cdot\varphi, \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; & x\in\partial\Omega. \end{array} \right. $

利用散度定理, 可得

$ I_{4} = \int_{\partial\Omega}\frac{\partial m_{\varepsilon}}{\partial \nu_{\varepsilon}}\cdot w_{\varepsilon}{\rm d}\sigma(x) = \int_{\partial\Omega}m_{\varepsilon}\cdot\frac{\partial w_{\varepsilon}}{\partial \nu_{\varepsilon}} {\rm d}\sigma(x). $

由Hölder不等式和迹定理, 就有

$ \begin{eqnarray*} \mid I_{4}\mid&\leq &C \parallel m_{\varepsilon}\parallel_{L^{p}(\partial\Omega)}\parallel \triangledown w_{\varepsilon}\parallel_{L^{q}(\partial\Omega)}\\ & \leq &C \left(\parallel m_{\varepsilon}\parallel_{L^{p}(\Omega)}+\parallel \triangledown m_{\varepsilon}\parallel_{L^{p}(\Omega)}\right)\cdot\parallel u_{0}\varphi\parallel_{W^{1, q}(\partial\Omega)}, \end{eqnarray*} $

其中$ 1/p+1/q = 1 $.

由Sobolev不等式, 有

$ \parallel u_{0}\varphi\parallel_{W^{1, q}(\partial\Omega)} \leq C\parallel u_{0}\parallel_{H^{3/2}(\partial\Omega)} \cdot \parallel \varphi\parallel_{H^{1}(\partial\Omega)}, $

其中$ 1/q = 1-3/2(n-1) $.

因此, 利用命题2.4, 可得对于任意的$ 0<\delta<1/p = 3/2(n-1) $, 可以得到

(3.6) $ \begin{eqnarray} \mid I_{4}\mid&\leq &C \left[\varepsilon^{\frac{3}{2(n-1)}}+\varepsilon^{\frac{3}{2(n-1)}-\delta}\right]\parallel u_{0}\parallel_{H^{3/2}(\partial\Omega)} \parallel \varphi\parallel_{H^{1}(\partial\Omega)}{}\\ & \leq& C \varepsilon^{\frac{3}{2(n-1)}-\delta}\parallel u_{0}\parallel_{H^{2}(\Omega)} \parallel \varphi\parallel_{H^{1}(\Omega)}. \end{eqnarray} $

最后, 结合(3.3)–(3.6)式, 则对于任意的$ \varphi\in H^{1}_{0}(\Omega;\Gamma_{1}) $, 有

(3.7) $ \begin{eqnarray} && \left| \int_{\Omega}A(x/\varepsilon)\triangledown M_{\varepsilon}\cdot \triangledown \varphi {\rm d}x+\int_{\partial\Omega} k(x/\varepsilon)M_{\varepsilon}\varphi {\rm d}\sigma(x) \right|{}\\ &\leq & C\varepsilon^{1/2}\parallel u_{0}\parallel_{H^{2}(\Omega)}\parallel \varphi \parallel_{H^{1}(\Omega)}+C\varepsilon^{\frac{3}{2(n-1)}-\delta}\parallel u_{0}\parallel_{H^{2}(\Omega)}\parallel \varphi \parallel_{H^{1}(\Omega)}{}\\ & \leq& C \varepsilon^{\frac{3}{2(n-1)}-\delta}\parallel u_{0}\parallel_{H^{2}(\Omega)}\parallel \varphi \parallel_{H^{1}(\Omega)}, \end{eqnarray} $

其中$ 0<\delta<3/2(n-1) $.

接下来, 假设

$ R_{\varepsilon} = \varepsilon\eta_{\varepsilon}\chi T_{\varepsilon}(\triangledown \widetilde{u}_{0}) $

和

$ S_{\varepsilon} = M_{\varepsilon}-R_{\varepsilon} = u_{\varepsilon}-u_{0}+\varepsilon(1-\eta_{\varepsilon})\chi T_{\varepsilon}(\triangledown \widetilde{u}_{0}), $

其中$ \eta_{\varepsilon} $是估计$ I_{2} $过程中出现的截断函数.

观察可得, $ S_{\varepsilon}\in H_{0}^{1}(\Omega;\Gamma_{1}) $.利用(3.7)式的结论可得

(3.8) $ \begin{equation} \left| \int_{\Omega}A(x/\varepsilon)\triangledown M_{\varepsilon}\cdot \triangledown S_{\varepsilon} {\rm d}x+\int_{\partial\Omega} k(x/\varepsilon)M_{\varepsilon}S_{\varepsilon} {\rm d}\sigma(x) \right|\leq C\varepsilon^{\frac{3}{2(n-1)}-\delta}\parallel u_{0}\parallel_{H^{2}(\Omega)}\parallel S_{\varepsilon}\parallel_{H^{1}(\Omega)}. \end{equation} $

利用命题2.2, 容易得出

$ \parallel R_{\varepsilon}\parallel_{H^{1}(\Omega)}\leq C\varepsilon^{1/2}\parallel u_{0}\parallel_{H^{2}(\Omega)}. $

因此, 利用(3.8)式和迹定理, 我们就有

(3.9) $ \begin{eqnarray} \parallel S_{\varepsilon}\parallel^{2}_{H^{1}(\Omega)}&\leq& C\parallel R_{\varepsilon}\parallel _{H^{1}(\Omega)}\parallel S_{\varepsilon}\parallel_{H^{1}(\Omega)}+C\varepsilon^{\frac{3}{2(n-1)} -\delta}\parallel u_{0}\parallel_{H^{2}(\Omega)}\parallel S_{\varepsilon} \parallel_{H^{1}(\Omega)}{}\\ &\leq& C\varepsilon^{1/2}\parallel u_{0}\parallel_{H^{2}(\Omega)}\parallel S_{\varepsilon}\parallel_{H^{1} (\Omega)}+C\varepsilon^{\frac{3}{2(n-1)}-\delta}\parallel u_{0}\parallel_{H^{2}(\Omega)}\parallel S_{\varepsilon} \parallel_{H^{1}(\Omega)}{}\\ &\leq &C\varepsilon^{\frac{3}{2(n-1)}-\delta}\parallel u_{0}\parallel_{H^{2}(\Omega)}\parallel S_{\varepsilon} \parallel_{H^{1}(\Omega)}. \end{eqnarray} $

于是, 就有

(3.10) $ \begin{eqnarray} \parallel M_{\varepsilon}\parallel_{H^{1}(\Omega)}&\leq& C\parallel S_{\varepsilon}\parallel_{H^{1}(\Omega)}+C\parallel R_{\varepsilon}\parallel_{H^{1}(\Omega)}{}\\ &\leq& C\varepsilon^{\frac{3}{2(n-1)}-\delta}\parallel u_{0}\parallel_{H^{2}(\Omega)}. \end{eqnarray} $

这就完成了定理1.1的证明.

由定理1.1和命题2.2, 再结合Minkowski不等式, 就有

(3.11) $ \begin{eqnarray} \parallel u_{\varepsilon}-u_{0}\parallel_{L^{2}(\Omega)}& \leq& C\varepsilon ^{\frac{3}{2(n-1)}-\delta}\parallel u_{0}\parallel_{H^{2}(\Omega)}+\parallel\varepsilon \chi T_{\varepsilon}(\bigtriangledown u_{0})\parallel_{L^{2}(\Omega)}{}\\ & \leq &C\varepsilon^{\frac{3}{2(n-1)}-\delta}\parallel u_{0}\parallel_{H^{2}(\Omega)}. \end{eqnarray} $

这就完成了定理1.2的证明.

本文中, 我们研究了Robin混合边值齐次化问题解的收敛率.通过齐次化技术、光滑算子和对偶方法, 我们解决了振荡积分和边界交叉项带来的困难, 最终得到了齐次化的一些定量结论.最终, 本文获得了解的$ H^{1} $和$ L^{2} $收敛率.我们的工作将对偶方法和光滑算子延拓到了混合边值问题的情形.

文中已经提到, Dirichlet和Neumann问题在诸多文献中已经有了相当多的理论结果.那么很自然的就会想到将Dirichlet或Neumann问题的相关结论推广到Robin问题上.众所周知, Robin问题在数学物理上有着很重要的作用, 如流体中多尺度材料的热传导问题, 再如工程材料表面的对流和相变问题等.

和经典齐次化中的结果相比, 文章获得的收敛率还不是最优的.我们知道, $ H^{1} $范数下收敛率的阶数应该是$ O({\varepsilon^{1/2}}) $.和一般的Dirichlet或Neumann边值问题不同的是, 本文研究的振荡Robin边值上出现了振荡因子, 这在一定程度上消耗了部分收敛速度.所以对此类问题, 如何得到最优的收敛率, 是下一步需要重点研究的问题.

虽然目前的工作还是纯理论上的, 但是仍然希望所得的结果能够更好地帮助理解Robin边值问题以及奇异振荡积分.