BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Lars Pastewka

Lars Pastewka contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
cond-mat.mtrl-scimath.NANumerical AnalysisComputational Engineering, Finance, and Sciencecond-mat.softphysics.comp-phphysics.flu-dyn

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
14works
0followers
7topics
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

14 published item(s)

preprint2022arXiv

Crack-path selection in phase-field models for brittle fracture

This work presents a critical overview of the effects of different aspects of model formulation on crack path selection in quasi-static phase field fracture. We consider different evolution methods, mechanics formulations, fracture dissipation energy formulations, and forms of the irreversibility condition. The different model variants are implemented with common numerical methods based on staggered solution of the phase-field and mechanics sub-problems via FFT-based solvers. These methods mix standard approaches with novel elements, such as the use of bound-constrained conjugate gradients for the phase field sub-problem and a heuristic method for near-equilibrium evolution. We examine differences in crack paths between model variants in simple model systems and microstructures with randomly heterogeneous Young's modulus. Our results indicate that near-equilibrium evolution methods are preferable for quasi-static fracture of heterogeneous microstructures compared to minimization and time-dependent methods. In examining mechanics formulations, we find distinct effects of crack driving force and the model for contact implicit in phase field fracture. Our results favor the use of a strain-spectral decomposition for the crack driving force but not the contact model. Irreversibility condition and fracture dissipation energy formulation were also found to affect crack path selection, but systematic effects were difficult to deduce due to the overall sensitivity of crack selection within the heterogeneous microstructures. Our findings support the use of the AT1 model over the AT2 model and irreversibility of the phase field within a crack set rather than the entire domain. Sensitivity to these differences in formulation was reduced but not eliminated by reducing the crack width parameter $\ell$ relative to the size scale of the random microstructures.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Efficient topology optimization using compatibility projection in micromechanical homogenization

The adjoint method allows efficient calculation of the gradient with respect to the design variables of a topology optimization problem. This method is almost exclusively used in combination with traditional Finite-Element-Analysis, whereas Fourier-based solvers have recently shown large efficiency gains for homogenization problems. In this paper, we derive the discrete adjoint method for Fourier-based solvers that employ compatibility projection. We demonstrate the method on the optimization of composite materials and auxetic metamaterials, where void regions are modelled with zero stiffness.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Optimal FFT-accelerated Finite Element Solver for Homogenization

We propose a matrix-free finite element (FE) homogenization scheme that is considerably more efficient than generic FE implementations. The efficiency of our scheme follows from a preconditioned well-scaled reformulation allowing for the use of the conjugate gradient or similar iterative solvers. The geometrically-optimal preconditioner -- a discretized Green's function of a periodic homogeneous reference problem -- has a block-diagonal structure in the Fourier space which permits its efficient inversion using the fast Fourier transform (FFT) techniques for generic regular meshes. This implies that the scheme scales as $\mathcal{O}(n \log(n))$ like FFT, rendering it equivalent to spectral solvers in terms of computational efficiency. However, in contrast to classical spectral solvers, the proposed scheme works with FE shape functions with local supports and is free of the Fourier ringing phenomenon. We showcase that the scheme achieves the number of iterations that are almost independent of spatial discretisation and scales mildly with the phase contrast. Additionally, we discuss the equivalence between our displacement-based scheme and the recently proposed strain-based homogenization technique with finite-element projection.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The surface-topography challenge: Problem definition

We present to the community a surface-definition problem, whose solution we consider to be critical for the proper description of contacts between nominally flat surfaces [1,2]. In 2015, Müser and Dapp issued the Contact Mechanics Challenge, which provided complete topography data for a fictional surface and asked theorists and modelers to compute the expected contact parameters for such a surface. This effort was a success, but exposed one glaring flaw in the community's understanding of the nature of contact: these models require as input a complete description of surface topography, which is rarely or never available for real-world surfaces [3-6]. The present challenge is to experimentalists: we will send you samples of two materials (one smoother and one rougher); you determine the surface topography of these materials. We call on you to measure such surfaces however you wish, using contact-based techniques, light scattering, microscopy, or other techniques. Examples of quantities of interest are: root-mean-square (RMS) parameters; the power spectral density (PSD); or the autocorrelation function (ACF). For the material, we have chosen chromium nitride, a wear- and corrosion-resistant coating used in industrial applications including automotive components, cutting tools, and die-casting. To participate, simply go to: https://contact.engineering/challenge to provide your shipping address and other information, then samples will be shipped out to you. The only requirement of participation is that your raw topography measurements are deposited on the free contact.engineering web app to facilitate data sharing. The purpose of this challenge is for our community to move towards: (a) better agreement on how to describe the multi-scale topography of experimental surfaces; and (b) better understanding of how to apply the well-developed models and theories to real-world surfaces.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Crack-front model for adhesion of soft elastic spheres with chemical heterogeneity

