• Nuclear Engineering and Technology
• /
• v.42 no.5
• /
• pp.500-519
• /
• 2010

AEGIS is a lattice physics code incorporating the latest advances in lattice physics computation, innovative calculation models and efficient numerical algorithms and is mainly used for light water reactor analyses. Though the primary objective of the AEGIS code is the preparation of a cross section set for SCOPE2 that is a three-dimensional pin-by-pin core analysis code, the AEGIS code can handle not only a fuel assembly but also multi-assemblies and a whole core geometry in two-dimensional geometry. The present paper summarizes the major calculation models and part of the verification/validation efforts related to the AEGIS code.

• International Journal of Contents
• /
• v.7 no.4
• /
• pp.1-5
• /
• 2011

Physics based fluid animation schemes need large computation cost due to tremendous degree of freedom. Many researchers tried to reduce the cost for solving the large linear system that is involved in grid-based schemes. GPU based algorithms and advanced numerical analysis methods are used to efficiently solve the system. Other groups studied local operation methods such as SPH (Smoothed Particle Hydrodynamics) and LBM (Lattice Boltzmann Method) for enhancing the efficiency. Our method investigates this efficiency problem thoroughly, and suggests novel paradigm in fluid animation field. Rather than physics based simulation, we propose a robust boundary handling technique for procedural fluid animation. Our method can be applied to arbitrary shaped objects and potential fields. Since only local operations are involved in our method, parallel computing can be easily implemented.

• Proceedings of the Korean Vacuum Society Conference
• /
• 2014.02a
• /
• pp.185.2-185.2
• /
• 2014

We have studied the atomic and electronic structure of graphene nanoribbons (GNRs) on a hexagonal boron nitride (h-BN) sheet with intercalated atoms using first-principles calculations. The h-BN sheet is an insulator with the band gap about 6 eV and then it may a good candidate as a supporting dielectric substrate for graphene-based nanodevices. Especially, the h-BN sheet has the similar bond structure as graphene with a slightly longer lattice constant. For the computation, we use the Vienna ab initio simulation package (VASP). The generalized gradient approximation (GGA) in the form of the PBE-type parameterization is employed. The ions are described via the projector augmented wave potentials, and the cutoff energy for the plane-wave basis is set to 400 eV. To include weak van der Waals (vdW) interactions, we adopt the Grimme's DFT-D2 vdW correction based on a semi-empirical GGA-type theory. Our calculations reveal that the localized states appear at the zigzag edge of the GNR on the h-BN sheet due to the flat band of the zigzag edge at the Fermi level and the localized states rapidly decay into the bulk. The open-edged graphene with a large corrugation allows some space between graphene and h-BN sheet. Therefore, atoms or molecules can be intercalated between them. We have considered various types of atoms for intercalation. The atoms are initially placed at the edge of the GNR or inserted in between GNR and h-BN sheet to find the effect of intercalated atoms on the atomic and electronic structure of graphene. We find that the impurity atoms at the edge of GNR are more stable than in between GNR and h-BN sheet for all cases considered. The nickel atom has the lowest energy difference of ~0.2 eV, which means that it is relatively easy to intercalate the Ni atom in this structure. Finally, the magnetic properties of intercalated atoms between GNR and h-BN sheet are investigated.

• Nuclear Engineering and Technology
• /
• v.54 no.7
• /
• pp.2670-2689
• /
• 2022

STREAM - a lattice transport calculation code with method of characteristics for the purpose of light water reactor analysis - has been developed by the Computational Reactor Physics and Experiment laboratory (CORE) of the Ulsan National Institute of Science and Technology (UNIST). Recently, efforts have been taken to develop a photon module in STREAM to assess photon heating and the influence of gamma photon transport on power distributions, as only neutron transport was considered in previous STREAM versions. A multi-group photon library is produced for STREAM based on the ENDF/B-VII.1 library with the use of the library-processing code NJOY. The developed photon solver for the computation of 2D and 3D distributions of photon flux and energy deposition is based on the method of characteristics like the neutron solver. The photon library and photon module produced and implemented for STREAM are verified on VERA pin and assembly problems by comparison with the Monte Carlo code MCS - also developed at UNIST. A short analysis of the impact of photon transport during depletion and thermal hydraulics feedback is presented for a 2D core also from the VERA benchmark.

• Nuclear Engineering and Technology
• /
• v.52 no.10
• /
• pp.2162-2172
• /
• 2020

Due to its computational efficiency and geometrical flexibility, the Method of Characteristics (MOC) has been widely used for light water reactor lattice physics analysis. Usually acceleration methods are necessary for MOC to achieve acceptable convergence on practical reactor physics problems. Among them, Coarse Mesh Finite Difference (CMFD) is very popular and can drastically reduce the number of transport iterations. In OpenMOC, CMFD acceleration was implemented but had the limitation of supporting only a uniform CMFD mesh, which would often lead to splitting MOC source regions, thus creating an unnecessary increase in computation and memory use. In this study, CMFD acceleration with a non-uniform Cartesian mesh is implemented into OpenMOC. We also propose a quadratic fit based CMFD prolongation method in the axial direction to further improve the acceleration when multiple MOC source regions are contained in one CMFD coarse mesh. Numerical results are presented to demonstrate the improvement of the CMFD acceleration capability in OpenMOC in terms of both efficiency and stability.

• Progress in Medical Physics
• /
• v.16 no.1
• /
• pp.24-31
• /
• 2005

The stereotactic radiosurgery (SRS) describes a method of delivering a high dose of radiation to a small tar-get volume in the brain, generally in a single fraction, while the dose delivered to the surrounding normal tissue should be minimized. To perform automatic plan of the SRS, a new method of multi-isocenter/shot linear accelerator (linac) and gamma knife (GK) radiosurgery treatment plan was developed, based on a physical lattice structure in target. The optimal radiosurgical plan had been constructed by many beam parameters in a linear accelerator or gamma knife-based radiation therapy. In this work, an isocenter/shot was modeled as a sphere, which is equal to the circular collimator/helmet hole size because the dimension of the 50% isodose level in the dose profile is similar to its size. In a computer-aided system, it accomplished first an automatic arrangement of multi-isocenter/shot considering two parameters such as positions and collimator/helmet sizes for each isocenter/shot. Simultaneously, an irregularly shaped target was approximated by cubic structures through computation of voxel units. The treatment planning method by the technique was evaluated as a dose distribution by dose volume histograms, dose conformity, and dose homogeneity to targets. For irregularly shaped targets, the new method performed optimal multi-isocenter packing, and it only took a few seconds in a computer-aided system. The targets were included in a more than 50% isodose curve. The dose conformity was ordinarily acceptable levels and the dose homogeneity was always less than 2.0, satisfying for various targets referred to Radiation Therapy Oncology Group (RTOG) SRS criteria. In conclusion, this approach by physical lattice structure could be a useful radiosurgical plan without restrictions in the various tumor shapes and the different modality techniques such as linac and GK for SRS.