Pith. sign in

REVIEW 3 major objections 7 minor 34 references

Phase-Field Modeling and Energy-Stable Schemes for Osmotic Flow through Semi-Permeable

T0 review · 3 major / 7 minor · reviewed 2026-08-07 · deepseek-v4-flash

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

Pith's one-line read A thermodynamically consistent phase-field model for osmotic flow across semi-permeable membranes, with discretizations that preserve the energy-dissipation law.

desk verdict A genuinely useful new phase-field model with a clean continuous energy law, but the discrete stability theorem does not cover the potential used in the code, and the SDC extension is claimed without proof. read the letter →

arxiv 2506.11374 v1 pith:5FIGPVWA submitted 2025-06-13 physics.flu-dyn cs.NAmath.NA

Ruihan Guo , Jie Shen , Shixin Xu , Xianmin Xu This is my paper · ORCID
classification physics.flu-dyncs.NAmath.NA MSC 65M6065H1076T0692C10
keywords phase-fieldmodelsemi-permeablemembraneosmoticpressureenergydissipationlawNavier-Stokes-Cahn-Hilliard-Allen-CahnlocaldiscontinuousGalerkinspectraldeferredcorrectiondouble-wellpotential
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 sets out to show that osmotic water transport across a semi-permeable membrane can be described by a phase-field model that is thermodynamically consistent: the total free energy, combining kinetic, entropic, and mixing contributions, can never increase along solutions. The model augments the Navier–Stokes–Cahn–Hilliard system with an Allen–Cahn-type transmembrane flux proportional to the chemical-potential imbalance, so water can cross the diffuse interface while solute cannot. The authors then construct numerical schemes, a first-order decoupled scheme and a spectral-deferred-correction lift to higher order, and prove discrete energy-dissipation laws for both. If correct, the framework gives a computational tool for osmotic effects in biology and industry where the simulated dynamics obey the second law of thermodynamics by construction.

What carries the argument

The machinery is the free-energy functional $E = E_{\mathrm{kin}} + E_{\mathrm{entropy}} + E_{\mathrm{mix}}$ with chemical potential $\tilde{\mu}_\phi = -\epsilon\Delta\phi + \frac{1}{\epsilon}F'(\phi) - \frac{\beta}{\epsilon}(\partial_\phi \zeta_+ C_+ + \partial_\phi \zeta_- C_-)$, and the Allen-Cahn-type transmembrane flux $S_\phi = -\frac{K}{Ca}\tilde{\mu}_\phi |\nabla\phi|$ that lets the interface move relative to the fluid. Its discrete support comes from a modified intermediate velocity $u_*$ that decouples the phase-field and concentration updates from the momentum update, and from an alternating-flux choice in the LDG scheme that makes the numerical operators adjoint-consistent so the telescoping energy argument survives discretization. The stabilization term $\frac{S}{\epsilon}(\phi^{n+1}-\phi^n)$ with $S \geq L/2$ is what makes the first-order semi-implicit treatment of $F(\phi)$ unconditionally energy-stable under condition (16).

What would settle it

Run the LDG scheme (30)-(32) with the untruncated potential and initial data whose phase field takes values outside $[-1,1]$, so that $|F''|$ exceeds any chosen $L$, and check whether the discrete energy $E_h^n$ is monotone non-increasing; a single time step where $E_h^{n+1} > E_h^n$ would refute the claimed unconditional stability as stated.

Watch

Extended reading notes

Core claim

On the paper's own terms, the central discovery is that the osmotic coupling can be encoded in a single extra flux term, $S_\phi = -\frac{K}{Ca}\tilde{\mu}_\phi |\nabla\phi|$, that preserves energy dissipation: Theorem 2.1 proves for the continuous system (5) that the time derivative of the total free energy equals a sum of negative semi-definite dissipation terms. The same dissipation structure is shown to survive time discretization (Theorem 3.1) and full local-discontinuous-Galerkin discretization (Theorem 3.2), provided the stabilizing constant $S$ is at least half the maximum of $|F''|$. The equilibrium condition $u=0$, $\tilde{\mu}_\phi=0$, $\mu_c^\pm=0$ then predicts that droplet equilibrium shape is determined by surface tension and solute loading, not by membrane permeability, which the simulations in Section 4 confirm.

