• Structural Engineering and Mechanics
• /
• 제39권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.

• 항공우주시스템공학회지
• /
• 제14권3호
• /
• pp.17-23
• /
• 2020

본 연구는 Chebyshev collocation operator를 지배 방정식의 시간 미분항에 적용하여 비정상 유동해석을 해석할 수 있는 기법을 개발한 논문이다. 시간적분으로 유속항은 내재적으로 처리하였으며 시간 미분항은 Chebyshev collocation operator을 적용하여 원천항 형태로 외재적으로 처리하여 부분 내재적 시간적분법을 적용하였다. 본 연구의 방법을 검증하기 위해 1차원 비정상 burgers 방정식과 2차원 진동하는 airfoil에 적용하였으며 기존의 비정상 유동 주파수 해석기법과 시험 결과를 비교하여 나타내었다. Chebyshev collocation operator는 주기적인 문제와 비주기적인 문제에 대해서 시간 미분항을 처리할 수 있으므로 추후 비주기적인 문제에 적용할 예정이다.

• Journal of Mechanical Science and Technology
• /
• 제19권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
• /
• 제17권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
• /
• 제21권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
• /
• 제55권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.

• 대한수학회보
• /
• 제53권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
• /
• 제61권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
• /
• 제68권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.

• 한국통신학회논문지
• /
• 제19권2호
• /
• pp.331-337
• /
• 1994

저항띠의 양끝에서 0( /square)으로 변하는 저항율을 가진 저항띠의 적자구조에 비스듬히 입사하는 E분극 평면파에 의한 전자파 산람눙제를 과수영역에서 모멘트 법을 이용하여 해석하였따. 이때 저항띠에 유도되는 전류일도는 2종 Chebyshev 다항식의 급수로 전개하였다 전개계수들은 과수영역에서 수치계산하였고, 본 논문에서의 변하는 저항율을 갖는 경우와 기존의 균일 저항율을 갖는 경우에 대해 기하광학적 반사계수의 수치계산 결과를 비교하였다. 그리고 기하광학적 반사계수의 크기에서 급변점들이 위치는 입사각과 스트립 주기를 변화시킴으로써 이동시킬 수 있었다. 이러한 급변점들은 전파모드와 감쇠모드 사이서 고차모드가 천이될 때 발생함을 알 수 있었다.