• Nuclear Engineering and Technology
• /
• 제55권12호
• /
• pp.4447-4464
• /
• 2023

This paper focuses on the development of a new computational model of the CNESTEN's TRIGA Mark II research reactor using the 3D continuous energy Monte-Carlo code TRIPOLI-4 (T4). This new model was developed to assess neutronic simulations and determine quantities of interest such as kinetic parameters of the reactor, control rods worth, power peaking factors and neutron flux distributions. This model is also a key tool used to accurately design new experiments in the TRIGA reactor, to analyze these experiments and to carry out sensitivity and uncertainty studies. The geometry and materials data, as part of the MCNP reference model, were used to build the T4 model. In this regard, the differences between the two models are mainly due to mathematical approaches of both codes. Indeed, the study presented in this article is divided into two parts: the first part deals with the development and the validation of the T4 model. The results obtained with the T4 model were compared to the existing MCNP reference model and to the experimental results from the Final Safety Analysis Report (FSAR). Different core configurations were investigated via simulations to test the computational model reliability in predicting the physical parameters of the reactor. As a fairly good agreement among the results was deduced, it seems reasonable to assume that the T4 model can accurately reproduce the MCNP calculated values. The second part of this study is devoted to the sensitivity and uncertainty (S/U) studies that were carried out to quantify the nuclear data uncertainty in the multiplication factor k_eff. For that purpose, the T4 model was used to calculate the sensitivity profiles of the k_eff to the nuclear data. The integrated-sensitivities were compared to the results obtained from the previous works that were carried out with MCNP and SCALE-6.2 simulation tools and differences of less than 5% were obtained for most of these quantities except for the C-graphite sensitivities. Moreover, the nuclear data uncertainties in the k_eff were derived using the COMAC-V2.1 covariance matrices library and the calculated sensitivities. The results have shown that the total nuclear data uncertainty in the k_eff is around 585 pcm using the COMAC-V2.1. This study also demonstrates that the contribution of zirconium isotopes to the nuclear data uncertainty in the k_eff is not negligible and should be taken into account when performing S/U analysis.

• Radiation Oncology Journal
• /
• 제28권3호
• /
• pp.155-165
• /
• 2010

목 적: 4-dimensional computed tomography (4DCT) 영상과 on board imaging (OBI) 및 real time position management (RPM) 장치로 매 회 치료 시마다 얻은 호흡연동 직각 kilovolt (KV) 준비 영상(gated orthogonal kilovolt setup image)을 이용해 간암 환자를 치료하는 동안 발생하는 종양 위치의 불확실성을 평가하고자 했다. 대상 및 방법: 3차원입체조형치료가 예정된 20명의 간암 환자를 대상으로 RPM과 전산화단층촬영모의치료기를 이용해 치료계획용 4DCT를 시행했다. 표적 근처에 위치한 간동맥화학색전술 후 집적된 리피오돌(lipiodol) 혹은 횡격막을 종양의 위치 변이를 측정하는 표지자로 선택했다. 표지자의 위치 차이를 이용해 온라인 분할간 및 분할중 내부 장기 변이와 움직임 진폭을 측정했다. 측정된 자료의 정량적 평가를 위해 통계 분석을 실시했다. 결 과: 20명 환자로부터 측정된 표지자의 분할간변이의 중앙값은 X (transaxial), Y (superior-inferior), Z (anterior-posterior) 축에서 각각 0.00 cm (범위, -0.50~0.90 cm), 0.00 cm (범위, -2.4~1.60 cm), 0.00 cm (범위, -1.10~0.50 cm) 였다. 4명의 환자에서 X, Y, Z축 중 하나 이상에서 0.5 cm를 초과하는 변이가 관찰되었다. 4DCT와 호흡연동 직각 준비 영상으로부터 얻은 표적의 움직임 진폭의 차이는 X, Y, Z 축에서 각각 중앙값이 -0.05 cm (범위, -0.83~0.60 cm), -0.15 cm (범위, -2.58~1.18 cm), -0.02 cm (범위, -1.37~0.59 cm) 였다. 두 영상간 표적의 움직임 진폭 차이가 1 cm를 초과하는 환자가 Y축 방향으로 3명 관찰되었으며, 0.5 cm 초과 1 cm 미만의 차이를 보이는 환자도 Y축과 Z축 방향을 합쳐 5명 관찰되었다. 분할중 표지자 위치 변이의 중앙값은 X, Y, Z축에서 각각 0.00 cm (범위, -0.30~0.40 cm), -0.03 cm (범위, -1.14~0.50 cm), 0.05 cm (범위, -0.30~0.50 cm)였으며 2명의 환자에서 1 cm를 초과하는 변이가 Y축 방향으로 관찰되었다. 결 론: 4DCT와 호흡연동 직각 KV 준비 영상으로 얻은 표지자의 분할간, 분할중 및 움직임 진폭에서 큰 변이가 관찰되었다.

