• Structural Engineering and Mechanics
• /
• v.11 no.4
• /
• pp.373-392
• /
• 2001

The main idea behind the paper is to present two alternative methods of homogenization of the heat conduction problem in composite materials, where the heat conductivity coefficients are assumed to be random variables. These two methods are the Monte-Carlo simulation (MCS) technique and the second order perturbation second probabilistic moment method, with its computational implementation known as the Stochastic Finite Element Method (SFEM). From the mathematical point of view, the deterministic homogenization method, being extended to probabilistic spaces, is based on the effective modules approach. Numerical results obtained in the paper allow to compare MCS against the SFEM and, on the other hand, to verify the sensitivity of effective heat conductivity probabilistic moments to the reinforcement ratio. These computational studies are provided in the range of up to fourth order probabilistic moments of effective conductivity coefficient and compared with probabilistic characteristics of the Voigt-Reuss bounds.

• The Bulletin of The Korean Astronomical Society
• /
• v.44 no.1
• /
• pp.46.1-46.1
• /
• 2019

We present an accurate and efficient method to calculate the gravitational potential of an isolated system in three-dimensional Cartesian and cylindrical coordinates subject to vacuum (open) boundary conditions. Our method consists of two parts: an interior solver and a boundary solver. The interior solver adopts an eigenfunction expansion method together with a tridiagonal matrix solver to solve the Poisson equation subject to the zero boundary condition. The boundary solver employs James's method to calculate the boundary potential due to the screening charges required to keep the zero boundary condition for the interior solver. A full computation of gravitational potential requires running the interior solver twice and the boundary solver once. We develop a method to compute the discrete Green's function in cylindrical coordinates, which is an integral part of the James algorithm to maintain second-order accuracy. We implement our method in the {\tt Athena++} magnetohydrodynamics code, and perform various tests to check that our solver is second-order accurate and exhibits good parallel performance.

• Bulletin of the Korean Chemical Society
• /
• v.7 no.4
• /
• pp.262-266
• /
• 1986

A multi-reference many-body perturbation theory (MRMBPT) method is critically tested in second order by comparing with the corresponding configuration interaction (CI) calculations. Excitation energies of the four-valence-electron states of transition metal atoms and ions are used for the comparison. The agreement between the second order MRMBPT and CI calculations is very reasonable, confirming the reliability of the second order MRMBPT method. The reliability of calculations with the present second order MRMBPT method was only been inferred empirically in the past since most results have been gauged by the agreement with experiment and/or with other MRMBPT calculations based upon different sets of orbitals and configuration spaces. The present MRMBPT method appears to be an efficient ab initio multi-reference method for the calculation of electron correlation effects in atoms and molecules, and it is shown how MRMBPT can be used to estimate core-core and core-valence correlation effects which are often omitted in CI calculations because too many configurations and correlating electrons are involved.

• Journal of the Korean Institute of Telematics and Electronics
• /
• v.23 no.5
• /
• pp.647-652
• /
• 1986

The controller tuning methods proposed by Yuwana-Seborg and Suh utilizes Pade first order approximation for the delay terms in the closed loop transfer function. In this paper, the use of a Pade second-order approximation method is investigated. The simulation results show that the new method is superior to pervious approaches such as Ziegler-Nichols and Cohen-Coon methods.

• Journal of the Korean Society for Precision Engineering
• /
• v.4 no.3
• /
• pp.35-43
• /
• 1987

A new analytic method for the identification of process control system is presented. A second and a third order transfer function are considered as the estimated model function. For the second order transfer function, the new method is compared with the exisiting ones, simulation results show that the new method is superior to the existing ones. And also, In case of the third order transfer function which is difficult to analyze mathematically, system identification is tried.

• The Pure and Applied Mathematics
• /
• v.17 no.4
• /
• pp.289-298
• /
• 2010

