• 한국융합학회논문지
• /
• 제9권3호
• /
• pp.257-263
• /
• 2018

고주파 열치료법은 1MHz 이상의 RF 고주파 전류를 전극을 통해 종양조직으로 공급, 가열하여 종양조직의 온도를 $42.0^{\circ}C$이상으로 상승시켜 열괴사시키는 치료법으로 알려져 있다. 인체 내에서 전자기에너지의 흡수와 전달을 위한 수학적으로 모델링과 생물학적인 신체조직의 온도 필드 분포의 해석과 평가를 통해 생체조직을 구성하는 분자들이 진동하면서 서로 마찰되어 열에너지로 전환되는 과정을 분석할 수 있다. 본 논문에서는 3차원 모델의 기하학적 형상의 성인남자 표준모델을 토대로 인체모형을 설정하고, 계산은 유한체적법 코드를 활용하였다. 인체에 전극을 장착하여 외부에서 유도되는 자기장을 모사한 쥴열이 공급되는 것으로 가정하였고 이에 대한 온도분포를 0-1,200초 동안 해석하였다. 시뮬레이션의 결과, 전달된 에너지는 전극의 가장자리로부터 전극의 안쪽으로, 피부 표면으로부터 피하층으로 점진적으로 침투되어 폐종양세포에 전달되어 종양이 열괴사하는 과정을 확인할 수 있었다.

• 한국산학기술학회논문지
• /
• 제18권8호
• /
• pp.25-33
• /
• 2017

본 연구에서는 간편하고 빠른 사출성형을 위하여 한쪽이 평면인 비구면 형상을 갖는 단일 fly-eye 렌즈형태의 2차 렌즈를 사용하는 MR-16 램프 시리즈의 설계를 보고하고자 한다. 인테리어 조명용 할로겐 램프를 대체하기 위해 효율이 높고 좁은 배광각을 갖도록 MR-16 램프용 광학계를 설계하여 이를 제작하고 성능을 확인하였다. 광학설계용 소프트웨어는 ORA사의 CODE V와 LightTools를 사용하였다. 그 결과, 제작된 MR-16 램프의 배광각은 $22.4^{\circ}$, $31.1^{\circ}$, $37.3^{\circ}$, $59.9^{\circ}$였으며, 광효율은 각각 76.5 lm/W, 75.2 lm/W, 72.0 lm/W, 77.8 lm/W 이고, 복사 균제도 81%에서 발산각은 $3^{\circ}$, $15^{\circ}$, $22^{\circ}$, $49^{\circ}$이다. 그리고 제작된 MR-16 램프에 확산시트를 2차 렌즈와 LED사이에 확산 시트를 적용하여 황색 띠를 제거한 후의 배광각은 각각 $20.8^{\circ}$, $31.5^{\circ}$, $37.8^{\circ}$, $68.7^{\circ}$로 측정되었다. 실험적으로 2차 렌즈를 같은 렌즈 4개를 붙이는 형식으로 fly-eye 형태로 함으로써 설계 및 제작이 용이하며 배광각과 균일한 조도 분포를 갖는데 있어서 문제가 없다는 것을 알 수 있었다.

• 대한원격탐사학회지
• /
• 제38권6_1호
• /
• pp.1385-1394
• /
• 2022

본 연구에서는 MODerate resolution Imaging Spectroradiometer (MODIS) 자료를 활용하여 20년(2002-2021)의 중분류 토지피복별 알베도 변화를 분석하고, Visible Infrared Imaging Radiometer Suite (VIIRS)의 10년(2012-2021) 자료를 활용해 MODIS 자료와의 차이를 분석하였다. MODIS와 VIIRS의 알베도 자료는 Bidirectional Reflectance Distribution Function (BRDF) 모델을 활용해 생산된 Sinusoidal Tile Grid 기반 500 m 공간해상도의 일단위 알베도 자료 MCD43A3와 VNP43IA3를 우리나라 범위에 대하여 구축하였다. Python 3.9 기반으로 작성된 코드를 활용하여 Reprojection을 하였으며, Resampling method는 Nearest neighbor를 적용하였다. 알베도 분석에는 단파 영역(Shortwave)의 White sky albedo와 Black sky albedo를 활용하였다. MODIS 자료를 활용한 20년의 알베도 분석 결과, 모든 토지이용에서 알베도가 상승하는 경향이 나타났다. 2000년대(2002-2011)에 비해 2010년대(2012-2021)의 평균 알베도가 산림 지역에서 0.0027의 가장 큰 상승값을 보였고, 그 다음으로 초지가 0.0024의 상승값을 보였다. VIIRS와 MODIS의 알베도를 비교한 결과, VIIRS의 알베도가 0.001에서 0.1 만큼 더 큰 것으로 나타났으며, 이는 영상의 촬영시기에 따른 지표면 반사도와 센서의 특성 차이에 의한 것으로 판단된다.