• Asia Marketing Journal
• /
• 제14권1호
• /
• pp.83-98
• /
• 2012

신차와 중고차가 함께 경쟁하는 시장에서 신차의 경쟁만을 모형화한다면 가격이나 기타 프로모션 탄력성의 추정이 왜곡될 수 있다. 그러나 자동차 시장을 연구대상으로 한 선행연구의 대부분이 신차 시장의 경쟁에만 관심을 기울였던 바, 합리적인 가격결정이나 프로모션 기획에 도움을 주기에 미흡한 점이 있었다. 본 연구는 신차의 가격결정 및 프로모션 기획이 향후 중고차 시장을 통해 리바운드되어 신차 매출에 다시 영향을 미친다는 점을 반영하여 모형을 설정하였다. 즉, 서로 다른 신차간의 (혹은 서로 다른 중고차간의) 교차탄력성보다, 동일 모델의 신차와 중고차간의 교차탄력성이 높다는 가정하에 모형을 설정하였다. 방법론적으로는 네스티드 로짓(Nested Logit) 모형을 설정하여 소비자의 자동차 선택은 단계적으로 이루어진다고 가정하였다. 즉, 1단계에서 자동차 모델을 선택하고, 모델이 정해지면 2단계에서 신차와 중고차 중 선택하는 구조를 가정하였다 실증분석은 미국 전역에서 2009년 1월부터 2009년 6월까지 판매된 모든 컴팩트 카 모델 중에서 시장점유율 상위 9개 모델의 신차와 중고차를 대상으로 하였다. 실증분석을 통하여 비교 대상 모형보다 제안된 모형이 모형 적합도 측면에서 우월하고 예측타당성도 높다는 것을 보여주었다. 제안된 모형으로 부터 추정된 모수를 사용하여 몇 가지 시나리오를 상정하여 시뮬레이션을 실시한 결과, 신차(중고차)가 점유율을 높이고자 리베이트를 실시할 경우 중고차(신차)는 현재의 시장점유율을 유지하기 위해 대응 가격할인을 실시하게 되는데 할인 폭은 반대의 경우에 비해 높다는(낮다는)점을 확인하였다. 또한 시뮬레이션 결과가 시사하는 바는 신차와 중고차가 함께 경쟁하는 시장에서 IIA(Independence of Irrelevant Alternatives)모형을 적용할 경우 동일모델의 신차와 중고차간의 교차 탄력성을 과소평가하게 되어 현상유지를 위한 가격할인을 실시할 경우 적정한 수준이하로 하게 된다는 것이다.

• Journal of Radiation Protection and Research
• /
• 제33권2호
• /
• pp.61-65
• /
• 2008

목 적: 전립선암 환자의 방사선 치료 시 전립선의 위치 변화 및 비뇨기계 부작용을 줄이기 위해 방광의 체적을 일정하게 유지하도록 한 후 초음파를 이용하여 방광의 체적 변화를 측정 하였다. 대상 및 방법: 2002년 12월부터 2007년 8월 까지 이대 목동병원에서 전립선 암으로 진단받고 근치적 목적의 방사선 치료를 받은 환자 26명을 대상으로 하였다. 13명의 실험군 환자는 방광의 체적을 일정하게 유지하기 위하여 모의조사, CT 스캔, 및 치료 1시간 전에 450 cc의 물을 마시도록 하였다. 나머지 13명의 대조군 환자에 대해서는 특별한 지시를 하지 않았다. 각 환자의 CT영상에서 방광의 체적을 측정하였다. 3차원 입체 조형 치료 계획이 설계되었고 처방선량 70.2 Gy 조사할 경우 방광과 직장의 유효체적을 계산 하였다. 실험군에 속해 있는 환자에 대해 방사선 치료 전 방광의 체적을 초음파를 이용하여 6-8주 간 매주 측정하였고 방광의 체적 변화를 분석하였다. 결과: CT상에서 측정 한 방광의 평균 체적 및 표준 오차는 실험군의 경우 $283.5{\pm}114.0\;cc$ (40%) 이고 대조군의 경우 $181.2{\pm}120.1\;cc$ (66%) 이다. 실험군의 경우 CT 상의 방광 체적과 CT scan 전에 획득한 초음파 영상상의 방광 체적은 통계적으로 유의하지 않으나 상관관계가 있는 것으로 보이고 초음파를 이용하여 측정한 방광의 체적은 CT상에서의 체적의 평균 62% 이다. 실험군에 있어서 초음파 시스템을 이용하여 주간 측정한 방광의 체적의 변화는 상당히 컸으며 최소 33%에서 최대 75%까지 변화 하였다. 결론: CT 스캔 전에 일정한 양의 물을 마시도록 지시한 실험군의 환자들의 경우 대조군대비 평균 방광의 체적을 크게 유지할 수 있었다. 그러나 일정한 양의 물을 마시도록 지시하더라도 6-8주간의 치료 기간 중에 방광의 체적을 일정히 유지 하는데는 어려움이 있었다.

