• Structural Engineering and Mechanics
• /
• v.17 no.5
• /
• pp.671-690
• /
• 2004

A meshless method with novel variation of point collocation by finite mixture approximation is developed in this paper, termed the meshless finite mixture (MFM) method. It is based on the finite mixture theorem and consists of two or more existing meshless techniques for exploitation of their respective merits for the numerical solution of partial differential boundary value (PDBV) problems. In this representation, the classical reproducing kernel particle and differential quadrature techniques are mixed in a point collocation framework. The least-square method is used to optimize the value of the weight coefficient to construct the final finite mixture approximation with higher accuracy and numerical stability. In order to validate the developed MFM method, several one- and two-dimensional PDBV problems are studied with different mixed boundary conditions. From the numerical results, it is observed that the optimized MFM weight coefficient can improve significantly the numerical stability and accuracy of the newly developed MFM method for the various PDBV problems.

• Journal of Ocean Engineering and Technology
• /
• v.28 no.2
• /
• pp.93-101
• /
• 2014

In general, particle simulation methods such as the MPS(Moving Particle Simulation) or SPH(Smoothed Particle Hydrodynamics) methods have some serious drawbacks for pressure solutions. The pressure field shows spurious high fluctuations both temporally and spatially. It is well known that pressure fluctuation primarily occurs because of the numerical approximation of the partial differential operators. The MPS and SPH methods employ a pre-defined kernel function in the approximation of the gradient and Laplacian operators. Because this kernel function is constructed artificially, an accurate solution cannot be guaranteed, especially when the distribution of particles is irregular. In this paper, we propose a particle simulation method based on the moving least-square technique for solving the partial differential operators using a Taylor-series expansion. The developed method was applied to the hydro-static pressure and dam-broken problems to validate it.

• Atmosphere
• /
• v.30 no.1
• /
• pp.75-89
• /
• 2020

Current in-situ airborne probes that measure the sizes of ice crystals smaller than 50 ㎛ are based on the concept that the measured intensity of light scattered by a particle in the forward and/or backward direction can be converted to particle size. The relationship between particle size and scattered light used in forward scattering probes is based on Mie theory, which assumes the refractive index of particle is known and all particles are spherical. Not only are small crystals not spherical, but also there are a wide variety of non-spherical shapes. Although it is well known that the scattering properties of non-spherical ice crystals differ from those of spherical shapes, the impacts of non-sphericity on derived in-situ particle size distributions are unknown. Thus, precise relationships between the intensity of scattered light and particle size and shape are required, as based on accurate calculations of scattering properties of ice crystals. In this study, single-scattering properties of ice crystals smaller than 50 ㎛ are calculated at a wavelength of 0.55 ㎛ using a numerically exact method (i.e., discrete dipole approximation). For these calculations, hexagonal ice crystals with varying aspect ratios are used to represent the shapes of natural small ice crystals to determine the errors caused by non-spherical ice crystals measured by forward scattering probes. It is shown that the calculated errors in sizing nonspherical ice crystals are at least 13% and 26% in forward (4~12°) and backward (168~176°) directions, respectively, and maximum errors are up to 120% and 132%.

• Journal of Aerospace System Engineering
• /
• v.7 no.1
• /
• pp.19-25
• /
• 2013

Smoothed Particle Hydrodynamics(SPH) method uses a grid of historical analysis and is not Lagrangian particles using the grid method. The Navier-Stokes equations were used to solve the viscous flow of the non-compressed. In this study, the numerical analysis of the three-dimensional Coin Drop Simulation using SPH method was performed, and the analysis results are compared with experimental results, and a similar behavior can be seen. The commercial program used was Abaqus/Explicit. SPH method to reduce the error by comparing the existing flow analysis or interpretation of the continuing research is needed in the future. That will enable real-time analysis of material obtained as a result of these numerical simulations similar to the actual flow phenomena, depending on the development of computer graphics technology to show visually. As a result, this method can be applied to the analysis fluid - structure interaction problems in a variety of fields.

• Journal of Advanced Marine Engineering and Technology
• /
• v.34 no.7
• /
• pp.981-990
• /
• 2010