• 한국의학물리학회지:의학물리
• /
• 제10권1호
• /
• pp.1-7
• /
• 1999

본 연구에서는 $\beta$선 방출 동위원소들 중에서 $^{32}$P으로부터 방출되는 $\beta$선에 의한 홉수선량 분포를 방사선원의 형태와 기하학적인 조건을 달리하여 컴퓨터를 이용한 모사실험을 통해 예측함으로써 balloon catheter 및 radioactive stent의 이용과 관련된 정보를 얻고자한다. $^{32}$p로부터 방출된 $\beta$선이 인체내에서 에너지를 전달하는 과정에 대한 모사실험은 EGS4 code system 을 이용하여 수행되었다. 인체내의 방사선 흡수선량은 선원의 형태와 위치를 고려하여 축방향과 반경방향으로 등간격으로 나누어 각 격자에서 계산되었다. $^{32}$P 에서 방출되는 $\beta$선 에너지는 Coulomb 포텐셜에 대한 Dirac방정식의 해를 이용하여 계산된 $\beta$선 스펙트럼의 결과를 사용하여 무작위로 선정되었다. 체적 선원과 표면선원에서 시료 표면으로부터 반경방향으로 깊이 0.5 mm내에 있는 표적체에서의 선량률은 각각 12.133 cGy/s per GBq (0.449 cGy/s per mCi, uncertainty: 1.51%)와 24.732 cGy/s per GBq (0.915 cGy/s per mCi, uncertainty: 1.01 %)이다. 선량률은 시료표면으로부터 축방향과 반경방향으로의 거리에 따라 감소한다. 본 연구 결과를 근거로하여 balloon catheter 및 radioactive stent에 $^{32}$P 핵종을 사용할 때 치료선량을 20 Gy로 할 경우 치료에 적합한 초기 방사능량은 각각 29.69 mCi(치료시간을 3분으로 제한할 때) 와 1.2278 $\mu$Ci (영구삽입)로 계산되었다. 또한 원통형 체적선원과 표면선원에 대하여 초기방사능의 크기를 1 mCi/ml의 방사능 체적 밀도와 0.1 mCi/$cm^2$의 방사능 면 밀도로 나타내었을 때 각 표적체에서의 흡수선량률을 계산하였다. 통일한 값의 방사능 체적 밀도와 방사능 면 밀도는 크기가 다른 모델에 대해서 비슷한 크기의 홉수선량을 유도하므로 $^{32}$p 방사선원의 초기 방사능 체적 밀도와 초기 방사능 면 밀도를 알고 있을 때 본 연구의 계산 결과를 이용하면 직경과 길이가 다른 $^{32}$P 핵종의 원통형 모델 주위의 홉수선량 분포를 쉽게 계산할 수 있다.

• 대한치과교정학회지
• /
• 제30권2호
• /
• pp.127-142
• /
• 2000

