• Journal of applied mathematics & informatics
• /
• v.25 no.1_2
• /
• pp.435-443
• /
• 2007

Computing the interior spectrum of large sparse generalized eigenvalue problems $Ax\;=\;{\lambda}Bx$, where A and b are large sparse and SPD(Symmetric Positive Definite), is often required in areas such as structural mechanics and quantum chemistry, to name a few. Recently, CG-type methods have been found useful and hence, very amenable to parallel computation for very large problems. Also, as in the case of linear systems proper choice of preconditioning is known to accelerate the rate of convergence. After the smallest eigenpair is found we use the orthogonal deflation technique to find the next m-1 eigenvalues, which is also suitable for parallelization. This offers advantages over Jacobi-Davidson methods with partial shifts, which requires re-computation of preconditioner matrx with new shifts. We consider as preconditioners Incomplete LU(ILU)(0) in two variants, ever-relaxation(SOR), and Point-symmetric SOR(SSOR). We set m to be 5. We conducted our experiments on matrices from discretizations of partial differential equations by finite difference method. The generated matrices has dimensions up to 4 million and total number of processors are 32. MPI(Message Passing Interface) library was used for interprocessor communications. Our results show that in general the Multi-Color ILU(0) gives the best performance.

• Smart Structures and Systems
• /
• v.25 no.2
• /
• pp.219-228
• /
• 2020

This work presents a modified continuum model to explore and investigate static and vibration behaviors of perforated piezoelectric NEMS structure. The perforated nanostructure is modeled as a thin perforated nanobeam element with Euler-Bernoulli kinematic assumptions. A size scale effect is considered by included a nonlocal constitutive equation of Eringen in differential form. Modifications of geometrical parameters of perforated nanobeams are presented in simplified forms. To satisfy the Maxwell's equation, the distribution of electric potential for the piezoelectric nanobeam model is assumed to be varied as a combination of a cosine and linear functions. Hamilton's principle is exploited to develop mathematical governing equations. Modified numerical finite model is adopted to solve the equation of motion and equilibrium equation. The proposed model is validated with previous respectable work. Numerical investigations are presented to illustrate effects of the number of perforated holes, perforation size, nonlocal parameter, boundary conditions, and external electric voltage on the electro-mechanical behaviors of piezoelectric nanobeams.

• Steel and Composite Structures
• /
• v.36 no.6
• /
• pp.643-656
• /
• 2020

This article presents a comprehensive static analysis of simply supported cross-ply carbon nanotubes reinforced composite (CNTRC) laminated nanobeams under various loading profiles. The nonlocal strain gradient constitutive relation is exploited to present the size-dependence of nano-scale. New higher shear deformation beam theory with hyperbolic function is proposed to satisfy the zero-shear effect at boundaries and parabolic variation through the thickness. Carbon nanotubes (CNTs), as the reinforced elements, are distributed through the beam thickness with different distribution functions, which are, uniform distribution (UD-CNTRC), V- distribution (FG-V CNTRC), O- distribution (FG-O CNTRC) and X- distribution (FG-X CNTRC). The equilibrium equations are derived, and Fourier series function are used to solve the obtained differential equation and get the response of nanobeam under uniform, linear or sinusoidal mechanical loadings. Numerical results are obtained to present influences of CNTs reinforcement patterns, composite laminate structure, nonlocal parameter, length scale parameter, geometric parameters on center deflection ad stresses of CNTRC laminated nanobeams. The proposed model is effective in analysis and design of composite structure ranging from macro-scale to nano-scale.

• Smart Structures and Systems
• /
• v.18 no.6
• /
• pp.1125-1143
• /
• 2016

Forced vibration analysis of a simple supported viscoelastic nanobeam is studied based on modified couple stress theory (MCST). The nanobeam is excited by a transverse triangular force impulse modulated by a harmonic motion. The elastic medium is considered as Winkler-Pasternak elastic foundation.The damping effect is considered by using the Kelvin-Voigt viscoelastic model. The inclusion of an additional material parameter enables the new beam model to capture the size effect. The new non-classical beam model reduces to the classical beam model when the length scale parameter is set to zero. The considered problem is investigated within the Timoshenko beam theory by using finite element method. The effects of the transverse shear deformation and rotary inertia are included according to the Timoshenko beam theory. The obtained system of differential equations is reduced to a linear algebraic equation system and solved in the time domain by using Newmark average acceleration method. Numerical results are presented to investigate the influences the material length scale parameter, the parameter of the elastic medium and aspect ratio on the dynamic response of the nanobeam. Also, the difference between the classical beam theory (CBT) and modified couple stress theory is investigated for forced vibration responses of nanobeams.

• Geomechanics and Engineering
• /
• v.13 no.4
• /
• pp.585-600
• /
• 2017

This study deals with simple solutions for a spherical or circular opening excavated in elastic-brittle plastic rock mass compatible with a linear Mohr-Coulomb (M-C) or a nonlinear Hoek-Brown (H-B) yield criterion. Based on total strain approach, the closed-form solutions of stresses and displacement are derived simultaneously for circular and spherical openings using original H-B and M-C yield criteria. Two simple numerical procedures are proposed for the solution of generalized H-B and M-C yield criteria. Based on incremental approach, the similarity solution is derived for circular and spherical openings using generalized H-B and M-C yield criteria. The classical Runge-Kutta method is used to integrate the first-order ordinary differential equations. Using three data sets for M-C and H-B models, the results of the radial displacements, the spreading of the plastic radius with decreasing pressure, and the radial and circumferential stresses in the plastic region are compared. Excellent agreement among the solutions is obtained for all cases of spherical and circular openings. The importance of the use of proper initial values in the similarity solution is discussed.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.18 no.9
• /
• pp.714-721
• /
• 2006

