Electrokinetic Flow Modeling Through Bone Scaffold

Debabrata DasGupta; Mohammad Osaid

REVIEW 4 major objections 5 minor 11 references

Electrokinetic Flow Modeling Through Bone Scaffold

T0 review · 4 major / 5 minor · reviewed 2026-08-10 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read Electrokinetic flow through a bone scaffold produces higher velocity and higher wall shear stress than pressure-driven flow at unit driving gradients, the paper argues.

desk verdict A legitimate modeling extension undone by an apples-to-oranges comparison and an abstract-vs-discussion contradiction; worth sending to reviewers, but not believable as-is. read the letter →

arxiv 2501.04025 v1 pith:KN4NKSR7 submitted 2024-12-30 physics.bio-ph

Mohammad Osaid , Debabrata Dasgupta This is my paper · ORCID

classification physics.bio-ph

keywords electrokineticflowbonescaffoldbioreactorelectro-osmosiswallshearstressoxygentransporthomogenizationperfusion

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

This paper models electrokinetic flow—flow driven by an electric field acting on the ionic double layer at the bone scaffold surface—through a trabecular bone scaffold in a perfusion bioreactor, and compares it with conventional pressure-driven flow. The authors aim to show that, at a unit electric field versus a unit pressure gradient, the electrokinetic mechanism produces higher fluid velocity, higher wall shear stress on the bone, and comparable oxygen concentration, and therefore could be a better transport mode for growing bone grafts. The study matters because nutrient and oxygen delivery inside scaffolds limits how large and viable tissue-engineered bone grafts can be, and the choice of how to pump medium is a design decision that affects both transport and mechanical stimulation of cells.

What carries the argument

The central mechanism is electro-osmotic flow: polar interactions between the bone surface and the nutrient medium create an electric double layer with a zeta potential, and an externally applied electric field accelerates the ions in that layer, dragging the surrounding fluid. The modeling pipeline solves the Poisson–Boltzmann equation for the double-layer potential and the Navier–Stokes equations with an electrical body-force term for the flow, on a reference volume element (RVE) of micro-CT-scanned trabecular bone, then upscales the result by homogenization. The argument turns on the difference in velocity profile shape: pressure-driven flow gives a parabolic profile, while electrokinetic flow gives an almost uniform plug profile, producing a higher velocity gradient and hence higher wall shear stress.

What would settle it

Measure the flow rate through the same trabecular bone scaffold under a controlled pressure drop and under a controlled applied voltage, and compute the delivered volumetric flux per watt of input power (hydraulic versus electrical). If pressure-driven flow delivers equal or higher flux per watt, or if the wall shear stress per unit power is not higher for electrokinetic flow, the paper's central claim of advantage would be undermined.

Watch

Extended reading notes

Core claim

On the paper's own terms, the central discovery is that electrokinetic transport through a bone scaffold outperforms pressure-driven perfusion when both are driven by unit gradients: for a unit electric field the velocity is considerably higher than for a unit pressure gradient, wall shear stress is higher because the velocity profile is nearly uniform (plug-like) rather than parabolic, and the oxygen concentration inside the scaffold is nearly the same in both cases. The higher velocity increases the flux of nutrient medium through the bone, which the authors argue delivers more glucose to cells and increases cellular activity, while the higher shear stress provides stronger mechanical stimulation. The conclusion states that an electrokinetics-based bioreactor therefore has several advantages over a pressure-based one for bone growth.

Load-bearing premise

The load-bearing premise is that a 1 V/m electric field and a 1 Pa/m pressure gradient are equally strong drivers, so the higher velocity and shear reported for electrokinetic flow may be an artifact of comparing two quantities that cannot be equated without an energy or cost normalization.

Editorial extensions

If this is right

Electrokinetic perfusion bioreactors could deliver nutrient medium through bone scaffolds at higher flux than pressure-driven ones operating at the same nominal unit gradient.
The higher wall shear stress from electrokinetic flow would give stronger mechanical stimulation to bone cells, which the paper links to increased cell growth and activity.
Oxygen concentration inside the scaffold would remain essentially the same as in pressure-driven flow, so the transport advantage comes without a loss of oxygen delivery.
Electrokinetic bioreactors offer simpler control, including easier actuation and manipulation of flow, compared with pressure-based systems.
The model's prediction of near-uniform velocity profiles suggests that electrokinetic flow may produce more uniform shear stress across the scaffold walls.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

