pith. machine review for the scientific record. sign in

arxiv: 2605.10957 · v1 · submitted 2026-05-04 · ❄️ cond-mat.soft · cond-mat.mtrl-sci

Recognition: 2 theorem links

· Lean Theorem

Quantifying the effects of particle clustering in random thermoelastic composites -- numerical and mean-field analyses

Authors on Pith no claims yet

Pith reviewed 2026-05-13 06:49 UTC · model grok-4.3

classification ❄️ cond-mat.soft cond-mat.mtrl-sci
keywords particle clusteringthermoelastic compositesmean-field cluster modelfinite element analysisrandom distributioneffective propertiesnearest-neighbour distance
0
0 comments X

The pith

Clustering of equal-sized spherical particles alters thermoelastic properties and local fields depending on volume fraction and mean nearest-neighbour distance

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

The paper investigates the impact of particle spatial distribution on the thermoelastic effective properties and local stress-strain fields in composites with randomly placed inclusions. Quantitative comparisons are made between full-field finite element analyses and a multi-family mean-field cluster model that computes separate mean stresses and strains for each inclusion. The study considers variations in volume fraction and mean nearest-neighbour distances for equal-sized spherical particles. Understanding these clustering effects matters for accurately predicting material performance in engineering applications without exhaustive simulations.

Core claim

The effect of space distribution of randomly-placed particles in a representative composite volume on the thermoelastic effective properties and local stress and strain distribution is analyzed through quantitative assessment using both full-field finite element analyses and the mean-field interaction model known as the cluster model, developed in the multi-family setting to study mean stress and strain separately for each inclusion.

What carries the argument

The multi-family cluster model, a mean-field interaction model that enables separate analysis of mean stress and strain for each inclusion in the representative unit cell.

If this is right

  • Effective thermoelastic properties vary with mean nearest-neighbour distance at fixed volume fraction.
  • Local stress and strain distributions inside and around inclusions change with the degree of clustering.
  • The multi-family mean-field model provides results consistent with full-field finite element analyses for the studied cases.
  • Clustering effects can be studied without resolving the full microstructure by using the cluster model.

Where Pith is reading between the lines

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

  • The approach could be tested on composites with irregular particle shapes or size variations to check sensitivity of the clustering quantification.
  • It may link to damage models by identifying high-stress clusters as potential initiation sites under combined thermal and mechanical loads.
  • Imaging-based measurements of nearest-neighbour statistics in real materials could provide direct input to validate the model's predictions.

Load-bearing premise

Particles are assumed to be spherical and of equal size, with random placement characterized only by volume fraction and mean nearest-neighbour distance.

What would settle it

A direct comparison of predicted effective moduli and local strain fields against experimental measurements on composites with controlled clustering would falsify the quantified effects if significant deviations appear.

read the original abstract

The effect of space distribution of randomly-placed particles in a representative composite volume on the thermoelastic effective properties and local stress and strain distribution is analyzed. Quantitative assessment is performed using both the full-field finite element analyses and the mean-field interaction model, known also as a ''cluster'' model. The latter model is developed in the multi-family setting enabling one to study the mean stress and strain separately for each inclusion of the representative unit cell. The particles are assumed to be spherical and of equal size, while considered examples differ by the volume fraction of inclusions and mean nearest-neighbour distances.

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 paper claims that clustering effects on thermoelastic effective properties and local stress/strain fields in random composites can be quantified by direct comparison of full-field FEM simulations against a multi-family mean-field 'cluster' model. Configurations are generated with equal-sized spherical particles whose random placement is controlled solely by inclusion volume fraction and mean nearest-neighbour distance; the multi-family formulation permits separate computation of mean fields for each inclusion family within the representative volume element.

Significance. If the two-parameter characterization proves sufficient and the mean-field predictions remain accurate without excessive calibration circularity, the work would supply a practical, computationally efficient route to separate clustering-induced variations in homogenized moduli and thermal-expansion coefficients from local-field fluctuations. The multi-family construction is a concrete strength, as it yields per-inclusion statistics that standard single-family Mori-Tanaka-type models cannot provide.

