边界层

不知道精确解时，怎样猜到层宽 ​

假设左端变化发生在宽度 $\delta$ 内，而函数改变量为数量级 1。于是这一带的导数量级约为

为了让被小系数压低的二阶项重新与一阶项抗衡，需要

令 $X=x/\epsilon$、$Y(X)=y(\epsilon X)$。原方程变成

首阶层方程 $Y_{0,XX}+Y_{0,X}=0$ 的解为 $A+B{e}^{-X}$，确有向区间内部衰减的分量。这同时检查了“导数足够大”与“能接回远处状态”两个要求。

主导平衡也会给出别的宽度。对 $\epsilon y^{″}-y=0$，若函数幅度仍为数量级 1，平衡是 $\epsilon /{\delta }^2\asymp 1$，所以 $\delta \asymp \sqrt{\epsilon }$。边界层没有一律等于参数本身的宽度。

为什么不是右端有层 ​

对原方程试取右端坐标 $X=(1-x)/\epsilon$，首阶方程变成

非恒定项为 $e^X$。向区间内部走对应 $X\to +\mathrm{\infty }$，它会增长而非消失，无法连接一个有界外解。因此本例应让外解满足右端条件，把左端条件交给层。

若一阶项的符号反过来，方程改为 $\epsilon y^{″}-y^{\prime }=1$，衰减分量出现在右端，层位置也随之反转。对变系数方程，若一阶系数变号或为零，还可能出现转向点与内部层；刚才的常系数判断不能直接覆盖这些位置。

层宽、容差与导数各自告诉什么 ​

以指数校正 $e^{-x/\epsilon }$ 为参照，要求它不超过 $\eta$ 时，需要 $x\ge \epsilon \log \left(1/\eta \right)$。固定容差下，这是 $O(\epsilon )$ 厚度；若容差为 $\epsilon$，厚度变为 $\epsilon \log \left(1/\epsilon \right)$。报告有效层宽时，应连同容差一起报告。

式(1)还给出

函数值在整段上保持有界，左端斜率却为 $O({\epsilon }^{-1})$。所以函数值近似准确，不代表导数同样准确；选取误差范数时必须说明是否包含导数。

从识别薄层到求解问题 ​

匹配渐近展开理路匹配渐近展开Matched asymptotic expansions · Matching asymptotic expansions在内外近似共同有效的重叠区匹配常数，扣除重复部分构造复合解，并单独验证其全区间误差。会用重叠区决定 $A,B$，再把内外近似合成一条整段公式。本页完成的则是更早的判断：层在哪里、需要放大多少，以及局部解是否朝正确方向衰减。

在数值计算中，局部加密网格理路网格与网格剖分Grid and mesh discretization用节点、单元及其几何映射离散计算区域，并以局部尺度、形状正则性和边界拟合质量描述一族网格。需使层内单元足以分辨这一尺度。固定参数时，式(1)是光滑函数；但当 $\epsilon$ 缩小时，导数界迅速增长，不能把固定函数的谱逼近常数当成参数一致的常数。谱配置法理路谱配置法Spectral collocation method在全局插值空间中求导，并用内部配点和边界方程组成微分方程的离散系统。因此也可能需要映射或分域。