CN110807246B - 一种Sub栅元尺度的反应堆热工水力控制体划分方法 - Google Patents

一种Sub栅元尺度的反应堆热工水力控制体划分方法 Download PDF

Info

Publication number
CN110807246B
CN110807246B CN201910930011.2A CN201910930011A CN110807246B CN 110807246 B CN110807246 B CN 110807246B CN 201910930011 A CN201910930011 A CN 201910930011A CN 110807246 B CN110807246 B CN 110807246B
Authority
CN
China
Prior art keywords
transverse
flow
vortex
sub
control body
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
Application number
CN201910930011.2A
Other languages
English (en)
Other versions
CN110807246A (zh
Inventor
陈广亮
田兆斐
张志俭
李磊
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Harbin Engineering University
Original Assignee
Harbin Engineering University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Harbin Engineering University filed Critical Harbin Engineering University
Priority to CN201910930011.2A priority Critical patent/CN110807246B/zh
Publication of CN110807246A publication Critical patent/CN110807246A/zh
Application granted granted Critical
Publication of CN110807246B publication Critical patent/CN110807246B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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

Abstract

本发明涉及核反应堆堆芯热工水力计算领域，具体涉及一种Sub栅元尺度的反应堆热工水力控制体划分方法。建立Sub栅元尺度流域特征及规律数据库；确定Sub栅元尺度横向控制体划分的区域；确定横流旋涡流域与非旋涡横流流域；开展Sub栅元尺度横向旋涡控制体区划分；开展横向非旋涡的横流基础控制体区划分；开展非旋涡横流控制体的细致划分。本发明能够模拟堆芯单个栅元以内参量的横向分布特性，提升核反应堆堆芯热工水力子通道分析空间分辨率，更适应于堆芯组件流域的实际热工水力过程。

Description

一种Sub栅元尺度的反应堆热工水力控制体划分方法
技术领域
本发明涉及核反应堆堆芯热工水力计算领域，具体涉及一种Sub栅元尺度的反应堆热工水力控制体划分方法。
背景技术
核动力的安全性与经济性受核反应堆堆芯热工水力状态预测技术水平的影响。精细化堆芯热工水力计算流体动力学(CFD)程序的计算能够在精细空间尺度上模拟堆芯的释热、传热与流动。然而堆芯精细化热工水力CFD计算分析的仿真网格数量巨大，计算资源占用量大，计算时间长。
由于控制体内物理参量的平均化处理，虽然增大CFD程序控制体的尺寸能够一定程度上提高计算效率，然而即使考虑全面的CFD控制方程在控制体尺寸增大到一定程度后将无法模拟细致的三维TH状态；且边界层内固体对液体的黏性影响需要通过小尺寸控制体进行计算；另外，针对宏观对象进行宏观量化总结的经验计算公式对大区域宏观平均参量的预测更为适用。
目前，通过对物理方程进行简化，减少控制体数量，考虑相邻横向大控制体间的质量、动量、能量输运，开发堆芯子通道计算程序，是堆芯热工水力细致状态模拟的主流工程化计算技术。然而，该类技术的研发始于20世纪60年代，受当时计算技术与堆芯热工水力机理认知的影响，程序控制体的最小横向范围为单个栅元级流道，不能模拟栅元流道内的参量细致分布、力与运动作用等细致热工水力特性。为提升堆芯热工水力计算分析技术在工程化应用中的精度，有必要开发新型堆芯子通道分析计算的控制体划分方法。
发明内容
本发明的目的在于提供一种Sub栅元尺度的反应堆热工水力控制体划分方法，以提升核反应堆堆芯热工水力子通道分析空间分辨率。
本发明实施例提供一种Sub栅元尺度的反应堆热工水力控制体划分方法，包括：
步骤一：建立Sub栅元尺度流域特征及规律数据库：根据不同工况时、不同空间位置处，辨识所研究堆芯组件不同流域Sub栅元尺度冷却剂的横流特性及其规律，建立数据库；
步骤二：确定Sub栅元尺度横向控制体划分的区域：根据研究域中具有复杂横流运动的区域，将该区域确定为待进行Sub栅元尺度横向控制体划分的区域；
步骤三：确定横流旋涡流域与非旋涡横流流域：根据步骤二所述的待进行Sub栅元尺度横向控制体划分的区域，基于堆芯计算分析域中Sub栅元尺度冷却剂的横流特性，通过对栅元级冷却剂通道内的横向流域进行分类，在横向上特征流域划分，得到横流旋涡流域和非旋涡横流流域；
步骤四：开展Sub栅元尺度横向旋涡控制体区划分：根据步骤三所述的横流旋涡流域，通过进行横向旋涡控制体划分，将每个横流旋涡流域划分，得到1个Sub栅元尺度的横流旋涡流域控制体；
步骤五：开展横向非旋涡的横流基础控制体区划分：根据步骤三所述的非旋涡横流流域，通过进行横向非旋涡的横流基础控制体区划分，基于每个栅元的横流出入口数量和位置，将每个非旋涡横流流域划分，得到4个横向非旋涡的横流基础控制体；
步骤六：开展非旋涡横流控制体的细致划分：根据计算效率与空间分辨率要求，基于步骤一所述的横流特性与规律，对步骤五所述的横向非旋涡的横流基础控制体区进行非旋涡横流控制体的细致划分，将每个横向非旋涡的横流基础控制体区划分，得到一个或多个非旋涡横流控制体；
本发明还包括这样一些结构特征：
所述步骤一中，所述Sub栅元尺度具体为：空间分辨率能够辨识栅元内部的横向二维流动特性，即辨识Sub栅元尺度的横向二维流动特性；
所述横流特性具体为：堆芯组件不同流域Sub栅元尺度冷却剂的横向流动特性，重点是含搅混翼定位格架下游附近流域的横向流动特性；其中，上述横向指与竖直方向垂直的方向；
所述建立数据库的具体内容为：针对不同运行工况时堆芯组件不同重要流域，记录Sub栅元尺度的横流特性及规律，重点记录各栅元流道内大尺寸旋涡以及所关心旋涡所在区域位置和范围，记录各栅元流道内非旋涡横流的横向流道与流动方向；其中，上述所关心旋涡通过研究者经验、精细化实验与计算研究经验、堆芯横流运动规律与机理研究的方式确定；
所述步骤五中，所述横向非旋涡的横流基础控制体的数目确定方法为：每个栅元的横向冷却剂流动出入口的数量为4，对应各横向出入口，将非旋涡横流流域分为4个横向非旋涡的横流基础控制体；
所述横向非旋涡的横流基础控制体的特征为：每个横向非旋涡的横流基础控制体区含有1个栅元横向流动出口或入口；
本发明的有益效果在于：
1.横向上小于单燃料栅元的Sub栅元尺度控制体的建立：传统反应堆子通道程序的横向最小控制体为单个栅元或更大的区域范围，无法模拟栅元内的：精细横流运动、旋涡运动、精细压力梯度、液体间的黏性力等。
本发明在子通道仿真中采用小于单燃料栅元的最小控制体，且控制体的划分形式是基于Sub栅元尺度精细化的热工水力特性，能够模拟栅元流道中横向Sub栅元尺度的热工水力过程，包括Sub栅元尺度的旋涡状态、栅元内与栅元间的横流状态、栅元内的压力分布、栅元内的液体间黏性力等细致的堆芯热工水力过程；
2.栅元及横截面积更大区域的参量非均匀化计算：若控制体的横向面积是单栅元横截面积或更大的横截面积，由于控制体内采用参量均匀化处理，参量在单栅元范围内的横向分布无法模拟，然而实际情况是单栅元范围压力、温度、速度等参量的横向非均匀分布特性显著，所以传统的子通道程序控制体划分方案的适用性受限。
本发明所设计的Sub栅元尺度的反应堆热工水力控制体划分方法由于采用Sub栅元级别的控制体，尺寸小，能够模拟堆芯单个栅元以内参量的横向分布特性，更适应于堆芯组件流域的实际热工水力过程。
附图说明
图1为典型栅元级流道的二维横截面示意图；
图2为本发明栅元级流道的横向流域划分示意图；
图3(a)为本发明栅元级流道中的流线及栅元中心旋转旋涡示意图；
图3(b)为本发明栅元级流道中的流线及栅元中心旋转旋涡示意图；
图4为本发明栅元级流道横截面的横向二维控制体划分示意图；
图5(a)为本发明栅元级流道横截面的横流旋涡流域控制体和横流基础控制体区划分示意图；
图5(b)为本发明栅元级流道横截面的横流旋涡流域控制体和横流基础控制体区划分示意图；
图6(a)为本发明栅元级流道横截面的非旋涡横流控制体划分示意图；每个横向非旋涡的横流基础控制体区划分成2个非旋涡横流控制体；
图6(b)为本发明栅元级流道横截面的非旋涡横流控制体划分示意图；每个横向非旋涡的横流基础控制体区划分成4个非旋涡横流控制体；
图7为一种Sub栅元尺度的反应堆热工水力控制体划分方法的流程图。
标号说明：1为横流旋涡流域；2为非旋涡横流流域；3为横流旋涡流域控制体；4为横向非旋涡的横流基础控制体。
具体实施方式
为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图对本发明做进一步描述：
本发明的技术方案是这样实现的：
步骤1：建立数据库，记录不同工况下，不同空间位置处，Sub栅元尺度流域的横流特性及规律；
步骤2：确定研究域中具有复杂横流运动的区域，并将该区域确定为进行Sub栅元尺度横向控制体划分的区域，该步骤可通过研究者经验、精细化实验与计算研究经验、堆芯横流运动规律与机理等方式完成。
步骤3：对所确定的Sub栅元尺度横向控制体划分区域进行流域分类，将流域在横向上划分为旋涡流域、非旋涡横流流域。
步骤4：进行横向上的旋涡控制体划分，将每个旋涡流域划分为一个Sub栅元尺度的旋涡流域控制体。
步骤5：进行横向上非旋涡的横流基础控制体区划分，根据每个栅元的横流出入口数量和位置，将非旋涡流域划分为与横流出入口数目相同的基础控制体区，且每个基础控制体区含有一个横流出口或入口。
步骤6：根据计算效率与空间分辨率要求，根据横流特性与规律，对每个横流基础控制体区进行非旋涡横流控制体的细致划分。
本发明还可以这样描述：
(1)辨识所研究堆芯组件不同流域Sub栅元尺度冷却剂的横流特性，尤其是辨识含搅混翼定位格架下游附近流域的横流特性；
(2)基于堆芯计算分析域中Sub栅元尺度的冷却剂横向流动特性，对栅元级冷却剂通道内的横向流域进行分类，分为横向旋涡流域、非旋涡横流流域；
(3)进行横向二维控制体划分，对栅元级冷却剂通道内的每个横向旋涡流域建立1个Sub栅元尺度的控制体，对栅元级冷却剂通道内的非旋涡横流流域分为4个横流基础控制体区；
(4)根据计算效率与空间分辨率要求，根据横流特性与规律，每个横流基础控制体区划分为1个或多个非旋涡横流控制体。
本发明还可以包括：
1.所述辨识所研究堆芯组件不同流域Sub栅元尺度冷却剂横流特性具体包括：
建立数据库，针对不同运行工况，为堆芯组件不同重要流域，记录Sub栅元尺度的横流特性及规律，重点记录各栅元流道内大尺寸旋涡或所关心旋涡所在的区域位置和范围，记录各栅元流道内非旋涡横流的横向流道与流动方向。
本发明所述的Sub栅元尺度是指空间分辨率能够辨识栅元内部的横向二维流动特性，即辨识Sub栅元尺度的横向二维流动特性，图1为典型栅元尺度的流道二维横截面。
本发明所述的横向指与竖直方向垂直的方向。
本发明所述所关心旋涡可通过研究者经验、精细化实验与计算研究经验、堆芯横流运动规律与机理研究等方式确定。
2.所述对研究流域中的每个栅元级冷却剂通道内的横向流域进行分类，分为横向旋涡流域、非旋涡横流流域，具体包括：
基于数据库中不同运行工况下不同流域Sub栅元尺度的横流特性及规律，针对研究域中每个栅元级冷却剂通道，在横向上进行特征流域划分。
图3中左右两图是横向研究域流域划分的两个示例，两个示例在横向上均包括2个相邻的栅元级流道，各栅元级流道中心处各有一个大尺寸的旋涡，图3中的两示例的旋涡旋转方向及相邻栅元流道间的横流流动方向存在差别，本发明将包含旋涡的中心流域划分为旋涡流域，栅元级流道中旋涡流域外的其它流域划分为非旋涡横流流域。
如图2所示，将各栅元流道内横向大尺寸旋涡或所关心旋涡的所在区域换分成1个横向旋涡流域，或分别为每个大尺寸旋涡或所关心旋涡所在的区域各自划分为横向旋涡流域。将每个栅元级冷却剂通道内除旋涡流域外的流域，划分为非旋涡横流流域。
3.所述进行横向二维控制体划分，对栅元级冷却剂通道内的每个横向旋涡流域建立1个Sub栅元尺度的旋涡流域控制体，对栅元级冷却剂通道内的非旋涡横流流域划分为4个横流基础控制体区，具体包括：
如图4所示，在横向上，将每个旋涡流域划分成1个横流旋涡控制体；对于非旋涡流域，由于每个栅元的横向冷却剂流动出入口的数量为4，对应各横向出入口，将非旋涡流域划分为4个横流基础控制体区，保证每个横流基础控制体区含有1个栅元横向流动出口或入口。
4.所述根据计算效率与计算空间分辨率要求，每个横流基础控制体区划分为1个或多个非旋涡横流控制体，具体包括：
如图6所示，栅元中每个横流基础控制体区根据计算效率与空间分辨率要求，根据横流特性与规律，可以继续划分为多个控制体；或不进行划分，即1个横流基础控制体区划分为1个非旋涡横流控制体。图6(a)为每个横流基础控制体区划分为2个非旋涡横流控制体的示例，图6(b)为每个横流基础控制体区划分为4个非旋涡横流控制体的示例。

Claims (2)

1.一种Sub栅元尺度的反应堆热工水力控制体划分方法，其特征在于，包括：
步骤一：建立Sub栅元尺度流域特征及规律数据库：根据不同工况时、不同空间位置处，辨识所研究堆芯组件不同流域Sub栅元尺度冷却剂的横流特性及其规律，建立数据库；所述Sub栅元尺度具体为：空间分辨率能够辨识栅元内部的横向二维流动特性，即辨识Sub栅元尺度的横向二维流动特性；
所述横流特性具体为：堆芯组件不同流域Sub栅元尺度冷却剂的横向流动特性，重点是含搅混翼定位格架下游附近流域的横向流动特性；其中，上述横向指与竖直方向垂直的方向；
所述建立数据库的具体内容为：针对不同运行工况时堆芯组件不同重要流域，记录Sub栅元尺度的横流特性及规律，重点记录各栅元流道内大尺寸旋涡以及所关心旋涡所在区域位置和范围，记录各栅元流道内非旋涡横流的横向流道与流动方向；其中，上述所关心旋涡通过研究者经验、精细化实验与计算研究经验、堆芯横流运动规律与机理研究的方式确定；
步骤二：确定Sub栅元尺度横向控制体划分的区域：根据所研究的堆芯组件流域中具有复杂横流运动的区域，将该区域确定为待进行Sub栅元尺度横向控制体划分的区域；
步骤三：确定横流旋涡流域与非旋涡横流流域：根据步骤二所述的待进行Sub栅元尺度横向控制体划分的区域，基于堆芯计算分析所研究堆芯组件流域中Sub栅元尺度冷却剂的横流特性，通过对栅元级冷却剂通道内的横向流域进行分类，在横向上进行特征流域划分，得到横流旋涡流域和非旋涡横流流域；
步骤四：开展Sub栅元尺度横向旋涡控制体区划分：根据步骤三所述的横流旋涡流域，通过进行横向旋涡控制体划分，将每个横流旋涡流域划分，得到1个Sub栅元尺度的横流旋涡流域控制体；
步骤五：开展横向非旋涡的横流基础控制体区划分：根据步骤三所述的非旋涡横流流域，通过进行横向非旋涡的横流基础控制体区划分，基于每个栅元的横流出入口数量和位置，将每个非旋涡横流流域划分，得到4个横向非旋涡的横流基础控制体；
步骤六：开展非旋涡横流控制体的细致划分：根据计算效率与空间分辨率要求，基于步骤一所述的横流特性与规律，对步骤五所述的横向非旋涡的横流基础控制体区进行非旋涡横流控制体的细致划分，将每个横向非旋涡的横流基础控制体区划分，得到一个或多个非旋涡横流控制体。
2.根据权利要求1所述的一种Sub栅元尺度的反应堆热工水力控制体划分方法，其特征在于：所述步骤五中，所述横向非旋涡的横流基础控制体的数目确定方法为：每个栅元的横向冷却剂流动出入口的数量为4，对应各横向出入口，将非旋涡横流流域分为4个横向非旋涡的横流基础控制体；
所述横向非旋涡的横流基础控制体的特征为：每个横向非旋涡的横流基础控制体区含有1个栅元横向流动出口或入口。
CN201910930011.2A 2019-09-29 2019-09-29 一种Sub栅元尺度的反应堆热工水力控制体划分方法 Active CN110807246B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910930011.2A CN110807246B (zh) 2019-09-29 2019-09-29 一种Sub栅元尺度的反应堆热工水力控制体划分方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910930011.2A CN110807246B (zh) 2019-09-29 2019-09-29 一种Sub栅元尺度的反应堆热工水力控制体划分方法

Publications (2)

Publication Number Publication Date
CN110807246A CN110807246A (zh) 2020-02-18
CN110807246B true CN110807246B (zh) 2022-07-08

Family

ID=69487913

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910930011.2A Active CN110807246B (zh) 2019-09-29 2019-09-29 一种Sub栅元尺度的反应堆热工水力控制体划分方法

Country Status (1)

Country Link
CN (1) CN110807246B (zh)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007232547A (ja) * 2006-02-28 2007-09-13 Hitachi Ltd 自然循環式沸騰水型原子炉における原子炉炉心性能計算方法と原子炉炉心性能計算装置
CN106897478A (zh) * 2016-12-19 2017-06-27 杨保文 基于阻力分布包含格架搅混效应的子通道分析方法

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4056091A (en) * 1975-04-23 1977-11-01 Moncrieff Yeates Alexander Joh Vortical flow aerothermodynamic fireplace unit
US8369475B2 (en) * 2009-07-01 2013-02-05 Westinghouse Electric Company Llc Nuclear fuel assembly support grid
WO2012005465A2 (ko) * 2010-07-08 2012-01-12 에스케이텔레콤 주식회사 무선랜 전파 환경 맵을 이용한 ａｐ 위치 추정 방법 및 장치
KR101471038B1 (ko) * 2010-07-08 2014-12-10 에스케이텔레콤 주식회사 무선랜 전파 환경 맵을 이용한 ａｐ 위치 추정 방법 및 장치
CN103050153B (zh) * 2012-12-20 2016-02-17 中国核动力研究设计院 用于核燃料组件的具有防勾挂以及交混作用的定位格架
CN108304620B (zh) * 2018-01-11 2020-07-24 西安交通大学 一种核反应堆堆芯棒束通道内硼扩散过程的计算方法
CN108897907A (zh) * 2018-05-25 2018-11-27 哈尔滨工程大学 一种用于反应堆堆芯热工水力分布并行计算的方法
CN108875212A (zh) * 2018-06-19 2018-11-23 哈尔滨工程大学 核动力装置主冷却剂系统热工水力多尺度分区仿真方法
CN108875213B (zh) * 2018-06-19 2022-04-12 哈尔滨工程大学 一种反应堆堆芯热工水力多尺度分析方法
CN109903870B (zh) * 2019-03-15 2020-06-26 西安交通大学 一种核动力系统跨维度耦合模拟方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007232547A (ja) * 2006-02-28 2007-09-13 Hitachi Ltd 自然循環式沸騰水型原子炉における原子炉炉心性能計算方法と原子炉炉心性能計算装置
CN106897478A (zh) * 2016-12-19 2017-06-27 杨保文 基于阻力分布包含格架搅混效应的子通道分析方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
板状燃料反应堆热工水力实时仿真程序研究与开发;张志俭等;《核动力工程》;20101231(第06期);全文 *

Also Published As

Publication number Publication date
CN110807246A (zh) 2020-02-18

Similar Documents

Publication Publication Date Title
Li et al. Large-eddy simulation of flow and pollutant dispersion in high-aspect-ratio urban street canyons with wall model
Liu et al. Population balance modeling of polydispersed bubbly flow in continuous-casting using multiple-size-group approach
Li et al. Modeling of mixing behavior in a combined blowing steelmaking converter with a filter-based Euler–Lagrange model
Zhang et al. CFD-PBM simulation of bubble columns: Effect of parameters in the class method for solving PBEs
Zhao et al. Numerical investigation of transverse flow and turbulent mixing in a helical cruciform fuel bundle
Zhang et al. Structural characteristics of a spiral symmetry stream anaerobic bioreactor based on CFD
Vahidifar et al. Numerical simulation of particle-laden flow in an industrial wastewater sedimentation tank
CN110807246B (zh) 一种Sub栅元尺度的反应堆热工水力控制体划分方法
Hassan et al. Computational investigation of magneto-hydrodynamic flow of newtonian fluid behavior over obstacles placed in rectangular cavity
Zhao et al. High throughput millifluidics mixer with double triangle baffle for improvement of mixing performance and reduction of flow resistance designed by grey relational analysis
Zhang et al. Numerical analysis on subcooled boiling in PWR coolant channel based on a modified multi-scale interface model
Wang et al. Enhanced particle focusing and sorting by multiple sheath stream in contraction–expansion microchannel
Hashemi et al. Simulation of turbulent flow around tandem piers
McIlroy Jr et al. Measurement of turbulent flow phenomena for the lower plenum of a prismatic gas-cooled reactor
Deng et al. Effect of upstream flow characteristics on the wake topology of a square-back truck
Sharma et al. Investigation of bio-inspired sawtooth riblets for boundary layer flow over a flat surface
Niceno et al. Multi-scale modeling and analysis of convective boiling: Towards the prediction of CHF in rod bundles
Fatehi-Nobarian et al. Investigation of the effect of velocity on secondary currents in semicircular channels on hydraulic jump parameters
Abolpour et al. A novel approach for obtaining optimum shape design with the minimum drag coefficient
Yu et al. Flow characteristics of dispersed two-phase flows in an 8× 8 rod bundle
Zhang et al. Three-Dimensional Numerical Simulation Research on Discharge Characteristics of Storm-Drain Inlet in Low-Lying Areas
Du et al. Numerical study of swirling effect on enhancing the bath dynamics in the top-submerged-lance furnace
Johnson CFD investigation of experimental data proposed to be a validation data set
Fan Counter Current Flow Studies Using Interface Capturing Simulations
Lee et al. The accuracy of computational simulation for naturally ventilated multi-span greenhouse

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant