• Laser Solutions
• /
• v.6 no.3
• /
• pp.37-45
• /
• 2003

Pulsed laser ablation is important in a variety of engineering applications involving precise removal of materials in laser micromachining and laser treatment of bio-materials. Particularly, detailed numerical simulation of complex laser ablation phenomena in air, taking the interaction between ablation plume and air into account, is required for many practical applications. In this paper, high-power pulsed laser ablation under atmospheric pressure is studied with emphasis on the vaporization model, especially recondensation ratio over the Knudsen layer. Furthermore, parametric studies are carried out to analyze the effect of laser fluence and background pressure on surface ablation and the dynamics of ablation plume. In the numerical calculation, the temperature, pressure, density, and vaporization flux on a solid substrate are obtained by a heat-transfer computation code based on the enthalpy method. The plume dynamics is calculated considering the effect of mass diffusion into the ambient air and plasma shielding. To verify the computation results, experiments for measuring the propagation of a laser induced shock wave are conducted as well.

• Journal of the Korean Society of Radiology
• /
• v.11 no.6
• /
• pp.423-428
• /
• 2017

In this study, the calculation of the effective spatial dose distribution of the diagnostic imaging laboratory of K university was performed by the Monte Carlo simulation. The radiation generator has a maximum tube voltage of 150 kVp and a maximum current of 700 mA. Using the results, we compared the spatial effective dose distributions of diagnostic imaging laboratory when the shielding door was closed and opened. In conclusion, it was found that the effective dose in the operating room of the diagnostic imaging laboratory does not exceed the annual dose limit (6 mSv/y) of the student (occasional visitor) even when the door is opened. However, since the effective dose when the door is open is about 16 times higher in front of the lead glass window and about 3,000 times higher in front of the doorway than the case when the door is closed, closing the shielding door at the time of the practical exercising reduces unnecessary radiation exposure by great extent.

• Journal of Radiation Protection and Research
• /
• v.39 no.4
• /
• pp.176-181
• /
• 2014

The object of this paper is to evaluate the fission product inventories and radiological doses in a non-LOCA event, based on the U.S. NRC's regulatory methodologies recommended by the TID-14844 and the RG 1.195. For choosing a non-LOCA event, one fuel assembly was assumed to be melted by a channel blockage accident. The Hanul nuclear power reactor unit 6 and the CE $16{\times}16$ fuel assembly were selected as the computational models. The burnup cross section library for depletion calculations was produced using the TRITON module in the SCALE6.1 computer code system. Based on the recently licensed values for fuel enrichment and burnup, the source term calculation was performed using the ORIGEN-ARP module. The fission product inventories released into the environment were obtained with the assumptions of the TID-14844 and the RG 1.195. With two kinds of source terms, the radiological doses of public in normal environment reflecting realistic circumstances were evaluated by applying the average condition of meteorology, inhalation rate, and shielding factor. The statistical analysis was first carried out using consecutive three year-meteorological data measured at the Hanul site. The annual-averaged atmospheric dispersion factors were evaluated at the shortest representative distance of 1,000 m, where the residents are actually able to live from the reactor core, according to the methodology recommended by the RG 1.111. The Korean characteristic-inhalation rate and shielding factor of a building were considered for a series of dose calculations.

• Nuclear Engineering and Technology
• /
• v.51 no.8
• /
• pp.2018-2025
• /
• 2019

This study systematically compared two hybrid deterministic/Monte Carlo transport codes, ADVANTG/MCNP and MAVRIC, in solving a difficult shielding problem for a real-world spent fuel storage cask. Both hybrid codes were developed based on the consistent adjoint driven importance sampling (CADIS) methodology but with different implementations. The dose rate distributions on the cask surface were of primary interest and their predicted results were compared with each other and with a straightforward MCNP calculation as a baseline case. Forward-Weighted CADIS was applied for optimization toward uniform statistical uncertainties for all tallies on the cask surface. Both ADVANTG/MCNP and MAVRIC achieved substantial improvements in overall computational efficiencies, especially for gamma-ray transport. Compared with the continuous-energy ADVANTG/MCNP calculations, the coarse-group MAVRIC calculations underestimated the neutron dose rates on the cask's side surface by an approximate factor of two and slightly overestimated the dose rates on the cask's top and side surfaces for fuel gamma and hardware gamma sources because of the impact of multigroup approximation. The fine-group MAVRIC calculations improved to a certain extent and the addition of continuous-energy treatment to the Monte Carlo code in the latest MAVRIC sequence greatly reduced these discrepancies. For the two continuous-energy calculations of ADVANTG/MCNP and MAVRIC, a remaining difference of approximately 30% between the neutron dose rates on the cask's side surface resulted from inconsistent use of thermal scattering treatment of hydrogen in concrete.

• Journal of Radiation Protection and Research
• /
• v.19 no.1
• /
• pp.13-22
• /
• 1994

The MCNP 4.2 code was used to calculate the thermal neutron flux distributions for $(n,\;{\gamma})$reaction in mainshell, annular plate, and subshell of the calandria of a CANDU 6 plant during operation. The thermal neutron flux distributions in calandria mainshell, annular plate, and subshell were in the range of $10^{11}{\sim}10^{13}\;neutrons/cm^2-sec$ which is somewhat higher than the previous estimates calculated by DOT 4.2 code. As an application to shielding analysis, photon dose rates outside the side and bottom shields were calculated. The resulting dose rates at the reactor accessible areas were below design target, $6 {\mu}Sv/h$. The methodology used in this study to evaluate the thermal neutron flux distribution for $(n,\;{\gamma})reaction$ can be applied to radiation shielding analysis of CANDU 6 type plants.

