Montmorillonite plays a key role in engineered barrier systems in the high-level radioactive waste repository because of its large sorption capacity and high swelling pressure. However, the sorption capacity of montmorillonite can be largely varied dependent on the surrounding environments. This study conducted the batch simulation for U(VI) sorption on Na-montmorillonite by utilizing the cation exchange and surface complexation coupled (2SP-NE-SC/CE) model and evaluated the effects of physicochemical properties (i.e., pH, temperature, competing cations, U(VI) concentration, and carbonate species) on U(VI) sorption. The simulation demonstrated that the U(VI) sorption was affected by physicochemical properties: the pH and temperature relate to aqueous U(VI) speciation, the competing cations relate to the cation exchange process and selectivity, the U(VI) concentration relates to saturation at sorption sites. For example, the Kd (L kg−1) of Na-montmorillonite represented the largest values of 2.7×105 L kg−1 at neutral pH condition and had significantly decreased at acidic pH<3, showing non-linear and diverse U(VI) sorption at the ranged pH from 2 to 11. Additionally, the U(VI) sorption on montmorillonite significantly decreased in presence of carbonate species. The U(VI) sorption for long-term in actual porewater chemistry and temperature of high-level radioactive waste repository represented that the sorption capacity of Na-montmorillonite was affected by various external properties such as concentration of competing cation, temperature, pH, and carbonate species. These results indicate that geochemical sorption capacity of bentonite should be evaluated by considering both geological and aquifer environments in the high-level radioactive waste repository.

To minimize the short-term thermal load on the repository facility, heat generating nuclides such as Cs-137 and Sr-90 should be separated from the spent nuclear fuel for efficiency of repository facility. In particular, Sr-90 must be separated because it generates high heat during the decay process. Recently, Korea Atomic Energy Research Institute (KEARI) is developing a waste burden minimization technology to reduce the environmental burden caused by the disposal of spent nuclear fuel and maximize the utilization of the disposal facility. The technology includes a nuclide management process that can maximize disposal efficiency by selectively separating and collecting major nuclides such as Cs, Sr, I, TRU/RE, and Tc/Se from spent nuclear fuel. Among the major nuclides, Sr nuclides dissolve in chloride phase during the chlorination process of spent nuclear fuel and recovered in the form of carbonate or oxide via reactive distillation. In this process, Ba nuclides are also recovered along with Sr nuclides due to their chemical similarity. In this study, we prepared group II nuclide ceramic waste form, Ba(x)Sr(1-x)TiO3 (x=0, 0.25, 0.5, 0.75, 1), using the solid-state reaction method by considering the various ratio of Sr/Ba nuclides generated from nuclide management process. The established waste form fabrication process was able to produce a stable waste form regardless of the ratio of Sr/Ba nuclides. To evaluate the stability of group II waste form, physicochemical properties such as leaching and thermal properties were evaluated. Also, the radiological properties of the Ba(x)Sr(1-x)TiO3 waste forms with various Sr/Ba ratios were evaluated, and the estimation of centerline temperature was carried out using the experimental thermal property data. These results provided fundamental data for long-term storage and management of group II nuclides waste form.

One of the promising candidates for heat transfer fluid is molten chloride salts. They have been studied in various fields such as the electrolyte of pyroprocessing, the molten salt reactor coolant, and the energy storage system media. Main considerations for utilizing molten chloride salts are the compatibility of salts with structural materials. The corrosion behavior of structural materials in molten chloride salts must be understood to identify suitable materials against the corrosive environment. In this study, the corrosion behavior of a candidate structural material, Hastelloy N, in molten LiCl- KCl salt at 500°C were investigated by the electrochemical impedance spectroscopy (EIS) method. The sheet type of Hastelloy N was utilized as the working electrode in LiCl-KCl to measure the EIS data for 100 hours with 5 hours of time intervals. The EIS data were measured in the frequency range from 104 Hz to 10-2 Hz with the AC signal (amplitude = 20 mV) at open circuit potential. The capacitance semicircle observed in Nyquist plots for all periods indicates that charge-transfer controlled reactions occur. As the immersion time progresses, the radius of the semicircle in Nyquist plots and the impedance and phase angle in Bode plots decrease. These behaviors suggest a decreasing reaction resistance and the corrosion reactions are accelerated with the immersion time. The EIS data were fitted using the equivalent circuit to achieve quantitative results. Two capacitor-resistor components were considered due to the overlapped shape of two valleys in phase angle. The depressed shape of the semicircle in Nyquist plots led to the use of the constant phase element(Q) instead of the capacitor. Therefore, R(Q(R(QR))) circuit was selected to fit the EIS data. Fitting results show that the charge transfer resistance decreases dramatically within 1 day and then converges. The film resistance shows no clear trends, but the increase of the film admittance value indicates the decreased film thickness. Consequently, the film appears to exist like the oxide layer but it does not act as a protective layer. The real-time EIS data were measured in molten salt and provides the corrosion behavior over time. The corrosion mitigation strategy should consider that the corrosion of Hastelloy N accelerates over time and its intrinsic film cannot act as the protective layer. The next steps of this study are to evaluate other candidate structural materials and to demonstrate the presence of the film.

