• 한국전자파학회논문지
• /
• 제18권7호
• /
• pp.819-827
• /
• 2007

본 논문은 논에 대한 전파 산란 특성을 측정하고, 수치적 계산 결과와 비교함으로써 벼의 수치적 모델링에 대한 전파 산란 알고리즘을 제시하고자 한다. 임피던스 경계 조건과 dyadic 그린함수로부터 임피던스 표면 위(수면)의 손실 있는 유전체(벼)에 대한 적분 방정식을 유도하였고, 모멘트 법을 이용하여 유전체의 체적 전류를 계산하였다. 또한, 몬테 카를로 방법을 적용하여 입사 각도 및 편파에 따른 후방 산란 계수를 수치적으로 계산하였다. 1.85 GHz의 측정 시스템을 이용하여 논의 후방 산란 계수를 측정하였으며, 본 논문에서 제시한 알고리즘의 계산 결과와 비교, 검증하였다.

• Nuclear Engineering and Technology
• /
• 제19권1호
• /
• pp.34-41
• /
• 1987

본 논문은 부하추종운전노심의 1차원 해석방법과 3차원 해석방법을 비교하기 위한 것이다. 이 목적으로 본 논문에서는 유한차분법에 의거한 1차원 HLOFO 코드와 수정형 Borresen 소격모형에 의한 3차원 CMSNAC 코드를 개발했다. CMSNAC 코드는 100-50-100% 출력변화를 수반하는 12-3-6-3 부하추종운전에 대한 고리 1호기 노심의 3차원 출려분포와 반응도 변수를 계산하는데 이용했다. 3차원 계산결과는 중성자속.체적 가증법으로 구한 단면적 및 버클링 입력에 의거한 1차원적 HLOFO 계산결과와 비교되었으며, 이 비교로부터 1차원 HLOFO 계산이 정착도 면에서 다소 뒤떨어지지만 3차원 계산에 가까운 결과를 얻을 수 있다는 것을 보였다.

• Journal of Radiation Protection and Research
• /
• 제26권2호
• /
• pp.93-99
• /
• 2001

해부학적으로 단순한 수학적인형팬텀의 한계를 극복하기 위한 voxel 머리팬텀을 제작하고 BNCT(Boron Neutron Capture Therapy) 시행 시 선량분포를 계산하였다. 일반목적 몬테칼로 코드인 MCNP4B의 반복구조 알고리즘을 이용하여 voxel 몬테칼로 계산체계를 수립하였고 두 가지 물질로 구성된 예시적 voxel 팬텀과 기하체조합팬텀의 계산값 비교를 통해 계산체계를 검증하였다. 미국 NLM(National Library of Medicine)에서 제공하는 VHP man 인체단층사진에 대한 분할 및 색인작업을 통해 voxel 머리팬텀을 제작하여 AP 및 PA 방향에서 입사하는 넓고 평행한 광자 및 중성자빔에 대한 선량값을 MIRD 팬텀의 계산값과 비교한 결과 중성자빔 AP 방향조사 시 MIRD 팬텀에서는 볼 수 없는 안구로 인한 중성자 감쇠현상을 확인할 수 있었다. 3차원 정밀계산이 필요한 BNCT 시술시 선량분포계산을 위해 뇌 중앙에 직경 5cm의 구형 뇌종양 체적을 정의하고 뇌와 종양의 붕소 함량을 조정하여 10keV 및 40keV 상부입사 중성자에 의한 장기별 흡수선량을 계산한 결과 종양에 $30{\mu}g/g$, 정상세포에 $3{\mu}g/g$의 붕소를 주입한 경우 붕소함량이 없을 때에 비해 2배 가량 큰 선량을 보였다. 본 연구를 통해 voxel몬테칼로기법을 이용한 선량평가체계를 수립하였고 정밀한 선량계산을 필요로 하는 치료방사선분야 선량계산에 실제 인체에 가까운 voxel팬텀의 응용가능성을 제시하였다.

• 한국의학물리학회지:의학물리
• /
• 제25권3호
• /
• pp.139-142
• /
• 2014

