BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Lars Pastewka

Lars Pastewka appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.mtrl-scicond-mat.softmath.NANumerical AnalysisComputational Engineering, Finance, and Sciencephysics.comp-phphysics.flu-dyn

Catalog footprint

What is connected

21works
7topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

21 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
preprint2016arXiv

Contact area of rough spheres: Large scale simulations and simple scaling laws

We use molecular simulations to study the nonadhesive and adhesive atomic-scale contact of rough spheres with radii ranging from nanometers to micrometers over more than ten orders of magnitude in applied normal load. At the lowest loads, the interfacial mechanics is governed by the contact mechanics of the first asperity that touches. The dependence of contact area on normal force becomes linear at intermediate loads and crosses over to Hertzian at the largest loads. By combining theories for the limiting cases of nominally flat rough surfaces and smooth spheres, we provide parameter-free analytical expressions for contact area over the whole range of loads. Our results establish a range of validity for common approximations that neglect curvature or roughness in modeling objects on scales from atomic force microscope tips to ball bearings.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Offset-corrected $Δ$-Kohn-Sham scheme for the prediction of X-ray photoelectron spectra of molecules and solids

Absolute binding energies of core electrons in molecules and bulk materials can be efficiently calculated by spin paired density-function theory employing a $Δ$ Kohn-Sham ($Δ$KS) scheme corrected by offsets that are highly transferable. These offsets depend on core level and atomic species and can be determined by comparing $Δ$KS energies to experimental molecular X-ray photoelectron spectra. We demonstrate the correct prediction of absolute and relative binding energies on a wide range of molecules, metals and insulators.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Finite-size scaling in the interfacial stiffness of rough elastic contacts

The total elastic stiffness of two contacting bodies with a microscopically rough interface has an interfacial contribution K that is entirely attributable to surface roughness. A quantitative understanding of K is important because it can dominate the total mechanical response and because it is proportional to the interfacial contributions to electrical and thermal conductivity in continuum theory. Numerical simulations of the dependence of K on the applied squeezing pressure p are presented for nominally flat elastic solids with a range of surface roughnesses. Over a wide range of p, K rises linearly with p. Sublinear power-law scaling is observed at small p, but the simulations reveal that this is a finite-size effect. We derive accurate, analytical expressions for the exponents and prefactors of this low-pressure scaling of K by extending the contact mechanics theory of Persson to systems of finite size. In agreement with our simulations, these expressions show that the onset of the low-pressure scaling regime moves to lower pressure as the system size increases.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

On the validity of the method of reduction of dimensionality: area of contact, average interfacial separation and contact stiffness

It has recently been suggested that many contact mechanics problems between solids can be accurately studied by mapping the problem on an effective one dimensional (1D) elastic foundation model. Using this 1D mapping we calculate the contact area and the average interfacial separation between elastic solids with nominally flat but randomly rough surfaces. We show, by comparison to exact numerical results, that the 1D mapping method fails even qualitatively. We also calculate the normal interfacial stiffness $K$ and compare it with the result of an analytical study. We attribute the failure of the elastic foundation model to the neglect of the long-range elastic coupling between the asperity contact regions.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Screened empirical bond-order potentials for Si-C

Typical empirical bond-order potentials are short ranged and give ductile instead of brittle behavior for materials such as crystalline silicon or diamond. Screening functions can be used to increase the range of these potentials. We outline a general procedure to combine screening functions with bond-order potentials that does not require to refit any of the potential's properties. We use this approach to modify Tersoff's [Phys. Rev. B 39, 5566 (1989)], Erhart & Albe's [Phys. Rev. B 71, 35211 (2005)] and Kumagai et al.'s [Comp. Mater. Sci. 39, 457 (2007)] Si, C and Si-C potentials. The resulting potential formulations correctly reproduce brittle materials response, and give an improved description of amorphous phases.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

When are rough surfaces sticky?

At the molecular scale there are strong attractive interactions between surfaces, yet few macroscopic surfaces are sticky. Extensive simulations of contact by adhesive surfaces with roughness on nanometer to micrometer scales are used to determine how roughness reduces the area where atoms contact and thus weakens adhesion. The material properties, adhesive strength and roughness parameters are varied by orders of magnitude. In all cases the area of atomic contact rises linearly with load, and the prefactor rises linearly with adhesive strength for weak interactions. Above a threshold adhesive strength, the prefactor changes sign, the surfaces become sticky and a finite force is required to separate them. A parameter-free analytic theory is presented that describes changes in these numerical results over up to five orders of magnitude in load. It relates the threshold strength to roughness and material properties, explaining why most macroscopic surfaces do not stick. The numerical results are qualitatively and quantitatively inconsistent with classical theories based on the Greenwood-Williamson approach that neglect the range of adhesion and do not include asperity interactions.

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