• Structural Engineering and Mechanics
• /
• v.78 no.4
• /
• pp.497-518
• /
• 2021

In this study, a new method for treating the wall boundary in smoothed particle hydrodynamics (SPH) is proposed to simulate free surface flows effectively. Unlike conventional methods of wall boundary treatment through boundary particles, in the proposed method, the wall boundary condition is directly imposed by adding boundary truncation terms to the mass and momentum conservation equations. Thus, boundary particles are not used in boundary modeling. Doing so, the wall boundary condition is accurately imposed, boundary modeling is simplified, and computation is made efficient without losing stability in SPH. Performance of the proposed method is demonstrated through several numerical examples: dam break, dam break with a wedge, sloshing, inclined bed, cross-lever rotation, pulsating tank and sloshing with a flexible baffle. These results are compared with available experimental results, analytical solutions, and results obtained using the boundary particle method.

• Journal of Mechanical Science and Technology
• /
• v.15 no.5
• /
• pp.671-680
• /
• 2001

Boundary condition is one of the major factors to influence the numerical stability and solution accuracy in numerical analysis. One of the most important physical boundary conditions in the flowfield analysis is the wall boundary condition imposed on the body surface. To solve a two-dimensional Euler equation, totally four numerical wall boundary conditions should be prescribed. Two of them are supplied by the flow tangency condition. The other two conditions, therefore, should be prepared additionally in a suitable way. In this paper, four different sets of wall boundary conditions are proposed and then applied to solve high-speed flowfields around a quarter circle geometry. A two-dimensional compressible Euler solver is prepared based on the finite volume method. This solver hires three different upwind schemes; Steger-Warmings flux vector splitting, Roes flux difference splitting, and Lious advection upstream splitting method. It is found that the way to specify the additional numerical wall boundary conditions strongly affects the overall stability and accuracy of the upwind schemes in high-speed flow calculation. The optimal wall boundary conditions should be also chosen very carefully depending on the numerical schemes used to solve the problem.

• Journal of Mechanical Science and Technology
• /
• v.19 no.11
• /
• pp.2077-2081
• /
• 2005

A symmetry breaking nonlinear fluid flow in a two-dimensional wall-driven square cavity taking symmetric boundary condition after some transients has been investigated numerically. It has been shown that the symmetry breaking critical Reynolds number is dependent on the time history of the boundary condition. The cavity has at least three stable steady state solutions for Re=300-375, and two stable solutions if Re>400. Also, it has also been showed that a particular solution among several possible solutions can be obtained by a controlled boundary condition.

• Proceedings of the KSME Conference
• /
• 2007.05b
• /
• pp.2797-2802
• /
• 2007

It has been confirmed that implementation of the no-slip boundary conditions for the lattice-Boltzmann method play an important role in the overall accuracy of the numerical solutions as well as the stability of the solution procedure. We in this paper propose a new algorithm, i.e. the method of the dynamic boundary condition for no-slip boundary condition. The distribution functions on the wall along each of the links across the physical boundary are assumed to be composed of equilibrium and nonequilibrium parts which inherit the idea of Guo's extrapolation method. In the proposed algorithm, we apply a dynamic equation to reflect the computational slip velocity error occurred on the actual wall boundary to the correction; the calculated slip velocity error dynamically corrects the fictitious velocity on the wall nodes which are subsequently employed to the computation of equilibrium distribution functions on the wall nodes. Along with the dynamic selfcorrecting process, the calculation efficiently approaches the steady state. Numerical results show that the dynamic boundary method is featured with high accuracy and simplicity.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2005.09b
• /
• pp.902-903
• /
• 2005

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.23 no.1
• /
• pp.33-42
• /
• 1999

The effects of wan thermal boundary condition were investigated for a uniform wall temperature, a uniform wall heat flux, and for coupled heat conduction In the channel wall with transverse rectangular ribs. Numerical investigations for steady laminar flow show behavior similar to that observed experimentally in the separated flow region for flow over a cylinder. Conjugate heat transfer with a low solid-fluid thermal conductivity ratio does not lead to the same results as for the uniform heat flux boundary condition, and heat transfer reversal is found on the back sides of the ribs.

