数值分析课程设计学生题目.docVIP

《数值分析》课程设计 本课程设计的内容为：每个小组的同学均应完成以下五个案例； 目标：能将数值分析课程中所学的算法知识熟练应用于实际问题中。 土木工程和环境工程师在设计一条排水渠道时必须考虑渠道的各种参数（如宽度，深度，渠道内壁光滑度）及水流速度、流量、水深等物理量之间的关系。 假设修一条横断面为矩形的水渠，其宽度为B，假定水流是定常的，也就是说水流速度不随时间而变化。 根据质量守恒定律可以得到 Q=UBH （1.1） 其中Q 是水的流量（），U是流速（），H是水的深度（）。 在水工学中应用的有关流速的公式是 （1.2） 这里n是Manning粗糙系数，它是一个与水渠内壁材料的光滑性有关的无量纲量；S是水渠的斜度系数，也是一个无量纲量，它代表水渠底每米内的落差。 把（1.2）代入（1.1）就得到 （1.3） 为了不同的工业目的（比如说要把污染物稀释到一定的浓度以下，或者为某工厂输入一定量的水），需要指定流量Q和B，求出水的深度。这样，就需要求解 （1.4） 一个具体的案例是 求出渠道中水的深度H。 所涉及的知识——非线性方程解法。 在化学工程中常常研究在一个封闭系统中同时进行的两种可逆反应 其中A，B，C和D代表不同的物质。反应达到平衡是有如下的平衡关系： 其中称为平衡常数，代表平衡状态时该物质的浓度。假定反应开始时各种物质的浓度为： 而且反应达到平衡时，由第一和第二种反应生成的C物质浓度分别为，于是平衡时满足的方程为： 用不同的数值方法求解上述方程。 所涉及的知识——非线性方程组解法。 湖水在夏天会出现分层现象，接近湖面温度较高，越往下温度变低。这种上热下冷的现象影响了水的对流和混合过程，使得下层水域缺氧，导致水生鱼类的死亡。如果把水温看成深度的函数T(x)，有某个湖的观测数据如下： T (°C) 22.8 22.8 22.8 20.6 13.9 11.7 11.1 11.1 x (m) 0 2.3 4.9 9.1 13.7 18.3 22.9 27.2 环境工程师希望： 用样条插值求出 T(x). 求在什么深度处达到最大（即） 所涉及的知识——插值、数值微分。 在排污管道设计中，工程师关心管道坡度、管子直径和污水流量之间的关系。对于圆截面管道这些量之间有如下经验公式： 其中Q代表流量（），S代表管道坡度（m/m），D代表圆管直径（m），是三个通过实验测定的经验参数。 有一组实验数据如下： 实验序号 S 1 2 3 4 5 6 7 8 9 0.302 0.604 0.906 0.302 0.604 0.902 0.302 0.604 0.906 0.001 0.001 0.001 0.01 0.01 0.01 0.05 0.05 0.05 0.0385 0.2283 0.6655 0.1293 0.7948 2.3100 0.3053 1.8975 5.5000 用适当的数值方法求出 所涉及的知识——最小二乘拟合。 在研究建筑物通过地板散失热量时，我们需要计算建筑物下方地基中的温度变化。假设建筑物是圆形的，其半径r=2m（如图所示） 假定： 室内温度恒定保持在25°C。 室外离开建筑物2m以外（即R≥4m）地基温度不受室内温度影响。 地层4m一下温度保持为10°C。 室外地表温度随昼夜温度变化而变化，其变化规律为 （12.1） 时间单位为小时。 我们再假设，地基是由均匀的黄土组成，其物性参数是 密度 导热系数 比热 我们要研究的是半径R=4m，高度H=4m的一块柱形地基中的温度变化问题。 由于几何上的对称性，我们可以沿对称轴做一个垂直剖面，并建立坐标系（如图所示） 在地基内P温度应当满足柱坐标下的热传导方程 （12.2） 其中 边界条件为 （12.3） （12.4） 根据地层传热学中的傅立叶定律可以得知 因为随着时间的流逝，在开始一瞬间对后来温度变化的影响逐渐消失，所以可以任意假设，不妨设其为10°C。 值得指出的是我们需要知道的是足够长的时间之后（t≥T），24小时地基温度的变化和由建筑物内P流失到地层中的热量。 所涉及的知识——数值微分，线性方程组求解，数值积分。 数值分析课程中一类算法的可视化设计 包括插值或数值微分，积分或解线性方程组或非线性方程的求根或常微分方程的数值解法  4m 4m y z z