• Nuclear Engineering and Technology
• /
• v.55 no.9
• /
• pp.3367-3382
• /
• 2023

Smoothed Particle Hydrodynamics (SPH) is a Lagrangian computational fluid dynamics method that has been widely used in the analysis of physical phenomena characterized by large deformation or multi-phase flow analysis, including free surface. Despite the recent implementation of eddy-viscosity models in SPH methodology, sophisticated turbulent analysis using Lagrangian methodology has been limited due to the lack of computational performance and numerical consistency. In this study, we implement the standard and dynamic Smagorinsky model and dynamic Vreman model as sub-particle scale models based on a weakly compressible SPH solver. The large eddy simulation method is numerically identical to the spatial discretization method of smoothed particle dynamics, enabling the intuitive implementation of the turbulence model. Furthermore, there is no additional filtering process required for physical variables since the sub-grid scale filtering is inherently processed in the kernel interpolation. We simulate lid-driven flow under transition and turbulent conditions as a benchmark. The simulation results show that the dynamic Vreman model produces consistent results with experimental and numerical research regarding Reynolds averaged physical quantities and flow structure. Spectral analysis also confirms that it is possible to analyze turbulent eddies with a smaller length scale using the dynamic Vreman model with the same particle size.

• Computers and Concrete
• /
• v.20 no.1
• /
• pp.1-10
• /
• 2017

This paper focuses on the mesh-size dependency in numerical simulations of reinforced concrete (RC) structures subjected to blast loading. A tensile failure criterion that can minimize the mesh-dependency of simulation results is introduced based on the fracture energy theory. In addition, conventional plasticity based damage models for concrete such as the CSC model and the HJC model, which are widely used for blast analyses of concrete structures, are compared with the orthotropic model that adopts the introduced tensile failure criterion in blast tests to verify the proposed criterion. The numerical predictions of the time-displacement relations at the mid-span of RC beams subjected to blast loading are compared with experimental results. The analytical results show that the numerical error according to the change in the finite element mesh size is substantially reduced and the accuracy of the numerical results is improved by applying a unique failure strain value determined by the proposed criterion.

• Journal of Korean Society for Atmospheric Environment
• /
• v.13 no.1
• /
• pp.65-77
• /
• 1997

A theoretical and numerical investigation on the boundary-layer flow over a two- or three-dimensional hill is presented. The numerical model is based on the finite volume method with boundary-fitted coordinates. The k-$\varepsilon$ turbulence model with modified wall function and the low-Reynolds-number model are employed. The hypothesis of Reynolds number independency for the atmospheric boundary-layer flow over aerodynamically rough terrain is confirmed by the numerical simulation. Comparisons of the mean velocity profiles and surface pressure distributions between the numerical predictions and the wind-tunnel experiments on the flow over a hill show good agreement. The linear theory provides generally good prediction of speed-up characteristics for the gentle-sloped hills. The flow separation occurs in the hill slope of 0.5 and the measured reattachment points are compared with the numerical prediction. It is found that the k- $\varepsilon$ turbulence model is reasonably accurate in predicting the attached flow, while the low- Reynolds-number model is more suitable to simulate the separated flows.ows.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.21 no.1
• /
• pp.63-78
• /
• 2009

This paper presents the development of dynamically combined Typhoon generated surge-tide-wave numerical model which is applicable to shallow water. The newly developed model is based on both POM (Princeton Ocean Model) for the surge and tide and WAM (WAve Model) for wind-generated waves, but is modified to be applicable to shallow water. In this paper which is the first paper of the two in a sequence, we verified the accuracy and numerical stability of the hydrodynamic part of the model which is responsible for the simulation of Typhoon generated surge and tide. In order to improve the accuracy and numerical stability of the combined model, we modified algorithms responsible for turbulent modeling as well as vertical velocity computation routine of POM. Verification of the model performance had been conducted by comparing numerical simulation results with analytic solutions as well as data obtained from field measurement. The modified POM is shown to be more accurate and numerically stable compare to the existing POM.

• Journal of Ocean Engineering and Technology
• /
• v.17 no.4
• /
• pp.1-7
• /
• 2003

