• Journal of applied mathematics & informatics
• /
• v.31 no.1_2
• /
• pp.299-309
• /
• 2013

In this paper, we introduce a new numerical technique which we call fractional Chebyshev finite difference method (FChFD). The algorithm is based on a combination of the useful properties of Chebyshev polynomials approximation and finite difference method. We tested this technique to solve numerically fractional BVPs. The proposed technique is based on using matrix operator expressions which applies to the differential terms. The operational matrix method is derived in our approach in order to approximate the fractional derivatives. This operational matrix method can be regarded as a non-uniform finite difference scheme. The error bound for the fractional derivatives is introduced. The fractional derivatives are presented in terms of Caputo sense. The application of the method to fractional BVPs leads to algebraic systems which can be solved by an appropriate method. Several numerical examples are provided to confirm the accuracy and the effectiveness of the proposed method.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.22 no.6
• /
• pp.387-396
• /
• 2010

A 3-D numeric model for the confined transient flow under sea ground have been developed. This is FDM model using Gauss-Seidel SOR (successive over-relaxation). This model shows the similar head distribution pattern to Theis analytic solution and MODFLOW simulation. The input flowrate to the aquifer and discharge of well have been compared. And it have been found that mass balance is influenced by the weight factor ${\alpha}$, i.e. fullyimplicit method (${\alpha}$=1) shows 5% error, but when ${\alpha}$ becomes to 0.5(Crank and Nicolson method) the mass balance becomes worse and the model result diverges. And the convergency of the model is not much different when $\lambda$ (over-relaxation factor)=0.8~1.5, but when $\lambda$>1.5, the model result diverges. The test-run shows that the well discharge becomes smaller when another well is near. This model can cover the isotropy$(Kx{\neq}Ky{\neq}Kz)$ and inhomogeneity, and can be used for the selection of well site, discharge calculation, and head prediction in case of the artificial recharge etc.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.19 no.11
• /
• pp.205-212
• /
• 2018

In order to overcome shortcomings of commercial analysis program for solving certain geotechnical problems, finite element method adopting isoparametric quadrilateral element was selected as a tool for analyzing soil behavior and calculating process was programmed. Two examples were considered in order to verify reliability of the developed program. One of the two examples is the case of acting isotropic confining pressure on finite element and the other is the case of acting shear stress on the sides of the finite element. Isoparametric quadrilateral element was considered as the finite element and displacements in the element can be expressed by node displacements and shape functions in the considered element. Calculating process for determining strain which is defined by derivatives using global coordinates was coded using the Jacobian and the natural coordinates. Four point Gauss rule was adopted to convert double integral which defines stiffness of the element into numerical integration. As a result of executing analysis of the finite element under isotropic confining pressure, calculated stress corresponding to four Gauss points and center of the element were equal to the confining pressure. In addition, according to the analyzed results for the element under shear stress, horizontal stresses and vertical stresses were varied with positions in the element and the magnitudes and distribution pattern of the stresses were thought to be rational.

• Structural Engineering and Mechanics
• /
• v.57 no.2
• /
• pp.281-294
• /
• 2016

The paper is devoted to study a mesh-free analysis method of structural elements of engineering structures based on B-spline Wavelet Basis Function. First, by employing the moving-least square method and the weighted residual method to solve the structural displacement field, the control equations and the stiffness equations are obtained. And then constructs the displacement field of the structure by using the m-order B-spline wavelet basis function as a weight function. In the end, the paper selects the plane beam structure and the structure with opening hole to carry out numerical analysis of deformation and stress. The Finite Element Method calculation results are compared with the results of the method proposed, and the calculation results of the relative error norm is compared with Gauss weight function as weight function. Therefore, the clarification verified the validity and accuracy of the proposed method.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.17 no.3
• /
• pp.279-293
• /
• 2004

