• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.29 no.1 s.232
• /
• pp.145-155
• /
• 2005

In this study, a finite element method program code which can be accomodate boundary conditions on the complex surfaces has been developed to simulate polymer extruder processes. The analysis method includes the fractional 4-step method for efficient computation time and compact usage of memory storage to solve the velocities and the pressure values from the Navier-Stokes equation. By using the developed program which was verified with simple Poiseuille flow mixture phenomena in single-and twin-screw extruder are analyzed. It is concluded that the proposed method resulte Poiseuille Poiseuille d in fair agreement with the exact solution of simple flow and the back flow near the entrance happens in single-screw model. It is identified that the location and values of maximum pressure in the twin screw extruder model. It is expected that the Velocity field found can be used to predict the degree of mixture in the extruder barrel.

• Structural Engineering and Mechanics
• /
• v.91 no.6
• /
• pp.627-642
• /
• 2024

In the present work, the deformation and stresses induced in a functionally graded disk have been reported for different loading conditions. The governing differential equation is solved using the classical method namely Navier's method by considering thermal and mechanical boundary conditions at the surface of the disk. To simplify solving the second-order differential equation, a plane stress condition was assumed. Following validation using a one-dimensional steady-state heat condition problem, temperature variations were computed for constant heat generation and varying conductivity. The research aims to investigate both the individual and combined effects of rotation, gravity, and temperature with constant heat generation on a hollow disk operating under complex loading conditions. The results demonstrated a high degree of accuracy when compared with those in existing literature. Material properties, such as Young's modulus, density, conductivity, and thermal expansion coefficient, were modeled using a power law variation along the disk's radius by considering aluminum as a base material. The proposed analytical method is straightforward, providing valuable insights into the behavior of disks under various loading conditions. This method is particularly useful for researchers and industries in selecting appropriate loading conditions and grading parameters for engineering applications, including aerospace components, energy systems, and rotary machinery parts.

• Journal of Mechanical Science and Technology
• /
• v.20 no.1
• /
• pp.167-175
• /
• 2006

A parametric study has been accomplished to figure out the effects of elliptic cylinder thickness, angle of attack, and Reynolds number on the unsteady lift and drag forces exerted on the elliptic cylinder. A two-dimensional incompressible Navier-Stokes flow solver is developed based on the SIMPLER method in the body-intrinsic coordinates system to analyze the unsteady viscous flow over elliptic cylinder. Thickness-to-chord ratios of 0.2, 0.4, and 0.6 elliptic cylinders are simulated at different Reynolds numbers of 400 and 600, and angles of attack of $10^{\circ},\;20^{\circ},\;and\;30^{\circ}$. Through this study, it is observed that the elliptic cylinder thickness, angle of attack, and Reynolds number are very important parameters to decide the lift and drag forces. All these parameters also affect significantly the frequencies of the unsteady force oscillations.

• Journal of Mechanical Science and Technology
• /
• v.17 no.3
• /
• pp.457-466
• /
• 2003

In this paper, the numerical finite element formulations were derived for the linearized Navier-Stokes' equations with assumptions of two-dimensional incompressible, homogeneous viscous fluid field, and small oscillation and the FAMD (Fluid Added Mass and Damping) code was developed for practical applications calculating the fluid added mass and damping. In formulations, a fluid domain is discretized with C$\^$0/-type quadratic quadrilateral elements containing eight nodes using a mixed interpolation method, i.e., the interpolation function for the velocity variable is approximated by a quadratic function based on all eight nodal points and the interpolation function for the pressure variable is approximated by a linear function based on the four nodal points at vertices. Using the developed code, the various characteristics of the fluid added mass and damping are investigated for the concentric cylindrical shell and the actual hexagon arrays of the liquid metal reactor cores.

• Proceedings of the KSME Conference
• /
• 2007.05b
• /
• pp.2666-2671
• /
• 2007

This work presents microchannel heat sink shape optimization procedure using Kriging method. Design variables relating to microchannel width, depth and fin width are selected, and thermal resistance has been taken as objective function. Design points are selected through a three-level fractional factorial design of sampling method. Navier-Stokes and energy equations for laminar flow and conjugate heat transfer are solved at these design points using a finite volume solver. Solutions are carefully validated with experimental results. Using the numerically evaluated objective function, a surrogate model (Kriging) is constructed and optimum point is searched by sequential quadratic programming. The process of shape optimization greatly improves the thermal performance of microchannel heat sink under constant pumping power.

