• Wind and Structures
• /
• v.30 no.1
• /
• pp.69-83
• /
• 2020

Modelling incompressible, neutrally stratified, barotropic, horizontally homogeneous and steady-state atmospheric boundary layer (ABL) is an important aspect in computational wind engineering (CWE) applications. The ABL flow can be viewed as a balance of the horizontal pressure gradient force, the Coriolis force and the turbulent stress divergence. While much research has focused on the increase of the wind velocity with height, the Ekman layer effects, entailing veering - the change of the wind velocity direction with height, are far less concerned in wind engineering. In this paper, a modified k-ε model is introduced for the ABL simulation considering wind veering. The self-sustainable method is discussed in detail including the precursor simulation, main simulation and near-ground physical quantities adjustment. Comparisons are presented among the simulation results, field measurement values and the wind profiles used in the conventional wind tunnel test. The studies show that the modified k-ε model simulation results are consistent with field measurement values. The self-sustainable method is effective to maintain the ABL physical quantities in an empty domain. The wind profiles used in the conventional wind tunnel test have deficiencies in the prediction of upper-level winds. The studies in this paper support future practical super high-rise buildings design in CWE.

• Journal of Navigation and Port Research
• /
• v.27 no.5
• /
• pp.473-478
• /
• 2003

Numerical study was made on the flow characteristics around a circular pipeline between parallel walls. The incompressible Navier-Stokes equations were solved by using a third-order upwind differential scheme. When the distance near a wall is small enough, the vortex shedding is almost completely suppressed because of the interaction with the wall boundary layer separation. This study aims to clarify the characteristics of the vortex shedding regime as the body approaches a wall as Reynolds number varies. The feature of separated vorticity dynamics is analyzed at different conditions with particular attention to the interaction between the pipeline wake and the induced separation on the plane walls.

• Journal of Ocean Engineering and Technology
• /
• v.16 no.5
• /
• pp.15-20
• /
• 2002

In this research, the finite difference lattice Boltzmann method(FDLBM) is used to analyze gravity currents in the lock exchange configuration that occur in many natural and man-made situations. At a lock those are seen when a gate is suddenly opened, and, in the atmosphere, when the thunderstorm outflows make a cold front. At estuaries in the ocean, the phenomenon is found between fresh water from a river and salt water in the sea. Since such interesting phenomena were recognized, pioneers have challenged to make them clear by conducing both experiments and analysis. Most of them were about the currents of liquid or Boussinesq fluids, which are assumed as incompressible. Otherwise, the difference in density of two fluids is small. The finite difference lattice Boltzmann method has been a powerful tool to simulate the flow of compressible fluids. Also, numerical predictions using FDLBM to clarify the gravity currents of compressible fluids exhibit all features, but typically observed in experimental flows near the gravity current head, including the lobe-and-cleft structure at the leading edge.

• Proceedings of the KSME Conference
• /
• 2002.08a
• /
• pp.705-706
• /
• 2002

Hemodynamics of the TPLS(Twin Pulse Life Support System) is numerically investigated to delineate the possibility of hemolysis in blood. Computational method employing finite element algorithm is utilized to solve the blood flow of the sac squeezed by moving actuator. We assume that the blood flow interacts with the sac material which is activated by the rigid body motion of the actuator. Valve dynamics at the ends of the sac is simplified as on/off type motion. We compute the transient viscous flow in the two-dimensional geometry of the blood sac. Incompressible laminar flow is simulated on the assumption of Newtonian fluid. Blood velocity has a step gradient near the throat of the sac formed by the moving actuator. According to the decrease of the gap size of blood passage, the magnitude of shear stress in the blood is dramatically increased. Numerical solutions show that the maximum value of shear stress in the blood flow in TPLS is relatively smaller than that of the roller type ECLS.

• Journal of Power System Engineering
• /
• v.21 no.2
• /
• pp.48-56
• /
• 2017

The steady-state, incompressible and three-dimensional numerical analysis was carried out to evaluate the velocity fields around the seabed tiller used for the improvement of the seabed soil and the pulling force and buoyancy generated by driving the seabed tiller. The turbulence model used in this study is a realizable $k-{\varepsilon}$ well known to be excellent for predicting the performance of the flow separation and recirculation flow as well as the boundary layer with rotation and strong back pressure gradient. As a results, a typical vortex pair appears near the adjacent rotor vane tip. When the current is stopped, there is no force when pulling the seabed tiller, but when the current flows at 1.2 knots, the force acts on the downstream side and the pulling force is much greater. In stationary currents, the buoyancy of the seabed tiller acts more strongly towards the seabed as the number of rotations of the rotor increases, but acts more strongly toward the sea surface at 1.2 knots of current.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.25 no.2
• /
• pp.156-163
• /
• 2001

