本文深入研究气泡动力学在空化模拟中的应用,参考现有空化机理研究和实验结果,对现有的空化模型进行考察。在此基础上对Singhal完全空化模型进行改进,在求解完整的Rayleigh-Plesset方程的基础上构造空化源项,以考虑空化气核发展的历史效应和非定常特性。使空化模型可以模拟剧烈变化的非定常流场,例如空泡溃灭过程等。本文应用改进模型,结合基于压力的全速算法对不同头型钝体绕流过程中的空化现象进行了数值模拟。分别计算了平头钝体,半球头钝体和锥头钝体三种情况。计算的压力系数分布与实验结果和已有研究的计算结果进行了比较,改进效果良好,验证了模型的有效性。该模型可以比较准确的捕捉空化区的位置,形态以及压力变化。同时分析了不同头型空化现象的区别,讨论了外形对空化区以及模型改进效果的影响。通过数值模拟,给出了半球头钝体和锥头钝体在一定水深重力梯度作用下的空化绕流流场,与无重力流场计算结果比较,分析了环境压力和重力对空化区的影响。得到结论认为较高的环境压力可以抑制空化的形成,而重力的作用导致空化区的不稳定。计算中捕捉到了该非定常过程中的气泡脱落和压力震荡现象,说明发展的空化模型可以应用于类似现象的模拟。同时结合界面捕捉算法,将运动界面捕捉和空化流场计算结合起来。计算了半球头钝体和锥头钝体出水流场空化区形态变化和压力分布,分析了空化泡的溃灭过程和产生的压力震荡特性。该过程体现了自由界面与空化区的相互作用,这种作用导致局部空泡发生溃灭,进而引起一系列随机的瞬时压力脉冲。计算中捕捉到的瞬时压力脉冲量级达100bar,与部分实验和计算结论进行比较,证实该结果具有一定合理性。同时也探讨了该模型的适用性和发展方向,提出结合能量方程求解,考虑温度变化的计算框架。该框架仍使用基于压力的算法,不耦合求解能量方程,而是在现有的流场求解基础上单独求解能量方程。在此基础上应用工业标准的水和蒸气热力学关系式计算温度。进而修正空化计算中与压力相关的各项参数,以此体现温度的影响。