pith. sign in

arxiv: 2411.16641 · v1 · submitted 2024-11-25 · 🧮 math.NA · cs.NA

On a Completely Discrete Discontinuous Galerkin Method for Incompressible Chemotaxis-Navier-Stokes Equations

Bikram Bir , Harsha Hutridurga , Amiya K. Pani This is my paper

Pith reviewed 2026-05-23 16:33 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords chemotaxisNavier-Stokesdiscontinuous Galerkinerror estimatesfinite elementKeller-Segelincompressible flowfully discrete scheme
0
0 comments X

The pith

A new projection operator enables optimal error estimates for a fully discrete discontinuous Galerkin scheme on the incompressible chemotaxis-Navier-Stokes system.

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

The paper constructs a completely discrete numerical method that pairs a discontinuous Galerkin finite element discretization in space with a semi-implicit first-order finite difference scheme in time for the coupled cell density, chemical concentration, and incompressible fluid equations. It introduces a specially chosen projection operator to close the error analysis and obtain optimal convergence rates in the L2 and H1 norms for the cell density, chemical signal, and fluid velocity, plus an optimal L2 bound for the pressure. These rates are then verified through numerical experiments. A reader would care because reliable, high-order accurate simulations are needed to study how chemical signaling drives collective cell motion inside flowing fluids without introducing uncontrolled numerical artifacts.

Core claim

The paper develops a fully discrete discontinuous Galerkin finite element scheme in space combined with a semi-implicit first-order finite difference method in time for the incompressible Chemotaxis(Keller-Segel)-Navier-Stokes system. With the help of a new projection, it derives optimal error estimates in L2 and H1-norms for the cell density, the concentration of chemical substances and the fluid velocity, together with an optimal error bound in L2-norm for the fluid pressure. Numerical simulations confirm the theoretical rates.

What carries the argument

A new projection operator defined on the discontinuous Galerkin finite element space that satisfies the approximation properties needed to close the error analysis for the coupled nonlinear system.

If this is right

  • The fully discrete solutions converge to the continuous solution at the optimal rates in the listed norms as the mesh size and time step approach zero.
  • The scheme remains stable and accurate for the nonlinear coupling between cell motion, chemical diffusion, and incompressible fluid flow.
  • An optimal L2 error bound holds for the computed fluid pressure without additional post-processing.
  • The method can be implemented directly for practical simulations whose results match the proven theoretical orders.

Where Pith is reading between the lines

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

  • The projection technique may simplify error analysis for other discontinuous Galerkin discretizations of nonlinear fluid-biology systems with similar structure.
  • Long-time simulations using this scheme could reveal whether the discrete solutions preserve qualitative features such as aggregation patterns observed in the continuous model.
  • The approach could be tested on three-dimensional domains or with variable fluid viscosity to check robustness beyond the two-dimensional analysis.

Load-bearing premise

The new projection operator is well-defined on the discontinuous Galerkin space and satisfies the approximation properties needed to close the error analysis for the coupled nonlinear system.

What would settle it

Numerical computations of observed convergence orders on successively refined meshes and time steps that systematically fail to attain the predicted optimal rates in the stated norms would falsify the error analysis.

Figures

Figures reproduced from arXiv: 2411.16641 by Amiya K. Pani, Bikram Bir, Harsha Hutridurga.

