• 지질공학
• /
• 제32권1호
• /
• pp.41-57
• /
• 2022

본 연구는 부산 기장지역에 분포하는 백악기 흑운모 화강암에 대한 사례연구를 통하여 DFN(discrete fracture network) 유동해석 기반의 REV(representative elementary volume) 크기 및 삼차원 수리전도텐서 추정을 위한 워크플로우를 제시하였다. 연구지역을 대표하는 이차원 DFN에 대하여 절리텐서의 일차불변량을 이용하여 절리의 빈도 및 길이가 보정되었다. 현장수압시험을 수행하여 선형의 등가유로관으로 취급된 현장의 절리에 대한 유효수리간극이 추정되었다. 연구지역의 이차원 DFN 블록에 대한 REV의 크기는 블록 크기에 따른 지향적 블록수리전도도의 변화를 분석하여 20 m × 20 m으로 결정하였다. 연구지역의 이차원 DFN에 대하여 산정된 지향적 블록수리전도도는 현장의 수리적 거동과 관련하여 강한 이방성을 지시한다. 다양한 방향으로의 이차원 DFN에서 산정된 지향적 블록수리전도도를 바탕으로 연구지역의 절리성 암반에 대한 삼차원 수리전도텐서가 추정되었다. 삼차원 수리전도텐서의 주성분 방향은 현장에서 구분된 주 절리군의 방향과 일치한다.

• 대한기계학회논문집B
• /
• 제20권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.

• 한국전산유체공학회지
• /
• 제6권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.

• 대한기계학회논문집
• /
• 제18권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.

• 한국정밀공학회:학술대회논문집
• /
• 한국정밀공학회 1995년도 추계학술대회 논문집
• /
• 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.

• 대한수학회지
• /
• 제47권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.

• 소성∙가공
• /
• 제1권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
• /
• 제17권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.

• 대한기계학회논문집B
• /
• 제28권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
• /
• 제2권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.