meshfree

How to deal with noisy data in a mesh free Galerkin method

I'm currently working with a particle tracking velocimetry method, specifically trying to extract shear and vorticity data from the non-gridded velocity results. Of several methods used, the most efficient and accurate was based on a non-Sibsonian element free method. Unlike a more typical Galerkin problem which finds displacements by solving a PDE, this method uses the displacements of natural neighbors to find local flow gradients. 
But here's the rub: this method is very sensitive to noisy data. Small errors in the displacement data cause large errors in the calculated gradients. Thus far smoothing the displacement data has been unsuccessful in correcting this problem. 


Hennadiy Netuzhylov's picture

CMMSE 2010 - Call for papers - Minisymposium on Sampling theory and meshfree methods

Please note the Minisymposium on Sampling theory and meshfree methods as a part of

the International Conference on Computational and Mathematical methods in Science and Engineering
CMMSE 2010, 26-30 June 2010, Almeria, Andalucia, Spain
--   http://gsii.usal.es/~CMMSE  --

 This mini-symposium aims to bring together such related areas of
mathematics and computational mechanics as Sampling Theory and Meshfree
Numerical Methods.


Stephane Bordas's picture

UK Newton Post-Doc Fellowship in Computational Mechanics

 
Dear All,

A new multi-million pound initiative to fund research collaborations and improve links between UK and overseas researchers has been launched.

The Newton International Fellowships aim to attract the most promising, early stage, post-doctoral researchers working overseas, who do not hold UK citizenship, in the fields of humanities, engineering, natural
and social sciences.

The scheme provides funding to successful candidates for up to 2 years to work with research groups at a UK research Institution and to establish long term international collaborations.


Ettore Barbieri's picture

The Future of Meshless Methods

I joined imechanica almost a year ago and I've been frequently following its interesting discussions, even the most animated ones. I think that a place like this is ideal to foster the exchange of ideas in the scientific community;

Moreover it is fantastic as a simple student like me can interact and easily ask questions to the most important researcher in the field of mechanics.

Hence, I thought it would have been the right place to pose a question which I believe is quite controversial. The debate I would like to open is about the future of meshless methods, are they still valid? It is worth to keep investigation in this area?


N. Sukumar's picture

Epi-convergence (max-ent bases), crack growth

In the attached paper, we have used Variational Analysis techniques (in particular, the theory of epi-convergence) to prove the continuity of maximum-entropy basis functions. In general, for non-smooth functionals, moving objectives and/or constraints, the tools of Newton-Leibniz calculus (gradient, point-convergence) prove to be insufficient; notions of set-valued mappings, set-convergence, etc., are required. Epi-convergence bears close affinity to Gamma- or Mosco-convergence (widely used in the mathematical treatment of martensitic phase transformations). The introductory material on convex analysis and epi-convergence had to be omitted in the revised version; hence the material is by no means self-contained. Here are a few more pointers that would prove to be helpful. Our main point of reference is Variational Analysis by RTR and RJBW; the Princeton Classic Convex Analysis by RTR provides the important tools in convex analysis. For convex optimization, the text Convex Optimization by SB and LV (available online) is excellent. The lecture slides provide a very nice (and gentle) introduction to some of the important concepts in convex analysis. The epigraphical landscape is very rich, and many of the applications would resonate with mechanicians.

On a different topic (non-planar crack growth), we have coupled the x-fem to a new fast marching algorithm. Here are couple of animations on growth of an inclined penny crack in tension (unstructured tetrahedral mesh with just over 12K nodes): larger `time' increment and smaller `time' increment. This is joint-work with Chopp, Bechet and Moes (NSF-OISE project). I will update this page as and when more relevant links are available.


N. Sukumar's picture

3. Meshfree methods and the patch test

We are in the framework of small-strain two-dimensional linear elasticity without any body forces. Consider a domain that is discretized by a union of triangles and/or quadrilaterals (`patch of elements').  For C0 conforming approximations such as triangular/quadrilateral finite elements, the finite element approximation can exactly reproduce an arbitrary linear displacement field. Hence, if the exact solution is linear, then the finite element solution must match (within machine precision) the exact solution. In simple terms, passing the patch test for linear elasticity with standard conforming finite elements provides verification of one's implementation and is used to assess the same when new elements are proposed. For conforming elements, it is a sufficient condition for convergence (2nd order PDEs), and hence is the first problem that is solved when a new element/method is proposed. To carry out the patch test, the following steps are performed:


N. Sukumar's picture

2. Is a mesh required in meshfree methods?

In meshfree (this is more in vogue than the term meshless) methods, two key steps need to be mentioned: (A) construction of the trial and test approximations; and (B) numerical evaluation of the weak form (Galerkin or Rayleigh-Ritz procedure) integrals, which lead to a linear system of equations (Kd = f). In meshfree Galerkin methods, the main departure from FEM is in (A): meshfree approximation schemes (linear combination of basis functions) are constructed independent of an underlying mesh (union of elements).

However, since a Galerkin method is typically used in solid mechanics applications, (B) arises and the weak form integrals need to be evaluated. Three main directions have been pursued to evaluate these integrals:


Zhigang Suo's picture

Where can I read about the basic ideas of the meshfree methods?

For someone with a background in solid mechanics and finite element methods, where should he go to read up on the elementary ideas of the meshfree methods?


Syndicate content