CONTAIN-LMR程序中池式钠火事故分析计算模型的验证

CONTAIN-LMR是针对以液态钠为冷却剂的反应堆而开发的安全壳事故一体化分析程序。我国目前的CONTAIN-LMR程序版本为2000年左右从法国引进,还未进行过面向工程设计的系统性地程序开发和验证。本文主要针对CONTAIN-LMR程序中模拟池式钠火事故的分析模型进行详细分析,并采用国际上的池式钠火实验进行验证,实验验证结果表明CONTAIN-LMR程序可以较准确地模拟池式钠火事故造成的钠工艺间内的温度、压力升高及放射性钠气溶胶行为。本文的研究结果初步表明CONTAIN-LMR程序可用于钠冷快堆的钠火事故分析。     

钠冷快堆混合钠火程序开发

钠管道泄漏继而发生钠的燃烧为钠冷快堆特有的事故。在喷雾钠火模型和池式钠火模型基础上,将钠喷雾燃烧和池式燃烧进行了耦合,并针对钠冷快堆钠工艺间的结构特点,最终开发了混合钠火计算程序COMSFIRE。使用该程序计算了FAUNA喷雾钠火试验和CADARACHE池式钠火试验,并与试验结果和部分程序计算结果进行了对比。同时设计了混合燃烧算例,并使用该程序与CONTAIN-LMR程序进行了对比。通过计算结果的对比和分析,初步验证了程序的正确性。     

钠冷快堆中池式钠火的计算分析

文章论述了根据池式钠火的特点建立了理论模型 ,编制了SPOOL程序。该程序模拟钠燃烧过程中钠和氧气的化学反应 ,钠燃烧热在各种介质中不同方式的传递 ,钠气溶胶的产生、沉积 ,以及在各种通风条件下多种介质的质量和能量交换等瞬态过程 ,描述了钠燃烧过程中各种特征参数随时间的变化。其主要的计算参数包括房间内气体的压力和温度、房间建筑结构的温度、钠气溶胶质量浓度等等。用俄罗斯别洛雅尔斯克核电站实验和法国卡桑德拉 3号实验的数据 ,对SPOOL程序进行验证的结果表明 ,该程序的计算结果可信。该程序为国内钠冷快堆中池式钠火事故的安全分析提供了分析方法     

快堆严重事故钠燃烧过程裂变产物释放模拟试验研究

钠燃烧过程产生的裂变产物及钠气溶胶迁移是快堆严重事故重要的源项之一。本研究对钠燃烧过程裂变产物随钠蒸汽和钠气溶胶迁移的行为进行分析,针对钠蒸发作用下裂变产物释放、钠燃烧作用下裂变产物释放以及气相空间气溶胶迁移行为分别提出了物理模型,并在确定计算方法的基础上通过CFD软件建模进行了仿真计算,最后通过开展小规模钠燃烧试验,获取了真实钠燃烧过程裂变产物沉降数据,对计算模型进行了修正和补充。试验数据与仿真计算结果表明,气溶胶迁移模型能够较好地表征裂变产物及钠气溶胶迁移行为,钠燃烧作用下裂变产物的释放系数为10-3时计算结果与试验结果较吻合。     

钠冷快增殖堆池式钠火事故分析计算

针对钠冷快堆严重事故下可能发生的池式钠火事故 ,描述了钠火现象 ,采用池式钠火程序SOFIREⅡ的“一腔室”模型 ,并在该模型基础上 ,更客观地模拟池式钠火过程 ,编制程序POOLFIRE。SOFIREⅡ认为池式钠火过程一开始就生成Na2 O和Na2 O2 的混合物 ,生成量也不确定。POOLFIRE认为早期生成Na2 O ,大约 30分钟后向Na2 O2 转变 ,最后计算得出池式钠火引起的安全壳内温度及压力响应 ,并与SOFIREⅡ结果进行了比较     

池式钠冷快堆系统分析程序开发

针对池式钠冷快堆的特点,在对快堆系统的水力模型、热工模型和中子动力学模型进行详细分类和建模的基础上,利用FORTRAN95语言开发了可用于池式钠冷快堆事故分析的系统分析程序（FASYS程序）。以中国实验快堆为计算对象对FASYS程序模型进行了初步验证,所获得的结果和试验值与其他系统程序计算值符合良好,证明了所开发的系统分析程序的正确性。     

8 CEFR放射性钠火风险分析

对CEFR的研究发现，一回路外的钠净化管道或阀门破裂将导致放射性钠火事故发生，其中309／1设备间的泄漏后果最严重。采用概率安全评价技术对该事故进行放射性风险分析，依据钠泄漏及钠火发生后相关系统的响应特性，确立了44个事故序列，建立了事件树．故障树模型。对模型进行定量分析后获得了该事件的厂内、外放射性风险，并通过对计算结果的分析，找出风险支配性因素，相应提出降低风险的措施。在整个分析过程中应用了概率安全分析软件、池式钠火分析程序及高架点源气溶胶扩散程序，获得了事故序列发生频率、钠气溶胶排放量、各种情况下高架烟囱周围3公里范围内的剂量分布等定量分析指标。     