임상에서 다양하게 사용되는 Multiloop Edgewise Arch Wire(MEAW)는 구치부 치아의 치축 직립을 위해 II급 고무줄이나 수직고무줄과 함께 사용된다. 이때 tip back bend의 정도는 각 환자의 구치치축의 근심경사의 정도에 따라 가감이 결정되나 이 장치가 실제 치근막에서 어떤 효과를 나타내는지 정확히 알려진 바는 없다. 이에 저자는 MEAW의 역학적 효과를 알아보기 위해 tip back bend의 정도를 변화시킨 후 II급 고무줄, 수직고무줄을 적용시켜 유한요소법을 통해 치근막에서 발생되는 응력의 분포와 크기를 정량적으로 비교, 관찰하여 장치의 역계를 이해하고 치료효과를 예측해 보고자 실험을 고안하였다. 발치나 교정치료의 경험이 없는 25세 남자 정상교합자를 대상으로 컴퓨터 단층촬영을 시행하여 얻은 방사선필림을 기초로 하여 유한요소 상악치아모델과, 치근막모델을 제작하였으며 여기에 beam element로 제작된 MEAW의 tip back bend를 5도, 10도, 15도로 달리하여 II급 고무줄 및 수직고무줄을 적용시켜 각 치근막에서의 압축, 인장력의 분포와 응력의 크기를 상호 비교한 결과 다음과 같은 결론을 얻었다. 1. MEAW에 tip back bend를 5도, 10도, 15도로 부여하고 II급 고무줄이나 수직고무줄을 적용시켰을때 전체 치근막에 나타나는 압축, 인장의 분포는 II급 고무줄, 수직고무줄 각각의 경우에서 tip back bend의 정도에 상관없이 유사한 양상을 나타냈다. 2. 각 치아의 치근막에 나타난 압축, 인장력의 크기는 5도,10도의 경우보다 15도 tip back bend 부여시 더 크게 나타났다. 3. 제1, 2대구치 원심측 치근막에서는 압축력이 나타났으며 특히 제2대구치 치근막의 압축부위가 제1대구치 치근막 압축부위보다 더 넓고 압축력의 크기도 더 크게 나타났다. 또한 tip back bend 정도가 클수록 압축력도 크고 압축부위도 더 넓게 나타났다. 4. 대구치 치근막의 압축력은 소구치 치근막의 압축력보다 상대적으로 더 작게 나타났다. 5. II급고무줄과 수직고무줄 적용 시를 비교해 보면 수직고무줄 장착시 전치부 치근막에 인장력이 더 넓게, 그 크기는 더 작게 나타났다. 반면에 구치부 치근막에 나타나는 응력의 분포와 크기는 별 차이를 보이지 않았다. 5. 전치부 치근막 인장부위에서 인장력은 견치에서 제일 컸다.

• 인터넷정보학회논문지
• /
• 제14권5호
• /
• pp.1-10
• /
• 2013

