Home  About us  Contact
 

Traditional Numerical Methods (traditional + numerical_methods)

Distribution by Scientific Domains
Engineering100%

Selected Abstracts

Solute transport in sand and chalk: a probabilistic approach

HYDROLOGICAL PROCESSES, Issue 5 2006
E. Carlier
Abstract A probabilistic approach is used to simulate particle tracking for two types of porous medium. The first is sand grains with a single intergranular porosity. Particle tracking is carried out by advection and dispersion. The second is chalk granulates with intergranular and matrix porosities. Sorption can occur with advection and dispersion during particle tracking. Particle tracking is modelled as the sum of elementary steps with independent random variables in the sand medium. An exponential distribution is obtained for each elementary step and shows that the whole process is Markovian. A Gamma distribution or probability density function is then deduced. The relationships between dispersivity and the elementary step are given using the central limit theorem. Particle tracking in the chalky medium is a non-Markovian process. The probability density function depends on a power of the distance. Experimental simulations by dye tracer tests on a column have been performed for different distances and discharges. The probabilistic approach computations are in good agreement with the experimental data. The probabilistic computation seems an interesting and complementary approach to simulate transfer phenomena in porous media with respect to the traditional numerical methods. Copyright © 2006 John Wiley & Sons, Ltd. [source]

Application of radial basis meshless methods to direct and inverse biharmonic boundary value problems

INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING, Issue 4 2005
Jichun Li
Abstract In this paper, we develop a non-iterative way to solve biharmonic boundary value problems by using a radial basis meshless method. This is an original application of meshless method to solving inverse problems without any iteration, since traditional numerical methods for inverse boundary value problems mainly are iterative and hence very time-consuming. Numerical examples are presented for inverse biharmonic boundary value problems and corresponding direct problems, since solving direct problems is a preliminary step for inverse problems. All our examples of direct and inverse problems are solved within seconds in CPU time on a standard PC, which makes our proposed technique a great potential candidate for wide-spread applications to other inverse problems. Copyright © 2004 John Wiley & Sons, Ltd. [source]

Lower-bound limit analysis by using the EFG method and non-linear programming

INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Issue 3 2008
Shenshen Chen
Abstract Intended to avoid the complicated computations of elasto-plastic incremental analysis, limit analysis is an appealing direct method for determining the load-carrying capacity of structures. On the basis of the static limit analysis theorem, a solution procedure for lower-bound limit analysis is presented firstly, making use of the element-free Galerkin (EFG) method rather than traditional numerical methods such as the finite element method and boundary element method. The numerical implementation is very simple and convenient because it is only necessary to construct an array of nodes in the domain under consideration. The reduced-basis technique is adopted to solve the mathematical programming iteratively in a sequence of reduced self-equilibrium stress subspaces with very low dimensions. The self-equilibrium stress field is expressed by a linear combination of several self-equilibrium stress basis vectors with parameters to be determined. These self-equilibrium stress basis vectors are generated by performing an equilibrium iteration procedure during elasto-plastic incremental analysis. The Complex method is used to solve these non-linear programming sub-problems and determine the maximal load amplifier. Numerical examples show that it is feasible and effective to solve the problems of limit analysis by using the EFG method and non-linear programming. Copyright © 2007 John Wiley & Sons, Ltd. [source]

Use of slopelimiter techniques in traditional numerical methods for multi-phase flow in pipelines and wells

INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 7 2005
R. J. Lorentzen
Abstract The aim of this paper is to show how simple and traditional methods for simulating multi-phase flow can be improved by introducing higher order accuracy. Numerical diffusion is reduced to a minimum by using slopelimiter techniques, and better predictions of flow rates and pressure are obtained. Slopelimiter techniques, originally developed to achieve higher order of accuracy in Godunov's method, is applied to a method following a finite element approach and a predictor,corrector shooting technique. These methods are tested and compared to a Godunov-type scheme recently developed for multi-phase flow. Implementation of Godunov-type schemes for multi-phase flow tends to be a complicated and challenging task. Introducing the slopelimiter techniques in the finite element approach and the predictor,corrector shooting technique is however simple, and provides an overall reduction of the numerical diffusion. The focus is on using these techniques to improve the mass transport description, since this is the main concern in the applications needed. The presented schemes represent different semi-implicit approaches for simulating multi-phase flow. An evaluation of higher order extensions, as well as a comparison by itself, is of large interest. We present a model for two-phase flow in pipelines and wells, and an outline of the numerical methods and the extensions to second order spatial accuracy. Several examples motivated by applications in underbalanced drilling are presented, and the advantages of using higher order schemes are illustrated. Copyright © 2005 John Wiley & Sons, Ltd. [source]

A novel finite point method for flow simulation

INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 12 2002
M. Cheng
Abstract A novel finite point method is developed to simulate flow problems. The mashes in the traditional numerical methods are supplanted by the distribution of points in the calculation domain. A local interpolation based on the properties of Taylor series expansion is used to construct an approximation for unknown functions and their derivatives. An upwind-dominated scheme is proposed to efficiently handle the non-linear convection. Comparison with the finite difference solutions for the two-dimensional driven cavity flow and the experimental results for flow around a cylinder shows that the present method is capable of satisfactorily predicting the flow separation characteristic. The present algorithm is simple and flexible for complex geometric boundary. The influence of the point distribution on computation time and accuracy of results is included. Copyright © 2002 John Wiley & Sons, Ltd. [source]