• 한국전산유체공학회지
• /
• 제14권4호
• /
• pp.101-108
• /
• 2009

We present a new finite-element scheme for direct numerical simulation of particle suspensions in simple shear flow of a viscoelastic fluid in 3D. The sliding tri-periodic representative cell concept has been combined with DEVSS/DG finite element scheme by introducing constraint equations along the domain boundary. Rigid body motion of the freely suspended particle is described by the rigid-shell description and implemented by Lagrangian multipliers on particle boundaries. We present the bulk rheology of suspensions through the numerical examples of single-, two- and many-particle problems, which represent a large number of such systems in simple shear flow. We report the steady bulk viscosity and the first normal stress coefficient, which show shear-thickening behavior for both properties.

• Korea-Australia Rheology Journal
• /
• 제17권4호
• /
• pp.171-180
• /
• 2005

We present numerical results for inertialess elliptic particle suspensions in a Newtonian fluid subject to simple shear flow, using the sliding bi-periodic frame concept of Hwang et al. (2004) such that a particulate system with a small number of particles could represent a suspension system containing a large number of particles. We report the motion and configurational change of elliptic particles in simple shear flow and discuss the inter-relationship with the bulk shear stress behaviors through several example problems of a single, two-interacting and ten particle problems in a sliding bi-periodic frame. The main objective is to check the feasibility of the direct simulation method for understanding the relationship between the microstructural evolution and the bulk material behaviors.

• Nuclear Engineering and Technology
• /
• 제37권6호
• /
• pp.511-524
• /
• 2005

Thermalhydraulic reactor simulation of tomorrow will require a new generation of codes combining at least three scales, the CFD scale in open medium, the component scale and the system scale. DNS will be used as a support for modelling more macroscopic models. NEPTUNE is such a new generation multi-scale platform developed jointly by CEA-DEN and EDF-R&D and also supported by IRSN and FRAMATOME-ANP. The major steps towards the next generation lie in new physical models and improved numerical methods. This paper presents the advances obtained so far in physical modelling for each scale. Macroscopic models of system and component scales include multi-field modelling, transport of interfacial area, and turbulence modelling. Two-phase CFD or CMFD was first applied to boiling bubbly flow for departure from nucleate boiling investigations and to stratified flow for pressurised thermal shock investigations. The main challenges of the project are presented, some selected results are shown for each scale, and the perspectives for future are also drawn. Direct Numerical Simulation tools with Interface Tracking Techniques are also developed for even smaller scale investigations leading to a better understanding of basic physical processes and allowing the development of closure relations for macroscopic and CFD models.

• 대한기계학회논문집B
• /
• 제30권5호
• /
• pp.422-429
• /
• 2006

Resolved simulations are presented fur gravity current flows aiming at studying their spreading rate. The simulations are performed for two extreme configurations such as planar and cylindrical and for 3 different Grashof numbers: $10^5,\;1{\times}10^6\;and\;10^7$. Varying the size of the heavy fluid release, the study is performed for several phases of spreading, namely acceleration, slumping and inertial phases. For the simulations, efficient spectral multi-domain code is used. From the simulations results it is concluded that 2-D results predicts well the mean front velocity during the slumping phase, but fails to predict it during the inertial phase of spreading. It is also observed that the vortex dynamics of the flow is not reproduced well by the 2-D simulation.

• 한국유변학회:학술대회논문집
• /
• 한국유변학회 2000년도 춘계 학술발표회 논문집
• /
• pp.15-18
• /
• 2000

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2008년도 12th Asian Congress of Fluid Mechanics
• /
• pp.29.1-29.1
• /
• 2008

• 유체기계공업학회:학술대회논문집
• /
• 유체기계공업학회 2006년 제4회 한국유체공학학술대회 논문집
• /
• pp.445-448
• /
• 2006

The effects of surface roughness on a spatially-developing turbulent boundary layer (TBL) were investigated by performing direct numerical simulations of TBLs over rough and smooth walls. The Reynolds number based on the momentum thickness was varied in the range $Re_{\theta}=300{\sim}1400$. The roughness elements used were periodically arranged two-dimensional spanwise rods, and the roughness height was $k=1.5{\theta}_{in}$, which corresponds to $k/{\delta}=0.045{\sim}0.125$. To avoid generating a rough wall inflow, which is prohibitively difficult, a step change from smooth to rough was placed $80{\theta}_{in}$ downstream from the inlet. The spatially-developing characteristics of the rough-wall TBL were examined. Along the streamwise direction, the friction velocity approached a constant value and a self-preserving form of the turbulent stress was obtained. Introduction of the roughness elements affected the turbulent stress not only in the roughness sublayer but also in the outer layer. Despite the roughness-induced increase of the turbulent stress in the outer layer, the roughness had only a relatively small effect on the anisotropic Reynolds stress tensor in the outer layer. Inspection of the triple products of the velocity fluctuations revealed that introducing the roughness elements onto the smooth wall had a marked effect on vertical turbulent transport across the whole TBL. By contrast, good surface similarity in the outer layer was obtained for the third-order moments of the velocity fluctuations.

• Journal of Mechanical Science and Technology
• /
• 제16권3호
• /
• pp.344-353
• /
• 2002

A reliable computational heat transfer model has been investigated to define the heat transfer characteristics of a spray column direct contact heat exchanger, which is often utilized in the process involving counterflows for heat and mass transfer operations. Most of the previous studies investigated are one-dimensional unsteady solutions based on rather fragmentary experimental data. Development of a multidimensional numerical model and a computational algorithm are essential to analyze the inherent multidimensional characteristics of a direct contact heat exchanger. The present study has been carried out numerically and establishes a solid simulation algorithm for the operation of a direct contact heat exchanger. Operational and system parameters such as the speed and direction of working fluid droplets at the injection point, and the effects of aspect ratio and void fraction of continuous fluid are examined thoroughly as well to assess their influence on the performance of a spray column.

• Wind and Structures
• /
• 제35권5호
• /
• pp.323-336
• /
• 2022

This study investigates the applicability of the direct identification of flutter derivatives in the time domain using Rational Function Approximation (RFA), where the extraction procedure requires either a combination of at least two wind speeds or one wind speed. In the frequency domain, flutter derivatives are identified at every wind speed. The ease of identifying flutter derivatives in the time domain creates a paradox because flutter derivative patterns sometimes change in higher-order polynomials. The first step involves a numerical study of RFA extractions for different deck shapes from existing bridges to verify the accurate wind speed combination for the extraction. The second step involves validating numerical simulation results through a wind tunnel experiment using the forced vibration method in one degree of freedom. The findings of the RFA extraction are compared to those obtained using the analytical solution. The numerical study and the wind tunnel experiment results are in good agreement. The results show that the evolution pattern of flutter derivatives determines the accuracy of the direct identification of RFA.

• 한국해양공학회지
• /
• 제33권5호
• /
• pp.387-393
• /
• 2019

In this study, the free motion of a riser due to vortex shedding was numerically simulated with Large Eddy Simulation (LES) and Detached Eddy Simulation (DES) turbulence models. A numerical simulation program was developed by applying the Rhie-Chow interpolation method to the pressure correction of the OpenFOAM standard solver pimpleDyMFoam. To verify the developed program, the vortex shedding around the fixed riser at Re = 3900 was calculated, and the results were compared with the existing experimental and numerical data. Moreover, the vortex-induced vibration of a riser supported by a linear spring was numerically simulated while varying the spring constant. The results are compared with published direct numerical simulation (DNS) results. The present calculation results show that the numerical method is appropriate for simulating the vortex-induced motion of a riser, including lock-in phenomena.