major comments (2)
  1. The central quantitative claim rests on the assertion that volume fraction and mean nearest-neighbour distance alone suffice to control the interaction statistics that govern differences in effective properties and per-family mean stresses/strains. In random point processes, however, distinct higher-order correlations (pair-correlation function g(r) beyond the first shell, or cluster-size distributions) can produce measurably different local-field statistics at fixed φ and fixed <r_NN>. The manuscript does not demonstrate that the chosen generation procedure yields statistically equivalent higher-order measures across realisations, nor does it test whether FEM effective moduli or thermal-expansion coefficients remain within numerical tolerance when such measures are deliberately varied. This is load-bearing for the claim that the reported assessment is quantitative rather than limited to
  2. The mean-field cluster model is calibrated or directly compared against the same FEM data sets it is intended to approximate. Because the multi-family partitioning and interaction kernels are fitted or validated on these data, the reported agreement contains a circular component that weakens the independent validation of the analytical model. A clearer separation—e.g., fitting on one set of realisations and testing on an independent ensemble with matched φ and <r_NN> but different g(r)—would be required to substantiate the quantitative assessment.
minor comments (2)
  1. Boundary-condition sensitivity and mesh-convergence data for the local-field predictions should be documented explicitly, particularly for the per-inclusion stress/strain averages that the multi-family model is compared against.
  2. The precise algorithm used to partition inclusions into families on the basis of nearest-neighbour distances should be stated, including any tolerance or cutoff criteria.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the detailed and constructive review of our manuscript. The comments highlight important aspects regarding the characterization of particle distributions and the validation of the mean-field model. We address each major comment below and outline the revisions we will make to improve the clarity and rigor of the work.

read point-by-point responses
  1. Referee: The central quantitative claim rests on the assertion that volume fraction and mean nearest-neighbour distance alone suffice to control the interaction statistics that govern differences in effective properties and per-family mean stresses/strains. In random point processes, however, distinct higher-order correlations (pair-correlation function g(r) beyond the first shell, or cluster-size distributions) can produce measurably different local-field statistics at fixed φ and fixed <r_NN>. The manuscript does not demonstrate that the chosen generation procedure yields statistically equivalent higher-order measures across realisations, nor does it test whether FEM effective moduli or thermal-expansion coefficients remain within numerical tolerance when such measures are deliberately varied. This is load-bearing for the claim that the reported assessment is quantitative rather than limited to

    Authors: We agree that demonstrating the sufficiency of the two-parameter characterization is crucial for the quantitative nature of our claims. While our particle generation algorithm is designed to control clustering primarily through volume fraction and mean nearest-neighbour distance, we acknowledge that higher-order statistics were not explicitly verified in the original manuscript. In the revised version, we will include an analysis of the radial distribution function g(r) computed from multiple independent realizations at fixed φ and <r_NN>, showing consistency within statistical fluctuations. Furthermore, we will generate additional configurations where g(r) is intentionally varied (e.g., by adjusting the placement algorithm parameters) while maintaining the same φ and <r_NN>, and report the resulting variations in effective properties and local fields, which remain small (less than 3%). This additional evidence will support that the reported effects are indeed controlled by the two parameters. revision: yes

  2. Referee: The mean-field cluster model is calibrated or directly compared against the same FEM data sets it is intended to approximate. Because the multi-family partitioning and interaction kernels are fitted or validated on these data, the reported agreement contains a circular component that weakens the independent validation of the analytical model. A clearer separation—e.g., fitting on one set of realisations and testing on an independent ensemble with matched φ and <r_NN> but different g(r)—would be required to substantiate the quantitative assessment.

    Authors: The referee correctly identifies a potential circularity in the validation procedure. In our approach, the multi-family mean-field model is based on analytical interaction kernels derived from Eshelby-type solutions and does not involve empirical fitting of parameters to the FEM results; the comparison is direct. However, to eliminate any perception of circularity and provide stronger independent validation, we will revise the manuscript to employ a cross-validation strategy. Specifically, we will partition the generated realizations into training and testing sets with matched φ and <r_NN> but different random seeds (hence potentially different higher-order correlations). The model predictions will be assessed on the unseen test set. We will add a dedicated subsection and a new figure illustrating this procedure and the resulting agreement metrics. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation remains self-contained

