• Proceedings of the Korean Society of Marine Engineers Conference
• /
• 2001.05a
• /
• pp.52-57
• /
• 2001

Laminar two-dimensional time-dependent flow past a circular cylinder is numerically investigated using direct numerical simulation for the low Reynolds number (Re=164∼280). The higher-order finite difference scheme is employed for the spatial distributions along with the second order Adams-Bashforth and the first order backward-Euler time integration. The convection term is applied by the 7th order up wind scheme and the pressure and viscosity terms are applied by the 4th order central difference. The grid system makes use of the regular grid system and it is generated by an equation. The calculated results of drag coefficients, lift coefficients, pressure distributions, and vorticity contours and other information are compared with experimental and numerical ones. These results obtained by the present DNS show good agreement with the previous studies.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 2001.05a
• /
• pp.205-208
• /
• 2001

Localized shear band is investigated through the analysis of one-dimensional model for simple shearing deformation of thermally rate dependent material. Generally mesh size or interval of nodes play an important role in determining the overall flow behavior of the material. In order to observe these size effects we adapted non-local theory by including higher order strain gradients of the equivalent strain into the constitutive equation for the flow stress. for the ease of convergence and numerical stability the inplicit finite difference scheme is employed.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.27 no.4
• /
• pp.585-592
• /
• 2003

Among various sensitivity evaluation techniques, semi-analytical method(SAM) is quite popular since this method is more advantageous than analytical method(AM) and global finite difference method(FDM). However, SAM reveals severe inaccuracy problem when relatively large rigid body motions are identified fur individual elements. Such errors result from the numerical differentiation of the pseudo load vector calculated by the finite difference scheme. In the present study, an iterative method combined with mode decomposition technique is proposed to compute reliable semi-analytical design sensitivities. The improvement of design sensitivities corresponding to the rigid body mode is evaluated by exact differentiation of the rigid body modes and the error of SAM caused by numerical difference scheme is alleviated by using a Von Neumann series approximation considering the higher order terms for the sensitivity derivatives.

• Journal of electromagnetic engineering and science
• /
• v.13 no.3
• /
• pp.158-164
• /
• 2013

This paper expands a previously proposed optimized higher order (2, 4) finite-difference time-domain scheme (H(2, 4) scheme) for use with lossy material. A low dispersion error is obtained by introducing a weighting factor and two scaling factors. The weighting factor creates isotropic dispersion, and the two scaling factors dramatically reduce the numerical dispersion error at an operating frequency. In addition, the results confirm that the proposed scheme performs better than the H(2, 4) scheme for wideband analysis. Lastly, the validity of the proposed scheme is verified by calculating a scattering problem of a lossy circular dielectric cylinder.

• Ocean Systems Engineering
• /
• v.10 no.2
• /
• pp.201-226
• /
• 2020

A systematic numerical comparative study of the performance of semicircular and rectangular submerged breakwaters interacting with solitary waves is the basis of this paper. To accomplish this task, Nwogu's extended Boussinesq model equations are employed to simulate the interaction of the wave with breakwaters. The finite difference technique has been used to discretize the spatial terms while a fourth-order predictor-corrector method is employed for time discretization in our numerical model. The proposed computational scheme uses a staggered-grid system where the first-order spatial derivatives have been discretized with fourth-order accuracy. For validation purposes, five test cases are considered and numerical results have been successfully compared with the existing analytical and experimental results. The performances of the rectangular and semicircular breakwaters have been examined in terms of the wave reflection, transmission, and dissipation coefficients (RTD coefficients) denoted by K_R, K_T, K_D. The latter coefficient K_D emerges due to the non-energy conserving K_R and K_T. Our computational results and graphical illustrations show that the rectangular breakwater has higher reflection coefficients than semicircular breakwater for a fixed crest height, but as the wave height increases, the two reflection coefficients approach each other. un the other hand, the rectangular breakwater has larger dissipation coefficients compared to that of the semicircular breakwater and the difference between them increases as the height of the crest increases. However, the transmission coefficient for the semicircular breakwater is greater than that of the rectangular breakwater and the difference in their transmission coefficients increases with the crest height. Quantitatively, for rectangular breakwaters the reflection coefficients K_R are 5-15% higher while the diffusion coefficients K_D are 3-23% higher than that for the semicircular breakwaters, respectively. The transmission coefficients K_T for rectangular breakwater shows the better performance up to 2.47% than that for the semicircular breakwaters. Based on our computational results, one may conclude that the rectangular breakwater has a better overall performance than the semicircular breakwater. Although the model equations are non-dissipative, the non-energy conserving transmission and reflection coefficients due to wave-breakwater interactions lead to dissipation type contribution.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.33 no.5
• /
• pp.279-286
• /
• 2020

