KR101406849B1 - 제어봉 삽입시간 측정을 위한 시뮬레이션 구축 방법 - Google Patents

제어봉 삽입시간 측정을 위한 시뮬레이션 구축 방법 Download PDF

Info

Publication number
KR101406849B1
KR101406849B1 KR1020120155251A KR20120155251A KR101406849B1 KR 101406849 B1 KR101406849 B1 KR 101406849B1 KR 1020120155251 A KR1020120155251 A KR 1020120155251A KR 20120155251 A KR20120155251 A KR 20120155251A KR 101406849 B1 KR101406849 B1 KR 101406849B1
Authority
KR
South Korea
Prior art keywords
control rod
insertion time
falls
reactor
cell
Prior art date
Application number
KR1020120155251A
Other languages
English (en)
Inventor
김경련
이철신
Original Assignee
한국전력기술 주식회사
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 한국전력기술 주식회사 filed Critical 한국전력기술 주식회사
Priority to KR1020120155251A priority Critical patent/KR101406849B1/ko
Priority to EP13866973.4A priority patent/EP2939246B1/en
Priority to PCT/KR2013/012064 priority patent/WO2014104696A1/en
Priority to CN201380035313.4A priority patent/CN104428841B/zh
Priority to US14/411,361 priority patent/US10192006B2/en
Application granted granted Critical
Publication of KR101406849B1 publication Critical patent/KR101406849B1/ko

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/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
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • 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
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/30—Circuit design
    • G06F30/36—Circuit design at the analogue level
    • G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
    • G—PHYSICS
    • G21—NUCLEAR PHYSICS; NUCLEAR ENGINEERING
    • G21C—NUCLEAR REACTORS
    • G21C17/00—Monitoring; Testing ; Maintaining
    • G—PHYSICS
    • G21—NUCLEAR PHYSICS; NUCLEAR ENGINEERING
    • G21C—NUCLEAR REACTORS
    • G21C7/00—Control of nuclear reaction
    • G21C7/06—Control of nuclear reaction by application of neutron-absorbing material, i.e. material with absorption cross-section very much in excess of reflection cross-section
    • G21C7/08—Control of nuclear reaction by application of neutron-absorbing material, i.e. material with absorption cross-section very much in excess of reflection cross-section by displacement of solid control elements, e.g. control rods
    • G—PHYSICS
    • G21—NUCLEAR PHYSICS; NUCLEAR ENGINEERING
    • G21C—NUCLEAR REACTORS
    • G21C7/00—Control of nuclear reaction
    • G21C7/06—Control of nuclear reaction by application of neutron-absorbing material, i.e. material with absorption cross-section very much in excess of reflection cross-section
    • G21C7/08—Control of nuclear reaction by application of neutron-absorbing material, i.e. material with absorption cross-section very much in excess of reflection cross-section by displacement of solid control elements, e.g. control rods
    • G21C7/12—Means for moving control elements to desired position
    • G—PHYSICS
    • G21—NUCLEAR PHYSICS; NUCLEAR ENGINEERING
    • G21D—NUCLEAR POWER PLANT
    • G21D3/00—Control of nuclear power plant
    • G21D3/001—Computer implemented control
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B01—PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
    • B01J—CHEMICAL OR PHYSICAL PROCESSES, e.g. CATALYSIS OR COLLOID CHEMISTRY; THEIR RELEVANT APPARATUS
    • B01J2208/00—Processes carried out in the presence of solid particles; Reactors therefor
    • B01J2208/00008—Controlling the process
    • B01J2208/00548—Flow
    • B01J2208/00557—Flow controlling the residence time inside the reactor vessel
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B01—PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
    • B01J—CHEMICAL OR PHYSICAL PROCESSES, e.g. CATALYSIS OR COLLOID CHEMISTRY; THEIR RELEVANT APPARATUS
    • B01J2208/00—Processes carried out in the presence of solid particles; Reactors therefor
    • B01J2208/00008—Controlling the process
    • B01J2208/00584—Controlling the density
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B01—PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
    • B01J—CHEMICAL OR PHYSICAL PROCESSES, e.g. CATALYSIS OR COLLOID CHEMISTRY; THEIR RELEVANT APPARATUS
    • B01J2208/00—Processes carried out in the presence of solid particles; Reactors therefor
    • B01J2208/00008—Controlling the process
    • B01J2208/00725—Mathematical modelling
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B01—PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
    • B01J—CHEMICAL OR PHYSICAL PROCESSES, e.g. CATALYSIS OR COLLOID CHEMISTRY; THEIR RELEVANT APPARATUS
    • B01J2219/00—Chemical, physical or physico-chemical processes in general; Their relevant apparatus
    • B01J2219/00049—Controlling or regulating processes
    • B01J2219/00051—Controlling the temperature
    • B01J2219/00074—Controlling the temperature by indirect heating or cooling employing heat exchange fluids
    • B01J2219/00076—Controlling the temperature by indirect heating or cooling employing heat exchange fluids with heat exchange elements inside the reactor
    • B01J2219/00081—Tubes
    • 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
    • G21—NUCLEAR PHYSICS; NUCLEAR ENGINEERING
    • G21D—NUCLEAR POWER PLANT
    • G21D3/00—Control of nuclear power plant
    • G21D3/001—Computer implemented control
    • G21D3/002—Core design; core simulations; core optimisation
    • G—PHYSICS
    • G21—NUCLEAR PHYSICS; NUCLEAR ENGINEERING
    • G21D—NUCLEAR POWER PLANT
    • G21D3/00—Control of nuclear power plant
    • G21D3/001—Computer implemented control
    • G21D3/005—Thermo-hydraulic simulations
    • 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
    • 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
    • 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
    • Y10—TECHNICAL SUBJECTS COVERED BY FORMER USPC
    • Y10S—TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y10S585/00—Chemistry of hydrocarbon compounds
    • Y10S585/919—Apparatus considerations
    • Y10S585/921—Apparatus considerations using recited apparatus structure
    • Y10S585/922—Reactor fluid manipulating device

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Plasma & Fusion (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • Mathematical Physics (AREA)
  • Chemical & Material Sciences (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • Data Mining & Analysis (AREA)
  • Algebra (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computational Mathematics (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • Microelectronics & Electronic Packaging (AREA)
  • Computing Systems (AREA)
  • Fluid Mechanics (AREA)
  • Monitoring And Testing Of Nuclear Reactors (AREA)

Abstract

본 발명은 제어봉 삽입시간 측정을 위한 시뮬레이션 구축 방법에 관한 것이다. 제어봉이 중력에 의해 원자로의 노심을 향하여 낙하할 때, 제어봉 삽입시간 측정을 위한 시뮬레이션 구축 방법에 있어서, 원자로 내벽, 상기 원자로에 수용되는 제어봉, 상기 제어봉이 안내되는 안내관, 상기 안내관의 하측에 마련되며 상기 안내관의 직경보다 작은 직경을 갖는 충격흡수관, 상기 안내관의 벽면에 마련되며 상기 원자로에 채워진 물이 유동하는 제1 유동구, 및 상기 충격흡수관의 저면에 마련되어 상기 물이 유동하는 제2 유동구를 3차원으로 모델링하는 단계; 상기 원자로 내부를 복수의 셀로 구획하여 유동장을 형성하는 유동장 형성단계, 여기서 상기 유동장은 상기 제어봉이 지나는 경로 상에 상기 제어봉의 위치가 변함에 따라 형상이 변화하는 복수의 가변셀로 이루어는 이동격자계와, 상기 제어봉의 위치 변화에 무관하게 그 형상이 유지되는 복수의 고정셀로 이루어지는 정렬격자계로 구분하여 형성되며; 3차원 전산유체역학 코드를 이용하여 열수력 현상을 계산하여 상기 제어봉이 상기 충격흡수관에 삽입될 때까지의 삽입시간에 대한 시뮬레이션 예측값을 산출하는 단계; 상기 예측값과 실제 발전소에서 제어봉이 낙하할 때의 삽입시간인 실제값과 비교하고, 상기 예측값과 상기 실제값의 오차범위가 임의의 기준 범위 내에 있는지 판별하며, 상기 오차범위를 벗어나는 경우 상기 가변셀의 크기 및/또는 상기 고정셀의 크기를 변경하는 셀변경 단계;를 포함하는 것을 특징으로 한다.

Description

제어봉 삽입시간 측정을 위한 시뮬레이션 구축 방법{Evaluation methodology for the measurement of control rod insertion time}
본 발명은 제어봉 삽입시간 측정을 위한 시뮬레이션 구축 방법에 관한 것으로, 특히 전산유체역학(Computational Fluid Dynamics)을 이용하여 중력에 의해 원자로 노심에 낙하하는 제어봉의 삽입시간을 측정하는 시뮬레이션을 구축하는 방법에 관한 것으로, 원자로의 3차원 열수력학적 영향을 고려하여 원자력 발전소에서 제어봉이 낙하하는 실제값에 유사한 값을 도출할 수 있도록 한 시뮬레이션을 구축하는 방법에 관한 것이다.
제어봉은 중성자흡수물질을 노심에 삽입, 인출하여 노심 내의 중성자 수를 변화시켜 노심의 반응도를 제어하는 역할을 하는 노심 구성품이다. 즉, 제어봉은 원자로 출력을 제어 및 정지시키는 목적으로 사용되는 부품으로서, 원자로 상단에 위치한 제어봉구동장치에 의해 상하로 이동된다. 제어봉은 제어봉구동장치로부터 낙하신호를 전달받으면 중력에 의해 원자로 노심에 자유낙하 된다. 제어봉 삽입시간은 원자로를 안전하게 운영하기 위해 반드시 평가되어야 하며, 허용기준을 만족하여야 한다.
종래의 제어봉 삽입시간은 일차원 코드를 사용하여 계산되었다. 일차원 코드는 유체의 수력저항, 마찰력 및 제어봉집합체의 무게 등을 고려하지만 제어봉집합체 안내관 내부에서 발생되는 3차원의 수력학적 영향을 전혀 고려하지 않기 때문에 계산의 오차가 큰 단점이 있다. 따라서, 종래의 제어봉 삽입시간은 보수적으로 평가하기 위해 많은 불확실도와 삽입시간의 영향을 미치는 인자에 대해 실제 설계값보다 큰 값을 적용하여 평가되었다. 이로 인해 제어봉 삽입시간은 실제 발전소를 운전하면서 측정된 값에 비해 크게 나타날 수 있다.
제어봉 삽입시간은 원자력발전소의 운전과 안전해석 평가를 수행하기 위해 사용되므로 과도한 보수성을 제거한 실제적인 값을 도출할 필요성이 있다. 제어봉 삽입시간은 실제 원자력발전소에서 측정된 값과 유사하게 평가하기 위해서는 제어봉, 원자로, 안내관 및 충격흡수관 등이 3차원의 공간적 분포를 갖는 열수력 현상을 고려하여야 한다.
한국등록특허공보 제10-1022161호
본 발명은 상술한 바와 같은 문제를 해결하기 위하여 안출된 것으로, 기존에 일차원 코드를 사용하여 제어봉 삽입시간의 오차가 큰 문제점을 해결하기 위하여 제어봉이 낙하할 때의 주변 구조물(원자로, 제어봉, 안내관, 충격흡수관 및 배수구 등)을 실제 형상으로 3차원 모델링하고, 제어봉의 삽입시간은 전산유체역학을 이용하여 3차원의 열수력학적 현상을 모사함으로서 실제 발전소에서 측정된 값과 유사한 값을 도출할 수 있도록 한 제어봉 삽입시간 측정을 위한 시뮬레이션 구축 방법을 제공함을 그 목적으로 한다.
본 발명에 따른 제어봉 삽입시간 측정을 위한 시뮬레이션 구축 방법은, 제어봉이 중력에 의해 원자로의 노심을 향하여 낙하할 때, 제어봉 삽입시간 측정을 위한 시뮬레이션 구축 방법에 있어서, 원자로 내벽, 상기 원자로에 수용되는 제어봉, 상기 제어봉이 안내되는 안내관, 상기 안내관의 하측에 마련되며 상기 안내관의 직경보다 작은 직경을 갖는 충격흡수관, 상기 안내관의 벽면에 마련되며 상기 원자로에 채워진 물이 유동하는 제1 유동구, 및 상기 충격흡수관의 저면에 마련되어 상기 물이 유동하는 제2 유동구를 3차원으로 모델링하는 단계; 상기 원자로 내부를 복수의 셀로 구획하여 유동장을 형성하는 유동장 형성단계, 여기서 상기 유동장은 상기 제어봉이 지나는 경로 상에 상기 제어봉의 위치가 변함에 따라 형상이 변화하는 복수의 가변셀로 이루어는 이동격자계와, 상기 제어봉의 위치 변화에 무관하게 그 형상이 유지되는 복수의 고정셀로 이루어지는 정렬격자계로 구분하여 형성되며; 3차원 전산유체역학 코드를 이용하여 열수력 현상을 계산하여 상기 제어봉이 상기 충격흡수관에 삽입될 때까지의 삽입시간에 대한 시뮬레이션 예측값을 산출하는 단계; 상기 예측값과 실제 발전소에서 제어봉이 낙하할 때의 삽입시간인 실제값과 비교하고, 상기 예측값과 상기 실제값의 오차가 임의의 기준범위 내에 있는지 판별하며, 상기 오차가 상기 기준범위를 벗어나는 경우 상기 가변셀의 크기 및/또는 상기 고정셀의 크기를 변경하는 셀변경 단계;를 포함하는 것을 특징으로 한다.
또한, 상기 제어봉이 낙하할 때, 상기 원자로에 채워진 물에 의한 상기 제어봉의 수력 저항을 계산하는 단계를 포함하는 것이 바람직하다.
또한, 상기 제어봉이 낙하할 때, 상기 물과 상기 제어봉의 마찰력 및 상기 제어봉의 무게를 고려하는 것이 바람직하다.
또한, 상기 기준범위는 상기 예측값과 상기 실제값이 5% 범위 내의 오차를 갖도록 설정되는 것이 바람직하다.
또한, 상기 열수력 현상으로는, 상기 제어봉이 낙하할 때 상기 제어봉 주변의 물의 압력변화, 온도변화, 및 밀도변화를 포함하는 것이 바람직하다.
본 발명에 따른 제어봉 삽입시간 측정을 위한 시뮬레이션 구축 방법은, 전산유체역학을 이용하여 제어봉이 낙하하는 과정을 3차원적 열수력학적 특성을 고려하여 모사하므로, 시뮬레이션에 의해 산출된 제어봉의 삽입시간(예측값)을 실제 발전소에서 제어봉이 낙하할 때의 실제값에 매우 근사한 값으로 확보할 수 있다.
또한, 제어봉의 삽입시간을 실제에 근사한 값으로 시뮬레이션하여, 제어봉의 낙하 속도, 안내관 및 충격흡수관의 압력분포, 제1,2 유동구를 통해 토출되는 유량에 관한 정보를 추가적으로 획득하여 원자로 설계 및 운전에 유용한 정보로 활용할 수 있다.
또한, 본 발명은 신규 원자력 발전소의 설계 및 운전은 물론, 현재 운전 중인 발전소에도 적용할 수 있으므로, 운전 중인 원자력 발전소의 운전에 유용한 정보를 제공할 수 있다.
도1은 본 발명에 따른 시뮬레이션 구축 방법의 흐름도,
도2는 제어봉 낙하 과정을 개략적으로 도시한 도면,
도3 및 도4는 제어봉의 위치에 따라 이동격자계와 정렬격자계를 개략적으로 도시한 도면,
제어봉은 제어봉 구동장치로부터 낙하신호가 수신되면 중력에 의해 원자로의 노심을 향하여 자유 낙하한다. 본 발명은 제어봉이 낙하하여 충격흡수관까지 삽입될 때까지의 시간을 시뮬레이션하는 방법에 관한 것이다.
이하, 본 발명에 따른 바람직한 실시예를 첨부된 도면을 참조하여 상세히 설명한다.
도1은 본 발명에 따른 시뮬레이션 구축 방법의 흐름도이고, 도2는 제어봉 낙하 과정을 개략적으로 도시한 도면이다. 도3 및 도4는 제어봉의 위치에 따라 이동격자계와 정렬격자계를 개략적으로 도시한 도면이다.
먼저, 도1를 참조하면, 본 발명의 실시예에 따른 제어봉 삽입시간 측정을 위한 시뮬레이션 구축 방법은, 3차원 모델링 단계, 유동장 형성 단계, 시뮬레이션 예측값 산출단계, 및 셀변경 단계를 포함한다.
상기 3차원 모델링 단계는, 제어봉(10)이 낙하하는 주변 구조물들을 3차원으로 모델링하는 단계이다. 본 발명의 실시예에 따르면, 원자로 내벽(20), 제어봉(10), 안내관(30), 충격흡수관(40), 제1 유동구(50) 및 제2 유동구(60)가 3차원으로 모델링된다.
상기 원자로는 제어봉(10)을 수용하며, 상기 안내관(30)은 상기 제어봉(10)이 자유 낙하할 때 상기 제어봉(10)을 안내한다. 상기 제어봉(10)은 상기 안내관(30)을 따라 하측으로 낙하한다. 상기 충격흡수관(40)은 상기 안내관(30)의 하측에 마련되며, 상기 안내관(30)의 직경보다 작은 직경을 갖는다.
상기 제어봉(10)은 상기 안내관(30)을 따라 낙하할 때 상기 충격흡수관(40)에 이르기까지는 급격하게 가속운동을 한다. 상기 제어봉(10)이 충격흡수관(40)에 진입하는 순간 제어봉(10)은 수력저항력이 증가하여 제어봉(10)이 낙하하는 속도가 감속된다. 즉, 안내관(30)의 직경보다 작은 직경을 갖는 충격흡수관(40)이 마련되어 상기 제어봉(10)이 충격흡수관(40)을 지날 때 하강 속도가 늦어지게 된다.
상기 제1 유동구(50)는 상기 안내관(30)의 벽면에 마련된다. 상기 제어봉(10)이 안내관(30)을 따라 낙하할 때 상기 제어봉(10)이 안내관(30)에 채워진 물을 밀어내게 되는데, 이때 상기 물은 상기 제1 유동구(50)를 통하여 안내관(30) 외측으로 유동한다. 즉, 상기 물은 원자로에 채워져 있고, 제어봉(10)이 낙하할 때 안내관(30)에 존재하는 물이 제1 유동구(50)를 통하여 안내관(30) 외측으로 유동하게 된다.
상기 제2 유동구(60)는 상기 충격흡수관(40)의 저면에 마련된다. 상기 충격흡수관(40)으로 진입한 제어봉(10)이 하측으로 더 낙하할 때 상기 충격흡수관(40)에 채워진 물은 제2 유동구(60)를 통하여 유동한다. 즉, 물이 제2 유동구(60)를 통하여 충격흡수관(40)의 외측으로 유동하게 된다. 이와 같이, 제어봉(10)이 낙하할 때, 제어봉(10)과 관련된 주변의 주요 구조물이 3차원으로 모델링된다.
이어서, 유동장을 형성한다. 상기 유동장은 이동격자계(80)와, 정렬격자계(70)를 포함한다.
상기 유동장이라함은 상기 원자로 내부에서 상기 제어봉(10)이 낙하할 때 제어봉(10)의 낙하에 따라 영향을 받는 영역을 의미하는 것으로, 원자로 내부를 복수의 셀로 구획하여 유동장을 형성한다.
제어봉(10)의 삽입 시간을 측정할 수 있도록, 노심 영역은 여러 개의 작은 셀(영역)로 분할된다. 이러한 작은 셀들을 격자라고 한다.
이동격자계(80)는 제어봉(10)이 지나는 경로 상에 상기 제어봉(10)의 위치가 변함에 따라 형상이 변화하는 복수의 셀의 집합이다. 이때, 제어봉(10)의 위치가 변함에 따라 형상이 변화하는 각 셀을 가변셀(81)이라 부르기로 한다. 즉, 이동격자계(80)는 제어봉(10)이 낙하할 때 제어봉(10)이 지나는 경로 상에 마련되는 것으로 가변셀(81)의 집합이다.
일반적으로, 이동격자계(80)는 유동장 내에 상대적으로 움직이는 물체가 있거나 유동장 내에서 어떤 값의 영향을 받아 유동장의 형상이 변화할 때 정해진 시간간격(time step)마다 격자를 변화시켜 계산을 수행하여 얻어진다.
즉, 본 실시예에 따르면, 제어봉(10)이 중력에 의해 노심으로 낙하하면서 제어봉(10)이 낙하하는 이동영역의 격자는 제어봉(10)과 함께 시간의 변화에 따라 움직이게 되므로, 제어봉(10)의 위치에 따라서 가변셀(81)이 변화하는 것으로 구현한 것이다.
본 실시예에 따르면, 도3 및 도4에 도시된 바와 같이 이동격자계(80)가 주어진다. 도3은 제어봉(10)의 상단부를 도시한 것으로, 낙하 전에 제어봉(10)의 상단부가 위치한 지점에서부터, 낙하 후 제어봉(10)의 상단부가 위치한 지점까지 제어봉(10)의 위치 변화에 따라 영향을 받는 이동격자계(80)가 형성되며, 이 부분에 가변셀(81)을 이용하여 이동격자계(80)를 구성한다. 또한, 도4는 제어봉(10)의 하단부를 도시한 것으로, 제어봉(10)이 낙하함에 따라 제어봉(10)의 하단부가 지나가는 경로를 이동격자계(80)로 구성한 것을 도시한 것이다.
상기 정렬격자계(70)는 상기 유동장 중에서 상기 이동격자계(80) 이외의 영역을 복수의 셀로 분할하여 형성한 것이다. 상기 정렬격자계(70)는 상기 제어봉(10)의 위치 변화에 무관하게 그 형상이 유지되는 복수의 고정셀(71)의 집합이다.
상기 정렬격자계(70)를 구성하는 각 고정셀(71)은 직사각형으로 이루어져 규칙적으로 분포된다. 참고적으로, 상기 정렬격자계(70)와 대비되는 용어인 비정렬격자계는, 직사각형 이외의 형상을 가지는 셀, 예컨대 직각인 변을 가지지 않는 사각형이나, 삼각형의 모양을 갖는 비정렬격자를 불규칙적으로 배치한 것을 일컫는다. 본 발명은 계산의 정확도 및 계산 시간의 단축을 위해서 정렬격자계(70)를 사용하며, 제어봉(10)이 이동하는 영역은 이동격자계(80)로 구성한다.
상기 시뮬레이션 예측값 산출단계는, 3차원 전산유체역학 코드를 이용하여 열수력 현상을 계산함으로써 상기 제어봉(10)이 상기 충격흡수관(40)에 삽입될 때까지의 삽입시간을 산출하는 단계이다.
상기 열수력 현상이라함은 임의의 물체가 떨어질 때 상기 물체의 주변에서 발생되는 물리적 변화를 의미하는 것으로, 본 실시예에서 제어봉(10)이 충격흡수관(40)까지 낙하할 때 제어봉(10) 주변에서의 물의 압력변화, 온도변화, 및 밀도변화 등을 의미한다.
따라서, 상기 3차원 전산유체역학 코드를 이용하여 열수력 현상을 계산할 때, 원자로 내부에 채워진 물의 압력, 온도, 및 밀도 등의 변화 양상이 고려되어 계산된다. 3차원 전산유체역학 코드를 이용하여 열수력 현상을 계산하는 방법은 공지된 방법을 사용하므로, 그 구체적인 설명은 생략한다.
상기 셀변경 단계는 상기 예측값이 실제 발전소에서 제어봉(10)이 낙하하여 충격흡수관(40)까지 삽입될 때까지 소요되는 실제값과 한계 오차 범위를 초과하는 경우에 가변셀(81) 및/또는 고정셀(71)의 크기를 변경하여 오차범위를 축소하기 위해서 수행된다.
먼저, 상기 예측값과 실제 발전소에서 제어봉(10)이 낙하할 때의 삽입시간인 실제값을 비교하고, 상기 예측값과 실제값의 오차가 임의의 기준 범위 내에 있는지 판별한다. 본 실시예에서, 기준범위는 상기 예측값과 상기 실제값이 5% 범위 내의 오차를 갖도록 설정된다. 기분범위를 5%를 초과하여 설정하면 예측값에 과하게 보수적으로 평가되는 결과가 되므로 바람직하지 않다.
상기 오차가 미리 설정된 기준범위를 벗어나는 경우 상기 가변셀(81)의 크기 및/또는 상기 고정셀(71)의 크기를 변경하는 셀변경 단계가 수행된다. 여기서, "및/또는"의 의미는, 가변셀(81)의 크기와 고정셀(71)의 크기를 동시에 변경하거나, 가변셀(81)의 크기 또는 고정셀(71)의 크기를 선택하여 어느 하나의 셀 크기를 변경할 수 있음을 의미한다.
제어봉(10)이 낙하할 때, 제어봉(10)의 주변에 배치되는 구조물은 매우 복잡하고, 정교한 구조물들이기 때문에 이동격자계(80)를 구성하는 가변셀(81) 및/또는 정렬격자계(70)를 구성하는 고정셀(71)의 조밀도를 최적화하도록 조절함으로써 정확한 예측값을 산출할 수 있게 된다.
또한, 본 발명의 실시예에 따른, 제어봉(10) 삽입시간 측정을 위한 시뮬레이션 구축 방법은, 수력 저항을 계산하는 단계 및 마찰력, 제어봉(10)의 무게를 고려하는 단계를 더 포함한다.
상기 수력저항 계산 단계는, 상기 제어봉(10)이 낙하할 때 상기 원자로에 채워진 물에 의한 상기 제어봉(10)의 수력 저항에 대한 영향을 고려하는 것으로, 3차원 열수력학적 현상을 모사하는 것과 더불어 실제 물에 의한 저항을 예측값의 산출에 파라미터로 고려함으로써 예측값을 더욱 정확하게 산출한다.
제어봉(10)이 중력에 의해 상기 충격흡수관(40)까지 낙하할 때, 물은 안내관(30) 및 충격흡수관(40)이 형성하는 좁은 유로를 따라 수직으로 상승한다. 이처럼 수직 방향으로 물이 유동하는 경우에 Darcy 방정식을 적용하여 수력 저항을 계산할 수 있으며, 구체적으로 Darcy 방정식을 적용하여 수력 저항을 계산해 내는 과정은 공지된 바에 의하며, 계산 과정 자체는 본 발명의 기술적 사상의 범위 밖에 있으므로 그 구체적인 설명은 생략한다.
상기 마찰력, 제어봉(10)의 무게를 고려하는 단계는, 제어봉(10)의 삽입 시간에 대한 예측값을 산출하는 또 다른 파라미터로서 작용한다.
상기 마찰력은, 제어봉(10)이 물이 채워진 안내관(30)과 충격흡수관(40)을 따라 낙하할 때, 물의 점성계수와 안내관(30) 벽면의 표면거칠기에 의해 발생되는 마찰력을 의미한다. 즉, 제어봉(10)이 낙하할 때 제어봉(10)과 접속하는 물의 표면에는 제어봉(10)이 낙하하는 방향과 반대방향으로 마찰력이 발생한다. 상기 마찰력은 제어봉(10)이 낙하하는 방향과 반대방향으로 작용하므로, 제어봉(10)의 삽입시간을 지연시키는 인자로 작용한다. 따라서, 상기 마찰력을 정확한 예측값 산출을 위해서 고려하는 것이다.
상기 제어봉(10)의 무게를 실제 발전소에서 사용되는 제어봉(10)의 무게를 파리미터로 적용하여 예측값 산출에 사용한다.
이처럼, 본 발명 실시예에 따른 제어봉 삽입시간 측정을 위한 시뮬레이션 구축 방법은, 제어봉(10)이 원자로의 상부에서 충격흡수관(40)까지 자유 낙하하는 과정을 실제 모습과 유사하게 시뮬레이션함으로써 제어봉(10) 삽입시간의 예측값을 실제 발전소에서 제어봉(10)이 낙하하는 실제 삽입시간과 유사하게 예측하는 효과를 제공한다.
또한, 예측값을 산출하는 파라미터로써, 수력저항, 마찰력, 및 제어봉(10)의 무게를 고려함으로써 예측값의 정확도를 더욱 향상시키는 효과를 제공한다.
또한 실제값에 근사한 예측값에 근거하여, 제어봉(10)의 낙하시 제어봉(10)의 낙하 속도, 안내관(30) 및 충격흡수관(40)의 압력분포, 제1,2 유동구(50,60)를 통해 토출되는 유량 정보를 추가적으로 획득하여 원자로 설계 및 운전에 유용한 정보를 제공할 수 있다.
이상, 본 발명을 바람직한 실시예를 들어 상세하게 설명하였으나, 본 발명은 상기 실시예에 한정되지 않으며, 본 발명의 범주를 벗어나지 않는 범위 내에서 여러 가지 많은 변형이 제공될 수 있다.
10... 제어봉 20... 원자로 내벽
30... 안내관 40... 충격흡수관
50... 제1 유동구 60... 제2 유동구
70... 정렬격자계 71... 고정셀
80... 이동격자계 81... 가변셀

Claims (5)

  1. 제어봉이 중력에 의해 원자로의 노심을 향하여 낙하할 때, 제어봉 삽입시간 측정을 위한 시뮬레이션 구축 방법에 있어서,
    원자로 내벽(20), 상기 원자로에 수용되는 제어봉(10), 상기 제어봉(10)이 안내되는 안내관(30), 상기 안내관(30)의 하측에 마련되며 상기 안내관(30)의 직경보다 작은 직경을 갖는 충격흡수관(40), 상기 안내관(30)의 벽면에 마련되며 상기 원자로에 채워진 물이 유동하는 제1 유동구(50), 및 상기 충격흡수관(40)의 저면에 마련되어 상기 물이 유동하는 제2 유동구(60)를 3차원으로 모델링하는 단계;
    상기 원자로 내부를 복수의 셀로 구획하여 유동장을 형성하는 유동장 형성단계, 여기서 상기 유동장은 상기 제어봉(10)이 지나는 경로 상에 상기 제어봉(10)의 위치가 변함에 따라 형상이 변화하는 복수의 가변셀(81)로 이루어는 이동격자계(80)와, 상기 제어봉(10)의 위치 변화에 무관하게 그 형상이 유지되는 복수의 고정셀(71)로 이루어지는 정렬격자계(70)로 구분하여 형성되며;
    3차원 전산유체역학 코드를 이용하여 열수력 현상을 계산하여 상기 제어봉(10)이 상기 충격흡수관(40)에 삽입될 때까지의 삽입시간에 대한 시뮬레이션 예측값을 산출하는 단계;
    상기 예측값과 실제 발전소에서 제어봉(10)이 낙하할 때의 삽입시간인 실제값과 비교하고, 상기 예측값과 상기 실제값의 오차가 임의의 기준범위 내에 있는지 판별하며, 상기 오차가 상기 기준범위를 벗어나는 경우 상기 가변셀(81)의 크기와 상기 고정셀(71)의 크기를 함께 변경하거나, 상기 가변셀(81) 또는 상기 고정셀(71) 중 어느 하나의 크기를 변경하는 셀변경 단계;를 포함하는 것을 특징으로 하는 제어봉 삽입시간 측정을 위한 시뮬레이션 구축 방법.
  2. 제1항에 있어서,
    상기 제어봉(10)이 낙하할 때, 상기 원자로에 채워진 물에 의한 상기 제어봉(10)의 수력 저항을 계산하는 단계를 포함하는 것을 특징으로 하는 제어봉 삽입시간 측정을 위한 시뮬레이션 구축 방법.
  3. 제1항에 있어서,
    상기 제어봉(10)이 낙하할 때, 상기 물과 상기 제어봉(10)의 마찰력 및 상기 제어봉(10)의 무게를 고려하는 것을 특징으로 하는 제어봉 삽입시간 측정을 위한 시뮬레이션 구축 방법.
  4. 제1항에 있어서,
    상기 기준범위는 상기 예측값과 상기 실제값이 5% 범위 내의 오차를 갖도록 설정되는 것을 특징으로 하는 제어봉 삽입시간 측정을 위한 시뮬레이션 구축 방법.
  5. 제1항에 있어서,
    상기 열수력 현상으로는, 상기 제어봉(10)이 낙하할 때 상기 제어봉(10) 주변의 물의 압력변화, 온도변화, 및 밀도변화를 포함하는 것을 특징으로 하는 제어봉 삽입시간 측정을 위한 시뮬레이션 구축 방법.
KR1020120155251A 2012-12-27 2012-12-27 제어봉 삽입시간 측정을 위한 시뮬레이션 구축 방법 KR101406849B1 (ko)

Priority Applications (5)

Application Number Priority Date Filing Date Title
KR1020120155251A KR101406849B1 (ko) 2012-12-27 2012-12-27 제어봉 삽입시간 측정을 위한 시뮬레이션 구축 방법
EP13866973.4A EP2939246B1 (en) 2012-12-27 2013-12-24 Simulation construction method for the measurement of control rod insertion time
PCT/KR2013/012064 WO2014104696A1 (en) 2012-12-27 2013-12-24 Simulation construction method for the measurement of control rod insertion time
CN201380035313.4A CN104428841B (zh) 2012-12-27 2013-12-24 测量控制棒插入时间的仿真构建方法
US14/411,361 US10192006B2 (en) 2012-12-27 2013-12-24 Simulation construction method for the measurement of control rod insertion time

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR1020120155251A KR101406849B1 (ko) 2012-12-27 2012-12-27 제어봉 삽입시간 측정을 위한 시뮬레이션 구축 방법

Publications (1)

Publication Number Publication Date
KR101406849B1 true KR101406849B1 (ko) 2014-06-13

Family

ID=51021669

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020120155251A KR101406849B1 (ko) 2012-12-27 2012-12-27 제어봉 삽입시간 측정을 위한 시뮬레이션 구축 방법

Country Status (5)

Country Link
US (1) US10192006B2 (ko)
EP (1) EP2939246B1 (ko)
KR (1) KR101406849B1 (ko)
CN (1) CN104428841B (ko)
WO (1) WO2014104696A1 (ko)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105045965B (zh) * 2015-06-26 2018-08-28 西安交通大学 基于仿真数据库的压力容器lpms质量估计方法
US11129906B1 (en) 2016-12-07 2021-09-28 David Gordon Bermudes Chimeric protein toxins for expression by therapeutic bacteria
WO2019015805A1 (en) * 2017-07-19 2019-01-24 Linde Aktiengesellschaft METHOD FOR DETERMINING STRAIN LEVELS IN MATERIAL OF PROCESS ENGINEERING APPARATUS
CN113177371A (zh) * 2021-04-20 2021-07-27 哈尔滨工程大学 一种反应堆堆芯组件流域流场顺次重构的cfd加速计算方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR100568562B1 (ko) 2003-11-07 2006-04-07 강문구 연속체 유체역학 해석기법과 볼륨렌더링 기법을 이용한 실시간 유체유동 시뮬레이션 및 렌더링 방법
KR100957061B1 (ko) 2009-11-27 2010-05-13 (주)에네시스 3차원 전산유체역학을 이용한 원자로 열수력 안전해석 방법
KR20110014399A (ko) * 2009-08-05 2011-02-11 한국전력공사 제어봉의 낙하 시간을 측정하기 위한 방법 및 컴퓨터 판독 가능한 기록 매체
KR101022161B1 (ko) 2010-08-30 2011-03-24 주식회사 리얼게인 제어봉 낙하시간 측정방법

Family Cites Families (22)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS54120388A (en) * 1978-03-10 1979-09-18 Hitachi Ltd Control rod driving mechanism of reactor
JPS5552998A (en) * 1978-10-16 1980-04-17 Hitachi Ltd Reactor recirculation flow rate control device
JPS5723893A (en) * 1980-07-18 1982-02-08 Tokyo Shibaura Electric Co Control rod insertion time judging device
JP3021100B2 (ja) * 1991-07-22 2000-03-15 株式会社日立製作所 ギャング制御棒制御装置及び原子炉の運転方法
EP0796495B1 (de) * 1994-12-08 1998-08-12 Siemens Aktiengesellschaft System und verfahren zur steuerung von steuerstäben einer kernkraftanlage
US5677938A (en) * 1995-03-13 1997-10-14 Peco Energy Company Method for fueling and operating a nuclear reactor core
JP3892572B2 (ja) 1998-02-20 2007-03-14 株式会社東芝 制御棒操作監視制御システムの試験方法および試験装置
JP3069558B1 (ja) * 1999-07-30 2000-07-24 株式会社日立製作所 制御棒案内管清掃装置
US6959269B1 (en) * 2001-04-20 2005-10-25 Lockheed Martin Corporation Method and system for simulating an unsteady flow field
US6824305B1 (en) * 2002-08-16 2004-11-30 The Texas A & M University System Local wall heat flux/temperature meter for convective flow and method of utilizing same
JP4197617B2 (ja) * 2003-01-23 2008-12-17 原子燃料工業株式会社 炉心制御棒反応度効果評価方法
US7620139B2 (en) 2005-06-30 2009-11-17 Global Nuclear Fuel - Americas, Llc Method of improving nuclear reactor performance during reactor core operation
US8010326B2 (en) * 2005-12-28 2011-08-30 Caterpillar Inc. Method and apparatus for automated grid formation in multi-cell system dynamics models
KR100927124B1 (ko) 2008-01-29 2009-11-18 부산대학교 산학협력단 원자로의 압력손실계수 계산방법
US8163158B2 (en) * 2008-05-12 2012-04-24 Enrg, Inc. Operation of an electrolysis cell
RU2553468C2 (ru) 2009-11-06 2015-06-20 ТерраПауэр, ЭлЭлСи Системы и способы регулирования реактивности в реакторе ядерного деления
KR101032867B1 (ko) * 2010-01-28 2011-05-06 한전원자력연료 주식회사 핵연료집합체용 안내관 마개
KR101250749B1 (ko) * 2011-01-18 2013-04-04 한국수력원자력 주식회사 제어봉제어계통 코일전류 감시장치를 이용한 제어봉 낙하시간 측정방법
US8797021B2 (en) * 2011-02-14 2014-08-05 Ge-Hitachi Nuclear Energy Americas Llc Electrochemical corrosion potential probe assembly
FR3004558B1 (fr) * 2013-04-10 2015-05-15 Areva Np Procedes de simulation de l'ecoulement d'un fluide dans une cuve d'un reacteur nucleaire et de calcul de la deformation mecanique d'assemblages d'un coeur de reacteur nucleaire, et produits programme d'ordinateur associes
US9122934B2 (en) * 2013-12-27 2015-09-01 Automotive Research & Testing Center Object detection method with a rising classifier effect and object detection device with the same
CA2891002C (en) * 2015-05-13 2022-09-06 Veronica Rose Zimmerman Modeling a bed plate and its use

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR100568562B1 (ko) 2003-11-07 2006-04-07 강문구 연속체 유체역학 해석기법과 볼륨렌더링 기법을 이용한 실시간 유체유동 시뮬레이션 및 렌더링 방법
KR20110014399A (ko) * 2009-08-05 2011-02-11 한국전력공사 제어봉의 낙하 시간을 측정하기 위한 방법 및 컴퓨터 판독 가능한 기록 매체
KR100957061B1 (ko) 2009-11-27 2010-05-13 (주)에네시스 3차원 전산유체역학을 이용한 원자로 열수력 안전해석 방법
KR101022161B1 (ko) 2010-08-30 2011-03-24 주식회사 리얼게인 제어봉 낙하시간 측정방법

Also Published As

Publication number Publication date
US10192006B2 (en) 2019-01-29
US20150294042A1 (en) 2015-10-15
CN104428841A (zh) 2015-03-18
EP2939246B1 (en) 2017-06-07
WO2014104696A1 (en) 2014-07-03
CN104428841B (zh) 2016-09-21
EP2939246A1 (en) 2015-11-04
EP2939246A4 (en) 2016-08-17

Similar Documents

Publication Publication Date Title
KR101406849B1 (ko) 제어봉 삽입시간 측정을 위한 시뮬레이션 구축 방법
Moonen et al. Performance assessment of Large Eddy Simulation (LES) for modeling dispersion in an urban street canyon with tree planting
Zeng et al. Applications of computational fluid dynamics to flow ratings at prototype spillways and weirs. I: Data generation and validation
Matsuura et al. MPS simulation of spreading behavior of molten material
Kaouachi et al. Numerical assessment of the inception point in different stepped spillway configurations
Chen et al. Comparing k-ε models on predictions of an impinging jet for ventilation of an office room
Manera et al. NEAMS IRP challenge problem 3: Mixing in large enclosures and thermal stratification
Al-Khatib et al. Development of empirical regression-based models for predicting mean velocities in asymmetric compound channels
CN116362001A (zh) 基于滑体质心高度势能计算高位滑坡入水速度及涌浪方法
Khassaf et al. Experimental investigation of compound side weir with modeling using computational fluid dynamic
Madarame et al. Gas entrainment inception at the border of a flow-swollen liquid surface
Stamataki et al. Study of dam break flow interaction with urban settlements over a sloping channel
Ho et al. Hydrodynamic modelling study of a fish exclusion system for a river diversion
Nakayama et al. Large-eddy simulation of turbulent winds during the Fukushima Daiichi Nuclear Power Plant accident by coupling with a meso-scale meteorological simulation model
Souli et al. Investigating supercritical bended flow using physical model and CFD
Ishay et al. Nuclear reactor containment flows-Modelling of stably stratified layer erosion by a turbulent jet
Hujala Quantification of large steam bubble oscillations and chugging using image analysis
Raiford et al. Turbulence schemes for modelling a submerged hydraulic jump
Khajehdezfuly et al. Dynamic Response of Inter-story Isolated High-rise Building Subjected to Wind Flow
Michalcova et al. Numerical and Experimental Models of Load on Buldings from the Effects of the Atmospheric Wind
Sarchami et al. Temperature fluctuations inside the CANDU reactor Moderator Test Facility (MTF)
Bosselut et al. Insertion and drop of control rod in assembly simulations and parametric analysis
Mochizuki et al. Development of ATR Type Steam Separator Performance Analysis Code and Its Validation,(I) Carryunder Characteristics
Seong et al. Development of a Partial Blockage Detection Algorithm through Temperature Fluctuation Analyses in a Liquid Metal Reactor
Hihara et al. Liquid Film Flow on Vertical and Successive Inclined Plate Modeling Flow Channel of Molten Control Rod in Severe Accident

Legal Events

Date Code Title Description
A201 Request for examination
E902 Notification of reason for refusal
E701 Decision to grant or registration of patent right
GRNT Written decision to grant
FPAY Annual fee payment

Payment date: 20170515

Year of fee payment: 4

FPAY Annual fee payment

Payment date: 20180515

Year of fee payment: 5

FPAY Annual fee payment

Payment date: 20190529

Year of fee payment: 6