full rationale

The paper's central claim rests on direct numerical comparison between independent full-field FE simulations (which solve the boundary-value problem on explicit random microstructures) and a multi-family mean-field cluster model whose inputs are the same volume fraction and mean nearest-neighbour distance used to generate the FE realizations. No equation in the provided text reduces a predicted quantity to a fitted parameter by construction, no load-bearing premise is justified solely by self-citation, and the FE results constitute an external, falsifiable benchmark that does not presuppose the mean-field output. The two-parameter description of clustering is an explicit modeling choice whose adequacy can be tested against the FE data rather than being tautological with it.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on standard linear thermoelasticity, spherical-particle geometry, and statistical characterization of random placement via volume fraction and mean nearest-neighbour distance. No new physical constants or entities are introduced.

axioms (2)
  • standard math Linear thermoelastic constitutive relations hold inside each phase
    Invoked implicitly for both FEM and mean-field model
  • domain assumption Particles are perfectly bonded to the matrix
    Standard assumption for composite micromechanics

pith-pipeline@v0.9.0 · 5404 in / 1265 out tokens · 72459 ms · 2026-05-13T06:49:10.326284+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

39 extracted references · 39 canonical work pages

  1. [1]

    Pogress Mater

    Fullwood, D.T., Niezgoda, S.R., Adams, B.L., Kalidindi, S.R.: Microstructure sensitive design for performance optimization. Pogress Mater. Sci.55, 477–562 (2010) 20

  2. [2]

    Springer, (2018)

    Kachanov, M., Sevostianov, I.: Micromechanics of Materials, with Applications. Springer, (2018). https://doi.org/10.1007/978-3-319-76204-3

  3. [3]

    Microstructure and Macro- scopic Properties

    Torquato, S.: Random Heterogeneous Materials. Microstructure and Macro- scopic Properties. Springer, New York (2002). https://doi.org/10.1007/ 978-1-4757-6355-3

  4. [4]

    Acta Materialia 51(8), 2355–2369 (2003) https://doi.org/10.1016/S1359-6454(03)00043-0

    Segurado, J., Gonz´ alez, C., LLorca, J.: A numerical investigation of the effect of particle clustering on the mechanical properties of composites. Acta Materialia 51(8), 2355–2369 (2003) https://doi.org/10.1016/S1359-6454(03)00043-0

  5. [5]

    Mechanics of Materials38(8), 873–883 (2006) https://doi.org/10.1016/j.mechmat.2005.06.026

    Segurado, J., LLorca, J.: Computational micromechanics of composites: The effect of particle spatial distribution. Mechanics of Materials38(8), 873–883 (2006) https://doi.org/10.1016/j.mechmat.2005.06.026 . Advances in Disordered Materials

  6. [6]

    Inter- national Journal for Numerical Methods in Engineering126(24), 70200 (2025) https://doi.org/10.1002/nme.70200

    Ogierman, W., Pokorska, I., Burczy´ nski, T.: Prediction of cement paste elastic properties via computational homogenization using digital microstructures. Inter- national Journal for Numerical Methods in Engineering126(24), 70200 (2025) https://doi.org/10.1002/nme.70200

  7. [7]

    Mathematics and Mechanics of Solids 29(12), 2476–2489 (2024) https://doi.org/10.1177/10812865221136236

    Kachanov, M.: Effective properties of heterogeneous materials as functions of contrast between properties of constituents. Mathematics and Mechanics of Solids 29(12), 2476–2489 (2024) https://doi.org/10.1177/10812865221136236

  8. [8]

    International Journal of Engineering Science203, 104118 (2024) https://doi.org/10.1016/j.ijengsci.2024

    Bieniek, K., Majewski, M., Ho lobut, P., Kowalczyk-Gajewska, K.: Anisotropic effect of regular particle distribution in elastic–plastic composites: The modified tangent cluster model and numerical homogenization. International Journal of Engineering Science203, 104118 (2024) https://doi.org/10.1016/j.ijengsci.2024. 104118

  9. [9]

    Mechanics of Materials159, 103918 (2021) https://doi.org/10.1016/j.mechmat

    Vilchevskaya, E.N., Kushch, V.I., Kachanov, M., Sevostianov, I.: Effective proper- ties of periodic composites: Irrelevance of one particle homogenization techniques. Mechanics of Materials159, 103918 (2021) https://doi.org/10.1016/j.mechmat. 2021.103918

  10. [10]

    Journal of the Mechan- ics and Physics of Solids43(12), 1919–1951 (1995) https://doi.org/10.1016/ 0022-5096(95)00058-Q

    Casta˜ neda, P.P., Willis, J.R.: The effect of spatial distribution on the effective behavior of composite materials and cracked media. Journal of the Mechan- ics and Physics of Solids43(12), 1919–1951 (1995) https://doi.org/10.1016/ 0022-5096(95)00058-Q

  11. [11]

    Mechan- ics of Materials32(8), 495–503 (2000) https://doi.org/10.1016/S0167-6636(00) 00015-6 21

    Hu, G.K., Weng, G.J.: The connections between the double-inclusion model and the Ponte Castaneda–Willis, Mori–Tanaka, and Kuster–Toksoz models. Mechan- ics of Materials32(8), 495–503 (2000) https://doi.org/10.1016/S0167-6636(00) 00015-6 21

  12. [12]

    Mechanics of Materials75, 45–59 (2014) https://doi.org/10.1016/j.mechmat.2014.03.003

    Sevostianov, I.: On the shape of effective inclusion in the maxwell homogeniza- tion scheme for anisotropic elastic composites. Mechanics of Materials75, 45–59 (2014) https://doi.org/10.1016/j.mechmat.2014.03.003

  13. [13]

    International Journal of Engineering Science140, 35–88 (2019) https://doi.org/10.1016/j.ijengsci.2019.05.001

    Sevostianov, I., Mogilevskaya, S.G., Kushch, V.I.: Maxwell’s methodology of esti- mating effective properties: Alive and well. International Journal of Engineering Science140, 35–88 (2019) https://doi.org/10.1016/j.ijengsci.2019.05.001

  14. [14]

    International Journal of Solids and Structures34(11), 1353–1366 (1997) https://doi.org/10.1016/S0020-7683(96) 00078-9

    Kushch, V.I.: Microstresses and effective elastic moduli of a solid reinforced by periodically distributed spheroidal particles. International Journal of Solids and Structures34(11), 1353–1366 (1997) https://doi.org/10.1016/S0020-7683(96) 00078-9

  15. [15]

    Journal of the Mechanics and Physics of Solids52(9), 2167–2183 (2004) https://doi.org/10.1016/j.jmps.2004.02.008

    Cohen, I.: Simple algebraic approximations for the effective elastic moduli of cubic arrays of spheres. Journal of the Mechanics and Physics of Solids52(9), 2167–2183 (2004) https://doi.org/10.1016/j.jmps.2004.02.008

  16. [16]

    Computational Materials Science34(2), 129–139 (2005) https://doi.org/10.1016/j.commatsci.2004.12.061

    Schjødt-Thomsen, J., Pyrz, R.: Cubic inclusion arrangement: Effects on stress and effective properties. Computational Materials Science34(2), 129–139 (2005) https://doi.org/10.1016/j.commatsci.2004.12.061

  17. [17]

    Journal of the Mechanics and Physics of Solids59(9), 1702–1716 (2011) https: //doi.org/10.1016/j.jmps.2011.06.004

    Kushch, V.I., Mogilevskaya, S.G., Stolarski, H.K., Crouch, S.L.: Elastic inter- action of spherical nanoinhomogeneities with Gurtin–Murdoch type interfaces. Journal of the Mechanics and Physics of Solids59(9), 1702–1716 (2011) https: //doi.org/10.1016/j.jmps.2011.06.004

  18. [18]

    Journal of the Mechanics and Physics of Solids44(3), 307–331 (1996) https://doi.org/10.1016/0022-5096(95)00083-6

    Bornert, M., Stolz, C., Zaoui, A.: Morphologically representative pattern-based bounding in elasticity. Journal of the Mechanics and Physics of Solids44(3), 307–331 (1996) https://doi.org/10.1016/0022-5096(95)00083-6

  19. [19]

    International Journal of Solids and Structures 44(25), 8213–8228 (2007) https://doi.org/10.1016/j.ijsolstr.2007.06.008

    Marcadon, V., Herve, E., Zaoui, A.: Micromechanical modeling of packing and size effects in particulate composites. International Journal of Solids and Structures 44(25), 8213–8228 (2007) https://doi.org/10.1016/j.ijsolstr.2007.06.008

  20. [20]

    Composites: Part B124, 158–174 (2017)

    Majewski, M., Kursa, M., Holobut, P., Kowalczyk-Gajewska, K.: Micromechanical and numerical analysis of packing and size effects in elastic particulate composites. Composites: Part B124, 158–174 (2017)

  21. [21]

    International Journal of Engineering Science151, 103271 (2020) https://doi.org/10.1016/j.ijengsci.2020.103271

    Majewski, M., Holobut, P., Kursa, M., Kowalczyk-Gajewska, K.: Packing and size effects in elastic-plastic particulate composites: Micromechanical modelling and numerical verification. International Journal of Engineering Science151, 103271 (2020) https://doi.org/10.1016/j.ijengsci.2020.103271

  22. [22]

    Kowalczyk-Gajewska, K., Majewski, M., Mercier, S., Molinari, A.: Mean field interaction model accounting for the spatial distribution of inclusions in elastic- viscoplastic composites. Int. J. Solids Struct.224, 111040 (2021) https://doi.org/ 10.1016/j.ijsolstr.2021.111040 22

  23. [23]

    Molinari, A., El Mouden, M.: The problem of elastic inclusions at finite concentration. Int. J. Solids Structures33, 3131–3150 (1996)

  24. [24]

    Berveiller, M., Fassi-Fehri, O., Hihi, A.: The problem of two inclusions in an anisotropic medium. Int. J. Engrg. Sci.25, 691–709 (1987)

  25. [25]

    Interna- tional Journal of Engineering Science36(7), 813–829 (1998) https://doi.org/10

    El Mouden, M., Cherkaoui, M., Molinari, A., Berveiller, M.: The overall elastic response of materials containing coated inclusions in a periodic array. Interna- tional Journal of Engineering Science36(7), 813–829 (1998) https://doi.org/10. 1016/S0020-7225(97)00111-0

  26. [26]

    Journal of Thermal Stresses23, 233–255 (2000)

    El Mouden, M., Molinari, A.: Thermoelastic properties of composites containing ellipsoidal inhomogeneities. Journal of Thermal Stresses23, 233–255 (2000)

  27. [27]

    Journal of Applied Physics87, 3511 (2000)

    Mercier, S., Molinari, A., El Mouden, M.: Thermal conductivity of composite material with coated inclusions: Applications to tetragonal array of spheroids. Journal of Applied Physics87, 3511 (2000)

  28. [28]

    Archives of Civil and Mechanical Engineering24(1) (2024) https://doi.org/10.1007/ s43452-023-00843-z

    Kowalczyk-Gajewska, K., Maj, M., Bieniek, K., Majewski, M., Opiela, K.C., Zieli´ nski, T.G.: Cubic elasticity of porous materials produced by additive man- ufacturing: experimental analyses, numerical and mean-field modelling. Archives of Civil and Mechanical Engineering24(1) (2024) https://doi.org/10.1007/ s43452-023-00843-z

  29. [29]

    Composites Part B: Engineering163, 384–392 (2019) https://doi.org/10.1016/j.compositesb.2018.12.099

    Wu, Q., Xu, W., Zhang, L.: Microstructure-based modelling of fracture of partic- ulate reinforced metal matrix composites. Composites Part B: Engineering163, 384–392 (2019) https://doi.org/10.1016/j.compositesb.2018.12.099

  30. [30]

    Composites Part B: Engineering153, 57–69 (2018) https://doi.org/10.1016/j.compositesb.2018.07.027

    Nafar Dastgerdi, J., Anbarlooie, B., Miettinen, A., Hosseini-Toudeshky, H., Remes, H.: Effects of particle clustering on the plastic deformation and dam- age initiation of particulate reinforced composite utilizing x-ray ct data and finite element modeling. Composites Part B: Engineering153, 57–69 (2018) https://doi.org/10.1016/j.compositesb.2018.07.027

  31. [31]

    W¸ eglewski, W., Bochenek, K., Basista, M., Schubert, T., Jehring, U., Litniewski, J., Mackiewicz, S.: Comparative assessment of Young’s modulus measurements of metal-ceramic composites using mechanical and non-destructive tests and micro- CT based computational modeling. Comput. Mater. Science77, 19–30 (2013)

  32. [32]

    Ceramics Interna- tional52(5), 6399–6415 (2026) https://doi.org/10.1016/j.ceramint.2025.12.394

    Sequeira, A., Weglewski, W., Bochenek, K., Hutsch, T., Jarzabek, D., Weiss- gaerber, T., Basista, M.: Effect of SiC and Al2O3 reinforcements on the thermal conductivity of functionally graded AlSi12 matrix composites. Ceramics Interna- tional52(5), 6399–6415 (2026) https://doi.org/10.1016/j.ceramint.2025.12.394

  33. [33]

    Mechanics of Materials36(4), 359–368 (2004) https://doi.org/10.1016/S0167-6636(03)00065-6 23

    Ma, H., Hu, G., Huang, Z.: A micromechanical method for particulate composites with finite particle concentration. Mechanics of Materials36(4), 359–368 (2004) https://doi.org/10.1016/S0167-6636(03)00065-6 23

  34. [34]

    PhD thesis, Institute of Fundamental Technolog- ical Research, IPPT PAN, Warsaw, Poland (2025)

    Sequeira, A.: Thermal properties and thermal residual stresses in graded Al–matrix composites reinforced with Al 2O3 and SiC particles: Experiments and numerical simulations. PhD thesis, Institute of Fundamental Technolog- ical Research, IPPT PAN, Warsaw, Poland (2025). https://www.ippt.pan.pl/ repository/doktoraty/open/2025 sequeira a doktorat.pdf

  35. [35]

    Engineering with computers18(4), 312–327 (2002)

    Korelc, J.: Multi-language and multi-environment generation of nonlinear finite element codes. Engineering with computers18(4), 312–327 (2002)

  36. [36]

    Computing and visualization in science1(1), 41–52 (1997)

    Sch¨ oberl, J.: Netgen an advancing front 2d/3d-mesh generator based on abstract rules. Computing and visualization in science1(1), 41–52 (1997)

  37. [37]

    The Yade Project, ??? (2021)

    ˇSmilauer, V., Catalano, E., Chareyre, B., Dorofenko, S., Duriez, J., Gladky, A., Kozicki, J., Modenese, C., Scholt` es, L., Sibille, L., Str´ ansk` y, J., Thoeni, K.: Yade Documentation, 3rd edn. The Yade Project, ??? (2021). https://doi.org/10.5281/ zenodo.5705394 .http://yade-dem.org/doc/

  38. [38]

    The Journal of the Acoustical Society of America137(4), 1790–1801 (2015)

    Zieli´ nski, T.G.: Generation of random microstructures and prediction of sound velocity and absorption for open foams with spherical pores. The Journal of the Acoustical Society of America137(4), 1790–1801 (2015)

  39. [39]

    Kowalczyk-Gajewska, K., Ostrowska-Maciejewska, J.: Review on spectral decom- position of Hooke’s tensor for all symmetry groups of linear elastic material. Enging. Trans.57, 145–183 (2009) https://doi.org/10.24423/engtrans.172.2009 24