• Journal of Environmental Science International
• /
• v.13 no.2
• /
• pp.175-180
• /
• 2004

A prompt gamma-ray neutron activation (PGNA) system was simulated by the Monte Carlo N-Particle transport code (MCNP-4A) to estimate the level at which the scattered photon fluence rate, the absolute efficiency of the HPGe-detector, the volume of the concrete sample and the $^{35}$ /Cl(n, ${\gamma}$) reaction rate in this sample contribute to the count rate in the NaCl concentration measurement. The n- ${\gamma}$ fluence rates at the ST-2 beam tube exit of the HANARO reactor were used as input data, and the GAMMA-X type HPGe detector was modeled to tally 1.1649 MeV ${\gamma}$ -rays emitted from the $^{35}$ Cl(n, ${\gamma}$) reaction in the concrete sample. For three cylindrical concrete samples of 13.8, 46.8 and 157.1 ㎤ volumes, respectively, the relations between the NaCl weight fractions of 0.1, 1, 2 and 5 % in each of the concrete samples and the 1.1 649 MeV pulses created in the HPGe detector model were studied. As a result, it was found that the count rate at the same NaCl concentration nearly depends on the volume of the samples in a simulated condition of the same NaCl concentration samples, and that the linearities of the NaCl concentration calibration curves were reasonable in the narrow range of the NaCl weight fraction.

• Journal of computational fluids engineering
• /
• v.21 no.2
• /
• pp.32-39
• /
• 2016

The Gradient Biome is a unique and large greenhouse(length 200 m, width 50 m, height:40 m) in which the elements of the weather, such as temperature and humidity, are controlled and reproduced in such a way as to create a continuous gradient from the tropical to frigid zones along specified longitudinal or transvers lines on the earth. One of the main purposes of the Gradient Biome is to observe the possible responses of the ecosystems (mainly plants), which are to be corresponding to each test climate and be introduced in the Biome, to the expected global warming. As one of the expected responses is the shift of the ecosystem(s) toward the region of suitable environment, there should be no artificial obstacles, which can prevent the shift, inside the facility. However, it is important but not so easy to find the ways of how the temperature and humidity in the Biome could be reproduced since the environmental variables tends to be homogeneous. In this paper, numerical simulations were carried out to find the effective control methods for air temperature and humidity inside the real scale Biome. One of the contributed solvers of OpenFOAM, which is an open source physics simulation code, was modified and used for the simulations.

• Journal of the Korean Society of Fisheries and Ocean Technology
• /
• v.46 no.1
• /
• pp.20-31
• /
• 2010

In this research, the submersible fish cage was designed to avoid structural and biological damage during harsh sea conditions. The submersible cage system consists of netting, mooring ropes, a floating collar, floats, sinkers and anchors. Whole elements of the cage were modeled on the mass-spring model. The computer simulations were carried out to investigate the dynamic behavior of the cage and to calculate mooring line tension subjected to tidal currents and waves. As expected, the tension values in the mooring line of the submerged position are 36% less compared to that of the surface cage under the same loading conditions. As the wave was used in combination with the current velocity of 1m/s, the average tensile load for the submerged cage showed 85% of the value for the floating cage. The simulation results provide an improved understanding of the behaviors of the structure and valuable information on the optimized design of the cage system exposed to open ocean environmental factors.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2011.11a
• /
• pp.274-278
• /
• 2011

We present an innovative method of multi-physics application involving energetic materials. We use an Eulerian methodology to address these problems. We have devised a new level set based tracking framework that can elegantly handle large gradients typically found in energetic response of high explosive and metals. Proper constitutive relations are employed to model the transient phases of gas, lliquid, and solid in the high strain rate regime. We use the confined and unconfined rate stick results to validate against the experimental data.

• Fisheries and Aquatic Sciences
• /
• v.14 no.2
• /
• pp.138-147
• /
• 2011

Modeling fishing-gear systems is essential to better understand the factors affecting their movement and for devising strategies to control movement. In this study, we present a generalized mathematical modeling methodology to analyze fishing gear and its various components. Fishing gear can be divided into a finite number of elements that are connected with flexible lines. We use an algorithm to develop a numerical method that calculates precisely the shape and movement of the gear. Fishinggear mathematical models have been used to develop software tools that can design and simulate dynamic movement of novel fishing-gear systems. The tool allowed us to predict the shape and motion of the gear based on changes in operation and gear design parameters. Furthermore, the tool accurately calculated the swept volume of towed gear and the surrounding volume of purse-seine gear. We analyzed the fished volume for trawl and purse-seine gear and proposed a new definition of fishing effort, incorporating the concept of fished space. This method may be useful for quantitative fishery research, which requires a good understanding of the selectivity and efficiency of fishing gear used in surveys.