基于FDS程序的钠火事故情景模拟技术研究

钠火是钠冷快堆的典型事故,钠火事故情景计算机模拟仿真是对钠火事故风险评价的有力工具。本文以常规火灾三维计算流体力学软件FDS为平台,增加钠火燃烧模型,包括燃烧热模型、燃烧速率模型、喷射液钠粒径分布模型等,完成了钠火情景模型的开发,并通过与SPHINCS钠火试验和计算结果的温度分布与氧气含量对比,验证了模拟技术和模型开发方案的可行性。本文的研究成果能为后续钠火仿真模拟程序的开发提供研究基础和经验参考。     

雾状钠火钠滴燃烧比率的模拟及应用

将雾状钠火中钠滴的燃烧分成预燃阶段和燃烧阶段,利用雾状钠火程序计算得到钠滴燃烧比率和时间的关系曲线,分别用幂函数、指数函数和线性函数对曲线进行拟合,拟合效果较好。拟合函数中包含钠滴下落时间和钠滴最大燃烧比率等参数,这些参数可通过钠滴下落燃烧试验或雾状钠火程序计算得到。通过推导得到了雾状钠火燃烧和单个钠滴燃烧的关系,钠滴燃烧比率的拟合函数被用来模拟雾状钠火燃烧的过程,包括用于计算已燃烧的钠质量、空气中未燃烧的钠质量、进入钠池的钠质量和雾状钠火的燃烧速率。当雾状钠火燃烧过程中钠泄漏流量恒定不变时,空气中未燃烧的钠质量和钠泄漏流量呈正比,雾状钠火的燃烧速率和钠泄漏流量呈正比。雾状钠火的燃烧速率和钠火造成的事故工艺间内的温度与压力变化直接相关。雾状钠火的燃烧速率被用来求解钠气溶胶的生成速率、钠燃烧火焰层和空气之间的传热、钠燃烧火焰层和墙壁之间的传热。总之,使用简单的函数模拟钠滴的燃烧比率曲线,将雾状钠火燃烧当成事故工艺间的热源和钠气溶胶源作为输入,便可模拟雾状钠火的整个燃烧过程,计算得到工艺间温度、压力和钠气溶胶浓度的变化。钠滴的燃烧比率曲线、雾状钠火的燃烧速率曲线还可与试验数据进行对比验证后作为雾状钠火模拟的输入,这种模拟方法可用于钠火事故安全分析中雾状钠火的模拟。     

钠冷快堆钠雾火事故三维数值模拟

针对钠冷快堆可能发生的钠雾火事故,开发和研制了一套用于钠火事故分析,能求解发生钠雾火事故后,事故空间3维的气体温度分布和化学成份分布的程序。本文详细介绍了分析求解钠雾火的燃烧模型和热传输模型,以及计算空间流场采用的计算方法和步骤,并将计算结果与试验值进行了比较：结果表明,二者符合较好。     

Response Surface Modeling of Aerosol Release Fraction in Sodium Pool Combustion

A response surface model has been proposed to evaluate an aerosol release fraction during sodium pool fire in a liquid metal fast reactor (LMFR). Air containing aerosols are radiative and they influence the allocation of combustion heat from the flame to atmospheric gas or sodium pool. Hence, the aerosol release fraction needs to be quantified based on the behavior of the aerosols and physics of mass and heat transfer. However, the aerosol release fraction is one of user-specified parameters of computer codes for the sodium fire safety analysis of the LMFR. In the present study, a response surface model of the aerosol release fraction has been developed based on numerical experiments of aerosol dynamics. For developing the model, aerosol dynamic equation has been solved coupled with thermal-hydraulics and chemical reaction. The authors obtained good agreement of the aerosol release fraction between the numerical experiments and the past experiments. Therefore, the aerosol behavior model has been validated with regard to the pool combustion phenomena and is reasonably applicable to the numerical experiment. Three influential variables on the release fraction are identified as pool temperature, gas temperature and oxygen molar fraction in the air. The proposed response surface model is a quadratic expression of the influential variables and can be easily employed in the sodium fire analysis code.     

Sodium Fire Accident Scenario Simulation Technology Study Based on FDS Code

