Locking-free hybrid high-order method for linear elasticity
2404.02768
0
0
Abstract
The hybrid-high order (HHO) scheme has many successful applications including linear elasticity as the first step towards computational solid mechanics. The striking advantage is the simplicity among other higher-order nonconforming schemes and its geometric flexibility as a polytopal method on the expanse of a parameter-free refined stabilization. The classical suggestion of a locking-free HHO discretization requires a split of the the reconstruction terms with an additional reconstruction of the divergence operator that might be motivated by the Stokes equations for the robust approximation in the incompressible limit, when one Lam'e parameter $lambdatoinfty$ becomes very large. This paper utilizes just one reconstruction operator for the linear Green strain and therefore does not rely on a split in deviatoric and spherical behavior. The a priori error analysis provides quasi-best approximation with $lambda$-independent equivalence constants. The reliable and (up to data oscillations) efficient a posteriori error estimates are stabilization-free and $lambda$-robust. The error analysis is carried out on simplicial meshes to allow conforming piecewise polynomials finite elements in the kernel of the stabilization terms. Numerical benchmarks provide empirical evidence for optimal convergence rates of the a posteriori error estimator in some associated adaptive mesh-refining algorithm also in the incompressible limit.
Create account to get full access
Overview
- This paper presents a new numerical method called the "locking-free hybrid high-order method" for solving linear elasticity problems.
- The method aims to overcome the "locking" issue that can arise in traditional numerical approaches, which can lead to inaccurate results when modeling materials with near-incompressible behavior.
- The proposed method combines elements from different numerical techniques, including the hybrid high-order (HHO) method and a stabilization strategy, to achieve a locking-free and robust solution.
Plain English Explanation
The paper describes a new way to solve a particular type of engineering problem called "linear elasticity." This problem deals with how materials, like metals or plastics, deform or change shape when force is applied to them. The traditional methods used to solve these problems can sometimes run into issues, particularly when the material being modeled is very rigid or "incompressible," meaning it's difficult to compress or squish.
The new method proposed in this paper tries to address this problem by combining several different numerical techniques in a smart way. The key idea is to use a approach called the "hybrid high-order (HHO) method," which allows for more flexibility in how the problem is discretized (broken down into smaller pieces for computation). This HHO method is then combined with an additional "stabilization" step that helps prevent the inaccuracies that can arise when modeling very rigid materials.
The end result is a new numerical method that can solve linear elasticity problems, even for materials that are nearly incompressible, without running into the "locking" issues that plague traditional approaches. This could be useful for engineers and scientists working on a variety of applications, from designing buildings and vehicles to modeling the behavior of biological tissues.
Technical Explanation
The paper introduces a new numerical method called the "locking-free hybrid high-order (LFHHO) method" for solving linear elasticity problems. The key motivation is to address the issue of "locking," which can occur when using traditional numerical techniques, such as the finite element method, to model materials with near-incompressible behavior.
The LFHHO method combines elements from the hybrid high-order (HHO) method and a stabilization strategy to achieve a locking-free and robust solution. The HHO method provides a flexible framework for discretizing the problem domain, allowing for the use of general polytopal meshes and high-order approximations. The stabilization strategy, inspired by the second-order correction method for loosely coupled discretizations, helps to mitigate the locking phenomenon by introducing an additional stabilization term in the formulation.
The authors present a thorough mathematical analysis of the LFHHO method, proving its well-posedness, stability, and optimal convergence rates. They also conduct numerical experiments to demonstrate the method's performance on a variety of test cases, including the classic Kirsch problem and a more complex problem involving a perforated plate.
Critical Analysis
The paper presents a well-designed and rigorously analyzed numerical method for solving linear elasticity problems with near-incompressible materials. The authors have clearly identified the locking issue as a significant challenge in traditional approaches and have developed a novel solution that combines established techniques in a thoughtful manner.
One potential limitation of the LFHHO method is that it may require more computational resources compared to lower-order methods, particularly for large-scale problems. The authors acknowledge this tradeoff and suggest that future work could focus on developing more efficient implementation strategies or adaptive refinement techniques to mitigate the increased computational cost.
Additionally, while the numerical experiments demonstrate the method's effectiveness on several benchmark problems, it would be valuable to see its performance tested on more realistic, industry-relevant applications. This could help to further validate the method's practical utility and identify any additional challenges that may arise in real-world scenarios.
Overall, the locking-free hybrid high-order method presented in this paper represents a significant contribution to the field of computational mechanics and could have important implications for a wide range of engineering and scientific applications.
Conclusion
This paper introduces a new numerical method called the "locking-free hybrid high-order (LFHHO) method" for solving linear elasticity problems, particularly those involving near-incompressible materials. The method combines elements from the hybrid high-order (HHO) approach and a stabilization strategy to overcome the "locking" issues that can arise in traditional numerical techniques.
The authors provide a thorough mathematical analysis of the LFHHO method, proving its well-posedness, stability, and optimal convergence rates. Numerical experiments on benchmark problems demonstrate the method's effectiveness in producing accurate and locking-free solutions.
The LFHHO method represents a significant advancement in the field of computational mechanics and could have important implications for a wide range of engineering and scientific applications, from designing robust structures to modeling the behavior of biological tissues. Future research could focus on developing more efficient implementation strategies and exploring the method's performance on real-world, industry-relevant problems.
This summary was produced with help from an AI and may contain inaccuracies - check out the links to read the original source documents!
Related Papers
First-Order Methods for Linearly Constrained Bilevel Optimization
Guy Kornowski, Swati Padmanabhan, Kai Wang, Zhe Zhang, Suvrit Sra
0
0
Algorithms for bilevel optimization often encounter Hessian computations, which are prohibitive in high dimensions. While recent works offer first-order methods for unconstrained bilevel problems, the constrained setting remains relatively underexplored. We present first-order linearly constrained optimization methods with finite-time hypergradient stationarity guarantees. For linear equality constraints, we attain $epsilon$-stationarity in $widetilde{O}(epsilon^{-2})$ gradient oracle calls, which is nearly-optimal. For linear inequality constraints, we attain $(delta,epsilon)$-Goldstein stationarity in $widetilde{O}(d{delta^{-1} epsilon^{-3}})$ gradient oracle calls, where $d$ is the upper-level dimension. Finally, we obtain for the linear inequality setting dimension-free rates of $widetilde{O}({delta^{-1} epsilon^{-4}})$ oracle complexity under the additional assumption of oracle access to the optimal dual variable. Along the way, we develop new nonsmooth nonconvex optimization methods with inexact oracles. We verify these guarantees with preliminary numerical experiments.
6/19/2024
🛠️
A simple and efficient convex optimization based bound-preserving high order accurate limiter for Cahn-Hilliard-Navier-Stokes system
Chen Liu, Beatrice Riviere, Jie Shen, Xiangxiong Zhang
0
0
For time-dependent PDEs, the numerical schemes can be rendered bound-preserving without losing conservation and accuracy, by a post processing procedure of solving a constrained minimization in each time step. Such a constrained optimization can be formulated as a nonsmooth convex minimization, which can be efficiently solved by first order optimization methods, if using the optimal algorithm parameters. By analyzing the asymptotic linear convergence rate of the generalized Douglas-Rachford splitting method, optimal algorithm parameters can be approximately expressed as a simple function of the number of out-of-bounds cells. We demonstrate the efficiency of this simple choice of algorithm parameters by applying such a limiter to cell averages of a discontinuous Galerkin scheme solving phase field equations for 3D demanding problems. Numerical tests on a sophisticated 3D Cahn-Hilliard-Navier-Stokes system indicate that the limiter is high order accurate, very efficient, and well-suited for large-scale simulations. For each time step, it takes at most $20$ iterations for the Douglas-Rachford splitting to enforce bounds and conservation up to the round-off error, for which the computational cost is at most $80N$ with $N$ being the total number of cells.
4/1/2024
✨
A global approach for the redefinition of higher-order flexibility and rigidity
Georg Nawratil
0
0
The famous example of the double-Watt mechanism given by Connelly and Servatius raises some problems concerning the classical definitions of higher-order flexibility and rigidity, respectively, as they attest the cusp configuration of the mechanism a third-order rigidity, which conflicts with its continuous flexion. Some attempts were done to resolve the dilemma but they could not settle the problem. As cusp mechanisms demonstrate the basic shortcoming of any local mobility analysis using higher-order constraints, we present a global approach inspired by Sabitov's finite algorithm for testing the bendability of a polyhedron, which allows us (a) to compute iteratively configurations with a higher-order flexion and (b) to come up with a proper redefinition of higher-order flexibility and rigidity. We also give algorithms for computing the flexion orders as well as the associated flexes. The presented approach is demonstrated on several examples (double-Watt mechanisms and Tarnai's Leonardo structure). Moreover, we determine all configurations of a given 3-RPR manipulator with a third-order flexion and present a corresponding joint-bar framework of flexion order 23.
5/27/2024
Physics-guided weak-form discovery of reduced-order models for trapped ultracold hydrodynamics
Reuben R. W. Wang, Daniel Messenger
0
0
We study the relaxation of a highly collisional, ultracold but nondegenerate gas of polar molecules. Confined within a harmonic trap, the gas is subject to fluid-gaseous coupled dynamics that lead to a breakdown of first-order hydrodynamics. An attempt to treat these higher-order hydrodynamic effects was previously made with a Gaussian ansatz and coarse-graining model parameter [R. R. W. Wang & J. L. Bohn, Phys. Rev. A 108, 013322 (2023)], leading to an approximate set of equations for a few collective observables accessible to experiments. Here we present substantially improved reduced-order models for these same observables, admissible beyond previous parameter regimes, discovered directly from particle simulations using the WSINDy algorithm (Weak-form Sparse Identification of Nonlinear Dynamics). The interpretable nature of the learning algorithm enables estimation of previously unknown physical quantities and discovery of model terms with candidate physical mechanisms, revealing new physics in mixed collisional regimes. Our approach constitutes a general framework for data-driven model identification leveraging known physics.
6/12/2024