The main objective of this study is to investigate the fluid flow and the heat transfer characteristics of nanofluids using multi-phase thermal LBM and to realize theenhancement of heat transfer characteristics considered in the Brownian motion. In multi-phase, fluid component($H_2O$) is driven by Boussinesq approximation, and nanoparticles component by the external force gravity and buoyancy. The effect of Brownian motion as a random movement is modified to the internal velocity of nanoparticles(Cu). Simultaneously, the particles of both the phases assume the local equilibrium temperature after each collision. It has been observed that when simulating $H_2O$-Cu nanoparticles, the heat transfer is the highest, at the particle volume fraction 0.5% of the particle diameter 10 nm. The average Nusselt number is increased approximately by 33% at the particle volume fraction 0.5% of the particle diameter 10 nm when compared with pure water.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.23 no.6 s.165
• /
• pp.968-978
• /
• 1999

The effective elastic properties of materials containing spherical inclusions were calculated by the elastic wave scattering theory. In the formulation additional scattering fields by the presence of random multiple scatterers that affects the effective properties were found by the single scattering approximation. In calculating the scattering fields the ensemble average on the displacements and strains inside the scatterer was found from the static approximation at long wavelength limit. The displacements were assumed to be equal to the incident field, while the strains were calculated by Eshelby's equivalent inclusion principle on the single inclusion problem. Four different models were considered and they reflected different degrees of multiple scattering effects based on the approximation introduced in the process of embedding the inclusion in the matrix. The expressions for the effective elastic constants were given in each model, and their relations to the results obtained from other scattering theory and elasticity theory were discussed. The theoretical predictions were compared with experimental results on the epoxy matrix composites containing tungsten particles of different sizes and volume fractions

• Journal of the Korean Society of Propulsion Engineers
• /
• v.19 no.4
• /
• pp.24-36
• /
• 2015

In the present study, numerical predictions of infrared spectrum of rocket plume with considering effect of particles based on approximation theories were performed by using a line-by-line radiation model with radiation databases. The high-resolution radiation databases were used to predict thermal emission spectra of gas molecules within the rocket plume regime. The particles were modeled as soot particles by using 1st term approximation of Mie theory and Rayleigh approximation. The reliability of modeled effect of soot particles using the two approximation theories was verified, and the spectral radiance of rocket plume was predicted based on the verification. The results were improved in the short wavelength range by considering the effect of soot particles.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.17 no.2
• /
• pp.67-85
• /
• 2013

Reproducing Polynomial Particle Method (RPPM) is one of meshless methods that use meshes minimally or do not use meshes at all. In this paper, the RPPM is employed for free vibration analysis of shear-deformable plates of the first order shear deformation model (FSDT), called Reissner-Mindlin plate. For numerical implementation, we use flat-top partition of unity functions, introduced by Oh et al, and patchwise RPPM in which approximation functions have high order polynomial reproducing property and the Kronecker delta property. Also, we demonstrate that our method is highly effective than other existing results for various aspect ratios and boundary conditions.

• Journal of the Korean Society of Manufacturing Process Engineers
• /
• v.8 no.3
• /
• pp.21-26
• /
• 2009

The finite element analysis of metal forming processes often fails because of severe mesh distortion at large deformation. As the concept of meshless methods, only nodal point data are used for modeling and solving. As the main feature of these methods, the domain of the problem is represented by a set of nodes, and a finite element mesh is unnecessary. This computational methods reduces time-consuming model generation and refinement effort. It provides a higher rate of convergence than the conventional finite element methods. The displacement shape functions are constructed by the reproducing kernel approximation that satisfies consistency conditions. In this research, A meshless method approach based on the reproducing kernel particle method (RKPM) is applied with metal forming analysis. Numerical examples are analyzed to verify the performance of meshless method for metal forming analysis.

• Structural Engineering and Mechanics
• /
• v.31 no.4
• /
• pp.407-421
• /
• 2009

The present paper addresses the optimal design of composite laminates with the aim of minimizing free-edge delamination stresses. A technique involving the application of particle swarm optimization (PSO) integrated with FEM was developed for the optimization. Optimization was also conducted with the zero-order method (ZOM) included in ANSYS. The semi-analytical method, which provides an approximation of the interlaminar normal stress of laminates under in-plane load, was used to partially validate the optimization results. It was found that optimal results based on ZOM are sensitive to the starting design points, and an unsuitable initial design set will lead to a result far from global solution. By contrast, the proposed method can find the global optimal solution regardless of initial designs, and the solutions were better than those obtained by ZOM in all the cases investigated.