2025年大学《数理基础科学》专业题库- 偏微分方程的数值解法

2025年大学《数理基础科学》专业题库——偏微分方程的数值解法考试时间：______分钟总分：______分姓名：______一、选择题（请将正确选项的代表字母填在题后括号内。每小题2分，共10分）1.对于一维热传导方程∂u/∂t=α∂²u/∂x²，以下差分格式中，具有CFL条件的是（）。A.显式格式：u_i^(n+1)=u_i^n+r(u_(i+1)^n-2u_i^n+u_(i-1)^n)，其中r=αΔt/Δx²B.隐式格式：u_i^(n+1)=u_i^n+r(u_(i+1)^n-2u_i^n+u_(i-1)^n)，其中r=αΔt/Δx²C.Crank-Nicolson格式：u_i^(n+1)=u_i^n+r/2[u_(i+1)^n-2u_i^n+u_(i-1)^n+u_(i+1)^(n+1)-2u_(i+1)^(n+1)+u_(i-1)^(n+1)]D.向后差分格式：u_i^(n+1)=u_i^n-r(u_(i+1)^n-2u_i^n+u_(i-1)^n)，其中r=αΔt/Δx²2.在求解二维拉普拉斯方程∇²u=0时，若采用中心差分格式，其离散形式为(u_(i+1,j)+u_(i-1,j)+u_(i,j+1)+u_(i,j-1)-4u_(i,j))/Δx²=0，则该格式是（）。A.一阶精度B.二阶精度C.三阶精度D.四阶精度3.对于有限元法，下列说法正确的是（）。A.空间离散是在整个求解区域上直接进行。B.形状函数必须满足分片常数条件。C.单元刚度矩阵的组装是基于物理方程在整个区域上的积分。D.有限元方法天然适用于求解非结构化网格问题。4.在求解线性方程组Ax=b，其中A是稀疏矩阵，下列迭代法中，通常收敛速度较快且适用于稀疏矩阵的是（）。A.Jacobi方法B.Gauss-Seidel方法C.SOR方法D.共轭梯度法5.下列哪种偏微分方程属于双曲型方程？（）A.Laplace方程∇²u=0B.热传导方程∂u/∂t=α∂²u/∂x²C.波动方程∂²u/∂t²=c²∂²u/∂x²D.扩散方程∂u/∂t=D∂²u/∂x²二、填空题（请将答案填在题后横线上。每小题3分，共15分）6.有限差分法中，为了使时间离散格式稳定，对于显式时间积分格式，其时间步长Δt与空间步长Δx的比值需要满足一定的条件，称为______条件。7.有限元法中，将求解区域划分为有限个单元是为了实现______，从而将复杂的区域问题简化为易于处理的小区域问题。8.对于二维泊松方程∇²u=f，在边界点i上，若采用中心差分格式离散内部点，则点i处的方程可以写为u_(i,j)=(u_(i+1,j)+u_(i-1,j)+u_(i,j+1)+u_(i,j-1)-Δx²f_i)/4。这种处理边界点的方法称为______。9.在求解线性方程组Ax=b时，若采用Jacobi迭代法，其迭代公式为x^(k+1)=Bx^(k)+g，其中B是由A的对角元构成的矩阵，则Jacobi迭代法收敛的必要条件是矩阵B的______范数小于1。10.有限体积法的基本思想是利用______对控制体积上的积分形式PDE进行离散，得到的离散方程天然满足守恒性。三、计算题（请按题目要求进行计算和推导。每小题10分，共30分）11.考虑一维热传导方程∂u/∂t=α∂²u/∂x²在区间[0,L]上，边界条件为u(0,t)=u(L,t)=0，初始条件为u(x,0)=f(x)。现采用向前差分格式对时间进行离散，中心差分格式对空间进行离散，步长分别为Δt和Δx。推导出节点i处的离散方程（即差分格式），并写出相应的迭代关系。12.推导二维拉普拉斯方程∇²u=0在笛卡尔坐标系下的中心差分格式。假设网格节点(i,j)的函数值为u_ij，空间步长为Δx=Δy。要求推导出u_ij的表达式，并指出其精度阶。13.考虑一个简单的有限元问题：在区域Ω=[0,1]上求解u''(x)=-1，边界条件为u(0)=0,u(1)=0。使用线性有限元方法，将区域划分为两个单元[0,0.5]和[0.5,1]，试推导出该问题的单元形状函数，并组装全局刚度矩阵和载荷向量（假设使用自然坐标系）。四、综合题（请结合所学知识进行分析和解答。每小题15分，共30分）14.对于二维波动方程∂²u/∂t²=c²(∂²u/∂x²+∂²u/∂y²)，边界条件为u|∂Ω=0，初始条件为u(x,y,0)=φ(x,y),∂u/∂t(x,y,0)=ψ(x,y)。(1)写出该方程的显式中心差分格式。(2)分析该格式的时间步长Δt需要满足的稳定性条件（CFL条件）。(3)与隐式格式相比，简述显式中心差分格式的优缺点。15.比较有限差分法（FDM）和有限元法（FEM）在求解偏微分方程时的主要异同点。(1)在离散思想上有何不同？(2)在处理边界条件和复杂几何区域方面有何差异？(3)在求解大型线性方程组方面有何不同考虑？(4)简述各自的主要优点。试卷答案一、选择题1.A2.B3.C4.D5.C二、填空题6.CFL(Courant-Friedrichs-Lewy)7.离散化(Discretization)8.强制边界条件(Essentialboundaryconditionimposition/Directimposition)9.谱(Spectral)10.控制体积(Controlvolume)三、计算题11.解析思路：首先应用向前差分格式离散时间导数∂u/∂t≈(u_i^(n+1)-u_i^n)/Δt。其次应用中心差分格式离散空间二阶导数∂²u/∂x²≈(u_(i+1,j)-2u_i,j+u_(i-1,j))/Δx²。将初始条件和边界条件代入，整理即可得到节点i处的离散方程。最后，将时间步长Δt用α和Δx表达，即可得到迭代关系。答案：差分格式为(u_i^(n+1)-u_i^n)/Δt=α(u_(i+1,j)-2u_i,j+u_(i-1,j))/Δx²。整理得u_i^(n+1)=u_i^n+r(u_(i+1,j)-2u_i,j+u_(i-1,j))，其中r=αΔt/Δx²。迭代关系为u_i^(n+1)=u_i^n+r(f_i^n+u_(i+1,j^n)-2u_i,j^n+u_(i-1,j^n))。12.解析思路：首先写出二维拉普拉斯方程的精确表达式∂²u/∂x²+∂²u/∂y²=0。然后分别用中心差分格式近似空间一阶导数∂u/∂x≈(u_(i+1,j)-u_(i-1,j))/(2Δx)和空间二阶导数∂²u/∂x²≈(u_(i+1,j)-2u_i,j+u_(i-1,j))/Δx²。同样，近似空间一阶导数∂u/∂y≈(u_(i,j+1)-u_(i,j-1))/(2Δy)和空间二阶导数∂²u/∂y²≈(u_(i,j+1)-2u_i,j+u_(i,j-1))/Δy²。将这两个近似式代入原方程，合并同类项即可得到u_ij的表达式。分析所得表达式中关于u_ij及其邻近点值的高阶项系数，判断其精度阶。答案：近似空间导数：∂u/∂x≈(u_(i+1,j)-u_(i-1,j))/(2Δx),∂²u/∂x²≈(u_(i+1,j)-2u_i,j+u_(i-1,j))/Δx²；∂u/∂y≈(u_(i,j+1)-u_(i,j-1))/(2Δy),∂²u/∂y²≈(u_(i,j+1)-2u_i,j+u_(i,j-1))/Δy²。代入方程得(u_(i+1,j)-2u_i,j+u_(i-1,j))/Δx²+(u_(i,j+1)-2u_i,j+u_(i,j-1))/Δy²=0。假设Δx=Δy，整理得u_ij=(u_(i+1,j)+u_(i-1,j)+u_(i,j+1)+u_(i,j-1))/4。此表达式仅涉及节点i及其直接相邻的四个节点，且近似表达式中u_ij的系数为1，高阶项系数为-1/Δx²和-1/Δy²。由于空间导数是用中心差分二阶格式近似，故该离散格式整体具有二阶空间精度。13.解析思路：首先，为线性有限元方法选取基函数。对于区间[0,1]上的线性有限元，常用基函数为φ_0(x)=1-x和φ_1(x)=x。节点0和1对应的基函数值分别为(φ_0(0),φ_0(0.5),φ_0(1))=(1,0.5,0)和(φ_1(0),φ_1(0.5),φ_1(1))=(0,0.5,1)。其次，写出单元刚度矩阵和载荷向量的积分表达式。单元刚度矩阵k_e=∫(Aφ_i'φ_j')dx，其中A=-1是导数φ_i'的系数。单元载荷向量f_e=∫(φ_j*(-1))dx。第三，对每个单元进行推导。对于单元[0,0.5]，积分区间为[0,0.5]，φ_0'=-1,φ_1'=1。k_e=∫((-1)*(-1)*φ_0+(-1)*(1)*φ_1)dx=∫(φ_0-φ_1)dx=[x-x]_0^0.5=0.5-0.5=0.25。f_e=∫(-φ_j)dx=-∫φ_jdx=-[xφ_j]_0^0.5=-0.5φ_j(0.5)。由于f_e是在单元内积的，对于单元[0,0.5]，f_e=-0.5φ_1(0.5)=-0.5*0.5=-0.25。对于单元[0.5,1]，积分区间为[0.5,1]，φ_0'=-1,φ_1'=1。k_e=∫((-1)*(-1)*φ_0+(-1)*(1)*φ_1)dx=∫(φ_0-φ_1)dx=[x-x]_0.5^1=1-0.5-(0.5-0.5)=0.25。f_e=∫(-φ_j)dx=-[xφ_j]_0.5^1=-(1*φ_1(1)-0.5*φ_1(0.5))=-(1*1-0.5*0.5)=-1+0.25=-0.75。注意，f_e需要考虑单元上的源项，此处源项为-1。第四，将两个单元的刚度矩阵和载荷向量组装成全局矩阵A和向量f。组装时需考虑单元节点与全局节点的对应关系。答案：单元[0,0.5]：形状函数φ_0(x)=1-x,φ_1(x)=x。单元刚度矩阵k_e1=0.25(2x2矩阵)。单元载荷向量f_e1=-0.25(长度为2向量，对应节点0和1，但仅单元内部贡献，此题f_e1应理解为单元内部积分结果，组装时需调整)。单元[0.5,1]：形状函数φ_0(x)=1-x,φ_1(x)=x。单元刚度矩阵k_e2=0.25(2x2矩阵)。单元载荷向量f_e2=-0.75(长度为2向量，对应节点0和1，此题f_e2应理解为单元内部积分结果，组装时需调整)。组装全局刚度矩阵A=[k_e1k_e2]=[0.25000.25]。组装全局载荷向量f=[f_e1f_e2]=[-0.25-0.75]。注意：此处的组装方式假设了全局节点编号为0,1,2...，而单元节点为0,1和0,1。实际组装需明确节点编号对应关系。更准确地说，对于线性元，单元刚度矩阵应乘以单元长度（这里是0.5），单元载荷向量应等于单元长度的源项平均值（这里是-0.5）。组装后A=[0.1250;00.125]，f=[-0.5-0.5]。如果题目要求使用自然坐标系推导，则需要将基函数表示为自然坐标λ的函数（φ_0(λ)=0.5(1-λ),φ_1(λ)=0.5(1+λ)），然后进行积分推导，结果会一致。四、综合题14.解析思路：(1)首先将时间导数∂²u/∂t²用向前差分近似：∂²u/∂t²≈(u(x,y,t+Δt)-2u(x,y,t)+u(x,y,t-Δt))/(Δt²)。将空间二阶导数∂²u/∂x²和∂²u/∂y²用中心差分近似：∂²u/∂x²≈(u_(i+1,j)-2u_i,j+u_(i-1,j))/Δx²，∂²u/∂y²≈(u_(i,j+1)-2u_i,j+u_(i,j-1))/Δy²。代入原方程得到近似的差分方程。(2)对该显式格式进行稳定性分析，通常采用傅里叶方法。假设解的形式为u(x,y,t)=e^(i(kx+ly-ωt))，代入差分方程。显式格式为u_i^(n+1)=u_i^n+(cΔt/Δx)²[u_(i+1,j)^(n)-2u_i,j^(n)+u_(i-1,j)^(n)]+(cΔt/Δy)²[u_(i,j+1)^(n)-2u_i,j^(n)+u_(i,j-1)^(n)]。代入试函数，得到关于amplificationfactorA的方程A=1+r_x²+r_y²，其中r_x=cΔt/Δx,r_y=cΔt/Δy。要求A的模小于等于1，即|1+r_x²+r_y²|≤1。由于r_x²≥0,r_y²≥0，故1+r_x²+r_y²≥1。因此，要满足稳定性条件，必须|1+r_x²+r_y²|≤1，这仅在r_x²+r_y²=0时成立，即Δt=0，这是无意义的。更常见的理解是对于波动方程的显式格式，需要满足Δt≤Δx/(c√2)或Δt≤Δy/(c√2)，这通常被称为广义CFL条件，确保相位速度的稳定性。但最严格的CFL条件是Δt≤Δx/(c)。此处直接写Δt≤cΔt/(c√2Δx)或类似形式，强调时间步长与空间步长、波速的关系。(3)比较优缺点：优点是计算简单，每个节点的新值只依赖于当前时刻的节点值和相邻节点值，可以显式求解，无需求解大型线性方程组。缺点是稳定性条件严格，时间步长受限，对于高波速或细网格问题可能非常小，导致计算效率低。答案：(1)显式中心差分格式为u(x,y,t+Δt)≈2u(x,y,t)-u(x,y,t-Δt)+(cΔt/Δx)²[u(x+Δx,y,t)-2u(x,y,t)+u(x-Δx,y,t)]+(cΔt/Δy)²[u(x,y+Δy,t)-2u(x,y,t)+u(x,y-Δy,t)]。(2)稳定性条件通常要求时间步长Δt满足Δt≤Δx/(c√2)或类似形式，确保相位速度的稳定性。更严格的解释是涉及r_x=cΔt/Δx和r_y=cΔt/Δy的关系，要求r_x²+r_y²≤1。(3)优点是计算格式简单，易于编程实现，属于显式求解，不需要迭代求解大型线性方程组。缺点是稳定性条件（CFL条件）较为严格，时间步长Δt受空间步长Δx,Δy和波速c的限制，可能导致计算效率不高，尤其是在需要高分辨率（小Δx,Δy）或高波速（大c）的情况下。15.解析思路：(1)离散思想：FDM直接在网格节点上近似求解，将偏微分方程转化为代数方程组。FEM将求解区域划分为单元，在每个单元上选择基函数进行插值，将偏微分方程转化为单元上的弱形式（积分形式），然后对整个区域积分，得到一个关于未知节点值的代数方程组。FVM基于控制体积概念，通过对积分形式的PDE在控制体积上进行离散，得到满足守恒性的代数方程。(2)边界条件和复杂几何：FDM需要显式地处理边界条件，对于复杂边界或非结构网格可能较困