The comparison uses 1 V/m and 1 Pa/m as if they were equivalent driving strengths; normalizing by energy input or pumping power could change or reverse the apparent advantage of electrokinetic flow, and the authors do not supply such a normalization.
If the plug-like velocity profile holds more generally, electrokinetic flow might reduce dead zones and give more spatially uniform nutrient delivery than parabolic pressure-driven flow, a consequence the paper does not explicitly develop.
The same electrokinetic model could be extended to other metabolites, such as glucose, and to waste-product removal, as well as to oscillating or pulsed electric fields, none of which are tested in the current steady-field setup.
A direct experimental test would be to measure the volumetric flow rate through a bone scaffold under a known pressure drop and a known applied voltage and compare performance per unit of electrical work versus hydraulic power.

Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Desk Editor's Note

read the letter

Colleague,

The useful question here is whether electrokinetic perfusion can replace pressure-driven flow in bone scaffolds. The paper builds a microCT-based geometry, borrows the homogenization framework from Bandopadhyay et al., and runs simulations for velocity, wall shear, and oxygen transport. That is a legitimate extension, and the microCT-derived RVE is real data, not a cartoon. I give them credit for that.

Where it falls apart is the comparison. The conclusion hinges on setting a unit pressure gradient against a unit electric field. Those are different physical quantities. The resulting velocity ratio depends on permeability, viscosity, zeta potential, permittivity, and ion concentration, none of which are pinned down in the paper. With typical values, 1 Pa/m can actually produce faster flow than 1 V/m, so the claimed advantage in velocity and shear may be an artifact of that arbitrary choice. A fair comparison needs a common operating condition: equal volumetric flow rate, equal wall shear, or equal power input. Without that, the central quantitative claim does not hold.

There is also an internal contradiction. The abstract touts improved oxygen transport, but the Discussion says oxygen concentration is almost the same in both cases. That needs to be reconciled. The paper is also missing mesh convergence, parameter values (zeta potential, ion concentration), and any code or data. Minor but annoying: two equations are numbered (8), which suggests the manuscript was assembled in a hurry.

The application itself is not silly. Electrokinetic flow in a trabecular bone scaffold could be worth modeling, and the homogenization approach is sensible. But as presented, the paper does not establish that electrokinetic transport is better. It would need a meaningful normalization, a parameter sweep, and resolution of the oxygen inconsistency.

I would send it to peer review rather than desk-reject, because the underlying question is legitimate and the scaffold geometry is real. A referee could push the authors to fix the comparison and add the missing details. But if I were the reviewer, my verdict would be major revision or reject, not accept. For a reading group, it is a good case study in how unit choices can drive a conclusion.

Best

Referee Report

4 major / 5 minor

Summary. The manuscript models electrokinetic flow through a bone scaffold RVE derived from micro-CT scans and compares it to pressure-driven flow. The stated central claims are that electrokinetic transport provides improved oxygen transport, higher fluid velocity, and higher wall shear stress than pressure-driven flow, and therefore is a preferable bioreactor mode for bone growth. The paper presents a homogenization-based theoretical framework, 2D and 3D simulations, and comparisons of velocity, shear stress, and oxygen concentration. The central quantitative conclusion rests on comparing a unit pressure gradient with a unit electric field.

Significance. If the comparison were physically meaningful, the paper could inform bioreactor design for bone tissue engineering by identifying electrokinetic driving as a way to control nutrient delivery and mechanical stimulation. The manuscript also attempts to combine homogenization theory with realistic micro-CT geometry, and it names concrete transport and shear-stress metrics. However, the central claim is undermined by the arbitrary normalization of the two driving forces and by an internal contradiction between the abstract and the Discussion. Because these issues affect the main conclusion rather than being local presentation problems, the current version does not provide reliable evidence for the claimed advantage.

major comments (4)

