• Korea-Australia Rheology Journal
• /
• v.19 no.4
• /
• pp.211-219
• /
• 2007

Here we demonstrate complex transient behavior of viscoelastic liquid described numerically with the Leonov model in straight and contraction channel flow domains. Finite element and implicit Euler time integration methods are employed for spatial discretization and time marching. In order to stabilize the computational procedure, the tensor-logarithmic formulation of the constitutive equation with SUPG and DEVSS algorithms is implemented. For completeness of numerical formulation, the so called traction boundaries are assigned for flow inlet and outlet boundaries. At the inlet, finite traction force in the flow direction with stress free condition is allocated whereas the traction free boundary is assigned at the outlet. The numerical result has illustrated severe forward-backward fluctuations of overall flow rate in inertial straight channel flow ultimately followed by steady state of forward flow. When the flow reversal occurs, the flow patterns exhibit quite complicated time variation of streamlines. In the inertialess flow, it takes much more time to reach the steady state in the contraction flow than in the straight pipe flow. Even in the inertialess case during startup contraction flow, quite distinctly altering flow patterns with the lapse of time have been observed such as appearing and vanishing of lip vortices, coexistence of multiple vortices at the contraction comer and their merging into one.

• Journal of computational fluids engineering
• /
• v.6 no.4
• /
• pp.26-34
• /
• 2001

The flow past a circular cylinder forced to vibrate transversely is numerically simulated by solving the two-dimensional Navier-Stokes equations modified by the vibration velocity of a circular cylinder at a Reynolds number of 164. The higher-order finite difference scheme is employed for the spatial discretization along with the second order Adams-Bashforth and the first order backward-Euler time integration. The calculated cylinder vibration frequency is between 0.60 and 1.30 times of the natural vortex-shedding frequency. The calculated oscillation amplitude extends to 25% of the cylinder diameter and in the case of the lock-in region it is 60%. It is made clear that the cylinder oscillation has influence on the wake pattern, the time histories of the drag and lift forces, power spectral density and phase diagrams, etc. It is found that these results include both the periodic (lock-in) and the quasi-periodic (non-lock-in) state. The vortex shedding frequency equals the driving frequency in the lock-in region but is independent in the non-lock-in region. The mean drag and the maximum lift coefficient increase with the increase of the forcing amplitude in the lock-in state. The lock-in boundaries are also established from the present direct numerical simulation.

• 한국전산유체공학회:학술대회논문집
• /
• 2001.05a
• /
• pp.181-188
• /
• 2001

The flow past a circular cylinder forced to vibrate transversely is numerically simulated by solving the two-dimensional Wavier-Stokes equations modified by the vibration velocity of a circular cylinder at a Reynolds number of 164. The higher-order finite difference scheme is employed for the spatial discretization along with the second order Adams-Bashforth and the first order backward-Euler time integration. The calculated cylinder vibration frequency is between 0.60 and 1.30 times of the natural vortex-shedding frequency. The calculated oscillation amplitude extends to $25\%$ of the cylinder diameter and in the case of the lock-in region it is $60\%$. It is made clear that the cylinder oscillation has influence on the wake pattern, the time histories of the drag and lift forces, power spectral density and phase diagrams, etc. It is found that these results include both the periodic (lock-in) and the quasi-periodic (non-lock-in) state. The vortex shedding frequency equals the driving frequency in the lock-in region but is independent in the non-lock-in region. The mean drag and the maximum lift coefficient increase with the increase of the forcing amplitude in the lock-in state. The lock-in boundaries are also established from the present direct numerical simulation.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.36 no.2
• /
• pp.115-122
• /
• 2008

The effects of characteristic condition number on the convergence of preconditioned Euler equations were investigated. The two-dimensional preconditioned Euler equations adopting Choi and Merkle's preconditioning and the temperature preconditioning are considered. Preconditioned Roe's FDS scheme was adopted for spatial discretization and preconditioned LU-SGS scheme was used for time integration. It is shown that the convergence characteristics of the Euler equations are strongly affected by the characteristic condition number, and there is an optimal characteristic condition number for a problem. The optimal characteristic condition numbers for the Choi and Merkle's preconditioning and temperature preconditioning are different.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.36 no.2
• /
• pp.123-130
• /
• 2008

The effects of characteristic condition number on the convergence of preconditioned Navier-Stokes equations were investigated. The two-dimensional preconditioned Navier-Stokes adopting Choi and Merkle's preconditioning and the temperature preconditioning are considered. Preconditioned Roe's FDS scheme was adopted for spatial discretization and preconditioned LU-SGS scheme was used for time integration. It is shown that the convergence characteristics of the Navier-Stokes equations are strongly affected by the characteristic condition number. Also it is shown that the optimal characteristic condition numbers for viscous flows are larger than that in inviscid flows.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.30 no.8
• /
• pp.37-45
• /
• 2002