For safeguarding dry storage facilities, a tomography system based on fast and thermal neutron detection was studied in Korea Institute of Nuclear Nonproliferation and Control. The study was conducted laboratory-scale experiments based on a custom built 1/10th scale model cask, He-4 gas scintillation detector array, and multiple 252Cf sources. A filtered back projection (FBP) was utilized to obtain the cask image via MATLAB. The Ram-Lak filter (ramp filter) was employed in FBP for improved the reconstructed image quality. The Ram-Lak filter is the increasing amplitude filter due to the increasing spatial frequency of the image. In spatial frequency, the frequency of brightness change in the low-frequency region is relatively low, and the frequency of brightness change in the highfrequency region is large. Thus, the high-frequency region in the neutron tomographic image is near the neutron sources and the cask, and the low-frequency region is outside of the cask and/or between the source and cask in the study. In order to apply the ramp filter, a Fourier transform is initially performed on projection data, and image reconstruction is performed with the corrected projection data. In this case, the filter is linearly changed. Therefore, a small filter value is applied at lower spatial frequencies to reduce the projection data, and a large filter value is applied at high spatial frequencies to reduce the projection data. The filter scale is a fraction of frequency amplitude, and the filter value applied to the projection data is determined according to the filter scale. This study was conducted for discussion of the image quality due to the effect of the filter scale used for image reconstruction of a neutron tomography system. The results show that in the experiment with one source, the source location was founded when we used the frequency scale of 0.5 and over. In the double or triple source experiment, the source locations and relative activities were found when we used a filter scale of 0.4 to 0.6. When the filter frequency scale of 0.7 to over, the relative activities are hard to know. It can be found that if the filter value is too large or too small, distortion may occur in the reconstruction results. Therefore, it seems reasonable to set a value between 0.4 and 0.6 as the scaling factor for the neutron tomography system. In the future, additional comparative studies will perform validation of the frequency scaling methods.

The parent and daughter nuclides in a radioactive decay chain arrive at secular equilibrium once they have a large half-life difference. The characteristics of this equilibrium state can be used to estimate the production time of nuclear materials. In this study, a mathematical model and algorithm that can be applied to radio-chronometry using the radioactive equilibrium relationship were investigated, reviewed, and implemented. A Bateman equation that can analyze the decay of radioactive materials over time was used for the mathematical model. To obtain a differential-based solution of the Bateman equation, an algebraic numerical solution approach and two different matrix exponential functions (Moral and Levy) were implemented. The obtained result was compared with those of commonly used algorithms, such as the Chebyshev rational approximation method and WISE Uranium. The experimental analysis confirmed the similarity of the results. However, the Moral method led to an increasing calculation uncertainty once there was a branching decay, so this aspect must be improved. The time period corresponding to the production of nuclear materials or nuclear activity can be estimated using the proposed algorithm when uranium or its daughter nuclides are included in the target materials for nuclear forensics.

Spent nuclear fuel is a very complex material because various elements such as fission products, transuranium elements and activation products are produced from initial fresh UO2 fuel after irradiation. These elements exist in UO2 with various forms and can change the structure and of physicochemical properties of UO2. These changes could provide the surface activation site that could enhance chemical reactions and corrosion processes, and would significantly affect the storage environment for long-term disposal of spent nuclear fuel. Therefore, it can be important to understand the characteristics of spent nuclear fuel to design reliable and safe geological repositories. However, it is too hard to study the characteristics of spent nuclear fuel, because it is a very complex material by itself and not easy to handle due to its radioactivity, and it is also difficult to independently understand the effects of each element. Therefore, a simulated spent nuclear fuel containing an element that forms a solid solution and epsilon particle was manufactured to understand the change in characteristics of each element. Most of the elements that form solid solutions are lanthanides or actinides and can change the structure of the UO2 lattice itself. The epsilon particles exist as metals at the grain boundaries of UO2. In this study, structural changes were measured using XRD, SEM, and Raman spectroscopy, and physical and chemical properties were also identified by measuring electrical conductivity and electrochemical properties. The results were summarized, and the effects of solid solution elements and epsilon particles on the structure and properties of UO2 matrix were compared and discussed.

