• 한국항공우주학회지
• /
• 제31권10호
• /
• pp.10-20
• /
• 2003

본 논문에서는 임의의 방향각을 가지고 적층된 평직 복합재료에 대하여 미세구조를 모델링하고 수치실험을 통하여 등가물성치를 계산하였다. 층간 상호이동에 의한 평직복합재료의 섬유다발의 배열의 경우의 수는 무한대에 이르므로 여기서는 몬테 카를로법을 이용하여 지정된 경우의 수에 대하여 모델링과 해석을 수행하는 전략을 채택하였다. 또한 평직복합재료의 미세구조를 효과적으로 고려하고 모델링 및 계산시간을 절감하기 위하여 수정된 분할영역 적분에 의한 단일변위장 마크로요소를 사용하였다. 계산결과, 평직복합재료 적층판의 등가물성치는 적층간에 따라서 편차가 크게 나타날 수 있음을 확인하였다. 또한 고전적층판 이론으로 계산한 값과 비교하였을 때 $0^{\circ}$, $45^{\circ}$ 이외의 적층각에서 값의 차이가 크게 났는데, 이로부터 평직복합재료 구조물의 경우 고전적층판이론으로 계산한 등가물성치를 사용할 때 주의가 필요함을 알 수 있었다.

• 대한화학회지
• /
• 제42권4호
• /
• pp.398-410
• /
• 1998

Boltzmann 방정식의 비선형 해법으로서 cumulant 모멘트 방법을 연구하였으며, Maxwell 분자모형 단원자분자 기체계의 정상충격파 문제에 대하여 적용하였다. 모멘트 방정식의 해는 Maxwell-Ikenberry-Truesdell(MIT) 반복법을 사용하였다. 원래의 MIT 반복법은 초기값을 평형분포함수로부터 구하지만, 본 연구에서는 반복계산의 초기값을 Mott-Smith의 두방식(bimodal)함수로부터 구하였다. 모멘트 계산은 2차 반복단계까지 수행하였으며, 강한 충격파에 대한 밀도, 온도, stress, heat flux 등의 윤곽과 충격파의 두께, 그리고 마하수 1.4 미만의 약한 충격파의 두께를 계산하였다. 1차 반복계산에서 충격파 윤곽에 대한 간단한 형태의 해석적 표현을 얻었으며, 이로부터 도출한 약한 충격파 두께에 대한 극한법칙은 Navier-Stokes 이론과 정확히 일치한다. 2차 반복계산에 의한 결과는 강한 충격파의 윤곽곡선 및 충격파 두께가 Monte Carlo 문헌값과 정량적으로 일치함을 보인다.

• 한국수학교육학회지시리즈E:수학교육논문집
• /
• 제24권1호
• /
• pp.177-193
• /
• 2010

본 연구는 고등학교 통계 영역의 확률분포에 제시되어 있는 정규분포를 이항분포에서 정규분포로의 근사, 정규분포곡선의 탐구, Monte Carlo 방법에 의한 정규분포곡선의 넓이 탐구, 정규분포곡선의 선형변환, 그리고 여러 형태의 정규분포곡선 탐구 등의 내용을 중심으로 CAS 계산기를 활용하여 탐구해보고자 한다. CAS 계산기의 도구적 기능인 사소화, 실험, 시각화, 집중의 측면에서 볼 때 지필로서는 교육과정에 제시된 확률분포의 목표를 달성하기 불가능하다고 판단된다. 따라서 본 연구에서는 CAS 계산기를 활용하여 정규분포곡선의 다양한 성질을 탐구하고 이러한 과정과 결과로부터 정규분포곡선에 대한 교수학적 시사점을 도출하고자 한다.

• 한국의학물리학회지:의학물리
• /
• 제15권2호
• /
• pp.100-104
• /
• 2004