Load-bearing premise

The load-bearing premise is that the double-well potential has bounded second derivative (condition (16)), but the actual potential $F(\phi)=\frac14(\phi^2-1)^2$ has $F''(\phi)=3\phi^2-1$, which is unbounded over all real $\phi$; the paper appeals to the common practice of using a truncated double-well potential for the proofs while the numerical scheme and tests use the untruncated one.

Editorial extensions

If this is right

  • At equilibrium the model predicts $u = 0$, $\tilde{\mu}_\phi = 0$, and $\mu_c^\pm = 0$, so droplet shape is set by surface tension and total solute mass, not by permeability; higher permeability only accelerates approach to equilibrium.
  • The first-order scheme and its SDC acceleration are decoupled: each step solves the phase-field and concentration system, then a Stokes-type system, then a pressure correction, which makes large three-dimensional simulations practical.
  • The discrete energy law includes a pressure-stabilization term $\frac{\Delta t^2}{2Re}\|\nabla p\|^2$, so the projection step contributes positively to the energy identity rather than spoiling it.
  • When osmotic effects are omitted, the model reduces to a degenerate Allen-Cahn-driven collapse in which droplets vanish; with solute present the droplet instead shrinks to a finite equilibrium size.
  • Permeability changes qualitative outcomes: in shear flow, impermeable droplets coalesce while permeable droplets shrink and avoid contact.

Reading between the lines

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

  • If the membrane permeability $K$ is made spatially varying or time-dependent, the same energy argument should yield a dissipation law with a time-dependent coefficient, because $K$ enters only the non-negative dissipative term; this extension is not stated in the paper.
  • The proof relies on replacing $F$ by a truncated potential, so a practical user should monitor $\max|\phi|$ over the run; a testable hypothesis is that the unconditionally stable behavior persists as long as $\phi$ stays in the interval where $F''$ is bounded by the chosen $S$.
  • The decoupled structure via $u_*$ suggests a natural splitting for adding other physics, such as electrodiffusion or active membrane stresses, because osmotic feedback enters only through the modified advection velocity.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

3 major / 7 minor

Summary. Section 2 derives, via the energetic variational approach, a Navier–Stokes–Cahn–Hilliard–Allen–Cahn model (5) for two-phase flow across semi-permeable membranes, with a transmembrane flux proportional to the chemical-potential imbalance through the K|∇φ|μ̃φ term; Theorem 2.1 verifies the continuous energy dissipation law. Section 3 presents a first-order, stabilized temporal scheme (11)–(14), proves discrete energy dissipation under the bounded-second-derivative assumption (16) in Theorem 3.1, extends the result to a fully discrete LDG scheme in Theorem 3.2, and constructs a semi-implicit SDC method for higher temporal accuracy. Section 4 reports convergence orders (Tables 1–2), energy decay (Figure 1), and numerical studies of membrane permeability, osmotic pressure, and shear flow on droplet equilibria (Figures 2–9). The paper's central claims are the thermodynamic consistency of the model, the energy stability of the first-order scheme and its LDG discretization, and the high-order accuracy of the SDC–LDG combination.

Significance. The continuous model is built transparently from an explicit energy and dissipation functional, so Theorem 2.1 is a genuine verification rather than a fitted identity, and the discrete stability proofs for the first-order and LDG schemes are carried out in detail with correct algebraic structure. The manuscript also provides a clear falsifiable prediction—the equilibrium droplet shape depends on solute content but not on permeability K—and the numerical experiments support it. If the gap described in the major comments (assumption (16) versus the implemented untruncated potential) is closed, the paper would be a sound technical contribution to phase-field modeling of osmotic flow, with value for both the model and the numerical framework. The significance is moderate: the ideas are largely assembled from well-established building blocks (stabilized semi-implicit treatment, alternating-flux LDG, SDC), but their combination for this NSCHAC system is new.