[Discussion, paragraph 3] The comparison between electrokinetic and pressure-driven transport is anchored on treating 'unit pressure gradient' and 'unit electric field' as equivalent driving strengths. These quantities have different units (Pa/m is a force per unit volume, V/m is a force per unit charge), so setting both to one is an arbitrary normalization. The resulting ratio of velocities depends on hidden material parameters such as permeability, viscosity, zeta potential, and permittivity, and no common operating constraint (equal volumetric flow rate, equal power input, equal maximum/average velocity, or equal wall shear) is imposed. This makes the reported higher velocity and higher wall shear for electrokinetic flow (Figures 6 and 7, Conclusion) artifacts of the chosen unit magnitudes rather than evidence of superiority.
[Abstract and Discussion, paragraph 2] The abstract claims 'improved oxygen transport' as an advantage of electrokinetic transport, but the Discussion states that 'the concentration of oxygen is almost the same in both the transport processes in spite of having completely different physics behind the transport.' Because oxygen concentration is the direct transport metric used in this study, the abstract's claim is contradicted by the manuscript's own results. This internal inconsistency prevents the reader from assessing what advantage, if any, electrokinetic flow actually provides for nutrient delivery.
[Governing Equations, Eqs. (1)-(10)] The theoretical framework is not sufficiently specified to support the simulation results. Eq. (8) appears twice: once for the diffusion equation and once for the dimensionless Poisson-Boltzmann equation, with the latter using symbols (A², ζ̄, ψ̄₀) that are not defined in the text. More importantly, no numerical values are given for the zeta potential ζ, bulk ion concentration n₀, or applied electric field magnitude used in the 2D and 3D simulations. Without these parameters, the electrokinetic velocity cannot be reproduced, and the claimed velocity advantage cannot be independently checked. The paper should state the parameter values and justify the chosen unit-field magnitude against any physical or experimental reference.
[Setup and Results] The manuscript provides no mesh convergence study and no validation of the homogenization assumption for the specific bone RVE. The RVE is a single 2×2×2 mm³ cube, and the paper asserts that it represents the larger bone sample, but no evidence is given that the RVE size is sufficient for scale separation or that the results are independent of the chosen sample location. These omissions are load-bearing because the central quantitative comparisons (velocity and shear) are derived from this geometry and upscaled to the whole scaffold.

minor comments (5)

[Governing Equations, Eq. (8)] The duplicate equation numbering should be corrected; the dimensionless Poisson-Boltzmann equation should have a distinct label and fully defined variables.
[Fig. 3, Fig. 5, Fig. 6, Fig. 7] The figures show color contours without visible axis labels, units, or color-bar scales. This makes it difficult to interpret the reported velocity and concentration values.
[Introduction, references] Reference [2] appears to be about numerical simulation of drops, unrelated to the sentence citing static vs. dynamic culture systems; please verify the citation.
[Results, 'Electrokinetic Transport in a 2D Domain'] The sentence 'The figure 5 illustrates that the velocity field gradient in this case is higher compared to pressure-driven flow' is not supported by a quantitative comparison in the text; the figure alone is insufficient.
[Discussion, paragraph 3] The phrase 'an interesting thing in this discussion is the velocity profile' is informal; consider a more precise statement about the shape of the velocity profile and its effect on wall shear.

Circularity Check

0 steps flagged · score 0.0 of 10

No circular derivation chain; the central comparison is a normalization/fairness issue, not an input-output equivalence.

full rationale

The paper's derivation chain is not circular. The governing equations (Poisson-Boltzmann, Navier-Stokes with an electrical body force, and convection-diffusion with Michaelis-Menten consumption) are standard physical models stated explicitly, and the reported velocity, shear stress, and oxygen concentration fields are numerical outputs of those equations under stated boundary conditions. Homogenization is imported from Bandopadhyay et al. (reference [11]), which is not a self-citation for the present authors and is an externally published method, so no load-bearing result reduces to a self-citation. No parameter is fitted to the reported outputs; zeta potential, permeability, diffusion coefficients, and consumption kinetics are all pre-specified inputs. The conclusion that electrokinetic flow gives higher velocity and wall shear than pressure-driven flow is a contingent result of those inputs, not an identity: the ratio of electroosmotic to pressure-driven velocity depends on material parameters (epsilon, zeta, permeability, viscosity) and could invert under different choices. The genuine weakness is that the comparison equates 1 V/m with 1 Pa/m, which are not physically commensurable driving strengths, and no common power, flow-rate, or shear constraint is imposed; this is an external-validity or fairness flaw, not circularity, because the output is not equivalent to the input by construction. The abstract's claim of improved oxygen transport is also inconsistent with the Discussion's statement that oxygen concentration is almost the same in both processes, but this is an internal inconsistency rather than a circular derivation.

Assumptions & free parameters 8 free parameters · 6 assumptions · 0 invented entities

The central comparison depends on a large set of assumed parameters and boundary conditions. The most consequential is the arbitrary unit-gradient normalization, which makes the velocity and shear advantages non-unique, followed by unspecified electrokinetic parameters such as zeta potential and bulk ion concentration.

free parameters (8)