In evaporative cooling applications, the evaporation water is supplied usually sufficiently larger than the amount evaporated to enlarge contact surface between the water and the air. Especially in indirect evaporative coolers, however, if the evaporation water flow rate is excessively large, the evaporative cooling effect is not used for heat absorption from the hot fluid but spent to the sensible cooling of the evaporation water itself. This would result in a decrease in the cooling performance of the indirect evaporative cooler. In this study, the effects of the evaporation water flow rate on the cooling performance are investigated theoretically. The cooling process in an indirect evaporative cooler is modeled into a set of linear differential equations and solved to obtain the exact solutions to the temperatures of the hot fluid, the moist air, and the evaporation water. Based on the exact solutions, it is analyzed how much the cooling performance is affected by the evaporation water flow rate. The results show that the decrease in the cooling effectiveness is substantial even for a small flow rate of the evaporation water and the relative decrease is more serious for a high-performance evaporative cooler.

• Proceedings of the KSME Conference
• /
• 2004.11a
• /
• pp.1152-1157
• /
• 2004

Temporal behavior of the laser induced incandescence (LII) signal is often used for soot particle sizing, which is possible because the cooling behavior of a laser heated particle is dependent on the particle size. In present study, LII signals of soot particles are modeled using two non-linear coupled differential equations deduced from the energy- and mass-balance of the process. The objective of this study is to see the effects of particle size, laser fluence on soot temperature characteristics and cooling behavior. Together with this, we focus on validating our simulation code by comparing with other previous results. Results of normalized LII signals obtained from various laser fluence conditions showed a good agreement with that of Dalzell and Sarofim's. It could be found that small particles cool faster at a constant laser fluence. And it also could be observed that vaporization is dominant process of heat loss during first 100ns after laser pulse, then heat conduction played most important role while thermal radiation had little influence all the time.

• Research in Mathematical Education
• /
• v.12 no.4
• /
• pp.283-291
• /
• 2008

The extension to a wide population of secondary education in many advanced countries seems to have led to a weakening of the mathematics curriculum. In response, many students have been classified as "gifted" so that they can access a stronger program. Apart from the difficulties that might arise in actually determining which students are gifted (Is it always clear what the term means?), there are dangers inherent in programs that might be devised even for those that are truly talented. Sometimes students are moved ahead to more advanced mathematics. Elementary students might be taught algebra or even subjects like trigonometry and vectors, and secondary students might be taught calculus, differential equations and linear algebra. It is my experience over thirty-five years of contact with bright students that acceleration to higher level mathematics is often not a good idea. In this paper, I will articulate some of the factors that have led me to this opinion and suggest alternatives. First, I would like to emphasize that in matters of education, almost every statement that can be made to admit counterexamples; my opinion on acceleration is no exception. Occasionally, a young Gauss or Euler walks in the door, and one has no choice but to offer the maximum encouragement and allow the student to go to the limit of his capabilities. A young genius can demonstrate an incredible amount of mathematical insight, maturity and mastery of technique. A classical example is probably the teen-age Euler, who in the 1720s was allowed regular audiences with Jean Bernoulli, the foremost mathematician of his day.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2008.11a
• /
• pp.256-263
• /
• 2008

The assumed modes method is developed to derive a set of linear differential equations describing the motion of a flexible wind turbine blade and to propose an approach to investigate the forced responses result from various wind excitations. In this work, we have adopted Euler beam theory and considered that the root of the blade is clamped at the rigid hub. And the aerodynamic parameters and forces are determined based on Blade Element Momentum (BEM) theory and quasi-steady airfoil aerodynamics. Numerical calculations show that this method gives good results and it can be used fur modeling and the forced vibration analysis including the coupling effect of wind-turbine blades, as well as turbo-machinery blades, aircraft propellers or helicopter rotor blades which may be considered as straight non-uniform beams with built-in pre-twist.

• 한국전산유체공학회:학술대회논문집
• /
• 2005.10a
• /
• pp.293-298
• /
• 2005

We calculate the coordinates of an axisymmetric nozzle with a central body. This nozzle ensures a transonic flow with a plane sound surface, which is orthogonal to the symmetry axis and has a wall kink at the sonic point, The Chaplygin transformation in the subsonic part of the flow leads the Dirichlet problem for a system of nonlinear equations. The definition domain of the solution in the velocity-hodograph plane is taken as a rectangle. This enables one to obtain the nozzle with a monotonic distribution of velocity along its subsonic part. In the nonlinear differential equation, the linear Chaplygin operator for plane flows is separated, which allows the iterative calculation of the solution. The supersonic part of the nozzle is calculated under the assumption that the flow at the nozzle exit is uniform and parallel to the symmetry axis; i.e., the supersonic jet outflows to the submerged space with the same pressure. The calculation is performed by the characteristic method. The exact solution of Tricomi equation for near-sonic flows with the straight sonic line is used to 'move away' the sound plane. The velocity distribution alone the supersonic part of the nozzle is also monotonic, which ensures the absence of the boundary-layer separation and, therefore, the adequacy of the ideal-gas model. calculations show that the flow in the supersonic part of the nozzle is continuous (compression shocks are absent)