pith. sign in

arxiv: 2607.00896 · v1 · pith:XEXFWEOKnew · submitted 2026-07-01 · 🧮 math.NA · cs.NA

Certification of PGD reduced-order models with separated spatial variables

Jean Ruel , Ludovic Chamoin , Fr\'ed\'eric Legoll , Arthur Leb\'ee This is my paper

Pith reviewed 2026-07-02 07:40 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords PGDmodel order reductionerror estimationConstitutive Relation Errordiffusion problemsplate-like domainscertificationadaptive strategy
0
0 comments X

The pith

A guaranteed global error estimate for PGD approximations of diffusion problems is derived using the Constitutive Relation Error method with a procedure for equilibrated fluxes.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper develops a certification method for Proper Generalized Decomposition (PGD) reduced-order models that use separation of spatial variables, suited to plate and shell geometries. For diffusion problems in plate-like domains, it provides guaranteed error bounds based on the Constitutive Relation Error (CRE) approach. The central challenge addressed is building equilibrated fluxes that align with the separated representation, leading to an adaptive strategy that controls both discretization error and the number of PGD modes. The method also extends to error control for specific quantities of interest.

Core claim

We introduce a guaranteed global error estimate associated with the PGD approximation. To this end, the error bounds are derived from the Constitutive Relation Error (CRE) method. The main difficulty of this approach lies in the construction of equilibrated fluxes, for which a dedicated procedure is proposed. Based on the resulting estimator, an adaptive strategy is developed to control both the discretization error and the number of PGD modes. This certification procedure is further extended to the error control in quantities of interest.

What carries the argument

The Constitutive Relation Error (CRE) method, with a dedicated procedure to construct equilibrated fluxes compatible with the separated spatial variables of the PGD solution.

If this is right

  • The procedure yields an adaptive strategy that simultaneously controls discretization error and the number of PGD modes.
  • The certification extends to controlling errors in quantities of interest.
  • Numerical examples demonstrate the reliability and efficiency of the error estimator for diffusion problems in plate-like domains.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The approach may allow reliable reduced-order modeling in engineering applications involving plate structures where computational speed and accuracy certification are both required.
  • Similar flux construction techniques could be explored for other model reduction methods that use variable separation.

Load-bearing premise

A dedicated procedure can be constructed to produce equilibrated fluxes that are compatible with the separated spatial representation of the PGD solution on plate-like domains.

What would settle it

Finding a case in which the proposed equilibrated flux construction fails to produce valid fluxes for a given PGD approximation, or where the computed error bound is exceeded by the actual error in a numerical test.

Figures

Figures reproduced from arXiv: 2607.00896 by Arthur Leb\'ee, Fr\'ed\'eric Legoll, Jean Ruel, Ludovic Chamoin.