Scaffold permeability K = 4.5 x 10^-10 (units not stated)
Assumed in the porous-media pressure-driven model; not measured on the actual bone sample.
Inlet velocity = 0.05 m/s
Chosen for the pressure-driven simplified geometry without justification.
Zeta potential zeta = unspecified constant
Central to electrokinetic velocity, but its value is never reported.
Bulk ion concentration n0 = unspecified
Needed for Debye length in Eq (10); not given.
Permittivity epsilon_r = unspecified
Appears in Poisson and Debye length; not given.
Applied electric field = 'unit' field, no units stated
Used as the electrokinetic driving force; without units the velocity comparison is ambiguous.
Oxygen diffusion coefficients = 3.29e-9 m^2/s (medium), 2.0e-9 m^2/s (cell layer)
Taken from literature, not measured for this scaffold; affect concentration results.
Inlet/outlet oxygen concentrations = 0.21 and 0.1 (or 2.1e-1 mol/m^3)
Boundary values chosen for the electrokinetic model; no sensitivity analysis.

assumptions (6)

domain assumption Incompressible Newtonian flow in the liquid domain with Navier-Stokes and electrical body force
Standard model for cell culture media; not experimentally verified here.
domain assumption Poisson-Boltzmann equilibrium for the electric double layer with constant zeta potential at walls
Invoked in Eqs (1)-(4); assumes equilibrium ion distribution near charged bone surfaces.
domain assumption Separation of scales epsilon = l/L << 1 and periodic microstructure for homogenization
Invoked in 'Theoretical Framework and Setup'; justifies upscaling RVE results, but no convergence test is shown.
domain assumption Michaelis-Menten kinetics describe oxygen consumption by bone cells
Assumed for the negative source term; parameters are from literature rather than measured for this system.
ad hoc to paper The single microCT-derived RVE is representative of the larger bone sample
No representativeness or statistical analysis is provided; one 2 mm cube is used.
domain assumption No-slip and no-flux boundary conditions at scaffold walls
Standard boundary conditions, but their applicability to cell-covered walls is not discussed.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Electrokinetic Flow Modeling Through Bone Scaffold." pith.science (2026). https://pith.science/paper/KN4NKSR7

@misc{pith2026250104025,
  author       = {Pith},
  title        = {Pith review of: Electrokinetic Flow Modeling Through Bone Scaffold},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/KN4NKSR7}},
  note         = {Machine review of arXiv:2501.04025}
}

read the original abstract

Fluid flow and mass transfer inside a bioreactor play a pivotal role in growing bone grafts, as cell proliferation is limited by the transport of nutrients and oxygen, as well as the removal of by-products from cells within the scaffold. Traditionally, perfusion bioreactors are used for tissue-engineered bone grafts. In this study, we modeled electrokinetic flow through the graft for tissue-engineered bone grafts and compared it to pressure-driven flows. The study highlights that electrokinetic transport offers several advantages over pressure-driven flow, such as improved oxygen transport and higher shear stress. Additionally, it provides generic benefits, including greater flexibility in control mechanisms, such as easier actuation and manipulation of flow.

Figures

Figures reproduced from arXiv: 2501.04025 by the authors.

**Figure 4.** Results of pressure-based flow on a simplified geometry. a) Variation in shear stress on the wall. b) Oxygen concentration in a different parts of the scaffold. • R = -(1/vₐₑₗₗ) × (Qₘc / (Cₘ + c)) Where: • vₐₑₗₗ is the cell volume (1.44 × 10⁻¹⁵ m³/cell), • Qₘ is the maximal oxygen consumption rate (1.86 × 10⁻¹⁸ mol/cell/s), • Cₘ is the oxygen concentration at half-maximal consumption (6.0 × 10⁻³ mol/m³). The dissolv… view at source ↗

**Figure 7.** Comparison of Oxygen (O2) concentraion and velocity profiles in electrokinetic and pressure driven flow. References 1. Martin, I., T. Smith, and D. Wendt, Bioreactor-based roadmap for the translation of tissue engineering strategies into clinical products. Trends in biotechnology, 2009. 27(9): p. 495- 502. 2. Dupont, J.-B. and D. Legendre, Numerical simulation of static and sliding drop with contact angle hysteresis… view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

11 extracted references · 11 canonical work pages

[1]

Smith, and D

Martin, I., T. Smith, and D. Wendt, Bioreactor-based roadmap for the translation of tissue engineering strategies into clinical products. Trends in biotechnology, 2009. 27(9): p. 495- 502

work page 2009
[2]

Dupont, J.-B. and D. Legendre, Numerical simulation of static and sliding drop with contact angle hysteresis. Journal of Computational Physics, 2010. 229(7): p. 2453-2478

work page 2010
[3]

Biotechnology and bioengineering, 2004

