稳压器波动管中热分层现象的数值模拟.pdfVIP

稳压器波动管中热分层现象的数值模拟 姚朝晖何枫孙喜明 (北京清华大学工程力学系，100084) 摘要：稳压器波动管内的热分层现象的研究是反应堆安全分析中的一个重要环节，对于这类有较大温差的低雷 诺数流动采用何种湍流模式进行数值模拟是当前反应堆热工流体力学中的～研究热点．本文通过两个热分屡流 动的典型算例考核了若干种湍流模式，在验证了Realizablek～占湍流模式计算热分层现象优越性的基础上。 采用该模式对稳压器波动管中的热分层现象进行了数值模拟，探讨了减少热分层现象的可能性，为反应堆安全 分析奠定了理论基础。 关键词；热分层现象；Realizablek～占湍流模式：安全分析 0引言 管线中的热分层现象会导致较大的热应力，从而使管线发生疲劳破坏和未预期的运动，因 No．88．08＆88．I 此美国核管会于1988年两次颁发公告(NRC。BulletinI)要求对核电站中有 温差存在的管线系统进行热分层研究，计算热应力，进而对管线的安全设计作出更合理的评估。 于在稳压器和压水堆(PWR)的热腿之间存在较大的温差，因此在连结稳压器和热腿的波动管 中有较大的热分层现象出现，所以对稳压器与热腿之间的波动管进行热分层现象的研究和计算 是核安全分析中所必不可少的。 热分层现象存在的必要条件是系统中存在温差和低的流速，温差导致介质的密度发生变 化，当管线中有带温差的介质低速注入时，热的介质由于热浮力的影响，将向上运动占据管线 的上方，这时两层流体之间被一薄的边界层分开，尽管在边界层上有热交换，但两层流体的总 体温度基本保持不变，形成熟分层现象。对于这种有较大温差的低雷诺数流动采用何种湍流模 式进行数值模拟是当前反应堆热工流体力学中的～研究重点口】。 本文探讨了k～s湍流模型、考虑热浮力k一占的湍流模式、各向异性湍流模式(ATM)、 11方程雷诺应力模式(RSM)和考虑热浮力的Realizablek～占湍流模式，通过两个热分层流 动的典型算例验证，雷诺应力模式和考虑热浮力的Realizable k～s湍流模式在一定程度可以 模拟热分层现象。但雷诺应力模式由于计算量很大，因此在应用上有一定的局限性。本文采用 考虑热浮力的Realizablek～s湍流模式，对波动管运行过程中温差较大，流速较小的危险工 况进行数值模拟，再现了热分层现象，得到了管线中温度场和速度场的分布，为热应力的计算 奠定了基础。 1基本方程和湍流模式 1．1基本方程 波动管中流动介质遵循如下的瞬态时均方程(f，，=1,2，3)啪 连续方程： 塑+一O-旦-(Ujp)：oot ％． ’ 动量方程 篙竽+毒c以㈨一言+毒c∥，等一p瓦，+偌，， cz， 焓方程： ㈥ i。h 型+堂竽=毒(寺差一俨，币)+一01)Ot Ox 0x Ox at城 J J、cp J。p 其中，U，，U，(f，J=1,2，3)是时均速度的分量：P为流体的时均压力：／2，为层流粘性 系数；P为流体密度；U，为速度U的扰动量；妒为温度r的扰动量；g，为重力加速度在3 个坐标方向上的分量；旯为热传导系数；Cp为比热；Sh为烩的源项：五万为雷诺应力：币 为湍流通量，通过模化或直接求解云石和u—j#可以建立不同的湍流模型。 izable 1．2考虑热浮力的ReaI k一￡湍流模式… Realizable k～占湍流模式中的“Realizable”意味着该湍流模式满足正应力的数学限制， 即与湍流的物理本质相符合． 在Realizable k～占湍流模式中湍动能七和湍动能耗散占的输运方程分别为： p瓦Ok+蚂考=只崛一胪+毒[(∥，+箦，考]