• Wind and Structures
• /
• v.16 no.5
• /
• pp.411-431
• /
• 2013

A computational study of vortex-induced transverse vibrations of a cylinder with low mass-damping is presented. An Arbitrary Lagrangian-Eulerian (ALE) formulation of the Unsteady Reynolds-Averaged Navier-Stokes equations (URANS), along with the Spalart-Allmaras (SA) one-equation turbulence model, are coupled conservatively with rigid body motion equations of the cylinder mounted on elastic supports in order to study the amplitude and frequency response of a freely vibrating cylinder, its flow-induced motion, Vortex Street, near-wake flow structure, and unsteady loading in a moderate range of Reynolds numbers. The time accurate response of the cylinder from rest to its limit cycle is studied to explore the effects of Reynolds number on the start of large displacements, motion amplitude, and frequency. The computational results are compared with published physical experiments and numerical studies. The maximum amplitudes of displacements computed for various Reynolds numbers are smaller than the experimental values; however, the overall agreement of the results is quite satisfactory, and the upper branch of the limit-cycle displacement amplitude vs. reduced velocity response is captured, a feature that was missed by other studies. Vortex shedding modes, lock-in phenomena, frequency response, and phase angles are also in agreement with experiments.

• Wind and Structures
• /
• v.39 no.3
• /
• pp.207-222
• /
• 2024

This paper presents a numerical simulation method for snow drift that takes into account the cohesion effect of snow particles. The critical state of free collapse accumulation of idealized snow particles is used to indirectly infer the effect of interparticle interactions on snow transport and re-accumulation. With the help of the Hertz-Mindlin with JKR cohesion contact model, the particle angle of repose is calibrated with a number of contact parameters through numerical experiment. The surface energy for a given property of snow particles is determined using the observed snow angle of repose, and a continuous-discrete snow drift two-way coupled numerical model incorporating these optimized contact parameters is developed. The snow redistribution pattern on a stepped flat roof structure is simulated, and the results are found to be consistent with those of the field measured in terms of phenomena and general laws, verifying the achievability and effectiveness of the presented method. To eliminate the influence of environmental conditions, wind tunnel tests are also conducted, and it is found that the reconstructed depth and reaccumulated angle of snowdrift resulting from the numerical simulation are in closer agreement with the experimental results, further confirming the enhancement achieved by introducing the contact effect.

• Clean Technology
• /
• v.25 no.3
• /
• pp.206-222
• /
• 2019

Recently, various multiphase flows have been developed, and among them some models have been commercialized. However, most of them have been developed based on a pressure-based approach; therefore, various numerical difficulties were involved inherently. Accordingly, in order to overcome these numerical difficulties, a multiphase flow model, MultiPhaSe flow (MPS), following a fractional-flow based approach was developed. In this study, by combining a contaminant transport module describing an enhanced dissolution effect of a surfactant with MPS, a MultiPhaSe flow and TranSport (MPSTS) model was developed. The developed model was verified using the analytical solution of Clement. The MPSTS model can simulate the process of surfactant enhanced aquifer remediation including interphase mass transfer and contaminant transport in multiphase flow by using the coupled particle tracking method and Lagrangian-Eulerian method. In this study, a surfactant was used in a non aqueous phase liquid (NAPL) contaminated area, and the effect of hydro-geological heterogeneity in the layered media on remediation efficiency was studied using the developed model. According to the numerical simulation, when hydraulic conductivity in a lower layer is 10 times, 20 times, and 50 times larger than that in an upper layer, the concentration of dissolved diesel in the lower layer is much higher than that in the upper layer because the surfactant moves faster along the lower layer owing to preferential flow; thus, the surfactant enhances dissolution of residual non aqueous phase liquid in the lower layer.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.40 no.6
• /
• pp.347-355
• /
• 2016

This study presents a Lattice Boltzmann Method (LBM) coupled with a momentum-exchange approach/fictitious domain (MEA/FD) method for the simulation of particle suspensions. The method combines the advantages of the LB and the FD methods by using two unrelated meshes, namely, a Eulerian mesh for the flow domain and a Lagrangian mesh for the solid domain. The rigid body conditions are enforced by the momentum-exchange scheme in which the desired value of velocity is imposed directly in the particle inner domain by introducing a pseudo body force to satisfy the constraint of rigid body motion, which is the key idea of a fictitious domain (FD) method. The LB-MEA/FD method has been validated by simulating two different cases, and the results have been compared with those through other methods. The numerical evidence illustrated the capability and robustness of the present method for simulating particle suspensions.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.26 no.12
• /
• pp.1738-1746
• /
• 2002

In the FSI (Fluid-Structure Interaction) problems, two different governing equations are to be solved together. One is fur the fluid and the other for the structure. Furthermore, a kinematic constraint should be imposed along the boundary between the fluid and the structure. We use the combined formulation, which incorporates both the fluid and structure equations of motion into a single coupled variational equation so that it is not necessary to calculate the fluid force on the surface of structure explicitly when solving the equations of motion of the structure. A two-dimensional channel flow divided by a Bernoulli-Euler beam is considered and the dynamic response of the beam under the influence of channel flow is studied. The Navier-Stokes equations are solved using a P2P1 Galerkin finite element method with ALE (Arbitrary Lagrangian-Eulerian) algorithm. The internal structural damping effect is not considered in this study and numerical results are compared with a previous work fer steady case. In addition to the Reynolds number, two non-dimensional parameters, which govern this fluid-structure system, are proposed. It is found that the larger the dynamic viscosity and density of the fluid are, the larger the damping of the beam is. Also, the added mass is found to be linearly proportional to the density of the fluid.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.44 no.4
• /
• pp.290-297
• /
• 2016

