海洋环境数值计算2013-1课件.pptVIP

海洋环境数值计算 海洋环境数值计算(1) 课程设置目的 在学习掌握数值计算相关知识的基础上，通过介绍和实际练习几个海洋动力学基本应用案例，掌握海洋环境数值计算中的基本方法，为今后应用数值计算方法研究海洋科学与工程技术问题打下基础。 一、前言 1、数值计算是科学研究的第三手段 第一手段——实验分析，形成理论 第二手段——理论分析 第三手段——数值计算或称科学计算 数值计算是伴随人类计算能力的增强而发展起来的新研究手段，计算机的高速发展使大规模科学计算有可能实现。 数值计算与实验分析、理论分析既相对独立，又互相补充、互相依赖。 优势： 模拟现实世界中的各种过程； 替代实验或作为实验的有效补充； 检验理论的正确性； 再现现实生活中无法重复或无法进行实验的现象或耗资巨大的实验工程。 2、数值计算促进交叉学科的形成 数值计算作为一种强有力的研究手段，被广泛应用于科学与工程技术领域，并发展成一系列的交叉学科分支：计算数学、计算物理学、计算化学、计算力学等等。 3、数值计算用于预测和工程设计的条件 对问题的物理本质了解得很清楚，并可用精确的数学模型来表达； 有足够计算能力的计算机，以便在实际可能的时间和费用下，用数值方法解出这些问题。 4、发展回顾 天体力学：轨道计算 空气动力学：航天、航空 天气动力学：数值天气预报 海洋动力学：海洋潮汐、环流、海气耦合、污染预测、生态系统等。 5、海洋环境数值计算 研究海洋中各种动力学过程、现象，进而建立相应的数值模型，最终实现对海洋中各种动力过程、现象的预测。 6、需要的基础 物理基础——需要理解动力过程的物理原理 数学基础——需要理解以偏微分方程为主的物理过程的数学描述和数值解法 编程基础——需要通过编程实现对数值模型的求解，建议使用FORTRAN语言 7、实验平台 建议采用WINDOWS XP 平台下的Compaq Visual FORTRAN 6.6，带编辑器和调试模块，比较方便。 或WINDOWS 7平台下的Intel Visual Fortran Compiler Professional(11.0.061) +MS Visual Studio 2008. 建议：不熟悉FORTRAN语言的同学，希望能在比较短的时间里，学习掌握编程语言 8、参考书： 数值计算类（电子版）： 海洋数值计算方法和数值模式 流体力学数值模拟 编程类的（电子版）： Fortran77和90/95编程入门 数学类的：偏微分方程数值解 解法类的：数值分析 二、数值计算的基本知识 （一）数值计算的主要方法 一个物理问题的数值解。方法有许多，从数学上表示可分为： 有限差分 有限元 有限基本解 有限分析 等几种基本解法，在这个基础上，学者们根据不同解题背景及需求提出了很多具体的解法。 在海洋研究中应用较多的有： 1、有限差分法 ——将微分方程转化为差分方程求解 优点：简单易用，计算量小 缺点：边界条件处理不是很理想；当方程中存在对流项时，计算结果会出现数值效应（频散或伪振荡）；Re数大时，空间步长被逼成很小，△x/L～1/Re 。 使用广泛。 2、有限元法 ——由试函数法演变而来，用一系列线性无关的基函数试解求基函数的系数 优点：非线性边界条件被包含在计算本身当中，不需另外列出。 缺点：计算量大，一般都需要逐时间步求解大型线性方程组；对流项存在时也会产生伪振荡；不可压缩条件（流体）较难满足。 主要变形方法有：加权余量法（配置法、子区域法、最小二乘法、矩量法、伽略金法Galerkin），变分法等。 五十年代用于飞机结构设计，六十年代中开始算流体，使用也较广泛。 3、谱方法 ——应当是有限差分法的一种变形 优点：将求解变量展开成三角级数（正交），接近于精确解。 缺点：只适合于规则网格 八十年代用于三维潮汐动力计算 4、有限分析法 ——利用局部分析解来建立数值离散格式 优点：可在较大雷诺数时得到精度较高的解 缺点：繁杂，工作量大 求解域＝》小子域＝》分析解 被用于不可压粘性流计算 5、有限基本解 ——绕流边界划为单元，每个单元布置基本解 优点：使问题减少一维 缺点：弯曲大的物体区域精度欠缺 用于解位势流 还有分步法、分布杂交法、多重网格法、流速分解法等，多重网格法结合有限差分可解区域关注度不同的问题。 国际上各类海洋模型 POM - Princeton Ocean Model MOM - GFDL Modular Ocean Model ROMS (Regional Ocean Modeling System) HYCOM (HYbrid Coordi