• Proceedings of the Korean Radioactive Waste Society Conference
• /
• 2009.06a
• /
• pp.187-187
• /
• 2009

Opening of fractures induced by shear dilation or normal deformation can be a significant source of fracture permeability change in fractured rock, which is important for the performance assessment of geological repositories for spent nuclear fuel. As the repository generates heat and later cools the fluid-carrying ability of the rocks becomes a dynamic variable during the lifespan of the repository. Heating causes expansion of the rock close to the repository and, at the same time, contraction close to the surface. During the cooling phase of the repository, the opposite takes place. Heating and cooling together with the, virgin stress can induce shear dilation of fractures and deformation zones and change the flow field around the repository. The objectives of this work are to examine the contribution of thermal stress to the shear slip of fracture in mid- and far-field around a KBS-3 type of repository and to investigate the effect of evolution of stress on the rock mass permeability. In the first part of this study, zones of fracture shear slip were examined by conducting a three-dimensional, thermo-mechanical analysis of a spent fuel repository model in the size of 2 km $\times$ 2 km $\times$ 800 m. Stress evolutions of importance for fracture shear slip are: (1) comparatively high horizontal compressive thermal stress at the repository level, (2) generation of vertical tensile thermal stress right above the repository, (3) horizontal tensile stress near the surface, which can induce tensile failure, and generation of shear stresses at the comers of the repository. In the second part of the study, fracture data from Forsmark, Sweden is used to establish fracture network models (DFN). Stress paths obtained from the thermo-mechanical analysis were used as boundary conditions in DFN-DEM (Discrete Element Method) analysis of six DFN models at the repository level. Increases of permeability up to a factor of four were observed during thermal loading history and shear dilation of fractures was not recovered after cooling of the repository. An understanding of the stress path and potential areas of slip induced shear dilation and related permeability changes during the lifetime of a repository for spent nuclear fuel is of utmost importance for analysing long-term safety. The result of this study will assist in identifying critical areas around a repository where fracture shear slip is likely to develop. The presentation also includes a brief introduction to the ongoing site investigation on two candidate sites for geological repository in Sweden.

• Tunnel and Underground Space
• /
• v.31 no.4
• /
• pp.270-288
• /
• 2021

We proposed a numerical method to simulate the hydro-mechanical behavior of rock fracture using a grain-based distinct element model (GBDEM) in the paper. As a part of DECOVALEX-2023 Task G, we verified the method via benchmarks with analytical solutions. DECOVALEX-2023 Task G aims to develop a numerical method to estimate the coupled thermo-hydro-mechanical processes within the crystalline rock fracture network. We represented the rock sample as a group of tetrahedral grains and calculated the interaction of the grains and their interfaces using 3DEC. The micro-parameters of the grains and interfaces were determined by a new methodology based on an equivalent continuum approach. In benchmark modeling, a single fracture embedded in the rock was examined for the effects of fracture inclination and roughness, the boundary stress condition and the applied pressure. The simulation results showed that the developed numerical model reasonably reproduced the fracture slip induced by boundary stress condition, the fracture opening induced by fluid injection, the stress distribution variation with fracture inclination, and the fracture roughness effect. In addition, the fracture displacements associated with the opening and slip showed good agreement with the analytical solutions. We expect the numerical model to be enhanced by continuing collaboration and interaction with other research teams of DECOVALEX-2023 Task G and validated in further study experiments.

• Advances in biomechanics and applications
• /
• v.1 no.2
• /
• pp.127-141
• /
• 2014