A two-phase Large Eddy Simulation(LES) has been conducted to investigate breakup and atomization of a liquid jet in a cross turbulent flow in a rectangular duct. Gas-droplet two-phase flow was solved by a coupled Eulerian-Lagrangian method which tracks every individual particles. Effects of liquid breakup models, sub-grid scale models, and a order of spatial discretization was investigated. The penetration depth in cross flow was comparable with experimental data by varying breakup model and LES scheme. SMD(Sauter Mean Diameter) distribution downstream of jet was analyzed.

• Journal of Ocean Engineering and Technology
• /
• v.32 no.5
• /
• pp.324-332
• /
• 2018

It is known that damages to the subsea cables used for electric power transmission between islands and countries, including renewable energy from offshore wind power, current, tides, etc., cost much to restore, which causes social and economic losses. Various types of fishing rigs and anchors have been reported to be the greatest hazards to subsea cables. It is possible to design and construct a suitable protection facility for a subsea cable by precisely estimating the underwater behavior of such hazardous apparatuses. In this study, numerical simulations of the underwater behaviors of various hazardous apparatuses were carried out using fluid-structure interaction (FSI) analysis as a basic study to simulate the actual behavior phenomena of hazardous apparatuses in relation to a subsea cable. In addition, the underwater drop characteristics according to the types of hazardous apparatuses were compared. In order to verify the accuracy of the FSI analysis method used in this study, we compared the test results for underwater drops of a steel ball bearing. Stock anchors, stockless anchors, and rocket piles, which were actually reported to be the cases of damage to subsea cables along the southwest coast of Korea, were considered as the hazardous apparatuses for the numerical simulations. Each hazardous apparatus was generated by a Lagrangian model and coupled with the fluid domain idealized by the Eulerian equation to construct the three-dimensional FSI analysis model. The accuracy of the numerical simulation results was verified by comparing them with the analytical solutions, and the underwater drop characteristics according to the types of hazard apparatuses were compared.

• Geomechanics and Engineering
• /
• v.34 no.2
• /
• pp.209-220
• /
• 2023

In this study, a series of three-dimensional numerical analyses were carried out to investigate the penetration performance of a dynamically installed Hall anchor. The advanced coupled Eulerian-Lagrangian (CEL) technique was adopted to accurately simulate the large soil deformation during the vertical penetration of a Hall anchor. In total, 52 numerical analyses were conducted to investigate the relationship between anchor penetration depth and the initial kinematic energy. Moreover, a sensitivity analysis was performed to investigate the effects of soil shear strength and soil type on the penetration mechanism of a drop anchor under self-weight. There is a monotonic increase in the penetration depth with an increasing anchor weight when the topsoil of the riverbed is not subjected to erosion. On the other hand, all the computed depths significantly increase when soil erosion is taken into consideration. This is mainly due to an enhanced initial kinematic energy from an increased dropping depth. Both depths increase exponentially with the initial kinematic energy. An enhanced shear strength can potentially increase the side resistance and end-bearing pressure around a drop anchor, thus significantly reducing the downward penetration of a hall anchor. Design charts are developed to directly estimate penetration depth and associated plastic zone due to dynamically installed anchor at arbitrary soil shear strength and anchor kinematic energy.

• Journal of the Korea institute for structural maintenance and inspection
• /
• v.24 no.6
• /
• pp.17-24
• /
• 2020

In this study, series of nonlinear seismic analysis were performed on a reinforced concrete intake tower surrounded by water. To consider the fluid effect around the structure, analysis models were composed using an added mass and CEL approach. At this time, the implicit method was used for the added mass model, and the explicit method was used for the fluid structure interaction model. The input motions were scaled to correspond to 500, 1000, and 2400 years return period of the same artificial earthquake. To estimate the counteractivity of the fluid coupled model, models without fluid effect were constructed and used as a reference. The material models of concrete and reinforcement were selected to consider the nonlinear behavior after yielding, and analysis were performed by ABAQUS. As results, in the acceleration response spectrum of the structure, it was found that the influence of the surrounding fluid reducing the peak frequency and magnitude corresponding to the fundamental frequency of the structure. However, the added mass model did not affect the peak value corresponding to the higher mode. The sectional moments were increased significantly in the case of the added mass model than those of the reference model. Especially, this amplification occurred largely for a small-sized earthquake response in which linear behavior is dominant. In the fluid structure interaction model, the sectional moment with a low frequency component amplifies compared to that of the reference model, but the sectional moment with a high requency component was not amplified. Based in these results, it was evaluated that the counteractivity of the additive mass model was greater than that of the fluid structure interaction model.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.40 no.3
• /
• pp.127-138
• /
• 2016

Diesel vehicles should be equipped with urea-selective catalytic reduction(SCR) system as a high-performance catalyst, in order to reduce harmful nitrogen oxide emissions. In this study, a three-dimensional Eulerian-Lagrangian CFD analysis was used to numerically predict the multiphase flow characteristics of the urea-SCR system, coupled with the chemical reactions of the system's transport phenomena. Then, the numerical spray structure was modified by comparing the results with the measured values from spray visualization, such as the injection velocity, penentration length, spray radius, and sauter mean diameter. In addition, the analysis results were verified by comparison with the removal efficiency of the nitrogen oxide emissions during engine and chassis tests, resulting in accuracy of the relative error of less than 5%. Finally, a verified CFD analysis was used to calculate the interanl flow of the urea-SCR system, thereby analyzing the characteristics of pressure drop and velocity increase, and predicting the uniformity index and overdistribution positions of ammonia.