Complexation of actinides and lanthanides with carboxylic organic ligands is known to facilitate migration of radionuclides from deep geological disposal systems of spent nuclear fuel. In order to examine the ligand-dependent structures of trivalent actinides and lanthanides, a series of Eu(III)-aliphatic dicarboxylate compounds, Eu2(oxalate)3(H2O)6, Eu2(malonate)3(H2O)6, and Eu2(succinate)3(H2O)2, were synthesized and characterized by using X-ray crystallography and time-resolved laser fluorescence spectroscopy. Powder X-ray diffraction results captured the transition of the coordination modes of aliphatic dicarboxylate ligands from side-on to end-on binding as the carbon chain length increases. This transition is illustrated in malonate bindings involving a combination of side-on and end-on modes. Strongly enhanced luminescence of these solid compounds, especially on the hypersensitive peak, indicates a low site symmetry of these solid compounds. Luminescence lifetimes of the compounds were measured to be increased, which is ascribed to the displacement of water molecules in the innersphere of Eu center upon bindings of the organic ligands. The numbers of remaining bound water molecules estimated from the increased luminescence lifetimes were in good agreement with crystal structures. The excitation-emission matrix spectra of these crystalline polymers suggest that oxalate ligands promote the sensitized luminescence of Eu(III), especially in the UV region. In the case of malonate and succinate ligands, charge transfer occurs in the opposite direction from Eu(III) to the ligands under UV excitation, resulting in weaker luminescence.

In this study, an aerosol process was introduced to produce CaCO3. The possibility of producing CaCO3 by the aerosol process was evaluated. The characteristics of CaCO3 prepared by the aerosol process were also evaluated. In the CaCO3 prepared in this study, as the heat treatment proceeded, the calcite phase disappeared. The portlandite phase and the lime phase were formed by the heat treatment. Even if the CO2 component is removed from the calcite phase, there is a possibility that the converted CO2 component could be adsorbed into the Ca component to form a calcite phase again. Therefore, in order to remove the calcite phase, carbon components should be removed first. The lime phase was formed when CO2 was removed from the calcite phase, while the portlandite phase was formed by the introducing of H2O to the lime phase. Therefore, the order in which each phase formed could be in the order of calcite, lime, and portlandite. The reason for the simultaneous presence of the portlandite phase and the lime phase is that the hydroxyl group (OH−) introduced by H2O was not removed completely due to low temperature and/or insufficient heating time. When the sufficient temperature (900°C) and heating time (60 min) were applied, the hydroxyl group (OH−) was removed to transform into lime phase. Since the precursor contained the hydrogen component, it could be possible that the moisture (H2O) and/or the hydroxyl group (OH−) were introduced during the heat treatment process.

Uranium-235, used for nuclear power generation, has brought radioactive waste. It could be released into the environment during reprocessing or recycling of the spent nuclear fuel. Among the radioactive waste nuclides, I-129 occurs problems due to its long half-life (1.57×107 y) with high mobility in the environment. Therefore, it should be captured and immobilized into a geological disposal system through a stable waste form. One of the methods to capture iodine in the off-gas treatment process is to use silver loaded zeolite filter. It converts radioactive iodine into AgI, one of the most stable iodine forms in the solid state. However, it is difficult to directly dispose of AgI itself in an underground repository because of its aqueous dissolution under reducing condition with Fe2+. It must be immobilized in the matrix materials to prevent release of iodine as a result of chemical reaction. Among the matrix glasses, silver tellurite glass has been proposed. In this study, additives including Al, Bi, Pb, V, Mo, and W were added into the silver tellurite glass. The thermal properties of each matrix for radioactive iodine immobilization were evaluated. The glasses were prepared by the melt-quenching method at 800°C for 1 h. Differential scanning calorimetry (DSC) was performed to evaluate the thermal properties of the glass samples. From the study, the glass transition temperature (Tg) was increased by adding additives such as V2O5, MoO3, or WO3 in the silver tellurite glass. The relative electro-static field (REF) values of V2O5, MoO3, and WO3 are about three times higher than that of the glass network former, TeO2. It could provide sufficient electro-static field (EF) to the TeO2 interacting with the non-bridging oxygen forming Te-O-M (M = V, Mo, W) links. Therefore, the addition of V2O5, MoO3, or WO3 reinforced the glass network cohesion to increase the Tg of the glass. The addition of MoO3or WO3 in the silver tellurite glass increased Tg and crystallization temperature (Tc) with remaining the glass stability.

