• Nuclear Engineering and Technology
• /
• v.39 no.3
• /
• pp.207-214
• /
• 2007

A new Monte Carlo method for solving heat conduction problems is developed in this study. Differing from other Monte Carlo methods, it is a transport approximation to the heat diffusion process. The method is meshless and thus can treat problems with complicated geometry easily. To minimize the boundary effect, a scaling factor is introduced and its effect is analyzed. A set of problems, particularly the heat transfer in the fuel sphere of PBMR, is calculated by this method and the solutions are compared with those of an analytical approach.

• Interaction and multiscale mechanics
• /
• v.5 no.2
• /
• pp.131-144
• /
• 2012

We study a radial basis function collocation method (RBFCM) to discretize a coupled nonlinear Schr$\ddot{o}$dinger equation (CNLSE) that governs a two dimensional rotating Bose-Einstein condensate (BEC) with an angular momentum rotation term. We exploit a RBFCM-continuation method (RBFCM-CM) to trace the solution curve of the CNLSE. We compare the performance of the RBFCM-CM with the FEM-CM. We observe that the RBFCM-CM is very robust in a coarse grid for resolving the ground state solution with many vortices when the angular momentum rotation is close to the limit. Numerical results demonstrate the efficiency and accuracy of the RBFCM-CM for computing the superfluid density of the ground level of the BEC.

• Steel and Composite Structures
• /
• v.32 no.3
• /
• pp.293-304
• /
• 2019

In the present paper, the element free Galerkin (EFG) method is developed for geometrically nonlinear analysis of deep beams considering small scale effect. To interpret the behavior of structure at the nano scale, the higher-order gradient elasticity nonlocal theory is taken into account. The radial point interpolation method with high order of continuity is used to construct the shape functions. The nonlinear equation of motion is derived using the principle of the minimization of total potential energy based on total Lagrangian approach. The Newmark method with the small time steps is used to solve the time dependent equations. At each time step, the iterative Newton-Raphson technique is applied to minimize the residential forces caused by the nonlinearity of the equations. The effects of nonlocal parameter and aspect ratio on stiffness and dynamic parameters are discussed by numerical examples. This paper furnishes a ground to develop the EFG method for large deformation analysis of structures considering small scale effects.

• Advances in nano research
• /
• v.12 no.6
• /
• pp.617-627
• /
• 2022

In the current study, the nonlinear impact of the Von-Kármán theory on the vibrational response of nonhomogeneous structures of functionally graded (FG) nano-scale tubes is investigated according to the nonlocal theory of strain gradient theory as well as high-order Reddy beam theory. The inhomogeneous distributions of temperature-dependent material consist of ceramic and metal phases in the radial direction of the tube structure, in which the thermal stresses are applied due to the temperature change in the thickness of the pipe structure. The general motion equations are derived based on the Hamilton principle, and eventually, the acquired equations are solved and modeled by the Meshless approach as well as a computer simulation via intelligent mathematical methodology. The attained results are helpful to dissect the stability of the MEMS and NEMS.

• Journal of the Korean GEO-environmental Society
• /
• v.25 no.4
• /
• pp.21-28
• /
• 2024

In this study, a polyhedral domain decomposition method for Smoothed Particle Hydrodynamics (SPH) analysis is introduced. SPH which is one of meshless methods is a numerical analysis method for fluid flow simulation. It can be useful for analyzing fluidic soil or fluid-structure interaction problems. SPH is a particle-based method, where increased particle count generally improves accuracy but diminishes numerical efficiency. To enhance numerical efficiency, parallel processing algorithms are commonly employed with the Cartesian coordinate-based domain decomposition method. However, for parallel analysis of complex geometric shapes or fluidic problems under dynamic boundary conditions, the Cartesian coordinate-based domain decomposition method may not be suitable. The introduced polyhedral domain decomposition technique offers advantages in enhancing parallel efficiency in such problems. It allows partitioning into various forms of 3D polyhedral elements to better fit the problem. Physical properties of SPH particles are calculated using information from neighboring particles within the smoothing length. Methods for sharing particle information physically separable at partitioning and sharing information at cross-points where parallel efficiency might diminish are presented. Through numerical analysis examples, the proposed method's parallel efficiency approached 95% for up to 12 cores. However, as the number of cores is increased, parallel efficiency is decreased due to increased information sharing among cores.

• Journal of the Korean GEO-environmental Society
• /
• v.25 no.4
• /
• pp.13-20
• /
• 2024