• International Journal of High-Rise Buildings
• /
• v.11 no.4
• /
• pp.331-346
• /
• 2022

Reproducing the horizontally homogeneous atmospheric boundary layer in computational wind engineering is essential for predicting the wind loads on structures. One of the important issues is to use fully developed inflow conditions, which will lead to the consistence problem between inflow condition and internal roughness. Thus, by analyzing the previous results of computational fluid dynamic modeling turbulent horizontally homogeneous atmospheric boundary layer, we modify the past hypotheses, detailly derive a new type of inflow condition for standard k-ε turbulence model. A group of remedial approaches including formulation for wall shear stress and fixing the values of turbulent kinetic energy and turbulent dissipation rate in first wall adjacent layer cells, are also derived to realize the consistence of inflow condition and internal roughness. By combing the approaches with four different sets of inflow conditions, the well-maintained atmospheric boundary layer flow verifies the feasibility and capability of the proposed inflow conditions and remedial approaches.

• Journal of computational fluids engineering
• /
• v.17 no.2
• /
• pp.42-52
• /
• 2012

Direct Numerical Simulation(DNS) of turbulent mass transfer in fully developed turbulent pipe flow has been performed to study the effect of wall boundary conditions on the concentration fields at $Re_{\tau}$=180 based on friction velocity and pipe radius. Fully developed turbulent pipe flows for Sc=0.71 are studied with two different wall boundary conditions, namely, constant mass flux and constant wall concentration. The mean concentration profiles and turbulent mass fluxes obtained from the present DNS are in good agreement with the previous numerical results currently available. To investigate the effects of wall boundary condition on the turbulent mass transfer, the mean concentration profile, root-mean-square of concentration fluctuation, turbulent mass fluxes and higher-order statistics(Skewness and Flatness factor) are compared for the two cases. Furthermore, the budgets of turbulent mass fluxes and concentration variance were computed and analyzed to elucidate the effects of wall boundary conditions on the turbulent mass transfer.

• Journal of computational fluids engineering
• /
• v.14 no.4
• /
• pp.23-30
• /
• 2009

The existing code which solves two-dimensional RANS(Reynolds Averaged Navier-Stokes) equations and 2-equation turbulence model equations was modified to enable numerical simulation of various supersonic flows. For this, various boundary conditions have been implemented to the code. Bleed boundary condition was incorporated into the code for calculating wall mean flow quantities. Furthermore, boundary conditions for the turbulence quantities along rough surfaces as well as porous walls were applied to the code. The code was verified and validated by comparing the computational results against the experimental data for the supersonic flows over bleed region on a flat plate. Furthermore, numerical simulations for supersonic shock boundary layer interaction with a bleed region were performed and their results were compared with the existing computational results.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.12
• /
• pp.1635-1643
• /
• 2000

The present work investigates a heat transfer phenomenon at the interface between a porous medium and an impermeable wall. In an effort to appropriately describe the heat transfer phenomenon at the interface, the heat transfer at the interface between the microchannel heat sink, which is an ideally organized porous medium, and the finite-thickness substrate is examined. From the examination, it is clarified that the he heat flux distribution at the interface is not uniform for the impermeable wall with finite thickness. On the other hand, the first approach, based on the energy balance for the representative elementary volume in the porous medium, is physically reason able. When the first approach is applied to the thermal boundary condition, and additional boundary condition based on the local thermal equilibrium assumption at the interface is used. This additional boundary condition is applicable except for the very th in impermeable wall. Hence, for practical situations, the first approach in combination with the local thermal equilibrium assumption at the interface is suggested as an appropriate thermal boundary condition. In order to confirm our suggestion, convective flows both in a microchannel heat sink and in a sintered porous channel subject to a constant heat flux condition are analyzed. The analytically obtained thermal resistance of the microchannel heat sink and the numerically obtained overall Nusselt number for the sintered porous channel are shown to be in close agreement with available experimental results when our suggestion for the thermal boundary conditions is applied.