Engineering Papers⌕ Search

Engineering topics

Seleson, Pablo

Publications and source records attributed to Seleson, Pablo.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Meshfree Methods

Meshfree methods have undergone substantial development and have received much attention in the last two decades. This new family of numerical methods is designed to inherit the main advantages of the finite element method such as compact supports of shape functions and good approximation properties while, at the same time, overcome the main disadvantages of the finite element method caused by the mesh dependence. The meshfree methods share a common feature that no mesh is needed and shape functions are constructed from sets of points, thus eliminating the need for time consuming mesh generation. The most significant advantage of meshfree methods is the flexibility in customizing approximation functions for desired regularity and for capturing essential physics and features of the particular problems of interest. Adaptivity formulation and multiple-scale solution strategies also can be implemented with relative ease. It has become clear that the meshfree methods provide considerable advantages over the conventional finite element methods in solving problems involving moving discontinuities, evolving material interfaces, multiple-scale phenomena, large material distortion and structural deformation, and fracture and damage processes. This Chapter gives an overview of many classes of meshfree methods, with more detailed discussions on Smoothed Particle Hydrodynamics (SPH), the Reproducing Kernel Particle Method (RKPM), Peridynamics (PD), the Material Point Method (MPM), as well as their applications in various challenging engineering problems.2

Chen, Jiun-Shyan↗
DOE OSTI2024

Co-design for Particle Applications at Exascale

Co-design across the Exascale Computing Project (ECP) has been critical for both enabling science applications and bringing disparate communities together. Developing and porting applications to the various high-performance computing (HPC) architectures on pre-exascale and exascale computers has been quite challenging due to the diversity of hardware features and software stacks. The Co-design Center for Particle Applications (CoPA) has developed and enhanced the Cabana and PROGRESS/BML libraries to facilitate the creation of new particle applications, make existing particle applications exascale capable, and allow teams to explore new capabilities. Particle methods from atomistic, mesoscale, continuum, through cosmological scales have been built with Cabana, along with new possibilities for application coupling. Similarly, the PROGRESS/BML library has enabled quantum particle applications with linear algebra solvers to use advanced hardware. Across these CoPA-developed libraries, the co-design abstraction layer combines performance portability with math library support to facilitate separation of concerns and directly support science runs.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

PDMATLAB2D: A Peridynamics MATLAB Two-dimensional Code

PDMATLAB2D is a meshfree peridynamics implementation in MATLAB suitable for simulation of two-dimensional fracture problems. The purpose of this code is twofold. First, it provides an entry-level peridynamics computational tool for educational and training purposes. Second, it serves as an accessible and easily modifiable computational tool for peridynamics researchers who would like to adapt the code for a multitude of peridynamics simulation scenarios. The current version of the code implements a bond-based brittle elastic peridynamic model and a critical stretch criterion for bond breaking. However, the code is designed to be extendable for other peridynamic models and computational features. In this paper, we provide an overview of the code structure and functions with illustrative examples. In conclusion, due to the integrated computation and postprocessing MATLAB capabilities, PDMATLAB2D can serve as an effective testbed for testing new constitutive models and advanced numerical features for peridynamics computations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Dirichlet-type absorbing boundary conditions for peridynamic scalar waves in two-dimensional viscous media

Construction of absorbing boundary conditions (ABCs) for nonlocal models is generally challenging, primarily due to the fact that nonlocal operators are commonly associated with volume constrained boundary conditions. Moreover, application of Fourier and Laplace transforms, which are essential for the majority of available methods for ABCs, to nonlocal models is complicated. In this paper, we propose a simple method to construct accurate ABCs for peridynamic scalar wave-type problems in viscous media. The proposed ABCs are constructed in the time and space domains and are of Dirichlet type. Consequently, their implementation is relatively simple, since no derivatives of the wave field are required. The proposed ABCs are derived at the continuum level, from a semi-analytical solution of the exterior domain using harmonic exponential basis functions in space and time (plane-wave modes). The numerical implementation is done using a meshfree collocation approach employed within a boundary layer adjacent to the interior domain boundary. The modes satisfy the peridynamic numerical dispersion relation, resulting in a compatible solution of the interior region (near-field) with that of the exterior region (far-field). The accuracy and stability of the proposed ABCs are demonstrated with several numerical examples in two-dimensional unbounded domains.

42 ENGINEERING↗
DOE OSTI2023

Peridynamic elastic waves in two-dimensional unbounded domains: Construction of nonlocal Dirichlet-type absorbing boundary conditions