In this study, we investigate the accuracy of higher order derivatives in the moving least square (MLS) difference method. An interpolation function is constructed by employing a Taylor series expansion via MLS approximation. The function is then applied to the mixed variational theorem in which the displacement and stress resultants are treated as independent variables. The higher order derivatives are evaluated by solving simply supported beams and cantilevers. The results are compared with the analytical solutions in terms of the order of polynomials, support size of the weighting function, and number of nodes. The accuracy of the higher order derivatives improves with the employment of the mean value theorem, especially for very high-order derivatives (e.g., above fourth-order derivatives), which are important in a classical asymptotic analysis.

• 한국전산유체공학회:학술대회논문집
• /
• 2001.10a
• /
• pp.84-90
• /
• 2001

Secondary instability of flow past a circular cylinder is examined using direct numerical simulation at Reynolds number 220 and 250. The higher-order finite difference scheme is employed for the spatial distributions along with the second order Adams-Bashforth and the first order backward-Euler time integration. In x-y plane, the convection term is applied by the 5th order upwind scheme, and the pressure and viscosity terms are applied by the 4th order central difference. In spanwise, Navier-Stokes equation is distributed using Spectral Method. The critical Reynolds number for this instability is found to be about Re=190. The secondary instability leads re three-dimensionality with a spanwise wavelength about 4 cylinder diameters at onset (A-mode). Results of three-dimensional effect in wake of a circular cylinder are represented with spanwise and streamwise vorticity contours as Reynolds numbers.

• Proceedings of the KSME Conference
• /
• 2001.11b
• /
• pp.570-577
• /
• 2001

Three-dimensional time-dependent flow past a circular cylinder is numerically investigated using direct numerical simulation for Reynolds number 280 and 300. The higher-order finite difference scheme is employed for the spatial distributions along with the second order Adams-Bashforth and the first order backward-Euler time integration. In x-y plane, the convection term is applied by the 5th order upwind scheme and the pressure and viscosity terms are applied by the 4th order central difference. And in spanwise, Navier-Stokes equation is distributed using of Spectral Method. At Reynolds number 259 the two-dimensional wake becomes linearly unstable to a second branch of modes with wavelength about 1.0 diameters at onset (B-mode). Present results of three-dimensional effects of in wake of a circular cylinder is represented with spanwise and streamwise vorticity contours as Reynolds numbers.

• Structural Engineering and Mechanics
• /
• v.48 no.5
• /
• pp.711-736
• /
• 2013

In this paper new implicit time integration called N-IHOA is presented for dynamic analysis of high damping systems. Here, current displacement and velocity are assumed to be functions of the velocities and accelerations of several previous time steps, respectively. This definition causes that only one set of weighted factors is calculated from the Taylor series expansion which leads to a simple approach and reduce the computational efforts. Moreover a comprehensive study on stability of the proposed method i.e., N-IHOA compared with IHOA integration which is performed based on amplification matrices proves the ability of the N-IHOA in high damping vibrations such as control systems. Also, wide range of numerical examples which contains single/multi degrees of freedom, damped/un-damped, free/forced vibrations from finite element/finite difference demonstrate that the accuracy and efficiency of the proposed time integration is more than the common approaches such as the IHOA, the Wilson-${\theta}$ and the Newmark-${\beta}$.

• Communications for Statistical Applications and Methods
• /
• v.21 no.3
• /
• pp.201-212
• /
• 2014

We provide a simple basic method to find bounds for higher order moments of unimodal distributions in terms of lower order moments when the random variable takes value in a given finite real interval. The bounds for moments in terms of the geometric mean of the distribution are also derived. Both continuous and discrete cases are considered. The bounds for the ratio and difference of moments are obtained. The special cases provide refinements of several well-known inequalities, such as Kantorovich inequality and Krasnosel'skii and Krein inequality.