In a recent preliminary inspection for disposal, the glass fiber waste (GFW), used as a pipe insulation, was judged as “pending evaluation” because some dust was found in drum opening tests. Therefore, additional inspection is required to ensure that the package corresponds with the acceptance criteria of the particulates. The dust was generated presumably due to GFW being used in a high-temperature environment for a long time, thus being easily degraded and crushed. For this reason, safety issues that may occur in the process of handling, transportation, and disposal are emerging. Therefore, in this study, a preliminary safety assessment of GFW disposal was performed, the exposure dose to the general public was derived, and compared with the dose limit. The evaluation was carried out in the following order: (1) evaluation of GFW radiation source term, (2) selection of accident scenario, (3) calculation of exposure dose, (4) comparison of evaluation results with dose limits, and confirmation of satisfaction. The average radioactivity of the GFW to be disposed of was used as the source term, and the main nuclides were identified as H-3, Fe-55, Co-60, Ni-63, and Pu-241. In general, the types of accidents that can occur at disposal facilities can be classified into falls, fires, collisions during transportation, off-site accidents, and nuclear criticality, and the accident scenarios are selected by analyzing and reviewing the probability of each accident. In this study, the accident analysis and scenarios presented in the safety assessment of the KORAD were reviewed, and the fire in the treatment facility, the fire in the storage facility, and the collision of the transport vehicle were selected as the evaluation scenarios. When an accident occurs, the radioactive material inside the container leaks out and diffuses into the atmosphere. In this evaluation, the internal and external exposure of the general public due to radioactive plume at the site boundary was evaluated and the dose conversion factors from ICRP-72 and FGR 12 were used. Based on the evaluation, general public was exposed to 0.004 mSv, 0.013 mSv, and 0.045 mSv, respectively, due to a fire at a treatment facility, at a storage facility, and in a transport vehicle. Most of the dose is due to internal exposure by Pu-241 nuclide, because the proportion of it in the waste is high, and when inhaled, the internal dose is high by emitting beta rays. It was confirmed that the result of dose was 0.4%, 1.3% and 4.5% of the annual dose limit, sufficiently satisfying the dose limit and safety.

Recently, concern regarding disposal of cellulosic material is growing as cellulose is known to produce complexing agent, isosaccharinic acid (ISA), upon degradation. ISA could enhance mobility of some radionuclides, thus increasing the amount of radionuclide released into the environment. Thus, evaluation on the possible impact of the cellulose degradation would be an important aspect in safety evaluation. In this paper, safety assessments conducted in Sweden and UK are studied, and the factors required to be considered for appropriate safety assessment of cellulose is analyzed. SKB (Sweden) conducted safety assessment of cellulose degradation as a part of long-term safety assessment of SFR. SKB determined that ISA would impact sorption of trivalent and tetravalent radionuclides (Eu, Am, Th, Np, Pa, Pu, U, Tc, Zr and Nb) at concentration higher than 10−4–10−3 M, and impact sorption of divalent radionuclides (Ni, Co, Fe, Be and Pb) at concentration higher than 10−2 M. Then, SKB conservatively set the upper limit of ISA concentration to be 10−4 M and conducted cellulose degradation evaluation on each waste package type, considering the expected disposal environment of SFR. Based on the calculated results, some of the waste packages showed concentration of ISA to be higher than 10−4 M, so SKB conservatively developed waste acceptance criteria to prevent ISA being produced to an extent of affecting the safety of the repository. SKB conducted safety assessment only for the repositories with pH above 12.5 and excluded 1BLA from the safety assessment as the expected pH of 1BLA is around 12, which is insufficient for cellulose to degrade. However, SKB set disposal limit for 1BLA as well, to minimize potential impact in future. Serco (UK) conducted safety assessment of cellulose degradation for the conceptual repository, which is a concrete vault with cementitious backfill. Serco estimated that the pH of repository would maintain around 12.4. Serco conservatively assumed that the pH would be sufficient for cellulose degradation to occur partially, and suggested application of appropriate degradation ratio for safety assessment of cellulose degradation. To conduct appropriate safety assessment of cellulose degradation, an appropriate ISA concentration limit based on radionuclide inventory list, and an appropriate cellulose degradation ratio based on the pH of disposal environment should be determined. As for guidance, below pH 12.5, cellulose degradation is not expected, and between pH 12.5–13, partial cellulose degradation is expected. In future, this study could be used as fundamental data to evaluate safety of the repository.

