• 천문학회보
• /
• 제36권2호
• /
• pp.114.2-114.2
• /
• 2011

We have developed a PDR code to reproduce the high rotational transitions of CO observed with Herschel-PACS. Part of these high-J CO line emission is produced by UV heated outflow walls around protostars. The local FUV radiation flux is calculated by using Monte Carlo method in (${\gamma}$, ${\alpha}$) grid taking anisotropic scattering into account. Kinetic temperature and Abundance of molecules were computed self-consistently. CO Line fluxes are calculated using RIG. We compare our PDR model with the results by Visser et al (2011) to show that the derived FUV radiation field strength can be affected by the grid resolution near the outflow wall and dust scattering.

• 한국정보디스플레이학회:학술대회논문집
• /
• 한국정보디스플레이학회 2005년도 International Meeting on Information Displayvol.I
• /
• pp.167-170
• /
• 2005

The mechanism of striation in the coplanar- and matrix-type plasma display panel (PDP) cells has been studied using the particle-in-cell Monte-Carlo Collision (PIC-MCC) model. The striation formation is related to the ionization energy of neutral atoms and the well-like deformation of space potential by space charge distribution. Negative wall charge accumulation by electrons on the MgO surface of the anode region is also one of the key factors for the formation of striation. The clearness of the striation phenomenon in PIC-MCC code in comparison with fluid code can be explained by using nonlocal electron kinetic effect.

• 한국추진공학회:학술대회논문집
• /
• 한국추진공학회 2004년도 제22회 춘계학술대회논문집
• /
• pp.245-250
• /
• 2004

The discharge current oscillation has been measured for various hollow anode widths and its axial positions using a 1㎾-class anode layer hall thruster. As a result, there were thresholds of magnetic flux density for stable discharge. The plasma structure inside the hollow anode was numerically analyzed using the fully kinetic 2D3V Particle-in-Cell (PIC) and Direct Simulation Monte Carlo (DSMC) methods. The results reproduced both stable and unstable operation modes. In the stable operation case, which corresponds to the case with low magnetic flux, the plasma penetrated into the hollow anode deeper than the case with higher magnetic flux density case. This suggests that comparably large substantial anode area should contribute to stable operation.

• 한국추진공학회:학술대회논문집
• /
• 한국추진공학회 2008년 영문 학술대회
• /
• pp.602-605
• /
• 2008

The anode sheath structure in the hollow anode of an anode-layer type Hall thruster was numerically computed using a fully kinetic 2D3V Particle-in-Cell and Direct Simulation Monte Carlo(PIC-DSMC) code. By treating both ions and electrons as particles, anode surface region, which is electrically non-neutral, was analyzed. In order to analyze in detail, the calculation code was parallelized using Message Passing Interface (MPI). The code successfully simulated the discharge current oscillation. In the low magnetic induction case, ion sheath appears in the anode surface because ionization is enough to maintain the plasma occurs in the anode hollow. As the magnetic induction increases, main ionization region move to outside of the anode. At the same time, anode sheath voltage decreases. In the high magnetic induction case, electron sheath appears on the anode surface periodically because the ionization occurs mainly in the discharge channel. This anode sheath condition shift can be explained using the simple sheath model.

• 한국방사선학회논문지
• /
• 제13권5호
• /
• pp.791-799
• /
• 2019

선량 증강 현상에서 발생하는 물리적 특성과 증강 물질과의 상호 작용으로부터 발생하는 이차입자 생성을 평가하였다. Geant 4, MIRD 두부 팬텀을 이용한 몬테카를로 전산 모사를 진행하였으며, 선형가속기에서 발생되는 4, 6, 10, 15, 18, 25 MV X선을 선원으로 적용하였다. 10, 20, 30 mg/g의 금(aurum), 가돌리늄(gadolinium) 증강 물질을 팬텀 내부 종양에 모사하였으며, 물리적 상호작용의 변화와 이차입자 발생에 따른 입자플루언스와 초기 에너지로부터 방사선가중인자를 고려하여 등가선량을 평가하였다. 방사선 선량 증강 물질에 의한 상호작용은 고 원자번호에서 기인하여 광전효과에 의한 에너지 흡수를 높이는 것으로 나타났으며, 10 MV 이상의 에너지에서는 광핵반응의 증가를 나타내었다. 이로 인해, 팬텀 내부에서 양성자, 중성자와 같은 이차입자 발생의 증가를 보였으며, 중성자에 의한 등가선량이 최대 424.2배 증가하는 것으로 나타났다. 본 연구는 선량 증강 현상에서의 에너지 전달, 흡수의 물리적 과정을 모사하여, 증강 현상에서 발생하는 물리적 특성을 분석하고자 하였다. 이러한 결과는 향후 in-vivo, in-vitro 선량 증강 실험을 위한 기초 자료로 활용될 수 있을 것으로 사료된다.

• 천문학회보
• /
• 제42권2호
• /
• pp.43.1-43.1
• /
• 2017

A three-dimensional (3D) $Ly{\alpha}$ radiative transfer code is developed to study the Wouthuysen-Field effect, which couples the 21 cm spin temperature of neutral hydrogen and the $Ly{\alpha}$ radiation field, and the escape fraction of $Ly{\alpha}$ from galaxies. The Monte Carlo code is capable of treating arbitrary 3D distributions of $Ly{\alpha}$ source, neutral hydrogen and dust densities, gas temperature, and velocity field. It is demonstrated that the resonance-line profile at the center approaches to the Boltzmann distribution with the gas temperature. A plane-parallel ISM model, which is appropriate for the neutral ISM of our Galaxy, is used to calculate the $Ly{\alpha}$ radiation field strength as a function of height above the galactic plane. We also use a two-phase, clumpy medium model which is composed of the cold and warm neutral media (WNM). It is found that the $Ly{\alpha}$ radiation field is strong enough to thermalize the 21 cm spin temperature in the WNM to the gas kinetic temperature. The escape fraction of $Ly{\alpha}$ is found to be a few percent, which is consistent with the $Ly{\alpha}$ observations of our Galaxy and external galaxies.

• 한국세라믹학회지
• /
• 제53권5호
• /
• pp.521-528
• /
• 2016

Kinetic Monte Carlo (KMC) and graph searches show that proton conduction limiting barriers and trajectories in $BaZr_{0.875}Y_{0.125}O_3$ are affected by the presence of other protons. At 1000 K, KMC limiting conduction barriers increase from 0.39 eV to 0.45 eV as the proton number is increased. The proton-proton radial distribution begins to rise at $2{\AA}$ and peaks at $4{\AA}$, which is half the distance expected, based on the proton concentration. Density functional theory (DFT) calculations find proton/proton distances of 2.60 and $2.16{\AA}$ in the lowest energy two-proton configurations. A simple average of the limiting barriers for 7-10 step periodic long range paths found via graph theory at 1100 K shows an increase in activation barrier from 0.32 eV to 0.37 eV when a proton is added. Both KMC and graph theory show that protons can affect each other's pathways and raise the overall conduction barriers.

• 한국터널지하공간학회 논문집
• /
• 제20권1호
• /
• pp.11-22
• /
• 2018

철도터널의 방재시설 계획 시 터널의 안전성을 정량적 위험도 평가에 의해서 정량화하여 방재시설의 적정성 여부를 판단하도록 하고 있다. 이에 본 연구에서는 터널화재에 대한 정량적 위험도 평가 시 결과에 크게 영향을 미치는 화재열차의 터널 내 정차확률을 예측하는 방법을 제시하였다. 이를 위해서 열차의 주행저항계수를 고려하여 타력운전거리를 계산하기 위한 모델을 개발하였으며, 타력운전특성과 비상제동거리를 고려하여 몬테카를로 시뮬레이션 기법에 의해 터널연장 및 경사도, 초기주행속도를 변수로 하여 화재열차가 터널에 정차할 확률을 예측하였다. 타력운전거리의 예측을 위한 운동방정식은 KTX II의 주행저항계수를 반영하여 분석하였다. KTX II 열차의 경우, 타력운전거리는 상향경사의 터널에서는 경사도가 증가할수록 감소하나 하향경사구간에서는 정지하지 않고 계속하여 주행이 가능하다. 화재열차의 터널 내 정차확률은 열차의 주행속도가 증가할수록, 경사도가 낮을수록 감소하며, 고속열차(주행속도 250 km/h 이상)는 인적오류를 고려하지 않는 경우, 경사도 0.5% 이하의 터널에서는 화재 시 열차가 터널에 정차할 확률은 0%이고, 경사도가 증가하고 터널연장이 증가하면 정차확률은 급격하게 증가하는 것으로 나타났다.

• Nuclear Engineering and Technology
• /
• 제56권4호
• /
• pp.1204-1212
• /
• 2024

To accurately calculate the heating distribution of the fast reactor, a neutron-photon library in MATXS format named Knight-B7.1-1968n × 94γ was processed based on the ENDF/B-VII.1 library for ultrafine groups. The neutron cross-section processing code MGGC2.0 was used to generate few-group neutron cross sections in ISOTXS format. Additionally, the self-developed photon cross-section processing code NGAMMA was utilized to generate photon libraries for neutron-photon coupled heating calculations, including photo-atom cross sections for the ISOTXS format, prompt photon production cross sections, and kinetic energy release in materials (KERMA) factors for neutrons and photons, and the self-shielding effect from the capture and fission cross sections of neutron to photon have been taken into account when the photon source generated by neutron is calculated. The interface code GSORCAL was developed to generate the photon source distribution and interface with the DIF3D code to calculate the neutron-photon coupling heating distribution of the fast reactor core. The neutron-photon coupled heating calculation route was verified using the ZPPR-9 benchmark and the RBEC-M benchmark, and the results of the coupled heating calculations were analyzed in comparison with those obtained from the Monte Carlo code MCNP. The calculations show that the library was accurately processed, and the results of the fast reactor neutron-photon coupled heating calculations agree well with those obtained from MCNP.

• 대한교통학회:학술대회논문집
• /
• 대한교통학회 1996년도 제29회 학술발표회
• /
• pp.59-86
• /
• 1996

Current methods for evaluating unsignalized intersections, and estimating level-of-service (LOS) is determined from efficiency-based criteria such as little or no delay to very long delays. At present, similar procedures to evaluate intersections using safety-based criteria do not exist. The improvement of sight distances at intersections is the most effective way of improving intersection safety. However, a set of procedures is necessary to account for the limitations in current methodology. Such an approach would build upon such methods, but also account for: deficiencies in the current deterministic solution for the determination of intersection sight distances; opportunity for an accident and severity of an accident; and cost-effectiveness of attaining various levels of sight distances. In this research, a model that estimates the degree of safety at two-way stop-controlled intersections is described. Only crossing maneuvers are considered in this study because accidents caused by the crossing maneuvers are the dominate type among intersection accidents. Monte Carlo methods are used to estimate the hazard at an intersection as a function of roadway features and traffic conditions. Driver`s minimum gap acceptance in the crossing vehicles and headway distribution on the major road are used in the crossing vehicles and headway distribution on the major road are used in the model to simulate the real intersectional maneuvers. Other random variables addressed in the model are: traffic speeds; preception-reaction times of both drivers in the crossing vehicles and drivers in oncoming vehicles on the major road; and vehicles on the major roads. The developed model produces the total number of conflicts per year per vehicle and total potential kinetic energy per year per vehicle dissipated during conflicts as measurements of safety at intersections. Based on the results from the developed simulation model, desirable sight distances for various speeds were determined as 350 feet, 450 feet and 550 feet for 40 mph, 50 mph and 60 mph prevailing speed on the major road, respectively. These values are seven to eight percent less than those values recommended by AASHTO. A safety based level-of-service (LOS) is also developed using the results of the simulation model. When the total number of conflicts per vehicle is less than 0.05 at an intersection, the LOS of the intersection is `A' and when the total number of conflicts per vehicle is larger than 0.25 at an intersection, the LOS is `F'. Similarly, when the total hazard per vehicle is less than 350, 000 1b-ft2/sec2, the LOS is `F'. Once evaluation of the current safety at the intersection is complete, a sensitivity analysis can be done by changing one or more input parameters. This will estimate the benefit in terms of time and budget of hazard reduction based upon improving geometric and traffic characteristics at the intersection. This method will also enable traffic engineers in local governments to generate a priority list of intersection improvement projects.