본 연구의 목적은 흡수선량 재구성, 방사선 치료간의 재구성된 선량의 등록, 선량-체적 히스토그램의 생산등을 수행하는 선량 재구성의 임상적 응용시스템을 만들고 그것을 변형된 전립선 팬텀에 적용하는 것이다. 이를 위해 변형가능한 전립선 팬텀을 20 cm 깊이와 40 cm너비의 물팬텀에 집에 넣었다. 이것의 영상을 얻고, 전립선, 정낭 및 항문의 윤곽을 그렸다. 동일 평면에서 네개의 조사문을 이용하여 세기 변조계획을 세웠다.항문에 20 ml의 물풍선을 삽입하여 장기를 변형시켰다. 영상을 다시 획득하여 위 장기의 윤곽을 그렸다. XVMC몬테칼로 코드를 사용하여 두 팬텀및 EPID내에서 선량반응 인자를 계산하였다. 세기변조계획에서 얻어진 방사선을 두팬텀에 조사하여 EPID에서 적분형 영상을 얻었다. Demons 방법을 사용하여 변형된 팬텀을 변형전 팬텀에 등록시켰다. 이를 통해 단위체적별 위치변이 정보를 얻었고 이를 이용해 두 팬텀의 재 구성된 선량을 합하여 변형전 팬텀에 생산해 냈다. 순방향으로 계산된 치료계획 선량을 합산된 재구성된 선량과 비교하였다. 200 cGy에서 전립선과 정낭이 받든 체적은 차이를 거의 보이지 않았으나, 210 cGy 이상에서는 3%가량 차이를 보였다. 항문에서는 150-200 cGy영역에서 재구성된 선량에 의하여 받은 체적은 치료 계획과 비교하여 3% 이상 적었다. 본 연구를 통하여 선량 재구성의 임상적 응용시스템이 성공적으로 만들어 졌다. 변형된 전립선 팬텀에 적용되어 작지 않은 선량의 차이를 목표장기와 보호 장기에 보였다.

• 한국광물학회지
• /
• 제18권4호통권46호
• /
• pp.249-257
• /
• 2005