Figure 1
Figure 1. Figure 1: Numerical errors with (P1, P0, P1, P1) pairs for Example 5.1. 0.0313 0.0625 0.125 0.25 0.5 h 10-5 10-4 10-3 10-2 10-1 100 Error e u e e c h 3 0.0313 0.0625 0.125 0.25 0.5 h 10-3 10-2 10-1 100 101 Error e u e e c e p h 2 [PITH_FULL_IMAGE:figures/full_fig_p024_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Numerical errors with (P2, P1, P2, P2) pairs for Example 5.1. 0.0313 0.0625 0.125 0.25 0.5 h 10-7 10-6 10-5 10-4 10-3 10-2 10-1 100 Error e u e e c h 4 0.0313 0.0625 0.125 0.25 0.5 h 10-5 10-4 10-3 10-2 10-1 100 101 Error e u e e c e p h 3 [PITH_FULL_IMAGE:figures/full_fig_p024_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Numerical errors with (P3, P2, P3, P3) pairs for Example 5.1. u2(x, y, z, t) = 10 cos(t)w(y) (wz(z)w(x) − w(z)wx(x)), u3(x, y, z, t) = 10 cos(t)w(x) (w(x)wy(y) − wx(x)wy(y)), ρ(x, y, z, t) = 10 cos(t)w(x)w(y)w(z), c(x, y, z, t) = 10 cos(t)wx(x)wy(y)wz(z), p(x, y, z, t) = 10 cos(t)(2x − 1)(2y − 1)(2z − 1), where w(x) = x 2 (1 − x) 2 and ws(s) is the partial derivatives of w(s) with respect to s. 24 [PITH_F… view at source ↗
Figure 4
Figure 4. Figure 4: Numerical errors with respect to time variable for Example [PITH_FULL_IMAGE:figures/full_fig_p025_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Numerical errors with (P1, P0, P1, P1) pairs for Example 5.2. 0.1 0.125 0.1667 0.25 0.5 h 10-7 10-6 10-5 10-4 10-3 10-2 10-1 100 Error e u e e c h 3 0.1 0.125 0.1667 0.25 0.5 h 10-5 10-4 10-3 10-2 10-1 100 101 Error e u e e c e p h 2 [PITH_FULL_IMAGE:figures/full_fig_p025_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Numerical errors with (P2, P1, P2, P2) pairs for Example 5.2. 25 [PITH_FULL_IMAGE:figures/full_fig_p025_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Numerical errors with (P3, P2, P3, P3) pairs for Example 5.2. 5.2 Movement of cells guided by the concentration We take an example where the numerical experiment is carried out for the system (3.31)-(3.34) (that is, for no source terms). This example is taken from [12]. Example 5.3. For this simulation, we consider a rectangular domain Ω = [0, 2] × [0, 1] and the initial data as: ρ0(x, y) = X 3 i=1 70 exp … view at source ↗
Figure 8
Figure 8. Figure 8: The evolution of the cell density (left) and the concentra [PITH_FULL_IMAGE:figures/full_fig_p027_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: The velocity vector (left) and the pressure contour (rig [PITH_FULL_IMAGE:figures/full_fig_p028_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: The evolution of the velocity components (left: primary c [PITH_FULL_IMAGE:figures/full_fig_p029_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: Deviation of mass of discrete cell density for linear(left) [PITH_FULL_IMAGE:figures/full_fig_p030_11.png] view at source ↗
read the original abstract

This paper deals with a fully discrete numerical scheme for the incompressible Chemotaxis(Keller-Segel)-Navier-Stokes system. Based on a discontinuous Galerkin finite element scheme in the spatial directions, a semi-implicit first-order finite difference method in the temporal direction is applied to derive a completely discrete scheme. With the help of a new projection, optimal error estimates in $L^2$ and $H^1$-norms for the cell density, the concentration of chemical substances and the fluid velocity are derived. Further, optimal error bound in $L^2$-norm for the fluid pressure is obtained. Finally, some numerical simulations are performed, whose results confirm the theoretical findings.

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 develops a fully discrete discontinuous Galerkin scheme for the incompressible chemotaxis-Navier-Stokes equations, combining DG spatial discretization with a semi-implicit first-order finite-difference time discretization. A new projection operator is introduced to derive optimal L² and H¹ error estimates for the cell density, chemical concentration, and fluid velocity, together with an optimal L² bound for the pressure; the analysis is supported by numerical experiments that confirm the predicted convergence rates.

Significance. If the new projection is shown to be well-defined on the DG space and to deliver the stated approximation properties for the coupled nonlinear system, the work supplies a rigorous optimal-rate error analysis for a challenging class of chemotaxis-fluid models. The combination of a completely discrete scheme, explicit projection-based estimates, and confirming numerical tests constitutes a solid contribution to numerical analysis of coupled PDE systems.

major comments (2)
  1. [Section 3] Section 3 (definition of the new projection): the manuscript must explicitly verify that the projection is well-defined on the chosen DG finite-element space and that its approximation properties hold uniformly for the nonlinear coupling terms, as these properties are load-bearing for closing the error estimates in L² and H¹.
  2. [Section 4] Section 4 (error analysis): the treatment of the convective and chemotactic nonlinearities in the error equations relies on the projection; any hidden dependence on mesh-dependent constants or inverse inequalities should be tracked explicitly to confirm the optimality of the rates.
minor comments (2)
  1. The notation for the projection operator and its properties could be collected in a single preliminary subsection for easier reference during the error analysis.
  2. A brief remark on the choice of the semi-implicit time discretization (why first-order rather than higher-order) would help readers assess the trade-off between stability and accuracy.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading of the manuscript and the constructive comments. We address each major comment below and will revise the manuscript to incorporate the requested clarifications.

read point-by-point responses

  1. Referee: [Section 3] Section 3 (definition of the new projection): the manuscript must explicitly verify that the projection is well-defined on the chosen DG finite-element space and that its approximation properties hold uniformly for the nonlinear coupling terms, as these properties are load-bearing for closing the error estimates in L² and H¹.

    Authors: We agree that an explicit verification strengthens the presentation. In the revised manuscript we will add a dedicated paragraph (or short subsection) in Section 3 that proves the projection is well-defined on the DG space by verifying the associated linear system is uniquely solvable, and that the approximation properties hold uniformly with respect to the nonlinear coupling terms appearing in the chemotaxis-Navier-Stokes system. revision: yes

  2. Referee: [Section 4] Section 4 (error analysis): the treatment of the convective and chemotactic nonlinearities in the error equations relies on the projection; any hidden dependence on mesh-dependent constants or inverse inequalities should be tracked explicitly to confirm the optimality of the rates.

    Authors: We accept the suggestion. In the revised Section 4 we will explicitly collect and bound every constant that arises when estimating the convective and chemotactic terms, including those originating from inverse inequalities. We will show that these constants remain independent of the mesh size in a manner that preserves the stated optimal L² and H¹ rates. revision: yes

Circularity Check

0 steps flagged

No significant circularity; error analysis self-contained

full rationale

The derivation introduces a new projection operator on the DG space, establishes its well-definedness and approximation properties independently, then applies these to obtain optimal L2/H1 error bounds for the fully discrete scheme on the coupled system. This is a standard, non-circular technique in numerical analysis of nonlinear PDEs; the estimates follow from the discrete scheme, projection properties, and standard inequalities rather than reducing to fitted inputs, self-definitions, or self-citation chains. No load-bearing steps match the enumerated circularity patterns.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on standard finite-element approximation theory plus the existence and properties of an unspecified new projection operator tailored to the discontinuous Galerkin space for the nonlinear chemotaxis-Navier-Stokes coupling.

axioms (1)
  • domain assumption The continuous problem admits a sufficiently regular solution so that the projection error terms can be bounded at the optimal rate.
    Required to convert consistency and stability estimates into the stated optimal convergence rates.

pith-pipeline@v0.9.0 · 5653 in / 1229 out tokens · 24344 ms · 2026-05-23T16:33:43.839955+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

39 extracted references · 39 canonical work pages

  1. [1]

    Adams and J

    R. Adams and J. F. Fournier. Sobolev spaces. Acad. Press, New York , 19(5), 1975

    work page 1975

  2. [2]

    D. N. Arnold. An interior penalty finite element method with discont inuous elements. SIAM J. Numer. Anal., 19(4):742–760, 1982

    work page 1982

  3. [3]

    Bajpai, D

    S. Bajpai, D. Goswami, and K. Ray. A priori error estimates of a d iscontinuous Galerkin method for the Navier-Stokes equations. Numer. Algorithms , 94(2):937–1002, 2023

    work page 2023

  4. [4]

    Cao and J

    X. Cao and J. Lankeit. Global classical small-data solutions for a t hree-dimensional chemotaxis Navier- Stokes system involving matrix-valued sensitivities. Calc. Var. Partial Differential Equations , 55(4):Art. 107, 39, 2016

    work page 2016

  5. [5]

    Caucao, E

    S. Caucao, E. Colmenares, G. N. Gatica, and C. Inzunza. A Bana ch spaces-based fully-mixed finite element method for the stationary chemotaxis-Navier-Stokes pr oblem. Comput. Math. Appl. , 145:65–89, 2023

    work page 2023

  6. [6]

    Chertock, K

    A. Chertock, K. Fellner, A. Kurganov, A. Lorz, and P. A. Marko wich. Sinking, merging and stationary plumes in a coupled chemotaxis-fluid model: a high-resolution numerica l approach. J. Fluid Mech. , 694:155–190, 2012

    work page 2012

  7. [7]

    P. G. Ciarlet. The Finite Element Method for Elliptic Problem . North-Holland Publ. Comp., Amsterdam, 1978

    work page 1978

  8. [8]

    Deleuze, C.-Y

    Y. Deleuze, C.-Y. Chiang, M. Thiriet, and T. W. H. Sheu. Numerical study of plume patterns in a chemotaxis-diffusion-convection coupling system. Comput. & Fluids , 126:58–70, 2016

    work page 2016

  9. [9]

    D. A. Di Pietro and A. Ern. Mathematical aspects of discontinuous Galerkin methods , volume 69. Springer Science & Business Media, 2011

    work page 2011

  10. [10]

    Douglas and T

    J. Douglas and T. Dupont. Interior penalty procedures for ellip tic and parabolic Galerkin methods. In Computing methods in applied sciences (Second Internat. Sy mpos., Versailles, 1975) , Lecture Notes in Phys., Vol. 58, pages 207–216. Springer, Berlin-New York, 1976

    work page 1975

  11. [11]

    R. Duan, A. Lorz, and P. Markowich. Global solutions to the cou pled chemotaxis-fluid equations. Comm. Partial Differential Equations , 35(9):1635–1673, 2010. 30

    work page 2010

  12. [12]

    Duarte-Rodr ´ ıguez, M

    A. Duarte-Rodr ´ ıguez, M. ´A. Rodr ´ ıguez-Bellido, D. A. Rueda-G´ omez, and ´E. J. Villamizar-Roa. Numerical analysis for a chemotaxis–Navier-Stokes system. ESAIM Math. Model. Numer. Anal. , 55(suppl.):S417–S445, 2021

    work page 2021

  13. [13]

    Epshteyn

    Y. Epshteyn. Discontinuous Galerkin methods for the chemota xis and haptotaxis models. J. Comput. Appl. Math. , 224(1):168–181, 2009

    work page 2009

  14. [14]

    Epshteyn and A

    Y. Epshteyn and A. Izmirlioglu. Fully discrete analysis of a discont inuous finite element method for the Keller-Segel chemotaxis model. J. Sci. Comput. , 40(1-3):211–256, 2009

    work page 2009

  15. [15]

    Epshteyn and A

    Y. Epshteyn and A. Kurganov. New interior penalty discontinuo us Galerkin methods for the Keller-Segel chemotaxis model. SIAM J. Numer. Anal. , 47(1):386–408, 2008/09

    work page 2008

  16. [16]

    Girault and P

    V. Girault and P. A. Raviart. Finite element approximation of the Navier-Stokes equatio ns, volume 749. Springer Berlin, 1979

    work page 1979

  17. [17]

    Girault, B

    V. Girault, B. Rivi` ere, and M. F. Wheeler. A discontinuous Galer kin method with nonoverlapping domain decomposition for the Stokes and Navier-Stokes problems. Math. Comp. , 74(249):53–84, 2005

    work page 2005

  18. [18]

    Girault, B

    V. Girault, B. Rivi` ere, and M. F. Wheeler. A splitting method usin g discontinuous Galerkin for the transient incompressible Navier-Stokes equations. M2AN Math. Model. Numer. Anal. , 39(6):1115–1147, 2005

    work page 2005

  19. [19]

    F. Hecht. New development in freefem++. J. Numer. Math. , 20(3-4):251–265, 2012

    work page 2012

  20. [20]

    A. T. Hill and E. S¨ uli. Approximation of the global attractor for the incompressible Navier-Stokes equations. IMA J. Numer. Anal. , 20(4):633–667, 2000

    work page 2000

  21. [21]

    D. Kay, V. Styles, and E. S¨ uli. Discontinuous Galerkin finite element approximation of the Cahn-Hilliard equation with convection. SIAM J. Numer. Anal. , 47(4):2660–2685, 2009

    work page 2009

  22. [22]

    Kaya and B

    S. Kaya and B. Rivi` ere. A discontinuous subgrid eddy viscosity method for the time-dependent Navier- Stokes equations. SIAM J. Numer. Anal. , 43(4):1572–1595, 2005

    work page 2005

  23. [23]

    H. G. Lee and J. Kim. Numerical investigation of falling bacterial p lumes caused by bioconvection in a three-dimensional chamber. Eur. J. Mech. B Fluids , 52:120–130, 2015

    work page 2015

  24. [24]

    X. H. Li, C.-W. Shu, and Y. Yang. Local discontinuous Galerkin me thod for the Keller-Segel chemotaxis model. J. Sci. Comput. , 73(2-3):943–967, 2017

    work page 2017

  25. [25]

    Z. Li, L. Xiao, M. Li, and H. Chen. Error estimates for the finite e lement method of the chemotaxis- Navier-Stokes equations. J. Appl. Math. Comput. , 69(4):3039–3065, 2023

    work page 2023

  26. [26]

    J. G. Liu and A. Lorz. A coupled chemotaxis-fluid model: global ex istence. Ann. Inst. H. Poincar´ e C Anal. Non Lin´ eaire, 28(5):643–652, 2011

    work page 2011

  27. [27]

    A. Lorz. Coupled chemotaxis fluid model. Math. Models Methods Appl. Sci. , 20(6):987–1004, 2010

    work page 2010

  28. [28]

    Masri, C

    R. Masri, C. Liu, and B. Riviere. Improved a priori error estimat es for a discontinuous Galerkin pressure correction scheme for the Navier-Stokes equations. Numer. Methods Partial Differential Equations , 39(4):3108–3144, 2023

    work page 2023

  29. [29]

    C. Qiu, Q. Liu, and J. Yan. Third order positivity-preserving dire ct discontinuous Galerkin method with interface correction for chemotaxis Keller-Segel equations. J. Comput. Phys. , 433:Paper No. 110191, 17, 2021

    work page 2021

  30. [30]

    W. H. Reed and T. R. Hill. Triangular mesh methods for the neutro n transport equation. Technical report, Los Alamos Scientific Lab., N. Mex.(USA), 1973

    work page 1973

  31. [31]

    Rivi` ere

    B. Rivi` ere. Discontinuous Galerkin methods for solving elliptic and pa rabolic equations: theory and implementation. SIAM, 2008

    work page 2008

  32. [32]

    Y. Tang, G. Zou, and J. Li. Unconditionally energy-stable finite e lement scheme for the chemotaxis-fluid system. J. Sci. Comput. , 95(1):Paper No. 1, 34, 2023. 31

    work page 2023

  33. [33]

    R. Temam. Navier-Stokes Equations: Theory and Numerical Analysis , volume 2. Amsterdam, North- Holland Publishing Co., Third edition, 1984

    work page 1984

  34. [34]

    Tuval, L

    I. Tuval, L. Cisneros, C. Dombrowski, C. W. Wolgemuth, J. O. Ke ssler, and R. E. Goldstein. Bacterial swimming and oxygen transport near contact lines. Proc. Natl. Acad. Sci. , 102(7):2277–2282, 2005

    work page 2005

  35. [35]

    M. Wang, G. Zou, B. Wang, and W. Zhao. Unconditionally energy- stable discontinuous Galerkin method for the chemo-repulsion-Navier-Stokes system. Comput. Math. Appl. , 150:132–155, 2023

    work page 2023

  36. [36]

    W. Wang. Global boundedness in a 2D chemotaxis-Navier-Stoke s system with flux limitation and nonlinear production. Math. Models Methods Appl. Sci. , 33(11):2361–2394, 2023

    work page 2023

  37. [37]

    M. F. Wheeler. An elliptic collocation-finite element method with inte rior penalties. SIAM J. Numer. Anal., 15(1):152–161, 1978

    work page 1978

  38. [38]

    M. Winkler. Global large-data solutions in a chemotaxis-(Navier- )Stokes system modeling cellular swimming in fluid drops. Comm. Partial Differential Equations , 37(2):319–351, 2012

    work page 2012

  39. [39]

    Zhang and X

    Q. Zhang and X. Zheng. Global well-posedness for the two-dime nsional incompressible chemotaxis- Navier-Stokes equations. SIAM J. Math. Anal. , 46(4):3078–3105, 2014. 32

    work page 2014