Sucosky, P., et al., Fluid mechanics of a spinner‐flask bioreactor. Biotechnology and bioengineering, 2004. 85(1): p. 34-46

work page 2004
[4]

Magrofuoco, and N

Flaibani, M., E. Magrofuoco, and N. Elvassore, Computational modeling of cell growth heterogeneity in a perfused 3D scaffold. Industrial & engineering chemistry research, 2010. 49(2): p. 859-869

work page 2010
[5]

Radisic, and B

Vunjak‐Novakovic, G., M. Radisic, and B. Obradovic, Cardiac tissue engineering: effects of bioreactor flow environment on tissue constructs. Journal of Chemical Technology & Biotechnology: International Research in Process, Environmental & Clean Technology, 2006. 81(4): p. 485-490

work page 2006
[6]

Biomaterials, 2001

Goldstein, A.S., et al., Effect of convection on osteoblastic cell growth and function in biodegradable polymer foam scaffolds. Biomaterials, 2001. 22(11): p. 1279-1288

work page 2001
[7]

Biotechnology and bioengineering, 2011

Grayson, W.L., et al., Optimizing the medium perfusion rate in bone tissue engineering bioreactors. Biotechnology and bioengineering, 2011. 108(5): p. 1159-1170

work page 2011
[8]

Journal of biomechanics, 2006

Boschetti, F., et al., Prediction of the micro-fluid dynamic environment imposed to three- dimensional engineered cell systems in bioreactors. Journal of biomechanics, 2006. 39(3): p. 418-425

work page 2006

Show all 11 references

[9]

Eur Cell Mater, 2016

Coughlin, T., et al., Primary cilia expression in bone marrow in response to mechanical stimulation in explant bioreactor culture. Eur Cell Mater, 2016. 32: p. 111-122

2016
[10]

Journal of Theoretical Biology, 2009

Whittaker, R.J., et al., Mathematical modelling of fibre-enhanced perfusion inside a tissue- engineering bioreactor. Journal of Theoretical Biology, 2009. 256(4): p. 533-546

2009
[11]

Physical Review E—Statistical, Nonlinear, and Soft Matter Physics, 2013

Bandopadhyay, A., et al., Electro-osmotic flows through topographically complicated porous media: Role of electropermeability tensor. Physical Review E—Statistical, Nonlinear, and Soft Matter Physics, 2013. 87(3): p. 033006

2013

Pith tools

Reviewed August 10, 2026 · model on record in the stance chip above.

[1] [1]

Smith, and D

Martin, I., T. Smith, and D. Wendt, Bioreactor-based roadmap for the translation of tissue engineering strategies into clinical products. Trends in biotechnology, 2009. 27(9): p. 495- 502

work page 2009

[2] [2]

Dupont, J.-B. and D. Legendre, Numerical simulation of static and sliding drop with contact angle hysteresis. Journal of Computational Physics, 2010. 229(7): p. 2453-2478

work page 2010

[3] [3]

Biotechnology and bioengineering, 2004

Sucosky, P., et al., Fluid mechanics of a spinner‐flask bioreactor. Biotechnology and bioengineering, 2004. 85(1): p. 34-46

work page 2004

[4] [4]

Magrofuoco, and N

Flaibani, M., E. Magrofuoco, and N. Elvassore, Computational modeling of cell growth heterogeneity in a perfused 3D scaffold. Industrial & engineering chemistry research, 2010. 49(2): p. 859-869

work page 2010

[5] [5]

Radisic, and B

Vunjak‐Novakovic, G., M. Radisic, and B. Obradovic, Cardiac tissue engineering: effects of bioreactor flow environment on tissue constructs. Journal of Chemical Technology & Biotechnology: International Research in Process, Environmental & Clean Technology, 2006. 81(4): p. 485-490

work page 2006

[6] [6]

Biomaterials, 2001

Goldstein, A.S., et al., Effect of convection on osteoblastic cell growth and function in biodegradable polymer foam scaffolds. Biomaterials, 2001. 22(11): p. 1279-1288

work page 2001

[7] [7]

Biotechnology and bioengineering, 2011

Grayson, W.L., et al., Optimizing the medium perfusion rate in bone tissue engineering bioreactors. Biotechnology and bioengineering, 2011. 108(5): p. 1159-1170

work page 2011

[8] [8]

Journal of biomechanics, 2006

Boschetti, F., et al., Prediction of the micro-fluid dynamic environment imposed to three- dimensional engineered cell systems in bioreactors. Journal of biomechanics, 2006. 39(3): p. 418-425

work page 2006