Sodium fire is a typical accident of sodium-cooled fast reactor. Simulation of sodium fire accident scenario by software is a powerful tool for risk assessment of sodium fire accident. In this paper, the conventional fire three-dimensional computational fluid dynamics software FDS was used as a platform to add a sodium fire combustion model, including combustion heat model, combustion rate model, spray liquid sodium particle size distribution model, etc., and complete the development of sodium fire scenario modeling analysis program. And through the comparison with SPHINCS sodium fire test and calculation results, the feasibility of the method and development plan was verified. The research results of this paper can provide the research basis and experience reference for the development of the subsequent sodium fire simulation program.     

Sodium pool combustion phenomena under natural convection airflow

Sodium pool fire is a design basis accident of sodium-cooled fast reactor. In this study, a numerical method for multi-dimensional modeling of sodium pool fire has been developed. It considers coupling of thermal-hydraulics, chemical reaction and aerosol dynamics equations. From the present multi-dimensional computation, phenomena of sodium pool fire are understood such as flow and temperature fields and aerosol mass distribution of various sizes. It has been found that the burning rate varies along the radial direction and the mass and heat transfer around the pool peripheral is maximum and most influential. The thermal-hydraulic phenomena in the near-surface region are very important to determine the sodium pool fire consequence such as the burning rate and aerosol emission. The area-averaged burning rate and aerosol release fraction calculated by the present numerical method are in agreement with experimental data.     

Sodium Fire Code SFIRE1C for Pool Fire Characteristics

Sodium pool fire code, SOFIRE II, written for the constant value of stoichiometric combustion ratio and heat of reaction is used to compute the buildup of pressure and temperature in a containment. In the SOFIRE II model, for the formation of a mixture of Na₂O and Na₂O₂ in the sodium pool, the input stoichiometric combustion ratio and heat of formation values need to be varied to corresponding values admissible for the mixture. In the present work, the SOFIRE II one-cell model is revised and the present version SFIRE1C (Sodium FIRE 1 Cell model) accounts for the formation of Na₂O in an early stage of the fire and shifts to the formation of Na₂O₂ at a later stage. Thus SFIRE1C computes in a more realistic manner the reaction products which are formed in the pool. The model for sodium oxide aerosol release is also modified in this version, by incorporating a more appropriate aerosol release rate equation. The calculated values using the SFIRE1C one-cell model are compared with sodium pool fire experimental results.     

A numerical study of radiation heat transfer in sodium pool combustion and response surface modeling of luminous flame emissivity

A response surface model of the luminous flame emissivity of sodium pool fire has been proposed for use in safety analysis computer codes of a liquid metal fast reactor. The liquid sodium burns in air resulting in not only heat generation but also release of sodium oxide aerosols of sub-micron diameters. Aerosols levitating in air are radiative and they influence the allocation of combustion heat from the flame to atmospheric gas or sodium pool. The emissivity of the flame needs to be quantified, as it is one of user-specified parameters of the computer codes for the sodium fire analysis. The response surface model of the flame emissivity is developed based on numerical experiments on the physics of mass and heat transfer and behavior of the aerosol. Thermal-hydraulic equations have been solved coupled with aerosol dynamics and chemical reaction. Three influential variables on the emissivity are identified as pool temperature, gas temperature and oxygen molar fraction in the air. It has been found that the emissivity is calculated reasonably as a function of the three variables. The proposed response surface model can be easily employed in the sodium fire analysis codes because it is a simple quadratic expression. For the safety evaluation of the sodium fire, combined use is recommended of the proposed model and the lumped-mass zone model code.     

Validation study of computer code sphincs for sodium fire safety evaluation of fast reactor

A computer code sphincs solves coupled phenomena of thermal hydraulics and sodium fire based on a multi-zone model. It deals with an arbitrary number of rooms, each of which is connected mutually by doorways and penetrations. With regard to the combustion phenomena, a flame sheet model and a liquid droplet combustion model are used for pool and spray fires, respectively, with the chemical equilibrium model based on the Gibbs free energy minimization method. The chemical reaction and mass and heat transfer are solved interactively. A specific feature of sphincs is detailed representation of thermalhydraulics of a sodium pool and a steel liner, which is placed on the floor to prevent sodium-concrete contact. The authors analyzed a series of pool combustion experiments, in which gas and liner temperatures are measured in detail. It has been found that good agreement is obtained and the sphincs code has been validated with regard to pool combustion phenomena. Further research needs are identified for pool spreading modeling considering thermal deformation of steel liner and measurement of pool fluidity property as a mixture of liquid sodium and reaction products. The sphincs code is to be used mainly in the safety evaluation of the consequence of a sodium fire accident in a liquid metal cooled fast reactor as well as fire safety analysis in general.     

钠冷快堆钠池火事故数值模拟

为了估计和预测钠火事故的后果，构建了以“有火焰薄层”为理论基础的燃烧模型和热传输模型，给出了程序计算结果与试验值的比较。比较结果证实，该计算结果可信、模型合理。程序可用来分析和预测钠池火事故。