pith. sign in

arxiv: 2507.01402 · v2 · submitted 2025-07-02 · 🧮 math.NA · cs.NA· math-ph· math.MP

Asymptotic Preserving and Accurate scheme for Multiscale Poisson-Nernst-Planck (MPNP) system

Clarissa Astuto , Giovanni Russo This is my paper

Pith reviewed 2026-05-19 06:58 UTC · model grok-4.3

classification 🧮 math.NA cs.NAmath-phmath.MP
keywords multiscale modelingPoisson-Nernst-Planck equationsasymptotic preserving schemestwo-species ion transportsurface trapsboundary condition reductionnumerical analysis for PDEs
0
0 comments X

The pith

A two-species multiscale Poisson-Nernst-Planck model replaces the small-range trap potential with a boundary condition and supplies an asymptotic-preserving numerical scheme that couples positive and negative ions.

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

The paper establishes a multiscale model for the Poisson-Nernst-Planck system that treats positive and negative ions together as they interact through the Coulomb field near a surface trap whose attraction range is much smaller than the overall domain. The detailed attractive potential is replaced by a boundary condition obtained from mass conservation and asymptotic analysis, and an asymptotic-preserving scheme is constructed to solve the resulting system accurately without resolving the small scale. A reader would care because most existing PNP treatments handle the two carriers independently, whereas here their coupled dynamics and the adsorption of one species at the trap surface are retained in the reduced description.

Core claim

We propose and validate a two-species Multiscale model for a Poisson-Nernst-Planck system that incorporates the simultaneous influence of both positive and negative ions interacting through the Poisson equation. The effect of the surface trap whose attraction range δ is much smaller than the macroscopic scale is replaced by a suitable boundary condition derived from mass conservation and asymptotic analysis. We also construct an asymptotic preserving and accurate scheme for the resulting MPNP system.

What carries the argument

The boundary condition for the surface trap, obtained from mass conservation and asymptotic analysis under the assumption that the attraction range δ is much smaller than the macroscopic scale, which replaces the detailed potential while preserving the essential ion dynamics and adsorption behavior.

If this is right

  • The numerical scheme remains stable and accurate even when the grid does not resolve the trap attraction range.
  • The model captures the correlated motion of both ion species and the selective adsorption of negative ions at the trap surface.
  • The asymptotic-preserving property ensures the discrete solution approaches the correct reduced model as the small scale vanishes.
  • The self-consistent Coulomb interaction between carriers is maintained through the Poisson equation at every step.

Where Pith is reading between the lines

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

  • The same boundary-condition reduction could be tested on other localized adsorption or reaction problems whose interaction length is small compared with the domain size.
  • Extending the scheme to time-dependent traps or moving boundaries would check whether the asymptotic preservation holds under more complex geometry.
  • Comparison with Monte-Carlo particle simulations of the same ion-trap setup could quantify how much information is retained by the continuum reduced model.

Load-bearing premise

The boundary condition is derived under the assumption that the attraction range δ is much smaller than the macroscopic scale, so the detailed potential can be replaced without losing key dynamics.

What would settle it

Numerical solutions of the full Poisson-Nernst-Planck system with a resolved but small finite δ compared against solutions of the reduced model with the proposed boundary condition, checking whether the ion densities and fluxes converge as δ approaches zero.

Figures

Figures reproduced from arXiv: 2507.01402 by Clarissa Astuto, Giovanni Russo.

