• 대한조선학회논문집
• /
• 제33권3호
• /
• pp.94-102
• /
• 1996

유한변형문제에서 변형구배텐서를 탄소성 성분으로 분해하기 위한 가산분해와 곱분해방법에 대해서 설명하고, 이 두 방법에서 파생되는 역학량들의 의미와 그 차이점을 보였다. 변형구배에 대한 기존의 곱분해와 가산분해로 얻어지는 변형속도구배는 가산적으로 표현되지 않으며, 소성변형속도구배는 탄성변형의 영향을 받고 있다. 본 연구에서는 공축소성 가정을 도입하고, 수정된 곱분해를 통하여 소성변형속도구배가 탄성변형에 영향을 받지 않는 가산적인 변형속도 구배를 얻었다.

• Journal of Mechanical Science and Technology
• /
• 제16권11호
• /
• pp.1522-1539
• /
• 2002

In the present study, an explicit algebraic stress model is shown to be the exact tensor representation of algebraic stress model by directly solving a set of algebraic equations without resort to tensor representation theory. This repeals the constraints on the Reynolds stress, which are based on the principle of material frame indifference and positive semi-definiteness. An a priori test of the explicit algebraic stress model is carried out by using the DNS database for a fully developed channel flow at Rer = 135. It is confirmed that two-point correlation function between the velocity fluctuation and the Laplacians of the pressure-gradient i s anisotropic and asymmetric in the wall-normal direction. Thus, a novel composite algebraic Reynolds stress model is proposed and applied to the channel flow calculation, which incorporates non-local effect in the algebraic framework to predict near-wall behavior correctly.

• 한국소성가공학회:학술대회논문집
• /
• 한국소성가공학회 2003년도 춘계학술대회논문집
• /
• pp.156-159
• /
• 2003

Surface normal vector and surface velocity are very important parameters to simulate rolling processes precisely. In this study, Local displacement functions are constructed for each node on the contact surface and parameters are found by the least square fitting of displacement on the neighbor nodes. Deformation gradient tensor is calculated from the displacement function and surface normal vector and velocity also can be derived. Flat rolling simulation model is presented on the basis of the suggested contact scheme. Series of rolling process simulation are carried out and the results are compared with the experiments.

• 한국가시화정보학회:학술대회논문집
• /
• 한국가시화정보학회 2005년도 추계학술대회 논문집
• /
• pp.58-63
• /
• 2005

Three-dimensional vortical structures in the wake behind a road vehicle has been visualized with the help of two-dimensional PIV measurement data. A three-dimensional velocity field has been reconstructed from several sectional measurement data in the x-y, y-z and z-x planes. Isovorticity surface observed by stacking only the sectional data in each plane, does not show the vortical structures within the recirculation region but represents only the strong shear flows. Thus, in the present study, the velocity component normal to the x-y plane is obtained by interpolating those velocities in the z-x plane. Then, a $\lambda_{2}$-definition which captures the local pressure minimum or vortex core, is applied to visualize the vortices in the recirculation region. The final results represent a successful configuration for the three-dimensional vortices.

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

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

• 대한기계학회논문집
• /
• 제15권1호
• /
• pp.262-273
• /
• 1991

본 연구에서는 벽면에서의 압력은 압력측정 정치를 제작하여 측정하였으며 평 균유동 및 난류량들은 열선 유속계(hot wire system)을 사용하여 측정하였다.

• 유변학
• /
• 제10권4호
• /
• pp.202-216
• /
• 1998

