• 제목/요약/키워드: Heat conduction model

검색결과 241건 처리시간 0.032초

스위스 Mont Terri rock laboratory에서 수행된 암반 히터시험(HE-D)에 대한 열-수리-역학적 복합거동 수치해석 (Numerical modelling of coupled thermo-hydro-mechanical behavior of Heater Experiment-D (HE-D) at Mont Terri rock laboratory in Switzerland)

  • 이창수;최희주;김건영
    • 터널과지하공간
    • /
    • 제30권3호
    • /
    • pp.242-255
    • /
    • 2020
  • 본 연구에서는 FLAC3D가 Opalinus Clay 암반의 열-수리-역학적 복합거동을 재현하고 이를 예측할 수 있는지 그 적용성을 검토하고자 국제공동연구 DECOVALEX-2015에서 참여하였으며, 그 일환으로 스위스 Mont Terri Rock Laboratory에서 수행된 Heater Experiment-D (HE-D)에 대한 모델링을 수행하였다. FLAC3D를 이용한 수치해석의 타당성을 평가하기 위해 현장시험에서 계측된 16 지점의 온도, 6 지점의 공극수압, 그리고 22 지점의 변형률 데이터와 비교하였다. 대상 암반의 열-수리-역학적인 이방성을 고려함으로써 Opalinus Clay 암반의 온도 변화 그리고 온도변화에 따른 공극수압의 변화와 같은 열-수리적 거동은 전반적으로 유사하게 나타났으나, 역학적 거동의 경우 변형률 데이터를 비교했을 때 온도와 공극수압과는 달리 계산된 변형률 일부만이 유사한 거동을 보였다.

Development and validation of multiphysics PWR core simulator KANT

  • Taesuk Oh;Yunseok Jeong;Husam Khalefih;Yonghee Kim
    • Nuclear Engineering and Technology
    • /
    • 제55권6호
    • /
    • pp.2230-2245
    • /
    • 2023
  • KANT (KAIST Advanced Nuclear Tachygraphy) is a PWR core simulator recently developed at Korea Advance Institute of Science and Technology, which solves three-dimensional steady-state and transient multigroup neutron diffusion equations under Cartesian geometries alongside the incorporation of thermal-hydraulics feedback effect for multi-physics calculation. It utilizes the standard Nodal Expansion Method (NEM) accelerated with various Coarse Mesh Finite Difference (CMFD) methods for neutronics calculation. For thermal-hydraulics (TH) calculation, a single-phase flow model and a one-dimensional cylindrical fuel rod heat conduction model are employed. The time-dependent neutronics and TH calculations are numerically solved through an implicit Euler scheme, where a detailed coupling strategy is presented in this paper alongside a description of nodal equivalence, macroscopic depletion, and pin power reconstruction. For validation of the steady, transient, and depletion calculation with pin power reconstruction capacity of KANT, solutions for various benchmark problems are presented. The IAEA 3-D PWR and 4-group KOEBERG problems were considered for the steady-state reactor benchmark problem. For transient calculations, LMW (Lagenbuch, Maurer and Werner) LWR and NEACRP 3-D PWR benchmarks were solved, where the latter problem includes thermal-hydraulics feedback. For macroscopic depletion with pin power reconstruction, a small PWR problem modified with KAIST benchmark model was solved. For validation of the multi-physics analysis capability of KANT concerning large-sized PWRs, the BEAVRS Cycle1 benchmark has been considered. It was found that KANT solutions are accurate and consistent compared to other published works.

Validation of a New Design of Tellurium Dioxide-Irradiated Target

  • Fllaoui, Aziz;Ghamad, Younes;Zoubir, Brahim;Ayaz, Zinel Abidine;Morabiti, Aissam El;Amayoud, Hafid;Chakir, El Mahjoub
    • Nuclear Engineering and Technology
    • /
    • 제48권5호
    • /
    • pp.1273-1279
    • /
    • 2016
  • Production of iodine-131 by neutron activation of tellurium in tellurium dioxide ($TeO_2$) material requires a target that meets the safety requirements. In a radiopharmaceutical production unit, a new lid for a can was designed, which permits tight sealing of the target by using tungsten inert gaswelding. The leakage rate of all prepared targets was assessed using a helium mass spectrometer. The accepted leakage rate is ${\leq}10^{-4}mbr.L/s$, according to the approved safety report related to iodine-131 production in the TRIGA Mark II research reactor (TRIGA: Training, Research, Isotopes, General Atomics). To confirm the resistance of the new design to the irradiation conditions in the TRIGA Mark II research reactor's central thimble, a study of heat effect on the sealed targets for 7 hours in an oven was conducted and the leakage rates were evaluated. The results show that the tightness of the targets is ensured up to $600^{\circ}C$ with the appearance of deformations on lids beyond $450^{\circ}C$. The study of heat transfer through the target was conducted by adopting a one-dimensional approximation, under consideration of the three transfer modes-convection, conduction, and radiation. The quantities of heat generated by gamma and neutron heating were calculated by a validated computational model for the neutronic simulation of the TRIGA Mark II research reactor using the Monte Carlo N-Particle transport code. Using the heat transfer equations according to the three modes of heat transfer, the thermal study of I-131 production by irradiation of the target in the central thimble showed that the temperatures of materials do not exceed the corresponding melting points. To validate this new design, several targets have been irradiated in the central thimble according to a preplanned irradiation program, going from4 hours of irradiation at a power level of 0.5MWup to 35 hours (7 h/d for 5 days a week) at 1.5MW. The results showthat the irradiated targets are tight because no iodine-131 was released in the atmosphere of the reactor building and in the reactor cooling water of the primary circuit.