Adhesion hysteresis can be caused by elastic instabilities that are triggered by surface roughness or chemical heterogeneity. However, the role of these instabilities in adhesion hysteresis remains poorly understood because we lack theoretical and numerical models accounting for realistic roughness. Our work focuses on the adhesion of soft elastic spheres with low roughness or weak heterogeneity, where the indentation process can be described as a Griffith-like propagation of a nearly circular external crack. We discuss how to describe the contact of spheres with chemical heterogeneity that leads to fluctuations in the local work of adhesion. We introduce a variational first-order crack-perturbation model and validate our approach using boundary-element simulations. The crack-perturbation model faithfully predicts contact shapes and hysteretic force-penetration curves, provided that the contact perimeter remains close to a circle and the contact area is simply connected. Computationally, the crack-perturbation model is orders of magnitude more efficient than the corresponding boundary element formulation, allowing for realistic heterogeneity fields. Furthermore, our crack-front formulation clarifies the connection of adhesion hysteresis to classic theories on pinning of elastic lines.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Elimination of ringing artifacts by finite-element projection in FFT-based homogenization

Micromechanical homogenization is often carried out with Fourier-accelerated methods that are prone to ringing artifacts. We here generalize the compatibility projection introduced by Vondřejc, Zeman & Marek [Comput. Math. Appl. 68, 156 (2014)] beyond the Fourier basis. In particular, we formulate the compatibility projection for linear finite elements while maintaining Fourier-acceleration and the fast convergence properties of the original method. We demonstrate that this eliminates ringing artifacts and yields an efficient computational homogenization scheme that is equivalent to canonical finite-element formulations on fully structured grids.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Height-averaged Navier-Stokes solver for hydrodynamic lubrication

Modelling hydrodynamic lubrication is crucial in the design of engineering components as well as for a fundamental understanding of friction mechanisms. The cornerstone of thin-film flow modelling is the Reynolds equation -- a lower-dimensional representation of the Stokes equation. However, the derivation of the Reynolds equation is based on assumptions and fixed form constitutive relations, that may not generally be valid, especially when studying systems under extreme conditions. Furthermore, these explicit assumptions about the constitutive behaviour of the fluid prohibit applications in a multiscale scenario based on measured or atomistically simulated data. Here, we present a method that considers the full compressible Navier-Stokes equation in a height-averaged sense for arbitrary constitutive relations. We perform numerical tests by using a reformulation of the viscous stress tensor for laminar flow to validate the presented method comparing to results from conventional Reynolds solutions. The versatility of the method is shown by incorporating models for mass-conserving cavitation, wall slip and non-Newtonian fluids. This allows testing of new constitutive relations that not necessarily need to take a fixed form, and may be obtained from experimental or simulation data.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Surface lattice Green's functions for high-entropy alloys

We study the surface elastic response of pure Ni, the random alloy FeNiCr and an average FeNiCr alloy in terms of the surface lattice Green's function. We propose a scheme for computing per-site Green's function and study their per-site variations. The average FeNiCr alloy accurately reproduces the mean Green's function of the full random alloy. Variation around this mean is largest near the edge of the surface Brillouin-zone and decays as $q^{-2}$ with wavevector $q$ towards the $Γ$-point. We also present expressions for the continuum surface Green's function of anisotropic solids of finite and infinite thickness and show that the atomistic Green's function approaches continuum near the $Γ$-point. Our results are a first step towards efficient contact calculations and Peierls-Nabarro type models for dislocations in high-entropy alloys.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Quantitative prediction of the fracture toughness of amorphous carbon from atomic-scale simulations