• Structural Engineering and Mechanics
• /
• v.21 no.5
• /
• pp.539-551
• /
• 2005

Finite elements based on isoparametric formulation are known to suffer spurious stiffness properties and corresponding stress oscillations, even when care is taken to ensure that completeness and continuity requirements are enforced. This occurs frequently when the physics of the problem requires multiple strain components to be defined. This kind of error, commonly known as locking, can be circumvented by using reduced integration techniques to evaluate the element stiffness matrices instead of the full integration that is mathematically prescribed. However, the reduced integration technique itself can have a further drawback - rank deficiency, which physically implies that spurious energy modes (e.g., hourglass modes) are introduced because of reduced integration. Such instability in an existing stiffness matrix is generally detected by means of an eigenvalue test. In this paper we show that a knowledge of the dimension of the solution space spanned by the column vectors of the strain-displacement matrix can be used to identify the instabilities arising in an element due to reduced/selective integration techniques a priori, without having to complete the element stiffness matrix formulation and then test for zero eigenvalues.

• Nuclear Engineering and Technology
• /
• v.52 no.8
• /
• pp.1603-1610
• /
• 2020

Due to ever-growing advancements in computers and relatively easy access to them, many efforts have been made to develop high-fidelity, high-performance, multi-physics tools, which play a crucial role in the design and operation of nuclear reactors. For this purpose in this study, the neutronic Monte Carlo and thermal-hydraulic sub-channel codes entitled MCNP and COBRA-EN, respectively, were applied for external coupling with each other. The coupled code was validated by code-to-code comparison with the internal couplings between MCNP5 and SUBCHANFLOW as well as MCNP6 and CTF. The simulation results of all code systems were in good agreement with each other. Then, as the second problem, the core of the VVER-1000 v446 reactor was simulated by the MCNP4C/COBRA-EN coupled code to measure the capability of the developed code to calculate the neutronic and thermohydraulic parameters of real and industrial cases. The simulation results of VVER-1000 core were compared with FSAR and another numerical solution of this benchmark. The obtained results showed that the ability of the MCNP4C/COBRA-EN code for estimating the neutronic and thermohydraulic parameters was very satisfactory.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.03a
• /
• pp.186-196
• /
• 2008

The dependence of the rheological properties of blood on shape, aggregation, and deformability of red blood cells (RBCs) has been investigated using hybrid systems by coupling fluid with solid models. We present a simple approach for simulating blood as a multi-component fluid, in which RBCs are modeled as droplets of acquired biconcave shape. We used lattice Boltzmann method (LBM) due to its excellent numerical stability as a simulation tool. The model enables us to control the droplet static shape by imposing non-isotropic surface tension force on the interface between the two components. The use of the proposed non-isotropic surface tension method is justified by the Norris hypothesis. This hypothesis states that the shape of the RBC is due to a non-uniform interfacial surface tension force acting on the RBC periphery. This force is caused by the unbalanced distribution of the lipid molecules on the surface of the RBC. We also used the same concept to investigate the dynamic shape change of the RBC while flowing through the microvasculature, and to explore the physics of the Fahraeus, and the Fahraeus-Lindqvist effects.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.10a
• /
• pp.186-196
• /
• 2008

The dependence of the rheological properties of blood on shape, aggregation, and deformability of red blood cells (RBCs) has been investigated using hybrid systems by coupling fluid with solid models. We present a simple approach for simulating blood as a multi-component fluid, in which RBCs are modeled as droplets of acquired biconcave shape. We used lattice Boltzmann method (LBM) due to its excellent numerical stability as a simulation tool. The model enables us to control the droplet static shape by imposing non-isotropic surface tension force on the interface between the two components. The use of the proposed non-isotropic surface tension method is justified by the Norris hypothesis. This hypothesis states that the shape of the RBC is due to a non-uniform interfacial surface tension force acting on the RBC periphery. This force is caused by the unbalanced distribution of the lipid molecules on the surface of the RBC. We also used the same concept to investigate the dynamic shape change of the RBC while flowing through the microvasculature, and to explore the physics of the Fahraeus, and the Fahraeus-Lindqvist effects.

• Journal of the Korean Institute of Telematics and Electronics
• /
• v.27 no.3
• /
• pp.63-71
• /
• 1990

The hydrodynamic model parameters for the submicron GaAs simulation are calculated using the Monte Carlo method. $\Gamma$, L-, and X-valleys are included in the conduction band of GaAs, and polar optic phonon, acoustic phonon, equivalent intervalley, non-equivalent intervalley, ionized impurity, and piezoelectric scattering are taken into account. The velocity-electric field strength curve obtained in this paper is in good agreement with experimental one. We present the results in tabular form so that other participants can make use of them to simulate the submicron GaAs devices by the hydrodynamic model.