RFID(Radio Frequency Identification)시스템은 하나의 RFDI리더, 다수의 RFID태그 장치들로 이루어진 비접촉방식의 근거리 무선 인식 기술이다. RFID태그는 자체적인 연산 수행이 가능한 능동형 태그와 이에 비해 성능은 떨어지지만 저렴한 가격으로 물류 유통에 적합한 수동형 태그로 나눌 수 있다. 데이터 처리 장치는 리더와 연결되어 리더가 전송받은 정보를 처리한다. RFID 시스템은 무선주파수를 이용해 다수의 태그를 빠른 시간에 인식할 수 있다. RFID시스템은 유통, 물류, 운송, 물품관리, 출입 통제, 금융 등 다양한 분야에서 응용되고 있다. 하지만 RFID시스템을 더욱 확산시키기 위해서는 가격, 크기, 전력소모, 보안 등 해결할 문제가 많다. 그 문제들 중에서 본 논문에서는 다수의 수동형 태그를 인식할 때 발생하는 충돌 문제를 해결하기 위한 알고리즘을 제안한다. RFID 시스템에서 다수의 태그를 인식하기 위한 충돌 방지 기법에는 확률적인 방식과 결정적인 방식 그리고 이를 혼합한 하이브리드 방식이 있다. 본 논문에서는 우선 기존에 있던 확률적 방식의 충돌방지기법인 알로하 기반 프로토콜과 결정적 방식의 충돌방지기법인 트리 기반 프로토콜에 대해 소개한다. 알로하 기반 프로토콜은 시간을 슬롯 단위로 나누고 태그들이 각자 임의로 슬롯을 선택하여 자신의 ID를 전송하는 방식이다. 하지만 알로하 기반 프로토콜은 태그가 슬롯을 선택하는 것이 확률적이기 때문에 모든 태그를 인식하는 것을 보장하지 못한다. 반면, 트리 기반의 프로토콜은 리더의 전송 범위 내에 있는 모든 태그를 인식하는 것을 보장한다. 트리 기반의 프로토콜은 리더가 태그에게 질의 하면 태그가 리더에게 응답하는 방식으로 태그를 인식한다. 리더가 질의 할 때, 두 개 이상의 태그가 응답 한다면 충돌이라고 한다. 충돌이 발생하면 리더는 새로운 질의를 만들어 태그에게 전송한다. 즉, 충돌이 자주 발생하면 새로운 질의를 자주 생성해야하기 때문에 속도가 저하된다. 그렇기 때문에 다수의 태그를 빠르게 인식하기 위해서는 충돌을 줄일 수 있는 효율적인 알고리즘이 필요하다. 모든 RFID태그는 96비트의 EPC(Electronic Product Code)의 태그ID를 가진다. 이렇게 제작된 다수의 태그들은 회사 또는 제조업체에 따라 동일한 프리픽스를 가진 유사한 태그ID를 가지게 된다. 이 경우 쿼리 트리 프로토콜을 이용하여 다수의 태그를 인식 하는 경우 충돌이 자주 일어나게 된다. 그 결과 질의-응답 수는 증가하고 유휴 노드가 발생하여 식별 효율 및 속도에 큰 영향을 미치게 된다. 이 문제를 해결하기 위해 충돌 트리 프로토콜과 M-ary 쿼리 트리 프로토콜이 제안되었다. 하지만 충돌 트리 프로토콜은 쿼리 트리 프로토콜과 마찬가지로 한번에 1비트씩 밖에 인식을 못한다는 단점이 있다. 그리고 유사한 태그ID들이 다수 존재할 경우, M-ary 쿼리 트리 프로토콜을 이용해 인식 하면, 불필요한 질의-응답이 증가한다. 본 논문에서는 이러한 문제를 해결하고자 M-ary 쿼리 트리 프로토콜의 매핑 함수를 이용한 m-비트 인식, 맨체스터 코딩을 이용한 태그 ID의 충돌정보, M-ary 쿼리 트리의 깊이를 하나 감소시킬 수 있는 예측 기법을 이용하여 성능을 향상시킨 적응형 M-ary 쿼리트리 프로토콜을 제안한다. 본 논문에서는 기존의 트리기반의 프로토콜과 제안하는 기법을 동일한 조건으로 실험하여 비교 분석 하였다. 그 결과 제안하는 기법은 식별시간, 식별효율 등에서 다른 기법들보다 성능이 우수하다.

• 한국의학물리학회지:의학물리
• /
• 제22권4호
• /
• pp.190-197
• /
• 2011

고선량률 근접치료에 사용되는 상업용 선원과 치료계획 시스템들은 AAPM TG 43에서 권고하는 점 및 선 선원에 의해 선량분포를 계산한다. 하지만, 근접치료용 선원에 대한 인체 내의 정확한 선량계산을 위해서 3차원 부피의 선원을 고려하는 MC 기반의 선량계산 방법이 필요하다. 본 연구에서는 microSelectron HDR Ir-192 선원을 작은 부분으로 분할하여 계산하는 미소선원 적분법을 이용하여 기하학적 인수를 계산하였다. 또한, 범용 방사선 수송코드인 MCNPX를 사용하여 30 cm 직경의 구형 물 팬텀 내에서 선원의 선량률을 계산하여 비등방성함수와 반경선량함수를 구하였다. 그 결과를 MC 기반 광자 수송코드인 MCPT를 사용하여 계산한 Williamson의 결과와 비교 및 분석하였다. 미소선원 적분법과 선 선원 근사법에 따른 기하학적 인수는 $r{\geq}0.5cm$에서는 0.2% 이내에서 일치하였고 r=0.1 cm일 때 1.33%의 차이를 보였다. 본 연구에서 계산된 비등방성함수와 반경선량함수가 Williamson의 계산된 결과의 차이는 비등방성함수의 경우 r=0.25 cm에 서 2.33%의 가장 큰 R-RMSE를 보였고 $r{\geq}0.5cm$에서는 1% 미만의 R-RMSE를 보였다. 반경선량함수의 경우는 r=0.1~14.0 cm에서 0.46%의 R-RMSE를 보였다. 미소선원 적분법과 선 선원 근사법으로 계산한 기하학적 인수는 $r{\geq}0.1cm$에서 잘 일치하지만 3차원의 Ir-192 선원을 적용하여 계산한 미소선원 적분법이 실제 기하학적 인수를 잘 반영할 것으로 생각된다. r=0.25 cm에서 비등방성함수를 제외하고는 MCPT와 MCNPX의 몬테칼로 코드를 이용하여 얻어진 비등방성함수와 반경선량함수는 각각의 몬테칼로 코드에 대한 불확실성 이내에서 잘 일치함을 확인하였다. 따라서 MCNPX 전산모사 결과를 통해 TG-43의 선량 계산식에 사용된 인자를 Williamson 등의 결과와 비교 및 검증함으로써, 추후 다른 종류의 선원에 대해서도 Monte Carlo 기반의 연구가 가능할 것으로 기대된다.

