• The Bulletin of The Korean Astronomical Society
• /
• v.36 no.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.07a
• /
• 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.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2004.03a
• /
• 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.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2008.03a
• /
• 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.

• Journal of the Korean Society of Radiology
• /
• v.13 no.5
• /
• pp.791-799
• /
• 2019

We evaluated the physical properties that occur to dose enhancement and changes from secondary particle production resulting from the interaction between enhancement material. Geant4 was used to perform a Monte Carlo simulation, and the medical internal radiation dose (MIRD) head phantom were employed. X-rays of 4, 6, 10, 15, 18, and 25 MV were used. Aurum (Au) and gadolinium (Gd) were applied within the tumor volume at 10, 20, and 30 mg/g, and an experiment using soft tissue exclusively was concomitantly performed for comparison. Also, particle fluence and initial kinetic energy of secondary particle of interaction were measured to calculate equivalent doses using the radiation weight factor. The properties of physical interaction by the radiation enhancement material showed the great increased in photoelectric effect as compared to the compton scattering and pair production, occurred with the highest, in aurum and gadolinium it is shown in common. The photonuclear effect frequency increased as the energy increased, thereby increasing secondary particle production, including alpha particles, protons, and neutrons. During dose enhancement using aurum, a maximum 424.25-fold increase in the equivalent dose due to neutrons was observed. This study was Monte Carlo simulation corresponds to the physical process of energy transmission in dose enhancement. Its results may be used as a basis for future in vivo and in vitro experiments aiming to improve effects of dose enhancement.

• The Bulletin of The Korean Astronomical Society
• /
• v.42 no.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.

• Journal of the Korean Ceramic Society
• /
• v.53 no.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.

• Journal of Korean Tunnelling and Underground Space Association
• /
• v.20 no.1
• /
• pp.11-22
• /
• 2018

The safety of tunnels is quantified by quantitative risk assessment when planning the disaster prevention facilities of railway tunnels, and it is decided whether they are appropriate. The purpose of this study is to estimate the probability of the train stopping in the tunnels at train fire, which has a significant effect on the results of quantitative risk assessment for tunnel fires. For this purpose, a model was developed to calculate the coasting distance of the train considering the coefficient of train running resistance. The probability of stopping in case of train fire in the tunnel is predicted by the Monte Carlo simulation method with the coasting distance and the emergency braking distance as parameters of the tunnel lengths and slopes, train initial driving speeds. The kinetic equations for predicting the coasting distance were analyzed by reflecting the coefficient train running resistance of KTX II. In the case of KTX II trains, the coasting distance is reduced as the slope increases in a tunnel with an upward slope, but it is possible to continue driving without stopping in a slope downward. The probability of the train stopping in the case of train fire in tunnel decreases as the train speed increases and the slope of the tunnel decreases. If human error is not taken into account, the probability that a high-speed train traveling at a speed of 250 km/h or above will stop in a tunnel due to a fire is 0% when the slope of the tunnel is 0.5% or less, and the probability of stopping increases rapidly as the tunnel slope increases and the tunnel length increases.

• Nuclear Engineering and Technology
• /
• v.56 no.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.

• Proceedings of the KOR-KST Conference
• /
• 1996.02a
• /
• 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.