• Nuclear Engineering and Technology
• /
• v.47 no.7
• /
• pp.791-803
• /
• 2015

The effective cross sections (XSs) in the direct whole core calculation code nTRACER are evaluated by the equivalence theory-based resonance-integral-table method using the WIMS-based library as an alternative to the subgroup method. The background XSs, as well as the Dancoff correction factors, were evaluated by the enhanced neutron-current method. A method, with pointwise microscopic XSs on a union-lethargy grid, was used for the generation of resonance-interference factors (RIFs) for mixed resonant absorbers. This method was modified by the intermediate-resonance approximation by replacing the potential XSs for the non-absorbing moderator nuclides with the background XSs and neglecting the resonance-elastic scattering. The resonance-escape probability was implemented to incorporate the energy self-shielding effect in the spectrum. The XSs were improved using the proposed method as compared to the narrow resonance infinite massbased method. The RIFs were improved by 1% in $^{235}U$, 7% in $^{239}Pu$, and >2% in $^{240}Pu$. To account for thermal feedback, a new feature was incorporated with the interpolation of pre-generated RIFs at the multigroup level and the results compared with the conventional resonance-interference model. This method provided adequate results in terms of XSs and k-eff. The results were verified first by the comparison of RIFs with the exact RIFs, and then comparing the XSs with the McCARD calculations for the homogeneous configurations, with burned fuel containing a mixture of resonant nuclides at different burnups and temperatures. The RIFs and XSs for the mixture showed good agreement, which verified the accuracy of the RIF evaluation using the proposed method. The method was then verified by comparing the XSs for the virtual environment for reactor applicationbenchmark pin-cell problem, as well as the heterogeneous pin cell containing burned fuel with McCARD. The method works well for homogeneous, as well as heterogeneous configurations.

• Journal of Environmental Science International
• /
• v.23 no.2
• /
• pp.207-218
• /
• 2014

A ventilation model was developed for predicting the air change per hour(ACH) in buildings and the airflow rates between zones of a multi-room building. In this model, the important parameters used in the calculation of airflow are wind velocity, wind direction, terrain effect, shielding effect by surrounding buildings, the effect of the window type and insect screening, etc. Also, the resulting set of mass balance equations required for the process of calculation of airflow rates are solved using a Conte-De Boor method. When this model was applied to the building which had been tested by Chandra et al.(1983), the comparison of predicted results by this study with measured results by Chandra et al. indicated that their variations were within -10%~+12%. Also, this model was applied to a building with five zones. As a result, when the wind velocity and direction did not change, terrain characteristics influenced the largest and window types influenced the least on building ventilation among terrain characteristics, local shieldings, and window types. Except for easterly and westerly winds, the ACH increased depending on wind velocity. The wind direction had influence on the airflow rates and directions through openings in building. Thus, this model can be available for predicting the airflow rates within buildings, and the results of this study can be useful for the quantification of airflow that is essential to the research of indoor air quality(temperature, humidity, or contaminant concentration) as well as to the design of building with high energy efficiency.

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

• Progress in Medical Physics
• /
• v.13 no.1
• /
• pp.21-26
• /
• 2002

A practical calculation algorithm which calculates the relative output factor(ROF) for irregular shaped electron field has been developed and evaluated the accuracy of the algorithm. The algorithm adapted two-source model, which assumes that the electron dose can be express as sum of the primary source component and the scattered component from the shielding block. Original two-source model has been modified in order to make the algorithm simpler and to reduce the number of parameters needed in the calculation, while the calculation error remains within clinical tolerance range. The primary source is assumed to have Gaussian distribution, while the scattered component follows the inverse square law. Depth and angular dependency of the primary and the scattered are ignored ROF can be calculated with three parameters such as, the effective source distance, the variance of primary source, and the scattering power of the block. The coefficients are obtained from the square shaped-block measurements and the algorithm is confirmed from the rectangular or irregular shaped-fields used in the clinic. The results showed less than 1.0 ％ difference between the calculation and measurements for most cases. None of cases which have bigger than 2.1 ％ have been found. By improving the algorithm for the aperture region which shows the largest error, the algorithm could be practically used in the clinic, since one can acquire the 1011 parameter's with minimum measurements(5∼6 measurements per cones) and generates accurate results within the clinically acceptable range.

• Journal of Radiation Protection and Research
• /
• v.6 no.1
• /
• pp.8-24
• /
• 1981

This paper presents flux-to-dose-rate conversion factors for neutrons and gamma rays based on the American National Standard Institute(ANSI) N666. These data are used to calculated the dose rate distribution of neutron and gamma ray in radiation fields. Neutron flux-to-dose-rate conversion factors for energies from $2.5{\times}10^{-8}$ to 20 MeV are presented; the corresponding energy range for gamma rays is 0.01 to 15 MeV. Flux-to-dose-rate conversion factors were calculated, under the assumption that radiation energy distribution has nonlinearity in the phantom, have different meaning from those values obtained by monoetiergetic radiation. Especially, these values were determined with the cross section library. The flux-to-dose-rate conversion factors obtained in this work were in a good agreement to the values presented by ANSI. Those data will be a useful for the radiation shielding analysis and the radiation dosimetry in the case of continuous energy distributions.