• 지적과 국토정보
• /
• 제47권2호
• /
• pp.107-120
• /
• 2017

산불은 한번 발생하면 기상, 지형 등 여러 악조건으로 인해 효과적인 진화가 어려워 넓은 면적으로 확대되는 경우가 많다. 따라서 산불의 예방이 중요하기 때문에 세계 각국에 다양한 산불위험지수와 예측시스템이 존재한다. 그러나 이러한 산불위험지수 및 지면건조지수가 우리나라의 산불발생에 적용가능한지에 대한 객관적인 평가는 이루어진 바 없다. 이에 본 연구에서는 1.5km 격자의 LDAPS(Local Analysis and Prediction System) 기상자료 및 1km 격자의 MODIS(Moderate-resolution Imaging Spectroradiometer) 위성자료를 활용하여 각종 산불위험지수와 지면건조지수의 우리나라 산불발생에 대한 민감도분석을 수행하고자 한다. 기상기반 산불위험지수로는 호주의 FFDI(forest fire danger index), 캐나다의 FFMC(fine fuel moisture code), 미국의 HI(Haines index), 그리고 학술연구에서 제시된 MNI(modified Nesterov index)를 산출하였고, 위성기반 지면건조지수인 NDDI(normalized difference drought index)와 TVDI(temperature vegetation dryness index)를 산출하여 우리나라 산불발생에 대한 적용가능성 실험을 수행하였다. 2013년 1월부터 2017년 5월까지 발생한 피해면적 1ha가 넘는 산불 120건과 6종류의 지수를 비교한 결과, FFDI는 피해면적 10ha가 넘는 모든 산불에 대하여 극도로 높은 CDF(cumulative density function) 값을 나타냈으며, FFDI와 FFMC는 피해면적 3ha가 넘는 산불에 대하여 평균 CDF 값이 0.95가 넘게 나타나는 등 매우 우수한 성능을 보였다. 반면, MNI는 이슬점온도와 기온의 차이가 크지 않은 우리나라의 계절적 특성 때문에 2월의 산불예측을 거의 하지 못하였고, TVDI는 전체적으로 산불발생에 대한 민감도가 낮은 것으로 나타났다. NDDI는 피해면적의 크기에 상관없이 평균 CDF 값이 안정적으로 높게 산출되어 위성기반 지면건조지수로서 보조적인 활용가능성이 있을 것으로 보인다. 이러한 산불위험지수와 지면건조지수를 취사선택 및 융합하여 활용한다면, 우리나라 산불예측에 일조할 수 있을 것으로 사료된다.

• 대한원격탐사학회지
• /
• 제13권2호
• /
• pp.99-120
• /
• 1997

