• The Journal of Engineering Geology
• /
• v.32 no.1
• /
• pp.41-57
• /
• 2022

A workflow is presented to estimate the size of a representative elementary volume and 3-D hydraulic conductivity tensor based on fluid flow analysis for a discrete fracture network (DFN). A case study is considered for a Cretaceous granitic rock mass at Gijang in Busan, Korea. The intensity and size of joints were calibrated using the first invariant of the fracture tensor for the 2-D DFN of the study area. Effective hydraulic apertures were obtained by analyzing the results of field packer tests. The representative elementary volume of the 2-D DFN was determined to be 20 m square by investigating the variations in the directional hydraulic conductivity for blocks of different sizes. The directional hydraulic conductivities calculated from the 2-D DFN exhibited strong anisotropy related to the hydraulic behavior of the study area. The 3-D hydraulic conductivity tensor for the fractured rock mass of the study area was estimated from the directional block conductivities of the 2-D DFN blocks generated for various directions in 3-D. The orientations of the principal components of the 3-D hydraulic conductivity tensor were found to be identical to those of delineated joint sets in the study area.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.20 no.8
• /
• pp.2593-2600
• /
• 1996

A modified k-$\varepsilon$ turbulence model having a generality is proposed in the present study, in which the constant $C_{\varepsilon2}$in the $\varepsilon$-equation is simply changed as a functional form of a new parameter both satisfying the tensor invariant condition and representing the extra straining effect on complex shear flows. With this model turbulent shear flows over two-dimensional obstacles placed in a channel are numerically studied for different blockage ratios and aspect ratios. Comparing with the available experimental data, the predicted results with the present model provide definite improvements over the standard model's results and work fairly well with the experimental data on the size of the recirculation zone, as well as mean velocity, wall static pressure, turbulent kinetic energy and Reynolds stresses.

• Journal of computational fluids engineering
• /
• v.6 no.4
• /
• pp.15-25
• /
• 2001

Numerical prediction of the diffusion controlled transition in a turbine gas pass is important because it can change the local heat transfer rate over a turbine blade as much as three times. In this study, the gas flow over turbine blade is simplified to the flat plate boundary layer, and an adaptive grid scheme redistributing grid points within the computation domain is proposed with a great emphasis on the construction of the grid control function. The function is sensitized to the second invariant of the mean strain tensor, its spatial gradient, and the interaction of pressure gradient and flow deformation. The transition process is assumed to be described with a κ-ε turbulence model. An elliptic solver is employed to integrate governing equations. Numerical results show that the proposed adaptive grid scheme is very effective in obtaining grid independent numerical solution with a very low grid number. It is expected that present scheme is helpful in predicting actual flow within a turbine to improve computation efficiency.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.18 no.11
• /
• pp.3039-3045
• /
• 1994

The k-$\varepsilon$ turbulence models by Launder et al.(1977, LPS) and Leschziner and Rodi(1981, LR) are modified to account for the secondary straining effect with having a generality in the present paper. The modified models are obtained by replacing the gradient Richardson number used to account for the secondary straining effect in the original models by a new parameter with a tensor-invariant correction form. These two modified models are used to predict the turbulent flow over a backward-facing step. In contrast to both standard and modified LR models, the modified LPS model is found to predict the reattachment point fairy well, as well as mean velocity, wall static pressure, turbulent kinetic energy and Reynolds shear stress in the recirculating region.

• Proceedings of the Korean Society of Precision Engineering Conference
• /
• 1995.10a
• /
• pp.736-739
• /
• 1995

The meshless adaptive method based on multiple scale analysis is developed to simulate large deformation problems. In the procedure, new particles are simply added to the orginal particle distribution because meshless methods do not require mesh structures in the formulations. The high scale component of the approximated solution detects the localized region where a refinement is needed. The high scale component of the second invariant od Green-Lagrangian strain tensor is suggested as the new high gradient detector for adaptive procedures. The feasibility of the proposed theory is demonstrated by a numerical experiment for the large deformation of hyperelastic materials.

• Journal of the Korean Mathematical Society
• /
• v.47 no.5
• /
• pp.1001-1015
• /
• 2010

