• Wind and Structures
• /
• 제30권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.

• 한국항해항만학회지
• /
• 제27권5호
• /
• pp.473-478
• /
• 2003

평행한 벽 사이에 원형배관(circular pipeline)을 놓고 그 주위의 유동특성에 대한 수치연구를 수행하였다. 비압축성 유체를 가지고, Navier-Stokes 방정식을 풀었고 3차 풍상(upwind) 차분의 수치해법을 이용하였다. 한쪽 벽과의 거리가 매우 작아질 때, 볼텍스 떨어짐이 상당히 억압되는데 이것은 벽 경계와의 상호 박리 작용 때문으로 간수된다. 본 연구는 레이놀드 수의 변함과 물체가 벽에 접근함에 따른 볼텍스 떨어짐의 특성을 규명하는데 있다. 원형배관 후류와 평행벽내 유기된 박리의 상호작용을 집중적으로 다루며 서로 다른 조건에서 박리와 와역한(vorticity dynamics)의 특성을 해석하였다.

• 한국해양공학회지
• /
• 제16권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.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2002년도 학술대회지
• /
• 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.

• 동력기계공학회지
• /
• 제21권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.

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

• 동력기계공학회지
• /
• 제18권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
• /
• 제17권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
• /
• 제4권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
• /
• 제6권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.