• 대한토목학회논문집
• /
• 제28권6B호
• /
• pp.731-741
• /
• 2008

기존의 한강 감조구간 수치모의에 대한 연구는 단면자료 획득의 어려움이나 유도의 조위자료가 없으므로 하류 경계단을 전류지점으로 하여 모의한 논문이 대부분이나 본 연구에서는 수치해도를 바탕으로 곡릉천 합류부 이후와 임진강 하구부의 단면을 가정하고 인천검조소의 조위자료를 유도지점으로 전이시켜 서해안 조위에 의한 한강하류부에서의 수리학적 거동을 해석하였다. 모형의 적용구간은 신곡수중보로부터 한강의 법적 하류단인 김포시 유도까지 총 36.8 km에 이르는 구간으로 흐름 및 혼합거동 해석은 RMA-2 모형과 서일원(2008)이 개발한 2차원 이송-분산 해석모형인 RAM4를 이용하였다. 전류지점에서의 수위 및 종횡방향 유속 실측자료와 수치모의 결과를 비교하여 흐름해석 결과를 검증하였다. 수위 관측소의 자료가 양호하고 인천검조소에서 높은 조차가 발생하는 2006년 6월 23일~25일 동안 모의한 결과 유도지점의 조위에 따라 총 5회의 역방향 흐름이 관찰되었고 최대 역류 길이는 장항 IC까지 총 32.9 km에 이르렀다. 최대 역방향 흐름의 발생 및 소멸 과정, 최고 유속선을 따른 수위 및 유속을 분석하였으며 이에 따른 비보존성 오염물질의 혼합거동을 해석하였다. 굴포천으로부터 유입된 오염물질은 하폭방향 퍼짐이 두드러지게 일어났지만 곡릉천에서 유입된 BOD는 곡릉천 합류부를 전후하여 중앙 및 좌안측 수심이 급격히 깊어지고 최대유속선이 곡릉천 방향인 우안에서 발생하고 있으므로 횡방향 혼합이 빠르게 완료되어 오염운의 반경이 상대적으로 좁은 것을 확인할 수 있었다. 또한 1차원 이송-분산방정식의 해석해를 적용해 유도지점의 염도 값을 인천검조소로부터 추산하여 한강 하류부에서의 수평 2차원 염수 혼합거동을 해석하였다.

• Journal of Animal Science and Technology
• /
• 제46권4호
• /
• pp.509-516
• /
• 2004

