CN115577583B - 一种铅基反应堆绕丝定位燃料棒流致振动的分析方法 - Google Patents
一种铅基反应堆绕丝定位燃料棒流致振动的分析方法 Download PDFInfo
- Publication number
- CN115577583B CN115577583B CN202211095591.6A CN202211095591A CN115577583B CN 115577583 B CN115577583 B CN 115577583B CN 202211095591 A CN202211095591 A CN 202211095591A CN 115577583 B CN115577583 B CN 115577583B
- Authority
- CN
- China
- Prior art keywords
- fluid
- fuel rod
- calculation
- domain
- lead
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/02—Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E30/00—Energy generation of nuclear origin
- Y02E30/30—Nuclear fission reactors
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Computing Systems (AREA)
- Fluid Mechanics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Monitoring And Testing Of Nuclear Reactors (AREA)
Abstract
一种铅基反应堆绕丝燃料棒流致振动分析方法,根据绕丝定位方式的燃料棒束组件进行建模,分别针对流体域以及固体域划分网格,并设置流固耦合交界面;流体域部分计算前调用铅铋合金物性关系式,在流固耦合交界面的轴向不同位置上设置监测点,监测不同位置的流体激振力随时间的变化;采用大涡模拟模型,并选择WALE模型作为亚格子应力模型,进行瞬态计算,直至收敛。固体域部分设置固定端面,将流体域压力场加载到流固耦合交界面上,直至固体域每个时间步的有限元计算达到收敛;对流体激振力以及振动位移进行频域功率谱密度函数分析。本发明为铅基反应堆绕丝定位燃料棒束的流致振动以及振动导致的微动磨损问题提供了评估与指导。
Description
技术领域
本发明涉及反应堆流致振动分析技术领域,具体涉及一种铅基反应堆堆芯绕丝定位燃料棒流致振动问题的分析方法。
背景技术
铅基反应堆是以液态金属铅或者铅铋合金(统称铅基材料)为冷却剂的一类反应堆。铅基快堆堆芯设计紧凑,燃料组件长期处于高辐照、高服役温度的严苛条件,冷却剂动态冲击下棒束振动造成疲劳损伤的同时,导致氧化层磨蚀剥落从而加剧腐蚀,长换料周期需求凸显了上述问题,验证影响铅基反应堆的设计与安全可靠运行。因此评估铅基快堆高密度冷却剂高速冲击下稠密燃料棒束的响应机制显得尤为重要。冷却剂冲刷下燃料组件的稳定性和完整性研究是铅基快堆的重要技术瓶颈,其可靠性是决定铅冷快堆设计与安全特性的关键因素之一。
压水堆中定位格架与燃料棒之间的振动磨损占压水堆燃料破损泄漏比例高达55%。铅基快堆燃料组件主要采用金属螺旋绕丝定位维持棒束径向间隙以紧凑堆芯并提高换热效率,液态铅基材料密度大且流速较快,堆芯稠密布置燃料棒束中存在类似振动失效现象。随着燃耗加深,芯块肿胀和包壳热膨胀使绕丝与相邻燃料棒间隙缩小甚至发生接触,高密度冷却剂冲刷燃料棒束产生显著附加作用力而诱发振动;定位绕丝扰动边界层造成强烈湍流搅混,燃料组件流场轴向以及周向各向异性更显著,从而加剧了振动。高频率高幅度振动响应导致燃料元件疲劳损伤,加速燃料元件包壳的腐蚀破损。
由于铅基合金冷却反应堆是四代反应堆堆型之一,各国对其研究还尚处探索阶段,对于铅基堆中绕丝定位燃料棒的流致振动的研究较少。考虑到液态铅基合金的不透明的物理特性,铅基堆中的流致振动主要是以数值模拟为主,分析方法主要集中在单根绕丝燃料棒的流体域与固体域的单向流固耦合,对于多根燃料棒形成的绕丝定位组件的分析还很缺乏。
发明内容
为了克服上述现有技术存在的问题,本发明的目的在于提供一种铅基反应堆绕丝定位燃料棒流致振动的分析方法,该方法针对铅基反应堆绕丝定位燃料棒的特点,通过构建反应堆绕丝定位燃料棒运行环境和设置铅基合金工质物性,开展瞬态热工水力计算,再将计算获得的流体激振力加载到燃料棒表面,通过有限元方法开展振动响应计算,并在时域和频域两方面分析流体激振力的波动和燃料棒的振动位移响应。
为了实现上述目的,本发明采取了以下技术方案予以实施:
一种铅基反应堆绕丝定位燃料组件流致振动的分析方法,步骤如下:
步骤1:绕丝棒束组件几何模型建模:根据绕丝定位方式的燃料棒组件进行建模,包括流体域和固体域两个部分,获得带绕丝棒束组件流动通道三维几何,以及绕丝棒束组件结构三维几何;分别对棒束组件流动通道三维几何和绕丝棒束组件结构三维几何两个部分进行网格划分,将燃料棒与绕丝表面设置为流固耦合交界面;
步骤2:绕丝棒束组件流动通道CFD计算:采用计算流体动力学方法进行流体域计算,获得瞬态计算时长内每个时间步的压力场以及速度场,具体分为以下步骤:
步骤2-1:输入燃料棒束通道入口初始速度和初始压力作为稳态计算初始值,并采用RANS模型作为稳态计算的湍流模型;
步骤2-2:调用铅铋合金物性关系式,由于不考虑温度场的影响,因此铅铋合金物性关系式主要包括铅铋合金的密度和动力粘度:
ρLBE=11096-1.3236T (1)
μLBE=(4.56-7.03×10-3T+3.61×10-6T2)×10-3 (2)
式中,
ρLBE——铅铋合金密度/kg·m-3
T——铅铋合金温度/K
μLBE——铅铋合金动力粘度/N·s·m-2
计算时温度取铅基反应堆正常工况下的温度500K进行稳态计算,直至稳态计算达到收敛,获得燃料棒束通道的稳态压力场以及速度场,作为瞬态计算的初始值;
步骤2-3:在流固耦合交界面的轴向不同位置设置监测点,输出不同位置处每个时间步的流体激振力的值;
步骤2-4:流体域瞬态计算时的湍流模型采用大涡模拟模型,过滤方式采取盒式过滤,亚格子应力模型选择WALE亚格子应力模型,进行瞬态计算,直至流体瞬态计算残差达到0.001的收敛条件;
步骤3:流固耦合交界面压力和剪切应力分布分析:将每个时间步的流体激振力的值通过周期图法进行频域分析,获得流体激振力的功率谱密度函数的结果,并且输出每个时刻燃料棒表面每个轴向位置的压力以及剪切力的值,分析压力和剪切力对绕丝定位燃料棒的影响;
步骤4:绕丝棒束组件流致振动时域分析:将流体域计算获得的流体随时间变化的压力场加载到流固耦合交界面上,并将燃料组件轴向两端固定,分别选取棒束的中间位置燃料棒以及周围其中一根燃料棒进行分析计算,直至每个时间步的有限元计算达到收敛,获得每个时间步的不同位置的燃料棒表面不同点的振动位移;
步骤5:绕丝棒束组件流致振动频域分析:将每个时间步的不同位置的燃料棒表面不同点的振动位移通过周期图法进行频域分析,获得振动位移的功率谱密度函数的结果;
步骤6:绕丝棒束组件固有频率分析:通过预应力模态分析,获得绕丝棒束组件各阶固有频率。并将步骤5的频域分析结果与固有频率结果相比较,分析不同位置燃料棒的振动与固有频率的关系。
和现有技术相比较,本发明具有如下优点:
1、本发明的分析方法能够对绕丝燃料棒束组件的流致振动进行分析,为组件中燃料棒之间的振动磨损计算提供了基础;
2、本发明流体域计算采用的是大涡模拟模型,能够更好地捕捉到湍流脉动特性,对于流体激振力的模拟更加精准;
3、本发明的分析方法通过对振动位移响应进行时域和频域两方面的分析,可以获得任意位置的燃料棒以及燃料棒上不同轴向高度的振动响应,增加了该方法的通用性。
附图说明
图1为本发明的计算流程框图。
图2为某种铅基反应堆绕丝定位燃料组件结构示意图。
具体实施方式
下面结合附图和具体实施方式对本发明方法进一步详细说明:
如图1所示,本发明一种铅基反应堆绕丝定位燃料棒流致振动的分析方法,步骤如下:
步骤1:绕丝棒束组件几何模型建模:根据绕丝定位方式的燃料棒组件进行建模,图2所示的模型包括流体域和固体域两个部分,分别为带绕丝棒束组件流动通道三维几何,以及绕丝棒束组件结构三维几何,分别对流动通道三维几何以及绕丝棒束组件结构三维几何两个部分进行网格划分,流体域部分采用多面体与棱柱层混合网格的形式进行网格划分,而固体域部分采用四面体网格的形式进行网格划分。并将燃料棒与绕丝表面设置为流固耦合交界面;
步骤2:绕丝棒束组件流动通道CFD计算:采用计算流体动力学方法进行流体域计算,获得瞬态计算时长内每个时间步的压力场以及速度场,具体分为以下步骤:
步骤2-1:输入燃料棒束通道入口初始速度和初始压力作为稳态计算初始值,并采用RANS模型作为稳态计算的湍流模型;
步骤2-2:调用铅铋合金物性关系式,由于不考虑温度场的影响,因此铅铋合金物性关系式主要包括铅铋合金的密度和动力粘度:
ρLBE=11096-1.3236T (1)
μLBE=(4.56-7.03×10-3T+3.61×10-6T2)×10-3 (2)
式中,
ρLBE——铅铋合金密度/kg·m-3
T——铅铋合金温度/K
μLBE——铅铋合金动力粘度/N·s·m-2
计算时温度取铅基反应堆正常工况下的温度500K进行稳态计算,直至稳态计算达到收敛,获得燃料棒束通道的稳态压力场以及速度场,作为瞬态计算的初始值;
步骤2-3:在流固耦合交界面的轴向不同位置设置监测点,监测点的位置具体设置在图2所示模型的中心位置和周围位置的燃料棒的轴向不同高度的点上,输出每根燃料棒的不同位置处每个时间步的流体激振力的值;
步骤2-4:流体域瞬态计算时的湍流模型采用大涡模拟模型,过滤方式采取盒式过滤,亚格子应力模型选择WALE亚格子应力模型,进行瞬态计算,直至流体瞬态计算残差达到0.001的收敛条件;
步骤3:流固耦合交界面压力和剪切应力分布分析。将每个时间步的流体激振力的值通过周期图法进行频域分析,获得流体激振力的功率谱密度函数的结果,并且输出每个时刻燃料棒表面每个轴向位置的压力以及剪切力的值,分析压力和剪切力对绕丝定位燃料棒的影响;
步骤4:绕丝棒束组件流致振动时域分析。将流体域计算得到的流体域随时间变化的压力场加载到流固耦合交界面上,并将燃料组件轴向两端固定,分别选取棒束的中间位置燃料棒以及周围其中一根燃料棒进行分析计算,直至每个时间步的有限元计算达到收敛,获得每个时间步的不同位置的燃料棒表面不同点的振动位移;
步骤5:绕丝棒束组件流致振动频域分析。将每个时间步的不同位置的燃料棒表面不同点的振动位移通过周期图法进行频域分析,获得振动位移的功率谱密度函数的结果。
步骤6:绕丝棒束组件固有频率分析。通过预应力模态分析,获得绕丝棒束组件各阶固有频率。并将步骤5的频域分析结果与固有频率结果相比较,分析不同位置燃料棒的振动与固有频率的关系。
Claims (4)
1.一种铅基反应堆绕丝定位燃料棒流致振动的分析方法,其特征在于:包括以下步骤:
步骤1:绕丝棒束组件几何模型建模:根据绕丝定位方式的燃料棒组件进行建模,包括流体域和固体域两个部分,获得带绕丝棒束组件流动通道三维几何,以及绕丝棒束组件结构三维几何;分别对绕丝棒束组件流动通道三维几何和绕丝棒束组件结构三维几何两个部分进行网格划分,将燃料棒与绕丝表面设置为流固耦合交界面;
步骤2:绕丝棒束组件流动通道CFD计算:采用计算流体动力学方法进行流体域计算,获得瞬态计算时长内每个时间步的压力场以及速度场,具体分为以下步骤:
步骤2-1:输入燃料棒束通道入口初始速度和初始压力作为稳态计算初始值,并采用RANS模型作为稳态计算的湍流模型;
步骤2-2:调用铅铋合金物性关系式,由于不考虑温度场的影响,因此铅铋合金物性关系式包括铅铋合金的密度和动力粘度:
ρLBE=11096-1.3236T (1)
μLBE=(4.56-7.03×10-3T+3.61×10-6T2)×10-3 (2)
式中,
ρLBE——铅铋合金密度/kg·m-3
T——铅铋合金温度/K
μLBE——铅铋合金动力粘度/N·s·m-2
计算时温度取铅基反应堆正常工况下的温度500K进行稳态计算,直至稳态计算达到收敛,获得燃料棒束通道的稳态压力场以及速度场,作为瞬态计算的初始值;
步骤2-3:在流固耦合交界面的轴向不同位置设置监测点,输出不同位置处每个时间步的流体激振力的值;
步骤2-4:流体域瞬态计算时的湍流模型采用大涡模拟模型,过滤方式采取盒式过滤,亚格子应力模型选择WALE亚格子应力模型,进行瞬态计算,直至流体瞬态计算残差达到0.001的收敛条件;
步骤3:流固耦合交界面压力和剪切应力分布分析:将每个时间步的流体激振力的值通过周期图法进行频域分析,获得流体激振力的功率谱密度函数的结果,并且输出每个时刻燃料棒表面每个轴向位置的压力以及剪切力的值,分析压力和剪切力对绕丝定位燃料棒的影响;
步骤4:绕丝棒束组件流致振动时域分析:将流体域计算获得的流体随时间变化的压力场加载到流固耦合交界面上,并将燃料组件轴向两端固定,分别选取棒束的中间位置燃料棒以及周围其中一根燃料棒进行分析计算,直至每个时间步的有限元计算达到收敛,获得每个时间步的不同位置的燃料棒表面不同点的振动位移;
步骤5:绕丝棒束组件流致振动频域分析:将每个时间步的不同位置的燃料棒表面不同点的振动位移通过周期图法进行频域分析,获得振动位移的功率谱密度函数的结果;
步骤6:绕丝棒束组件固有频率分析:通过预应力模态分析,获得绕丝棒束组件各阶固有频率;并将步骤5的频域分析结果与固有频率结果相比较,分析不同位置燃料棒的振动与固有频率的关系。
2.如权利要求1所述的一种铅基反应堆绕丝定位燃料棒流致振动的分析方法,其特征在于:步骤1中,所获得的带绕丝棒束组件流动通道三维几何为包含绕丝几何特征的完整三维几何模型。
3.如权利要求1所述的一种铅基反应堆绕丝定位燃料棒流致振动的分析方法,其特征在于:步骤1中,流体域采用六面体和棱柱层混合网格的方法进行网格划分,固体域采用四面体网格的方法进行网格划分。
4.如权利要求1所述的一种铅基反应堆绕丝定位燃料棒流致振动的分析方法,其特征在于:步骤4中,固体域的有限元计算的时间步长要与流体域的流体动力学计算输出保存的时间步长相同,保证每个时间步加载一次流体激振力。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211095591.6A CN115577583B (zh) | 2022-09-06 | 2022-09-06 | 一种铅基反应堆绕丝定位燃料棒流致振动的分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211095591.6A CN115577583B (zh) | 2022-09-06 | 2022-09-06 | 一种铅基反应堆绕丝定位燃料棒流致振动的分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115577583A CN115577583A (zh) | 2023-01-06 |
CN115577583B true CN115577583B (zh) | 2023-07-04 |
Family
ID=84581302
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211095591.6A Active CN115577583B (zh) | 2022-09-06 | 2022-09-06 | 一种铅基反应堆绕丝定位燃料棒流致振动的分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115577583B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112699620A (zh) * | 2021-01-08 | 2021-04-23 | 西安交通大学 | 基于计算流体力学的反应堆堆芯热工水力特性分析方法 |
CN113782105A (zh) * | 2021-08-31 | 2021-12-10 | 西安交通大学 | 一种摇摆条件下液态金属铅铋合金流动传热特性分析方法 |
CN114997081A (zh) * | 2022-05-31 | 2022-09-02 | 中山大学 | 一种大型压水堆堆芯结构流致振动的数字反应堆仿真分析方法 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110020479B (zh) * | 2019-04-09 | 2022-07-08 | 中国核动力研究设计院 | 一种圆筒结构随机湍流激励诱发振动的分析方法 |
EP3896703A1 (en) * | 2020-04-14 | 2021-10-20 | Westinghouse Electric Sweden AB | Computer-based simulation methods for boiling water reactors (bwr) |
CN111859752B (zh) * | 2020-07-18 | 2021-11-16 | 西安交通大学 | 一种核反应堆蒸汽发生器流致振动计算方法 |
CN113642132B (zh) * | 2021-08-25 | 2023-07-14 | 中国核动力研究设计院 | 一种基于cfd的管束结构流致振动评估临界流速分析方法 |
CN114036871B (zh) * | 2021-11-26 | 2023-10-20 | 中国核动力研究设计院 | 基于瞬态分析的堆内棒状结构湍流激振分析方法及装置 |
CN114282460B (zh) * | 2021-12-31 | 2023-05-16 | 西安交通大学 | 一种铅铋快堆堆芯热工水力特性分析方法 |
-
2022
- 2022-09-06 CN CN202211095591.6A patent/CN115577583B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112699620A (zh) * | 2021-01-08 | 2021-04-23 | 西安交通大学 | 基于计算流体力学的反应堆堆芯热工水力特性分析方法 |
CN113782105A (zh) * | 2021-08-31 | 2021-12-10 | 西安交通大学 | 一种摇摆条件下液态金属铅铋合金流动传热特性分析方法 |
CN114997081A (zh) * | 2022-05-31 | 2022-09-02 | 中山大学 | 一种大型压水堆堆芯结构流致振动的数字反应堆仿真分析方法 |
Also Published As
Publication number | Publication date |
---|---|
CN115577583A (zh) | 2023-01-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8149983B2 (en) | Method for analysis of pellet-cladding interaction | |
de Stordeur | Drag coefficients for fuel-element spacers | |
CN115577583B (zh) | 一种铅基反应堆绕丝定位燃料棒流致振动的分析方法 | |
Massih et al. | Analyses of pellet-cladding mechanical interaction behaviour of different ABB Atom fuel rod designs | |
Fisher et al. | Measurement of dynamic interaction between a vibrating fuel element and its support | |
Lee et al. | Preliminary Multi-Physics Analysis of a 2x2 Rod Array Using CUPID/GIFT Coupled Code | |
Page et al. | The Performance of Zirconium Alloy Clad UO2 Fuel for Canadian Pressurized and Boiling Water Power Reactors... | |
Torres et al. | RMB experimental program on the hydrodynamical behavior of fuel assemblies | |
Li et al. | Thermal hydraulic analysis of a deformed SFR fuel assembly using CFD | |
Revankar et al. | Assessment of steam generator tube flaw size and leak rate models | |
JP4323901B2 (ja) | 沸騰水型原子炉のボイド率 | |
Shin et al. | Effect of Rod Internal Pressure on Simulation of Halden Test IFA-650.9 with FE-based Fuel Analysis Code MERCURY | |
Mazurkiewicz et al. | Meeting the industry challenge to achieve leak free performance | |
Uddin et al. | The Huge Missing Factor in Leak-Before-Break Analysis: How a Circumferential Through-Wall-Crack in a Pipe System Changes the Flexibility and Reduces the Applied Moments | |
Merckx | DISTORTION OF HELICAL WIRE WRAPPINGS | |
Pickman | Interactions between fuel pins and assembly components | |
Kim | Plan for Structural Analysis of Fuel Assembly for Seismic and Loss of Coolant Accident Loading Considering End-Of-Life Condition for APR1400 NRC Design Certification | |
Jo et al. | Life estimation of PWR steam generator U-tubes subjected to foreign object-induced fretting wear | |
Scuro et al. | A CFD analysis of blockage length on a partially blocked fuel rod | |
Du et al. | Benchmark Study of RBHT Experiment for the Effect of Spacer Grids on Reflood Heat Transfer With LOCUST Code | |
Ibrahim et al. | Turbulence prediction in two-dimensional bundle flows using large eddy simulation | |
Han et al. | Flow-induced vibration analysis of steam generators and fuel assemblies with the VIBIC computer code | |
Demick | Using SA508/533 for the HTGR Vessel Material | |
Mukhopadhyay et al. | A methodology for on-line fatigue life monitoring of Indian nuclear power plant components | |
Scott | PRTR MECHANICAL SEAL PUMP OPERATING EXPERIENCE FOR SEPTEMBER 1959-AUGUST 1960 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |