BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

John Lowengrub

John Lowengrub appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.NANumerical Analysiscond-mat.softphysics.flu-dynQuantitative MethodsBiological PhysicsCell Behaviorcond-mat.stat-mechMachine Learningmath-phmath.APmath.MPnlin.PSphysics.comp-phPopulations and Evolution

Catalog footprint

What is connected

15works
15topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

15 published item(s)

preprint2026arXiv

An Efficient Constant-Coefficient MSAV Scheme for Computing Vesicle Growth and Shrinkage

We present a fast, unconditionally energy-stable numerical scheme for simulating vesicle deformation under osmotic pressure using a phase-field approach. The model couples an Allen-Cahn equation for the biomembrane interface with a variable-mobility Cahn-Hilliard equation governing mass exchange across the membrane. Classical approaches, including nonlinear multigrid and Multiple Scalar Auxiliary Variable (MSAV) methods, require iterative solution of variable-coefficient systems at each time step, resulting in substantial computational cost. We introduce a constant-coefficient MSAV (CC-MSAV) scheme that incorporates stabilization directly into the Cahn-Hilliard evolution equation rather than the chemical potential. This reformulation yields fully decoupled constant-coefficient elliptic problems solvable via fast discrete cosine transform (DCT), eliminating iterative solvers entirely. The method achieves O(N^2 log N) complexity per time step while preserving unconditional energy stability and discrete mass conservation. Numerical experiments verify second-order temporal and spatial accuracy, mass conservation to relative errors below 5 x 10^-11, and close agreement with nonlinear multigrid benchmarks. On grids with N >= 2048, CC-MSAV achieves 6-15x overall speedup compared to classical MSAV with optimized preconditioning, while the dominant Cahn-Hilliard subsystem is accelerated by up to two orders of magnitude. These efficiency gains, achieved without sacrificing accuracy, make CC-MSAV particularly well suited for large-scale simulations of vesicle dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Nonlinear simulation of vascular tumor growth with chemotaxis and the control of necrosis

In this paper, we develop a sharp interface tumor growth model in two dimensions to study the effect of both the intratumoral structure using a controlled necrotic core and the extratumoral nutrient supply from vasculature on tumor morphology. We first show that our model extends the benchmark results in the literature using linear stability analysis. Then we solve this generalized model numerically using a spectrally accurate boundary integral method in an evolving annular domain, not only with a Robin boundary condition on the outer boundary for the nutrient field which models tumor vasculature, but also with a static boundary condition on the inner boundary for pressure field which models the control of tumor necrosis. The discretized linear systems for both pressure and nutrient fields are shown to be well-conditioned through tracing GMRES iteration numbers. Our nonlinear simulations reveal the stabilizing effects of angiogenesis and the destabilizing ones of chemotaxis and necrosis in the development of tumor morphological instabilities if the necrotic core is fixed in a circular shape. When the necrotic core is controlled in a non-circular shape, the stabilizing effects of proliferation and the destabilizing ones of apoptosis are observed. Finally, the values of the nutrient concentration with its fluxes and the pressure level with its normal derivatives, which are solved accurately at the boundaries, help us to characterize the corresponding tumor morphology and the level of the biophysical quantities on interfaces required in keeping various shapes of the necrotic region of the tumor. Interestingly, we notice that when the necrotic region is fixed in a 3-fold non-circular shape, even if the initial shape of the tumor is circular, the tumor will evolve into a shape corresponding to the 3-fold symmetry of the shape of the fixed necrotic region.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A multi-scale model of virus pandemic: Heterogeneous interactive entities in a globally connected world

This paper is devoted to the multidisciplinary modelling of a pandemic initiated by an aggressive virus, specifically the so-called \textit{SARS--CoV--2 Severe Acute Respiratory Syndrome, corona virus n.2}. The study is developed within a multiscale framework accounting for the interaction of different spatial scales, from the small scale of the virus itself and cells, to the large scale of individuals and further up to the collective behaviour of populations. An interdisciplinary vision is developed thanks to the contributions of epidemiologists, immunologists and economists as well as those of mathematical modellers. The first part of the contents is devoted to understanding the complex features of the system and to the design of a modelling rationale. The modelling approach is treated in the second part of the paper by showing both how the virus propagates into infected individuals, successfully and not successfully recovered, and also the spatial patterns, which are subsequently studied by kinetic and lattice models. The third part reports the contribution of research in the fields of virology, epidemiology, immune competition, and economy focused also on social behaviours. Finally, a critical analysis is proposed looking ahead to research perspectives.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Complex Far-Field Geometries Determine the Stability of Solid Tumor Growth with Chemotaxis

In this paper, we develop a sharp interface tumor growth model to study the effect of the tumor microenvironment using a complex far-field geometry that mimics a heterogeneous distribution of vasculature. Together with different nutrient uptake rates inside and outside the tumor, this introduces variability in spatial diffusion gradients. Linear stability analysis suggests that the uptake rate in the tumor microenvironment, together with chemotaxis, may induce unstable growth, especially when the nutrient gradients are large. We investigate the fully nonlinear dynamics using a spectrally accurate boundary integral method. Our nonlinear simulations reveal that vascular heterogeneity plays an important role in the development of morphological instabilities that range from fingering and chain-like morphologies to compact, plate-like shapes in two-dimensions.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Nonlinear limiting dynamics of a shrinking interface in a Hele-Shaw cell

The flow in a Hele-Shaw cell with a time-increasing gap poses a unique shrinking interface problem. When the upper plate of the cell is lifted perpendicularly at a prescribed speed, the exterior less viscous fluid penetrates the interior more viscous fluid, which generates complex, time-dependent interfacial patterns through the Saffman-Taylor instability. The pattern formation process sensitively depends on the lifting speed and is still not fully understood. For some lifting speeds, such as linear or exponential speed, the instability is transient and the interface eventually shrinks as a circle. However, linear stability analysis suggests there exist shape invariant shrinking patterns if the gap $b(t)$ is increased more rapidly: $\displaystyle b(t)=\left(1-\frac{7}{2}τ\mathcal{C} t\right)^{-{2}/{7}}$, where $τ$ is the surface tension and $\mathcal{C}$ is a function of the interface perturbation mode $k$. Here, we use a spectrally accurate boundary integral method together with an efficient time adaptive rescaling scheme, which for the first time makes it possible to explore the nonlinear limiting dynamical behavior of a vanishing interface. When the gap is increased at a constant rate, our numerical results quantitatively agree with experimental observations (Nase et al., Phys. Fluids, vol. 23, 2011, pp. 123101). When we use the shape invariant gap $b(t)$, our nonlinear results reveal the existence of $k$-fold dominant, one-dimensional, web-like networks, where the fractal dimension is reduced to almost one at late times. We conclude by constructing a morphology diagram for pattern selection that relates the dominant mode $k$ of the vanishing interface and the control parameter $\mathcal{C}$.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

An Energy Stable Finite-Difference Scheme for Functionalized Cahn-Hilliard Equation and its Convergence Analysis

We present and analyze an unconditionally energy stable and convergent finite difference scheme for the Functionalized Cahn-Hilliard equation. One key difficulty associated with the energy stability is based on the fact that one nonlinear energy functional term in the expansion appears as non-convex, non-concave. To overcome this subtle difficulty, we add two auxiliary terms to make the combined term convex, which in turns yields a convex-concave decomposition of the physical energy. As a result, an application of the convex splitting methodology assures both the unique solvability and the unconditional energy stability of the proposed numerical scheme. To deal with a 4-Laplacian solver in an $H^{-1}$ gradient flow at each time step, we apply an efficient preconditioned steepest descent algorithm to solve the corresponding nonlinear systems. In addition, a global in time $H_{\rm per}^2$ stability of the numerical scheme is established at a theoretical level, which in turn ensures the full order convergence analysis of the scheme. A few numerical results are presented, which confirm the stability and accuracy of the proposed numerical scheme.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Dynamics of a multicomponent vesicle in shear flow