Fracture is the ultimate source of failure of amorphous carbon (a-C) films, however it is challenging to measure fracture properties of a-C from nano-indentation tests and results of reported experiments are not consistent. Here, we use atomic-scale simulations to make quantitative and mechanistic predictions on fracture of a-C. Systematic large-scale K-field controlled atomic-scale simulations of crack propagation are performed for a-C samples with densities of $ρ=2.5, \, 3.0 \, \text{ and } 3.5~\text{g/cm}^{3}$ created by liquid quenches for a range of quench rates $\dot{T}_q = 10 - 1000~\text{K/ps}$. The simulations show that the crack propagates by nucleation, growth, and coalescence of voids. Distances of $ \approx 1\, \text{nm}$ between nucleated voids result in a brittle-like fracture toughness. We use a crack growth criterion proposed by Drugan, Rice \& Sham to estimate steady-state fracture toughness based on our short crack-length fracture simulations. Fracture toughness values of $2.4-6.0\,\text{MPa}\sqrt{\text{m}}$ for initiation and $3-10\,\text{MPa}\sqrt{\text{m}}$ for the steady-state crack growth are within the experimentally reported range. These findings demonstrate that atomic-scale simulations can provide quantitatively predictive results even for fracture of materials with a ductile crack propagation mechanism.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Pattern formation during deformation of metallic nanolaminates

We used nonequilibrium molecular dynamics simulations to study the shear deformation of metallic composites composed of alternating layers of Cu and Au. Our simulations reveal the formation of "vortices" or "swirls" if the bimaterial interfaces are atomically rough and if none of the {111} planes that accommodate slip in fcc materials is exactly parallel to this interface. We trace the formation of these patterns back to grain rotation, induced by hindering dislocations from crossing the bimaterial interface. The instability is accompanied by shear-softening of the material. These calculations shed new light on recent observations of pattern formation in plastic flow, mechanical mixing of materials and the common formation of a tribomutation layer in tribologically loaded systems.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Structural and elastic properties of amorphous carbon from simulated quenching at low rates

We generate representative structural models of amorphous carbon (a-C) from constant-volume quenching from the liquid with subsequent relaxation of internal stresses in molecular dynamics simulations using empirical and machine-learning interatomic potentials. By varying volume and quench rate we generate structures with a range of density and amorphous morphologies. We find that all a-C samples show a universal relationship between hybridization, bulk modulus and density despite having distinct cohesive energies. Differences in cohesive energy are traced back to slight changes in the distribution of bond-angles that will likely affect thermal stability of these structures.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Surface flaws control strain localization in the deformation of Cu$\vert$Au nanolaminates

We carried out matched experiments and molecular dynamics simulations of the compression of nanopillars prepared from Cu$\vert$Au nanolaminates with 25 nm layer thickness. The stress-strain behavior obtained from both techniques are in excellent agreement. Variation of the layer thickness in simulations reveals an increase of the strength with decreasing layer thickness. Pillars fail through the formation of shear bands whose nucleation we trace back to the existence of surface flaws. Our combined approach demonstrates the crucial role of contact geometry in controlling the deformation mode and suggests that modulus-matched nanolaminates should be able to suppress strain localization while maintaining controllable strength.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

The universal emergence of self-affine roughness from deformation

Most natural and man-made surfaces appear to be rough on many length scales. There is presently no unifying theory of the origin of roughness or the self-affine nature of surface topography. One likely contributor to the formation of roughness is deformation, which underlies many processes that shape surfaces such as machining, fracture, and wear. Using molecular dynamics, we simulate the bi-axial compression of single-crystal Au, the high-entropy alloy Ni$_{36.67}$Co$_{30}$Fe$_{16.67}$Ti$_{16.67}$, and amorphous Cu$_{50}$Zr$_{50}$, and show that even surfaces of homogeneous materials develop a self-affine structure. By characterizing subsurface deformation, we connect the self-affinity of the surface to the spatial correlation of deformation events occurring within the bulk and present scaling relations for the evolution of roughness with strain.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Revised Periodic Boundary Conditions: Fundamentals, Electrostatics, and the Tight-Binding Approximation

Many nanostructures today are low-dimensional and flimsy, and therefore get easily distorted. Distortion-induced symmetry-breaking makes conventional, translation-periodic simulations invalid, which has triggered developments for new methods. Revised periodic boundary conditions (RPBC) is a simple method that enables simulations of complex material distortions, either classically or quantum-mechanically. The mathematical details of this easy-to-implement approach, however, have not been discussed before. Therefore, in this paper we summarize the underlying theory, present the practical details of RPBC, especially related to a non-orthogonal tight-binding formulation, discuss selected features, electrostatics in particular, and suggest some examples of usage. We hope this article to give more insight to RPBC, to help and inspire new software implementations capable of exploring the physics and chemistry of distorted nanomaterials.

0 citations0 reviewsNo code linkedOpen access