• Communications of the Korean Mathematical Society
• /
• v.10 no.1
• /
• pp.39-48
• /
• 1995

The ADI method was introduced by Peaceman and Rachford [6] in 1955, to solve the discretized boundary value problems for elliptic and parabolic PDEs. The finite difference discretization of the model elliptic problem $$ (1) -\Delta u = f, \Omega = [0, 1] \times [0, 1] $$ $$ u = 0 on \delta \Omega $$ with 5-point centered finite difference discretization, with n +2 mesh-points in the x - direction and m + 2 points in the y direction, leads to the solution of a linear system of equations of the form $$ (2) Au = b $$ where A is a matrix of dimension $N = n \times m$. Without loss of generality and for the sake of simplicity, we will assume for the remainder of this paper that m = n, so that $N = n^2$.

• Journal of the Korean Mathematical Society
• /
• v.52 no.5
• /
• pp.1069-1096
• /
• 2015

In this article we introduce a numerical method, named Gegenbauer wavelets method, which is derived from conventional Gegenbauer polynomials, for solving fractional initial and boundary value problems. The operational matrices are derived and utilized to reduce the linear fractional differential equation to a system of algebraic equations. We perform the convergence analysis for the Gegenbauer wavelets method. We also combine Gegenbauer wavelets operational matrix method with quasilinearization technique for solving fractional nonlinear differential equation. Quasilinearization technique is used to discretize the nonlinear fractional ordinary differential equation and then the Gegenbauer wavelet method is applied to discretized fractional ordinary differential equations. In each iteration of quasilinearization technique, solution is updated by the Gegenbauer wavelet method. Numerical examples are provided to illustrate the efficiency and accuracy of the methods.

• Honam Mathematical Journal
• /
• v.33 no.2
• /
• pp.223-229
• /
• 2011

In investigation of boundary-value problems for certain partial differential equations arising in applied mathematics, we often need to study the solution of system of partial differential equations satisfied by hypergeometric functions and find explicit linearly independent solutions for the system. Here we choose the Exton functions $X_1$ and $X_2$ among his twenty functions to show how to find the linearly independent solutions of partial differential equations satisfied by these functions $X_1$ and $X_2$.

• Proceedings of the Korean Society For Composite Materials Conference
• /
• 2003.10a
• /
• pp.75-78
• /
• 2003

In order to describe the Bauschinger and transient behavior of orthotropic fiber-reinforced composites, a combined isotropic-kinematic hardening law based on the non-linear kinematic hardening rule was considered here, in particular, based on the Chaboche type law. In this modified constitutive law, the anisotropic evolution of the back-stress was properly accounted for. Also, to represent the orthotropy of composite materials, Hill's 1948 quadratic yield function and the orthotropic elasticity constitutive equations were utilized. Furthermore, the numerical formulation to update the stresses was also developed based on the incremental deformation theory for the boundary value problems. Numerical examples confirmed that the new law based on the anisotropic evolution of the back-stress complies well with the constitutive behavior of highly anisotropic materials such as fiber-reinforced composites.

• Journal of the Korean Society of Manufacturing Process Engineers
• /
• v.6 no.4
• /
• pp.81-85
• /
• 2007

A safety valve is used for protecting the pressure vessel and facilities by discharging the operating fluid into the valve from the accident when the pressure is over the designated value. The fluid is sulfurous acid and nitric acid. etc. in the semi-conductor assembly line. Thus the valve elements material must be acid resistance. Teflon, which is used generally as inner parts of a valve, tends to easily sticks to sliding surface by thermal expansion under high temperature. Some studies are performed to change teflon to another material and shape to have a better fluidity under the condition. The analysis of the thermal expansion is conducted by commercial FEM software to improve the problems. Boundary conditions were temperature and load in this study. From the analysis, the thermal expansion of stainless steel is verified to be lower than that of teflon under high temperature. Thus coupled teflon/stainless steel-made valve is applied to assembly line without danger due to thermal expansion.

