🔬 materials science

Solving the Dissipation Inequality not as a constitutive restriction

This paper formulates and solves a constrained optimization procedure that treats the nonlinear Dissipation Inequality as a constraint equation rather than a constitutive restriction, enabling the automatic correction of faulty material specifications to ensure compliance with continuum mechanics postulates, as demonstrated through the elastoplastic response of a rate-dependent bar.

Original authors: Maximiliano Larrain Silva, Amit Acharya

Published 2026-08-04

📖 4 min read☕ Coffee break read

CC BY 4.0

Original authors: Maximiliano Larrain Silva, Amit Acharya

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine you are a chef trying to bake the perfect cake. You have a recipe (the laws of physics) that says, "You must never use more sugar than the bowl can hold, or the cake will collapse." But sometimes, you might accidentally grab a handful of sugar that is way too big. In the world of materials science, scientists try to predict how things like metal bars or rubber bands stretch and squish. They use "constitutive equations," which are fancy recipes describing how a material behaves. However, there is one unbreakable rule in the universe: the Second Law of Thermodynamics. In simple terms, this law says that whenever you deform a material, it must always "pay a tax" in the form of energy loss (called dissipation). You can't get energy back out of nowhere; the material must always lose at least a tiny bit of energy as heat. If a scientist's recipe accidentally predicts that a material gains energy while being stretched, the recipe is broken, and the physics doesn't make sense. The big question is: what do you do when your best guess for a material's behavior breaks this fundamental rule? Do you throw the whole recipe away, or is there a way to fix it just enough to make it work?

This paper by Maximiliano Larrain Silva and Amit Acharya tackles exactly that problem. They propose a clever new way to treat the "energy tax" rule not as a strict limit on how we write our recipes, but as a safety net that automatically corrects our mistakes. They set up a test case with a metal bar that is being stretched. They intentionally wrote a "bad" recipe for how the metal's internal structure (plastic strain) changes over time—a recipe that, if followed blindly, would violate the Second Law by suggesting the bar creates energy out of thin air during certain moments. Instead of discarding this flawed recipe, their method introduces a "correction variable," which acts like a tiny, invisible hand that nudges the material's behavior just enough to keep the energy tax positive.

Think of it like driving a car with a slightly broken GPS. The GPS (the material's constitutive equation) tells you to turn left into a lake (violating the laws of physics). The driver (the scientist's model) knows this is impossible. Usually, you'd have to rewrite the GPS map entirely. But this paper suggests a different approach: keep the GPS, but add a "correction steering wheel" that automatically turns the car slightly right whenever the GPS tries to drive you into the lake. The goal is to make the smallest possible turn to stay safe. The authors found that when they applied this method, the "correction steering wheel" kicked in exactly when the bad recipe tried to break the rules. It forced the material to stop deforming plastically for a moment, creating what they call an "elastic gap"—a brief pause where the material acts purely like a spring, refusing to flow until the rules are safe again.

The paper demonstrates this using both mathematical formulas and computer simulations. They tested two types of materials: one that reacts quickly to how fast it's stretched (rate-sensitive) and one that reacts slowly (rate-insensitive). In both cases, their computer code successfully "fixed" the broken recipe. When the bad recipe tried to make the energy tax negative, the correction variable stepped in, driving the dissipation to exactly zero (the bare minimum allowed) rather than letting it go negative. The result was a stress-strain curve that looked almost identical to the "perfect" mathematical solution they calculated by hand, with errors as low as 0.05% in some cases. The only time the computer struggled slightly was when the bad recipe changed its mind very suddenly, causing a tiny spike in error, but even then, the fix held up.

The authors show that this approach doesn't just patch the problem; it finds the most logical solution among infinite possibilities. By trying to keep the correction as small as possible, the method naturally selects the solution that stays closest to the original, flawed recipe while still obeying the laws of physics. It's a bit like a self-correcting autopilot that doesn't just crash the plane when the navigation fails, but gently steers it back to a safe path, ensuring the flight continues without breaking the laws of aerodynamics. This work suggests that even when our understanding of a material's behavior is incomplete or slightly wrong, we can still get reliable predictions by letting the fundamental laws of thermodynamics act as a gentle, automatic referee.

Technical Summary: Solving the Dissipation Inequality as a Constraint Equation

Problem Statement
The paper addresses a fundamental challenge in continuum mechanics: the difficulty of accommodating the Second Law of Thermodynamics (the Dissipation Inequality) as a strict constraint when the constitutive behavior of a material is not fully known or is specified in a manner that may inadvertently violate thermodynamic principles. Specifically, the authors investigate a scenario where a prescribed constitutive equation for a rate-dependent, elastic-plastic bar is intentionally chosen to violate the Dissipation Inequality over a specific time interval. The goal is to develop a self-consistent procedure that accepts such a "faulty" specification as input but produces a solution that strictly adheres to the Second Law with minimal deviation from the prescribed material response.