선량 측정기의 공간적인 반응특성 때문에 나타나는 detector의 크기효과는 임상적인 선량측정을 부정확하게 만드는 중요한 원인이기에 많은 연구의 대상이 되어왔다. 관례적으로 detector response kernel은 detector 자체의 크기가 측정한 방사선의 선량분포에 대해 미친 영향에 대한 정보를 포함하고 있다. 이 kernel에 대해 다양한 수학적 모델들이 제안되었고 실험적으로 이론적으로 연구되어왔다. 이 논문은 convolution이론과 Monte Carlo simulation만을 이용하여 detector의 kernel을 결정하는 방법을 제시한다. 이 수치해석적인 방법을 사용하여 물 phantom에 잠긴 Farmer형 ion chamber의 detector response kernel을 계산하였다. 계산된 kernel은 기존의 parabolic 모델의 특성과 Gaussian 모델의 특성을 동시에 나타내고 있다. 이 kernel과 deconvolution 방법을 사용하여 측정된 6MV, 10${\times}$10 $\textrm{cm}^2$, 0.5${\times}$10 $\textrm{cm}^2$ 광자선으로부터 크기효과를 제거하였다. 크기효과가 제거된 방사선의 선량분포는 꼬리부분을 제외하고는 film이나 pin-point ion chamber에 의해 측정된 결과와 유사한 선량분포를 나타냈다.

• 대한의용생체공학회:의공학회지
• /
• 제25권2호
• /
• pp.89-95
• /
• 2004

세기조절방사선치료(IMRT)에서는 일반적인 방사선 치료에서 사용되는 조사면에 비해 비교적 작은 크기의 빔조각(beamlet)을 사용하여 방사선의 세기를 조절하는 새로운 치료법으로 이에 대한 비균질 효과는 많은 연구가 필요하다. 우리는 기하학적으로 일정한 비균질 팬텀들에서 몬테카를로 시뮬레이션에 의한 선량값을 라디오크로믹 필름에 의한 선량값과 회선/중첩 방법에 의한 선량 계산 값과 서로 비교하였다. 몬테 카를로 모사를 위하여 EGS4 코드 기반의 BEAM 코드를 사용하였으며 이를 이용하여 Varian 2300C/D 선형가속기의 두부를 호사하였다. 측정과 모사에 사용된 조사면은 1${\times}$1$\textrm{cm}^2$, 2${\times}$2$\textrm{cm}^2$, 그리고 5${\times}$5$\textrm{cm}^2$이었다. 또한 팬텀의 물질은 솔리드 워터, 폐 등가 물질, 뼈 등가 물질을 사용하여 세 경우의 비극질 팬텀들을 설정하여 방사선을 조사하였다. 회선/중첩 방법과 몬테 카를로 방법에 의한 선량 계산치는 광자 측면선량의 경우 $\pm$1 mm, 깊이선량의 경우 $\pm$2% 이내로 선량측정치와 잘 일치함을 볼 수 있었다. 결론적으로 회선/중첩 방법과 몬테 카를로 방법이 소조사면에서도 필름 측정 데이터와 잘 일치함을 확인할 수 있었다.

• 지질공학
• /
• 제27권4호
• /
• pp.429-438
• /
• 2017

지질학, 자원공학의 다양한 분야에서 적용되는 중성자검층의 활용성과 측정값의 신뢰성을 높이기 위하여 Monte Carlo 알고리듬에 기초한 MCNP 모델링을 수행하였다. 본 연구에서는 중성자 존데와 지층의 수치모형화를 기본으로 존데 고유의 교정곡선과 MCNP 모델링으로 계산한 교정곡선과의 비교를 통해 모델링의 적정성을 확인하고 암상변화, 공극 유체 특성, 시추공 지름 변화, 케이싱 영향, 공내수 영향을 모델링 결과를 이용하여 고찰하였다. 모델링 결과, 암상 변화에 따른 중성자 계수율 비율의 변화를 정량적으로 파악하였다. 시추공 지름이 존데와 비슷한 3인치의 경우, 지름이 큰 경우보다 계수율의 비가 예상보다 높게 나타났다. 이와 같은 결과는 공내수 영향이 작은 영향으로 해석되었다. 케이싱 내에서의 반응과 나공에서의 반응을 비교할 때 전반적으로 차이가 작았으며 특히 지층의 공극률이 증가하면 케이싱 영향이 감소하여 구분이 어려웠다. 지하수위 상부에 대한 모델링 결과는 지하수위 하부와는 반대의 경향을 나타냈으며 지하수위 파악에도 정성적으로 이용이 가능할 것으로 예상된다. 다양한 시추공 환경에 대한 모델링 결과는 중성자검층 현장자료의 자료처리와 해석이 유용하게 이용될 것으로 예상된다.

• 한국의학물리학회지:의학물리
• /
• 제16권2호
• /
• pp.53-61
• /
• 2005

