• Title/Summary/Keyword: Numerical solution.

Search Result 4,022, Processing Time 0.031 seconds

Numerical Solution of Nonlinear Diffusion in One Dimensional Porous Medium Using Hybrid SOR Method

  • Jackel Vui Lung, Chew;Elayaraja, Aruchunan;Andang, Sunarto;Jumat, Sulaiman
    • Kyungpook Mathematical Journal
    • /
    • v.62 no.4
    • /
    • pp.699-713
    • /
    • 2022
  • This paper proposes a hybrid successive over-relaxation iterative method for the numerical solution of a nonlinear diffusion in a one-dimensional porous medium. The considered mathematical model is discretized using a computational complexity reduction scheme called half-sweep finite differences. The local truncation error and the analysis of the stability of the scheme are discussed. The proposed iterative method, which uses explicit group technique and modified successive over-relaxation, is formulated systematically. This method improves the efficiency of obtaining the solution in terms of total iterations and program elapsed time. The accuracy of the proposed method, which is measured using the magnitude of absolute errors, is promising. Numerical convergence tests of the proposed method are also provided. Some numerical experiments are delivered using initial-boundary value problems to show the superiority of the proposed method against some existing numerical methods.

Accurate semi-analytical solution for nonlinear vibration of conservative mechanical problems

  • Bayat, Mahmoud;Pakar, Iman
    • Structural Engineering and Mechanics
    • /
    • v.61 no.5
    • /
    • pp.657-661
    • /
    • 2017
  • In this paper, it has been tried to propose a new semi analytical approach for solving nonlinear vibration of conservative systems. Hamiltonian approach is presented and applied to high nonlinear vibration systems. Hamiltonian approach leads us to high accurate solution using only one iteration. The method doesn't need any small perturbation and sufficiently accurate to both linear and nonlinear problems in engineering. The results are compared with numerical solution using Runge-Kutta-algorithm. The procedure of numerical solution are presented in detail. Hamiltonian approach could be simply apply to other powerfully non-natural oscillations and it could be found widely feasible in engineering and science.

The finite element model of pre-twisted Euler beam based on general displacement solution

  • Huang, Ying;Chen, Changhong;Zou, Haoran;Yao, Yao
    • Structural Engineering and Mechanics
    • /
    • v.69 no.5
    • /
    • pp.479-486
    • /
    • 2019
  • Based on the displacement general solution of a pre-twisted Euler-Bernoulli beam, the shape function and stiffness matrix are deduced, and a new finite element model is proposed. Comparison analyses are made between the new proposed numerical model based on displacement general solution and the ANSYS solution by Beam188 element based on infinite approach. The results show that developed numerical model is available for the pre-twisted Euler-Bernoulli beam, and that also provide an accuracy finite element model for the numerical analysis. The effects of pre-twisted angle and flexural stiffness ratio on the mechanical property are also investigated.

Combined Wave Reflection and Diffraction near the Upright Breakwater (직립 방파제 주위에서 파랑의 반사 및 회절의 혼합)

  • Shin, Seung Ho;Gug, Seung Gi;Yeom, Won Gi;Lee, Joong Woo
    • Journal of Korean Port Research
    • /
    • v.5 no.1
    • /
    • pp.71-81
    • /
    • 1991
  • This study deals with the analytical and numerical solution for the combined wave reflection and diffraction near the impermeable rigid upright breakwater, subject to the excitation of a plane simple harmonic wave coming from infinity. Three cases are presented : a) the analytical solution near a thin semi-infinite breakwater, b) the analytical solution near the semi-infinite breakwaters of arbitrary edge angles, $30^{\circ},\;45^{\circ},\;and\;90^{\circ}$, c) the numerical solution near a detached thin breakwater the results are presented in amplification factor and wave height diagrams. Moreover, the amplification factors near the structure(2 wavelength before and behind the structure) are compared for the given cases. A finite difference technique for the numerical solution was applied to the integral equation obtained for the wave potential.

  • PDF

The Solution of Mild-Slope Equation using Power Series (멱급수를 이용한 완경사 방정식의 해)

  • Jung, Tae-Hwa;Lee, Seung-Oh;Park, Jin-Ho;Cho, Yong-Sik
    • Journal of the Korean Society of Hazard Mitigation
    • /
    • v.8 no.1
    • /
    • pp.133-138
    • /
    • 2008
  • To analyze incident waves traveling from the deep ocean is very important in that it is based on resolving problems occurred in coastal areas. In general, numerical models and analytical solutions are used to analyze wave transformation. Although a numerical model can be applied to various bottoms and wave conditions, it may have some cumbersome numerical errors. On the other hand, an analytical solution has an advantage of obtaining the solution quickly and accurately without numerical errors. The analytical solution can, however, be utilized only for specific conditions. In this study, the analytical solution of the mild-slope equation has been developed. It can be applied to various conditions combing a numerical technique and an analytical approach while minimizing the numerical errors. As a result of comparing the obtained solutions in this study with those of the previously developed numerical model, A good agreement was observed.