기존의 'Orthotropic' 근사모델의 개선된 형태인 ORW를 새로운 유동 자료를 이용하여 수치적으로 구하였다. 기존의 'Orthotropic' 근사모델인 ORF나 ORL은 특히 전단유동 하에서 상호작용상수 $C_1$<0.001인 경우 비물리적 진동특성을 나타낸다. 물론 center-gated disk와 같은 비균일 유동하에서도 비물리적 진동특성을 나타내고 'Distribution Function Calculation'과 비교하여 배향 상태를 낮게 예측한다. 이런 현상들은 바로 least-square 최적화 시 사용된 유동 자료에 기인한 것을 알 수 있었다. 작은 상호작용계수의 균일 유동 자료를 이용하여 최적화를 한 ORW의 경우 비물리적 진동특성도 나타나지 않았고 균일 및 비균일 유동하에서 모두 정성적으로 잘 일치함을 확인할 수 있었다. 최적화 시 사용된 함수의 선택은 근사모델을 발전시키는데 그다지 영향을 미치지 못하였다. 하지만, 모든 배향 텐서의 eigenvalue들을 고려하면 보다 정량적으로 발전시킬 수 있지만 이들의 함수모양 선택은 중요하고 어려운 문제다. 비교를 위하여 ORW와 다른 여러 가지 근사모델을 이용하여 Film-gated strip과 Center-gated disk에 대한 연계효과 및 평면속도구배를 포함한 사출성형 충전공정의 수치모사를 수행하였다. ORW가 'Distribution Function Calculation' 과 비교하여 정량적으로도 거의 비슷한 결과를 예측함을 보여주지만 실제 실험자료와 비교하였을 때 약간의 차이가 있음을 확인하였다. 따라서 좀 더 정확히 섬유의 배향도를 예측하기 위해서는 섬유들의 상호작용을 나타내는 항의 모델링의 변화가 요구된다.

• Korea-Australia Rheology Journal
• /
• 제15권3호
• /
• pp.131-150
• /
• 2003

A new numerical algorithm of finite element methods is presented to solve high Deborah number flow problems with geometric singularities. The steady inertialess planar 4 : 1 contraction flow is chosen for its test. As a viscoelastic constitutive equation, we have applied the globally stable (dissipative and Hadamard stable) Leonov model that can also properly accommodate important nonlinear viscoelastic phenomena. The streamline upwinding method with discrete elastic-viscous stress splitting is incorporated. New interpolation functions classified as rational interpolation, an alternative formalism to enhance numerical convergence at high Deborah number, are implemented not for the whole set of finite elements but for a few elements attached to the entrance comer, where stress singularity seems to exist. The rational interpolation scheme contains one arbitrary parameter b that controls the singular behavior of the rational functions, and its value is specified to yield the best stabilization effect. The new interpolation method raises the limit of Deborah number by 2∼5 times. Therefore on average, we can obtain convergent solution up to the Deborah number of 200 for which the comer vortex size reaches 1.6 times of the half width of the upstream reservoir. Examining spatial violation of the positive definiteness of the elastic strain tensor, we conjecture that the stabilization effect results from the peculiar behavior of rational functions identified as steep gradient on one domain boundary and linear slope on the other. Whereas the rational interpolation of both elastic strain and velocity distorts solutions significantly, it is shown that the variation of solutions incurred by rational interpolation only of the elastic strain is almost negligible. It is also verified that the rational interpolation deteriorates speed of convergence with respect to mesh refinement.

• 대한기계학회논문집B
• /
• 제29권11호
• /
• pp.1265-1276
• /
• 2005

The elliptic conceptual second moment models for turbulent heat fluxes, which are proposed on the basis of elliptic-blending and elliptic-relaxation equations, are applied to calculate the combined forced and natural turbulent convection in a vertical plane channel. The models satisfy the near-wall balance between viscous diffusion, viscous dissipation and temperature-pressure gradient correlation, and also have the characteristics of approaching its respective conventional high Reynolds number model far away from the wall. Also the models are closely linked to the elliptic blending model which is used for the prediction of Reynolds stress. In order to calibrate the heat flux models, firstly, the distributions of mean temperature and scala flux in fully developed channel flow with constant wall difference temperature are solved by the present models. The buoyancy effect on the turbulent characteristics including the mean velocity and temperature, the Reynolds stress tensor, and the turbulent heat flux vector are examined. In the opposing flow, the turbulent transport is greatly enhanced with both the Reynolds stresses and the turbulent heat fluxes being remarkably increased; whereas, in the aiding flow, the opposite change is observed. The results of prediction are directly compared to the DNS to assess the performance of the model predictions and show that the behaviors of the turbulent heat transfer in the whole flow region are well captured by the present models.