-
Residual stresses couple microscopic and macroscopic scales
Authors:
Sebastian Steinhäuser,
Timm Treskatis,
Stefan Turek,
Thomas Voigtmann
Abstract:
We show how residual stresses emerge in a visco-elastic material as a signature of its past flow history, through an interplay between flow-modified microscopic relaxation and macroscopic features of the flow. Long-lasting temporal-history dependence of the microscopic dynamics and nonlinear rheology are incorporated through the mode-coupling theory of the glass transition (MCT). The theory's inte…
▽ More
We show how residual stresses emerge in a visco-elastic material as a signature of its past flow history, through an interplay between flow-modified microscopic relaxation and macroscopic features of the flow. Long-lasting temporal-history dependence of the microscopic dynamics and nonlinear rheology are incorporated through the mode-coupling theory of the glass transition (MCT). The theory's integral constitutive equation (ICE) is coupled to continuum mechanics in a finite-element method (FEM) scheme that tracks the flow history through the Finger tensor. The method is suitable for a calculation of residual stresses from a "first-principles" starting point following well-understood approximations. As an example, we calculate within a schematic version of MCT the stress-induced optical birefringence pattern of an amorphous solid cast into the shape of a slab with a cylindrical obstacle and demonstrate how FEM-MCT can predict the dependence of material properties on the material's processing history.
△ Less
Submitted 24 July, 2023;
originally announced July 2023.
-
An Accelerated Dual Proximal Gradient Method for Applications in Viscoplasticity
Authors:
Timm Treskatis,
Miguel A. Moyers-Gonzalez,
Chris J. Price
Abstract:
We present a very simple and fast algorithm for the numerical solution of viscoplastic flow problems without prior regularisation. Compared to the widespread alternating direction method of multipliers (ADMM / ALG2), the new method features three key advantages: firstly, it accelerates the worst-case convergence rate from $O(1/\sqrt{k})$ to $O(1/k)$, where $k$ is the iteration counter. Secondly, e…
▽ More
We present a very simple and fast algorithm for the numerical solution of viscoplastic flow problems without prior regularisation. Compared to the widespread alternating direction method of multipliers (ADMM / ALG2), the new method features three key advantages: firstly, it accelerates the worst-case convergence rate from $O(1/\sqrt{k})$ to $O(1/k)$, where $k$ is the iteration counter. Secondly, even for nonlinear constitutive models like those of Casson or Herschel-Bulkley, no nonlinear systems of equations have to be solved in the subproblems of the algorithm. Thirdly, there is no need to augment the Lagrangian, which eliminates the difficulty of choosing a penalty parameter heuristically.
In this paper, we transform the usual velocity-based formulation of viscoplastic flow problems to a dual formulation in terms of the stress. For the numerical solution of this dual problem we apply FISTA, an accelerated first-order optimisation algorithm from the class of so-called proximal gradient methods. Finally, we conduct a series of numerical experiments, focussing on stationary flow in two-dimensional square cavities.
Our results confirm that Algorithm FISTA*, the new dual-based FISTA, outperforms state-of-the-art algorithms such as ADMM / ALG2 by several orders of magnitude. We demonstrate how this speedup can be exploited to identify the free boundary between yielded and unyielded regions with previously unknown accuracy. Since the accelerated algorithm relies solely on Stokes-type subproblems and nonlinear function evaluations, existing code based on augmented Lagrangians would require only few minor adaptations to obtain an implementation of FISTA*.
△ Less
Submitted 24 September, 2016; v1 submitted 16 September, 2015;
originally announced September 2015.
-
A trust-region SQP method for the numerical approximation of viscoplastic fluid flow
Authors:
Timm Treskatis,
Miguel A. Moyers-Gonzalez,
Chris J. Price
Abstract:
We present a new approach to the problem of stationary viscoplastic duct flow as modelled by the Herschel-Bulkley model, with Bingham fluids included as a special case.
While the mathematical formulation of this problem is conventionally based on a variational inequality, or equivalently, on a nonsmooth minimisation problem for the flow velocity, we suggest an alternative approach. Considering t…
▽ More
We present a new approach to the problem of stationary viscoplastic duct flow as modelled by the Herschel-Bulkley model, with Bingham fluids included as a special case.
While the mathematical formulation of this problem is conventionally based on a variational inequality, or equivalently, on a nonsmooth minimisation problem for the flow velocity, we suggest an alternative approach. Considering the Lagrangian dual in terms of the stress, rather than the velocity, turns out to be advantageous in numerous ways. The objective functional possesses higher regularity, which ensures applicability of second order methods. Our numerical experiments with a trust-region SQP algorithm also demonstrate clearly superior performance compared to the widely used augmented Lagrangian method, although no artificial regularisation is introduced into the problem.
Hence, besides providing a new theoretical angle to a classical problem, our results also pave the way for an entirely new class of numerical approaches to simulating flows of viscoplastic fluids.
△ Less
Submitted 14 October, 2015; v1 submitted 29 April, 2015;
originally announced April 2015.