Researcher profile

Peter Wriggers

Peter Wriggers contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
6works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

6 published item(s)

preprint2022arXiv

Parallel Multiphysics Simulation for the Stabilized Optimal Transportation Meshfree (OTM)

This paper presents a parallel \PG{implementation} for the Optimal Transportation Meshfree (OTM) method on large CPU clusters. Communications are handled with the Message Passing Interface (MPI). The Recursive Coordinate Bisection (RCB) algorithm is utilized for domain decomposition and for implementing dynamic load-balancing strategy. This work involves three new concepts to reduce the computational efforts: Dynamic halo regions, Efficient data management strategies for ease of addition and deletion of nodes and material points using advanced STL container, and nearest neighborhood communication for detection of neighbors and communication. Also, Linked Cell approach has been implemented to further reduce the computational efforts. Parallel performance analysis is investigated for challenging multiphysics applications like Taylor rod impact and serrated chip formation process. Adequate scalability of parallel implementation for these applications is reported.

preprint2021arXiv

A matrix-free isogeometric Galerkin method for Karhunen-Loève approximation of random fields using tensor product splines, tensor contraction and interpolation based quadrature

The Karhunen-Loève series expansion (KLE) decomposes a stochastic process into an infinite series of pairwise uncorrelated random variables and pairwise $L^2$-orthogonal functions. For any given truncation order of the infinite series the basis is optimal in the sense that the total mean squared error is minimized. The orthogonal basis functions are determined as the solution of an eigenvalue problem corresponding to the homogeneous Fredholm integral equation of the second kind, which is computationally challenging for several reasons. Firstly, a Galerkin discretization requires numerical integration over a $2d$ dimensional domain, where $d$, in this work, denotes the spatial dimension. Secondly, the main system matrix of the discretized weak-form is dense. Consequently, the computational complexity of classical finite element formation and assembly procedures as well as the memory requirements of direct solution techniques become quickly computationally intractable with increasing polynomial degree, number of elements and degrees of freedom. The objective of this work is to significantly reduce several of the computational bottlenecks associated with numerical solution of the KLE. We present a matrix-free solution strategy, which is embarrassingly parallel and scales favorably with problem size and polynomial degree. Our approach is based on (1) an interpolation based quadrature that minimizes the required number of quadrature points; (2) an inexpensive reformulation of the generalized eigenvalue problem into a standard eigenvalue problem; and (3) a matrix-free and parallel matrix-vector product for iterative eigenvalue solvers. Two higher-order three-dimensional benchmarks illustrate exceptional computational performance combined with high accuracy and robustness.

preprint2020arXiv

A global-local approach for hydraulic phase-field fracture in poroelastic media

In this work, phase-field modeling of hydraulic fractures in porous media is extended towards a global-local approach. Therein, the failure behavior is solely analyzed in a (small) local domain. In the surrounding medium, a simplified and linearized system of equations is solved. Both domains are coupled by Robin-type interface conditions. The fracture(s) inside the local domain are allowed to propagate and consequently both subdomains change within time. Here, a predictor-corrector strategy is adopted in which the local domain is dynamically adjusted to the current fracture pattern. The resulting framework is algorithmically described in detail and substantiated with some numerical tests.

preprint2020arXiv

A machine learning based plasticity model using proper orthogonal decomposition

Data-driven material models have many advantages over classical numerical approaches, such as the direct utilization of experimental data and the possibility to improve performance of predictions when additional data is available. One approach to develop a data-driven material model is to use machine learning tools. These can be trained offline to fit an observed material behaviour and then be applied in online applications. However, learning and predicting history dependent material models, such as plasticity, is still challenging. In this work, a machine learning based material modelling framework is proposed for both elasticity and plasticity. The machine learning based hyperelasticity model is developed with the Feed forward Neural Network (FNN) directly whereas the machine learning based plasticity model is developed by using of a novel method called Proper Orthogonal Decomposition Feed forward Neural Network (PODFNN). In order to account for the loading history, the accumulated absolute strain is proposed to be the history variable of the plasticity model. Additionally, the strain-stress sequence data for plasticity is collected from different loading-unloading paths based on the concept of sequence for plasticity. By means of the POD, the multi-dimensional stress sequence is decoupled leading to independent one dimensional coefficient sequences. In this case, the neural network with multiple output is replaced by multiple independent neural networks each possessing a one-dimensional output, which leads to less training time and better training performance. To apply the machine learning based material model in finite element analysis, the tangent matrix is derived by the automatic symbolic differentiation tool AceGen. The effectiveness and generalization of the presented models are investigated by a series of numerical examples using both 2D and 3D finite element analysis.

preprint2020arXiv

The Neural Particle Method -- An Updated Lagrangian Physics Informed Neural Network for Computational Fluid Dynamics

Numerical simulation is indispensable in industrial design processes. It can replace expensive experiments and even reduce the need for prototypes. While products designed with the aid of numerical simulation undergo continuous improvement, this must also be true for numerical simulation itself. Up to date, no general purpose numerical method is available which can accurately resolve a variety of physics ranging from fluid to solid mechanics including large deformations and free surface flow phenomena. These complex multi-physics problems occur for example in Additive Manufacturing processes. In this sense, the recent developments in Machine Learning display promise for numerical simulation. It has recently been shown that instead of solving a system of equations as in standard numerical methods, a neural network can be trained solely based on initial and boundary conditions. Neural networks are smooth, differentiable functions that can be used as a global ansatz for Partial Differential Equations (PDEs). While this idea dates back to more than 20 years ago [Lagaris et al., 1998], it is only recently that an approach for the solution of time dependent problems has been developed [Raissi et al., 2019]. With the latter, implicit Runge Kutta schemes with unprecedented high order have been constructed to solve scalar-valued PDEs. We build on the aforementioned work in order to develop an Updated Lagrangian method for the solution of incompressible free surface flow subject to the inviscid Euler equations. The method is easy to implement and does not require any specific algorithmic treatment which is usually necessary to accurately resolve the incompressibility constraint. Due to its meshfree character, we will name it the Neural Particle Method (NPM). It will be demonstrated that the NPM remains stable and accurate even if the location of discretization points is highly irregular.

preprint2019arXiv

An Adaptive Global-Local Approach for Phase-Field Modeling of Anisotropic Brittle Fracture

This work addresses an efficient Global-Local approach supplemented with predictor-corrector adaptivity applied to anisotropic phase-field brittle fracture. The phase-field formulation is used to resolve the sharp crack surface topology on the anisotropic/non-uniform local state in the regularized concept. To resolve the crack phase-field by a given single preferred direction, second-order structural tensors are imposed to both the bulk and crack surface density functions. Accordingly, a split in tension and compression modes in anisotropic materials is considered. A Global-Local formulation is proposed, in which the full displacement/phase-field problem is solved on a lower (local) scale, while dealing with a purely linear elastic problem on an upper (global) scale. Robin-type boundary conditions are introduced to relax the stiff local response at the global scale and enhancing its stabilization. Another important aspect of this contribution is the development of an adaptive Global-Local approach, where a predictor-corrector scheme is designed in which the local domains are dynamically updated during the computation. To cope with different finite element discretizations at the interface between the two nested scales, a non-matching dual mortar method is formulated. Hence, more regularity is achieved on the interface. Several numerical results substantiate our developments.