本文采用大涡模拟方法和小球运动方程及内置边界条件等方法，进行了三维湍流风沙跃移运动和三维沙波纹演化过程的数值模拟，主要研究工作和结果如下： 我们研制的三维拟谱法与内置边界方法和风沙运动模型相结合的大涡模拟算法和程序适合于模拟三维风沙跃移运动和三维沙波纹的形成。通过对单颗沙粒的跃移轨迹和单宽输沙率的数值模拟，验证了程序的可靠性。 在研究单颗沙粒的跃移运动中，使用单向模型对单颗沙粒在不同受力情况下的运动轨迹进行分析后，发现重力和空气阻力是沙粒所受到的最重要的两个作用力。此外，静电力对沙粒运动轨迹的影响要明显大于空气升力，Magnus力的影响主要由沙粒角速度决定，当角速度很小时，其影响可以忽略，当角速度很大时，其影响非常显著。并且对单向和双向模型得到的单颗沙粒跃移轨迹进行比较后，发现在双向耦合模型情况下沙粒跃移的距离近一些，这是由于沙粒对风有了阻碍作用。 在沙粒群体跃移运动中，静电力对沙粒的平均速度以及风沙流结构都有明显影响，空气升力则影响微小；沙粒沿Y轴的运动速度不易受到Z轴方向受力的影响，即使这个力很大。计算得到风沙流结构分为三层：线性增加层；饱和层；当高度超过饱和层时，输沙量很快下降的下降层。这与实验结果一致。并且静电力对饱和层的高度有明显影响。 在湍流风沙跃移运动和蠕移运动的共同作用下得到三维沙波纹的演化过程，最终收敛的沙波纹迎风坡平缓，背风坡陡峭，这与自然形成的沙波纹的形态是一致的。在数值模拟中还再现了三维沙波纹的自我修复过程，这从另一个方面说明了我们的数值模拟方法和所采用的跃移、蠕移模型等是合理的，可以反映沙波纹演化过程中的典型自组织特征。分析了沙波纹对三维湍流风场的影响，发现沙波纹地形的出现和发展使得湍流风场在法向有着更强烈的能量交换，揭示了湍流起沙的机理。沙波纹与地表湍流脉动涡量间有着明显的对应关系，这有助于我们解释沙波纹的形成机理。