In this paper, we propose a second-order prediction/correction (SPC) domain decomposition method for solving one dimensional linear hyperbolic partial differential equation $u_{tt}+a(x,t)u_t+b(x,t)u=c(x,t)u_{xx}+{\int}(x,t)$. The method can be applied to variable coefficients problems and singular problems. Unconditional stability and error analysis of the method have been carried out. Numerical results support stability and efficiency of the method.

• Journal of applied mathematics & informatics
• /
• v.26 no.5_6
• /
• pp.1057-1069
• /
• 2008

We consider singularly perturbed Boundary Value Problems (BVPs) for third and fourth order Ordinary Differential Equations(ODEs) of convection-diffusion type with discontinuous source term and a small positive parameter multiplying the highest derivative. Because of the type of Boundary Conditions(BCs) imposed on these equations these problems can be transformed into weakly coupled systems. In this system, the first equation does not have the small parameter but the second contains it. In this paper a computational method named as 'An asymptotic finite element method' for solving these systems is presented. In this method we first find an zero order asymptotic approximation to the solution and then the system is decoupled by replacing the first component of the solution by this approximation in the second equation. Then the second equation is independently solved by a fitted mesh Finite Element Method (FEM). Numerical experiments support our theoritical results.

• Journal of Ocean Engineering and Technology
• /
• v.12 no.2 s.28
• /
• pp.33-42
• /
• 1998

In the presence of incident waves with different frequencies, the second order sum and difference frequency waves due to the nonlinearity of the incident waves come into existence. Although the magnitudes of the forces produced on a Tension Leg Platform(TLP) by these nonlinear waves are small, they act on the TLP at sum and difference frequencies away from those of the incident waves. So, the second order sum and difference frequency wave loads produced close to the natural frequencies of TLPs often give greater contributions to high and low frequency resonant responses. The second order wave exciting forces and moments have been obtained by the method based on direct integration of pressure acting on the submerged surface of a TLP. The components of the second order forces which depend on first order quantities have been evaluated using the three dimensional source distribution method. The numerical results of time domain analysis for the nonlinear wave exciting forces in regular waves are compared with the numerical ones of frequency domain analysis. The results of comparison confirmed the validity of the proposed approach.

• Journal of applied mathematics & informatics
• /
• v.36 no.1_2
• /
• pp.135-154
• /
• 2018

In this paper, an almost second order overlapping Schwarz method for singularly perturbed third order convection-diffusion type problem is constructed. The method splits the original domain into two overlapping subdomains. A hybrid difference scheme is proposed in which on the boundary layer region we use the combination of classical finite difference scheme and central finite difference scheme on a uniform mesh while on the non-layer region we use the midpoint difference scheme on a uniform mesh. It is shown that the numerical approximations which converge in the maximum norm to the exact solution. We proved that, when appropriate subdomains are used, the method produces convergence of second order. Furthermore, it is shown that, two iterations are sufficient to achieve the expected accuracy. Numerical examples are presented to support the theoretical results. The main advantages of this method used with the proposed scheme are it reduce iteration counts very much and easily identifies in which iteration the Schwarz iterate terminates.

• Journal of Korean Society of Steel Construction
• /
• v.21 no.6
• /
• pp.627-636
• /
• 2009

The purpose of this study was to evaluate the validity of the second-order analysis using the load amplification factor suggested by design codes. For this purpose, the first-order analysis with the B1 and B2 factors suggested by KBC 2005 and the direct analysis with the load amplification factor suggested by KBC 2009 (draft) were performed for three-story -one-bay and five-story-three-bay unbraced steel frames. The results of the analyses were compared with the results of the second-order inelastic analysis to evaluate the validity of the suggested methods. The main parameters of the analysis were the scale of the frame, the axial load ratio of the column, and the methods of analysis. The research results showedthat the method suggested by KBC 2005 does not properly consider the second-order effect under the high axial load ratio, but the direct analysis method suggested by KBC 2009 (draft) properly estimates the second-order effect without any serious problem.