• Journal of applied mathematics & informatics
• /
• v.23 no.1_2
• /
• pp.141-152
• /
• 2007

In this paper a numerical method is presented to solve singularly perturbed two points boundary value problems for second order ordinary differential equations consisting a discontinuous source term. First, in this method, an asymptotic expansion approximation of the solution of the boundary value problem is constructed using the basic ideas of a well known perturbation method WKB. Then some initial value problems and terminal value problems are constructed such that their solutions are the terms of this asymptotic expansion. These initial value problems are happened to be singularly perturbed problems and therefore fitted mesh method (Shishkin mesh) are used to solve these problems. Necessary error estimates are derived and examples provided to illustrate the method.

• Journal of applied mathematics & informatics
• /
• v.7 no.2
• /
• pp.391-408
• /
• 2000

In the curved geometries, from the solution of the classical Riemann problem in the plane, the asymptotic solutions of the compressible Euler equation are presented. The explicit formulae are derived for the third order approximation of the generalized Riemann problem form the conventional setting of a planar shock-interface interaction.

• Journal of applied mathematics & informatics
• /
• v.35 no.3_4
• /
• pp.277-302
• /
• 2017

In this paper, we have proposed a numerical method for Singularly Perturbed Boundary Value Problems (SPBVPs) of reaction-diffusion type of third order Ordinary Differential Equations (ODEs). The SPBVP is reduced into a weakly coupled system of one first order and one second order ODEs, one without the parameter and the other with the parameter ${\varepsilon}$ multiplying the highest derivative subject to suitable initial and boundary conditions, respectively. The numerical method combines boundary value technique, asymptotic expansion approximation, shooting method and finite difference scheme. The weakly coupled system is decoupled by replacing one of the unknowns by its zero-order asymptotic expansion. Finally the present numerical method is applied to the decoupled system. In order to get a numerical solution for the derivative of the solution, the domain is divided into three regions namely two inner regions and one outer region. The Shooting method is applied to two inner regions whereas for the outer region, standard finite difference (FD) scheme is applied. Necessary error estimates are derived for the method. Computational efficiency and accuracy are verified through numerical examples. The method is easy to implement and suitable for parallel computing. The main advantage of this method is that due to decoupling the system, the computation time is very much reduced.

• Coupled systems mechanics
• /
• v.7 no.1
• /
• pp.61-77
• /
• 2018

This paper establish a high concentration ratio approximation of linear elastic properties of materials with periodic microstructure with cubic inclusions. The approximation is derived using first few terms of power series expansion of the solution of the equivalent eigenstrain problem with a homogeneous eigenstrain approximation. Viability of the approximation at high concentration ratios is proved by comparison with a numerical solution of the homogenization problem. To this end some theoretical result of symmetry properties of the homogenization problem are given. Using these results efficient numerical computation on a reduced computational domain is presented.

• Kyungpook Mathematical Journal
• /
• v.59 no.4
• /
• pp.735-769
• /
• 2019

In this paper, we consider the Robin-Dirichlet problem for a nonlinear wave equation of Kirchhoff-Carrier type with a viscoelastic term. Using the Faedo-Galerkin method and the linearization method for nonlinear terms, the existence and uniqueness of a weak solution are proved. An asymptotic expansion of high order in a small parameter of a weak solution is also discussed.

• Journal of applied mathematics & informatics
• /
• v.33 no.5_6
• /
• pp.485-502
• /
• 2015

A computational method for solving singularly perturbed boundary value problem of differential equation with shift arguments of mixed type is presented. When shift arguments are sufficiently small (o(ε)), most of the existing method in the literature used Taylor's expansion to approximate the shift term. This procedure may lead to a bad approximation when the delay argument is of O(ε). The main idea for this work is to deal with constant shift arguments, which are independent of ε. In the present method, we construct the formally asymptotic solution of the problem using the method of composite expansion. The reduced problem is solved numerically by using operator compact implicit method, and the second problem is solved analytically. Error estimate is derived by using the maximum norm. Numerical examples are provided to support the theoretical results and to show the efficiency of the proposed 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.

• Communications for Statistical Applications and Methods
• /
• v.7 no.1
• /
• pp.55-64
• /
• 2000

We propose a criterion for testing homogeneity of matrix intraclass covariance matrices of K multivariate normal populations, It is based on a variable transformation intended to propose and develop a likelihood ratio criterion that makes use of properties of eigen structures of the matrix intraclass covariance matrices. The criterion then leads to a simple test that uses an asymptotic distribution obtained from Box's (1949) theorem for the general asymptotic expansion of random variables.

• The Journal of the Acoustical Society of Korea
• /
• v.33 no.4
• /
• pp.232-237
• /
• 2014

A passive sonar of warship is composed of several directional or omni-directional sensors. In order to model the acoustic signal received into a warship sonar, the wave propagation modeling is usually required from arbitrary noise source to all sensors equipped to the sonar. However, the full calculation for all sensors is time-consuming and the performance of sonar simulator deteriorates. In this study, we suggest an asymptotic method to estimate the sonar signal arrived to sensors adjacent to the reference sensor, where it is assumed that all information of eigenrays is known. This method is developed using Taylor series for the time delay of eigenray and similar to Fraunhofer and Fresnel approximation for sonar aperture. To validate the proposed method, some numerical experiments are performed for the passive sonar. The approximation when the second-order term is kept is vastly superior. In addition, the error criterion for each approximation is provided with a practical example.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.13 no.1
• /
• pp.51-61
• /
• 2000

This paper discusses the homogenization method to determine effective average elastic constants of a linear structure by considering its microstructure. A detailed description on the homogenization method is given for the linear elastic material and then the finite element approximation is performed for an investigation of elastic properties. An asymptotic expansion is carried out in the cross-section area, or in the unit cell. Two and three lay-up structures made up of individual isotropic constituents are chosen for numerical examples to check discrepancies between results generated by this theoretical development and the conventional approach. Asymptotic characteristics of the process in extracting the stiffness of structure locally formed by spatial repetitions yield underestimated values of stiffness. These discrepancies are detected by the asymptotic corrective term which is ascribed to considerations of microscopic perturbations and proved in the finite element formulation. The asymptotic analysis is the more reasonable in analysing the composite material, rather than the conventional approach to calculate the macroscopic average for elastic properties.