偏微分方程数值解试卷及答案.docxVIP

偏微分方程数值解试卷及答案

考试时间：______分钟总分：______分姓名：______

一、

简述椭圆型、抛物型、双曲型偏微分方程的典型特征（至少包括定义、阶数、典型方程形式、守恒性等方面），并举例说明。

二、

推导一维热传导方程$u_t=a^2u_{xx}$在区间$[0,L]$上，采用中心差分格式（步长分别为$\Deltax$和$\Deltat$）的离散格式。请给出差分格式，并指出其截断误差阶。

三、

对于二维拉普拉斯方程$\Deltau=0$，在矩形区域$[0,L]\times[0,H]$上，边界条件为$u|_{x=0}=u|_{x=L}=0$，$u|_{y=0}=f(x)$，$u|_{y=H}=g(x)$。请推导其基于五点差分格式的数值解法，并说明如何处理非齐次边界条件$u|_{y=0}=f(x)$。

四、

说明什么是CFL条件（Courant-Friedrichs-Lewycondition），并解释其在求解一维波动方程$u_{tt}=a^2u_{xx}$的显式差分格式（例如，采用中心差分）中的意义。若方程中波速$a=10$，空间步长$\Deltax=0.1$，为了保证数值格式的稳定性，时间步长$\Deltat$应满足什么条件？

五、

有限差分法求解偏微分方程时，如何定义并估计其数值解的误差？（请从截断误差和离散误差两个方面进行说明）。

六、

有限元法求解稳态问题（如泊松方程$\Deltau=f$）时，基于变分原理，其基本步骤通常包括哪些？请简述之。

七、

考虑在单位正方形区域$\Omega=[0,1]\times[0,1]$上求解拉普拉斯方程$\Deltau=0$，假设在边界$x=0$处$u=0$，在$y=0$处$u=0$，在$x=1$处$u=1$，在$y=1$处$u=2$。请使用有限差分法（采用合适的网格，例如4x4网格），给出该问题的数值格式（即系数矩阵和右端向量），并写出求解该线性方程组的矩阵形式。

八、

比较有限差分法、有限元法和有限体积法在求解偏微分方程时的主要异同点。在哪些方面它们各有优劣？

九、

对于一维波动方程$u_{tt}=a^2u_{xx}$，其定解问题为：$u(x,0)=\phi(x)$，$u_t(x,0)=\psi(x)$，$u(0,t)=u(L,t)=0$。请推导其基于差分格式的显式求解公式（即Yee格式或类似格式），并简述求解步骤。

试卷答案

一、

*解析思路：要求区分三种典型PDE类型。首先回顾定义：椭圆型方程对任何非零测试函数，其积分形式$\int_\OmegaL(u)v\,dx\equiv0$对所有满足$v|_{\partial\Omega}=0$的$v$都成立；抛物型方程含有时间导数项，通常与能量估计相关，如热方程；双曲型方程通常描述波传播，其特征线决定了解的结构。然后分别列举典型方程形式（如$-u_{xx}+u_{yy}=f$为椭圆，$u_t-u_{xx}=0$为双曲，$u_t=u_{xx}$为抛物），说明阶数（通常二阶），并强调其物理或数学上的主要特征，如椭圆型方程的调和函数性质、抛物型的扩散/演化性质、双曲型的波前传播性质。守恒性是抛物型和双曲型方程的重要特征，椭圆型通常不直接关联守恒律。举例时需给出具体方程并说明其类型。

二、

*解析思路：推导热传导方程的显式中心差分格式。首先写出连续方程$u_t=a^2u_{xx}$。然后在网格点$(x_i,t_n)$处，使用中心差分近似时间导数：$u_t(x_i,t_n)\approx\frac{u(x_i,t_{n+1})-u(x_i,t_n)}{\Deltat}$。使用中心差分近似空间二阶导数：$u_{xx}(x_i,t_n)\approx\frac{u(x_{i+1},t_n)-2u(x_i,t_n)+u(x_{i-1},t_n)}{(\Deltax)^2}$。将这两个近似代入原方程，得到：

$\frac{u(x_i,t_{n+1})-u(x_i,t_n)}{\Deltat}=a^2\frac{u(x_{i+1},t_n)-2u(x_i,t_n)+u(x_{i-1},t_n)}{(\Deltax)^2}$

整理得到：

$u(x_i,t_{n+1})=u(x_i,t_n)+\frac{a^2\Deltat}{(\Deltax)^2}[u(x_{i+1},