• Journal of the Society of Naval Architects of Korea
• /
• v.33 no.3
• /
• pp.94-102
• /
• 1996

An interpretation for the additive and multiplicative decomposition theory of the deformation gradient tensor in finite deformation problems is presented. the conventional methods have not provided the additive deformation velocity gradient. Moreover the plastic deformation velocity gradients are not free from elastic deformations. In this paper, a modified multiplicative decomposition is introduced with the assumption of coaxial plastic deformation velocity gradient. This strategy well gives the additive deformation velocity gradient in which the plastic deformation velocity gradient is not affect4d by the elastic deformation.

• Journal of Mechanical Science and Technology
• /
• v.16 no.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.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 2003.05a
• /
• 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.12a
• /
• 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.

• 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.

• 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.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.15 no.1
• /
• pp.262-273
• /
• 1991

In the present study, the steady, incompressible, isothermal, developing flow in a 90.deg. rectangular cross sectional strongly curved duct with aspect ratio 1:1.5 and Reynolds number of 9.4*10$^{4}$ has been investigated. Measurements of components of mean velocities, pressures, and corresponding components of the Reynolds stress tensor are obtained with a hot-wire anemometer and pitot tube. In general, flow in a curved duct is characterized by the secondary vortices which are driven mainly by centrifugal force-radial pressure gradient imbalance, and the stress field stabilizing effects near the convex wall and destablizing effects close to the concave wall. It was found that the secondary mean velocities attain values up to 39% of the bulk velocity and are largely responsible for the convections of Reynolds stress in the cross stream plane. Therefor upstream of the bend the Reynolds stress are low. Corresponding to the small boundary layer thickness. At successive planes, large values of Reynolds stress were observed near the concave surface and the side wall.

• The Korean Journal of Rheology
• /
• v.10 no.4
• /
• pp.202-216
• /
• 1998

Improved version of previous 'Orthotropic' closure approximation, termed 'ORW' has been numerically developed using new homogeneous flow data. Previous 'Orthotropic' closure approximation, i.e., ORF or ORL showed non-physical oscillation for interaction coefficient $C_1$<0.001 at simple shear flow. It also shows non-physcial oscillation and under-prediction compared with 'Distribution Function Calculation' at non-homogeneous flow of center-gated disk. These phenomena are mainly due to the flow data of 'Distribution Function Calculation' which were used for least-square optimization. ORW obtained by fitting flow data of low interaction coefficient does not show non-physical oscillation and results in reasonably good behaviors at non-homogeneous flows as well as homogeneous flows. Fitting function forms have not been found to improve overall behaviors. It has been found that considering all the eigenvalues of orientation tensor (including the third eigenvalues) might end up with a better closure approximation than just considering the first and second eigenvalues. It is, however, very important and yet difficult to select appropriate function forms of eigenvalues. Numerical simulation including coupling and in-plane velocity gradient effects were performed for injection mold filing process with a film-gated strip and a center-gated disk using ORW and various other closure approximations for comparisons. Although ORW is in excellent agreement with 'Distribution Function Calculation', the predicted results seem to have consistent error in comparison with experimental data. The diffusivity term with constant interaction coefficient might have to be further investigated in order to accurately describe orientation states.

• Korea-Australia Rheology Journal
• /
• v.15 no.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.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.29 no.11 s.242
• /
• 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.