A preliminary study of a new pulsatile pump that will work to a frequency greater than 1 Hz, is presented. The fluid-structure interaction between a Newtonian blood flow and a piston drive that moves with periodic speed is simulated. The mechanism is of double effect and has four valves, two at the input flow and two at the output flow; the valves are simulated with specified velocity of closing and reopening. The simulation is made with finite elements software named COMSOL Multiphysics 3.3 to resolve the flow in a preliminary planar configuration. The geometry is 2D to determine areas of high speeds and high shear stresses that can cause hemolysis and platelet aggregation. The opening and closing valves are modelled by solid structure interacting with flow, the rhythmic opening and closing are synchronized with the piston harmonic movement. The boundary conditions at the input and output areas are only normal traction with reference pressure. On the other hand, the fluid structure interactions are manifested due to the non-slip boundary conditions over the piston moving surfaces, moving valve contours and fix pump walls. The non-physiologic frequency pulsatile pump, from the viewpoint of fluid flow analysis, is predicted feasible and with characteristic of low hemolysis and low thrombogenesis, because the stress tension and resident time are smaller than the limit and the vortices are destroyed for the periodic flow.

• Journal of Mechanical Science and Technology
• /
• v.14 no.1
• /
• pp.84-92
• /
• 2000

A new bounce back boundary method of the second order in error is proposed for the lattice Boltzmann fluid simulation. This new method can be used for the arbitrarily irregular lattice geometry of a non-slip boundary. The traditional bounce back boundary condition for the lattice Boltzmann simulation is of the first order in error. Since the lattice Boltzmann method is the second order scheme by itself, a boundary technique of the second order has been desired to replace the first order bounce back method. This study shows that, contrary to the common belief that the bounce back boundary condition is unilaterally of the first order, the second order bounce back boundary condition can be realized. This study also shows that there exists a generalized bounce back technique that can be characterized by a single interpolation parameter. The second order bounce back method can be obtained by proper selection of this parameter in accordance with the detailed lattice geometry of the boundary. For an illustrative purpose, the transient Couette and the plane Poiseuille flows are solved by the lattice Boltzmann simulation with various boundary conditions. The results show that the generalized bounce back method yields the second order behavior in the error of the solution, provided that the interpolation parameter is properly selected. Coupled with its intuitive nature and the ease of implementation, the bounce back method can be as good as any second order boundary method.

• International Journal of Fluid Machinery and Systems
• /
• v.2 no.1
• /
• pp.80-91
• /
• 2009

In this study, a prediction theory for specific noise that is the overall characteristic of the fan has been proposed. This theory is based on total pressure prediction and broadband noise prediction. The specific noises of two forward curved fans with different number of blades were predicted. The flow around the impeller having 120 blades (MF120) was more biased at a certain positions than the impeller with 40 blades (MF40). An effective domain of the energy conversion of MF40 has extended overall than MF120. The total pressure was affected by the slip factor and pressure loss caused by the vortex flow. The suppression of a major pressure drop by the vortex flow and expansion of the effective domain for energy conversion contributed to an increase in the total pressure of MF40 at the design point. The position of maximum relative velocity was different for each fan. The relative velocity of MF120 was less than that of MF40 due to the deviation angle. The specific noise of MF120 was 2.7 dB less than that of MF40 due to the difference in internal flow. It has been quantitatively estimated that the deceleration in the relative velocity contributed to the improvement in the overall performance.

• The KSFM Journal of Fluid Machinery
• /
• v.10 no.1 s.40
• /
• pp.56-63
• /
• 2007

Because the mixing efficiency is influenced remarkably by varying the geometrical configurations, the study of flow characteristics inside the mechanical agitator is very important to improve the performances. The draught tube in the agitator makes intermixing between the screw and tube by interrupting radial flow, and it makes circulation region in a mixing chamber. In general, the helical screw agitator with a draught tube (HSA) is proved more efficient to mix than the others. Consequently, such as the shapes of helical screw, number of pitches and the variation of angular velocity are the main parameters for improving the capacity of HSA. And also the suspension of the solid particles in the agitator can be determined these parameters. The rate of solids suspension in the mixing chamber was quantified with a statistical average value, of. Numerical analyses were carried out, using a commercial CFD code, Fluent, to obtain the velocity, pressure and particle distributions under steady, laminar flow and no-slip conditions. Results are graphically depicted with various parameters.

• Journal of the Korean Society for Precision Engineering
• /
• v.2 no.2
• /
• pp.69-75
• /
• 1985