major comments (3)
  1. [3.1, condition (16); Theorems 3.1, 3.2; scheme (11); Section 4] The unconditional stability claim (restated in Section 5) is not established for the implemented scheme. Theorems 3.1 and 3.2 assume condition (16), max|F''(φ)| ≤ L, but the temporal scheme (11) and all tests in Section 4 use F(φ)=1/4(φ²−1)², whose second derivative 3φ²−1 is unbounded on ℝ, so no finite L exists. The paragraph after (16) acknowledges this and cites the common practice of truncation, yet the manuscript never states that a truncated potential is actually used in the experiments, and no L∞ bound on the discrete φ that would keep F''(φ) ≤ 2S is established; the discrete energy laws (27) and (53) require (S − f'(ξ^n)/2)(φ^{n+1}−φ^n)² ≥ 0, which forces f'(ξ^n) ≤ 2S at every step. Since the Cahn-Hilliard equation admits no general maximum principle, the proofs as written apply only to a potential different from the one implemented. The authors should either (a) implement and explicitly define the truncated potential, giving the resulting L and S, or (b) prove or verify that the computed solutions stay in a region where (16) holds, and in either case state the stabilization parameter S used in Section 4.
  2. [4.4, Example 4.4; boundary conditions (4)] The shear-flow experiment imposes velocity boundary data u|y=0=(−1,0)^T and u|y=6.4=(1,0)^T, whereas the model (5) and the stability theorems (Theorem 2.1 and Theorem 3.2) assume the homogeneous condition u|∂Ω=0 from (4); the LDG flux convention (33) specifies bun+1_p=0 and inherits the remaining boundary fluxes from the interior, so it is not stated how non-homogeneous Dirichlet velocity data is enforced. As a result, the energy-stability guarantee does not cover Example 4.4, and the text should say so explicitly, along with a description of how the moving-wall boundary conditions are implemented in the LDG scheme.
  3. [4, Examples 4.1–4.4] The stabilization parameter S, which the theorems require to satisfy S ≥ L/2, is never reported in Section 4, so the stability hypothesis cannot be checked for the computations of Tables 1–2 and Figures 1–9; the value of S (and, if a truncated potential is used, its definition) must be given for each experiment.
minor comments (7)
  1. [A, Appendix] In the derivation of I2 + I3, the terms 'ζ− ln C+ c∞ ∂C−/∂t' contain typos (C+ should be C− in the ζ− terms); please correct the logarithms and make the notation consistent with the objective functional (2).
  2. [Table 1] In Table 1, the C+ row for m=3 reports an L2 error of 5.82E-04, which is larger than the m=2 value of 1.15E-04, yet the order column lists 0.98; this looks like a typo (probably 5.82E-05) and should be corrected, since the current entry contradicts the claimed first-order convergence.
  3. [3.1, scheme (11)] The quotient (ζ^{n+1}_±−ζ^n_±)/(φ^{n+1}−φ^n) is undefined when φ^{n+1}=φ^n; the paper should state the convention used (e.g., the limiting value of the secant derivative).
  4. [3.1] The scheme (11)–(14) is called 'decoupled,' but Step 1 is a coupled system for φ^{n+1} and C^{n+1}_± through the ζ-derivative term and through u^n_*; a sentence clarifying that the decoupling is with respect to the pressure-velocity system (Steps 2–3) would be helpful.
  5. [Abstract and Section 5] The abstract and conclusion describe 'high-order, energy-stable numerical schemes' as a whole, but the energy-stability proofs (Theorems 3.1, 3.2) apply only to the first-order scheme and its LDG version; no stability proof is given for the SDC scheme, whose energy decay is shown only numerically in Figure 1(b). Please rephrase to avoid implying a stability proof for the SDC method.
  6. [4.3] The stated equilibrium conditions 'u = 0, μ̃φ = 0, μ±c = 0' do not follow from the dissipation law of Theorem 2.1: the last dissipation term vanishes already when ∇μ±c = 0 on each phase subdomain, and the equilibrium concentrations are fixed by total solute mass; the simulations' equilibrium concentrations are not close to C± = 1 implied by μ±c = ln C± = 0. The text should state the gradient condition rather than μ±c = 0.
  7. [2, dimensionless form] In the dimensionless equations (5), the diffusivities D± appear as bare ratios, but their normalization is not defined after the scaling; please define the dimensionless D± (e.g., D±/D∗).

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the model is derived from an explicit energy/dissipation variational principle, the discrete stability theorems are proven from the scheme, and self-citations are tool citations rather than load-bearing premises.

full rationale

The central energy-dissipation law (Theorem 2.1) is not circular: the model (5) is obtained by specifying the fluxes and stresses so that the time derivative of the energy (2) matches the dissipation functional (3), and the theorem algebraically verifies this constitutive closure. This is the standard energetic variational construction, and the paper makes the construction transparent in Section 2 and Appendix A. The discrete energy-stability results (Theorems 3.1 and 3.2) are proved directly from the semi-implicit scheme using explicit Taylor expansions; condition (16) is stated as a hypothesis of the theorems. The paper explicitly acknowledges that F(phi)=1/4(phi^2-1)^2 does not satisfy (16) and cites [3] for the common truncated-potential practice; whether the reported numerical tests actually use the truncated potential is a correctness gap, not a circular argument. References [15,16,17] are self-citations for the LDG and SDC methods used as tools, and no load-bearing uniqueness theorem or fitted parameter is imported from them. The convergence studies are performed against manufactured exact solutions, so the numerical claims are independently checkable. No quoted step reduces to its own input; the only notable concern is the unverified applicability of the stability theorems to the untruncated potential in the experiments, which is a rigor gap rather than circularity.

Assumptions & free parameters 5 free parameters · 5 assumptions · 0 invented entities

The model relies on an explicitly chosen energy and dissipation functional, which are constitutive assumptions. The numerical stability proofs add the truncated-potential assumption. No new physical entities are introduced. The remaining items are standard mathematical and domain assumptions for phase-field modeling.

free parameters (5)
  • K (dimensionless membrane permeability) = 0, 0.1, 0.3
    Chosen by hand in Examples 4.3 and 4.4 to illustrate the effect of permeability; not fitted to data.
  • M (mobility) = 1 (accuracy test), 0.03 (Examples 4.3-4.4)
    Selected as test values; not fitted.
  • beta (osmotic coupling) = 1 (accuracy test), 2*epsilon (Examples 4.3-4.4)
    Chosen to define the test cases; in Examples 4.3-4.4 beta is set to 2*epsilon, giving a specific osmotic strength.
  • epsilon (interface thickness) = 1 (accuracy test), 0.03 (Examples 4.3-4.4)
    Chosen as the diffuse interface width; in the accuracy test epsilon=1 is not a sharp-interface limit but serves the manufactured solution test.
  • S (stabilization parameter) = S >= L/2
    Appears in the scheme (11) and (30); must be chosen large enough to satisfy the stability condition. Its value is not specified in the numerical tests, which is a reproducibility gap.
assumptions (5)
  • domain assumption The total energy functional (2) with kinetic, entropy, and mixing terms is the correct free energy for the system.
    This is the starting point of the energy variational approach; it encodes the physics of the system but is not derived.
  • ad hoc to paper The dissipative functional (3) with the interface permeability term K|grad phi| |mu|^2 is the correct dissipation.
    The specific form, especially the K term proportional to |grad phi|, is assumed to model transmembrane flux; it is the key constitutive assumption of the model.
  • ad hoc to paper Condition (16): there exists L such that max|F''| <= L, and the double-well potential is truncated to satisfy it.
    Required for the unconditional energy stability proofs (Theorems 3.1, 3.2). The untruncated F does not satisfy (16) globally.
  • domain assumption The membrane is impermeable to solute: there is no flux of C across the interface, and water flux follows the Allen-Cahn-type law S_phi = -K mu |grad phi|.
    The model assumes the solute is confined to each phase while water crosses the membrane. This is a strong biophysical assumption for some but not all semi-permeable membranes.
  • standard math Solutions are sufficiently smooth for integration by parts and the LDG numerical fluxes; boundary conditions (4) hold.
    Needed for the energy dissipation derivations in the appendix.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Phase-Field Modeling and Energy-Stable Schemes for Osmotic Flow through Semi-Permeable." pith.science (2026). https://pith.science/paper/5FIGPVWA

@misc{pith2026250611374,
  author       = {Pith},
  title        = {Pith review of: Phase-Field Modeling and Energy-Stable Schemes for Osmotic Flow through Semi-Permeable},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/5FIGPVWA}},
  note         = {Machine review of arXiv:2506.11374}
}
read the original abstract