위성 탐사 지구 관측 자료, 특히 가시광선 영역의 자료는 태양광이 지구 대기계와 복잡 한 상호 작용을 거친후 위성센서에 의해서 감지되어 수집된 것이다. 따라서 가시광선 영역의 지 구 관측 위성 자료를 정량적으로 분석하기 위해서는 대기의 산란과 흡수에 의한 대기효과에 대한 정량적 보정이 필요하다. 본 연구에서는 대기효과에 대한 정량적 보정을 위해서 복사 전달 모델 의 활용 가능성을 조사하였다. 이를 위해 위성 원격탐사 응용 분야에 이용되고 있는 복사 전달 모델 중의 하나인 LOWTRAN7의 복사 전달 모의 결과 중 단파 자료를 복사 관측 자료인 CAGEX 자료와 비교하였다. CAGEX 관측 자료는 대기복사 모델의 검정을 위하여 NASA Langley Research Center에서 수집한 자료로써, 1)대기 sounding, aerosols, 구름의 특성 자료, 2) 지구복사계, 직달일사계, 차광전천일사계에 의해 측정된 ARM(Aerosol Radiation Measurement) 자료, 3)Fu-Luio 모델에 의해 모의된 장파, 단파 flux 등의 자료로 구성되어 있다. 대기복사 전달 모의를 위하여 에어로졸의 광학적 특성은 CAGEX의 column optical depth, Spinhime의 산란 계 수 수직 분포와 D' Almeida의 에어로졸 복사 특성 값으로부터 도출되었다. LOWTRAN7의 복사 모의는 완전히 맑은 날에 해당하는 31개의 경우에 대하여 수행되었으며, 이 모의 결과 중 단파영 역에서의 지표면에서의 상향복사와 지표면으로의 하향직달복사 및 하향 확산 복사 그리고 대기 상단에서의 상향 복사를 CAGEX 관측 자료와 각각 비교 하였다. 비교 결과 CAGEX 자료에 대한 LOWTRAN7 결과의 표준 오차는 지표면에서의 하향 확산 복사(6.9%)를 제외한 모든 복사 항목 들이 5%이내였다. 이 결과로 보아 하향 확산 복사항의 오차가 가장 크며 이 오차가 관련된 나머 지 항목의 오차를 일으키는 역할을 하는 것으로 추측할 수 있다. 결론적으로 지표면에서의 하향 확산 복사 항목은 에어로졸에 의해 생기는 항목이므로 향후 복사 모델을 고려할 때 에어로졸의 산란에 의한 부분이 더 고려 된다면 복사 전달 모델을 이용하여 정량적으로 대기 효과를 보정하 는데 오차를 줄일 수 있을 것으로 기대된다.

• 천문학논총
• /
• 제12권1호
• /
• pp.111-143
• /
• 1997

We have developed a spherical FCT code in order to simulate the interaction of supernova remnants with stellar wind bubbles. We assume that the density profile of the supernova ejecta follows the Chevalier mode1(1982) where the outer portion has a power-law density distribution($\rho{\propto}\gamma^{-n}$) and the SN ejecta has a kinetic energy of $10^{51}$ ergs. The structure of wind bubble has been calculated with the stellar mass loss rate $\dot{M}=5\times10^{-6}M_{\odot}/yr$ and the wind velocity $\upsilon=2\times10^3$ km/s We have simulated seven models with different initial conditions In the first two models we computed the evolution of SNRs with n=7 and n=14 in the uniform medium The numerical results agree with the Chevalier's similarity solution at early times. When all of the power-law portion of the ejecta is swept up by the reverse shock, the evolution slowly converges to the Sedov-Taylor stage. There is not much difference between the two cases with different n's The other five models simulate SNRs produced inside wind bubbles. In model III, we consider the SN ejecta of 1.4 $M_{\odot}$ and the radius of bubble ~2.76 pc so that ratio of the mass $\alpha(=M_{W.S}/M_{ej}$ is 2. We follow the complex hydrodynamic flows produced by the interaction of SN shocks with stellar shocks and with the contact discontinuities, In the model III, the time scale for the SN shock to cross the wind shell $\tau_{cross}$ is similar to the time scale for the reverse shock to sweep the power-law density profile $\tau_{bend}$. Hence the SN shock crosses the wind shell. At late times SN shock produces another shell in the ambient medium so that we have a SNR with double shell structure. From the numerical results of the remaining models, we have found that when $\tau_{cross}/\tau_{bend}\leq2$, or equivalently when $\alpha\leq50$, the SNRs produced inside wind bubbles have double shell structure. Otherwise, either the SN shock does not cross the wind shell or even if it crosses at one time, the reverse shock reflected at the center accelerates the wind shell to merge into the SN shock Our results confirm the conclusion of Tenorio-Tagle et a1(1990).