Let M be a $C^{\infty}$ real hypersurface in $\mathbb{C}^{n+1}$, $n\;{\geq}\;1$, locally given as the zero locus of a $C^{\infty}$ real valued function r that is defined on a neighborhood of the reference point $P\;{\in}\;M$. For each k = 1,..., n we present a necessary and sufficient condition for there to exist a complex manifold of dimension k through P that is contained in M, assuming the Levi form has rank n - k at P. The problem is to find an integral manifold of the real 1-form $i{\partial}r$ on M whose tangent bundle is invariant under the complex structure tensor J. We present generalized versions of the Frobenius theorem and make use of them to prove the existence of complex submanifolds.

• Transactions of Materials Processing
• /
• v.1 no.1
• /
• pp.33-41
• /
• 1992

Sintered P/M connecting rod is forged to increase density and to satisfy dimensional specifications. Flow of the materials is different form that of wrought materials due to pores in the preform. The Mises yield function was modified to. include the first invariant of stress tensor, and the associated flow rule was derived by applying the normality rule to the yield function. Axisymmetric and plane-strain finite element analyes were carried out for the ring and beam portions of the connecting rod, respectively. The flow of the preform and density change of the analysis are presented in this paper. A load-stroke curve was also presented by superimposing analysis results for the ring and beam portions.

• Geomechanics and Engineering
• /
• v.17 no.4
• /
• pp.385-392
• /
• 2019

The present fractional-order plasticity models for granular soil are mainly established under the triaxial compression condition, due to its difficult in analytically solving the fractional differentiation of the third stress invariant, e.g., Lode's angle. To solve this problem, a three dimensional fractional-order elastoplastic model based on the transformed stress method, which does not rely on the analytical solution of the Lode's angle, is proposed. A nonassociated plastic flow rule is derived by conducting the fractional derivative of the yielding function with respect to the stress tensor in the transformed stress space. All the model parameters can be easily determined by using laboratory test. The performance of this 3D model is then verified by simulating multi series of true triaxial test results of rockfill.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.28 no.6
• /
• pp.672-681
• /
• 2004

Direct numerical simulation was carried out to study the vortical structures of the flow around a wall-mounted cube in a channel at Re=1,000 and Re=3,500 based on cubic height and bulk mean velocity. The cubic obstacle is situated in the entrance region of the channel flow where the boundary layers are developing. Upstream of the obstacle, steady and unsteady laminar horseshoe vortex systems are observed at Re=1,000 and Re=3,500, respectively; the near-wake flow is turbulent in both cases. The flow separates at each leading sharp edge of the cube, and subsequent vortex roll-up is noticed in the corresponding free-shear layer. The vortex shedding from the upper leading edge (upper vortices) and that from the two lateral leading edges (lateral vortices) are both quasi-periodic and their frequencies are computed. The upper and lateral vortices further develop into hairpin and Λ vortices, respectively. A series of instantaneous contours of the second invariant of velocity gradient tensor helps us identify spatial and temporal behaviors of the vortices in detail. The results indicate that the length and time scales of the vortical structures at Re=3,500 are much shorter than those at Re:1,000. Correlations between the upper and lateral vortices are also reported.

• International Journal of Aeronautical and Space Sciences
• /
• v.2 no.2
• /
• pp.65-72
• /
• 2001

This paper demonstrates an explicit-implicit, finite element analysis for linear as well as nonlinear hygrothermal stress problems. Additional features, such as moisture diffusion equation, crack element and virtual crack extension(VCE) method for evaluating J-integral are implemented in this program. The Linear Elastic Fracture Mechanics(LEFM) Theory is employed to estimate the crack driving force under the transient condition for an existing crack. Pores in materials are assumed to be saturated with moisture in the liquid form at the room temperature, which may vaporize as the temperature increases. The vaporization effects on the crack driving force are also studied. The ideal gas equation is employed to estimate the thermodynamic pressure due to vaporization at each time step after solving basic nodal values. A set of field equations governing the time dependent response of porous media are derived from balance laws based on the mixture theory. Darcy's law is assumed for the fluid flow through the porous media. Perzyna's viscoplastic model incorporating the Von-Mises yield criterion are implemented. The Green-Naghdi stress rate is used for the invariant of stress tensor under superposed rigid body motion. Isotropic elements are used for the spatial discretization and an iterative scheme based on the full Newton-Raphson method is used for solving the nonlinear governing equations.