A Study of the Effect of Grouting Region on the Solution of Line Source Analysis (그라우팅 영역이 선형열원 해석에 미치는 영향에 관한 연구)

  • Lee, Se-Kyoun;Woo, Joung-Son;Ro, Jeong-Geun
    • Journal of Energy Engineering
    • /
    • v.19 no.3
    • /
    • pp.143-150
    • /
    • 2010
  • Line source method of borehole system assumes the entire surrounding medium is uniform. However, thermal properties of grouting region are considerably different from those of surrounding soil. In this study we investigate the effect of grouting materials on the solution of line source method with the aid of numerical analysis. This numerical model generates the temperature of borehole fluid with which line source solution can be obtained. Then this solution can be compared with input condition of numerical model. The results of this comparison show that thermal conductivity and borehole thermal resistance of line source solution are approximately 86% and 91% of the input condition of numerical model. Chart method is developed in this study to find the numerical input conditions (thermal conductivity and borehole thermal resistance) from the line source solution. Thermal response test of test borehole is conducted, the results of which are approximately consistent with the Chart method. Thermal property changes of grouting materials on the line source solution are also examined.

Numerical Simulation of Solution Droplets and Falling Films in Horizontal Tube Absorbers

  • Phan Thanh-Tong;Lee Ho-Saeng;Yoon Jung-In;Kim Eun-Pil
    • Journal of Advanced Marine Engineering and Technology
    • /
    • v.30 no.5
    • /
    • pp.597-607
    • /
    • 2006
  • This paper presents a numerical simulation of the behavior of the LiBr solution droplets and falling films in horizontal tube banks of absorber. The model developed here accounts for the details of the droplets formation and impact process for absorption on horizontal tubes including the heat transfer from solution film to the tube wall. Especially. the characteristic of unsteady behavior of solution flow has been investigated. Flow visualization studies shown that the solution droplets and falling films have some of the complex characteristics. It is found that. with the numerical conditions similar to the operating condition of an actual absorption chiller/heater, the outlet solution temperature and heat flux from solution film to the tube wall have a stable periodic behavior with time. The solution droplets and falling films in horizontal tube banks of absorber is a periodic unsteady flow. The results from this model are compared with previous experimental observation taken with a high-speed digital video camera and shown good agreement.

A PARAMETRIC SCHEME FOR THE NUMERICAL SOLUTION OF THE BOUSSINESQ EQUATION

  • Bratsos, A.G.
    • Journal of applied mathematics & informatics
    • /
    • v.8 no.1
    • /
    • pp.45-57
    • /
    • 2001
  • A parametric scheme is proposed for the numerical solution of the nonlinear Boussinesq equation. The numerical method is developed by approximating the time and the space partical derivatives by finite-difference re placements and the nonlinear term by an appropriate linearized scheme. The resulting finite-difference method is analyzed for local truncation error and stability. The results of a number of numerical experiments are given for both the single and the double-soliton wave. AMS Mathematics Subject Classification : 65J15, 47H17, 49D15.

Numerical solution of linear elasticity by preconditioning cubic spline collocation

  • Lee, Yong-Hun
    • Communications of the Korean Mathematical Society
    • /
    • v.11 no.3
    • /
    • pp.867-880
    • /
    • 1996
  • Numerical approximations to the linear elasticity are traditionally based on the finite element method. In this paper we propose a new formulation based on the cubic spline collocation method for linear elastic problem on the unit square. We present several numerical results for the eigenvalues of the matrix represented by cubic collocation method and preconditioner matrix which is preconditioned by FEM and FDM. Finally we present the numerical solution for some example equation.

  • PDF

HOPF BIFURCATION IN NUMERICAL APPROXIMATION OF THE SUNFLOWER EQUATION

  • Zhang Chunrui;Zheng Baodong
    • Journal of applied mathematics & informatics
    • /
    • v.22 no.1_2
    • /
    • pp.113-124
    • /
    • 2006
  • In this paper we consider the numerical solution of the sunflower equation. We prove that if the sunflower equation has a Hopf bifurcation point at a = ao, then the numerical solution with the Euler-method of the equation has a Hopf bifurcation point at ah = ao + O(h).