论文首先对所使用的数值方法进行了论述：简要介绍了高精度分块耦合算法，此算法是针对不可压缩粘性流体的，采用任意曲线坐标系下原始变量形式的N-S方程进行描述。对于运动边界问题，本文引入了任意拉各朗日-欧拉法(ALE)进行处理。差分离散过程中，粘性项采用四阶精度的Pade紧致格式，对流项采用三阶迎风紧致格式离散，时间上采用四阶显式Runge-Kutta格式离散。作者所做的贡献是结合分块耦合方法的天然粗粒度并行特性，利用MPI函数库实现了并行，并且将基于分块耦合算法子域划分的任务分配策略与基于负载平衡的任务分配策略进行了比较，发现虽然分块耦合算法能够提供天然的任务分配方案，但却很难保证各进程上的负载平衡。另一个计算方法是沉浸边界法(Immersed Boundary Method)，该方法是基于笛卡儿网格系统的，因此在网格生成上可节省大量时间，而且可以大大提高计算效率，尤其适合含运动物体的流场求解。本文先归纳了前人在沉浸边界法上所做的努力，随后发展了自己的求解器，先后在交错网格和非交错网格上实现了沉浸边界法，并结合大涡模拟(LES)对高雷诺数下的湍流数值模拟进行了尝试。用高精度分块耦合方法结合ALE对弹性支撑方柱绕流进行了数值模拟。模型简化为“弹簧-阻尼-质量”系统，柱体运动采用经典龙格-库塔方法求解。同时模拟了方柱的驰振和涡致振动，捕捉到了驰振到涡致振动的转变，捕获了涡致振动的“锁定”、“拍”、“相位开关”等现象。对单圆柱绕流尾流三维特性进行了研究，找到了尾流向三维转捩的临界雷诺数Re=195，低于这个临界雷诺数时，尾流保持二维状态，此时的流场数值模拟可以近似为二维，但当超过这个临界雷诺数时，二维求解不再适用。当雷诺数超过临界值后，随着雷诺数的进一步提高，尾流先后经过两种过渡的三维模态，这两种模态分别被称作模态A(Mode A)和模态B(Mode B)。这两种转捩区的模态有各自不同的特征和物理机制，文中分别进行了研究。在此基础上，还研究了圆柱的旋转振荡对尾涡的影响。此研究是在Re=260下进行的，这个雷诺数下的圆柱尾流处在Mode B的状态。研究发现，当旋转振荡的振幅足够大时，能将尾流带回二维状态。对串列双圆柱绕流三维特性进行了研究。从两个不同方面进行了分析：首先是固定雷诺数Re=220而改变上下游圆柱之间的间距(1.5≤L／D≤8)，研究表明，当L／D≤3.5时，尾流呈现二维状态，而当L／D≥4时，尾流向三维转捩。据此可以判断临界间距存在于区间3.5≤L／D≤4；其次，将间距固定在L／D=3.5，来研究此间距下出现尾流三维特性的临界雷诺数，通过计算得到临界雷诺数为Re=250。在这部分的研究中，首次发现：在一个圆柱尾流场放置另一圆柱可以一定程度上抑制尾流向三维的转捩，从而推迟向湍流的转捩。对置于方形管道里的垂直交叉双圆柱结构进行了研究，本算例采用高精度分块耦合方法的并行程序，首先给出了详细的分块处理程序，以及内部拟边界、角线和角点的处理。圆柱端部紧贴管道壁面，两圆柱间距为圆柱直径的5倍，计算雷诺数为200。研究发现，交叉结构的中心流场较复杂，上游圆柱尾流的规则卡门涡街被下游圆柱截断成复杂的缠绕结构。两圆柱间的相互影响，以及管壁的影响，形成了流场特殊的涡结构。对无限长垂直交叉双圆柱绕流进行了研究，包括两圆柱都静止的情况以及下游圆柱强迫振动或涡致振动的情况。对于两圆柱都静止的情况，当L／D＞3时，下游圆柱对上游圆柱尾流的影响只限定在下游圆柱的尾流所扫过的范围之内；当L／D＜3时，下游圆柱对上游圆柱尾流的影响扩大，下游圆柱尾流扫过区上下出现两排三维流向二次涡结构。当下游圆柱做强迫振动或涡致振动时，尾流的涡结构和宏观量都发生了改变，尤其当下游圆柱弹性固定作涡致振动时，与独立的圆柱绕流相比，振幅响应的幅度降低，共振区域的带宽增加。最后一章提出了用计算流体力学的方法来研究鱼群游动的流体动力学性质，对一些实验根本无法实现的内容进行了探索性的研究，例如：巡游状态下，特定队形对游动推进效率的影响。从前人观察得到的菱形编队中抽出三条鱼组成的基本单元，当一条鱼位于另外两条鱼中间的下游区域时，如果处于加速巡游状态，则上游鱼尾流产生的反卡门涡街使下游鱼体两侧流体产生了向前的运动，从而提高了下游鱼游推进效率，减少了功率消耗。本文首次对前人关于鱼群减阻和提效的假设进行了验证。