pith. sign in

arxiv: 2510.07552 · v1 · pith:ICUHHRAEnew · submitted 2025-10-08 · ❄️ cond-mat.stat-mech · physics.flu-dyn

Thermodynamically Consistent Continuum Theory of Magnetic Particles in High-Gradient Fields

Marko Tesanovic , Daniel M. Markiewitz , Marcus L. Popp , Martin Z. Bazant , Sonja Berensmeier This is my paper

Pith reviewed 2026-05-18 08:37 UTC · model grok-4.3

classification ❄️ cond-mat.stat-mech physics.flu-dyn
keywords magnetic particleshigh-gradient magnetic fieldscontinuum theoryfree-energy functionalparticle capturehomogenizationMason numbermagnetic separation
0
0 comments X

The pith

A thermodynamically consistent continuum theory derives magnetic particle transport and capture directly from a free-energy functional.

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

This paper builds a single set of equations for magnetic particles in strong field gradients by beginning with a free-energy expression that accounts for magnetic energy, particle mixing entropy, and steric repulsions. Homogenization then supplies a concentration-dependent susceptibility that feeds into coupled equations for magnetism, mass transport, and fluid momentum. No separate rules are needed to halt particle motion; instead, shielding, shaped deposits, and thin layers appear on their own. A phase diagram in Mason number organizes the behavior into thermodynamically controlled, transitional, and dynamically controlled regimes. The work supplies a predictive route for designing magnetic separation processes without patchwork empirical fixes.

Core claim

The central claim is that a free-energy functional coupling magnetic energy, entropic mixing, and steric interactions, together with homogenization theory for a concentration-dependent susceptibility, produces unified equations for magnetism, mass transport, and momentum balance. These equations allow field shielding, anisotropic deposition, and boundary-layer confinement to emerge naturally. Simulations recover canonical capture morphologies such as axially aligned plumes and crescent-shaped deposits, and organize captured mass data into three distinct regimes on a dimensionless Mason-number phase diagram.

What carries the argument

The free-energy functional that couples magnetic energy, entropic mixing, and steric interactions, combined with homogenization theory to obtain concentration-dependent susceptibility, which then closes the unified transport and momentum equations.

If this is right

  • Field shielding and anisotropic deposition appear without any added shut-off criteria.
  • Canonical morphologies such as axial plumes and crescent deposits arise across ranges of field strength and flow.
  • Captured particle mass organizes into three regimes on a Mason-number phase diagram.
  • The equations provide a platform for in-silico optimization of high-gradient magnetic separation.

Where Pith is reading between the lines

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

  • The same free-energy construction could be adapted to three-dimensional geometries for more complex industrial separators.
  • Digital-twin models of separation equipment might incorporate this continuum description as their core physics module.
  • Analogous thermodynamic closures could be developed for other field-driven colloidal systems, such as those using electric fields.

Load-bearing premise

The free-energy functional together with homogenization theory correctly captures the dominant interactions among magnetic particles, fluid, and field.

What would settle it

Systematic mismatch between the predicted capture morphologies or the locations of the three Mason-number regimes and direct experimental measurements in high-gradient magnetic fields would falsify the theory.

Figures

Figures reproduced from arXiv: 2510.07552 by Daniel M. Markiewitz, Marcus L. Popp, Marko Tesanovic, Martin Z. Bazant, Sonja Berensmeier.