The focus of this paper is on application of peridynamics (PD) to propagation of elastic waves in unbounded domains. We construct absorbing boundary conditions (ABCs) derived from a semi-analytical solution of the PD governing equation at the exterior region. This solution is made up of a finite series of plane waves, as fundamental solutions (modes), which satisfy the PD dispersion relations. The modes are adjusted to transmit the energy from the interior region (near field) to the exterior region (far field). The corresponding unknown coefficients of the series are found in terms of the displacement field at a layer of points adjacent to the absorbing boundary. This is accomplished through a collocation procedure at subregions (clouds) around each absorbing point. The proposed ABCs offer appealing advantages, which facilitate their application to PD. They are of Dirichlet-type, hence their implementation is relatively simple as no derivatives of the field variables are required. They are constructed in the time and space domains and thus application of Fourier and Laplace transforms, cumbersome for nonlocal models, is not required. At the discrete level, the modes satisfy the same numerical dispersion relations of the near field, which makes the far-field solution compatible with that of the near field. We scrutinize the performance of the proposed ABCs through several examples. So our investigation shows that the proposed ABCs perform stably in time with an appropriate level of accuracy even in problems characterized by highly-dispersive propagating waves, including crack propagation in semi-unbounded brittle solids.

42 ENGINEERING↗
DOE OSTI2022

PDMATLAB2D

The PDMATLAB2D code is a meshfree peridynamics implementation in MATLAB suitable for simulation of two-dimensional fracture problems. The initial version implements bond-based brittle elastic peridynamics models and a critical stretch criterion for bond breaking. PDMATLAB2D provides an entry-level peridynamics computational tool for educational and training purposes. It also serves as an accessible and easily modifiable computational tool for peridynamics researchers who would like to adapt the code for a multitude of peridynamics simulations.

Seleson, Pablo↗
DOE OSTI2022

CabanaPD

CabanaPD is a meshfree peridynamics application built with Cabana and Kokkos. Kokkos enables performance portability across hardware architectures and Cabana provides particle capabilities including multi-node MPI support. The main components of CabanaPD are particle initialization, neighbor list generation, force and energy computation, time integration, and multi-node particle communication. In addition, options for creating pre-cracked regions and particle boundary conditions are available. CabanaPD currently enables two common force models: prototype microelastic brittle (PMB) and linear peridynamic solid (LPS). For both of these models, versions with and without fracture as well as linearized model options are available. CabanaPD is designed to be extensible for addition of more complex force models, boundary conditions, etc.

Reeve, Sam↗
DOE OSTI2022

A New Proof That the Number of Linear Elastic Symmetries in Two Dimensions Is Four

In this work, we present an elementary and self-contained proof that there are exactly four symmetry classes of the elasticity tensor in two dimensions: oblique, rectangular, square, and isotropic. In two dimensions, orthogonal transformations are either reflections or rotations. The proof is based on identification of constraints imposed by reflections and rotations on the elasticity tensor, and it simply employs elementary tools from trigonometry, making the proof accessible to a broad audience. For completeness, we identify the sets of transformations (rotations and reflections) for each symmetry class and report the corresponding equations of motions in classical linear elasticity.

42 ENGINEERING↗
DOE OSTI2022

Data-driven surrogate modeling of hPIC ion energy-angle distributions for high-dimensional sensitivity analysis of plasma parameters' uncertainty

In this work, we present a data-driven strategy for effective construction of a surrogate model in high-dimensional parameter space for the ion energy-angle distribution (IEAD) output of hPIC simulations of plasma-surface interactions. The methodology is based on a bin-by-bin least-squares fitting of the IEAD in the parameter space. The fitting is performed in a transformed coordinate system to normalize the IEAD, and it employs sparse grids for sampling the parameter space to overcome sampling challenges in high dimensions. The surrogate model is significantly cheaper computationally than direct hPIC simulations yet maintains high fidelity to them, providing a fast emulator for hPIC simulations. Sensitivity analysis based on the surrogate model is utilized to characterize the dependence of the ion impact angle and energy moments on the physical parameters.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Particle-in-cell (hPIC) data

Data generated from the hPIC particle-in-cell code to construct a surrogate ion energy-angle distribution (IEAD) model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

An Advanced Meso-Scale Peridynamic Modeling Technology using High-Performance Computing for Cost-Effective Product Design and Testing of Carbon Fiber Reinforced Polymer Composites in Light-weight Vehicles

We study a peridynamic composite modeling technology based on the discontinuous Galerkin finite element method, implemented in the commercial LS-DYNA software, for modeling and prediction of failure in carbon fiber reinforced polymer composites. The proposed technology is developed for the material failure analysis at the meso-scale, which provides the prevailing fiber-matrix interaction mechanism, without adoption of the representative volume element method and thus avoiding complicated numerical calibration procedures. Three types of experimental tests—in-plane coupon test, out-of-plane coupon test, and a component crash test—are simulated in a high-performance computing environment to assess the performance of the proposed peridynamic composite modeling technology.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A hybrid meshfree discretization to improve the numerical performance of peridynamic models

