• Journal of the Korea Computer Graphics Society
• /
• v.23 no.1
• /
• pp.9-16
• /
• 2017

Simulating the clumping and stiffness of wet hair or fur is a challenging problem. The dynamics of wet hair or fur is characterized by the clumping and stiffness at the tip, which is easily seen in running animals or headbanging scenes. Existing methods address these phenomenon within pre-set scenarios. But there is no consensus on the method of depicting the details of wet hair. Hence, the present paper proposes a new method of modeling the clumping and stiffness of wet hair or fur. Previous studies focused on modeling the absorption of water into hair or fur, whereas this paper highlights a realistic simulation of wet hair. Unlike dry hair strands, wet hair strands adjacent to one another are subjected to the clumping force and gather together, while at the same time becoming stiff as the saturation of water increases. The proposed method builds on the surface tension model based on SPH (smoothed particle hydrodynamics) to simulate the clumping force and to adjust the hair elasticity by giving stiffness constraints. The present method enables a realistic simulation of wet hair by maintaining the clumping force of the wet hair even in dynamic motions, and by simulating the stiffness of hair in line with water saturation.

• The Bulletin of The Korean Astronomical Society
• /
• v.38 no.2
• /
• pp.101.2-101.2
• /
• 2013

We present a new computational fluid dynamic (CFD) simulation code. The code employs the moving and polyhedral unstructured mesh scheme, which is known as a superior approach to the conventional SPH (smoothed particle hydrodynamics) and AMR (adaptive mesh refinement) schemes. The code first generates unstructured meshes by the Voronoi tessellation at every time step, and then solves the Riemann problem for surfaces of every Voronoi cell to update the hydrodynamic states as well as to move former generated meshes. For the second-order accuracy, the MUSCL-Hancock scheme is implemented. To increase efficiency for generating Voronoi tessellation we also develop the incremental expanding method, by which the CPU time is turned out to be just proportional to the number of particles, i.e., O(N). We will discuss the applications of our code in the context of cosmological simulations as well as numerical experiments for galaxy formation.

• The Bulletin of The Korean Astronomical Society
• /
• v.39 no.2
• /
• pp.65.2-65.2
• /
• 2014

We present a new hydrodynamic simulation code based on the Voronoi tessellation for estimating the density precisely. The code employs both of Lagrangian and Eulerian description by adopting the movable mesh scheme, which is superior to the conventional SPH (smoothed particle hydrodynamics) and AMR (adaptive mesh refinement) schemes. The code first generates unstructured meshes by the Voronoi tessellation at every time step, and then solves the Riemann problem for all surfaces of each Voronoi cell so as to update the hydrodynamic states as well as to move current meshes. Besides, the IEM (incremental expanding method) is devised to compute the Voronoi tessellation to desired degree of speed, thereby the CPU time is turned out to be just proportional to the number of particles, i.e., O(N). We discuss the applications of our code in the context of cosmological simulations as well as numerical experiments for galaxy formation.

• The Bulletin of The Korean Astronomical Society
• /
• v.36 no.2
• /
• pp.57.2-57.2
• /
• 2011

To study the formation and evolution of sub-galactic scale structures, we have added SPH (Smoothed Particle Hydrodynamics) method into an existing cosmological PMTree code, GOTPM. To follow the evolution of gas particles, we consider heating/cooling processes, star formation, and energy & metal feedback by supernova explosion. We have performed various tests for the new code and found that the results reproduce observed quantities or follow the known analytic solutions. We present a test simulation of isolated disk galaxy with a focus on whether the star formation reproduces the observed features.

• Journal of The Korean Astronomical Society
• /
• v.36 no.3
• /
• pp.241-248
• /
• 2003

Owing to several observational evidences and theoretical predictions for morphological evolution of galaxies, it is now widely accepted that galaxies do evolve from late types to early ones along the Hubble sequence. It is also well established that non-axisymmetric potentials of bar-like or oval mass distributions can change the morphology of galaxies significantly during the Hubble time. Here, we review the observational and theoretical grounds of the secular evolution driven by bar-like potentials, and present the results of SPH simulations for the response of the gaseous disks to the imposed potentials to explore the secular evolution in the central regions of barred galaxies.

• 한국전산유체공학회:학술대회논문집
• /
• 2009.04a
• /
• pp.239-244
• /
• 2009