A virtual reality technology for multipurpose numerical simulation is developed to reproduce and investigate a variety of ocean environmental problems in a 3D Numerical Wave Tank(NWT). The governing equations for solving incompressible fluid motion are Navier-Stokes equation and continuity equation. The Marker-Density function technique is adopted to implement the fully nonlinear freesurface kinematic condition. The marine environmental situations, i.e., waves, currents, etc., are reproduced by use of multi-segmented wavemakers on the basis of the so-called ″snake-principle″. In this paper, some numerical reproduction techniques for regular, and irregular waves, multi-directional waves, Bull's-eye wave. wave-current, and solitary wave are presented, and a model test in motion with large amplitude of roll angle is conducted in the developed 3D-NWT, using a overlaid grid system.

• Nuclear Engineering and Technology
• /
• v.44 no.8
• /
• pp.847-854
• /
• 2012

Despite extensive research efforts, the mechanism of the nucleate boiling phenomena is still not clear. A direct numerical simulation of the boiling phenomena is one of the promising approaches in order to clarify its heat transfer characteristics and discuss their mechanism. Therefore, many DNS procedures have been developed based on recent highly advancing computer technologies. This brief review focuses on the state of the art in direct numerical simulation of the pool boiling phenomena over the past two decades. In this review, the fundamentals of the boiling phenomena and the bubble departure and micro-layer models are briefly introduced, and then the numerical procedures for tracking or capturing interface/surface shape such as the front tracking method, level set method, volume of fluid treatments, and other methods (Lattice Boltzmann method, phase-field method and so on) are briefly reviewed.

• Journal of the Korean Society of Safety
• /
• v.16 no.4
• /
• pp.74-81
• /
• 2001

A direct numerical simulation (DNS) was applied to nonpremixed counter-flow diffusion flames between oxidizer and fuel ducts. The objective of this study is to evaluate the numerical method for simulation of axisymmetric counterflow diffusion flames. Effects of computational domain size and grid size were scrutinized, and then the method was applied to air-methane diffusion flames. The results at zero gravity conditions were in good agreement with those obtained by the one-dimension flame code OPPDIF. It was confirmed thai the numerical method is applicable to the diffusion flames at the normal gravity conditions since the results clearly showed the effects of buoyancy and velocity ratio.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.18 no.4
• /
• pp.383-389
• /
• 2006

In three dimensional wave field, a direct numerical simulation model, which is able to handle free-fall and slope-fall of a waterdrop due to overtopping, is proposed to evaluate the overtopping rate on the coastal structures like an absorbing revetment. A comparison between the numerical model and existing experimental results for overtopping rate was made to validate the proposed numerical model's accuracy, and showed fairly good agreement between them. It is confirmed in numerical and hydraulic(existing) model test that the overtopping quantity on a absorbing revetment becomes larger with an increase in Ursell number. Also, the overtopping rate estimated by 3-D numerical model is compared with it obtained by 2-D numerical model.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.17 no.2
• /
• pp.86-97
• /
• 2005

Accurate estimation of the wave-induced pore water pressure in the seabed is key factor in studying the stability of the seabed in the vicinity of coastal structure. Most of the existing numerical models for wave structure seabed interaction have been linked through applying hybrid numerical technique which is analysis method separating the wave field and seabed regime. Therefore, it is necessary to develope a numerical model f3r simulating accurately wave$\cdot$structure$\cdot$ seabed interaction under wave loadings by the single domain approach for wave field and seabed regime together. In this study, direct numerical simulation is newly proposed. In this model, modeled fluid drag has been used to detect the hydraulic properties according to the varied geometrical shape inside the porous media by considering the turbulence resistance as well as laminar resistance. Contrary to hybrid numerical technique, direct numerical simulation avoids the explicit formulation of the boundary conditions at the fluid/porous media interface. A good agreement has been obtained by the comparison between existed experimental results by hydraulic model test and direct numerical simulation results far wave $\cdot$structure$\cdot$seabed interaction. Therefore, the newly proposed numerical model is a powerful tool for estimating the nonlinear dynamic responses among a structure, its seabed foundation and water waves.

• Journal of the korean Society of Automotive Engineers
• /
• v.11 no.3
• /
• pp.29-36
• /
• 1989

A multidimensional numerical simulation for flows in an engine with axisymmetric geometry was performed. Three kinds of differencing schemes, namely, skew upwind differencing scheme (SUDS), interpolated upwind differencing scheme (IUDS), upwind differencing scheme (UDS), are used in a comparative study. Simultaneously, the effects of the artificial dampings and the grids on numerical results are estimated. Compared with the measurements, the calculations with SUDS and proper artificial damping show very similar qualitative tendency with observed results. But there are some discrepancies due to numerical errors and unclear boundary conditions.