한국잔디식재 경량박층형 옥상녹화의 열수지 해석 (Heat Budget Analysis of Light Thin Layer Green Roof Planted with Zoysia japonica)

  • 김세창;이현정;박봉주
    • 한국조경학회지
    • /
    • 제40권6호
    • /
    • pp.190-197
    • /
    • 2012
  • 본 연구는 저관리 경량형 옥상녹화의 열수지를 정량화하기 위하여 모의실험을 수행하였다. 콘크리트블록을 이용하여 가로 $1.2m{\times}$세로 $1.2m{\times}$높이 1m의 모의실험구 2개를 제작하여 하나는 한국잔디를 식재한 토심 10cm의 옥상녹화 실험구로 하였으며, 옥상녹화를 실시하지 않은 나머지 하나는 대조구로 하였다. 2012년 6월 6일~7일에 기온, 상대습도, 풍향, 풍속 등 외부 기상환경, 증발산량, 옥상녹화와 대조구인 콘크리트의 표면과 천장면의 온도, 열류량, 열수지를 측정하였다. 6일과 7일의 일최고기온은 각각 $29.4^{\circ}C$$30^{\circ}C$로 초여름 날씨로는 매우 무더웠다. 6일과 7일의 일증발산량은 각각 $2,686.1g/m^2$$3,312.8g/m^2$이었으며, 일사량이 증가함에 따라 증발산량도 증가하였다. 외부 기온이 가장 높은 시간일 때의 옥상녹화 표면과 대조구 표면의 온도차이를 비교해 보면 대조구 표면이 $5.5^{\circ}C$(6일)와 $2.2^{\circ}C$(7일) 더 높은 것으로 나타났다. 옥상녹화 천장과 대조구 천장의 온도차이를 비교해 보면 $13.1^{\circ}C$(6일)와 $14.2^{\circ}C$(7일)로 나타나 옥상녹화에 의한 실내온도저감효과가 매우 큰 것으로 나타났다. 옥상녹화와 대조구인 콘크리트 표면의 열수지를 분석한 결과, 옥상녹화는 순복사량과 잠열량이 많은 비율을 차지하고 있었으나 대조구는 현열량과 전도열량의 비율이 높았다. 즉, 옥상녹화는 한국잔디와 식재토양에 의한 증발산으로 열을 식혀주는 잠열의 비율이 높은 반면에 대조구인 콘크리트 표면은 식물이 식재되어 있지 않기 때문에 열을 식혀주는 잠열은 $0W/m^2$인 반면에 주위 온도를 증가시키는 현열량과 건물내로 전달되는 전도열량의 비율이 높아 표면과 천장면의 온도를 상승시켰다. 본 연구는 옥상녹화의 온도저감과 열수지 해석을 시도했다는 점에서는 매우 의의가 크다고 할 수 있지만, 옥상녹화는 토심과 토양배합비, 식물종, 계절 등 여러 변수에 따라 온도저감효과가 다르게 나타날 수 있기 때문에 이에 대한 후속연구가 필요하다.