한우 종모우 선발을 위한 유전능력 평가에서 고려되는 형질들 중 이산형 형태로 조사되는 근내지방도의 유전변이가 추정방법에 따라 어느 정도 차이가 있는지 알아보기 위한 모의실험을 실시하였다. 모의실험 자료는 연속변량으로 간주되는 도체중 및 배장근단면적과 근내지방도의 잠재변수를 다변량 정규분포함수에서 생성하였고 근내지방도의 잠재변수를 이용하여 특정 임계값을 중심으로 순서화된 근내지방도 점수로 변화 하였따. 근내지방도의 점수 부여방법으로써 비거세우에서 조사된 근내지방도의 점수 1${\sim}$5점 사이에 정규분포에서 크게 어긋나는 분포특성을 갖도록 자료(DSI)를 생성하였고 또한 한우 거세우에서 현재 조사되고 있는 점수 1${\sim}$7점 사이에 정규 분포에 좀더 접근한 분포특성을 갖는 모의 자료(DS2)를 생성하였다. 분석방법간에 유전변이 추정의 정확도를 알아보기 위하여 1) 생성된 이들 자료를 선형으로 간주하고 다형질 혼합 선형 개체모형에서 REML 분석방법으로 유전변이를 추정하였고 2) 특정 임계치를 중심으로 잠재변수가 존재한다는 가정하에 다형질 임계 개체 혼합모형을 설정하여 Gibbs sampling 방법으로 유전변이를 추정하였다. 여기서 추정된 유전변이(유전력, 유전상관 및 잔차상관)에 대하여 모수와의 차이를 검정함으로써 편의되는 정도를 알아보았다. 모의실험은 각 자료에 대하여 10회 실시하였다. 분석결과, 근내지방도의 유전력 추정치는 DS1에서는 다형질 임계개체혼합모형을 설정하여 Gibbs sampling 방법으로 모수에 대한 사후분포의 평균으로 계산한 결과 참값과 유의적인 차이가 없는 것으로 분석되었다. 반면에 근내지방도를 선형으로 간주하고 다형질 선형 개체혼합모형에 의한 유전력 추정치는 모수보다 매우 낮은 유전력을 보였다(0.500 vs 0.315). 유전상관 추정치는 선형모형에서의 REML 방법 또는 임계모형에서의Gibbs sampling 방법에서 모두 모수와 유의적인 차이가 없는 것으로 분석되었으나 근내지방도의 잔차상관에 있어서 REML 방법으로 분석하였을 경우에 모수보다 낮게 추정되었다. 반면에 범주형 모형에서는 모수와 추정치 간에 유의적인 차이가 없는 것으로 분석되었다. 또한 7개의 범주형으로 조사된 자료(DS2)에서 이들 추정치는 DS1에서와 동일한 경향을 보였는데 그 편의 정도는 다소 적어지는 경향을 보였다. 따라서 이산형으로 조사되는 근내지방도에 대한 유전변이를 추정하기 위해서는 범주형 임계모형이 선형모형 보다 사소 정확한 추정을 할 수 있을 것으로 판단 되었다.

• 한국도로학회논문집
• /
• 제10권1호
• /
• pp.123-134
• /
• 2008

국내에 설치되고 있는 가드레일의 단부처리시설(end treatments)은 차량의 관통이나 전복을 유발할 수 있는 형상을 하고 있어 사고 시 탑승자가치명상을 입을 가능성이 매우 높다. 실제로 차량이 가드레일의 단부에 충돌하였을 경우, 가드레일이 차량내부를 관통하여 탑승자에게 직접적인 상해를 입히는 경우가 발생되고 있다. 본 연구에서는 미국과 유럽의 성능평가기준(performance evaluation criteria)에 근거하여 국내 실정에 적합한 단부처리시설에 관한 성능평가기준(안)을 제시하였다. 또한 외국의 설치사례 및 연구사례를 조사하여 단부처리시설의 mechanism을 파악하고 국내 도로환경 에 적합한 단부처리시설을 제안하였다. LS-DYNA 프로그램을 이용한 컴퓨터 시뮬레이션을 통해 제안된 단부처 리시설에 대한 성능을 검증한 결과 성능평가기준(안)에 모두 만족한 결과를 보였으며, 개발 단부처리시설에 대한 실물차량 충돌시험을 수행한 결과 성능평가기준(안)의 탑승자 보호성능, 단부처리시설의 거동, 충돌 후 차량의 거동을 모두 만족하였다. 국내 최초로 탑승자의 안전을 고려하여 개발된 단부처리시설이 고속도로나 국도에 설치된다면 단부충돌사고의 심각성을 감소시켜 도로의 안전성을 향상시킬 수 있을 것으로 사료된다.

• 대한전자공학회논문지SP
• /
• 제37권2호
• /
• pp.1-14
• /
• 2000

로젯 주사 탐색기는 적외선 유도 미사일에 장착되어 표적을 추적하는 장치이다. 단소자 검출기가 로젯 패턴의 형태로 공간을 주사함으로써 표적의 2차원 영상을 획득할 수 있다. 검출된 영상은 시계내의 위치에 따라서 형태가 변하고 대상 물체의 수가 고정되어 있지 않기 때문에 unsupervised clustering 방법을 이용하여 이들을 구분한다. 기존의 ISODATA 방식은 씨앗점(seed point)과 대상 화소간의 거리를 이용하여 clustering하기 때문에 물체의 모양이 복잡하거나 병합 및 분리 파라미터 값이 변하면 clustering 결과가 실제와 다르게 나타난다. 본 논문에서는 이러한 단점을 개선한 새로운 clustering 방법인 ALCA (Arrav Linkage Clustering Algorithm)을 제안한다. 이 방식은 화소가 저장된 메모리 번호의 연속성을 이용하여 clustering하기 때문에 초기 씨앗점과 병합 및 분리 파라미터를 필요로 하지 않는다. 따라서 대상 물체의 모양과 관계없이 clustering을 할 수 있다. 대상 물체의 clustering를 기존 방식과 비교 평가함으로써 제안된 방식의 우수성을 확인한다. 또한 제안된 ALCA을 로젯 주사 탐색기의 반대응 능력으로 이용하여 3차원 시뮬레이터상에서 추적 실험을 행한다. 기존 방식과 비교 평가를 통하여 제안된 ALCA 방식이 로젯 주사 탐색기의 반대응 능력으로서 우수한 성능을 가지고 있음을 확인한다.

