• 한국전산유체공학회:학술대회논문집
• /
• 1998.05a
• /
• pp.138-144
• /
• 1998

The BGK equation, which is the kinetic model equation of Boltzmann equation, is solved using FDDO(finite difference with the discrete-ordinate method) to compute the rarefied flow of monatomic gas. Using reduced velocity distribution and discrete ordinate method, the scalar equation is transformed into a system of hyperbolic equations. High resolution ENO(Essentially Non-Oscillatory) scheme based on Harten-Yee's MFA(Modified Flux Approach) method with Strang-type explicit time integration is applied to solve the system equations. The calculated results are well compared with the experimental density field of NACA0012 airfoil, validating the developed computer code. Next. the computed results of circular cylinder flow for various Knudsen numbers are compared with the DSMC(Direct Simulation Monte Carlo) results by Vogenitz et al. The present scheme is found to be useful and efficient far the analysis of two-dimensional rarefied gas flows, especially in the transitional flow regime, when compared with the DSMC method.

• 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.