Figure 1
Figure 1. Figure 1: Scheme of the experimental setup. On the top left there is a zoom in of the anions and cations behavior at the air-surface of the bubble: the cations (blue) are composed by hydrophilic heads; the anions (red) have hydrophobic tails inside the air bubble, and hydrophilic heads on the surface. In this work, we consider the PNP model for the diffusion of the two carriers [26, 41], in presence of an external t… view at source ↗
Figure 2
Figure 2. Figure 2: Example of the potentials V±(x) and, after a change of variable, the corresponding U±(ξ) for δ = 10−2 . Now we write the system for the interval Ωδ f = [δL, 1] ∂c± ∂t = − ∂J± ∂x , in Ωδ f (11a) J± = −D±  ∂c± ∂x ± c± ∂Φ ∂x  , in Ωδ f (11b) −ε ∂ 2Φ ∂x2 = c+ m+ − c− m− , in Ωδ f (11c) J±|x=1 = 0 (11d) ∂Φ ∂x [PITH_FULL_IMAGE:figures/full_fig_p006_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Space and time accuracy orders of the δ-model at final time t = 0.1, for different values of δ and ε. Simulations details are provided in Section 4. The discretization of (65a) at the left boundary reads M 2  dc−,0 dt + dc−,1 dt  = D−  c−,1 − c−,0 h − c−,1 + c−,0 2 Φ1 − Φ0 h  . (66) Now we discretize the boundary condition for Φ, that at the left boundary becomes ε Φ1 − Φ0 h = M c−,0 + c−,1 2 . (67) Eq… view at source ↗
Figure 4
Figure 4. Figure 4: Space and time accuracy orders of the 0-model at final time t = 0.1, for different values of ε. Simulations details are provided in Section 4. for different values of N and show the results in [PITH_FULL_IMAGE:figures/full_fig_p017_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Difference between the solutions c δ − and c 0 − of the δ-model and 0-model, respectively. The quantity is computed is Eq. (77). 17 [PITH_FULL_IMAGE:figures/full_fig_p017_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Qualitative comparison between the solutions c δ ± and c 0 ± of the δ-model and the 0-model, respec￾tively, for different values of δ. Zoom-in in panels (e),(d). 18 [PITH_FULL_IMAGE:figures/full_fig_p018_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Space and time accuracy orders of the QNL system (see Eqs. (78-79)) at final time t = 0.1, for different values of ε. Simulations details are provided in Section 4. In these tests ∆t = 0.1∆x. with boundary conditions M 2 ∂C ∂t = De ∂C ∂n + DbC ∂Φ ∂n on ΓB (81a) − M 2 ∂C ∂t = Db ∂C ∂n + DeC ∂Φ ∂n on ΓB (81b) ∂Φ ∂n = M 2ε C − M 2 Q on ΓB (81c) ∂C ∂n = 0 on ΓS (81d) ∂Φ ∂n = 0 on ΓS (81e) Numerical tests in on… view at source ↗
Figure 8
Figure 8. Figure 8: Discretization of the computational domain. Ω is the green region inside the unit square R. (a): classification of the grid points: the blue points are the internal ones while the red circles denote the ghost points. (b): points of intersection between the grid and the boundary Γ (see the definition of A and B in Algorithm 1). Here, we define the set of ghost points G, which are grid points that belong to … view at source ↗
Figure 9
Figure 9. Figure 9: Grid before and after snapping technique. (a): representation of the cell related to the internal point P (blue points), whose distance from Γ is less than h 2 ; (b): zoom-in of the shape of the domain, after the grid point P has changed its classification, from internal to ghost point (red circles). where, for each K, φi |K is the restriction of φi in cell K, which is a bilinear function and takes value 1… view at source ↗
Figure 10
Figure 10. Figure 10: Scheme of the three quadrature points (circles) for each edge li , i = 0, · · · , 4. The squared points represent the vertices Pi , i = 0, · · · , 4, of the polygon P. Let us consider a general integrable function f(x, y) defined in Ω, with F(x, y) = R f(x, y)dx (in our strategy we consider the primitive in x direction; analogue results can be obtained integrating in y direction). 24 [PITH_FULL_IMAGE:fig… view at source ↗
Figure 11
Figure 11. Figure 11: The AP diagram. P ε is the original problem and P ε h its numerical approximation characterized by a discretization parameter h. The AP property corresponds to the request that P 0 h is consistent with P 0 as ε → 0, independently of h. To reformulate system (89) using the computational matrices for the spatial derivatives, we express it as follows B[vh] ∂Ch ∂t −  M 2 ∂Ch ∂t − ε M 2 ∂Qh ∂t , vh  L2(ΓB,h)… view at source ↗
Figure 12
Figure 12. Figure 12: Charge conservation test: we plot the difference diffcons defined in Eq. (104), as a function of the number of cells of the space discretization (a). In panel (b), we show the same quantity in function of time, considering an explicit discretization in time. In (b) ∆t = 0.1h 2 . The initial conditions are given by C in(x, t = 0) = c in +(x, t = 0) m+ + c in −(x, t = 0) m− (102a) Q in(x, t = 0) = 1 ε  c i… view at source ↗
Figure 13
Figure 13. Figure 13: Space and time accuracy of the QNL system in two dimensions (see Eqs. (98–101)) at final time t = 0.3125, for different values of ε. In this test the order is calculated with Richardson extrapolation technique. Simulation details are provided in Section 7. Charge conservation If zero flux boundary conditions are adopted on the external boundary, system (59), (60) conserves the total charges of both anions… view at source ↗
Figure 14
Figure 14. Figure 14: Time accuracy orders of the QNL system in two dimensions (see Eqs. (98–101)) at final time t = 0.1, for different values of ε. Simulation details are provided in Section 7. In panel (d) the error in Φ is not plotted, however for vanishing ε, the value of the potential Φ remains constant within machine precision. 31 [PITH_FULL_IMAGE:figures/full_fig_p031_14.png] view at source ↗
Figure 15
Figure 15. Figure 15: Profiles of the concentrations c+ (solid lines) and c− (dashed lines) in two dimensions, at x = 0.5. We show the solutions at different times t = 5, 10, 15, 20, and for different ε. In panels (a)-(d) the initial volume is v0 = 10−6 while in panel (e) v0 = 10−11 and ∆t = 0.01h. 32 [PITH_FULL_IMAGE:figures/full_fig_p032_15.png] view at source ↗
Figure 16
Figure 16. Figure 16: Contour plot of the difference defined in Eq. (105), for different values of ε. in regimes where the Debye length is small but not negligible. While the Quasi-Neutral limit provides a simplified model in the asymptotic regime ε = 0, the case ε ≪ 1 introduces significant numerical challenges: the system becomes stiff and the condition number of the discretized matrix increases, leading to potential loss of… view at source ↗
Figure 17
Figure 17. Figure 17: Time evolution of anion concentration c− at different times t = 5, 10, 15, and 20. We mark in red the concentration values at the boundary of the bubble ΓB,h. 40 [PITH_FULL_IMAGE:figures/full_fig_p040_17.png] view at source ↗
Figure 18
Figure 18. Figure 18: Positivity of the solution for different values of ε. (a) (b) [PITH_FULL_IMAGE:figures/full_fig_p041_18.png] view at source ↗
Figure 19
Figure 19. Figure 19: Positivity of the solution for different initial conditions. 41 [PITH_FULL_IMAGE:figures/full_fig_p041_19.png] view at source ↗
read the original abstract

In this paper, we propose and validate a two-species Multiscale model for a Poisson-Nernst-Planck (PNP) system, focusing on the correlated motion of positive and negative ions under the influence of a trap. Specifically, we aim to model surface traps whose attraction range, of length $\delta$, is much smaller then the scale of the problem. The physical setup we refer to is an anchored gas drop (bubble) surrounded by a flow of charged surfactants {(composed by positive and negative ions) that diffuses in water. When the diffusing surfactants reach the surface of the trap, the negative ions are adsorbed because of their hydrophobic tail that is attracted by the air bubble}. As in our previous works, the effect of the attractive potential is replaced by a suitable boundary condition derived by mass conservation and asymptotic analysis. The novelty of this work is the extension of the model proposed in \cite{astuto2023multiscale}, now incorporating the influence of both carriers -- positive and negative ions -- simultaneously, which is often neglected in traditional approaches that treat ion species independently. The two carriers interact through the Coulomb potential, that is computed by a Poisson equation. [...]

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

1 major / 2 minor

Summary. The paper proposes and validates a two-species multiscale Poisson-Nernst-Planck (MPNP) model for charged surfactants near a surface trap (e.g., an anchored gas bubble) whose short-range attractive potential has length δ much smaller than the macroscopic scale. The attractive potential is replaced by an effective boundary condition derived from mass conservation and asymptotic analysis; this extends the single-species model of Astuto et al. (2023) by retaining both positive and negative ions coupled through the self-consistent Poisson potential. An asymptotic-preserving numerical scheme is constructed and tested for accuracy and preservation properties.

Significance. If the two-species effective boundary condition is shown to capture the correct leading-order adsorption rates without missing cross-coupling terms, the work would supply a practical, asymptotically consistent tool for simulating multiscale PNP systems that avoids resolving the microscopic attraction length δ. The explicit treatment of both carriers addresses a frequent modeling simplification in the literature.

major comments (1)
  1. [§2.3] §2.3 (Derivation of the effective boundary condition): The single-species predecessor replaces the short-range potential by a local flux condition obtained from a separable boundary-layer problem. In the two-species extension the Poisson equation couples the densities of both carriers, so the leading-order boundary-layer problem for the normal fluxes is no longer separable. The manuscript does not exhibit the explicit two-species boundary-layer calculation; it is therefore unclear whether additional cross-terms of the same order appear in the effective adsorption rates for the negative ions.
minor comments (2)
  1. The statement that the scheme is 'parameter-free' should be qualified by listing any numerical parameters (e.g., stabilization constants or mesh-size ratios) that remain after the asymptotic analysis.
  2. Figure 4 (error plots) would benefit from a log-log inset confirming the expected convergence rate under mesh refinement.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the careful reading of our manuscript and the constructive comment on the derivation of the effective boundary condition. We address the point below and will incorporate clarifications in the revision.

read point-by-point responses

  1. Referee: [§2.3] §2.3 (Derivation of the effective boundary condition): The single-species predecessor replaces the short-range potential by a local flux condition obtained from a separable boundary-layer problem. In the two-species extension the Poisson equation couples the densities of both carriers, so the leading-order boundary-layer problem for the normal fluxes is no longer separable. The manuscript does not exhibit the explicit two-species boundary-layer calculation; it is therefore unclear whether additional cross-terms of the same order appear in the effective adsorption rates for the negative ions.

    Authors: We appreciate the referee pointing out the need for explicit details on the two-species boundary-layer analysis. In our derivation, the leading-order inner problem is formulated as a coupled system for the densities and the potential, solved subject to the short-range potential and far-field matching. Because the Poisson equation is linear in the leading-order correction and the boundary-layer coordinate is stretched by δ, the cross-coupling enters only through the self-consistent potential at higher order; the leading-order flux conditions for each species remain independent of additional cross-terms of O(1). To make this transparent, we will add the full two-species boundary-layer calculation (including the explicit solution of the inner problem and the resulting effective adsorption rates) to Section 2.3 in the revised manuscript. revision: yes

Circularity Check

0 steps flagged

Minor self-citation to single-species predecessor; central BC derivation and scheme remain independent

full rationale

The paper extends the single-species model from the authors' prior work via citation but grounds the two-species boundary condition explicitly in mass conservation plus asymptotic analysis of the short-range trap potential. No fitted parameter is renamed as a prediction, no self-referential definition appears in the equations, and the numerical scheme is constructed and validated separately from the cited predecessor. The self-citation supplies context rather than load-bearing justification for the core result.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the validity of the asymptotic reduction of the trap potential to a boundary condition and on standard assumptions of the PNP system; no new free parameters or invented entities are introduced beyond the model extension.

axioms (1)
  • domain assumption The attraction range δ is much smaller than the macroscopic scale, justifying replacement of the potential by a boundary condition via asymptotic analysis and mass conservation.
    Invoked in the abstract to derive the trap model from the physical setup.

pith-pipeline@v0.9.0 · 5748 in / 1311 out tokens · 28532 ms · 2026-05-19T06:58:10.981718+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

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.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Standard versus Asymptotic Preserving Time Discretizations for the Poisson-Nernst-Planck System in the Quasi-Neutral Limit

    math.NA 2025-11 unverdicted novelty 4.0

    IMEX time discretizations for the PNP system remain asymptotically stable for all Debye lengths and require no special assumptions on initial conditions, unlike standard schemes.

Reference graph

Works this paper leans on

63 extracted references · 63 canonical work pages · cited by 1 Pith paper

  1. [1]

    A cut finite element method for non-newtonian free surface flows in 2d-application to glacier modelling

    Josefin Ahlkrona and Daniel Elfverson. A cut finite element method for non-newtonian free surface flows in 2d-application to glacier modelling. Journal of Computational Physics: X , 11:100090, 2021

    work page 2021

  2. [2]

    Zero-electron-mass and quasi-neutral limits for bipolar euler–poisson systems

    Nuno J Alves and Athanasios E Tzavaras. Zero-electron-mass and quasi-neutral limits for bipolar euler–poisson systems. Zeitschrift f¨ ur angewandte Mathematik und Physik, 75(1):17, 2024

    work page 2024

  3. [3]

    Implicit-explicit runge-kutta methods for time-dependent partial differential equations

    Uri M Ascher, Steven J Ruuth, and Raymond J Spiteri. Implicit-explicit runge-kutta methods for time-dependent partial differential equations. Applied Numerical Mathematics, 25(2-3):151–167, 1997

    work page 1997

  4. [4]

    High order multiscale methods for advection-diffusion equation in highly oscillatory regimes: application to surfactant diffusion and generalization to arbitrary domains

    Clarissa Astuto. High order multiscale methods for advection-diffusion equation in highly oscillatory regimes: application to surfactant diffusion and generalization to arbitrary domains. Communications in Computational Physics , 37(5):1327–1357, 2025

    work page 2025

  5. [5]

    A nodal ghost method based on variational formulation and regular square grid for elliptic problems on arbitrary domains in two space dimensions

    Clarissa Astuto, Daniele Boffi, Giovanni Russo, and Umberto Zerbinati. A nodal ghost method based on variational formulation and regular square grid for elliptic problems on arbitrary domains in two space dimensions. Computer Methods in Applied Mechanics and Engineering , 443:118041, 2025

    work page 2025

  6. [6]

    A finite-difference ghost-point multigrid method for multi-scale modelling of sorption kinetics of a surfactant past an oscillating bubble

    Clarissa Astuto, Armando Coco, and Giovanni Russo. A finite-difference ghost-point multigrid method for multi-scale modelling of sorption kinetics of a surfactant past an oscillating bubble. Journal of Computational Physics, 476:111880, 2023

    work page 2023

  7. [7]

    A comparison of the coco-russo scheme and ghost-fem for elliptic equations in arbitrary domains

    Clarissa Astuto, Armando Coco, and Umberto Zerbinati. A comparison of the coco-russo scheme and ghost-fem for elliptic equations in arbitrary domains. arXiv preprint arXiv:2405.16582 , 2024

    work page arXiv 2024

  8. [8]

    Self-regulated biological trans- portation structures with general entropy dissipations, part I: The 1D case

    Clarissa Astuto, Jan Haskovec, Peter Markowich, and Simone Portaro. Self-regulated biological trans- portation structures with general entropy dissipations, part I: The 1D case. Journal of Dynamics and Games, 2023

    work page 2023

  9. [9]

    Time multiscale modeling of sorption kinetics i: uniformly accurate schemes for highly oscillatory advection-diffusion equation

    Clarissa Astuto, Mohammed Lemou, and Giovanni Russo. Time multiscale modeling of sorption kinetics i: uniformly accurate schemes for highly oscillatory advection-diffusion equation. Multiscale Modeling & Simulation, 23(2):753–776, 2025

    work page 2025

  10. [10]

    Self-regulated biological transportation structures with general entropy dissipation: 2D case and leaf-shaped domain

    Clarissa Astuto, Peter Markowich, Simone Portaro, and Giovanni Russo. Self-regulated biological transportation structures with general entropy dissipation: 2D case and leaf-shaped domain. arXiv preprint arXiv:2408.15680, 2024

    work page arXiv 2024

  11. [11]

    Multiscale modeling of sorption kinetics

    Clarissa Astuto, Antonio Raudino, and Giovanni Russo. Multiscale modeling of sorption kinetics. Multiscale Modeling & Simulation , 21(1):374–399, 2023. 35

    work page 2023

  12. [12]

    Space-time unfitted finite element methods for time-dependent problems on moving domains

    Santiago Badia, Hridya Dilip, and Francesc Verdugo. Space-time unfitted finite element methods for time-dependent problems on moving domains. Computers & Mathematics with Applications, 135:60–76, 2023

    work page 2023

  13. [13]

    An asymptotically stable semi-lagrangian scheme in the quasi-neutral limit

    Radoin Belaouar, Nicolas Crouseilles, Pierre Degond, and Eric Sonnendr¨ ucker. An asymptotically stable semi-lagrangian scheme in the quasi-neutral limit. Journal of Scientific Computing , 41:341–365, 2009

    work page 2009

  14. [14]

    Study of a finite volume scheme for the drift-diffusion system

    Marianne Bessemoulin-Chatard, Claire Chainais-Hillairet, and M-H Vignal. Study of a finite volume scheme for the drift-diffusion system. asymptotic behavior in the quasi-neutral limit. SIAM Journal on Numerical Analysis, 52(4):1666–1691, 2014

    work page 2014

  15. [15]

    High order semi-implicit schemes for time dependent partial differential equations

    Sebastiano Boscarino, Francis Filbet, and Giovanni Russo. High order semi-implicit schemes for time dependent partial differential equations. Journal of Scientific Computing , 68:975–1001, 2016

    work page 2016

  16. [16]

    Implicit-explicit methods for evolutionary partial differential equations

    Sebastiano Boscarino, Lorenzo Pareschi, and Giovanni Russo. Implicit-explicit methods for evolutionary partial differential equations. SIAM, 2024

    work page 2024

  17. [17]

    Foam fractionation of globular proteins

    Lynn Brown, Ganesan Narsimhan, and PC Wankat. Foam fractionation of globular proteins. Biotech- nology and Bioengineering, 36(9):947–959, 1990

    work page 1990

  18. [18]

    Asymptotic-preserving scheme for a bi-fluid euler-lorentz model

    St´ ephane Brull, Pierre Degond, Fabrice Deluzet, and Alexandre Mouton. Asymptotic-preserving scheme for a bi-fluid euler-lorentz model. Kinetic and Related Models , 4(4):991–1023, 2011

    work page 2011

  19. [19]

    Cutfem: discretizing geometry and partial differential equations

    Erik Burman, Susanne Claus, Peter Hansbo, Mats G Larson, and Andr´ e Massing. Cutfem: discretizing geometry and partial differential equations. International Journal for Numerical Methods in Engineer- ing, 104(7):472–501, 2015

    work page 2015

  20. [20]

    Fictitious domain finite element methods using cut elements: Ii

    Erik Burman and Peter Hansbo. Fictitious domain finite element methods using cut elements: Ii. a stabilized nitsche method. Applied Numerical Mathematics, 62(4):328–341, 2012

    work page 2012

  21. [21]

    Cutfem based on extended finite element spaces

    Erik Burman, Peter Hansbo, and Mats G Larson. Cutfem based on extended finite element spaces. Numerische Mathematik, 152(2):331–369, 2022

    work page 2022

  22. [22]

    High order asymptotic preserving penalized numerical schemes for the euler-poisson system in the quasi-neutral limit

    Nicolas Crouseilles, Giacomo Dimarco, and Saurav Samantaray. High order asymptotic preserving penalized numerical schemes for the euler-poisson system in the quasi-neutral limit. arXiv preprint arXiv:2409.04807, 2024

    work page arXiv 2024

  23. [23]

    An asymptotic preserving scheme for the euler equations in a strong magnetic field.Journal of Computational Physics, 228(10):3540–3558, 2009

    Pierre Degond, Fabrice Deluzet, Afeintou Sangam, and M-H Vignal. An asymptotic preserving scheme for the euler equations in a strong magnetic field.Journal of Computational Physics, 228(10):3540–3558, 2009

    work page 2009

  24. [24]

    Multigrid methods for the ghost finite element approximation of elliptic problems

    Hridya Dilip and Armando Coco. Multigrid methods for the ghost finite element approximation of elliptic problems. arXiv preprint arXiv:2505.05105 , 2025

    work page arXiv 2025

  25. [25]

    Asymptotic preserving implicit-explicit runge–kutta methods for nonlinear kinetic equations

    Giacomo Dimarco and Lorenzo Pareschi. Asymptotic preserving implicit-explicit runge–kutta methods for nonlinear kinetic equations. SIAM Journal on Numerical Analysis , 51(2):1064–1087, 2013

    work page 2013

  26. [26]

    Poisson-nernst-planck systems for ion channels with permanent charges

    Bob Eisenberg and Weishi Liu. Poisson-nernst-planck systems for ion channels with permanent charges. SIAM journal on mathematical analysis , 2007. 36

    work page 2007

  27. [27]

    Adsorption in slit-like and cylindrical micropores in the henry’s law region

    Douglas H Everett and John C Powl. Adsorption in slit-like and cylindrical micropores in the henry’s law region. a model for the microporosity of carbons. Journal of the Chemical Society, Faraday Transactions 1: Physical Chemistry in Condensed Phases , 72:619–636, 1976

    work page 1976

  28. [28]

    Existence and uniqueness results for a kinetic model in bulk-surface surfactant dynamics

    JR Fern´ andez, Piotr Kalita, Stanislaw Migo´ orski, M Carmen Mu˜ niz, and C N´ u˜ nez. Existence and uniqueness results for a kinetic model in bulk-surface surfactant dynamics. SIAM Journal on Mathe- matical Analysis, 48(5):3065–3089, 2016

    work page 2016

  29. [29]

    Existence and unique- ness results for a kinetic model in bulk-surface surfactant dynamics

    JR Fern´ andez, Piotr Kalita, Stanislaw Mig´ orski, M Carmen Mu˜ niz, and C N´ u˜ nez. Existence and unique- ness results for a kinetic model in bulk-surface surfactant dynamics. SIAM Journal on Mathematical Analysis, 48(5):3065–3089, 2016

    work page 2016

  30. [30]

    A class of asymptotic-preserving schemes for kinetic equations and related problems with stiff sources

    Francis Filbet and Shi Jin. A class of asymptotic-preserving schemes for kinetic equations and related problems with stiff sources. Journal of Computational Physics , 229(20):7625–7648, 2010

    work page 2010

  31. [31]

    An asymptotic preserving scheme for the es-bgk model of the boltzmann equation

    Francis Filbet and Shi Jin. An asymptotic preserving scheme for the es-bgk model of the boltzmann equation. Journal of Scientific Computing , 46(2):204–224, 2011

    work page 2011

  32. [32]

    Arbitrary lagrangian–eulerian finite-element method for com- putation of two-phase flows with soluble surfactants

    Sashikumaar Ganesan and Lutz Tobiska. Arbitrary lagrangian–eulerian finite-element method for com- putation of two-phase flows with soluble surfactants. Journal of Computational Physics , 231(9):3685– 3702, 2012

    work page 2012

  33. [33]

    The initial time layer problem and the quasineutral limit in the semiconductor drift-diffusion model

    Ingenuin Gasser, C David Levermore, Peter A Markowich, and Christian Schmeiser. The initial time layer problem and the quasineutral limit in the semiconductor drift-diffusion model. European Journal of Applied Mathematics , 12(4):497–512, 2001

    work page 2001

  34. [34]

    A cut finite element method for coupled bulk-surface problems on time-dependent domains

    Peter Hansbo, Mats G Larson, and Sara Zahedi. A cut finite element method for coupled bulk-surface problems on time-dependent domains. Computer Methods in Applied Mechanics and Engineering , 307:96–116, 2016

    work page 2016

  35. [35]

    Efficient asymptotic-preserving (ap) schemes for some multiscale kinetic equations

    Shi Jin. Efficient asymptotic-preserving (ap) schemes for some multiscale kinetic equations. SIAM Journal on Scientific Computing , 21(2):441–454, 1999

    work page 1999

  36. [36]

    A hierarchy of hydrodynamic models for plasmas

    Ansgar Jungel and Yue-Jun Peng. A hierarchy of hydrodynamic models for plasmas. quasi-neutral limits in the drift-diffusion equations. Asymptotic Analysis, 28, 2000

    work page 2000

  37. [37]

    Performance and stability of the organosilicon surfactant l-77: effect of ph, concentration, and temperature

    Moritz Knoche, Hiroto Tamura, and Martin J Bukovac. Performance and stability of the organosilicon surfactant l-77: effect of ph, concentration, and temperature. Journal of agricultural and food chemistry, 39(1):202–206, 1991

    work page 1991

  38. [38]

    High order unfitted finite element methods on level set domains using isopara- metric mappings

    Christoph Lehrenfeld. High order unfitted finite element methods on level set domains using isopara- metric mappings. Computer Methods in Applied Mechanics and Engineering , 300:716–733, 2016

    work page 2016

  39. [39]

    Analysis of a high-order unfitted finite element method for elliptic interface problems

    Christoph Lehrenfeld and Arnold Reusken. Analysis of a high-order unfitted finite element method for elliptic interface problems. IMA Journal of Numerical Analysis , 38(3):1351–1387, 2018. 37

    work page 2018

  40. [40]

    Basaran, and Elias I

    Ying-Chih Liao, Osman A. Basaran, and Elias I. Franses. Effects of dynamic surface tension and fluid flow on the oscillations of a supported bubble. Colloids and Surfaces A: Physicochemical and Engineering Aspects, 282-283:183–202, 2006. A Collection of Papers in Honor of Professor Ivan B. Ivanov (Laboratory of Chemical Physics and Engineering, University...

    work page 2006

  41. [41]

    Benzhuo Lu and Y.C. Zhou. Poisson-nernst-planck equations for simulating biomolecular diffusion- reaction processes ii: Size effects on ionic distributions and diffusion-reaction rates. Biophysical Journal, 2011

    work page 2011

  42. [42]

    Shape oscillations of drops in the presence of surfactants

    Hui-Lan Lu and Robert E Apfel. Shape oscillations of drops in the presence of surfactants. Journal of fluid mechanics, 222:351–368, 1991

    work page 1991

  43. [43]

    Mathematical modelling of surfactant self-assembly at interfaces

    CE Morgan, CJW Breward, Ian M Griffiths, and Peter D Howell. Mathematical modelling of surfactant self-assembly at interfaces. SIAM Journal on Applied Mathematics , 75(2):836–860, 2015

    work page 2015

  44. [44]

    New notions on a fundamental principle of respiratory mechanics: the retractile force of the lung, dependent on the surface tension in the alveoli

    V Neergaard. New notions on a fundamental principle of respiratory mechanics: the retractile force of the lung, dependent on the surface tension in the alveoli. Zeitschrift fur Gesundh. und Exp. Medizin , 66:373–394, 1929

    work page 1929

  45. [45]

    Pulmonary surfactant: a surface chemistry viewpoint

    RH Notter and PE Morrow. Pulmonary surfactant: a surface chemistry viewpoint. Annals of biomedical engineering, 3:119–159, 1975

    work page 1975

  46. [46]

    Diffusion on a curved surface coupled to diffusion in the volume: Application to cell biology

    Igor L Novak, Fei Gao, Yung-Sze Choi, Diana Resasco, James C Schaff, and Boris M Slepchenko. Diffusion on a curved surface coupled to diffusion in the volume: Application to cell biology. Journal of computational physics , 226(2):1271–1290, 2007

    work page 2007

  47. [47]

    Osher and Ronald P

    Stanley. Osher and Ronald P. Fedkiw. Level set methods and dynamic implicit surfaces . Applied mathematical sciences ; 153. Springer, New York ;, 2002

    work page 2002

  48. [48]

    Asymptotic preserving monte carlo methods for the boltzmann equation

    Lorenzo Pareschi and Giovanni Russo. Asymptotic preserving monte carlo methods for the boltzmann equation. Transport Theory and Statistical Physics , 29(3-5):415–430, 2000

    work page 2000

  49. [49]

    Implicit-explicit runge-kutta schemes for stiff systems of differ- ential equations

    Lorenzo Pareschi and Giovanni Russo. Implicit-explicit runge-kutta schemes for stiff systems of differ- ential equations. Recent trends in numerical analysis , 3:269–289, 2000

    work page 2000

  50. [50]

    High order asymptotically strong-stability-preserving methods for hyperbolic systems with stiff relaxation

    Lorenzo Pareschi and Giovanni Russo. High order asymptotically strong-stability-preserving methods for hyperbolic systems with stiff relaxation. In Hyperbolic Problems: Theory, Numerics, Applications: Proceedings of the Ninth International Conference on Hyperbolic Problems held in CalTech, Pasadena, March 25–29, 2002 , pages 241–251. Springer, 2003

    work page 2002

  51. [51]

    Implicit–explicit runge–kutta schemes and applications to hy- perbolic systems with relaxation

    Lorenzo Pareschi and Giovanni Russo. Implicit–explicit runge–kutta schemes and applications to hy- perbolic systems with relaxation. Journal of Scientific computing , 25:129–155, 2005

    work page 2005

  52. [52]

    Numerical simulations of a rising drop with shape oscillations in the presence of surfactants

    Antoine Piedfert, Benjamin Lalanne, Olivier Masbernat, and Fr´ ed´ eric Risso. Numerical simulations of a rising drop with shape oscillations in the presence of surfactants. Physical Review Fluids, 3(10):103605, 2018. 38

    work page 2018

  53. [53]

    Anomalous sorption kinetics of self-interacting particles by a spherical trap

    Antonio Raudino, Antonio Grassi, Giuseppe Lombardo, Giovanni Russo, Clarissa Astuto, and Mario Corti. Anomalous sorption kinetics of self-interacting particles by a spherical trap. Communications in Computational Physics, 31(3):707–738, 2022

    work page 2022

  54. [54]

    Lewis Fry Richardson. IX. the approximate arithmetical solution by finite differences of physical prob- lems involving differential equations, with an application to the stresses in a masonry dam.Philosophical Transactions of the Royal Society of London. Series A, Containing Papers of a Mathematical or Physical Character, 210(459-470):307–357, 1911

    work page 1911

  55. [55]

    Code verification by the method of manufactured solutions

    Patrick J Roache. Code verification by the method of manufactured solutions. J. Fluids Eng., 124(1):4– 10, 2002

    work page 2002

  56. [56]

    A remark on computing distance functions

    Giovanni Russo and Peter Smereka. A remark on computing distance functions. Journal of computa- tional physics, 163(1):51–67, 2000

    work page 2000

  57. [57]

    J. A. Sethian. Level set methods and fast marching methods: evolving interfaces in computational geometry, fluid mechanics, computer vision, and materials science . Cambridge monographs on applied and computational mathematics 3. Cambridge University Press, 2nd ed edition, 1999

    work page 1999

  58. [58]

    A level set approach for computing solutions to incompressible two-phase flow

    Mark Sussman, Peter Smereka, and Stanley Osher. A level set approach for computing solutions to incompressible two-phase flow. Journal of computational physics , 114(1):146–159, 1994

    work page 1994

  59. [59]

    Unsteady evolution of slip and drag in surfactant-contaminated superhydrophobic channels

    Samuel Tomlinson, Fr´ ed´ eric Gibou, Paolo Luzzatto-Fegiz, Fernando Temprano-Coleto, Oliver Jensen, and Julien R Landel. Unsteady evolution of slip and drag in surfactant-contaminated superhydrophobic channels. Journal of Fluid Mechanics , 2023

    work page 2023

  60. [60]

    Laminar drag reduction in surfactant-contaminated superhydrophobic channels

    Samuel D Tomlinson, Fr´ ed´ eric Gibou, Paolo Luzzatto-Fegiz, Fernando Temprano-Coleto, Oliver E Jensen, and Julien R Landel. Laminar drag reduction in surfactant-contaminated superhydrophobic channels. Journal of Fluid Mechanics , 963:A10, 2023

    work page 2023

  61. [61]

    Role of dynamic surface tension in slide coating

    Jose E Valentini, William R Thomas, Paul Sevenhuysen, Tsung S Jiang, Hae O Lee, Yi Liu, and Shi Chern Yen. Role of dynamic surface tension in slide coating. Industrial & engineering chemistry research, 30(3):453–461, 1991

    work page 1991

  62. [62]

    Wesseling

    P. Wesseling. Principles of computational fluid dynamics . Springer, 2001

    work page 2001

  63. [63]

    Surfactant solutions new methods of investigations

    R Zana. Surfactant solutions new methods of investigations. Surfactant science series , 1 1986. 39 (a) (b) (c) (d) Figure 17: Time evolution of anion concentration c− at different times t = 5, 10, 15, and 20. We mark in red the concentration values at the boundary of the bubble ΓB,h. 40 (a) (b) Figure 18: Positivity of the solution for different values of...

    work page 1986