Near-surface disposal facility is more susceptible to intrusion than underground repository, resulting in more possible pathways for contaminant release. Alike human intrusion, animals (e.g. Ants, Moles, etc.) could intrude into the disposal site to excavate burrows, which could cause direct release of contaminants to biosphere. In this paper, animal intrusion is demonstrated using GoldSim’s commercial contaminant transport module and impact on the integrity of the near-surface disposal facility is evaluated in terms of fractional release rate of the contaminants. In this study, the near-surface disposal facility is modelled with a single concrete vault to contain radionuclide according to LLW concentration limit stated in NSSC notice No.2020-6. The release of contaminants is modelled to occur directly after the institutional control period, and the contaminants are mostly transported from the concrete vault to cover layers via diffusion. To produce mathematical model of the release of the contaminants due to animal intrusion, firstly, the fraction of burrow volume for each cover layer is calculated separately for each animal species, based on their maximum possible intrusion depth. In this study, fractions of burrow volume for ants and moles are calculated based on their maximum possible intrusion depths, where for ants is 2–3 m, and for moles is 0.1–0.135 m. Then, assuming that the contaminants are distributed homogeneously throughout each cover layers by diffusion, fraction of contaminants transported into the uppermost layer via excavation of the burrow is calculated for each layer based on burrow volume, and fraction of contaminants removed from the uppermost layer to the layers below via collapse of the burrow is also calculated based on the burrow volume. Lastly, the net transportation of contaminants into and out of the burrow via excavation and collapse, respectively, is calculated and demonstrated using direct transfer rate function of the GoldSim. Based on the simulated result, the maximum mass flux is too minor to cause a meaningful impact on the safety. The peak mass flux of the most sensitive radionuclide, I-129, is witnessed at around year 1,470, with a flux value of 5.36×10−6 g·yr−1. This minor release of the contaminants could be due to cover layers being much thicker than the maximum possible intrusion depth of the animals, preventing the animal intrusion into the deeper layers of higher radionuclide concentration. In future, this study can be used to provide a guidance and fundamental data for scenario development and safety evaluation of the near-surface disposal facility.

Recently, concern regarding disposal of cellulosic material is growing as cellulose is known to produce complexing agent, isosaccharinic acid (ISA), upon degradation. ISA could enhance mobility of some radionuclides, thus increasing the amount of radionuclide released into the environment. Evaluation on the possible impact of the cellulose degradation would be an important aspect in safety evaluation. In this paper, the maximum safe disposal amount cellulose is evaluated considering the disposal environment of silos of 1st phase disposal facility. The key factor governing the impact of cellulose degradation is pH of disposal environment, as cellulose is known to degrade partially at pH above 12.5, and completely at pH above 13. Thus, disposal environment should be analyzed as to determine the extent of degradation. As silos are constructed with large amount of cement, porewater within concrete walls would be of very high pH. However, for high pH porewater to be released into the pores of crushed rock, which is filling up the silos, lower pH groundwater (commonly pH 7) should flow into the silos through the concrete walls. This causes dilution of the high pH concrete porewater, resulting in a lower pH as the silos are filled, reaching to expected pH of 11.8–12.3, which is below cellulose degradation condition. Thus, cellulose degradation is not expected, but to quantitatively evaluate safe disposal amount of cellulose, partial degradation is assumed. Upon literature review, the most conservative ISA concentration, enhancing radionuclide mobility, is determined to be 1.0×10−4 M and to reach this concentration, cellulose mass equivalent to 6wt% of cement of the repository, is required to be degraded. However, this ratio is derived based on complete degradation of cellulose into ISA, so for partial degradation, degradation ratio and yield ratio of ISA should be considered. Commonly, cellulosic material (e.g. cotton, paper, etc.) has degree of polymerization (DP) between 1,000–2,000, and with this DP, degradation ratio is estimated to be about 10%. Furthermore, yield ratio of ISA is known to be 80%. Considering all these aspects, about 1.79×107 kg of cellulose could be disposed, which if converted into number of drums, considering cellulose content of dry active waste, more than 100,000 drums (200 L) could be disposed with negligible impact on safety. Based on the result, negligible impact of cellulose degradation is expected for safety of 1st phase disposal facility. In future, this study could be used as fundamental data for revising waste acceptance criteria.