The gridless (or meshfree) methods, such as MPS, SPH, FPM an so forth, are feasible and robust for the problems with moving boundary and/or complicated boundary shapes, because these methods do not need to generate a grid system. In this study, a gridless solver, which is based on the combination of moving least square interpolations on a cloud of points with point collocation for evaluating the derivatives of governing equations, is presented for two-dimensional unsteady incompressible Navier-Stokes problem in the low Reynolds number. A MAC-type algorithm was adopted and the Poission equation for the pressure was solved by successively in the moving least square sense. Some weighing functions were tested in order to investigate the up-winding effect for the convection term. Some typical problems were solved by the presented solver for the validation and the results obtained were compared with analytic solutions and the numerical results by conventional CFD methods, such as FVM.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.31 no.1A
• /
• pp.25-33
• /
• 2011

The impact damage behavior of steel-strengthened concrete panels exposed to explosive loading is investigated. Since real explosion experiments require the vast costs to facilities as well as the blast and impact damage mechanisms are too complicated, numerical analysis has lately become a subject of special attention. However, for engineering problems involving blast wave and fragment impact, there is no single numerical method that is appropriate to the various problems. In order to evaluate the retrofit performance of a steel-strengthened concrete panel subject to blast wave and fragment impact loading, an explicit analysis program, AUTODYN is used in this work. The multi-solver coupling methods such as Euler-Lagrange and SPH-Lagrange coupling method in order to improve efficiency and accuracy of numerical analysis is implemented. The simplified and idealized two dimensional and axisymmetric models are used in order to obtain a reasonable computation running time. As a result of the analysis, concrete panels subject to either blast wave or fragment impact loading without the steel plate are shown the scabbing and perforation. The perforation can be prevented by concrete panels reinforced with steel plate. The numerical results show good agreement with the results of the experiments.

• Nuclear Engineering and Technology
• /
• v.53 no.5
• /
• pp.1403-1415
• /
• 2021

The transient capability of a SP₃ based pin-wise core analysis code SPHINCS is developed and verified through the analyses of the C5G7-TD benchmark. Spatial discretization is done by the fine mesh finite difference method (FDM) within the framework of the coarse mesh finite difference (CMFD) formulation. Pin size fine meshes are used in the radial fine mesh kernels. The time derivatives of the odd moments in the time-dependent SP₃ equations are neglected. The pin homogenized group constants and Super Homogenization (SPH) factors generated from the 2D single assembly calculations at the unrodded and rodded conditions are used in the transient calculations via proper interpolation involving the approximate flux weighting method for the cases that involve control rod movement. The simplifications and approximations introduced in SPHINCS are assessed and verified by solving all the problems of C5G7-TD and then by comparing with the results of the direct whole core calculation code nTRACER. It is demonstrated that SPHINCS yields accurate solutions in the transient behaviors of core power and reactivity.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.41 no.11
• /
• pp.1065-1075
• /
• 2017

In this study, high velocity impact tests along with modeling of material behavior and numerical analyses were conducted to predict the dispersion behavior of the debris resulting from a high velocity impact fracture. For the impact tests, two different materials were employed for both the projectile and the target plate - the first setup employed aluminum alloy while the second employed steel. The projectile impacts the target plate with a velocity of approximately 1 km/s were enforced to generate the impact damages in the aluminum witness plate through the fracture debris. It was confirmed that, depending on the material employed, the debris dispersion behavior as well as the dispersion radii on the witness plate varied. A numerical analysis was conducted for the same impact test conditions. The smoothed particle hydrodynamics (SPH)-finite element (FE) coupled technique was then applied to model the fracture and damage upon the debris. The experimental and numerical results for the diameters of the perforation holes in the target plate and the debris dispersion radii on the witness plate were in agreement within a 5％ error. In addition, the impact test using steel was found to be more threatening as proven by the larger debris dispersion radius.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.50 no.12
• /
• pp.829-838
• /
• 2022

This study conducts Hypervelocity Impact(HVI) simulations considering space objects with various shapes and different impact angles. A commercial nonlinear structural dynamics analysis code, LS-DYNA, is used for the present simulation study. The Smoothed Particle Hydrodynamic(SPH) method is applied to represent the impact phenomena with hypervelocity. Mie-Grüneisen Equation of State and Johnson-Cook material model are used to consider nonlinear structural behaviors of metallic materials. The space objects with various shapes are modeled as a sphere, cube, cylinder, and cone, respectively. The space structure is modeled as a thin plate(200 mm×200 mm×2 mm). HVI simulations are conducted when space objects with various shapes with 4.119 km/s collide with the space structures, and the impact phenomena such as a debris cloud are analyzed considering the space objects with various shapes having the same mass at the different impact angles of 0°, 30° and 45° between the space object and space structure. Although space objects have the same kinetic energy, different debris clouds are generated due to different shapes. In addition, it is investigated that the size of the debris cloud is decreased by impact angles.