Search results

The purpose of this paper is to develop a new local radial basis function collocation method (LRBFCM) for one‐domain solving of the non‐linear convection‐diffusion equation, as it appears in mixture continuum formulation of the energy transport in solid‐liquid phase change systems.

The method is structured on multiquadrics radial basis functions. The collocation is made locally over a set of overlapping domains of influence and the time stepping is performed in an explicit way. Only small systems of linear equations with the dimension of the number of nodes in the domain of influence have to be solved for each node. The method does not require polygonisation (meshing). The solution is found only on a set of nodes.

The computational effort grows roughly linearly with the number of the nodes. Results are compared with the existing steady analytical solutions for one‐dimensional convective‐diffusive problem with and without phase change. Regular and randomly displaced node arrangements have been employed. The solution is compared with the results of the classical finite volume method. Excellent agreement with analytical solution and reference numerical method has been found.

A realistic two‐dimensional non‐linear industrial test associated with direct‐chill, continuously cast aluminium alloy slab is presented.

A new meshless method is presented which is simple, efficient, accurate, and applicable in industrial convective‐diffusive solid‐liquid phase‐change problems with non‐linear material properties.

This paper gives a bibliographical review of the finite element and boundary element parallel processing techniques from the theoretical and application points of view. Topics include: theory – domain decomposition/partitioning, load balancing, parallel solvers/algorithms, parallel mesh generation, adaptive methods, and visualization/graphics; applications – structural mechanics problems, dynamic problems, material/geometrical non‐linear problems, contact problems, fracture mechanics, field problems, coupled problems, sensitivity and optimization, and other problems; hardware and software environments – hardware environments, programming techniques, and software development and presentations. The bibliography at the end of this paper contains 850 references to papers, conference proceedings and theses/dissertations dealing with presented subjects that were published between 1996 and 2002.

The purpose of this study is to provide a robust numerical method for a two parameter singularly perturbed delay parabolic initial boundary value problem (IBVP).

To solve the problem, the authors have used a hybrid scheme combining the midpoint scheme, the upwind scheme and the second-order central difference scheme for the spatial derivatives. The backward Euler scheme on a uniform mesh is used to approximate the time derivative. Here, the authors have used Shishkin type meshes for spatial discretization.

It is observed that the proposed method converges uniformly with almost second-order spatial accuracy with respect to the discrete maximum norm.

This paper deals with the numerical study of a two parameter singularly perturbed delay parabolic IBVP. To solve the problem, the authors have used a hybrid scheme combining the midpoint scheme, the upwind scheme and the second-order central difference scheme for the spatial derivatives. The backward Euler scheme on a uniform mesh is used to approximate the time derivative. The convergence analysis is carried out. It is observed that the proposed method converges uniformly with almost second-order spatial accuracy with respect to the discrete maximum norm. Numerical experiments illustrate the efficiency of the proposed scheme.

Gives a bibliographical review of the finite element meshing and remeshing from the theoretical as well as practical points of view. Topics such as adaptive techniques for meshing and remeshing, parallel processing in the finite element modelling, etc. are also included. The bibliography at the end of this paper contains 1,727 references to papers, conference proceedings and theses/dissertations dealing with presented subjects that were published between 1990 and 2001.

This paper is concerned with the numerical solution of multi‐dimensional convection dominated convection‐diffusion problems. These problems are characterized by a large parameter, K, multiplying the convection terms. The goal of this work is the development and analysis of effective preconditioners for iteratively solving the large system of linear equations arising from various finite element and finite difference discretizations with grid size h. When centered finite difference schemes and standard Galerkin finite element methods are used, h must be related to K by the stability constraint, Kh ≤ C₀, where the constant C₀ is sufficiently small. A class of preconditioners is developed that significantly reduces the condition number for large K and small h. Furthermore, these preconditioners are inexpensive to implement and well suited for parallel computation. It is shown that under suitable assumptions, the number of iterations remains bounded as h ↓0 with K fixed and, at worst, grows slowly as K ↓ ∞. Numerical results are presented illustrating the theory. It is also shown how to apply the theoretical results to more general convection‐diffusion problems and alternative discretizations (including streamline diffusion methods) that remain stable as Kh ↓ ∞.

This paper presents a boundary element formulation for transient convection‐diffusion problems employing the fundamental solution of the corresponding steady‐state equation with constant coefficients and a dual reciprocity approximation. The formulation allows the mathematical problem to be described in terms of boundary values only. Numerical results show that the BEM does not present oscillations or damping of the wave front as appear in other numerical techniques.

The paper presents the numerical performance of the preconditioned generalized conjugate gradient (PGCG) methods in solving non‐linear convection — diffusion equations. Three non‐linear systems which describe a non‐isothermal chemical reactor, the chemically driven convection in a porous medium and the incompressible steady flow past a sphere are the test problems. The standard second order accurate centred finite difference scheme is used to discretize the models equations. The discrete approximations are solved with a double iterative process using the Newton method as outer iteration and the PGCG algorithm as inner iteration. Three PGCG techniques, which emerge to be the best performing, are tested. Laplace‐type operators are employed for preconditioning. The results show that the convergence of the PGCG methods depends strongly on the convection—diffusion ratio. The most robust algorithm is GMRES. But even with GMRES non‐convergence occurs when the convection—diffusion ratio exceeds a limit value. This value seems to be influenced by the non‐linearity type.

A hybrid numerical‐analytical approach, based on recent developments in the generalized integral transform technique, is presented for the solution of a class of non‐linear transient convection‐diffusion problems. The original partial differential equation is integral transformed into a denumerable system of coupled non‐linear ordinary differential equations, which is numerically solved for the transformed potentials. The hybrid analysis convergence is illustrated by considering the one‐dimensional non‐linear Burgers equation and numerical results are presented for increasing truncation orders of the infinite ODE system.

In this paper we investigate the performance of CGS, BCGSTAB and GMRES with ILU preconditioner for solving convection‐diffusion problems. Numerical experiments indicate that BCGSTAB appears to be an efficient and stable method. CGS sometimes suffers from severe numerical instability. GMRES shows a higher suitability and stability but the overall convergence rate may be lower.

In this paper new cubic v‐splines monotonic one‐dimensional profiles are presented, for the finite volume solution of convection‐diffusion problems. By studying the profile in normalized variables, some weight functions have been determined for the profile. Being free of the requirement that the volumes be equal, the volume size can be reduced where needed. Numerical properties of the proposed method were formally analysed and are confirmed by numerical examples included here.