Recognition: 2 theorem links
· Lean TheoremQuantifying the effects of particle clustering in random thermoelastic composites -- numerical and mean-field analyses
Pith reviewed 2026-05-13 06:49 UTC · model grok-4.3
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.
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
- 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.
Referee Report
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)
- 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
- 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)
- 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.
- 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
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
-
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
-
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
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
axioms (2)
- standard math Linear thermoelastic constitutive relations hold inside each phase
- domain assumption Particles are perfectly bonded to the matrix
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
The mean-field interaction model (also known as a cluster model) was extended to the thermoelastic composite materials in [26]. ... tensors Γ_Ij ... polarisation tensor P0 ...
-
IndisputableMonolith/Foundation/DimensionForcing.leanalexander_duality_circle_linking unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
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
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
-
[1]
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
work page 2010
-
[2]
Kachanov, M., Sevostianov, I.: Micromechanics of Materials, with Applications. Springer, (2018). https://doi.org/10.1007/978-3-319-76204-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
work page 2002
-
[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]
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]
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]
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]
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]
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]
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
work page 1919
-
[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]
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]
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]
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]
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]
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]
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]
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]
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]
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)
work page 2017
-
[21]
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]
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]
Molinari, A., El Mouden, M.: The problem of elastic inclusions at finite concentration. Int. J. Solids Structures33, 3131–3150 (1996)
work page 1996
-
[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)
work page 1987
-
[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
work page 1998
-
[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)
work page 2000
-
[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)
work page 2000
-
[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
work page 2024
-
[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]
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]
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)
work page 2013
-
[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]
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]
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
work page 2025
-
[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)
work page 2002
-
[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)
work page 1997
-
[37]
ˇ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/
work page 2021
-
[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)
work page 2015
-
[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
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.