• Structural Engineering and Mechanics
• /
• v.39 no.5
• /
• pp.669-682
• /
• 2011

A Chebyshev spectral method (CSM) for the dynamic analysis of non-uniform Timoshenko beams under various boundary conditions and concentrated masses at their ends is proposed. The matrix-based Chebyshev spectral approach was used to construct the spectral differentiation matrix of the governing differential operator and its boundary conditions. A matrix condensation approach is crucially presented to impose boundary conditions involving the homogeneous Cauchy conditions and boundary conditions containing eigenvalues. By taking advantage of the standard powerful algorithms for solving matrix eigenvalue and generalized eigenvalue problems that are embodied in the MATLAB commands, chebfun and eigs, the modal parameters of non-uniform Timoshenko beams under various boundary conditions can be obtained from the eigensolutions of the corresponding linear differential operators. Some numerical examples are presented to compare the results herein with those obtained elsewhere, and to illustrate the accuracy and effectiveness of this method.

• Journal of Aerospace System Engineering
• /
• v.14 no.3
• /
• pp.17-23
• /
• 2020

In this paper, the efficient periodic unsteady flow analysis is developed by using a Chebyshev collocation operator applied to the time differential term of the governing equations. The partial implicit time integration method was also applied in the governing equation for a fluid, which means flux terms were implicitly processed for a time integration and the time derivative terms were applied explicitly in the form of the source term by applying the Chebyshev collocation operator. To verify this method, we applied the 1D unsteady Burgers equation and the 2D oscillating airfoil. The results were compared with the existing unsteady flow frequency analysis technique, the Harmonic Balance Method, and the experimental data. The Chebyshev collocation operator can manage time derivatives for periodic and non-periodic problems, so it can be applied to non-periodic problems later.

• Journal of Mechanical Science and Technology
• /
• v.19 no.9
• /
• pp.1753-1760
• /
• 2005

The pseudo spectral method is applied to the free vibration analysis of double-span Timoshenko beams. The analysis is based on the Chebyshev polynomials. Each section of the double-span beam has its own basis functions, and the continuity conditions at the intermediate support as well as the boundary conditions are treated separately as the constraints of the basis functions. Natural frequencies are provided for different thickness-to-length ratios and for different span ratios, which agree with those of Euler-Bernoulli beams when the thickness-to-length ratio is small but deviate considerably as the thickness-to-length ratio grows larger.

• Journal of Mechanical Science and Technology
• /
• v.17 no.10
• /
• pp.1458-1465
• /
• 2003

An eigenvalue analysis of the circular Mindlin plates with free boundary conditions is presented. The analysis is based on the Chebyshev-Fourier pseudospectral method. Even though the eigenvalues of lower vibration modes tend to convergence more slowly than those of higher vibration modes, the eigenvalues converge for sufficiently fine pseudospectral grid resolutions. The eigenvalues of the axisymmetric modes are computed separately. Numerical results are provided for different grid resolutions and for different thickness-to-radius ratios.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.21 no.3
• /
• pp.109-124
• /
• 2017

The numerical solution of the Stokes equation with discontinuous viscosity and singular force term is challenging, due to the discontinuity of pressure, non-smoothness of velocity, and coupled discontinuities along interface.In this paper, we give an efficient algorithm to solve this problem by employing spectral Legendre and Chebyshev approximations.First, we present the algorithm for a problem defined in rectangular domain with straight line interface. Then it is generalized to a domain with smooth curve boundary and interface by employing spectral element method. Numerical experiments demonstrate the accuracy and efficiency of our algorithm and its spectral convergence.

• Nuclear Engineering and Technology
• /
• v.55 no.7
• /
• pp.2516-2533
• /
• 2023