두부의 전열물성 및 유효열전도도의 추정 2. 대두단상질의 고유열전도도 측정과 두부의 유효열전도도의 추정 (Thermophysical Properties of the Soybean Curd and Prediction of its Thermal Conductivity 2. The 'intrinsic' thermal conductivity of soybean protein and prediction of the thermal conductivity of soybean curd)

  • 공재열
    • 한국수산과학회지
    • /
    • 제15권3호
    • /
    • pp.219-225
    • /
    • 1982
  • 1) 2성분계 도지두부의 유효열전도도($\lambda_e$)를 $0\sim20^{\circ}C$$-5\sim-20^{\circ}C$의 범위내에서 측정한 결과, 수분함양의 증감에 따라 $\lambda_e$값도 증감했으며, 동결점이상의 온도대에서는 온도의 상승과 함께 $\lambda_e$값도 커지고, 동결점 이하에서는 온도의 하강과 함께 $\lambda_e$값이 커졌다. 2) 3성분계 두부의 유효열전도는 수분함양이 증가함에 따라 $\lambda_e$값은 커져가나 지질의 함양이 증가함에 따라 $\lambda_e$값은 작아졌다. $\lambda_e$값의 온도의존성은 2성분계도지두부의 경우와 유사했다. 3)대대단백질의 고유열전도도는 미동결상태에서 $\lambda_p=0.300[W/m{\bullet}k](0\sim-20^{\circ}C)$이고, 동결상태에서 $\lambda_p=0.488[W/m{\bullet}k](-5\sim-20^{\circ}C)$였으며 이들은 실험온도 범위내에서 온도의존성을 나타내지 않았다. 4) 대두단상질의 고유열전도도의 추정치 $\lambda_p=0.300$$\lambda_p=0.488$의 타당성을 지질의 함양이 상이한 2종류의 3성분계 두부에 적용한 결과, 실측치와 계산치는 잘 일치했다.

  • PDF

루트부 갭이 있는 양면 필릿용접 이음부의 용접잔류응력 분포 (Distribution of Welding Residual Stresses in T-joint Weld with Root Gap)

  • 방한서;김성환;김영표;이창우
    • 대한조선학회논문집
    • /
    • 제39권3호
    • /
    • pp.81-88
    • /
    • 2002
  • 용접구조물의 루트부는 외력에 의한 응력 집중에 의해 파손되기 쉽다. 따라서 구조물의 안전성 및 신뢰성 측면에서 홈 가공한 그루브 용접에 의한 완전용입 용접이 일반적으로 요구되어진다. 하지만 필릿 T-이음부 용접은 루트부의 갭과 같은 불완전 용입부를 만들어내기 쉬움에도 불구하고 홈 가공 시간 및 용접봉 소모량을 줄이기 위해 이러한 필릿용접이 자주 행해지고 있다. 따라서, 본 연구에서는 필릿 용접구조물의 플래이트(또는 플랜지)와 웨브 부분에 발생하는 용접잔류 응력과, 특히 불완전 용입에 의한 루트부 갭을 갖는 양면 T-이음부의 노치부분에 발생하는 잔류응력 분포를 해석하고자 하였다. 해석을 위해서 서브머지드 아-크 용접에 의한 단층 및 다층패스용접 모델을 선정하였으며, 열전도 및 열탄소성 이론을 고려한 유한요소 프로그램을 사용하였다.

냉동 물류 창고 내 도크시스템을 통한 에너지 손실량 분석 (Analysis of Amount of Energy Loss for a Dock System in the Cold Distribution Center)

  • 양성준;김영주;허준;김태성
    • 설비공학논문집
    • /
    • 제29권8호
    • /
    • pp.419-428
    • /
    • 2017
  • In this study, energy loss due to ventilation load in the dock system was analyzed through simulation. Also, flow generated in the dock system of the warehouse was measured using manufactured measuring devices. Numerical simulation was conducted by simulating the most common picking tasks by examining the actual working environment. Incompressible and unsteady turbulent flows were assumed, and the turbulence model was the k-e standard model. Proper grid was selected through grid dependency test. Measurement was conducted using Honeywell and Vaisala sensors, and flow and temperature inside the warehouse were measured and compared with simulation results to validate simulation. When comparing amount of loss occurring in two hours and amount of loss occurring in 15 minutes, docking time of the former was eight times longer but energy loss was 3.8 times lower. Ventilation load occurring during the initial period after opening docking system accounted for a large proportion of total ventilation load. Also, comparing the load when the dock was closed and the load when the truck was parked, ventilation load was significantly higher than load due to heat conduction from the wall. Therefore, in improving the docking system, it is effective to reduce the gap by improving compatibility of the docking system and truck, rather than wall material.

외기대류계수에 관한 실험적 연구 (Experimental Study on Coefficient of Air Convection)

  • 전상은;김진근
    • 콘크리트학회논문집
    • /
    • 제15권2호
    • /
    • pp.305-313
    • /
    • 2003
  • 본 연구에서는 외기와의 열전달을 나타내는 외기대류계수에 관한 실험을 실시하였다. 외기대류계수에 관한 기존의 모델에서 나타났던 문제점을 해결하기 위해 실험변수로 풍속외에 거푸집의 종류(목재, 철재)를 선정하였다. 실험결과를 이용하여 외기대류계수를 산정하고자 열평형 방정식을 이용한 수치해법을 사용하였으며, 이론적인 고찰을 통해 각 거푸집별로 풍속에 따른 외기대류계수의 변화를 예측할 수 있는 모델식을 제안하였다. 제안된 모델식에 의하면, 모든 경우에 풍속에 따라 외기대류계수가 증가하는 경향을 보였으나 거푸집의 사용여부나 거푸집 재료에 따라 다른 양상을 보이는 것을 알 수 있었다. 이러한 양상의 차이는 거푸집 재료의 열특성에 의해 결정되는 것으로 외기대류계수는 거푸집이 없는 경우, 철재 거푸집을 사용한 경우, 목재 거푸집을 사용한 경우의 순으로 풍속의 영향을 받는 것으로 나타났다. 제안된 모델식을 이용하면 수화열에 의한 콘크리트 구조물의 온도해석시 보다 정확한 결과를 얻을 수 있을 것으로 사료되며, 향후 이러한 열특성계수에 대한 연구가 필요할 것으로 판단된다.