Efficient and accurate calculation of spatial integrals is of major interest in the numerical implementation of peridynamics (PD). The standard way to perform this calculation is a particle-based approach that discretizes the strong form of the PD governing equation. This approach has rapidly been adopted by the PD community since it offers some advantages. Additionally, it is computationally cheaper than other available schemes, can conveniently handle material separation, and effectively deals with nonlinear PD models. Nevertheless, PD models are still computationally very expensive compared with those based on the classical continuum mechanics theory, particularly for large-scale problems in three dimensions. This results from the nonlocal nature of the PD theory which leads to interactions of each node of a discretized body with multiple surrounding nodes. Here, we propose a new approach to significantly boost the numerical efficiency of PD models. We propose a discretization scheme that employs a simple collocation procedure and is truly meshfree; i.e., it does not depend on any background integration cells. In contrast to the standard scheme, the proposed scheme requires a much smaller set of neighboring nodes (keeping the same physical length scale) to achieve a specific accuracy and is thus computationally more efficient. Our new scheme is applicable to the case of linear PD models and within neighborhoods where the solution can be approximated by smooth basis functions. Therefore, to fully exploit the advantages of both the standard and the proposed schemes, a hybrid discretization is presented that combines both approaches within an adaptive framework. The high performance of the developed framework is illustrated by several numerical examples, including brittle fracture and corrosion problems in two and three dimensions.

42 ENGINEERING↗
DOE OSTI2021

A Review of Local-to-Nonlocal Coupling Methods in Nonlocal Diffusion and Nonlocal Mechanics

Local-to-nonlocal (LtN) coupling refers to a class of methods aimed at combining nonlocal and local modeling descriptions of a given system into a unified coupled representation. This allows to consolidate the accuracy of nonlocal models with the computational expediency of their local counterparts, while often simultaneously removing nonlocal modeling issues such as surface effects. The number and variety of proposed LtN coupling approaches have significantly grown in recent years, yet the field of LtN coupling continues to grow and still has open challenges. This review provides an overview of the state of the art of LtN coupling in the context of nonlocal diffusion and nonlocal mechanics, specifically peridynamics. Furthermore, we present a classification of LtN coupling methods and discuss common features and challenges. The goal of this review is not to provide a preferred way to address LtN coupling but to present a broad perspective of the field, which would serve as guidance for practitioners in the selection of appropriate LtN coupling methods based on the characteristics and needs of the problem under consideration.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

An FEM-Based Peridynamic Model for Failure Analysis of Unidirectional Fiber-Reinforced Laminates

To predict the mixed damage modes of unidirectional fiber-reinforced polymer (FRP) laminates under dynamic loading, an FEM-based peridynamic model is introduced in this paper. Based on its geometric structure and material composition, a long fiber lamina is considered a transversely isotropic medium as a result of homogenization at the meso-scale. The laminated structure is modeled by stacking surface mesh layers with arbitrary fiber angles along the thickness direction. The peridynamic bonds between Gauss points connect the separated elements. These bonds are classified as inner-layer bonds and inter-layer bonds. To represent the anisotropy of a laminate, the micro-elastic modulus of the inner-layer and inter-layer bonds is calculated from the anisotropic engineering material constants separately. To capture complex failure behaviors of laminate structures, an empirical damage model is proposed for the tension/ compression breakage of peridynamic bonds. This damage model can control the in-plane and delamination failure process. Finally, benchmark tests are conducted to validate the elastic response of laminates under dynamic loading. In terms of damage analysis, the proposed model can capture the complex damage modes and resistive force of laminate structures.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Overall equilibrium in the coupling of peridynamics and classical continuum mechanics

Coupling peridynamics based computational tools with those using classical continuum mechanics can be very beneficial, because it can provide a means to generate a computational method that combines the efficiency of classical continuum mechanics with the capability to simulate crack propagation, typical of peridynamics. This paper presents an overlooked issue in this type of coupled computational methods: the lack of overall equilibrium. This can be the case even if the coupling strategy satisfies the usual numerical tests involving rigid body motions as well as uniform and linear strain distributions. We focus our investigation on the lack of overall equilibrium in an approach to couple peridynamics and classical continuum mechanics recently proposed by the authors. In our examples, the magnitude of the out-of-balance forces is a fraction of a percent of the applied forces, but it cannot be assumed to be a numerical round-off error. We show analytically and numerically that the main reason for the existence of out-of-balance forces is a lack of balance between the local and nonlocal tractions at the coupling interface. This usually results from the presence of high-order derivatives of displacements in the coupling zone.

97 MATHEMATICS AND COMPUTING↗