• 한국전산유체공학회:학술대회논문집
• /
• 2009.11a
• /
• pp.177-184
• /
• 2009

Cavitating flow simulation is of practical importance for many engineering systems, such as marine propellers, pump impellers, nozzles, injectors, torpedoes, etc. The present work has focused on the simulation of cavitating flow past cylinders with strong side flows. The governing equation is the Navier-Stokes equation based on the homogeneous mixture model. The momentum and energy equation is in the mixture phase while the continuity equation is solved liquid and vapor phase, separately. An implicit dual time and preconditioning method are employed for computational analysis. For the code validation, the results from the present solver have been compared with experiments and other numerical results. A fairly good agreement with the experimental data and other numerical results have been obtained. After the code validation, the strong side flow was applied to include the wake flow effects of the submarine or ocean tide.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.28 no.8 s.227
• /
• pp.929-936
• /
• 2004

In the present work, shape of the mixing vane in Plus7 fuel assembly has been optimized numerically using three-dimensional Reynolds-averaged Navier-Stokes analysis of flow and heat transfer. Standard $k-{\epsilon}$ model is used as a turbulence closure. The Response surface method is employed as an optimization technique. The objective function is defined as a combination of heat transfer rate and inverse of friction loss. Bend angle and base length of mixing vane are selected as design variables. Thermal-hydraulic performances for different shapes of mixing vane have been discussed, and optimum shape has been obtained as a function of weighting factor in the objective function.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.26 no.8
• /
• pp.1088-1094
• /
• 2002

A channel flow with a high Reynolds number but coarse grids is numerically studied to investigate the prediction possibility of its turbulence which is three-dimensional and time-dependent. In the present paper, a Reynolds-Averaged Navier-Stokes (RANS) model, a Large Eddy Simulation (LES) and a Navier-Stokes equation with no model are tested with a new approach of hybrid RANS/LES, which reduces to RANS model in the boundary layers and at separation, and to Smagorinsky-like LES downstream of separation, and then compared with each other. It is found that the simulations of hybrid RANS/LES method sustain turbulence like those of LES and with no model, and the results are stable and fairly accurate. This indicates strongly that gradual improvements could lead to a simple, stable, and accurate approach to predict turbulence phenomena of wall-bounded flow.

• Journal of Mechanical Science and Technology
• /
• v.14 no.9
• /
• pp.1005-1012
• /
• 2000

Numerical optimization techniques combined with a three-dimensional thin-layer Navier-Stokes solver are presented to find an optimum shape of a stator blade in an axial compressor through calculations of single stage rotor-stator flow. Governing differential equations are discretized using an explicit finite difference method and solved by a multi-stage Runge-Kutta scheme. Baldwin-Lomax model is chosen to describe turbulence. A spatially-varying time-step and an implicit residual smoothing are used to accelerate convergence. A steady mixing approach is used to pass information between stator and rotor blades. For numerical optimization, searching direction is found by the steepest decent and conjugate direction methods, and the golden section method is used to determine optimum moving distance along the searching direction. The object of present optimization is to maximize efficiency. An optimum stacking line is found to design a custom-tailored 3-dimensional blade for maximum efficiency with the other parameters fixed.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.3 no.2
• /
• pp.161-171
• /
• 1999

Aerodynamic sensitivity analysis codes are developed via the hand-differentiation using a direct differentiation method and an adjoint method respectively from discrete two-dimensional compressible Navier-Stokes equations. Unlike previous other researches, Baldwin-Lomax algebraic turbulence model is also differentiated by hand to obtain design sensitivities with respect to design variables of interest in turbulent flows. Discrete direct sensitivity equations and adjoint equations are efficiently solved by the same time integration scheme adopted in the flow solver routine. The required memory for the adjoint sensitivity code is greatly reduced at the cost of the computational time by allowing the large banded flux jacobian matrix unassembled. Direct sensitivity code results are found to be exactly coincident with sensitivity derivatives obtained by the finite difference. Adjoint code results of a turbulent flow case show slight deviations from the exact results due to the limitation of the algebraic turbulence model in implementing the adjoint formulation. However, current adjoint sensitivity code yields much more accurate sensitivity derivatives than the adjoint code with the turbulence eddy viscosity being kept constant, which is a usual assumption for the prior researches.