We present a thermodynamically consistent phase-field model for simulating fluid transport across semi-permeable membranes, with a particular focus on osmotic pressure effects. The model extends the classical Navier-Stokes-Cahn-Hilliard (NSCH) system by introducing an Allen-Cahn-type transmembrane flux governed by chemical potential imbalances, resulting in a strongly coupled system involving fluid motion, solute transport, and interface dynamics. To solve this system efficiently and accurately, we develop high-order, energy-stable numerical schemes. The local discontinuous Galerkin (LDG) method is employed for spatial discretization, offering high-order accuracy and geometric flexibility. For temporal integration, we first construct a first-order decoupled scheme with rigorous energy stability, and then improve temporal accuracy via a semi-implicit spectral deferred correction (SDC) method. Numerical experiments confirm the theoretical properties of the proposed scheme and demonstrate the influence of osmotic pressure and membrane permeability on droplet morphology at equilibrium. The framework offers a robust and versatile tool for modeling transmembrane fluid transport in both biological and industrial applications.

Figures

Figures reproduced from arXiv: 2506.11374 by the authors.

Figure 1
Figure 1. Energy evolution of the first order scheme and the second order SDC scheme. [PITH_FULL_IMAGE:figures/full_fig_p018_1.png] view at source ↗
Figure 2
Figure 2. Initial condition for Example 4.3 with c0 = 0.15. profiles of concentration and interface are given as follows (see [PITH_FULL_IMAGE:figures/full_fig_p019_2.png] view at source ↗
Figure 3
Figure 3. Time evolution of the free energy functional with different permeability: [PITH_FULL_IMAGE:figures/full_fig_p020_3.png] view at source ↗
Figures from the paper (6 more)
Figure 4
Figure 4. Figure 4: Time evolution of the mass (R Ω ζ+C+dx) with different permeability: K = 0, 0.1 and 0.3 with c0 = 0.15. 20 [PITH_FULL_IMAGE:figures/full_fig_p020_4.png]
Figure 5
Figure 5. Figure 5: Snapshots of interface ϕ = 0 at different time slots with different permeability: K = 0, 0.1 and 0.3 with c0 = 0.15. 1 2 3 4 5 6 1 2 3 4 5 0.3 0.275 0.25 0.225 0.2 0.175 0.15 0.125 0.1 0.075 0.05 (a) K = 0 1 2 3 4 5 6 1 2 3 4 5 0.5 0.455 0.41 0.365 0.32 0.275 0.23 0.18…
Figure 6
Figure 6. Figure 6: Snapshots of ζ+C+ at equilibrium with different permeability: K = 0, 0.1 and 0.3 with c0 = 0.15. equilibrium, the concentrations C ± inside and outside the droplet become spatially uniform, with their values determined by the total initial solute mass and the geometry …
Figure 7
Figure 7. Figure 7: The profiles of the interface when the osmotic pressure is neglected at time slots [PITH_FULL_IMAGE:figures/full_fig_p022_7.png]
Figure 8
Figure 8. Figure 8: The interface profiles of the droplet at equilibrium for different initial concentrations [PITH_FULL_IMAGE:figures/full_fig_p022_8.png]
Figure 9
Figure 9. Figure 9: Snapshots of ϕ with K = 0 (first row) and K = 0.3 (second row) at different time t = 1.5 (first column), t = 6 (second column), t = 15 (third column) and t = 25 (fourth column). C+(x, y, 0) = 0.1 + exp − (x − 2.35)2 + (y − 4.05)2 0.15 ! + exp − (x − 4.05)2 + (y − 2.35)…

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

34 extracted references · 34 canonical work pages

  1. [1]

    An immersed membrane method for mass transfer across flex- ible semipermeable membranes in flows

    Farhad A Amiri and Junfeng Zhang. An immersed membrane method for mass transfer across flex- ible semipermeable membranes in flows. International Communications in Heat and Mass Transfer, 128:105601, 2021

    work page 2021

  2. [2]

    Diffuse-interface methods in fluid mechanics

    Donald M Anderson, Geoffrey B McFadden, and A A Wheeler. Diffuse-interface methods in fluid mechanics. Annual Review of Fluid Mechanics, 30(1):139–165, 1998

    work page 1998

  3. [3]

    AnL∞ bound for solutions of the Cahn-Hilliard equation.Archive for Rational Mechanics and Analysis, 133(2):129–144, 1995

    Luis A Caffarelli and Nora E Muler. AnL∞ bound for solutions of the Cahn-Hilliard equation.Archive for Rational Mechanics and Analysis, 133(2):129–144, 1995

    work page 1995

  4. [4]

    Periodic migration in a physical model of cells on micropatterns.Physical Review Letters, 111(15):158102, 2013

    Brian A Camley, Yanxiang Zhao, Bo Li, Herbert Levine, and Wouter-Jan Rappel. Periodic migration in a physical model of cells on micropatterns.Physical Review Letters, 111(15):158102, 2013. 24

    work page 2013

  5. [5]

    Decoupled energy stable schemes for phase-field vesicle membrane model.Journal of Computational Physics, 302:509–523, 2015

    Rui Chen, Guanghua Ji, Xiaofeng Yang, and Hui Zhang. Decoupled energy stable schemes for phase-field vesicle membrane model.Journal of Computational Physics, 302:509–523, 2015

    work page 2015

  6. [6]

    Convergence analysis of a fully discrete finite difference scheme for the Cahn-Hilliard-Hele-Shaw equation

    Wenbin Chen, Yuan Liu, Cheng Wang, and Steven Wise. Convergence analysis of a fully discrete finite difference scheme for the Cahn-Hilliard-Hele-Shaw equation. Mathematics of Computation, 85(301):2231–2257, 2016

    work page 2016

  7. [7]

    An energy stable fourth order finite difference scheme for the Cahn–Hilliard equation.Journal of Computational and Applied Mathematics, 362:574–595, 2019

    Kelong Cheng, Wenqiang Feng, Cheng Wang, and Steven M Wise. An energy stable fourth order finite difference scheme for the Cahn–Hilliard equation.Journal of Computational and Applied Mathematics, 362:574–595, 2019

    work page 2019

  8. [8]

    Water homeostasis and cell volume maintenance and regulation

    Eric Delpire and Kenneth B Gagnon. Water homeostasis and cell volume maintenance and regulation. Current Topics in Membranes, 81:3–52, 2018

    work page 2018

Show all 34 references
  1. [9]

    A finite element method for a phase field model of nematic liquid crystal droplets.Communications in Computational Physics, 25:155–188, 2019

    Amanda E Diegel and Shawn W Walker. A finite element method for a phase field model of nematic liquid crystal droplets.Communications in Computational Physics, 25:155–188, 2019

    2019

  2. [10]

    Energetic variational approaches in modeling vesicle and fluid interactions.Physica D: Nonlinear Phenomena, 238(9-10):923–930, 2009

    Qiang Du, Chun Liu, Rolf Ryham, and Xiaoqiang Wang. Energetic variational approaches in modeling vesicle and fluid interactions.Physica D: Nonlinear Phenomena, 238(9-10):923–930, 2009

    2009

  3. [11]

    A phase field approach in the numerical study of the elastic bending energy for vesicle membranes.Journal of Computational Physics, 198(2):450–468, 2004

    Qiang Du, Chun Liu, and Xiaoqiang Wang. A phase field approach in the numerical study of the elastic bending energy for vesicle membranes.Journal of Computational Physics, 198(2):450–468, 2004

    2004

  4. [12]

    Energy variational analysis of ions in water and channels: Field theory for primitive models of complex ionic fluids.The Journal of Chemical Physics, 133(10):104104, 2010

    Bob Eisenberg, Yunkyong Hyon, and Chun Liu. Energy variational analysis of ions in water and channels: Field theory for primitive models of complex ionic fluids.The Journal of Chemical Physics, 133(10):104104, 2010

    2010

  5. [13]

    Diffuse interface modelling of soluble surfactants in two-phase flow.Communications in Mathematical Sciences, 12:1475–1522, 2013

    Harald Garcke, Kei Fong Lam, and Björn Stinner. Diffuse interface modelling of soluble surfactants in two-phase flow.Communications in Mathematical Sciences, 12:1475–1522, 2013

    2013

  6. [14]

    An immersed boundary method for mass transfer across permeable moving interfaces.Journal of Computational Physics, 278:148–168, 2014

    Xiaobo Gong, Zhaoxin Gong, and Huaxiong Huang. An immersed boundary method for mass transfer across permeable moving interfaces.Journal of Computational Physics, 278:148–168, 2014

    2014

  7. [15]

    Ruihan Guo and Yinhua Xia. Arbitrary high-order fully-decoupled numerical schemes for phase-field models of two-phase incompressible flows.Communications on Applied Mathematics and Computation, 6:625–657, 2024

    2024

  8. [16]

    Ruihan Guo and Yan Xu. High order numerical simulations for the binary fluid-surfactant system using the discontinuous Galerkin and spectral deferred correction methods.SIAM Journal on Scientific Computing, 42:B353–B378, 2020

    2020

  9. [17]

    Semi-implicit spectral deferred correction methods for highly nonlinear partial differential equations.Journal of Computational Physics, 338:269–284, 2017

    Ruihan Guo, Yan Xu, and Yinhua Xia. Semi-implicit spectral deferred correction methods for highly nonlinear partial differential equations.Journal of Computational Physics, 338:269–284, 2017

    2017

  10. [18]

    Zhenlin Guo, Ping Lin, and John S Lowengrub. A numerical method for the quasi-incompressible Cahn–Hilliard–Navier–Stokes equations for variable density flows with a discrete energy law.Journal of Computational Physics, 276:486–507, 2014

    2014

  11. [19]

    Intracellular water homeostasis and the mammalian cellular osmotic stress response

    Steffan N Ho. Intracellular water homeostasis and the mammalian cellular osmotic stress response. Journal of Cellular Physiology, 206(1):9–15, 2006

    2006

  12. [20]

    Cellular pressure and volume regulation and implications for cell mechanics

    Hongyuan Jiang and Sean X Sun. Cellular pressure and volume regulation and implications for cell mechanics. Biophysical Journal, 105(3):609–619, 2013

    2013

  13. [21]

    Modeling water transport across elastic boundaries using an explicit jump method

    Anita T Layton. Modeling water transport across elastic boundaries using an explicit jump method. SIAM Journal on Scientific Computing, 28(6):2189–2207, 2006. 25

    2006

  14. [22]

    Lowengrub and L

    J. Lowengrub and L. Truskinovsky. Quasi-incompressible Cah-Hilliard fluids and topological transitions. Proceedings of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences, 454(1978):2617–2654, October 1998

    1978

  15. [23]

    John Lowengrub, Jun Allard, and Sebastian Aland. Numerical simulation of endocytosis: Viscous flow driven by membranes with non-uniformly distributed curvature-inducing molecules.Journal of Computational physics, 309:112–128, 2016

    2016

  16. [24]

    A phase field model for mass transport with semi-permeable interfaces.Journal of Computational Physics, 464:111334, 2022

    Yuzhe Qin, Huaxiong Huang, Yi Zhu, Chun Liu, and Shixin Xu. A phase field model for mass transport with semi-permeable interfaces.Journal of Computational Physics, 464:111334, 2022

    2022

  17. [25]

    Diffuse-interface approximations of osmosis free boundary problems.SIAM Journal on Applied Mathematics, 76(3):910–929, 2016

    Andreas Rätz. Diffuse-interface approximations of osmosis free boundary problems.SIAM Journal on Applied Mathematics, 76(3):910–929, 2016

    2016

  18. [26]

    An energy stableC 0 finite element scheme for a quasi-incompressible phase-field model of moving contact line with variable density

    Lingyue Shen, Huaxiong Huang, Ping Lin, Zilong Song, and Shixin Xu. An energy stableC 0 finite element scheme for a quasi-incompressible phase-field model of moving contact line with variable density. Journal of Computational Physics, 405:109179, 2020

    2020

  19. [27]

    An energy stableC 0 finite element scheme for a phase-field model of vesicle motion and deformation.SIAM Journal on Scientific Computing, 44(1):B122–B145, 2022

    Lingyue Shen, Zhiliang Xu, Ping Lin, Huaxiong Huang, and Shixin Xu. An energy stableC 0 finite element scheme for a phase-field model of vesicle motion and deformation.SIAM Journal on Scientific Computing, 44(1):B122–B145, 2022

    2022

  20. [28]

    Phase field modeling and compu- tation of vesicle growth or shrinkage.Journal of Mathematical Biology, 86(6):97, 2023

    Xiaoxia Tang, Shuwang Li, John S Lowengrub, and Steven M Wise. Phase field modeling and compu- tation of vesicle growth or shrinkage.Journal of Mathematical Biology, 86(6):97, 2023

    2023

  21. [29]

    An immersed boundary method for mass transfer through porous biomembranes under large deformations.Journal of Computational Physics, 413:109444, 2020

    Xiaolong Wang, Xiaobo Gong, Kazuyasu Sugiyama, Shu Takagi, and Huaxiong Huang. An immersed boundary method for mass transfer through porous biomembranes under large deformations.Journal of Computational Physics, 413:109444, 2020

    2020

  22. [30]

    LocaldiscontinuousGalerkinmethodsforhigh-ordertime-dependentpartial differential equations

    YanXuandChi-WangShu. LocaldiscontinuousGalerkinmethodsforhigh-ordertime-dependentpartial differential equations. Comminations in Computational Physics, 7:1–46, 2010

    2010

  23. [31]

    A second-order energy stable BDF numerical schemefortheCahn-Hilliardequation

    Yue Yan, Wenbin Chen, Cheng Wang, and Steven M Wise. A second-order energy stable BDF numerical schemefortheCahn-Hilliardequation. Communications in Computational Physics, 23(2):572–602, 2018

    2018

  24. [32]

    Convergence analysis for the Invariant Energy Quadratization (IEQ) schemes for solving the Cahn–Hilliard and Allen–Cahn equations with general nonlinear potential

    Xiaofeng Yang and Guo-Dong Zhang. Convergence analysis for the Invariant Energy Quadratization (IEQ) schemes for solving the Cahn–Hilliard and Allen–Cahn equations with general nonlinear potential. Journal of Scientific Computing, 82(3):1–28, 2020

    2020

  25. [33]

    Numerical approximations for the molecular beam epitaxial growth model based on the invariant energy quadratization method.Journal of Computational Physics, 333:104–127, 2017

    Xiaofeng Yang, Jia Zhao, and Qi Wang. Numerical approximations for the molecular beam epitaxial growth model based on the invariant energy quadratization method.Journal of Computational Physics, 333:104–127, 2017

    2017

  26. [34]

    A numerical method for osmotic water flow and solute diffusion with deformable membrane boundaries in two spatial dimension.Journal of Computational Physics, 350:728–746, 2017

    Lingxing Yao and Yoichiro Mori. A numerical method for osmotic water flow and solute diffusion with deformable membrane boundaries in two spatial dimension.Journal of Computational Physics, 350:728–746, 2017. 26 A Model Derivation Details We calculate the time derivative of th...

    2017

Pith tools

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