pith. sign in

arxiv: 2604.19077 · v1 · submitted 2026-04-21 · 🧮 math.NA · cs.NA

High-Order Multi-Scale Method and Its Convergence Analysis for Nonlinear Thermo-Electro-Mechanical Coupling Problems of Composite Structures

Hao Dong This is my paper

Pith reviewed 2026-05-10 02:26 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords multi-scale asymptotic methodthermo-electro-mechanical couplingnonlinear problemscomposite structuresconvergence analysisTaylor serieserror estimationperiodic heterogeneity
0
0 comments X

The pith

High-order multi-scale asymptotic model with Taylor corrections simulates nonlinear thermo-electro-mechanical coupling in periodic composites with explicit convergence rates.

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

The paper develops a high-order multi-scale method for time-dependent nonlinear problems coupling thermal, electrical, and mechanical effects in composite materials that have repeating patterns along with temperature-dependent properties and electrical heating. It combines asymptotic expansions with Taylor series to generate correction terms that capture fine-scale behavior without resolving every detail directly. Local error analysis shows that balances of heat, charge, and stress hold pointwise, while global estimates give explicit rates at which the approximation improves. An efficient algorithm splits work into offline preparation and online evaluation to cut computational expense. Tests on problems with strongly varying coefficients confirm gains in accuracy and speed over conventional approaches.

Core claim

By employing the multi-scale asymptotic approach and the Taylor series technique, a high-accuracy multi-scale asymptotic model featuring novel high-order correction terms is established for nonlinear multi-physics simulation of periodic solid structures. A local point-wise error analysis is derived to theoretically and physically illustrate the local balance preserving of heat quantity, electric charge and stress, thereby enabling high-accuracy multi-scale computation. Moreover, a global error estimation is obtained that provides an explicit convergence rate for high-order multi-scale solutions. An efficient numerical algorithm featuring off-line and on-line stages is presented with a full a

What carries the argument

Multi-scale asymptotic expansion combined with Taylor series to generate high-order correction terms that handle nonlinear thermo-electro-mechanical couplings and spatial heterogeneity in periodic structures.

If this is right

  • Local point-wise error analysis confirms preservation of physical balances for heat, electric charge, and stress.
  • Global error estimates supply explicit convergence rates as the small-scale period shrinks.
  • The off-line and on-line algorithm lowers cost for problems whose material coefficients oscillate strongly.
  • Numerical experiments verify higher accuracy and lower run times than standard methods for the targeted nonlinear time-dependent cases.

Where Pith is reading between the lines

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

  • The same expansion strategy could be carried over to other nonlinear multi-physics models that involve periodic microstructures.
  • Design optimization loops for composites under combined loads would become cheaper because the online stage re-uses precomputed data.
  • Performance when periodicity holds only approximately could be checked by adding controlled perturbations to the test geometries.

Load-bearing premise

The composite structures are strictly periodic so that the multi-scale asymptotic expansion and Taylor series remain valid for the time-dependent nonlinear couplings that include temperature-dependent properties and Joule heating.

What would settle it

Compare the high-order multi-scale solution against a highly resolved direct numerical simulation on a periodic test domain and check whether the observed error decreases at the exact rate given by the global error estimation when the period length is reduced toward zero.

read the original abstract

This study proposes a high-order multi-scale method tailored for time-dependent nonlinear thermo-electro-mechanical coupling problems of composite structures with highly spatial heterogeneity, which incorporate temperature-dependent material properties and Joule heating effect. By employing the multi-scale asymptotic approach and the Taylor series technique, a high-accuracy multi-scale asymptotic model featuring novel high-order correction terms is established for nonlinear multi-physics simulation of periodic solid structures. A local point-wise error analysis is derived to theoretically and physically illustrate the local balance preserving of heat quantity, electric charge and stress,thereby enabling high-accuracy multi-scale computation. Moreover, a global error estimation is obtained that provides an explicit convergence rate for high-order multi-scale solutions. Furthermore, an efficient numerical algorithm featuring with off-line and on-line stages is presented meticulously, accompanied by a corresponding error analysis. Numerical experiments are conducted to showcase the competitive advantages of the proposed method for simulating the time-dependent nonlinear thermo-electro-mechanical coupling problems with highly oscillatory and discontinuous coefficients, demonstrating superior numerical accuracy and reduced computational cost.

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.