The analysis of linear static and free vibration problems of isotropic and laminated composite plates and shells is performed by the improved 9-node shell element with the new strain displacement relationship. In that relationship, the effect of new additional terms between the bending strain and displacement has been investigated in the warping problem. Natural co ordinate based strains, stresses and constitutive equations are used. The assumed natural strain method is used to alleviate both membrane and shear locking behavior from the element. The Lanczos method is employed in the calculation of the eigenvalues of laminated composite structures and the Gauss integration rule is adopted to evaluate the mass matrix. The numerical examples are compared with the analytical solutions to validate the current formulation and the results presented could be useful for the understanding of the behaviour of laminates under free vibration conditions.

• Membrane Journal
• /
• v.22 no.6
• /
• pp.470-480
• /
• 2012

In this study, for a linear tubular membrane with constant diameter the mass balance considering permeate velocity of the each unit length was established. On this basis, mathematical modelling of flows in a pipe was solved using nonlinear second order differential equations as well as steady-state equation. Since this equation is nonlinear, Gauss-Seidel method or another iteration method were used to solve the differential equations. Simulation algorithm for numerical solutions was presented. Also since the permeate flow is varied as operating condition, the solution of equations at each conditions using numerical integrations such as Simpson's rules was used. In order to analyze and compare simulation results, we have performed experiments using a hollow fiber membrane with almost identical tubular membrane. Comparison of theoretical and experimental results, pressure drop, flow rate, and permeate flow in a hollow fiber membrane, were illustrated.

• Journal of the Society of Naval Architects of Korea
• /
• v.38 no.4
• /
• pp.30-38
• /
• 2001

The current is considered in the conventional manoeuvering equation. This equation is represented as the nonlinear state and measurement equations in which external forces and the direction and the velocity of current are augmented as that variables. The external forces are modeled as the third-order Gauss-Markov processes and the direction and the velocity of current are assumed to be constant. The augmented state variables are estimated with extended Kalman-Bucy filter and the fixed-interval smoother. While Hwang estimated motion state variables, hydrodynamic coefficients and the current variables simultaneously by using extended Kalman filter, external forces of surge, sway and yaw and the direction and the velocity of current are the only parameters to be estimated in the estimation-before-modeling method. The current variables are satisfactorily estimated in simulation process where the measurement noise is present.

• Journal of KSNVE
• /
• v.9 no.5
• /
• pp.914-921
• /
• 1999

There are many theoretical investigations to analyze higher order mode of reactive type single expansion chambers with offset inlet/outlet locations. But the conventional method has the restriction that the ratio between the area of inlet(or outlet) pipe and that of chamber must be natural number. In the paper, a new method was suggested to apply the Kim's method to silencer with circular cross-section. Not only the offset location but also the magnitude of cross-section area of inlet/outlet pipe can be considered by the suggested method. The predictions by this new method also compared with those by the finite element method and Munjal's method in order to verify the accuracy of the suggested method presented here.

• Proceedings of the Korean Statistical Society Conference
• /
• 2004.11a
• /
• pp.1-8
• /
• 2004

We address the problem of parameter estimation in multivariate distributions under ignorable non-monotone missing data. The factoring likelihood method for monotone missing data, termed by Robin (1974), is extended to a more general case of non-monotone missing data. The proposed method is algebraically equivalent to the Newton-Raphson method for the observed likelihood, but avoids the burden of computing the first and the second partial derivatives of the observed likelihood Instead, the maximum likelihood estimates and their information matrices for each partition of the data set are computed separately and combined naturally using the generalized least squares method. A numerical example is also presented to illustrate the method.

• 한국연소학회:학술대회논문집
• /
• 1999.10a
• /
• pp.219-230
• /
• 1999

The Equations of Chemical kinetics are very stiff, which forces the use of an implicit scheme. The problem of implicit scheme, however, is that the jacobian must be solved at each time step. In this paper, we examined the methodology that can be stable without full chemical jacobian, This method is derived by applying the different time steps to the chemical source term. And the lower triangular chemical jacobian is derived. This is called the preconditioned time differencing method and represents partial implicit method. We show that this method is more stable in chemical kinetics than the full implicit method and that this is more efficient in supersonic combustion problem than the full jacobian method with same accuracy.