Figure 1
Figure 1. Figure 1: FIG. 1. Hierarchical structure of the magnetic nanoparticle (MNP) suspension. ( [PITH_FULL_IMAGE:figures/full_fig_p008_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Steady-state magnetic particle concentration fields at [PITH_FULL_IMAGE:figures/full_fig_p014_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Dimensionless phase diagram for field–induced particle capture: Captured mass versus [PITH_FULL_IMAGE:figures/full_fig_p017_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Temporal evolution of magnetic particle accumulation (left) and magnetic field magnitude [PITH_FULL_IMAGE:figures/full_fig_p019_4.png] view at source ↗
read the original abstract

Magnetic particles underpin a broad range of technologies, from water purification and mineral processing to bioseparations and targeted drug delivery. The dynamics of magnetic particles in high-gradient magnetic fields-encompassing both their transport and eventual capture-arise from the coupled interplay of field-driven drift, fluid advection, and particle-field feedback. These processes remain poorly captured by existing models relying on empirical closures or discrete particle tracking. Here, we present a thermodynamically consistent continuum theory for collective magnetic particle transport and capture in high-gradient fields. The framework derives from a free-energy functional that couples magnetic energy, entropic mixing, and steric interactions, yielding a concentration-dependent susceptibility via homogenization theory. The resulting equations unify magnetism, mass transport, and momentum balances without ad hoc shut-off criteria, allowing field shielding, anisotropic deposition, and boundary-layer confinement to emerge naturally. Simulations predict canonical capture morphologies-axially aligned plumes, crescent-shaped deposits, and nonlinear shielding-across field strengths and flow regimes, consistent with trends reported in prior experimental and modeling studies. By organizing captured particle mass data into a dimensionless phase diagram based on the Mason number, we reveal three distinct regimes-thermodynamically controlled, transitional, and dynamically controlled. This perspective provides a predictive platform for in silico optimization and extension to three-dimensional geometries, and informing digital twin development for industrial-scale high-gradient magnetic separation processes.

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 presents a thermodynamically consistent continuum theory for the collective transport and capture of magnetic particles in high-gradient fields. Starting from a free-energy functional that incorporates magnetic energy, entropic mixing, and steric interactions, homogenization theory is used to obtain a concentration-dependent susceptibility. The resulting unified set of equations for magnetism, mass transport, and momentum balance is claimed to eliminate ad-hoc shut-off criteria, allowing field shielding, anisotropic deposition, and boundary-layer confinement to emerge naturally. Simulations are reported to reproduce canonical capture morphologies (axially aligned plumes, crescent-shaped deposits, nonlinear shielding) and to organize captured mass data into a Mason-number phase diagram that identifies three regimes: thermodynamically controlled, transitional, and dynamically controlled.

Significance. If the central derivation and its numerical implementation hold, the work supplies a predictive, closure-free framework for high-gradient magnetic separation that unifies multiple physical balances from a single free-energy principle. This could enable in-silico optimization of processes in water purification, mineral processing, bioseparations, and targeted drug delivery, and support digital-twin development for industrial-scale systems. The explicit avoidance of empirical shut-off criteria and the natural emergence of shielding and anisotropic effects constitute a clear methodological advance over existing discrete-particle or empirically closed models.

major comments (2)
  1. [Homogenization / effective susceptibility derivation] Homogenization section (derivation of concentration-dependent susceptibility): The effective susceptibility is obtained by homogenization that presupposes a clear separation between the macroscopic field-gradient length scale and the microscopic particle diameter or inter-particle spacing. In the high-gradient regime targeted by the paper, this separation can become marginal or absent, raising the possibility that the local averaging step introduces uncontrolled errors into the force term and the predicted shielding behavior. A quantitative assessment of the validity of the scale-separation assumption (e.g., via a comparison of homogenization predictions against direct micro-scale simulations at representative gradient strengths) is needed to substantiate the central claim of thermodynamic consistency and predictive accuracy.
  2. [Simulation results / Mason-number phase diagram] Numerical results and phase-diagram section: The manuscript states that simulations reproduce observed morphologies and that the Mason-number diagram reveals three distinct regimes, yet the support rests on qualitative trend agreement rather than quantitative metrics (capture efficiency, deposit mass versus field strength or flow rate, or direct comparison with experimental data sets). Without such benchmarks or an explicit statement of the numerical scheme, boundary conditions, and mesh convergence, it is difficult to judge whether the reported morphologies and regime boundaries are robust predictions or artifacts of discretization.
minor comments (2)
  1. [Model formulation] The free-energy functional is introduced without an explicit listing of all coefficients and their physical interpretations; adding a compact table of symbols and parameter values would improve readability.
  2. [Figures] Figure captions for the morphology and phase-diagram plots should state the precise values of the Mason number, susceptibility, and flow parameters used in each panel.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their insightful and constructive comments, which have helped clarify several aspects of our work. We respond to each major comment below and indicate the changes planned for the revised manuscript.

read point-by-point responses

  1. Referee: [Homogenization / effective susceptibility derivation] Homogenization section (derivation of concentration-dependent susceptibility): The effective susceptibility is obtained by homogenization that presupposes a clear separation between the macroscopic field-gradient length scale and the microscopic particle diameter or inter-particle spacing. In the high-gradient regime targeted by the paper, this separation can become marginal or absent, raising the possibility that the local averaging step introduces uncontrolled errors into the force term and the predicted shielding behavior. A quantitative assessment of the validity of the scale-separation assumption (e.g., via a comparison of homogenization predictions against direct micro-scale simulations at representative gradient strengths) is needed to substantiate the central claim of thermodynamic consistency and predictive accuracy.

    Authors: We thank the referee for this important observation on the scale-separation assumption underlying the homogenization step. The derivation follows the standard two-scale homogenization procedure for effective constitutive relations in heterogeneous media, which by construction yields a thermodynamically consistent concentration-dependent susceptibility from the underlying free-energy functional. While we acknowledge that the separation of scales can become marginal under the strongest gradients considered, the resulting continuum model still enforces thermodynamic consistency and permits the natural emergence of shielding and anisotropic effects without ad-hoc closures. In the revised manuscript we will add an explicit paragraph in the homogenization section that quantifies the scale-separation parameter for the simulated regimes and discusses its implications for the force term. A direct, quantitative benchmark against fully resolved micro-scale simulations would be a valuable extension but lies outside the scope of the present study; we will note this as a recommended direction for future work. revision: partial

  2. Referee: [Simulation results / Mason-number phase diagram] Numerical results and phase-diagram section: The manuscript states that simulations reproduce observed morphologies and that the Mason-number diagram reveals three distinct regimes, yet the support rests on qualitative trend agreement rather than quantitative metrics (capture efficiency, deposit mass versus field strength or flow rate, or direct comparison with experimental data sets). Without such benchmarks or an explicit statement of the numerical scheme, boundary conditions, and mesh convergence, it is difficult to judge whether the reported morphologies and regime boundaries are robust predictions or artifacts of discretization.

    Authors: We agree that additional quantitative support and numerical documentation will strengthen the presentation. In the revised manuscript we will expand the numerical-methods section to include a full description of the finite-element discretization, time-stepping scheme, boundary conditions, and mesh-convergence studies performed for representative cases. We will also add quantitative figures showing captured mass and capture efficiency as functions of Mason number, together with direct comparisons to experimental trends reported in the high-gradient magnetic separation literature. These additions will allow readers to assess the robustness of the three identified regimes more rigorously. revision: yes

Circularity Check

0 steps flagged

Derivation from free-energy functional via homogenization is self-contained

full rationale

The paper constructs its continuum equations from an explicit free-energy functional that incorporates magnetic energy, entropic mixing, and steric interactions, then applies homogenization theory to obtain a concentration-dependent susceptibility that enters the coupled transport and momentum balances. No equation or step in the provided description reduces a claimed prediction or result to a fitted parameter or prior self-citation by construction. The Mason number used for the phase diagram is a standard external dimensionless group, and the emergence of shielding and deposition morphologies is presented as a consequence of the derived equations rather than an input. The framework therefore remains independent of the target outputs and does not exhibit any of the enumerated circularity patterns.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

Information is limited to the abstract; the central construction rests on a free-energy functional whose explicit form and coefficients are not given, plus homogenization theory whose assumptions are not detailed.

free parameters (1)
  • coefficients in the free-energy functional
    Magnetic, entropic, and steric terms are invoked but their specific numerical values or fitting procedures are not stated in the abstract.
axioms (1)
  • domain assumption Homogenization theory yields a concentration-dependent magnetic susceptibility from the free-energy functional
    Invoked to close the magnetic response without further derivation shown in the abstract.

pith-pipeline@v0.9.0 · 5793 in / 1362 out tokens · 47136 ms · 2026-05-18T08:37:37.608851+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]

    This 6 hierarchical organization—spanning single nanoparticles, aggregates, and ultimately the ho- mogenized slurry—is depicted in Fig

    Macroscopic magnetic susceptibility In equilibrium, MNPs assemble into mesoscale aggregates containing roughly 5 % solid phase (Φ M N P ≈0.05), with the remaining volume occupied by entrained fluid [17]. This 6 hierarchical organization—spanning single nanoparticles, aggregates, and ultimately the ho- mogenized slurry—is depicted in Fig. 1. For modeling p...

    work page 2000

  2. [2]

    This formulation allows the framework to be augmented for phase-field models

    Magnetochemical thermodynamics We propose the following free-energy functional for an inhomogeneous magnetic suspen- sion: G= Z V dr g(˜c,∇ψ)− 1 2 µ0∇ψ· 1 +χ ef f a (˜c,∇ψ) ∇ψ (14) The first term g(˜c,∇ψ) is the enthalpy density, which depends on the concentration of mag- netic aggregates and theH-field. This formulation allows the framework to be augment...

    work page

  3. [3]

    Numerical simulation Numerical simulations were carried out in COMSOL Multiphysics ® 6.2 using built-in physics interfaces to solve the coupled magnetic, flow, and transport equations derived in the theory section. Three modules were combined: (i) Magnetic Fields, No Currents for magnetostatics, (ii) Laminar Flow for incompressible Stokes flow, and (iii) ...

    work page

  4. [4]

    Let Ω f be the fluid domain andc(x, t) the volumetric particle concentration [m −3]

    Captured MNP Mass In 2D we compute capture as anexcess inventory per unit depth. Let Ω f be the fluid domain andc(x, t) the volumetric particle concentration [m −3]. The instantaneous inventory is N(t) = Z Ωf c(x, t) dA[m −1]. Two otherwise identical simulations are run (same mesh,v 0,B 0, numerics): afield-oncase N on(t) and afield-offbaselineN off(t). T...

    work page

  5. [5]

    M. Lee, A. Shelke, S. Singh, J. Fan, P. Zaleski, and S. Afkhami, Physical Review E106 (2022), 10.1103/physreve.106.015104, publisher: American Physical Society (APS)

    work page doi:10.1103/physreve.106.015104 2022

  6. [6]

    Nithya, A

    R. Nithya, A. Thirunavukkarasu, A. B. Sathya, and R. Sivashankar, Environmental Chemistry Letters19, 1275 (2021), publisher: Springer Science and Business Media LLC

    work page 2021

  7. [7]

    Srivastava, G

    A. Srivastava, G. L. Devnani, and P. Gupta, Materials Advances6, 3043 (2025), publisher: Royal Society of Chemistry (RSC)

    work page 2025

  8. [8]

    S. P. Schwaminger, P. Fraga-Garc´ ıa, M. Eigenfeld, T. M. Becker, and S. Berensmeier, Frontiers in Bioengineering and Biotechnology7(2019), 10.3389/fbioe.2019.00233, publisher: Frontiers

    work page doi:10.3389/fbioe.2019.00233 2019

  9. [9]

    F. Chen, K. A. Smith, and T. A. Hatton, AIChE Journal58, 2865 (2012)

    work page 2012

  10. [10]

    Kolm, IEEE Transactions on Magnetics12, 450 (1976), conference Name: IEEE Transac- tions on Magnetics

    H. Kolm, IEEE Transactions on Magnetics12, 450 (1976), conference Name: IEEE Transac- tions on Magnetics. 24

    work page 1976

  11. [11]

    Tesanovic, T

    M. Tesanovic, T. Bardel, R. Karl, and S. Berensmeier, Current Research in Biotechnology 10, 100324 (2025)

    work page 2025

  12. [12]

    M. Z. Bazant, M. S. Kilic, B. D. Storey, and A. Ajdari, Advances in Colloid and Interface Science152, 48 (2009)

    work page 2009

  13. [13]

    Ebeler, F

    M. Ebeler, F. Pilgram, K. Wolz, G. Grim, and M. Franzreb, Biotechnology Journal13, 1700448 (2018), eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1002/biot.201700448

    work page doi:10.1002/biot.201700448 2018

  14. [14]

    Ebeler, O

    M. Ebeler, O. Lind, N. Norrman, R. Palmgren, and M. Franzreb, New Biotechnology42, 48 (2018)

    work page 2018

  15. [15]

    Krolitzki, S

    E. Krolitzki, S. Steck, A. Nazifi, M. Abt, S. P. Schwaminger, and S. Berensmeier, Separation and Purification Technology323, 124429 (2023)

    work page 2023

  16. [16]

    N. A. Brechmann, P.-O. Eriksson, K. Eriksson, S. Oscarsson, J. Buijs, A. Shokri, G. Hj¨ alm, and V. Chotteau, Biotechnology Progress35, e2775 (2019), eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1002/btpr.2775

    work page doi:10.1002/btpr.2775 2019

  17. [17]

    Tesanovic, I

    M. Tesanovic, I. Zimmermann, M. Schneider, and S. Berensmeier, Separation and Purification Technology369, 132791 (2025), publisher: Elsevier BV

    work page 2025

  18. [18]

    Turco, A

    F. Turco, A. Wegelius, O. Lind, N. Norrman, A.-C. Magnusson, C. Sund-Lundstr¨ om, B. Nor´ en, J. Hedberg, and R. Palmgren, Molecular Therapy - Methods & Clinical Development30, 394 (2023), publisher: Elsevier

    work page 2023

  19. [19]

    Z. Hu, J. Liu, L. Han, Y. Wang, D. Lu, X. Zheng, and Z. Xue, Minerals Engineering183, 107609 (2022), publisher: Elsevier BV

    work page 2022

  20. [20]

    Choomphon-anomakhun, A

    N. Choomphon-anomakhun, A. D. Ebner, M. Natenapit, and J. A. Ritter, Journal of Mag- netism and Magnetic Materials428, 493 (2017)

    work page 2017

  21. [21]

    Tesanovic, J

    M. Tesanovic, J. P. de Souza, M. Z. Bazant, and S. Berensmeier, AIChE Journal71, e18733 (2025), eprint: https://aiche.onlinelibrary.wiley.com/doi/pdf/10.1002/aic.18733

    work page doi:10.1002/aic.18733 2025

  22. [22]

    Wiener, Sachsischen Gesellshaft der Wissenschaften32, 509 (1912)

    O. Wiener, Sachsischen Gesellshaft der Wissenschaften32, 509 (1912)

    work page 1912

  23. [23]

    Torquato,Random Heterogeneous Materials, edited by S

    S. Torquato,Random Heterogeneous Materials, edited by S. S. Antman, L. Sirovich, J. E. Marsden, and S. Wiggins, Interdisciplinary Applied Mathematics, Vol. 16 (Springer, New York, NY, 2002)

    work page 2002

  24. [24]

    E. P. Furlani and K. C. Ng, Physical Review E73, 061919 (2006), publisher: American Physical Society

    work page 2006

  25. [25]

    Hashin and S

    Z. Hashin and S. Shtrikman, Journal of Applied Physics33, 3125 (1962). 25

    work page 1962

  26. [26]

    Hashin and S

    Z. Hashin and S. Shtrikman, Journal of the Mechanics and Physics of Solids10, 343 (1962)

    work page 1962

  27. [27]

    Fraggedakis, M

    D. Fraggedakis, M. Mirzadeh, T. Zhou, and M. Z. Bazant, Journal of The Electrochemical Society167, 113504 (2020), publisher: IOP Publishing

    work page 2020

  28. [28]

    M. Z. Bazant, B. D. Storey, and A. A. Kornyshev, Physical Review Letters106, 046102 (2011), publisher: American Physical Society

    work page 2011

  29. [29]

    R. E. Garc´ ıa, C. M. Bishop, and W. C. Carter, Acta Materialia52, 11 (2004)

    work page 2004

  30. [30]

    N. F. Carnahan and K. E. Starling, The Journal of Chemical Physics51, 635 (1969)

    work page 1969

  31. [31]

    S. R. d. Groot and P. Mazur,Non-Equilibrium Thermodynamics(Amsterdam, North-Holland Pub. Co.; New York, Interscience Publishers, 1962)

    work page 1962

  32. [32]

    G. K. Batchelor, Journal of Fluid Mechanics74, 1 (1976)

    work page 1976

  33. [33]

    M. R. Maxey and J. J. Riley, The Physics of Fluids26, 883 (1983)

    work page 1983

  34. [34]

    Bagchi and S

    P. Bagchi and S. Balachandar, Physics of Fluids15, 3496 (2003)

    work page 2003

  35. [35]

    M. H. Kasbaoui, D. L. Koch, and O. Desjardins, Journal of Fluid Mechanics859, 174 (2019)

    work page 2019

  36. [36]

    J. F. Brady, The Journal of Chemical Physics99, 567 (1993)

    work page 1993

  37. [37]

    Mewis and N

    J. Mewis and N. J. Wagner,Colloidal suspension rheology, 1st ed., Cambridge series in chemical engineering (Cambridge University Press, Cambridge, 2013)

    work page 2013

  38. [38]

    M. Z. Bazant, Accounts of Chemical Research46, 1144 (2013)

    work page 2013

  39. [39]

    T. Gao, Y. Han, D. Fraggedakis, S. Das, T. Zhou, C.-N. Yeh, S. Xu, W. C. Chueh, J. Li, and M. Z. Bazant, Joule5, 393 (2021), publisher: Elsevier

    work page 2021

  40. [40]

    Magnettechnologie in der Verfahrenstechnik w¨ assriger Medien,

    M. Franzreb, “Magnettechnologie in der Verfahrenstechnik w¨ assriger Medien,” (2003). 26

    work page 2003