明尼苏达大学Fortran实现的SIMPLE算法

"该资源是一个关于Simple算法的FORTRAN语言编写的程序代码，源自于明尼苏达大学的研究生课程，经过改编以适应个人计算机环境，并修正了一些不适当的表达。这个程序仅供教学使用，禁止出版，并可以作为研究目的的代码框架进行二次开发。" 在FORTRAN编程中，Simple算法通常用于解决流体力学中的数值模拟问题，特别是求解 Navier-Stokes 方程的简化版本。这个程序可能是一个基础的流体动力学模拟软件，它包含了以下几个关键部分： 1. **主程序** (MAINPROGRAM): 这是程序的入口点，负责初始化、调用各个子程序以及控制程序流程。在示例代码中，`LOGICALLSTOP` 是一个逻辑变量，用于判断是否停止程序执行。`COMMON/CNTL/LSTOP` 声明了一个公共变量 `LSTOP`，它在多个子程序间共享，用于控制循环。 2. **文件操作** (`OPEN` 和 `CLOSE`): `OPEN(08,FILE='teresul')` 打开了名为 'teresul' 的文件，可能是用于存储计算结果。在程序结束或特定条件下，通过 `CLOSE(08)` 关闭文件。 3. **子程序调用**: - `SETUP0`, `SETUP1`, `SETUP2`: 这些可能是设置初始条件、参数或网格的子程序。`SETUP` 通常用于准备计算前的必要工作。 - `GRID`: 可能是定义网格结构的子程序，网格化是流体模拟中的重要步骤，用于将连续域离散化为可计算的网格点。 - `START`: 可能是初始化计算过程的子程序，如设置初值、边界条件等。 - `DENSE`: 这个名字暗示了它可能涉及求解系统的稠密矩阵，即进行迭代求解流体方程的过程。 - `BOUND`: 处理边界条件，这是数值方法中的关键部分，确保流体在物理边界上的行为符合实际要求。 - `OUTPUT`: 可能是输出计算结果到文件的子程序，供后续分析或可视化。 4. **循环与控制**: 主程序中的 `IF(.NOT.LSTOP) GOTO 15` 结构形成了一个循环，只要 `LSTOP` 为 `.NOT.` (非真)，就会继续执行循环，否则跳出并执行后续的关闭文件和停止程序的命令。 5. **注释**: 注释中提到，此程序由西安交通大学的讲师在1995年11月提供，表明这是一份有一定历史背景的教学材料。 这个FORTRAN程序展示了如何用编程实现Simple算法来解决流体动力学问题，涉及到网格生成、初始化、迭代求解、边界处理和结果输出等多个步骤。对于学习FORTRAN编程和流体模拟的人来说，这是一个宝贵的实践案例。