In this paper, a parallel analysis algorithm for Smoothed Particle Hydrodynamics (SPH), one of the numerical methods for fluidic materials, is introduced. SPH, which is a meshless method, can represent the behavior of a continuum using a particle-based approach, but it demands substantial computational resources. Therefore, parallel analysis algorithms are essential for SPH simulations. The domain decomposition algorithm, which divides the computational domain into partitions to be independently analyzed, is the most representative method among parallel analysis algorithms. In Discrete Element Method (DEM) and Molecular Dynamics (MD), the Cartesian coordinate-based domain decomposition method is popularly used because it offers advantages in quickly and conveniently accessing particle positions. However, in SPH, it is important to share particle information among partitioned domains because SPH particles are defined based on information from nearby particles within the smoothing length. Additionally, maintaining CPU load balance is crucial. In this study, a highly parallel efficient algorithm is proposed to dynamically minimize the size of orthogonal domain partitions to prevent excess CPU utilization. The efficiency of the proposed method was validated through numerical analysis models. The parallel efficiency of the proposed method is evaluated for up to 30 CPUs for fluidic models, achieving 90% parallel efficiency for up to 28 physical cores.

• Structural Engineering and Mechanics
• /
• v.14 no.6
• /
• pp.713-732
• /
• 2002

A local radial point interpolation method (LRPIM) based on local residual formulation is presented and applied to solid mechanics in this paper. In LRPIM, the trial function is constructed by the radial point interpolation method (PIM) and establishes discrete equations through a local residual formulation, which can be carried out nodes by nodes. Therefore, element connectivity for trial function and background mesh for integration is not necessary. Radial PIM is used for interpolation so that singularity in polynomial PIM may be avoided. Essential boundary conditions can be imposed by a straightforward and effective manner due to its Delta properties. Moreover, the approximation quality of the radial PIM is evaluated by the surface fitting of given functions. Numerical performance for this LRPIM method is further studied through several numerical examples of solid mechanics.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.31 no.4
• /
• pp.35-43
• /
• 2003

In this paper, an efficient computational algorithm for the implementation of the newly proposed node activation technique is presented, and its computational aspects are thoroughly investigated. To verify the validity, convergence, and efficiency of the node activation technique, various numerical examples are worked out including the problems of Poisson equation, 2D elasticity problems, and 3D elasticity problems. From the numerical tests, it is verified that one can arbitrarily activate and handle the nodal points of interest in finite element model with very little loss of the numerical accuracy.

• Journal of Advanced Research in Ocean Engineering
• /
• v.1 no.1
• /
• pp.1-13
• /
• 2015

Post-Katrina investigations revealed that most earthen levee damage occurred on the levee crest and landward-side slope as a result of either wave overtopping, storm surge overflow, or a combination of both. In this paper, combined wave overtopping and storm surge overflow of a levee embankment strengthened with high performance turf reinforcement mat (HPTRM) system was studied in a purely Lagrangian and meshless approach, two-dimensional smoothed particle hydrodynamics (SPH) model. After the SPH model is calibrated with full-scale overtopping test results, the overtopping discharge, flow thickness, flow velocity, average overtopping velocity, shear stress, and soil erosion rate are calculated. New equations are developed for average overtopping discharge. The shear stresses on landward-side slope are calculated and the characteristics of soil loss are given. Equations are also provided to estimate soil loss rate. The range of the application of these equations is discussed.

• Structural Engineering and Mechanics
• /
• v.21 no.2
• /
• pp.151-168
• /
• 2005

In this paper an energy method is presented for the linear buckling analysis of first order shear deformable plates. The displacement fields are defined in terms of the shape functions, which correspond to a set of predefined points and are composed of significantly high order polynomials. The locations of these points are found by mapping the geometry using the naturalized coordinates and bilinear shape functions. In order to evaluate the method, fully clamped and simply supported rectangular plates subjected to uniform uniaxial compressive loading on two opposite edges of the plate are investigated thoroughly and the results are compared with the exact solution given in the monograph of Timoshenko and Gere (1961). The method is extended to the analysis of perforated plates, wherein the negative stiffness computed over the opening area from in-plane and out-of-plane deformation modes is superimposed to the stiffness of the full plate. Numerical results are then favorably compared with those obtained by finite element methods. Other cases such as; rectangular plates with eccentrically located openings of different shapes are studied and reported in this paper with regards to the effect of aspect ratio, hole size, and hole position on the buckling. For a square plate with a large circular opening at the center, diameter being 80 percent of the length, the present method yields buckling coefficient 12.5 percent higher than the one from the FEM.