Several practical methods for accelerating the depletion calculation in a GPU-based Monte Carlo (MC) code PRAGMA are presented including the multilevel spectral collapse method and the vectorized Chebyshev rational approximation method (CRAM). Since the generation of microscopic reaction rates for each nuclide needed for the construction of the depletion matrix of the Bateman equation requires either enormous memory access or tremendous physical memory, both of which are quite burdensome on GPUs, a new method called multilevel spectral collapse is proposed which combines two types of spectra to generate microscopic reaction rates: an ultrafine spectrum for an entire fuel pin and coarser spectra for each depletion region. Errors in reaction rates introduced by this method are mitigated by a hybrid usage of direct online reaction rate tallies for several important fissile nuclides. The linear system to appear in the solution process adopting the CRAM is solved by the Gauss-Seidel method which can be easily vectorized on GPUs. With the accelerated depletion methods, only about 10% of MC calculation time is consumed for depletion, so an accurate full core cycle depletion calculation for a commercial power reactor (BEAVRS) can be done in 16 h with 24 consumer-grade GPUs.

• Bulletin of the Korean Mathematical Society
• /
• v.53 no.2
• /
• pp.423-440
• /
• 2016

The bubble stabilization technique of Chebyshev-Legendre high-order element methods for one dimensional advection-diffusion equation is analyzed for the proposed scheme by Canuto and Puppo in [8]. We also analyze the finite element lower-order preconditioner for the proposed stabilized linear system. Further, the numerical results are provided to support the developed theories for the convergence and preconditioning.

• Kyungpook Mathematical Journal
• /
• v.61 no.4
• /
• pp.831-843
• /
• 2021

In this paper we look at the three dimensional stagnation point flow problem over a rough rotating disk. We study the theoretical behaviour of the stagnation point flow, or forced flow, in the presence of a slip factor in which convective instability stationary modes appear. We make a numerical investigation of the effects of slip on the behaviour of the flow components of the stagnation point flow where the disk is rough. We provide, for the first time in the literature, a complete convective instability analysis and an energy analysis. Suitable similarity transformations are used to reduce the Navier-Stokes equations and the continuity equation into a system of highly non-linear coupled ordinary differential equations, and these are solved numerically subject to suitable boundary conditions using the bvp4c function of MATLAB. The convective instability analysis and the energy analysis are performed using the Chebyshev spectral method in order to obtain the neutral curves and the energy bars. We observe that the roughness of the disk has a destabilising effect on both Type-I and Type-II instability modes. The results obtained will be prominently treated as benchmarks for our future studies on stagnation flow.

• Structural Engineering and Mechanics
• /
• v.68 no.2
• /
• pp.171-182
• /
• 2018

For the first time, nonlocal damped vibration and buckling analyses of arbitrary tapered bidirectional functionally graded solid circular nano-plate (BDFGSCNP) are presented by employing modified spectral Ritz method. The energy method based on Love-Kirchhoff plate theory assumptions is applied to derive neutral equilibrium equation. The Eringen's nonlocal continuum theory is taken into account to capture small-scale effects. The characteristic equations and corresponding first mode shapes are calculated by using a novel modified basis in spectral Ritz method. The modified basis is in terms of orthogonal shifted Chebyshev polynomials of the first kind to avoid employing adhesive functions in the spectral Ritz method. The fast convergence and compatibility with various conditions are advantages of the modified spectral Ritz method. A more accurate multivariable function is used to model two-directional variations of elasticity modulus and mass density. The effects of nanoscale, in-plane pre-load, distributed dashpot, arbitrary tapering, pinned and clamped boundary conditions on natural frequencies and buckling loads are investigated. Observing an excellent agreement between results of current work and outcomes of previously published works in literature, indicates the results' accuracy in current work.

• The Journal of Korean Institute of Communications and Information Sciences
• /
• v.19 no.2
• /
• pp.331-337
• /
• 1994

The scatting problem by E-polarized plane wave with obique incidence on a tapered resistive strip grating with zero resistivity(perfectly conducting) at strip-edges is analyzed by the method of moments in the spectral domain. Then the induced surface current density on the strip is expanded in a series of Chebyshev polynomials of the second kind. The expasion coefficients are calculated numerically in the spectral domain, the numerical results of the geometric-optical reflection coefficient for the tapered resistivity in this paper are compared with those for the existing uniform resistivity. And the position of sharp variation points in the magnitude of the geometric-optical reflection coefficient can be moved by varying the incident angle and the strip spacing, It is found out that these sparp variation points are due to the transition of higher mode between the propagation mode and the evanescent mode.