• Journal of applied mathematics & informatics
• /
• v.13 no.1_2
• /
• pp.183-194
• /
• 2003

We will explain a new method for obtaining the nearly optimal domain for optimal shape design problems associated with the solution of a nonlinear wave equation. Taking into account the boundary and terminal conditions of the system, a new approach is applied to determine the optimal domain and its related optimal control function with respect to the integral performance criteria, by use of positive Radon measures. The approach, say shape-measure, consists of two steps; first for a fixed domain, the optimal control will be identified by the use of measures. This function and the optimal value of the objective function depend on the geometrical variables of the domain. In the second step, based on the results of the previous one and by applying some convenient optimization techniques, the optimal domain and its related optimal control function will be identified at the same time. The existence of the optimal solution is considered and a numerical example is also given.

• Proceedings of the KSME Conference
• /
• 2002.08a
• /
• pp.711-712
• /
• 2002

A numerical study of a hemodynamical model for the tumor angiogenesis is carried out. The tumor angiogenesis process is comprised of a sequence of events; secretion of tumor angiogenesis factor(TAF) from the solid tumor, degradation of the basement membrane of nearby blood vessels, migration and proliferation of the endothelial cells. The model takes into account the effect of TAF concentration and endothelial cell density, and their conservation equations are represented as a set of one-dimensional initial boundary value problems. These equations are discretized by using a finite difference method in which the second order schemes both in time and in space are used. The effects of the parameters contained in the model are Investigated extensively through the numerical simulation of the discretized model. The result for the typical case compares very well with the known result.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.24 no.10 s.181
• /
• pp.2469-2476
• /
• 2000

The meshing of the domain has long been the major bottleneck in performing the finite element analysis. Research efforts which are so-called meshfree methods have recently been directed towards eliminating or at least easing the requirement for meshing of the domain. In this paper, a new meshfree method for solving nonlinear boundary value problem, based on the bubble packing technique and Delaunay triangle is proposed. The method can be efficiently implemented to the problems with singularity by using formly distributed nodes.

• Journal of Aerospace System Engineering
• /
• v.1 no.1
• /
• pp.25-35
• /
• 2007

In this paper, the interlaminar crack problems of a fiber metal laminate (FML) under generalized plane deformation are studied using the theory of anisotropic elasticity. The crack is considered to be embedded in the matrix interlaminar region (including adhesive zone and resin rich zone) of the FML. Based on Fourier integral transformation and the stress matrix formulation, the current mixed boundary value problem is reduced to solving a system of Cauchy-type singular integral equations of the 1st kind. Within the theory of linear fracture mechanics, the stress intensity factors are defined on terms of the solutions of integral equations and numerical results are obtained for in-plane normal (mode I) crack surface loading. The effects of location and length of crack in the 3/2 and 2/1 ARALL, GLARE or CARE type FML's on the stress intensity factors are illustrated.

• Structural Engineering and Mechanics
• /
• v.26 no.1
• /
• pp.69-86
• /
• 2007

A four-noded plate bending quadrilateral (PBQ4) and an eight-noded plate bending quadrilateral (PBQ8) element based on Mindlin plate theory have been adopted for modeling the thick plates on elastic foundations using Winkler model. Transverse shear deformations have been included, and the stiffness matrices of the plate elements and the Winkler foundation stiffness matrices are developed using Finite Element Method based on thick plate theory. A computer program is coded for this purpose. Various loading and boundary conditions are considered, and examples from the literature are solved for comparison. Shear locking problem in the PBQ4 element is observed for small value of subgrade reaction and plate thickness. It is noted that prevention of shear locking problem in the analysis of the thin plate is generally possible by using element PBQ8. It can be concluded that, the element PBQ8 is more effective and reliable than element PBQ4 for solving problems of thin and thick plates on elastic foundations.