• International Journal of High-Rise Buildings
• /
• v.8 no.4
• /
• pp.255-264
• /
• 2019

This paper describes large eddy simulation of wind pressures on a square cylinder in a uniform flow and a high-rise building immersed in an atmospheric turbulent boundary layer. For the atmospheric boundary layer case, the inflow turbulence is generated by a numerical wind tunnel. In the numerical simulation, particular attention is devoted to the performance of an auto hexahedral non-structural mesh. Both simulations are performed for three grid systems: an auto hexahedral non-structured grid, a structured Cartesian grid and a non-structured triangular prism grid, and for three grid numbers. The present study shows that the auto hexahedral unstructured mesh achieves the best simulation results for wind pressures on the square cylinder and the high-rise building. When the grid number is sufficiently large, the differences among the results obtained from the three investigated grid systems are not significant. However, the advantage of the auto hexahedral unstructured mesh becomes clear when the grid number decreases, because it enables a balanced distribution of orthogonal grids. The results described in this paper demonstrate that the auto hexahedral non-structured mesh has good potential applicability to simulation of urban flows.

• Journal of computational fluids engineering
• /
• v.9 no.3
• /
• pp.73-80
• /
• 2004

A LS(Level Set) formulation is developed for computing two-phase flows on non- orthogonal meshes. Compared with the VOF(Volume-of-Fluid) method based on a non-smooth volume-fraction function, the LS method can calculate an interfacial curvature more accurately by using a smooth distance function. Also, it is quite straightforward to implement for 3-D irregular meshes compared with the VOF method requiring much more complicated geometric calculations. The LS formulation is implemented into a general purpose program for 3-D flows and verified through several test problems.

• Nuclear Engineering and Technology
• /
• v.51 no.5
• /
• pp.1181-1194
• /
• 2019

A new moving-mesh Finite Volume Method (FVM) for the efficient solution of the two-dimensional neutron diffusion equation is introduced. Many other moving-mesh methods developed to solve the neutron diffusion problems use a relatively large number of sophisticated mathematical equations, and so suffer from a significant complexity of mathematical calculations. In this study, the proposed method is formulated based on simple mathematical algebraic equations that enable an efficient mesh movement and CV deformation for using in practical nuclear reactor applications. Accordingly, a computational framework relying on a new moving-mesh FVM is introduced to efficiently distribute the meshes and deform the CVs in regions with high gradient variations of reactor power. These regions of interest are very important in the neutronic assessment of the nuclear reactors and accordingly, a higher accuracy of the power densities is required to be obtained. The accuracy, execution time and finally visual comparison of the proposed method comprehensively investigated and discussed for three different benchmark problems. The results all indicated a higher accuracy of the proposed method in comparison with the conventional fixed-mesh FVM.

• Proceedings of the Korea Association of Crystal Growth Conference
• /
• 1996.06b
• /
• pp.285-309
• /
• 1996

Non oxide ceramics such as nitrides of transition metals have shown significant potential for future economic impact, in diverse applications in ceramic, aerospace and electronic industries, as refractory products, abrasives and cutting tools, aircraft components, and semi-conductor substrates amid others. Combustion synthesis has become an attractive alternative to the conventional furnace technology to produce these materials cheaply, faster and at a higher level of purity. However he process os highly exothermic and manifests complex dynamics due to its strongly non-linear nature. In order to develop an understanding of this process and to study the effect of operational parameters on the final outcome, numerical modeling is necessary, which would generated essential knowledge to help scale-up the process. the model is based on a system of parabolic-hyperbolic partial differential equations representing the heat, mass and momentum conservation relations. The model also takes into account structural change due to sintering and volumetric expansion, and their effect on the transport properties of the system. The solutions of these equations exhibit steep moving spatial gradients in the form of reaction fronts, propagating in space with variable velocity, which gives rise to varying time scales. To cope with the possibility of extremely abrupt changes in the values of the solution over very short distances, adaptive mesh techniques can be applied to resolve the high activity regions by ordering grid points in appropriate places. To avoid a control volume formulation of the solution of partial differential equations, a simple orthogonal, adaptive-mesh technique is employed. This involves separate adaptation in the x and y directions. Through simple analysis and numerical examples, the adaptive mesh is shown to give significant increase in accuracy in the computations.

• Journal of Mechanical Science and Technology
• /
• v.16 no.2
• /
• pp.182-191
• /
• 2002

The reliable stress intensity factor analysis is required for fracture mechanics design or safety evaluation of mechanical joints at which cracks often initiate and grow. It has been reported that cracks in mechanical joints usually nucleate as corner cracks at the faying surface of joints and grow as elliptical arc through cracks. In this paper, three dimensional finite element analyses are performed for elliptical arc through cracks in mechanical joints. Thereafter stress intensity factors along elliptical crack front including two surface points are determined by the virtual crack closure technique. Virtual crack closure technique is a method to calculate stress intensity factor using the finite element analysis and can be applied to non-orthogonal mesh. As a result, the effects of clearance on the stress intensity factor are investigated and crack shape are then predicted.