856 research outputs found
Addressing Integration Error for Polygonal Finite Elements Through Polynomial Projections: A Patch Test Connection
Polygonal finite elements generally do not pass the patch test as a result of
quadrature error in the evaluation of weak form integrals. In this work, we
examine the consequences of lack of polynomial consistency and show that it can
lead to a deterioration of convergence of the finite element solutions. We
propose a general remedy, inspired by techniques in the recent literature of
mimetic finite differences, for restoring consistency and thereby ensuring the
satisfaction of the patch test and recovering optimal rates of convergence. The
proposed approach, based on polynomial projections of the basis functions,
allows for the use of moderate number of integration points and brings the
computational cost of polygonal finite elements closer to that of the commonly
used linear triangles and bilinear quadrilaterals. Numerical studies of a
two-dimensional scalar diffusion problem accompany the theoretical
considerations
Trefftz Difference Schemes on Irregular Stencils
The recently developed Flexible Local Approximation MEthod (FLAME) produces
accurate difference schemes by replacing the usual Taylor expansion with
Trefftz functions -- local solutions of the underlying differential equation.
This paper advances and casts in a general form a significant modification of
FLAME proposed recently by Pinheiro & Webb: a least-squares fit instead of the
exact match of the approximate solution at the stencil nodes. As a consequence
of that, FLAME schemes can now be generated on irregular stencils with the
number of nodes substantially greater than the number of approximating
functions. The accuracy of the method is preserved but its robustness is
improved. For demonstration, the paper presents a number of numerical examples
in 2D and 3D: electrostatic (magnetostatic) particle interactions, scattering
of electromagnetic (acoustic) waves, and wave propagation in a photonic
crystal. The examples explore the role of the grid and stencil size, of the
number of approximating functions, and of the irregularity of the stencils.Comment: 28 pages, 12 figures; to be published in J Comp Phy
Application of Meshless Methods for Thermal Analysis
Many numerical and analytical schemes exist for solving heat transfer problems. The meshless method is a particularly attractive method that is receiving attention in the engineering and scientific modeling communities. The meshless method is simple, accurate, and requires no polygonalisation. In this study, we focus on the application of meshless methods using radial basis functions (RBFs) — which are simple to implement — for thermal problems. Radial basis functions are the natural generalization of univariate polynomial splines to a multivariate setting that work for arbitrary geometry with high dimensions. RBF functions depend only on the distance from some center point. Using distance functions, RBFs can be easily implemented to model heat transfer in arbitrary dimension or symmetry
Adaptive meshless refinement schemes for RBF-PUM collocation
In this paper we present an adaptive discretization technique for solving
elliptic partial differential equations via a collocation radial basis function
partition of unity method. In particular, we propose a new adaptive scheme
based on the construction of an error indicator and a refinement algorithm,
which used together turn out to be ad-hoc strategies within this framework. The
performance of the adaptive meshless refinement scheme is assessed by numerical
tests
Inverse heat conduction problems by using particular solutions
Based on the method of fundamental solutions, we develop in this paper a new computational method to solve two-dimensional transient heat conduction inverse problems. The main idea is to use particular solutions as radial basis functions (PSRBF) for approximation of the solutions to the inverse heat conduction problems. The heat conduction equations are first analyzed in the Laplace transformed domain and the Durbin inversion method is then used to determine the solutions in the time domain. Least-square and singular value decomposition (SVD) techniques are adopted to solve the ill-conditioned linear system of algebraic equations obtained from the proposed PSRBF method. To demonstrate the effectiveness and simplicity of this approach, several numerical examples are given with satisfactory accuracy and stability.Peer reviewe
Moving-boundary problems solved by adaptive radial basis functions
The objective of this paper is to present an alternative approach to the conventional level set methods for solving two-dimensional moving-boundary problems known as the passive transport. Moving boundaries are associated with time-dependent problems and the position of the boundaries need to be determined as a function of time and space. The level set method has become an attractive design tool for tracking, modeling and simulating the motion of free boundaries in fluid mechanics, combustion, computer animation and image processing. Recent research on the numerical method has focused on the idea of using a meshless methodology for the numerical solution of partial differential equations. In the present approach, the moving interface is captured by the level set method at all time with the zero contour of a smooth function known as the level set function. A new approach is used to solve a convective transport equation for advancing the level set function in time. This new approach is based on the asymmetric meshless collocation method and the adaptive greedy algorithm for trial subspaces selection. Numerical simulations are performed to verify the accuracy and stability of the new numerical scheme which is then applied to simulate a bubble that is moving, stretching and circulating in an ambient flow to demonstrate the performance of the new meshless approach. (C) 2010 Elsevier Ltd. All rights reserved
- …