Voronoi 입자기반 개별요소모델을 이용한 암석 균열의 열에 의한 미끄러짐 해석: 국제공동연구 DECOVALEX-2023 Task G(Benchmark simulation) (Voronoi Grain-Based Distinct Element Modeling of Thermally Induced Fracture Slip: DECOVALEX-2023 Task G (Benchmark Simulation))

  • 박정욱;박찬희;이창수
    • 터널과지하공간
    • /
    • 제31권6호
    • /
    • pp.593-609
    • /
    • 2021
  • 본 연구에서는 입자기반 개별요소모델(grain-based distinct element model, GBDEM)을 이용하여 결정질 암석 내 포함된 균열의 열-역학적 거동을 평가할 수 있는 수치해석기법을 제시하고 열에 의한 균열의 미끄러짐 거동을 해석하였다. 이는 DECOVALEX-2023 프로젝트 Task G의 일환으로 수행된 벤치마크 모델링 연구로, Task G는 결정질 암반 내 균열의 열-수리-역학적 복합거동을 해석하기 위한 수치해석기법을 개발하는 데에 목표가 있다. 여기에서는 Voronoi diagram을 이용하여 다면체 개별입자의 집합체로서 해석모델을 생성하고, 입자 및 입자간 접촉에서 발생하는 열-역학적 거동을 개별요소프로그램인 3DEC을 통해 해석하였다. 암석 시험편의 탄성거동을 재현하기 위하여 등가연속체 개념을 적용하여 입자와 접촉의 미시물성을 산정하였으며, 균열에 상응하는 접촉에는 Coulomb slip model을 부여하여 인장강도와 전단강도를 갖는 불연속면을 모사하였다. 경계응력과 열응력에 의한 균열의 거동을 수치적으로 모델링하였으며, 경계조건에 따라 균열의 미끄러짐이 발생하는 열-역학적 메커니즘을 정량적으로 분석하였다. 해석 결과, 본 연구에서 제시한 해석모델이 암석 내 열팽창과 열응력의 증가, 균열 응력과 변위, 경계조건의 영향 등을 합리적으로 재현하고 있음을 확인하였다. 본 연구의 해석모델은 Task G에 참여하는 국외 연구팀들과의 의견 교류와 워크숍을 통해 지속적으로 개선하는 한편, 향후 실내실험에 적용하여 타당성을 검증할 예정이다.

차세대 인공위성 전기저항제트 가스추력기의 다물리 수치모사 (MULTI-PHYSICAL SIMULATION FOR THE DESIGN OF AN ELECTRIC RESISTOJET GAS THRUSTER IN THE NEXTSAT-1)

  • 장세명;최진철;한조영;신구환
    • 한국전산유체공학회지
    • /
    • 제21권2호
    • /
    • pp.112-119
    • /
    • 2016
  • NEXTSat-1 is the next-generation small-size artificial satellite system planed by the Satellite Technology Research Center(SatTReC) in Korea Advanced Institute of Science and Technology(KAIST). For the control of attitude and transition of the orbit, the system has adopted a RHM(Resisto-jet Head Module), which has a very simple geometry with a reasonable efficiency. An axisymmetric model is devised with two coil-resistance heaters using xenon(Xe) gas, and the minimum required specific impulse is 60 seconds under the thrust more than 30 milli-Newton. To design the module, seven basic parameters should be decided: the nozzle shape, the power distribution of heater, the pressure drop of filter, the diameter of nozzle throat, the slant length and the angle of nozzle, and the size of reservoir, etc. After quasi one-dimensional analysis, a theoretical value of specific impulse is calculated, and the optima of parameters are found out from the baseline with a series of multi-physical numerical simulations based on the compressible Navier-Stokes equations for gas and the heat conduction energy equation for solid. A commercial code, COMSOL Multiphysics is used for the computation with a FEM (finite element method) based numerical scheme. The final values of design parameters indicate 5.8% better performance than those of baseline design after the verification with all the tuned parameters. The present method should be effective to reduce the time cost of trial and error in the development of RHM, the thruster of NEXTSat-1.