iMechanica - Comments for "rate of convergence for Peridynamics"
https://imechanica.org/node/8667
Comments for "rate of convergence for Peridynamics"enSecond order accuracy
https://imechanica.org/comment/15283#comment-15283
<a id="comment-15283"></a>
<p><em>In reply to <a href="https://imechanica.org/node/8667">rate of convergence for Peridynamics</a></em></p>
<div class="field field-name-comment-body field-type-text-long field-label-hidden"><div class="field-items"><div class="field-item even"><p>
Phil,
</p>
<p>
Yes, second order accuracy on smooth problems does occur for linear peridynamics.
</p>
<p>
My colleague Max Gunzburger at FSU and his student have a paper under review where a 1d finite peridynamic bar is discretized with continuous and discontinuous FE and convergence rates verified to both smooth and discontinuous manufactured solutions. Another colleague Qiang Du at PSU and his student show the approximation error for both a 1d finite bar and 2d problem on a square; this paper has been accepted for publication in the SIAM J. of Num. Anal. You might contact both Max and Qiang for copies.
</p>
<p>
Any implementation has details, such as those you list, but mass lumping is not an issue if discontinuous FE are used. This may (or maybe not) seem obvious but is a point that deserves a bit of discussion. Discontinuous Galerkin methods for conventional weak formulations of linear elasticity are nonconforming. Asymptoptic convergence to discontinuous problems only occurs if the element boundaries are aligned with the discontinuities. On the other hand, a discontinuous FE for peridynamics is conforming and any discontunuities in the solutiuon only affect the rate of convergence.
</p>
<p>
Thanks for your interest.
</p>
<p>
Rich
</p>
</div></div></div><ul class="links inline"><li class="comment_forbidden first last"><span><a href="/user/login?destination=node/8667%23comment-form">Log in</a> or <a href="/user/register?destination=node/8667%23comment-form">register</a> to post comments</span></li>
</ul>Thu, 12 Aug 2010 03:07:32 +0000Rich Lehoucqcomment 15283 at https://imechanica.orgRich, thank you for the
https://imechanica.org/comment/15269#comment-15269
<a id="comment-15269"></a>
<p><em>In reply to <a href="https://imechanica.org/comment/15238#comment-15238">Rate of convergence for peridynamics</a></em></p>
<div class="field field-name-comment-body field-type-text-long field-label-hidden"><div class="field-items"><div class="field-item even"><p>
Rich, thank you for the pointer to the SAND report; the refs from chapter one seem to go in the right direction and several of them were new to me. I wanted to follow up on a few points you made, and make sure I'm heading in the right direction. To start things off we make the necessary assumptions of existence, uniqueness, and smoothness (no cracks) in space and time. So then we can make comparisons to typical FE behavior and to manufactured solutions. In most methods the details can get a bit tricky: the way boundary conditions are applied, the way time stepping is initialized, mass lumping, order of reproducibility for elements or MLS approximations, the definition of error, etc. etc. What I ultimately hope to learn is the specifics of all those picky details that are necessary to get second order convergence with, say, a typical L-2 norm. If Peridynamics can demonstrate second order accuracy for smooth problems that would be fantastic, and lower order for discontinuous problems is to be expected. Cheers,
</p>
<p>
Phil
</p>
</div></div></div><ul class="links inline"><li class="comment_forbidden first last"><span><a href="/user/login?destination=node/8667%23comment-form">Log in</a> or <a href="/user/register?destination=node/8667%23comment-form">register</a> to post comments</span></li>
</ul>Tue, 10 Aug 2010 21:27:17 +0000wallstedtcomment 15269 at https://imechanica.orgRate of convergence for peridynamics
https://imechanica.org/comment/15238#comment-15238
<a id="comment-15238"></a>
<p><em>In reply to <a href="https://imechanica.org/node/8667">rate of convergence for Peridynamics</a></em></p>
<div class="field field-name-comment-body field-type-text-long field-label-hidden"><div class="field-items"><div class="field-item even"><p>
Phil,
</p>
<p>
There are several published papers that consider the rate of convergence as the mesh size is refined. Please see chapter 1 of <strong>Peridynamic Theory of Solid Mechanics</strong> by <a href="http://www.sandia.gov/%7Esasilli/">Stewart Silling</a> and <a href="http://www.cs.sandia.gov/~rblehou/">I</a>, to appear in <em><strong>Advances in Applied Mechanics</strong></em>, and available as <a href="http://www.sandia.gov/%7Erblehou/snl-sand2010-1233J.pdf">Technical report SAND 2010-1233J</a>, for pointers to the literature. Keep in mind that before one can talk about reducing error, it must be demonstrated that there is a solution and that it's unique. These sort of mathematical questions are being addressed and chapter 1 cites the relevant papers and work in progress.
</p>
<p>
Note that first or second order convergence rates depend upon the norm used and the <em>regularity</em>, i.e. smoothness, of the solution and the quality of the approximation. One of the interesting aspects of peridynamics is that the equations are valid even if the deformation is discontinuous and the location of the discontunuities are not required to be known in advance. The resulting mathematical theory and numerical analysis then explain that the best possible convergence rate may only be square root. However, away from the discontunity, the convergence rate will be first or second order, depending how the error is measured and the quality of the approximation.
</p>
<p>
Regards,
</p>
<p>
Rich Lehoucq</p>
</div></div></div><ul class="links inline"><li class="comment_forbidden first last"><span><a href="/user/login?destination=node/8667%23comment-form">Log in</a> or <a href="/user/register?destination=node/8667%23comment-form">register</a> to post comments</span></li>
</ul>Sat, 07 Aug 2010 04:14:59 +0000Rich Lehoucqcomment 15238 at https://imechanica.org