Referee Report

2 major / 2 minor

Summary. The manuscript proposes a high-order multi-scale method for time-dependent nonlinear thermo-electro-mechanical coupling problems in periodic composite structures featuring temperature-dependent properties and Joule heating. It combines multi-scale asymptotic analysis with Taylor series expansions to derive a model with novel high-order correction terms, establishes local point-wise error analysis to demonstrate preservation of heat, charge, and stress balances, obtains global error estimates providing explicit convergence rates, presents an efficient offline-online numerical algorithm with accompanying error analysis, and reports numerical experiments on problems with highly oscillatory and discontinuous coefficients.

Significance. If the derivations hold, particularly the applicability of the Taylor-augmented expansion and the supporting error bounds, the work would advance multi-scale modeling for nonlinear multi-physics problems by supplying higher-order accuracy, explicit convergence guarantees, and computational efficiency via the two-stage algorithm.

major comments (2)
  1. [§3] §3 (multi-scale asymptotic expansion with Taylor series): The Taylor expansion around the slow-scale variable to linearize nonlinearities (temperature-dependent moduli and quadratic Joule heating) implicitly requires C^2 regularity in the macroscopic variable. With highly oscillatory and discontinuous coefficients across cell interfaces, this differentiability is not guaranteed, which directly affects the validity of the high-order correctors and the global error estimate.
  2. [§4] §4 (global error estimation, Theorem 4.1): The explicit convergence rate for the high-order multi-scale solution is stated, but the analysis does not address how discontinuities in the coefficients propagate into the Taylor remainder terms or affect the local balance preservation; this is load-bearing for the central claim of superior accuracy and rate.
minor comments (2)
  1. [Abstract] Abstract: The local point-wise error analysis is described as illustrating physical balance preservation, but the specific norms or function spaces are not indicated, reducing clarity.
  2. [§5] §5 (numerical experiments): Error plots versus scale parameter or mesh size should overlay the predicted theoretical rates to allow direct visual verification of the claimed convergence.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the thorough review and valuable comments on our manuscript. We address each major comment point by point below, providing clarifications and indicating where revisions will be made to improve the presentation and rigor of the analysis.

read point-by-point responses

  1. Referee: [§3] §3 (multi-scale asymptotic expansion with Taylor series): The Taylor expansion around the slow-scale variable to linearize nonlinearities (temperature-dependent moduli and quadratic Joule heating) implicitly requires C^2 regularity in the macroscopic variable. With highly oscillatory and discontinuous coefficients across cell interfaces, this differentiability is not guaranteed, which directly affects the validity of the high-order correctors and the global error estimate.

    Authors: The Taylor series expansion is applied to the nonlinear constitutive relations with respect to the slow-scale macroscopic variables (e.g., the macroscopic temperature), while the fast-scale periodic oscillations and discontinuities are accounted for through the solution of the cell problems on the reference cell. The required C^2 regularity is assumed for the macroscopic solution, which follows from the well-posedness of the homogenized nonlinear system under standard assumptions on the averaged coefficients. The discontinuities at microscale interfaces do not violate the differentiability in the slow variable, as the expansion is performed after averaging. We will add a detailed discussion of the regularity assumptions in Section 3 of the revised manuscript to clarify this point. The high-order correctors remain valid under these assumptions, and the local error analysis holds in the distributional sense or away from interfaces. revision: partial

  2. Referee: [§4] §4 (global error estimation, Theorem 4.1): The explicit convergence rate for the high-order multi-scale solution is stated, but the analysis does not address how discontinuities in the coefficients propagate into the Taylor remainder terms or affect the local balance preservation; this is load-bearing for the central claim of superior accuracy and rate.

    Authors: In the proof of Theorem 4.1, the Taylor remainder terms are estimated using the boundedness of the second derivatives of the nonlinear functions and the periodicity of the coefficients. The propagation of discontinuities is controlled by considering the weak formulation and integrating over subdomains where the coefficients are smooth, with jump conditions handled at the interfaces. The local balance preservation is demonstrated pointwise in the sense of distributions, preserving the integral balances of heat, charge, and stress. To strengthen the manuscript, we will expand the error analysis section to include more explicit bounds on the remainder terms accounting for the discontinuous coefficients and provide additional details on the local preservation properties. revision: partial

Circularity Check

0 steps flagged

No circularity: derivation proceeds from asymptotic expansion and Taylor series without reduction to inputs

full rationale

The paper establishes its high-order multi-scale model by applying the multi-scale asymptotic approach combined with Taylor series expansion to handle nonlinearities (temperature-dependent properties and Joule heating) in periodic structures, then derives local point-wise error analysis for balance preservation and global error estimation for convergence rates. These steps are presented as direct consequences of the expansion technique and standard error bounding arguments rather than any fitted parameter renamed as prediction, self-definitional loop, or load-bearing self-citation chain. The numerical algorithm and experiments serve as validation but do not form the core derivation. No quoted equation or claim reduces the output to the input by construction, confirming the analysis is self-contained.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Abstract provides insufficient detail to identify concrete free parameters, axioms, or invented entities; the approach appears to rest on standard assumptions of periodicity and asymptotic validity for nonlinear problems.

pith-pipeline@v0.9.0 · 5477 in / 1123 out tokens · 32213 ms · 2026-05-10T02:26:21.242724+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

40 extracted references · 40 canonical work pages

  1. [1]

    Abdulle, E

    A. Abdulle, E. Weinan, B. Engquist, and E. Vanden-Eijnden , The heterogeneous multi- scale method, Acta Numerica, 21 (2012), pp. 1–87

    work page 2012

  2. [2]

    I. V. Andrianov, H. Topol, and V. V. Danishevskyy , Asymptotic analysis of heat transfer in composite materials with nonlinear thermal properties , International Journal of Heat and Mass Transfer, 111 (2017), pp. 736–754

    work page 2017

  3. [3]

    Bensoussan, J.-L

    A. Bensoussan, J.-L. Lions, and G. Papanicolaou, Asymptotic Analysis for Periodic Struc- tures, vol. 5 of Studies in Mathematics and its Applications, North-Holland, Amsterdam, 38 H. DONG 1978

    work page 1978

  4. [4]

    H. G. C and S.-V. T, Boundary integral formulations for transient linear thermoelasticity with combined-type boundary conditions, SIAM Journal on Mathematical Analysis, 53 (2021), pp. 3888–3911

    work page 2021

  5. [5]

    L.-Q. Cao, Multiscale asymptotic expansion and finite element methods for the mixed bound- ary value problems of second order elliptic equation in perforated domains , Numerische Mathematik, 103 (2006), pp. 11–45

    work page 2006

  6. [6]

    Cao, J.-Z

    L.-Q. Cao, J.-Z. Cui, and D.-C. Zhu, Multiscale asymptotic analysis and numerical simulation for the second order helmholtz equations with rapidly oscillating coefficients over general convex domains, SIAM Journal on Numerical Analysis, 40 (2002), pp. 543–577

    work page 2002

  7. [7]

    P. W. Chung, K. K. Tamma, and R. R. Namburu, Homogenization of temperature-dependent thermal conductivity in composite materials, Journal of Thermophysics and Heat Transfer, 15 (2001), pp. 10–17

    work page 2001

  8. [8]

    Cioranescu and P

    D. Cioranescu and P. Donato, An introduction to homogenization, Oxford university press, 1999

    work page 1999

  9. [9]

    H. Dong, J. Cui, Y. Nie, R. Ma, K. Jin, and D. Huang, Multi-scale computational method for nonlinear dynamic thermo-mechanical problems of composite materials with temperature- dependent properties, Communications in Nonlinear Science and Numerical Simulation, 118 (2023), p. 107000

    work page 2023

  10. [10]

    H. Dong, X. Guan, and Y. Nie , Multiscale method and convergence analysis for coupled nonlinear thermomechanical problems in heterogeneous shells , SIAM Journal on Scientific Computing, 47 (2025), pp. B190–B219

    work page 2025

  11. [11]

    H. Dong, Z. Yang, and Y. Nie , Higher-order multiscale method and its convergence analysis for nonlinear thermo-electric coupling problems of composite structures , Multiscale Mod- eling & Simulation, 23 (2025), pp. 1744–1772

    work page 2025

  12. [12]

    H. Dong, X. Zheng, J. Cui, Y. Nie, Z. Yang, and Q. Ma , Multi-scale computational method for dynamic thermo-mechanical performance of heterogeneous shell structures with orthog- onal periodic configurations, Computer Methods in Applied Mechanics and Engineering, 354 (2019), pp. 143–180

    work page 2019

  13. [13]

    H. Dong, X. Zheng, J. Cui, Y. Nie, Z. Yang, and Z. Yang , High-order three-scale com- putational method for dynamic thermo-mechanical problems of composite structures with multiple spatial scales, International Journal of Solids and Structures, 169 (2019), pp. 95– 121

    work page 2019

  14. [14]

    Dong and L.-Q

    Q.-L. Dong and L.-Q. Cao , Multiscale asymptotic expansions and numerical algorithms for the wave equations of second order with rapidly oscillating coefficients , Applied Numerical Mathematics, 59 (2009), pp. 3008–3032

    work page 2009

  15. [15]

    Efendiev, J

    Y. Efendiev, J. Galvis, G. Li, and M. Presho , Generalized multiscale finite element meth- ods. nonlinear elliptic equations , Communications in Computational Physics, 15 (2014), pp. 733–755

    work page 2014

  16. [16]

    Efendiev and T

    Y. Efendiev and T. Y. Hou , Multiscale finite element methods: theory and applications , vol. 4, Springer Science & Business Media, 2009

    work page 2009

  17. [17]

    Fern´andez and K

    J. Fern´andez and K. Kuttler , A dynamic thermoviscoelastic problem: numerical analysis and computational experiments, Quarterly journal of mechanics and applied mathematics, 63 (2010), pp. 295–314

    work page 2010

  18. [18]

    Henning and A

    P. Henning and A. M ˚alqvist, Localized orthogonal decomposition techniques for boundary value problems, SIAM Journal on Scientific Computing, 36 (2014), pp. A1609–A1634

    work page 2014

  19. [19]

    Hou, X.-H

    T. Hou, X.-H. Wu, and Z. Cai, Convergence of a multiscale finite element method for elliptic problems with rapidly oscillating coefficients , Mathematics of Computation, 68 (1999), pp. 913–943

    work page 1999

  20. [20]

    G. C. Hsiao and T. S´anchez-Vizuet, Time-domain boundary integral methods in linear ther- moelasticity, SIAM Journal on Mathematical Analysis, 52 (2020), pp. 2463–2490

    work page 2020

  21. [21]

    T. J. Hughes , Multiscale phenomena: Green’s functions, the dirichlet-to-neumann formula- tion, subgrid scale models, bubbles and the origins of stabilized methods , Computer Meth- ods in Applied Mechanics and Engineering, 127 (1995), pp. 387–401

    work page 1995

  22. [22]

    Kuttler, M

    K. Kuttler, M. Shillor, and J. Fern´andez, Existence for the thermoviscoelastic thermistor problem, Differential Equations and Dynamical Systems, 17 (2009), pp. 217–233

    work page 2009

  23. [23]

    B. Li, H. Gao, and W. Sun , Unconditionally optimal error estimates of a crank–nicolson galerkin method for the nonlinear thermistor equations, SIAM Journal on Numerical Analy- sis, 52 (2014), pp. 933–954

    work page 2014

  24. [24]

    P. Li, Y. Dong, M. Tang, J. Mao, L. J. Jiang, and H. Ba ˘gcı, Transient thermal analysis of 3-d integrated circuits packages by the dgtd method , IEEE Transactions on Components, Packaging and Manufacturing Technology, 7 (2017), pp. 862–871. MULTI-SCALE METHOD FOR THERMO-ELECTRO-MECHANICAL COUPLING PROBLEMS39

    work page 2017

  25. [25]

    Lin and Q

    Q. Lin and Q. Zhu , The Preprocessing and Preprocessing for the Finite Element Method , Shanghai Scientific & Technical Publishers, 1994

    work page 1994

  26. [26]

    Lunati and P

    I. Lunati and P. Jenny , Multiscale finite-volume method for compressible multiphase flow in porous media, Journal of Computational Physics, 216 (2006), pp. 616–636

    work page 2006

  27. [27]

    C. Ma, C. Alber, and R. Scheichl , Wavenumber explicit convergence of a multiscale gen- eralized finite element method for heterogeneous helmholtz problems , SIAM Journal on Numerical Analysis, 61 (2023), pp. 1546–1584

    work page 2023

  28. [28]

    M˚alqvist and T

    A. M˚alqvist and T. Stillfjord, Finite element convergence analysis for the thermoviscoelas- tic joule heating problem , BIT Numerical Mathematics, 57 (2017), pp. 787–810

    work page 2017

  29. [29]

    A. H. Muliana and J. S. Kim , A two-scale homogenization framework for nonlinear effective thermal conductivity of laminated composites , Acta Mechanica, 212 (2010), pp. 319–347

    work page 2010

  30. [30]

    O. A. Ole¨ınik, A. Shamaev, and G. Yosifian , Mathematical problems in elasticity and ho- mogenization, vol. 26, Elsevier, 1992

    work page 1992

  31. [31]

    A. A. Pankov, G-convergence and homogenization of nonlinear partial differential operators , vol. 422, Springer Science & Business Media, 2013

    work page 2013

  32. [32]

    Qin and Q.-S

    Q.-H. Qin and Q.-S. Yang , Macro-micro theory on multifield coupling behavior of heteroge- neous materials, Springer Berlin, Germany:, 2008

    work page 2008

  33. [33]

    X. Wang, L. Cao, and Y. Wong , Multiscale computation and convergence for coupled ther- moelastic system in composite materials , Multiscale Modeling & Simulation, 13 (2015), pp. 661–690

    work page 2015

  34. [34]

    Wu and X

    X. Wu and X. Xu, Existence for the thermoelastic thermistor problem, Journal of mathematical analysis and applications, 319 (2006), pp. 124–138

    work page 2006

  35. [35]

    Xie and M

    J. Xie and M. Swaminathan, Electrical–thermal cosimulation with nonconformal domain de- composition method for multiscale 3-d integrated systems , IEEE Transactions on Compo- nents, Packaging and Manufacturing Technology, 4 (2013), pp. 588–601

    work page 2013

  36. [36]

    Y. Xing, Y. Yang, and X. Wang , A multiscale eigenelement method and its application to periodical composite structures, Composite Structures, 92 (2010), pp. 2265–2275

    work page 2010

  37. [37]

    Yang and Y.-L

    Y.-B. Yang and Y.-L. Jiang , Optimal error estimates of a lowest-order galerkin-mixed fem for the thermoviscoelastic joule heating equations , Applied Numerical Mathematics, 183 (2023), pp. 86–107

    work page 2023

  38. [38]

    Z. Yang, J. Cui, Y. Sun, J. Liang, and Z. Yang , Multiscale analysis method for thermo- mechanical performance of periodic porous materials with interior surface radiation , In- ternational Journal for Numerical Methods in Engineering, 105 (2016), pp. 323–350

    work page 2016

  39. [39]

    W. Yu, S. Cheng, Z. Li, L. Liu, Z. Zhang, Y. Zhao, Y. Guo, and S. Liu , The application of multi-scale simulation in advanced electronic packaging , Fundamental Research, (2024)

    work page 2024

  40. [40]

    H. Zhai, Q. Wu, N. Yoshikawa, K. Xiong, and C. Chen , Space-time asymptotic expan- sion method for transient thermal conduction in the periodic composite with temperature- dependent thermal properties, Computational Materials Science, 194 (2021), p. 110470

    work page 2021