ExaCMech: GPU-Native Crystal Plasticity Constitutive Library

Tue, 15 Apr 2025 00:00:00 +0000

The Problem ExaCMech Solves
#

Crystal plasticity finite element codes spend a large fraction of their time doing one thing: the constitutive update. Given a material’s current state (its crystal orientation, internal hardening variables, elastic strain) and a prescribed deformation over a time step, compute the resulting stress and update the material state. This has to be done at every quadrature point in the mesh, which in a production micromechanics simulation means evaluating physically complex, iterative nonlinear equations simultaneously at tens of millions of points. For the problem to be tractable at scale, those evaluations need to run on the GPU, and the models need to be structured in a way that maps naturally to how GPUs actually execute work.

SNLS: A Small Nonlinear Solver Library That Punches Above Its Weight

Tue, 15 Apr 2025 00:00:00 +0000

Where It Started
#

SNLS predates both ExaCMech and ExaConstit. It was developed at LLNL specifically to solve the problem of running material model constitutive updates on the GPU, at a time when no existing nonlinear solver library was designed to do that. General-purpose solvers like MINPACK or PETSc are built for problems that live on the host. They carry significant overhead per call and are not structured for batched execution across millions of independent points simultaneously. In crystal plasticity specifically, every quadrature point in the mesh needs its own nonlinear solve to update stress and internal state at each time step, and those systems are dense, stiff, and need to run as fast as possible without any per-point failure being an option. That was the gap SNLS was built to fill.

Raja on Robert Carson

ExaCMech: GPU-Native Crystal Plasticity Constitutive Library

The Problem ExaCMech Solves #

SNLS: A Small Nonlinear Solver Library That Punches Above Its Weight

Where It Started #

The Problem ExaCMech Solves
#

Where It Started
#