• 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.

• Journal of The Korean Astronomical Society
• /
• v.47 no.3
• /
• pp.87-98
• /
• 2014

We develop a parallel cosmological hydrodynamic simulation code designed for the study of formation and evolution of cosmological structures. The gravitational force is calculated using the TreePM method and the hydrodynamics is implemented based on the smoothed particle hydrodynamics. The initial displacement and velocity of simulation particles are calculated according to second-order Lagrangian perturbation theory using the power spectra of dark matter and baryonic matter. The initial background temperature is given by Recfast and the temperature uctuations at the initial particle position are assigned according to the adiabatic model. We use a time-limiter scheme over the individual time steps to capture shock-fronts and to ease the time-step tension between the shock and preshock particles. We also include the astrophysical gas processes of radiative heating/cooling, star formation, metal enrichment, and supernova feedback. We test the code in several standard cases such as one-dimensional Riemann problems, Kelvin-Helmholtz, and Sedov blast wave instability. Star formation on the galactic disk is investigated to check whether the Schmidt-Kennicutt relation is properly recovered. We also study global star formation history at different simulation resolutions and compare them with observations.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.15 no.1
• /
• pp.1-10
• /
• 1991

Two-dimensional explicit finite element code was developed. The transient dynamics code can analyse large deformations of non-linear materials subjected to extremely high strain rates. The Lagrangian finite element program uses an explicit time integration operator to integrate the equations of motion, thus the stiffness matrix is not introduced. Cylinder upsetting and ring compression problems are simulated to check the effects of friction and inertia force. It is shown that (1) calculated results agree very well with experimental results, (2) constant shear friction method overestimates the decrease of inner ring radius and then underestimates after on in comparison with the Coulomb friction method, and (3) the effect of the increase in initial strain rate is similar to the effect of higher frictional coefficient.

• Journal of the Society of Naval Architects of Korea
• /
• v.42 no.4 s.142
• /
• pp.330-340
• /
• 2005

Compressible two-phase flow is analyzed based on the arbitrary Lagrangian-Eulerian (ALE) formulation. For water, Tamman type stiffened equation of state is used. Numerical fluxes are calculated using the ALE two-phase Godunov scheme which assumes only that the speed of sound and pressure can be provided whenever density and internal energy are given. Effects of the approximations of a material interface speed are Investigated h method Is suggested to assign a rigid body boundary condition effectively To validate the developed code, several well-known problems are calculated and the results are compared with analytic or other numerical solutions including a single material Sod shock tube problem and a gas/water shock tube problem The code is applied to analyze the refraction and transmission of shock waves which are impacting on a water-gas interface from gas or water medium.

• 한국전산유체공학회:학술대회논문집
• /
• 1998.05a
• /
• pp.1-14
• /
• 1998

Underwater explosion properties for TNT, an ideal explosive, and DXD-04, a nonideal explosive, were numerically modeled with a one-dimensional Lagrangian hydrodynamic code. The equation of state parameters for detonation products for TNT and DXD-04 were obtained from the BKW code, assuming complete reaction. Burn of TNT was modeled by using the Chapman-Jouguet(CJ) volume burn technique, a programmed-burn technique, assuming instantaneous detonation reaction. Burn of DXD-04 was modeled by using the same technique and by using the reaction rate calibrated from two-dimensional steady-state detonation experiments. The calculations for TNT reproduced the experimental peak pressure of the shock wave propagating through water with an error of $3.0\%$ and the experimental oscillation period of the bubble formed of detonation products with an error of $2.3\%$. For DXD-04, the CJ volume burn technique could not reproduce the experimental observations. When the reaction rate calibrated from two-dimensional steady-state detonation experimental data, the calculated peak pressure was slightly higher by $7.3\%$ than the experimental data, but the calculated shock profile was in good agreement. The bubble period was reproduced with an error of $1.8\%$. These results demonstrated that underwater explosion properties for an ideal explosive can be predicted by using a programmed burn technique, and that, however, those for a nonideal explosive can be predicted only when a well-calibrated reaction rate is used.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.43 no.1
• /
• pp.8-22
• /
• 2015

Most jet spray and atomization simulations are done with the Eulerian method which has inherent disadvantage in representing jet breakups and droplets. Full Lagrangian particles method called Smoothed Particle Hydrodynamics(SPH) is used in this work. We develop the SPH code and perform validations that confirm the suitability of our SPH method for simulating liquid jet atomization problem. Then, we conduct the simulation of liquid-liquid swirl coaxial injector for comparison against the experimental data.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2004.04a
• /
• pp.111-118
• /
• 2004

For the precise analysis of high velocity impact problem though FEM with element erosive method, the adequate mesh size and critical equivalent plastic strain(EQPS) is chosen prior to the simulation. In this research, it is strongly required from a standpoint that critical EQPS is used to decide whether perforation occurs or not. The optimization of dual armor plate consisting of 4340 steel and 2024 aluminium against a die steel sphere with high-velocity has been suggested using Lagrangian explicit time-integration code, NET2D. The response surface method based on the design of experiment is utilized for the size optimization. The optimized thickness of each layer, in which perforation does not occur, the strength of multi-layer is maximized and total weight is minimized, is obtained at a constant velocity of a pellet with a designated total thickness.

• Journal of the Korean Society of Manufacturing Technology Engineers
• /
• v.8 no.1
• /
• pp.60-66
• /
• 1999

This study is to develop design sensitivity analysis method based on continuum theory for the actual buckling load of vehicle passenger compartment with respect to sizing design variables. For nonlinear structural analysis, both geometric and material nonlinear effects are considered. The total Lagrangian formulation for incremental equilibrium analysis and one-point linear eigenvalue problem for buckling analysis are utilized. Numerical methods are presented to evaluate design sensitivity expressions, using structural analysis results from FEM code. Optical design of vehicle passenger compartment with buckling constraint solved using Gradient projection method.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.10 no.6
• /
• pp.702-713
• /
• 1998

The characteristics of flow and particle collection for an electrocyclone with a central wire inside a high efficiency Stairmand cyclone was numerically analysed. Turbulent flow field was modeled by the Reynolds stress model and solved with an FVM code FLUENT. Particle motion and in-situ charging were simultaneously solved by a Lagrangian integration with time. The flow field obtained was in good agreement with experiments in the outer region. The characteristics of collection enhancement due to electric force were well manifested and well explained based on first principles. The effect of the in-situ charging process was very similar to the case of a simplified assumption of saturated charging, and the effect of the hopper was proved negligible.

• 한국전산유체공학회:학술대회논문집
• /
• 2010.05a
• /
• pp.460-463
• /
• 2010

Ice accretion on aircraft surface can greatly impair the aerodynamic performance of aircraft. As an alternative to the traditional Lagrangian particle tracking approach, an Eulerian-based droplet impingement and ice accretion code for air flows containing water droplets was developed A CFD solver was also developed to solve the clean airflow. The results of present method were compared with experimental data and previous icing codes such as LEWICE and FENSAP-ICE and were confirmed to show good agreement each other in qualitative and quantitative ways.