We study the fully nonlinear, nonlocal dynamics of two-dimensional multicomponent vesicles in a shear flow with matched viscosity of the inner and outer fluids. Using a nonstiff, pseudo-spectral boundary integral method, we investigate dynamical patterns induced by inhomogeneous bending for a two phase system. Numerical results reveal that there exist novel phase-treading and tumbling mechanisms that cannot be observed for a homogeneous vesicle. In particular, unlike the well-known steady tank-treading dynamics characterized by a fixed inclination angle, here the phase-treading mechanism leads to unsteady periodic dynamics with an oscillatory inclination angle. When the average phase concentration is around 1/2, we observe tumbling dynamics even for very low shear rate, and the excess length required for tumbling is significantly smaller than the value for the single phase case. We summarize our results in phase diagrams in terms of the excess length, shear rate, and concentration of the soft phase. These findings go beyond the well known dynamical regimes of a homogeneous vesicle and highlight the level of complexity of vesicle dynamics in a fluid due to heterogeneous material properties.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Analysis of the diffuse-domain method for solving PDEs in complex geometries

In recent work, Li et al.\ (Comm.\ Math.\ Sci., 7:81-107, 2009) developed a diffuse-domain method (DDM) for solving partial differential equations in complex, dynamic geometries with Dirichlet, Neumann, and Robin boundary conditions. The diffuse-domain method uses an implicit representation of the geometry where the sharp boundary is replaced by a diffuse layer with thickness $ε$ that is typically proportional to the minimum grid size. The original equations are reformulated on a larger regular domain and the boundary conditions are incorporated via singular source terms. The resulting equations can be solved with standard finite difference and finite element software packages. Here, we present a matched asymptotic analysis of general diffuse-domain methods for Neumann and Robin boundary conditions. Our analysis shows that for certain choices of the boundary condition approximations, the DDM is second-order accurate in $ε$. However, for other choices the DDM is only first-order accurate. This helps to explain why the choice of boundary-condition approximation is important for rapid global convergence and high accuracy. Our analysis also suggests correction terms that may be added to yield more accurate diffuse-domain methods. Simple modifications of first-order boundary condition approximations are proposed to achieve asymptotically second-order accurate schemes. Our analytic results are confirmed numerically in the $L^2$ and $L^\infty$ norms for selected test problems.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Collective Properties of a Transcription Initiation Model under Varying Environment

The dynamics of gene transcription is tightly regulated in eukaryotes. Recent experiments have revealed various kinds of transcriptional dynamics, such as RNA polymerase II pausing, that involves regulation at the transcription initiation stage, and the choice of different regulation pattern is closely related to the physiological functions of the target gene. Here we consider a simplified model of transcription initiation, a process including the assembly of transcription complex and the pausing and releasing of the RNA polymerase II. Focusing on the collective behaviors on a population level, we explore potential regulatory functions this model can offer. These functions include fast and synchronized response to environmental change, or long-term memory about the transcriptional status. As a proof of concept we also show that, by selecting different control mechanisms cells can adapt to different environments. These findings may help us better understand the design principles of transcriptional regulation.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Energy Stable Multigrid Method for Local and Non-local Hydrodynamic Models for Freezing

In this paper we present a numerical method for hydrodynamic models that arise from time dependent density functional theories of freezing. The models take the form of compressible Navier-Stokes equations whose pressure is determined by the variational derivative of a free energy, which is a functional of the density field. We present unconditionally energy stable and mass conserving implicit finite difference methods for the models. The methods are based on a convex splitting of the free energy and that ensures that a discrete energy is non-increasing for any choice of time and space step. The methods are applicable to a large class of models, including both local and non-local free energy functionals. The theoretical basis for the numerical method is presented in a general context. The method is applied to problems using two specific free energy functionals: one local and one non-local functional. A nonlinear multigrid method is used to solve the numerical method, which is nonlinear at the implicit time step. The non-local functional, which is a convolution operator, is approximated using the Discrete Fourier Transform. Numerical simulations that confirm the stability and accuracy of the numerical method are presented.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Kinetic Density Functional Theory of Freezing

