[论文解读] Scalar boundary value problems on junctions of thin rods and plates. I. Asymptotic analysis and error estimates

Renata Bunoiu, Giuseppe Cardone+1 more|arXiv (Cornell University)|Mar 24, 2014

Advanced Mathematical Modeling in Engineering参考文献 37被引用 10

一句话总结

本文针对具有对比材料特性的薄板与细杆连接结构上的标量边值问题，发展了渐近分析与误差估计，采用小厚度参数 $ h $ 和导热率对比参数 $ \gamma_j(h) = \gamma_j h^{-\alpha} $ 进行建模。研究推导了在 Neumann 与 Dirichlet 边界条件下解在各向异性加权 Sobolev 范数下的收敛速率，表明当 $ \alpha = 1 $ 时收敛率为 $ O(|\ln h|^{-1}) $，当 $ \alpha = 0 $ 时为 $ O(h^{3/2}) $，并给出了包含边界层与降维问题的显式渐近展开式。

ABSTRACT

We derive asymptotic formulas for the solutions of the mixed boundary value problem for the Poisson equation on the union of a thin cylindrical plate and several thin cylindrical rods. One of the ends of each rod is set into a hole in the plate and the other one is supplied with the Dirichlet condition. The Neumann conditions are imposed on the whole remaining part of the boundary. Elements of the junction are assumed to have contrasting properties so that the small parameter, i.e. the relative thickness, appears in the differential equation, too, while the asymptotic structures crucially depend on the contrastness ratio. Asymptotic error estimates are derived in anisotropic weighted Sobolev norms.

研究动机与目标

• 分析具有对比材料特性的薄板与细杆连接结构上 Poisson 方程解的渐近行为。
• 在各向异性加权 Sobolev 范数下，为解的渐近逼近推导严格的误差估计。
• 表征对比参数 $ \alpha $（0 或 1）对渐近结构与收敛速率的影响。
• 建立并确定在连接处与杆端附近边界层的存在性与形式，尤其针对 $ \alpha = 1 $ 的情形。
• 通过渐近展开与传输条件，证明从三维模型向一维与二维模型降维的合理性。

提出的方法

• 在薄板 $ \Omega_\bullet(h) $ 与细杆 $ \Omega_j(h) $ 的连接结构上，建立 Poisson 方程的混合边值问题，其中大部分边界采用 Neumann 条件，杆端采用 Dirichlet 条件。
• 引入表示相对厚度的小参数 $ h $，以及用于建模不同材料导热率的对比参数 $ \gamma_j(h) = \gamma_j h^{-\alpha} $。
• 通过匹配渐近展开法构造解 $ u(h,x) $ 的渐近展开式，对 $ \alpha = 1 $ 的情形引入边界层校正项 $ A_j(\ln h) \widetilde{\mathbf{W}}^j $。
• 推导极限问题：在板区域 $ \omega_0 $ 上为二维 Dirichlet 问题，在每根杆 $ (0,l_j) $ 上为一维问题，其解分别为 $ U_0^\# $ 与 $ U_j^\# $。
• 利用加权 Sobolev 范数与 Kondratiev 理论控制奇异性，并在 $ L^2 $ 与 $ H^1 $ 类型范数下推导误差估计。
• 通过能量估计与各向异性不等式，证明三维全解向降维的二维与一维模型的收敛性。

实验结果

研究问题

• RQ1 Poisson 方程在薄板与细杆连接结构上的解的渐近行为如何依赖于对比参数 $ \alpha $？
• RQ2 在连接处与杆端附近，解的结构是怎样的，尤其对 $ \alpha = 1 $ 的情形，边界层是如何形成的？
• RQ3 渐近逼近在 $ h $ 上的收敛速率是多少，$ \alpha = 0 $ 与 $ \alpha = 1 $ 两种情形下的收敛速率有何不同？
• RQ4 在各向异性加权 Sobolev 范数下的误差估计如何反映几何奇点与材料对比的影响？
• RQ5 三维问题能否有效降维为耦合的二维-一维模型，该降维的精度如何？

主要发现

• 当 $ \alpha = 1 $ 时，解在 $ L^2 $ 与 $ H^1 $ 范数下强收敛于降维的二维-一维模型，收敛率为 $ O(|\ln h|^{-1}) $，误差估计为 $ \| \nabla_x (u_0 - U_0^0 - \sum_j \chi_{0j} A_j \widetilde{\mathbf{W}}^j) \|_{L^2(\Omega_\bullet)} + h^{-1} \sum_j \| \nabla_x (u_j - U_j^0) \|_{L^2(\Omega_j)} \leq c h \mathcal{N} $。
• 当 $ \alpha = 0 $ 时，收敛率为 $ O(h^{3/2}) $，误差估计为 $ \| \nabla_x (u_0 - U_0^\#) \|_{L^2(\Omega_\bullet)} + \sum_j \| \nabla_x (u_j - U_j^\#) \|_{L^2(\Omega_j)} \leq c h^{3/2} \mathcal{N} $。
• 杆端的 Dirichlet 条件导致 Green 函数出现对数奇点，从而在板与杆上的 $ H^1 $-范数中产生 $ O(|\ln h|^{-1}) $ 的收敛率。
• 边界层校正项 $ A_j(\ln h) \widetilde{\mathbf{W}}^j $ 对 $ \alpha = 1 $ 至关重要，能够捕捉连接处与杆端附近的奇异行为。
• 极限问题为在 $ \omega_0 $ 上的二维 Dirichlet 问题与在每根杆 $ (0,l_j) $ 上的一维混合问题，其解满足 $ U_0^\# \in H^2(\omega_0) $ 与 $ U_j^\# \in H^2(0,l_j) $。
• 在板区域 $ \omega_0 \times (0,1) $ 上的 $ H^1 $-范数收敛率为 $ O(|\ln h|^{-1/2}) $，反映出梯度中对数奇点的影响。