Source-linked AI summary
A linear mass-lumped finite element method for the Landau-Lifshitz-Gilbert equation: unconditional energy dissipation and length preservation
Qiumei Huang, Binghong Li, Xiaoli Li, Cheng Wang, Jiang Yang
TL;DR
Standard finite element discretizations struggle to preserve both the LLG equation’s unit-length constraint and original energy dissipation. The paper develops a linear mass-lumped finite element scheme using interpolation, renormalization, and a hybrid finite element–finite difference framework, obtaining nodewise length preservation, unconditional energy dissipation, and an optimal-order error estimate.
Problem
Standard finite element methods typically fail to maintain the exact unit-length constraint, while simultaneously preserving energy dissipation remains challenging.
Method
The paper proposes a fully discrete linear mass-lumped finite element scheme combining interpolation operators, mass lumping, renormalization, and finite-difference-inspired nodewise analysis.
Results
The scheme preserves nodewise unit length and unconditional energy dissipation while supporting an optimal-order error estimate.
Takeaways & Limitations
The method integrates finite element geometric flexibility with finite-difference pointwise constraint preservation for constrained dissipative systems.
Takeaways & Limitations
The error analysis relies heavily on nodewise finite-difference properties, creating an essential difficulty for finite element spatial approximation.
Abstract
from arXiv · showhide
We develop a linear, unconditionally energy-dissipative, mass-lumped finite element method for the highly nonlinear Landau--Lifshitz--Gilbert (LLG) equation on quasi-uniform triangular meshes. The method is built on a projection strategy for enforcing the nonconvex pointwise constraint $|\mathbf{m}| = 1$, whose simultaneous preservation with unconditional energy stability remains challenging for standard finite element discretizations. The key innovation is a unified hybrid finite element-finite difference framework that underlies both the design and the analysis of the proposed method. In the scheme construction, we exploit the weak formulation and nodal structure of mass-lumped finite element method, while incorporating suitable interpolation operators and a node-wise length-preserving mechanism inspired by finite difference discretizations. This combination yields a linear scheme that preserves the node-wise unit-length constraint and satisfies a discrete energy dissipation law. The same hybrid framework also plays a central role in the error analysis, where the weak formulation and quasi-uniform mesh structure of finite elements are combined with interpolation-based and nodewise finite difference method to control the strongly nonlinear damping term and to establish an optimal-order error estimate. More importantly, the proposed method provides a unified framework that systematically integrates the geometric flexibility of finite element method with the pointwise constraint-preserving property of finite difference method, and thus offers a general strategy for designing and analyzing structure-preserving discretizations of constrained dissipative systems. Numerical experiments, including a classical blow-up simulation, confirm the predicted accuracy, energy dissipation, and robustness of the method.
1. Introduction.
The paper addresses the difficulty of simultaneously preserving the LLG equation’s unit-length constraint and energy dissipation, proposing a linear finite element scheme with a hybrid finite element–finite difference design and analysis.
- The LLG equation models magnetization dynamics in ferromagnetic materials as a time-dependent nonlinear equation.
- Standard finite element methods typically fail to preserve the exact non-convex unit-length constraint |m| = 1.
- Renormalization enforces unit length by normalizing intermediate magnetization, but rigorous error analysis remains challenging.
- Existing normalization methods are not known to be unconditionally energy stable, while some alternatives require modified energy laws or new finite element spaces at each time step.
- The paper proposes a fully discrete linear finite element scheme combining interpolation, mass lumping, and renormalization to enforce |m| = 1 and unconditional original energy dissipation.
- Its hybrid framework combines finite element structure with finite-difference nodewise properties for scheme construction and error analysis.
2. Some preliminary notations and results.
This section establishes the notation, mesh assumptions, grid operators, mass-lumped finite element framework, interpolation estimates, and discrete inequalities used later in the analysis.
- Notation and mesh assumptions: The paper fixes Sobolev-space notation, time-step notation, quasi-uniform triangular partitions, and primarily periodic or homogeneous Neumann boundary conditions.The mesh consists mainly of adjacent triangle pairs forming approximate parallelograms.
- Mass-lumped finite elements: The mass-lumped finite element setting uses piecewise linear functions, vertex-based Lagrange interpolation, and a discrete inner product defined through nodal values.The section also records an estimate for the discrepancy between lumped and standard inner products.
- Grid construction: On rectangular domains, right-triangular meshes generated from quasi-uniform grids support the finite-difference operators, with general triangulations handled under corresponding mesh conditions.The extension to general domains is described as natural under the stated conditions.
- Grid operators: The computational grid, grid-function space, discrete gradient, and componentwise matrix operations are introduced for scalar and vector grid functions.The discrete gradient inner product combines corresponding matrix entries through multiplication and summation.
- Analytical tools: Interpolation and inverse inequalities, matched central interpolation identities, and a discrete Gronwall lemma provide analytical tools for subsequent estimates.The stated inequalities apply under dimension- and exponent-dependent conditions, while the constants in the inverse estimate depend only on the domain.
3. An energy-dissipative and length preserving finite element scheme.
The section introduces a projection-based finite element scheme that preserves node-wise unit length and original energy dissipation, then establishes stability, uniqueness, and mesh extensions.
- Stability and preservation: The constructed scheme preserves the node-wise unit-length constraint and the original energy dissipation law.The analysis introduces a weak Laplacian and derives the energy law from the scheme identities and a node-wise estimate.
- Scheme construction: The scheme combines time projection with mass-lumped finite element spatial discretization to construct a fully discrete method.The projection includes a node-wise normalization operation that is simple to implement.
- Scheme construction: The method is connected to the original LLG equation through mass-lumped integration by parts and vector identities at mesh nodes.The resulting reformulation uses the approximation mt = −βm × ∆m −γm × (m × ∆m).
- Stability and preservation: The scheme admits unconditional stability, including a uniform H1 bound for the numerical solution.The stability result also yields existence and uniqueness of the scheme solution.
- Mesh considerations: The stability analysis extends to general triangular mesh partitions under a Delaunay-type condition.Under this condition, the stated inequality keeps the analysis valid on general unstructured triangular meshes.
4. Theoretical analysis.
The analysis reformulates the nonlinear scheme to avoid a difficult damping term, then combines mass-lumped finite element identities, interpolation estimates, and nodewise finite-difference expansions to establish optimal-order convergence under stated regularity and mesh conditions.
- 4. Theoretical analysis.: The equivalent reformulation avoids the complicated nonlinear damping term and enables optimal-rate error analysis through linearized stability estimates.Summation by parts is used to derive the reformulation before the convergence argument.
- 4. Theoretical analysis.: Mass-lumping and interpolation operators express finite element inner products and nonlinear gradients through finite-difference approximations on triangular cells.The construction uses node values and constant cellwise gradients to obtain the required expansions.
- 4. Theoretical analysis.: Nodewise identities and interpolation values are used to transform nonlinear inner-product terms into forms needed for the convergence proof.These identities support the discrete expansions and consistency estimates used later in the analysis.
- 4. Theoretical analysis.: The consistency analysis combines interpolation and nodewise finite-difference estimates, with super-convergence and auxiliary interpolation results supplying precise truncation-error bounds.The argument assumes a smooth exact solution and uses estimates for the Lagrange interpolation operator.
- 4. Theoretical analysis.: Under the theorem’s regularity, smallness, and scaling assumptions, the fully discrete scheme satisfies an optimal error estimate with L2 error bounded by C0(h2+∆t).The proof closes through discrete Gronwall and comparison estimates between discrete and continuous norms.
5. Numerical results.
Numerical experiments verify the method’s predicted convergence, energy dissipation, quasi-uniform mesh performance, and behavior in a classical blowup simulation.
- 5.1. Convergence test.: Second-order spatial and first-order temporal convergence are observed in the discrete L2_h norm, matching Theorem 4.5.The tests use Δt = h^2 for spatial convergence and Δt = h for temporal convergence.
- 5.1. Convergence test.: The convergence tables on quasi-uniform grids satisfy the mesh condition hmax/hmin ≤ 3.This supports testing the method beyond uniform meshes.
- 5.2. Energy dissipation.: The original energy decreases monotonically across different spatial mesh sizes, confirming dissipative behavior.The simulation uses time step Δt = 0.01.
- 5.3. Phenomenon of blowup.: In the blowup simulation, the magnetization preserves (0, 0, 1)^T at the origin while turning toward (0, 0, −1)^T nearby over time.Snapshots are taken at t = 0, 0.06, 0.15, 0.30, 0.32, and 0.35.
- 5.2. Energy dissipation.: Centrally refined quasi-uniform 128 × 128 meshes produce monotonic energy decay and match the uniform 256 × 256 reference considerably well.The comparison includes a standard uniform 128 × 128 mesh and the reference solution from [26].
6. Conclusion.
The paper proposes and analyzes a linear mass-lumped finite element method that simultaneously preserves node-wise unit length and unconditional original energy dissipation. A hybrid finite element–finite difference analysis yields optimal-order convergence on triangular meshes, while higher-order extensions remain an open challenge.
- 6. Conclusion.: The scheme enforces |m_h| = 1 at every finite element node and preserves unconditional original energy dissipation.The method is presented as a classical finite element scheme distinct from the tangent-plane finite element method.
- 6. Conclusion.: Interpolation operators combine finite element and finite difference analysis to control the nonlinear term and derive optimal-order error estimates.The analysis uses an equivalent weak formulation and triangular meshes.
- 6. Conclusion.: Developing higher-order finite element methods that retain both unit-length preservation and unconditional original energy dissipation remains a significant challenge.