A theory of freezing of a dense hard sphere gas is presented. Starting from a revised Enskog theory, hydrodynamic equations that account for non-local variations in the density but local variations in the flow field are derived using a modified Chapman Enskog procedure. These hydrodynamic equations, which retain structural correlations, are shown to be effectively a time dependent density functional theory. The ability of this theory to capture the solid liquid phase transition is established through analysis and numerical simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

POPE: Post Optimization Posterior Evaluation of Likelihood Free Models

In many domains, scientists build complex simulators of natural phenomena that encode their hypotheses about the underlying processes. These simulators can be deterministic or stochastic, fast or slow, constrained or unconstrained, and so on. Optimizing the simulators with respect to a set of parameter values is common practice, resulting in a single parameter setting that minimizes an objective subject to constraints. We propose a post optimization posterior analysis that computes and visualizes all the models that can generate equally good or better simulation results, subject to constraints. These optimization posteriors are desirable for a number of reasons among which easy interpretability, automatic parameter sensitivity and correlation analysis and posterior predictive analysis. We develop a new sampling framework based on approximate Bayesian computation (ABC) with one-sided kernels. In collaboration with two groups of scientists we applied POPE to two important biological simulators: a fast and stochastic simulator of stem-cell cycling and a slow and deterministic simulator of tumor growth patterns.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Diffuse interface models of locally inextensible vesicles in a viscous fluid

We present a new diffuse interface model for the dynamics of inextensible vesicles in a viscous fluid. A new feature of this work is the implementation of the local inextensibility condition in the diffuse interface context. Local inextensibility is enforced by using a local Lagrange multiplier, which provides the necessary tension force at the interface. To solve for the local Lagrange multiplier, we introduce a new equation whose solution essentially provides a harmonic extension of the local Lagrange multiplier off the interface while maintaining the local inextensibility constraint near the interface. To make the method more robust, we develop a local relaxation scheme that dynamically corrects local stretching/compression errors thereby preventing their accumulation. Asymptotic analysis is presented that shows that our new system converges to a relaxed version of the inextensible sharp interface model. This is also verified numerically. Although the model does not depend on dimension, we present numerical simulations only in 2D. To solve the 2D equations numerically, we develop an efficient algorithm combining an operator splitting approach with adaptive finite elements where the Navier-Stokes equations are implicitly coupled to the diffuse interface inextensibility equation. Numerical simulations of a single vesicle in a shear flow at different Reynolds numbers demonstrate that errors in enforcing local inextensibility may accumulate and lead to large differences in the dynamics in the tumbling regime and differences in the inclination angle of vesicles in the tank-treading regime. The local relaxation algorithm is shown to effectively prevent this accumulation by driving the system back to its equilibrium state when errors in local inextensibility arise.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Analysis of a Mixture Model of Tumor Growth

We study an initial-boundary value problem (IBVP) for a coupled Cahn-Hilliard-Hele-Shaw system that models tumor growth. For large initial data with finite energy, we prove global (local resp.) existence, uniqueness, higher order spatial regularity and Gevrey spatial regularity of strong solutions to the IBVP in 2D (3D resp.). Asymptotically in time, we show that the solution converges to a constant state exponentially fast as time tends to infinity under certain assumptions.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Convergence Analysis of a Second Order Convex Splitting Scheme for the Modified Phase Field Crystal Equation

In this paper we provide a detailed convergence analysis for an unconditionally energy stable, second-order accurate convex splitting scheme for the Modified Phase Field Crystal equation, a generalized damped wave equation for which the usual Phase Field Crystal equation is a special degenerate case. The fully discrete, fully second-order finite difference scheme in question was derived in a recent work [2]. An introduction of a new variable ψ, corresponding to the temporal derivative of the phase variable ϕ, could bring an accuracy reduction in the formal consistency estimate, because of the hyperbolic nature of the equation. A higher order consistency analysis by an asymptotic expansion is performed to overcome this difficulty. In turn, second order convergence in both time and space is established in a discrete L^\infty (0,T; H^3) norm.

0 citations0 reviewsNo code linkedOpen access