The object of this study is to discuss the effect of cutting fluid on the mechanism of chip formation in orthogonal cutting. Rehbinder effect has been known as a phenomenon, the reduction of mechanical strength, when the metal is exposed in a polar organic environment or the surface of metal is coated with some polar organic substances. About the cause of Rehbinder effect there have been many different ideas by Rehbinder, Merchant, Shaw, Sakakida and etc. In this report, the effects of surface active medium (magic ink) upon the mechanism of chip formation on the orthogonal cutting of copper and the mechanical properties of the work material are experimentally discussed with constant rake angle. Under the condition of polar organic environment the experimental results are as follows; 1) The chip thickness becomes thinner and slip line pitch on the free surface of chip becomes smaller than that of dried cutting area. 2) The order of alternation of cutting ratio was changed. 3) The friction angle on the tool face is not affected by the depth of cut. 4) The cutting force and shear strain on the shear plane decrease remarkably, therefore the work material must be embrittled under polar organic environment.

• Proceedings of the Korean Society of Precision Engineering Conference
• /
• 2004.10a
• /
• pp.496-500
• /
• 2004

A new, commercially available polishing process called magnetorheological finishing is used to polish and figure precision optics. To understand and model this process correctly it is important to determine the mechanical properties of the fluid under the influence of the magnetic field. Magnetorheological (MR) fluids are commonly modeled as Bingham fluids, so one of the essential properties to measure is the yield stress. Since MR fluids are inherently anisotropic, the yield stress will depend on the mutual orientation of the magnetic field and the direction of deformation. The relative orientation of the field and deformation in polishing does not coincide with common rheological setups, so a new rheometer has been designed and tested. This new magnetorheometer design has been shown to give correct stresses during calibration experiments using Newtonian fluids with a known viscosity. The measured stress has also been shown to have a magnitude consistent with published finite element approximations for magnetic fluids. The design of the instrument was complicated because of the requirements imposed upon the magnetic field, and the difficulty in satisfying the no slip boundary condition. Our results show the importance of having a homogeneous field in the test region during measurements. The solutions to these problems and discussion of the measurements on nonmagnetic and magnetic fluids are given.

• The KSFM Journal of Fluid Machinery
• /
• v.17 no.6
• /
• pp.13-20
• /
• 2014

It is necessary to design a compressor with high pressure ratio that satisfies the IMO(international maritime organization) NOx emission regulation for the marine diesel engine. Impeller was designed using the modified slip factor with the flow coefficient. The main purpose of this study is to investigate the sensitivity of the compressor performance by the vaned diffuser geometries. The first vaned diffuser type was based on a NACA airfoil, the second was channel diffuser, and the third was conformally transformated configuration of a NACA65(4A10)06 airfoil. The sensitivity of the performance was calculated using a commercial CFD program for three different diffuser geometries. The channel diffuser showed the wide range of operation and higher pressure characteristics, comparing with the others. This is attributed to the flow stability at diffuser. Combined with this results with impeller design, the optimized compressor was designed and verified by the test results.

• 유체기계공업학회:학술대회논문집
• /
• 2005.12a
• /
• pp.448-453
• /
• 2005

Temperature Control Valve (TCV) is one of the useful temperature control devices, which is used to control constant temperature of working fluid in power and chemical plants and domestic water supply systems. TCV is composed of body, cylinder and piston, and the body shape has a symmetrical H-type. In general, it has several inlet and outlet holes, and its shape is like as tubular sleeve. The piston has three rings two rings of the end of piston have the function of controlling inlet flow rate with hot and cold working fluids, the center ring has the function of preventing hot and cold water from intermixing. Consequently, the shapes of piston and cylinder are the main design parameters in the performance of TCV. In this study, numerical analyses were carried out with two different piston and cylinder shapes to investigate the functions as a temperature control valve and the flow characteristics according to piston opening grade in TCV. Using a commercial code, FLUENT, velocity and pressure fields in TCV are obtained under steady, standard $k -{\epsilon}$ turbulence model and no-slip condition.