석류석은 하부맨틀의 최상부인 약 660 km깊이까지 올리빈, 휘석과 함께 주요한 구성광물 중의 하나이다. 석류석은 약 660 km를 지나 하부맨틀에서 페롭스카이트와 코런덤으로 상변이를 하는 것으로 알려져 있으나, 실험방법 및 상변이 깊이와 상변이 경로에 대해서 아직까지 논쟁의 대상이 되고 있다. 실험은 천연산 알만딘(($Fe_{2.52}Ca_{0.21}Mg_{0.18}Mn_{0.12})Al_{2.23}Si_{2.97}O_{12}$)에 마오-벨 타입의 다이아몬드 앤빌기기를 이용하여 실온에서 62 GPa까지 압축실험을 시행하여 체적탄성계수를 결정하였다. 실험결과는 다음과 같다 : 62 GPa에서 격자상수 = $10.775\;{\AA}$, 체적 = $1251.16\;{\AA}^{3}$, X-선밀도 = $5.265\;g/cm^{3}$, 버치-머내한 상태방정식을 이용하여 계산한 체적탄성계수는 156 GPa이다(이때 $K_{0}\;'$는 4로 가정함). 본 연구는 포항가속기연구소의 방사광을 이용하여 시행한 국내 첫 고압실험결과이다.

• 대한토목학회논문집
• /
• 제29권5B호
• /
• pp.429-439
• /
• 2009

댐 붕괴로 인한 극한홍수가 발생하였을 경우, 홍수경보에 대한 대응시간은 일반적인 홍수의 경우보다 훨씬 짧다. 수치모형은 홍수파의 전파양상을 예측하고, 범람지역, 홍수파 도달시간 그리고 침수심 등에 관한 정보를 제공하는데 있어 강력한 도구가 될 수 있다. 그러나 댐 붕괴로 인한 홍수파의 전파는 불연속 흐름이나 마른하도의 전파를 포함하고 있으므로, 수학적으로 표현하기 어려운 경우가 많다. 그럼에도 불구하고 최근에 유한체적기법을 이용하여 댐 붕괴로 인한 홍수범람을 모의하기 위한 수치모형의 개발이 많이 이루어졌다. 유한체적기법은 적분보존형 방정식을 기본으로 하고 있으므로, 불연속 흐름이나 충격파의 해석에 용이하다. 따라서, 본 연구에서는 2차원 보존형 천수방정식의 해석을 위해 유한체적기법과 Riemann 근사해법을 이용한 수치모형을 개발하였다. 그리고 예측단계와 수정단계에서 연속방정식과 운동량 방정식의 보존변수 재구성을 위해 수면경사법과 연계한 MUSCL 기법을 적용하여 시간과 공간에서 2차정확도를 얻었다. 개발한 유한체적모형을 2차원 부분적 댐 붕괴 해석 및 삼각형 융기를 가진 하도에 대한 댐 붕괴 해석에 적용하고, 적용결과를 실험자료 및 기존 연구자의 계산결과와 비교하여 개발모형을 검증하였다.

• 한국의학물리학회지:의학물리
• /
• 제19권4호
• /
• pp.263-268
• /
• 2008

방사선치료 시 환자 호흡에 의한 종양의 움직임 영향을 줄이기 위해 사용하는 호흡연동 방사선치료 방법을 폐암 환자의 방사선수술에 적용하였을 경우와 종양의 움직임 영역을 모두 고려하여 계획용표적체적(PTV)을 설정한 일반적인 방사선 수술 경우를 비교, 평가하고자, 각 치료계획에서 계산한 종양 주변의 중요 장기에 조사되는 선량분포의 변화를 분석하고, 그 영향을 고찰하였다. 총 5명의 방사선수술 대상 비소세포성 폐암 환자의 4차원 전산화단층촬영(CT) 모의치료 영상을 획득한 후, 호흡주기 내 50% 위상에서 구성한 CT 영상을 기반으로 호흡연동 방사선치료 방법의 방사선수술 계획을 수립하고, 호흡주기 내 모든 위상에서의 CT 영상들에서 생성한 내부표적체적(ITV)을 기반으로 한 방사선수술 계획을 수립하여 각 계획에서 계산된 중요장기들의 선량체적히스토그램(DVH)을 비교, 분석하였다. 호흡연동 방법을 적용하지 않은 ITV를 기반으로 한 치료계획과 비교한 결과, 호흡연동 방사선치료 방법을 적용하였을 경우의 중요장기 평균선량이 기관지는 평균 $16.88{\pm}9.97%$, 척수는 평균 $34.13{\pm}19.15%$, 흉벽은 평균 $28.42{\pm}18.49%$, 폐의 경우는 평균 $32.48{\pm}16.66%$의 상대적인 감소를 보여, 폐암 환자의 방사선수술 시 호흡연동 방사선치료 방법 적용의 유용성 및 타당성을 확인할 수 있었다.

• 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%$ 정도의 측정오차가 발생할 수 있다. 그러므로 인형 모의피폭체 내의 선량을 정확하게 측정하기 위해서는 선량보정인자를 각 선량계에 필히 적용해주어야 한다.

• Nuclear Engineering and Technology
• /
• 제16권2호
• /
• pp.70-79
• /
• 1984

대형 고속증식로용 핵계산 체제(KAERI-26군 단면적 library/1DX/2DB)를 이용하여 SUPER-PHE-NIX I 평형 노심의 초기상태에서의 diagrid 및 Pad팽창에 대한 반응도 계수를 계산, 검토하였다. 노심은 R-Z등가 model로 묘사하고, 2차원 다군 확산 이론 전산코드인 2DB를 사용하여 임계도를 계산하였다. 기초계산으로서 반경방향 및 축방향 균일 괭창에 대한 반응도 계산과 노심구성물질의 원자수 밀도 변화와 노심체적 변화에 대한 반응도 변화량계산을 수행하였다. 균일 팽창으로 고려한 diagrid팽창에 대한 반응도 계수는 -0.553pcm/mil로 계산되었다. 한편 반응도의 온도계수는 -1.0766pcm/$^{\circ}C$로 환산되어 프랑스 발표치 -1.09pcm/$^{\circ}C$와 1.2%오차내로 일치하였다. Diagrid 팽창효과 졔산방법을 활용하여 노심의 불균일 팽창을 유발하는 pad팽창에 대한 반응도 계수 계산을 수행한 바 매우 유용함을 알았다. Pad팽창에 대한 반응도 계수는 평균팽창 model의 경우 -0.2743pcm/mi1, pancake집적 model의 경우 -0.2786pcm/mi1로 각각 계산되었다. 또한 자유팽창 노심에서의 온도변화에 따른 pad 팽창에 대한 반응도 계수는 각 model에 대하여 -0.5766pcm/$^{\circ}C$, -0.5858pcm/$^{\circ}C$로 계산되어 프랑스 발표치 -0.574pcm/$^{\circ}C$와 2% 오차내로 일치하였다.

• 한국추진공학회지
• /
• 제3권4호
• /
• pp.58-66
• /
• 1999

초음속 이젝터 유동특성은 압축성, 비정상성, 충격파 둥으로 인하여 충분히 알려져 있지 않다. 종래의 이론적/실험적 연구결과들은 대부분 1차원 가정 하에서 얻어진 것들이며, 현재까지 수치계산법을 이용하여 이젝터 내부 유동장을 해석한 실례가 많지 않다. 뿐만 아니라 기존의 수치계산 결과들은 2차유동의 질량유량이 매우 작거나 없는 단순한 유동장에 대한 것들이었다. 본 연구에서는 초음속 이젝터-펌프 유동장을 수치적으로 해석하기 위하여 축대칭 수치계산 모델을 이용하였다. 수치계산은 축대칭 압축성 Navier-Stokes 방정식에 유한체적법을 적용하였으며, 표준형 $\kappa$ -$\varepsilon$ 난류모델을 이용하였고, 2차유동의 질량유량을 변화시켜 1차유동과 2차유동의 질량비의 변화에 따른 이젝터 내부 유동장의 변화를 조사하였다. 그 결과, 2차목을 가지지 않는 일정 단면적의 혼합부를 가지는 이젝터의 경우, 이젝터 내부의 유동장 특성은 2차유동의 질량유량의 변화와 거의 무관하게 나타났다. 그러나 2차목을 가지는 경우는 2차유동의 질량유량의 증가에 따른 노즐출구 주위에서 배압으로 인하여 이젝터 내부 유동장이 크게 변화하였다.