A two-dimensional unsteady, inviscid flow solver has been developed for the simulation of airfoil-vortex interactions on unstructured dynamically adapted meshes. The Euler solver is based on a second-order accurate implicit time integration using a point Gauss-Seidel relaxation scheme and a dual time-step subiteration. A vertex-centered, finite-volume discretization is used in conjunction with the Roe's flux-difference splitting. An unsteady solution-adaptive dynamic mesh scheme is used by adding and deleting mesh points to take account of both spatial and temporal variations of the flow field. The effect of vortex interaction on the aeroelastic displacement of an airfoil attached to the idealized two degree-of-freedom spring system is investigated.

• Wind and Structures
• /
• v.26 no.1
• /
• pp.35-43
• /
• 2018

This paper describes the effect of different testing parameters (configuration of infrastructure and truck position on road) on truck aerodynamic coefficients under cross wind conditions, by means of a numerical approach known as Large Eddy Simulation (LES). In order to estimate the air flow behaviour around both the infrastructure and the truck, the filtered continuity and momentum equations along with the Smagorinsky-Lilly model were solved. A solution for these non-linear equations was approached through the finite volume method (FVM) and using temporal and spatial discretization schemes. As for the results, the aerodynamic coefficients acting on the truck model exhibited nearly constant values regardless of the Reynolds number. The flat ground is the infrastructure where the rollover coefficient acting on the truck model showed lowest values under cross wind conditions (yaw angle of $90^{\circ}$), while the worst infrastructure studied for vehicle stability was an embankment with downward-slope on the leeward side. The position of the truck on the road and the value of embankment slope angle that minimizes the rollover coefficient were determined by successfully applying the Response Surface Methodology.

• Structural Engineering and Mechanics
• /
• v.34 no.6
• /
• pp.685-702
• /
• 2010

In this paper, a nonlinear finite element procedure is presented for the dynamic analysis of reinforced concrete shell structures. A computer program, named RCAHEST (Reinforced Concrete Analysis in Higher Evaluation System Technology), was used. A 4-node flat shell element with drilling rotational stiffness was used for spatial discretization. The layered approach was used to discretize the behavior of concrete and reinforcement in the thickness direction. Material nonlinearity was taken into account by using tensile, compressive and shear models of cracked concrete and a model of reinforcing steel. The smeared crack approach was incorporated. The low-cycle fatigue of both concrete and reinforcing bars was also considered to predict a reliable dynamic behavior. The solution to the dynamic response of reinforced concrete shell structures was obtained by numerical integration of the nonlinear equations of motion using Hilber-Hughes-Taylor (HHT) algorithm. The proposed numerical method for the nonlinear dynamic analysis of reinforced concrete shell structures was verified by comparison of its results with reliable experimental and analytical results.

• 한국전산유체공학회:학술대회논문집
• /
• 1998.05a
• /
• pp.145-155
• /
• 1998

A computer code has been developed for the computation of the viscous flow around a ship model with the free surface. In this code, the incompressible Reynolds-averaged Navier-Stokes equations are solved numerically by a finite difference method which employes second-order finite differences for the spatial discretization and a four-stage Runge-Kutta scheme for the temporal integration of the governing equations. For the turbulence closure, a modified version of the Baldwin-Lomax model is exploited. The location of the free surface is determined by solving the equation of the kinematic free-surface condition using the Lax-Wendroff scheme and the boundary-fitted grid is generated at each time step so that one of the grid surfaces always coincides with the free surface. An inviscid approximation of the dynamic free-surface boundary condition is applied as the boundary conditions for the velocity and pressure on the free surface. To validate the computational method and the computer code developed in the present study, the numerical computations are carried out for both Wigley parabolic hull and Series 60 $C_B=0.6$ ship model and the computational results are compared with the experimental data.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.25 no.10
• /
• pp.1384-1391
• /
• 2001

Three-dimensional flow analysis is implemented to investigate the flow through transonic axial-flow compressor rotor(NASA R67) and to evaluate the performances of Abid's low-Reynolds-number k-$\omega$ and Baldwin-Lomax turbulence models. A finite volume method is used fur spatial discretization. The equations are solved implicitly in time by the use of approximate factorization. The upwind difference scheme is used for inviscid terms and viscous terms are approximated with central difference. The flux-difference-splitting method of Roe is used to obtain fluxes at the cell faces. Numerical analysis is performed near peak efficiency and near stall. The results are compared with the experimental data for NASA R67 rotor. Blade-to-Blade Mach number distributions are compared to confirm the accuracy of the code. From the results, it is concluded that Abid'k-$\omega$ model is better for the calculation of flow rate and efficiency than Baldwin-Lomax model. But, the predictions for Mach number and shock structure are almost the same.