Deep geological disposal (DGD) of spent nuclear fuels (SNF) at 500 m–1 km depth has been the mainly researched as SNF disposal method, but with the recent drilling technology development, interest in deep borehole disposal (DBD) at 5 km depth is increasing. In DBD, up to 40SNF canisters are disposed of in a borehole with a diameter of about 50 cm, and SNF is disposed of at a depth of 2–5 km underground. DBD has the advantage of minimizing the disposal area and safely isolating highlevel waste from the ecosystem. Recently, due to an increasing necessity of developing an efficient alternative disposal system compared to DGD domestically, technological development for DBD has begun. In this paper, the research status of canister operation technology and plans for DBD demonstration tests, which subjects are being studied in the project of developing a safety-enhancing high-efficiency disposal system, are introduced. The canister operation technology for DBD can be divided into connection device development and operation technology. The developing connection device, emplacing and retrieving canisters in borehole, adopted the concept of a wedge thus making replacement equipment at the surface unnecessary. The new connection device has the advantage of being well applied with emplacement facilities only by simple mechanical operation. The technology of operating a connection device in DBD can be divided into drill pipe, coiled tubing, free-drop, and wireline. The drill pipe is a proven method in the oil industry, but requiring huge surface equipment. The coiled tubing method uses a flexible tube and shares disadvantages as the drill pipe. The free-drop is a convenient method of dropping canister into a borehole, but has a weakness in irretrievability in an accident. Finally, the wireline method can be operational on a small scale using hydraulic cranes, but the number of operated canisters at once is limited. The test facility through which the connection device is to be tested consists of dummy canister, borehole, lifting part, monitoring part, and connecting device. The canister weight is determined according to the emplacement operation unit. The lifting part will be composed following wireline consisting of a crane, a wire and a winding system. The monitoring part will consist of an external monitoring system for hoists and trolleys, and an internal monitoring system for the connection device’s location, pressure, and speed. In this project, a demonstration test will be conducted in a borehole with 1km depth, 10 cm diameter provided by KAERI to verify operation in the actual drilling environment after design improvement of the connecting device. If a problem is found through the demonstration test, the problem will be improved, and an improved connection device will be tested to an extended borehole with a 2 km depth, 40 cm diameter.

The design of the high-level radioactive waste (HLW) repository is made for isolating the HLW from the groundwater system by using artificial and natural barriers. Granite is usually considered to be a great natural barrier for the HLW repository in various countries including Sweden, Canada, and Korea due to its low hydraulic permeability. However, many fractures that can act as conduits for groundwater and radionuclides exist in granite. Furthermore, the decay heat generated by the HLW can induce groundwater acceleration through the fracture. Since the direction, magnitude, and lasting time of the heat-induced groundwater flow can be differed depending on the fracture geometry, the effect of fracture geometry on the groundwater flow around the repository should be carefully analyzed. In this study, groundwater models were conducted with various fracture geometries to quantify the effect of various properties of fractures (or fracture networks) on the heat-induced groundwater flow. In all models, the pressure around the repository only lasted for a short period after it peaked at 0.1 years. In contrast, the temperature lasted for 10,000 years after the disposal inducing the convective groundwater flow. Single fracture models with different orientations were conducted to evaluate the variations in groundwater velocities around the repository depending on the fracture slope. According to the results, the groundwater velocity on the fracture was the fastest when the regional groundwater flow direction and the fracture direction coincided. In double fracture models, various inclined fractures were added to the horizontal fracture. Due to the intersecting, the groundwater flow velocity showed a discontinuous change at the intersecting point. Lastly, the discrete fracture network models were conducted with different fracture densities, length distributions, and orientations. According to the modeling results, the groundwater flow was significantly accelerated when the fracture network density increased, or the average fracture length increased. However, the effect of the fracture orientation was not significant compared to the other two network properties.