微分方程数值解课程设（2019年）.docxVIP

微分方程数值解课程设计 （五点差分格式测试及溢油事故的扩散行为模拟预测） 2019年11月 内容： 1. 数值测试：椭圆型方程五点差分格式理论结论验证（大作业） 包括： 方法理论结果验证：断误差分析（讨论是否收敛及收敛阶，通过选取不同的步长，观测步长与误差的关系）。 简单的可视化等。 2. 应用实例: 溢油事故的扩散行为模拟预测 船舶溢油事故给海洋环境带来了严重的危害。溢油事故的发生在危害海洋生态环境的同时，也给当地的渔业及旅游业等造成了不可估量的经济损失。随着时间的移动，在外界条件的作用下，进入海洋中的溢油可能会发生风化、扩散及漂移等一系列的变化，致使处理它的难度增加。需要更加高效的溢油回收、处置和清除的技术。处理溢油的常见方法有微生物分解法、化学处理法以及物理回收法微生物。分解法主要是依赖于能够对石油进行分解的微生物；化学处理法主要是通过消油剂及特定条件下的燃烧来完成的；物理回收方法主要采用机械装置，如专业收油器等强行回收海面溢油。从环境保护和资源回收再利用的方面考虑，由于消油剂具有毒性及微生物的分解缓慢，采用机械设备回收海面溢油是最理想的处理方式。 下面我们讨论用围油栏结合收油器装置的处理窄河道溢油问题。 当紧急溢油事故发生时，我们通常采用围油栏和收油装置结合的方法对溢油污染进行紧急处理。溢油发生时我们先用围油栏对溢油进行拦截使溢油范围缩小到一定程度，再利用收油装置对汇聚起来的溢油进行收集。然而我们在使用收油器时，需要注意收油器的额定收油量。在收油的过程中，单位时间收油量应与其额定收油量相同，当比额定收油量大时， 将会导致收油装置的堵塞、失效；当比额定收油量小时，收油装置就不能达到最佳的工作状态。通过数值模拟，我们可以进一步确定围油栏内的油浓度变化，进而根据收油装置在单位时间内额定收油量来选择合适的收油速度。 当溢油发生后，围油栏会立刻围住溢油，假设河道较窄，因此可认为流动的溢油在Y方向上的浓度分布是一致的，故只需考虑X方向上的溢油变化，即该物理问题可通过一维浓度扩散方程来进行模拟（我们采取无量纲的形式，只考虑理想状态下的溢油扩散问题进行模拟）。 （1） 其中水中的对流速度为u（假设恒定）, 扩散系数为D，表示某空间某时刻油的浓度。并假设初始时刻的溢油浓度为 （2） 现采用对溢油区域进行围油栏铺设的方法对溢油污染进行紧急处理，试模拟以下几种情形： 假设围油栏装置不动，在其上、下游处铺设收油装置，当溢油扩散到收油装置处时，收油装置会立刻吸收溢油，即上、下游收油装置处溢油始终为零，用数学式描述（边值为零）： （3a） 即讨论模型： 假设围油栏装置不动，在其上游处铺设围油栏，可把围油栏边界用数学表达式描述为；而在下游处铺设收油装置，当溢油在水流及围油栏的作用下向下游运动到收油装置处时，收油装置会立刻吸收溢油，即下游收油装置处溢油始终为零，用数学式描述为 。此时边界条件为： (3b） 即讨论模型： （选做）在上游处铺设围油栏并且围油栏从上游向下游对溢油进行拦截，围油栏不断移动且速度为（比如取），可把围油栏边界用数学表达式描述为,表示围油栏移动的距离；下游处铺设收油装置，当溢油在水流及围油栏的作用下向下游运动到收油装置处时，收油装置会立刻吸收溢油，即下游收油装置处溢油始终为零，用数学式描述为。此时边界条件为： (3c） 即讨论模型： 试分析（可取，D=1）： 数值模拟不同装置下，不同水流速度（水流速度分别为时）的影响对溢油问题的处理； 试着用不同的数值方法分析问题（1）； （3） 比较两种（三种）装置的处理效果。 可具体试着分析以下等问题：当围油栏由上游向下游进行围油作业时，随着围油栏的移动，围油栏边界处的溢油量变化规律？同一水流速度下，围油栏运动到不同位置时油浓度值变化规律？水流速度不同时溢油浓度值的变化规律？ 课程设计须知： 严禁抄袭！若发现雷同，不区分抄与被抄，一并作成绩不及格处理。 一、格式要求具体见附件(课程设计报告模块)，基本框架包括： （1）题目 （2）摘要：（150字左右） 摘要是对论文（设计报告）内容不加注释和评论的简短陈述，要求扼要说明研究工作的目的、主要方法、研究结果、结论、科学意义或应用价值等，是一篇具有独立性和完整性的短文。摘要中不宜使用公式、图表以及非公知公用的符号和术语， 不标注引用文献编号。 （3）关键词：四个（按重要性、频率排序，分号隔开） （4）引言：（背景、问题的提出等） （5）正文主体 正文基本内容包括（可