본 연구에서는 kV X-선에 대하여 원통형전리함의 선질인자를 몬테칼로 계산과 측정으로 결정하였다. 대상 X-선은 KFDA에 설치되어 있으며 에너지스펙트럼과 선질이 명시된 $60\~300kV$ 범위의 X-선(ISO-4037)이었으며, 전리함은 PTW N23333과 N30001이었다. 이론적 계산과 몬테칼로 계산으로 구한 공기커마와 공동의 흡수선량으로부터 IAEA프로토콜의 $R_{\mu}$와 $R_{Q,Q_{0}}$인자를 결정하였다. 두 전리함의 반응 범위는 $80\~250kV$ 영역에서 모두 $\pm3\~4\%$를 보였으나, N30001 전리함의 경우에 $110\~300kV$ 영역에서 $\pm1.7\%$로서 평탄한 특성을 보였다. 본 연구로부터 kV X-선에 대한 선량평가에 대하여 고찰할 수 있었으며, 프로토콜의 적용을 위하여 에너지의존성의 평가가 중요함을 알 수 있었다.

• 한국의학물리학회지:의학물리
• /
• 제15권2호
• /
• pp.77-83
• /
• 2004

모델 기반으로 광자선 선량을 계산하는 전산화치료계획장치를 취역검사(commissioning)하기 위해서는 많은 파라미터들을 조절해가면서 측정된 심부선량 및 가로방향선량분포 등을 맞추어야한다. 우리는 몬테카를로 전산 모사를 이용하여 Pinnacle$^3$ 시스템의 취역검사에 필요한 광자선의 에너지 스펙트럼, 오염 전자(contaminant electron), 축외선질연화(off-axis softening) 및 입자 유량 증가 등을 기술하는 파라미터들을 구하였다. 몬테칼로 계산을 통해 실험으로는 측정이 쉽지 않은 이러한 양들의 변화량을 알 수 있었으나 축외선질연화 및 입자 유량 증가 변수의 경우에는 Pinnacle$^3$ 시스템을 이용한 계산과 측정값에 불일치를 보였다. 전산화치료계획장치의 취역검사에 몬테카를로 방식으로부터 얻은 파라미터 값을 그대로 이용하는 문제는 추가 연구가 필요하다.

• 한국의학물리학회지:의학물리
• /
• 제26권1호
• /
• pp.36-41
• /
• 2015

정위체부방사선수술(SBRT) 환자의 선량분포를 계산하기 위해 전산모사 방식을 이용한 응용프로그램을 개발 하였다. 본 소프트웨어는 최근 이용이 활발하게 증가하고 있는 Geant4를 기반으로 개발 하였다. 환자에 조사하기 위한 광자선 스펙트럼은 이전 연구에서 구한 선형가속기 스펙트럼 자료를 사용하였다. 치료계획시스템과 유사한 조사면을 구현하기 위하여 PrimaryGeneratorAction 클래스에서 MLC 조사면 형태를 반영하도록 하였다. 본 연구에서는 8개 조사면에 대한 계산을 수행하였으며 이 때 갠트리의 각도는 PrimaryGeneratorAction 클래스에서 회전 매트릭스를 사용하여 선원의 위치를 변경하는 방법을 사용하였다. 환자에 대한 물질 자료는 CT의 dicom 파일에서 픽셀 크기, 매트릭스 크기 등의 정보와 픽셀의 HU를 밀도로 변환한 파일을 생성한 다음 이 파일을 이용 환자의 모델링에 이용 하였다. 환자의 물질 구성과 기하학적 자료의 입력에 있어 EGSnrc 코드와의 비교를 통하여 계산의 효율성을 비교하였다.

• Nuclear Engineering and Technology
• /
• 제14권4호
• /
• pp.196-203
• /
• 1982

본 논문에서는 고리 1호기의 기사용 핵연료 집합체를 수송하기 위한 Cask를 설계하였다. 이를 위하여 고리 1호기의 기사용 핵연료 집합체로부터 방출되는 감마선과 중성자를 계산하여 MORSE 및 ANISN전산 코드로써 차폐 계산을 수행하였다. 그 결과, 9개의 집합체를 동시에 수송할 수 있는 Steel Cask가 가장 적합하다는 것을 밝혔다. 이 Steel Cask에 대한 안전성을 평가하기 위하여 연료봉의 중심 온도와 복재온도를 계산하여 핵연료의 용융점보다 훨씬 낮음을 증명하였다. 또한 KENO와 MORSE전산 코드를 사용하여 임계도 계산을 수행하여 미임계 상태임을 증명하였다. 이로써 9개의 기사용 핵연료 집합체를 동시에 수송할 수 있는 Steel Cask를 간단히 설계하였다.