The experiment is conducted on the rapidly rotating incompressible flow within a confined cylinder using LDV(Laser Doppler Velocimetry). The configurations of interest are the flows between a rotating upper disk with a bar and a stationary lower disk enclosed within a cylinder. The flow is considered to be an axisymmetric undisturbed basic flow. The results show that the flow is strongly dependent on the radius and the shape of bar but is negligibly affected by the Reynolds number in turbulent flow. It is observed that in the lid-driven case the main forms near the wall as the Reynolds number increases. The thin bar causes the second axial flow due to the suction effect and the thick bar causes the main flow to be pulled toward the surface of the bar. The step bar shows the dual effect of the two. 1:2 tilt bar shows that the main flow distributes wider than the other cases in which interference occurs due step bar.

• Journal of Power System Engineering
• /
• v.18 no.3
• /
• pp.59-65
• /
• 2014

The computational fluid dynamics was carried out to evaluate turbulent models on the swirling flow of a gun-type gas burner(GTGB) according to the mesh size. The commercial SC/Tetra software was used for a steady-state, incompressible and three-dimensional numerical analysis. In consequence, the velocity magnitude from the exit of a GTGB and the flowrate predicted by the turbulent models of MP k-${\varepsilon}$, Realizable k-${\varepsilon}$ and RNG k-${\varepsilon}$ agree with the results measured by an experiment very well. Moreover, the turbulent kinetic energy predicted by the turbulent model of standard k-${\varepsilon}$ with mesh type C only agrees with the experimental result very well along the radial distance. On the other hand, the detailed prediction of the information of swirling flow field near the exit of a GTGB at least needs a CFD analysis using a fairly large-sized mesh such as a mesh type C.

• Korea-Australia Rheology Journal
• /
• v.17 no.3
• /
• pp.99-110
• /
• 2005

This work presents results of finite element analysis of isothermal incompressible creeping viscoelastic flows with the tensor-logarithmic formulation of the Leonov model especially for the planar geometry with singular comers in the domain. In the case of 4:1 contraction flow, for all 5 meshes we have obtained solutions over the Deborah number of 100, even though there exists slight decrease of convergence limit as the mesh becomes finer. From this analysis, singular behavior of the comer vortex has been clearly seen and proper interpolation of variables in terms of the logarithmic transformation is demonstrated. Solutions of 4:1:4 contraction/expansion flow are also presented, where there exists 2 singular comers. 5 different types spatial resolutions are also employed, in which convergent solutions are obtained over the Deborah number of 10. Although the convergence limit is rather low in comparison with the result of the contraction flow, the results presented herein seem to be the only numerical outcome available for this flow type. As the flow rate increases, the upstream vortex increases, but the downstream vortex decreases in their size. In addition, peculiar deflection of the streamlines near the exit comer has been found. When the spatial resolution is fine enough and the Deborah number is high, small lip vortex just before the exit comer has been observed. It seems to occur due to abrupt expansion of the elastic liquid through the constriction exit that accompanies sudden relaxation of elastic deformation.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.4 no.4
• /
• pp.477-487
• /
• 2012

We investigated the aerodynamic characteristics of a three-dimensional (3D) wing with an endplate in the vicinity of the free surface by solving incompressible Navier-Stokes equations with the turbulence closure model. The endplate causes a blockage effect on the flow, and an additional viscous effect especially near the endplate. These combined effects of the endplate significantly reduce the magnitudes of the velocities under the lower surface of the wing, thereby enhancing aerodynamic performance in terms of the force coefficients. The maximum lift-to-drag ratio of a wing with an endplate is increased 46% compared to that of wing without an endplate at the lowest clearance. The tip vortex of a wing-with-endplate (WWE) moved laterally to a greater extent than that of a wing-without-endplate (WOE). This causes a decrease in the induced drag, resulting in a reduction in the total drag.

• Interaction and multiscale mechanics
• /
• v.6 no.2
• /
• pp.237-255
• /
• 2013

In this paper, a meshfree-enriched finite element method (ME-FEM) is introduced for the large deformation analysis of nonlinear path-dependent problems involving contact. In linear ME-FEM, the element formulation is established by introducing a meshfree convex approximation into the linear triangular element in 2D and linear tetrahedron element in 3D along with an enriched meshfree node. In nonlinear formulation, the area-weighted smoothing scheme for deformation gradient is then developed in conjunction with the meshfree-enriched element interpolation functions to yield a discrete divergence-free property at the integration points, which is essential to enhance the stress calculation in the stage of plastic deformation. A modified variational formulation using the smoothed deformation gradient is developed for path-dependent material analysis. In the industrial metal forming problems, the mortar contact algorithm is implemented in the explicit formulation. Since the meshfree-enriched element shape functions are constructed using the meshfree convex approximation, they pose the desired Kronecker-delta property at the element edge thus requires no special treatments in the enforcement of essential boundary condition as well as the contact conditions. As a result, this approach can be easily incorporated into a conventional displacement-based finite element code. Two elasto-plastic problems are studied and the numerical results indicated that ME-FEM is capable of delivering a volumetric locking-free and pressure oscillation-free solutions for the large deformation problems in metal forming analysis.