Engineering topics
Publications and source records attributed to Seleson, Pablo.
This report describes the objectives, activities, and outcomes of the Workshop on Experimental and Computational Fracture Mechanics, which was held at Louisiana State University in Baton Rouge, Louisiana on February 26–28, 2020.
Diffusion-type problems in (nearly) unbounded domains play important roles in various fields of fluid dynamics, biology, and materials science. The aim of this paper is to construct accurate absorbing boundary conditions (ABCs) suitable for classical (local) as well as nonlocal peridynamic (PD) diffusion models. The main focus of the present study is on the PD diffusion formulation. The majority of the PD diffusion models proposed so far are applied to bounded domains only. In this study, we propose an effective way to handle unbounded domains both with PD and classical diffusion models. For the former, we employ a meshfree discretization, whereas for the latter the finite element method (FEM) is employed. The proposed ABCs are time-dependent and Dirichlet-type, making the approach easy to implement in the available models. The performance of the approach, in terms of accuracy and stability, is illustrated by numerical examples in 1D, 2D, and 3D.
This paper explores the restrictions imposed by bond-based peridynamics, particularly with respect to plane strain and plane stress models. We begin with a review of the derivations in Gerstle et al. (2005) wherein for isotropic materials a Poisson’s ratio restriction of $\frac{1}{4}$ for plane strain and $\frac{1}{3}$ for plane stress is deduced. Next, we show Cauchy’s relations are an intrinsic limitation of bond-based peridynamics and specialize these restrictions to plane strain and plane stress models. This generalizes the results from Gerstle et al. (2005) and demonstrates that the Poisson’s ratio restrictions described in Gerstle et al. (2005) are merely a consequence of Cauchy’s relations for isotropic materials. We conclude with a discussion of the validity of peridynamic plane strain and plane stress models formulated from two-dimensional bond-based peridynamic models.
Herein a general state-based peridynamics model is developed to simulate transport of fluids in an arbitrary heterogeneous porous medium. The generality encompasses modeling of multiphase, multi-component flow of non-Newtonian and compressible fluids, which is often encountered in but not limited to subsurface reservoirs. Peridynamic model is especially useful for solving non-local problems, such as crack propagation, since it does not assume spatial continuity of field variables. Thus, the formulation presented here, combined with peridynamics-based damage model, can be used to simulate hydraulic fracturing with complex fluids. To demonstrate its capability to simulate multi-phase flow in porous media, the derived model is verified against the analytical Buckley-Leverett solution for immiscible Newtonian two-phase flow. Further, the non-Newtonian two-phase fluid flow in porous media is verified by simulating the polymer flood process involving immiscible displacement of a Newtonian fluid by a non-Newtonian fluid against a generalized solution obtained by Wu et al. The non-local solutions are shown to be consistent with the corresponding local solutions in limiting cases. Moreover, mass conservation of all the phases is satisfied, irrespective of discretization and extent of non-locality.
Abstract not provided.
The peridynamic theory of solid mechanics is a nonlocal reformulation of the classical continuum mechanics theory. At the continuum level, it has been demonstrated that classical (local) elasticity is a special case of peridynamics. Such a connection between these theories has not been extensively explored at the discrete level. This paper investigates the consistency between nearest-neighbor discretizations of linear elastic peridynamic models and finite difference discretizations of the Navier–Cauchy equation of classical elasticity. While nearest-neighbor discretizations in peridynamics have been numerically observed to present grid-dependent crack paths or spurious microcracks, this paper focuses on a different, analytical aspect of such discretizations. We demonstrate that, even in the absence of cracks, such discretizations may be problematic unless a proper selection of weights is used. Specifically, we demonstrate that using the standard meshfree approach in peridynamics, nearest-neighbor discretizations do not reduce, in general, to discretizations of corresponding classical models. We study nodal-based quadratures for the discretization of peridynamic models, and we derive quadrature weights that result in consistency between nearest-neighbor discretizations of peridynamic models and discretized classical models. The quadrature weights that lead to such consistency are, however, model-/discretization-dependent. We motivate the choice of those quadrature weights through a quadratic approximation of displacement fields. The stability of nearest-neighbor peridynamic schemes is demonstrated through a Fourier mode analysis. Finally, an approach based on a normalization of peridynamic constitutive constants at the discrete level is explored. This approach results in the desired consistency for one-dimensional models, but does not work in higher dimensions. The results of the work presented in this paper suggest that even though nearest-neighbor discretizations should be avoided in peridynamic simulations involving cracks, such discretizations are viable, for example for verification or validation purposes, in problems characterized by smooth deformations. Furthermore, we demonstrate that better quadrature rules in peridynamics can be obtained based on the functional form of solutions.
Here, a notion of material homogeneity is proposed for peridynamic bodies with variable horizon but constant bulk properties. A relation is derived that scales the force state according to the position-dependent horizon while keeping the bulk properties unchanged. Using this scaling relation, if the horizon depends on position, artifacts called ghost forces may arise in a body under a homogeneous deformation. These artifacts depend on the second derivative of the horizon and can be reduced by employing a modified equilibrium equation using a new quantity called the partial stress . Bodies with piecewise constant horizon can be modeled without ghost forces by using a simpler technique called a splice . As a limiting case of zero horizon, both the partial stress and splice techniques can be used to achieve local-nonlocal coupling. Computational examples, including dynamic fracture in a one-dimensional model with local-nonlocal coupling, illustrate the methods.
Abstract not provided.
A notion of material homogeneity is proposed for peridynamic bodies with vari- able horizon but constant bulk properties. A relation is derived that scales the force state according to the position-dependent horizon while keeping the bulk properties un- changed. Using this scaling relation, if the horizon depends on position, artifacts called ghost forces may arise in a body under homogeneous deformation. These artifacts de- pend on the second derivative of horizon and can be reduced by use of a modified equilibrium equation using a new quantity called the partial stress . Bodies with piece- wise constant horizon can be modeled without ghost forces by using a technique called a splice between the regions. As a limiting case of zero horizon, both partial stress and splice techniques can be used to achieve local-nonlocal coupling. Computational examples, including dynamic fracture in a one-dimensional model with local-nonlocal coupling, illustrate the methods.