Methodology
The authors formulate the problem as a nonlinear constrained optimization task. Instead of treating the Dissipation Inequality as a restriction on the choice of constitutive functions (as is traditional), they treat it as a constraint equation within the solution process.

Mathematical Formulation:
- The plastic strain rate evolution is decomposed into a prescribed part, $f_c(\sigma, t)$ , and a corrective field, $a(x, t)$ .
- The system seeks to determine the fields for displacement ( $u$ ), plastic strain ( $p$ ), the corrective field ( $a$ ), and a dissipation density field ( $s$ ).
- The objective is to minimize the integral of a functional $H(a, s) = \frac{1}{2}c_a a^2 + \frac{1}{2}c_s s^2$ over time. This minimization ensures that the corrective field $a$ and the dissipation $s^2$ are as small as possible, reflecting a "faith" in the prescribed constitutive model while strictly enforcing non-negative dissipation.
- The constraints include force equilibrium, the modified constitutive relation, and the Dissipation Inequality.
Analytical Solution:
- A closed-form solution is derived for the 1D bar problem. The optimization yields a specific selection of the corrective field $a(t)$ and dissipation $s(t)$ that satisfies the constraints.
- In the limit where the ratio of material constants approaches zero, the solution naturally selects a state where the plastic strain rate becomes zero whenever the prescribed evolution would violate the Second Law, effectively creating an "elastic gap."
Computational Scheme:
- The authors employ a Dual-to-Primal (DtP) variational formulation based on a convex dual variational principle.
- The method utilizes a sequence of convex optimization problems to solve the governing equations.
- The numerical algorithm proceeds in two phases:
  - Gradient Flow: A gradient descent scheme in a fictitious time variable drives the initial approach to the solution.
  - Newton-Raphson: A Newton-Raphson scheme with step-size control refines the solution.
- Crucially, the algorithm restricts the solution domain (the "DtP zone") to ensure the convexity of the dual functional, preventing the solution from violating the Second Law.
- The method uses finite element discretization in time for the dual fields.

Key Results
The scheme was tested on a 1D bar under quasi-static stretching with a prescribed plastic evolution function that becomes negative (violating the Second Law) during a specific time interval.

Correction of Constitutive Violation: The computational scheme successfully "corrected" the intentionally faulty constitutive specification. When the prescribed evolution $f_c$ threatened to produce negative dissipation, the corrective field $a(t)$ activated to ensure the total dissipation remained non-negative.
Elastic Gaps: The activation of the corrective field resulted in the plastic strain rate becoming zero during the violation interval. Consequently, the stress-strain response exhibited an "elastic gap" where the material behaved purely elastically ( $d\sigma/d\epsilon = E$ ) despite the increasing external load. This behavior emerged naturally from the optimization rather than being prescribed.
Accuracy: The numerical solutions (Finite Element Method) showed excellent agreement with the analytical closed-form solutions.
- For the rate-sensitive case ( $m=1$ ), the error in the stress-strain relation was below 0.05%.
- For the rate-insensitive case ( $m=0.1$ ), the error was below 0.8%.
- Errors were localized to regions where the modulating function $g(t)$ underwent rapid transitions, suggesting the errors stem from mesh resolution rather than the algorithmic scheme itself.
Role of Initial Conditions: The study highlighted that the initial base state in the numerical algorithm acts as a selection parameter. A small initial value for the dissipation field guided the numerical solution to the specific closed-form solution derived in the analytical section. Larger initial values converged to different, yet physically valid, solutions within the infinite family of solutions permitted by the constraints.

Significance and Claims
The paper claims to demonstrate a viable scheme for treating the Dissipation Inequality not as a pre-condition for constitutive modeling, but as a constraint to be satisfied during the solution of boundary value problems.

Robustness to Uncertainty: The approach allows for the use of constitutive models that may be incomplete or imperfect, automatically correcting them to satisfy fundamental thermodynamic postulates with "minimal" deviation.
Computational Feasibility: The authors provide a first numerical confirmation that the dual variational principle and the associated DtP scheme can enforce the Second Law in scenarios where the law might otherwise be violated by the specified process.
Application Scope: The authors suggest the technique is particularly relevant for complex mechanical models where Second Law restrictions on constitutive functions are intricate, or for coupling disparate phenomena where standard constitutive restrictions are difficult to formulate. They also note the emergence of "elastic gaps" in the stress-strain response, which bears resemblance to behaviors observed in strain-gradient plasticity theories, though this is presented as an emergent feature of the specific example rather than a new physical theory.

The work does not propose new experimental setups or claim to solve all constitutive modeling issues, but rather offers a specific mathematical framework to handle the constraint of thermodynamic admissibility in the presence of uncertain or potentially flawed material specifications.

Technical Summary: Solving the Dissipation Inequality as a Constraint Equation

More like this