Figure 1
Figure 1. Figure 1: 2D Poisson problem with a uniform source term: exact solution [PITH_FULL_IMAGE:figures/full_fig_p010_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: 2D Poisson problem with a uniform source term: exact flux [PITH_FULL_IMAGE:figures/full_fig_p011_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: 2D Poisson problem with a prescribed flux: exact solution [PITH_FULL_IMAGE:figures/full_fig_p013_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: 2D Poisson problem with a prescribed flux: second component of the exact flux [PITH_FULL_IMAGE:figures/full_fig_p014_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: 2D Poisson problem with a prescribed flux: exact relative errors and error estimate [PITH_FULL_IMAGE:figures/full_fig_p015_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: 2D Poisson problem with a prescribed flux: initial values of [PITH_FULL_IMAGE:figures/full_fig_p020_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: 2D Poisson problem with a prescribed flux: meshes obtained at the end of the adaptive PGD procedure (left) and error convergence for [PITH_FULL_IMAGE:figures/full_fig_p021_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: 2D Poisson problem with a prescribed flux: history of convergence for the adaptive PGD strategy. [PITH_FULL_IMAGE:figures/full_fig_p021_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: 2D diffusion problem with discontinuous coefficient: PGD solution with 4 modes obtained at the end of the adaptive PGD procedure (left) and history of convergence (right) [PITH_FULL_IMAGE:figures/full_fig_p022_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: 2D diffusion problem with discontinuous coefficient: components of the SA PGD flux qˆ h n with 5 modes at the end of the adaptive PGD procedure. The final meshes (with Nω = 36 and NI = 31) are plotted in [PITH_FULL_IMAGE:figures/full_fig_p022_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: 2D di [PITH_FULL_IMAGE:figures/full_fig_p023_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: 3D diffusion problem in a laminated plate: PGD solution with 3 modes obtained at the end of the adaptive PGD procedure (left) and history of convergence (right) [PITH_FULL_IMAGE:figures/full_fig_p024_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: 3D diffusion problem in a laminated plate: components x and y of the SA PGD flux qˆ h n with 4 modes at the end of the adaptive PGD procedure. 5. Goal-oriented error estimation for the PGD approximation In this section, we show how to extend the previous tools to the framework of goal-oriented error estimation. It is indeed well known that controlling and adapting a numerical solution by measuring the err… view at source ↗
Figure 14
Figure 14. Figure 14: 3D diffusion problem in a laminated plate: initial values of η ω K (left) and η I L (right) when considering m = 3 modes. 101 102 103 104 10-2 10-1 100 Uniform refinement Adaptive refinement [PITH_FULL_IMAGE:figures/full_fig_p025_14.png] view at source ↗
Figure 15
Figure 15. Figure 15: 3D diffusion problem in a laminated plate: meshes obtained at the end of the adaptive PGD procedure (left) and error convergence for uniform and adaptive mesh refinement (right). The objective of the following is therefore to estimate the PGD error in quantity of interest Q (ePGD) = Q(u) − Q  u h m  . 5.1. Adjoint problem and upper bound on the error in quantity of interest To this end, the adjoint prob… view at source ↗
Figure 16
Figure 16. Figure 16: 2D homogenization problem: PGD solution u h m (left) and meshes (right) obtained at the end of the goal-oriented adaptive PGD procedure. Eventually, the convergence of the error estimator on Q2 is shown in [PITH_FULL_IMAGE:figures/full_fig_p028_16.png] view at source ↗
Figure 17
Figure 17. Figure 17: 2D homogenization problem: components of the PGD flux [PITH_FULL_IMAGE:figures/full_fig_p029_17.png] view at source ↗
Figure 18
Figure 18. Figure 18: 2D homogenization problem: convergence of the PGD error estimate for [PITH_FULL_IMAGE:figures/full_fig_p030_18.png] view at source ↗
read the original abstract

Model order reduction techniques have become an attractive approach for obtaining fast approximations of multidimensional problems. Besides computational efficiency, ensuring the reliability of the resulting approximations is of primary importance. This work focuses on the certification of PGD-based reduced-order models based on the separation of spatial variables, which are particularly well suited to plate and shell geometries. Considering diffusion problems defined in plate-like domains, we introduce a guaranteed global error estimate associated with the PGD approximation. To this end, the error bounds are derived from the Constitutive Relation Error (CRE) method. The main difficulty of this approach lies in the construction of equilibrated fluxes, for which a dedicated procedure is proposed. Based on the resulting estimator, an adaptive strategy is developed to control both the discretization error and the number of PGD modes. This certification procedure is further extended to the error control in quantities of interest. We provide several numerical examples illustrating the reliability and efficiency of our procedure.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Circularity Check

0 steps flagged

No significant circularity detected

full rationale

The paper derives guaranteed CRE-based error bounds for PGD approximations on plate-like domains by proposing a dedicated construction of equilibrated fluxes that respects the separated spatial representation. This construction is presented as an original technical contribution, not as a renaming, self-definition, or tautological fit. The abstract and context indicate the method extends standard CRE without reducing the central claim to its own inputs or to unverified self-citations. Numerical examples are used to illustrate reliability, preserving independent content in the derivation chain.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

Abstract-only view yields minimal ledger entries; the central claim rests on the domain assumption that plate-like geometry permits clean variable separation and on the existence of a flux-construction algorithm whose details are not supplied.

axioms (1)
  • domain assumption Diffusion problems are posed on plate-like domains where spatial variables can be separated in the PGD representation.
    Explicitly stated in the abstract as the setting for the method.

pith-pipeline@v0.9.1-grok · 5697 in / 1156 out tokens · 21851 ms · 2026-07-02T07:40:49.846272+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

27 extracted references · 23 canonical work pages

  1. [1]

    A. Nouy, Generalized spectral decomposition method for solving stochastic finite element equations: Invariant subspace problem and dedicated algorithms, Computer Methods in Applied Mechanics and Engineering 197 (51) (2008) 4718–4736.doi:10.1016/j.cma.2008.06.012

    work page doi:10.1016/j.cma.2008.06.012 2008

  2. [2]

    A. Nouy, A priori model reduction through Proper Generalized Decomposition for solving time-dependent par- tial differential equations, Computer Methods in Applied Mechanics and Engineering 199 (23) (2010) 1603– 1626.doi:10.1016/j.cma.2010.01.009

    work page doi:10.1016/j.cma.2010.01.009 2010

  3. [3]

    Chinesta, P

    F. Chinesta, P. Ladeveze, E. Cueto, A Short Review on Model Order Reduction Based on Proper Generalized Decomposition, Archives of Computational Methods in Engineering 18 (4) (2011) 395–404.doi:10.1007/ s11831-011-9064-7

    2011

  4. [4]

    Chinesta, R

    F. Chinesta, R. Keunings, A. Leygue, The Proper Generalized Decomposition for Advanced Numerical Simula- tions: A Primer, SpringerBriefs in Applied Sciences and Technology, Springer Cham, 2014

    2014

  5. [5]

    P. Ladevèze, Nonlinear Computational Structural Mechanics: New Approaches and Non-Incremental Meth- ods of Calculation, Mechanical Engineering Series, Springer, New York, NY , 1999.doi:10.1007/ 978-1-4612-1432-8

    1999

  6. [6]

    Ammar, F

    A. Ammar, F. Chinesta, P. Diez, A. Huerta, An error estimator for separated representations of highly multi- 30 dimensional models, Computer Methods in Applied Mechanics and Engineering 199 (25) (2010) 1872–1880. doi:10.1016/j.cma.2010.02.012

    work page doi:10.1016/j.cma.2010.02.012 2010

  7. [7]

    Ladevèze, L

    P. Ladevèze, L. Chamoin, On the verification of model reduction methods based on the Proper Generalized Decomposition, Computer Methods in Applied Mechanics and Engineering 200 (23) (2011) 2032–2047.doi: 10.1016/j.cma.2011.02.019

    work page doi:10.1016/j.cma.2011.02.019 2011

  8. [8]

    J. P. Moitinho de Almeida, A basis for bounding the errors of Proper Generalised Decomposition solutions in solid mechanics, International Journal for Numerical Methods in Engineering 94 (10) (2013) 961–984.doi: 10.1002/nme.4490

    work page doi:10.1002/nme.4490 2013

  9. [9]

    Allier, L

    P.-E. Allier, L. Chamoin, P. Ladevèze, Proper Generalized Decomposition computational methods on a bench- mark problem: introducing a new strategy based on Constitutive Relation Error minimization, Advanced Mod- eling and Simulation in Engineering Sciences 2 (1) (2015) 1–25.doi:10.1186/s40323-015-0038-4

    work page doi:10.1186/s40323-015-0038-4 2015

  10. [10]

    Chamoin, F

    L. Chamoin, F. Pled, P.-E. Allier, P. Ladevèze, A posteriori error estimation and adaptive strategy for PGD model reduction applied to parametrized linear parabolic problems, Computer Methods in Applied Mechanics and Engineering 327 (2017) 118–146.doi:10.1016/j.cma.2017.08.047

    work page doi:10.1016/j.cma.2017.08.047 2017

  11. [11]

    J. Reis, J. P. Moitinho de Almeida, P. Díez, S. Zlotnik, Error estimation for Proper Generalized Decomposition solutions: A dual approach, International Journal for Numerical Methods in Engineering 121 (23) (2020) 5275– 5294.doi:10.1002/nme.6452

    work page doi:10.1002/nme.6452 2020

  12. [12]

    J. Reis, J. P. Moitinho de Almeida, P. Díez, S. Zlotnik, Error estimation and adaptivity for PGD based on complementary solutions applied to a simple 1D problem, Advanced Modeling and Simulation in Engineering Sciences 7 (1) (2020) 1–22.doi:10.1186/s40323-020-00180-3

    work page doi:10.1186/s40323-020-00180-3 2020

  13. [13]

    J. Reis, J. P. Moitinho de Almeida, P. Díez, S. Zlotnik, Error estimation for Proper Generalized Decomposition solutions: Dual analysis and adaptivity for quantities of interest, International Journal for Numerical Methods in Engineering 122 (3) (2021) 752–776.doi:10.1002/nme.6559

    work page doi:10.1002/nme.6559 2021

  14. [14]

    Chamoin, P

    L. Chamoin, P. Ladevèze, Robust control of PGD-based numerical simulations, European Journal of Computa- tional Mechanics 21 (3-6) (2012) 195–207.doi:10.1080/17797179.2012.721499

    work page doi:10.1080/17797179.2012.721499 2012

  15. [15]

    Bognet, F

    B. Bognet, F. Bordeu, F. Chinesta, A. Leygue, A. Poitou, Advanced simulation of models defined in plate geometries: 3D solutions with 2D computational complexity, Computer Methods in Applied Mechanics and Engineering 201-204 (2012) 1–12.doi:10.1016/j.cma.2011.08.025

    work page doi:10.1016/j.cma.2011.08.025 2012

  16. [16]

    Bognet, A

    B. Bognet, A. Leygue, F. Chinesta, Separated representations of 3D elastic solutions in shell geometries, Ad- vanced Modeling and Simulation in Engineering Sciences 1 (1) (2014) 4.doi:10.1186/2213-7467-1-4. 31

    work page doi:10.1186/2213-7467-1-4 2014

  17. [17]

    Giner, B

    E. Giner, B. Bognet, J. J. Ródenas, A. Leygue, F. J. Fuenmayor, F. Chinesta, The Proper Generalized Decomposi- tion (PGD) as a numerical procedure to solve 3D cracked plates in linear elastic fracture mechanics, International Journal of Solids and Structures 50 (10) (2013) 1710–1720.doi:10.1016/j.ijsolstr.2013.01.039

    work page doi:10.1016/j.ijsolstr.2013.01.039 2013

  18. [18]

    Vidal, L

    P. Vidal, L. Gallimard, O. Polit, Proper Generalized Decomposition and layer-wise approach for the modeling of composite plate structures, International Journal of Solids and Structures 50 (14-15) (2013) 2239–2250.doi: 10.1016/j.ijsolstr.2013.03.034

    work page doi:10.1016/j.ijsolstr.2013.03.034 2013

  19. [19]

    Nadal, A

    E. Nadal, A. Leygue, F. Chinesta, M. Beringhier, J. J. Ródenas, F. J. Fuenmayor, A separated representation of an error indicator for the mesh refinement process under the Proper Generalized Decomposition framework, Computational Mechanics 55 (2) (2015) 251–266.doi:10.1007/s00466-014-1097-y

    work page doi:10.1007/s00466-014-1097-y 2015

  20. [20]

    Perotto, A

    S. Perotto, A. Veneziani, Coupled model and grid adaptivity in hierarchical reduction of elliptic problems, Jour- nal of Scientific Computing 60 (2014) 505–536.doi:10.1007/s10915-013-9804-y

    work page doi:10.1007/s10915-013-9804-y 2014

  21. [21]

    Chamoin, F

    L. Chamoin, F. Legoll, An introductory review on a posteriori error estimation in Finite Element computations, SIAM Review 65 (4) (2023) 963–1028.doi:10.1137/21M1464841

    work page doi:10.1137/21m1464841 2023

  22. [22]

    Becker, R

    R. Becker, R. Rannacher, An optimal control approach to a posteriori error estimation in finite element methods, Acta Numerica 10 (2001) 1–102.doi:10.1017/S0962492901000010

    work page doi:10.1017/s0962492901000010 2001

  23. [23]

    P. Ladevèze, Strict upper error bounds on computed outputs of interest in computational structural mechanics, Computational Mechanics 42 (2) (2008) 271–286.doi:10.1007/s00466-007-0201-y

    work page doi:10.1007/s00466-007-0201-y 2008

  24. [24]

    Chamoin, P

    L. Chamoin, P. Ladevèze, A non-intrusive approach of goal-oriented error estimation for evolution problems solved by the finite element method, European Journal of Computational Mechanics 17 (5-7) (2008) 981–992. doi:10.3166/remn.17.981-992

    work page doi:10.3166/remn.17.981-992 2008

  25. [25]

    Giraldi, A

    L. Giraldi, A. Nouy, G. Legrain, P. Cartraud, Tensor-based methods for numerical homogenization from high- resolution images, Computer Methods in Applied Mechanics and Engineering 254 (14-15) (2013) 154–169. doi:10.1016/j.cma.2012.10.012

    work page doi:10.1016/j.cma.2012.10.012 2013

  26. [26]

    Le Bris, F

    C. Le Bris, F. Legoll, W. Minvielle, Special quasirandom structures: A selection approach for stochastic homog- enization, Monte Carlo Methods and Applications 22 (1) (2016) 25–54.doi:doi:10.1515/mcma-2016-0101

    work page doi:10.1515/mcma-2016-0101 2016

  27. [27]

    Bourgeat, A

    A. Bourgeat, A. Piatnitski, Approximations of effective coefficients in stochastic homogenization, Annales de l’Institut Henri Poincaré (B) Probability and Statistics 40 (2) (2004) 153–165.doi:https://doi.org/10. 1016/j.anihpb.2003.07.003. 32

    2004