• 한국지구과학회지
• /
• 제35권4호
• /
• pp.249-258
• /
• 2014

대기 중 이산화탄소 등의 농도가 지속적으로 증가하고 있음에도 최근 10여 년 동안(2002-현재) 전지구 지표 온도는 거의 답보상태에 머물러 있다. 이처럼 온실기체 강제력에도 불구하고, 지구 온난화 경향이 사라진 듯 보이는 현상을 지구 온난화 멈춤(hiatus)이라 한다. 이 연구는 HadGEM2-AO가 모의한 RCP8.5 시나리오 실험(95년간) 자료를 분석하여, 온난화 멈춤 시기의 특징을 분석하였다. 온난화 멈춤 기간을 나타내는 시계열은 동서 평균한 연직 해수 온도 분포를 EOF 분석하여 구한 두 번째 PC (PC2)로 정의하였다. PC2를 이용하여 온난화 멈춤과 엔소와의 관련성, 기후시스템의 변화 등을 분석하였다. 라니냐 지수(NINO3지수에 -1을 곱하여 정의)가 PC2를 약 11개월 앞서는 것으로 보아 라니냐 발생이 온난화 멈춤을 유도할 수 있음을 발견하였다. 또한 기후시스템의 냉각은 해수 표층의 열이 해양 내부로 침강으로 나타남을 보였다. 이는 해양의 열흡수에 의해 전지구 온도 상승률이 약화되었음을 의미한다. 온난화 멈춤 시기에 북태평양과 북반구 극지는 양의 온도 편차가 나타났으며, 열대 해양에서는 무역풍이 강화되었다.

• Journal of Radiation Protection and Research
• /
• 제31권2호
• /
• pp.97-104
• /
• 2006

최근 방사선 치료 및 진단 분야에서 선량 측정을 위하여 다양하게 사용되고 있는 MOSFET 선량계는 검출부위가 실리콘으로 이루어져 있으며 다른 검출기들과 마찬가지로 어느 정도의 에너지 의존도와 방향 의존도를 보인다. 따라서 MOSFET 선량제가 공기 중이 아닌 모의피폭체 내에서 선량 측정에 사용될 경우 낮은 에너지를 갖는 산란 광자 등 이차 광자들로 인하여 선량을 실제보다. 높게 평가하게 된다. 본 연구에서는 MOSFET 선량계의 에너지 의존도와 방향 의존도로 인하여 발생하는 오차를 보정하기 위한 선량보정인자를 몬테카를로 전산모사 기법을 이용하여 계산하였다. 먼저 사용되는 인형 모의 피폭체의 체적소 모의피폭체(Voxel Phanom)를 CT 영상을 이용하여 제작하였으며 제작된 체적소 모의 피폭체를 몬테카를로 전산코드로 구현한 후, 모의피폭체 내 각 선량계 지점에서 입사하는 광자의 에너지 및 방향별 에너지 스펙트럼을 계산하였다. 각각의 선량계 지점에서 0.662 MeV와 1.25 MeV의 광자빔을 고려하였으며 또한 MOSFET 선량계의 방향은 실리콘 베이스 방향과 에폭시 방향을 고려하였다. 주어진 선량제 지점에서의 선량보정인자는 계산된 에너지 의존도들의 중간간을 이용하여 결정하였으며 이렇게 결정된 각 선량계 지점에서 선량보정인자는 0.89-0.97 범위의 값들을 나타내었다. 본 연구결과에 따르면 MOSFET 선량계를 이용하여 인형 모의피폭체 내에 선량을 측정할 때 에너지 의존도와 방향 의존도를 고려하지 않을 경우 측정 위치에 따라 $3{\sim}11%$ 정도의 측정오차가 발생할 수 있다. 그러므로 인형 모의피폭체 내의 선량을 정확하게 측정하기 위해서는 선량보정인자를 각 선량계에 필히 적용해주어야 한다.