REVIEW 2 major objections 51 references
Reviewed by Pith at T0; open to challenge.
T0 means a machine referee read the full paper against a public rubric. The mark states how deep the mechanical check went, never who wrote it. the ladder, T0–T4 →
T0 review · grok-4.3
FUSE couples independent FEBio models through structured field exchange without modifying solvers.
2026-07-03 01:20 UTC pith:UJVBT3K3
load-bearing objection FUSE is a practical plugin that reuses FEBio's data maps for partitioned multiphysics coupling in time-decoupled problems, with examples in cartilage and bone. the 2 major comments →
FUSE: A Partitioned Field-Exchange Framework for Coupling Physics Simulations in FEBio
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
FUSE is a partitioned coupling plugin for FEBio that enables communication between independently defined models through bidirectional field transfer, spatial mapping, and filtered exchange using existing data maps and outputs, following a time-decoupled strategy with a primary model on longer scales and secondary models on shorter horizons, demonstrated to reproduce reference solutions in coupled mechanics-chemistry and mechanics-biology problems.
What carries the argument
The FUSE plugin implementing time-decoupled partitioned coupling via FEBio's data maps, output fields, and user-specified filters.
Load-bearing premise
Problems are best addressed by time-decoupled partitioned coupling where fast mechanics influence slower biology or chemistry, relying on existing FEBio data maps without solver changes.
What would settle it
Running a test case with tightly coupled fast bidirectional interactions and observing if the partitioned FUSE solution deviates from a monolithic reference beyond acceptable error.
If this is right
- Coupled multiphysics models can be assembled from independent components without custom pipelines.
- Bidirectional transfer and spatial field mapping are handled for problems spanning mechanical and slower processes.
- Workflows remain maintainable by keeping physics implementations separate from coupling logic.
- Applications include mechanical-chemical coupling in cartilage and mechanical-biological feedback in bone healing.
Where Pith is reading between the lines
- The framework could extend to other simulation environments if similar data map interfaces exist.
- Strong coupling on overlapping timescales might require additional synchronization beyond the current decoupled approach.
- Standardizing the field exchange could facilitate community-shared model components for biomechanics.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript presents FUSE, a partitioned coupling plugin for FEBio that enables separately defined models to exchange fields (bidirectional, spatial mapping, filtered) via existing data maps and output fields without solver modifications. It employs a time-decoupled strategy with a primary model on the longer timescale repeatedly initializing and updating secondary models on shorter horizons. The central claim is that this reproduces reference coupled solutions, demonstrated on cartilage mechanical-chemical degradation and bone healing with tissue-mechanical feedback.
Significance. If the reproduction claims hold with supporting evidence, FUSE offers a maintainable mechanism for multiphysics workflows in FEBio by isolating coupling logic from physics implementations. This addresses a practical barrier to model reuse in computational biomechanics involving mechanics, transport, chemistry, and biology, leveraging pre-existing FEBio infrastructure rather than requiring new solver code.
major comments (2)
- [Abstract] Abstract: the central claim that 'the framework was able to reproduce reference coupled solutions' while handling bidirectional transfer, spatial field mapping, and filtered exchange supplies no quantitative metrics, error analysis, convergence rates, or specific test-case comparisons, which is load-bearing for validating the partitioned exchange approach.
- [Abstract] Abstract: the time-decoupled primary/secondary initialization pattern is explicitly scoped to problems where fast mechanics drive slower evolution, but no analysis or test is provided on stability, accuracy loss, or failure modes when this separation of timescales does not hold.
Simulated Author's Rebuttal
We thank the referee for the constructive comments. We address each major comment below and indicate the revisions we will make to strengthen the manuscript.
read point-by-point responses
-
Referee: [Abstract] Abstract: the central claim that 'the framework was able to reproduce reference coupled solutions' while handling bidirectional transfer, spatial field mapping, and filtered exchange supplies no quantitative metrics, error analysis, convergence rates, or specific test-case comparisons, which is load-bearing for validating the partitioned exchange approach.
Authors: We agree that the abstract would be strengthened by including quantitative support for the reproduction claim. The full manuscript contains L2-norm error comparisons and maximum deviation metrics for the cartilage degradation and bone-healing cases, but these are not summarized in the abstract. In the revised version we will add a concise statement reporting the observed error levels (e.g., relative L2 errors below 2 % for the primary fields) and the specific test-case comparisons performed. revision: yes
-
Referee: [Abstract] Abstract: the time-decoupled primary/secondary initialization pattern is explicitly scoped to problems where fast mechanics drive slower evolution, but no analysis or test is provided on stability, accuracy loss, or failure modes when this separation of timescales does not hold.
Authors: The framework is deliberately scoped to problems possessing a clear separation of timescales, as stated in the abstract, introduction, and methods. We do not claim applicability or provide validation for cases lacking this separation, because those situations would require a different (e.g., tightly coupled) strategy. We will add a short limitations paragraph in the discussion section that explicitly notes this scope restriction and advises users to verify timescale assumptions before applying the method. revision: partial
Circularity Check
No significant circularity
full rationale
This is a software-framework description paper with no derivations, equations, fitted parameters, or predictions that could reduce to inputs by construction. The central claims rest on empirical reproduction of reference solutions using pre-existing FEBio data maps, output fields, and filters; the time-decoupled partitioned strategy is explicitly scoped to a class of problems rather than derived from self-referential premises. No load-bearing self-citations, ansatzes, or uniqueness theorems appear in the construction.
Axiom & Free-Parameter Ledger
read the original abstract
Computational biomechanics increasingly requires models that combine mechanics, transport, chemistry, and biological regulation across different spatial and temporal scales. The FEBio simulation software provides extensive open-source capabilities for modeling these processes using monolithic approaches. However, assembling independently developed physics models into reproducible coupled workflows remains challenging. Existing approaches often require custom scripts or external software pipelines, which can limit model reuse and complicate development. We present FUSE, the FEBio Unified Simulation and Exchange framework, a partitioned coupling plugin that enables separately defined FEBio models to communicate through structured field exchange. FUSE is designed for problems that are best solved independently, particularly when fast mechanical responses influence slower biological or chemical evolution. The framework uses a time-decoupled strategy in which a primary model advances on the longer time scale, while one or more secondary models are repeatedly initialized, supplied with updated fields, solved over shorter time horizons, with results returned to the primary model. Field exchange utilizes existing FEBio data maps, output fields, and user-specified filters, allowing coupled workflows to be constructed without modifying the underlying solvers. The framework was able to reproduce reference coupled solutions while handling bidirectional transfer, spatial field mapping, and filtered exchange of model variables. Example applications demonstrated coupling between mechanical loading and chemical degradation in injured cartilage and interaction between biological tissue formation and mechanical feedback during bone healing. By separating coupling logic from physics implementation, FUSE provides a practical mechanism for building maintainable multiphysics workflows within FEBio.
Reference graph
Works this paper leans on
-
[1]
D. Ambrosi, G.A. Ateshian, E.M. Arruda, S.C. Cowin, J. Dumais, A. Goriely, G.A. Holzapfel, J.D. Humphrey, R. Kemkemer, E. Kuhl, J.E. Olberding, L.A. Taber, K. Garikipati, Perspectives on biological growth and remodeling, J Mech Phys Solids, 59 (2011) 863-883
work page 2011
-
[2]
G.A. Ateshian, J.D. Humphrey, Continuum mixture models of biological growth and remodeling: past successes and future opportunities, Annu Rev Biomed Eng, 14 (2012) 97-111
work page 2012
- [3]
-
[4]
Y .F. Guo, M.R.K. Mofrad, A.B. Tepole, On modeling the multiscale mechanobiology of soft tissues: Challenges and progress, Biophys Rev-Us, 3 (2022)
work page 2022
-
[5]
G.A. Ateshian, On the theory of reactive mixtures for modeling biological growth, Biomech Model Mechanobiol, 6 (2007) 423-445
work page 2007
-
[6]
G.A. Ateshian, R.J. Nims, S. Maas, J.A. Weiss, Computational modeling of chemical reactions and interstitial growth and remodeling involving charged solutes and solid-bound molecules, Biomech Model Mechanobiol, 13 (2014) 1105-1120
work page 2014
-
[7]
M.M. Rahman, P.N. Watton, C.P. Neu, D.M. Pierce, A chemo-mechano-biological modeling framework for cartilage evolving in health, disease, injury, and treatment, Comput Meth Prog Bio, 231 (2023)
work page 2023
-
[8]
L.T. Edgar, S.A. Maas, J.E. Guilkey, J.A. Weiss, A coupled model of neovessel growth and matrix mechanics describes and predicts angiogenesis in vitro, Biomech Model Mechanobiol, 14 (2015) 767-782
work page 2015
-
[9]
M.S. Ghiasi, J.E. Chen, E.K. Rodriguez, A. Vaziri, A. Nazarian, Computational modeling of human bone fracture healing affected by different conditions of initial healing stage, Bmc Musculoskel Dis, 20 (2019)
work page 2019
-
[10]
D.E. Keyes, L.C. McInnes, C. Woodward, W. Gropp, E. Myra, M. Pernice, J. Bell, J. Brown, A. Clo, J. Connors, E. Constantinescu, D. Estep, K. Evans, C. Farhat, A. Hakim, G. Hammond, G. Hansen, J. Hill, T. Isaac, X.M. Jiao, K. Jordan, D. Kaushik, E. Kaxiras, A. Koniges, K. Lee, A. Lott, Q.M. Lu, J. Magerlein, R. Maxwell, M. McCourt, M. Mehl, R. Pawlowski, A...
work page 2013
-
[11]
H.J. Bungartz, F. Lindner, B. Gatzhammer, M. Mehl, K. Scheufele, A. Shukaev, B. Uekermann, preCICE - A fully parallel library for multi-physics surface coupling, Comput Fluids, 141 (2016) 250-258
work page 2016
-
[12]
C.A. Felippa, K.C. Park, C. Farhat, Partitioned analysis of coupled mechanical systems, Comput Method Appl M, 190 (2001) 3247-3270
work page 2001
- [13]
-
[14]
G.A. Ateshian, M.R. Herron, S.A. Maas, J.A. Weiss, Impact of FEBio, 2026
work page 2026
-
[15]
G.A. Ateshian, S. Maas, J.A. Weiss, Multiphasic finite element framework for modeling hydrated mixtures with multiple neutral and charged solutes, J Biomech Eng, 135 (2013) 111001
work page 2013
-
[16]
G.A. Ateshian, J.J. Shim, S.A. Maas, J.A. Weiss, Finite Element Framework for Computational Fluid Dynamics in FEBio, J Biomech Eng, 140 (2018) 0210011-02100117
work page 2018
- [17]
- [18]
- [19]
-
[20]
S.A. Labelle, A.M. Poulson, S.A. Maas, A. Rauff, G.A. Ateshian, J.A. Weiss, Spatial Configurations of 3D Extracellular Matrix Collagen Density and Anisotropy Simultaneously Guide Angiogenesis, Plos Comput Biol, 19 (2023)
work page 2023
-
[21]
S.A. Labelle, M.S. Sadrabadi, S. Baek, M.R.K. Mofrad, J.A. Weiss, A. Arzani, Multiscale Kinematic Growth Coupled With Mechanosensitive Systems Biology in Open-Source Software, J Biomech Eng-T Asme, 147 (2025)
work page 2025
-
[22]
(!!! INV ALID CITATION !!! [22, 23])
-
[23]
S. Roberts, J. Loffeld, A. Sarshar, C.S. Woodward, A. Sandu, Implicit Multirate GARK Methods, J Sci Comput, 87 (2021)
work page 2021
- [24]
-
[25]
G. Chourdakis, K. Davis, B. Rodenberg, M. Schulte, F. Simonis, B. Uekermann, G. Abrams, H.J. Bungartz, L. Cheung Yau, I. Desai, K. Eder, R. Hertrich, F. Lindner, A. Rusch, D. Sashko, D. Schneider, A. Totounferoush, D. V olland, P. V ollmer, O.Z. Koseomur, preCICE v2: A sustainable and user-friendly coupling library, Open Res Eur, 2 (2022) 51
work page 2022
-
[26]
A.S.A. Eskelinen, P. Tanska, C. Florea, G.A. Orozco, P. Julkunen, A.J. Grodzinsky, R.K. Korhonen, Mechanobiological model for simulation of injured cartilage degradation via pro- inflammatory cytokines and mechanical stimulus, Plos Comput Biol, 16 (2020) e1007998
work page 2020
-
[27]
S.M. Finley, D.S. Brodke, N.T. Spina, C.A. DeDen, B.J. Ellis, FEBio finite element models of the human lumbar spine, Comput Methods Biomech Biomed Engin, 21 (2018) 444-452
work page 2018
-
[28]
A.R. MacLeod, H. Rose, H.S. Gill, A Validated Open-Source Multisolver Fourth-Generation Composite Femur Model, J Biomech Eng, 138 (2016)
work page 2016
- [29]
-
[30]
S. Maas, F. Muhib, J.A. Weiss, FEBio Fuse Github Repository, 2026
work page 2026
-
[31]
Maas, FEBio Plugins, FEBio Developer's Manual, 2026
S. Maas, FEBio Plugins, FEBio Developer's Manual, 2026
work page 2026
-
[32]
Maas, FEBio Task Plugins, FEBio Developer's Manual, 2026
S. Maas, FEBio Task Plugins, FEBio Developer's Manual, 2026
work page 2026
-
[33]
S. Maas, The FEBioChem plugin for FEBio: solving the reaction-diffusion equation in the context of a (non-deformable) mixture with multiple chemical species., 2026
work page 2026
-
[34]
H. Weinans, R. Huiskes, H.J. Grootenboer, The Behavior of Adaptive Bone-Remodeling Simulation-Models, J Biomech, 25 (1992) 1425-1441
work page 1992
-
[35]
D.R. Carter, W.C. Hayes, Bone compressive strength: the influence of density and strain rate, Science, 194 (1976) 1174-1176
work page 1976
-
[36]
D.R. Carter, W.C. Hayes, The compressive behavior of bone as a two-phase porous structure, J Bone Joint Surg Am, 59 (1977) 954-962
work page 1977
-
[37]
M.G. Mullender, R. Huiskes, H. Weinans, A Physiological Approach to the Simulation of Bone Remodeling as a Self-Organizational Control Process, J Biomech, 27 (1994) 1389-1394. 31
work page 1994
- [38]
-
[39]
A. Paz, G.A. Orozco, P. Tanska, J.J. Garcia, R.K. Korhonen, M.E. Mononen, A novel knee joint model in FEBio with inhomogeneous fibril-reinforced biphasic cartilage simulating tissue mechanical responses during gait: data from the osteoarthritis initiative, Comput Methods Biomech Biomed Engin, 26 (2023) 1353-1367
work page 2023
-
[40]
A.S.A. Eskelinen, M.E. Mononen, M.S. Venäläinen, R.K. Korhonen, P. Tanska, Maximum shear strain-based algorithm can predict proteoglycan loss in damaged articular cartilage, Biomech Model Mechan, 18 (2019) 753-778
work page 2019
-
[41]
A. Bailon-Plaza, M.C. van der Meulen, A mathematical framework to study the effects of growth factor influences on fracture healing, J Theor Biol, 212 (2001) 191-209
work page 2001
-
[42]
A. Bailon-Plaza, M.C. van der Meulen, Beneficial effects of moderate, early loading and adverse effects of delayed or excessive loading on bone healing, J Biomech, 36 (2003) 1069-1077
work page 2003
-
[43]
S. Ghimire, S. Miramini, M. Richardson, P. Mendis, L.H. Zhang, Role of Dynamic Loading on Early Stage of Bone Fracture Healing, Ann Biomed Eng, 46 (2018) 1768-1784
work page 2018
- [44]
-
[45]
E. Cory, A. Nazarian, V . Entezari, V . Vartanians, R. Müller, B.D. Snyder, Compressive axial mechanical properties of rat bone as functions of bone volume fraction, apparent density and micro-ct based mineral density, J Biomech, 43 (2010) 953-960
work page 2010
-
[46]
M.J. Gardner, M.C.H. van der Meulen, D. Demetrakopoulos, T.M. Wright, E.R. Myers, M.P. Bostrom, In vivo cyclic axial compression affects bone healing in the mouse tibia, J Orthop Res, 24 (2006) 1679-1686
work page 2006
-
[47]
K.E. Williams, F. Muhib, E. Dinh, K.E. Leguineche, A. Hajarizadeh, J.W. Rosenthal, T. Guyer, T. Seah, N.J. Willett, J.A. Weiss, R.E. Guldberg, Subject-specific multivariate modeling for regenerative rehabilitation of bone healing, Apl Bioeng, 9 (2025)
work page 2025
-
[48]
E.C. Jeffery, T.L.A. Mann, J.A. Pool, Z.Y . Zhao, S.J. Morrison, Bone marrow and periosteal skeletal stem/progenitor cells make distinct contributions to bone maintenance and repair, Cell Stem Cell, 29 (2022) 1547-+
work page 2022
-
[49]
O.D. de Lageneste, A. Julien, R. Abou-Khalil, G. Frangi, C. Carvalho, N. Cagnard, C. Cordier, S.J. Conway, C. Colnot, Periosteum contains skeletal stem cells with high bone regenerative potential controlled by Periostin, Nat Commun, 9 (2018)
work page 2018
-
[50]
E.R. Moore, Y .X. Zhu, H.S. Ryu, C.R. Jacobs, Periosteal progenitors contribute to load- induced bone formation in adult mice and require primary cilia to sense mechanical stimulation (vol 9, 190, 2018), Stem Cell Res Ther, 9 (2018)
work page 2018
-
[51]
Q. He, J.W. Lu, Q.S. Liang, L.T. Yao, T.F. Sun, H. Wang, M. Duffy, X. Jiang, Y .W. Lin, J.H. Lee, J. Ahn, N.A. Dyment, F. Mourkioti, J.D. Boerckel, L. Qin, fibroadipogenic progenitors in muscle are crucial for bone fracture repair, P Natl Acad Sci USA, 122 (2025)
work page 2025
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.