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

David M. Bortz

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

ResearcherUnclaimed source record
math.NAmath.APmath.DSmath.OCPopulations and Evolutioncond-mat.softMachine LearningNumerical AnalysisQuantitative MethodsBiological PhysicsBiomoleculesCell BehaviorComputational Engineering, Finance, and Sciencemath.CAmath.FAmath.PRmath.STphysics.flu-dynStatistics TheorySubcellular ProcessesTissues and Organs

Catalog footprint

What is connected

17works
21topics
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

17 published item(s)

preprint2022arXiv

Asymptotic Analysis of a General Multi-Structured Population Model

Structured populations are ubiquitous across the biological sciences. Mathematical models of these populations allow us to understand how individual physiological traits drive the overall dynamics in aggregate. For example, linear age- or age-and-size-structured models establish constraints on individual growth under which the age- or age-and-size- distribution stabilizes, even as the population continues to grow without bound. However, individuals in real-world populations exhibit far more structural features than simply age and size. Notably, cyanobacteria contain carboxysome organelles which are central to carbon fixation and can be older (if inherited from parent cells) or younger (if created after division) than the enveloping cell. Motivated by a desire to understand how carboxysome age impacts growth at the colony level, we develop a multi-structured model which allows for an arbitrary (but finite) number of structure variables. We then derive and solve the renewal equation for cell division to obtain an asymptotic solution, and show that, under certain conditions, a stable structural distribution is reached.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Learning diffusion coefficients, kinetic parameters, and the number of underlying states from a multi-state diffusion process: robustness results and application to PDK1/PKC$α$, dynamics

Systems driven by Brownian motion are ubiquitous. A prevailing challenge is inferring, from data, the diffusion and kinetic parameters that describe these stochastic processes. In this work, we investigate a multi-state diffusion process that arises in the context of single particle tracking (SPT), wherein the motion of a particle is governed by a discrete set of diffusive states, and the tendency of the particle to switch between these states is modeled as a random process. We consider two models for this behavior: a mixture model and a hidden Markov model (HMM). For both, we adopt a Bayesian approach to sample the distributions of the underlying parameters and implement a Markov Chain Monte Carlo (MCMC) scheme to compute the posterior distributions, as in Das, Cairo, Coombs (2009). The primary contribution of this work is a study of the robustness of this method to infer parameters of a three-state HMM, and a discussion of the challenges and degeneracies that arise from considering three states. Finally, we investigate the problem of determining the number of diffusive states using model selection criteria. We present results from simulated data that demonstrate proof of concept, as well as apply our method to experimentally measured single molecule diffusion trajectories of monomeric phosphoinositide-dependent kinase-1 (PDK1) on a synthetic target membrane where it can associate with its binding partner protein kinase C alpha isoform (PKC$α$) to form a heterodimer detected by its significantly lower diffusivity. All matlab software is available here: \url{https://github.com/MathBioCU/SingleMolecule}

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Online Weak-form Sparse Identification of Partial Differential Equations

This paper presents an online algorithm for identification of partial differential equations (PDEs) based on the weak-form sparse identification of nonlinear dynamics algorithm (WSINDy). The algorithm is online in a sense that if performs the identification task by processing solution snapshots that arrive sequentially. The core of the method combines a weak-form discretization of candidate PDEs with an online proximal gradient descent approach to the sparse regression problem. In particular, we do not regularize the $\ell_0$-pseudo-norm, instead finding that directly applying its proximal operator (which corresponds to a hard thresholding) leads to efficient online system identification from noisy data. We demonstrate the success of the method on the Kuramoto-Sivashinsky equation, the nonlinear wave equation with time-varying wavespeed, and the linear wave equation, in one, two, and three spatial dimensions, respectively. In particular, our examples show that the method is capable of identifying and tracking systems with coefficients that vary abruptly in time, and offers a streaming alternative to problems in higher dimensions.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Learning Mean-Field Equations from Particle Data Using WSINDy

We develop a weak-form sparse identification method for interacting particle systems (IPS) with the primary goals of reducing computational complexity for large particle number $N$ and offering robustness to either intrinsic or extrinsic noise. In particular, we use concepts from mean-field theory of IPS in combination with the weak-form sparse identification of nonlinear dynamics algorithm (WSINDy) to provide a fast and reliable system identification scheme for recovering the governing stochastic differential equations for an IPS when the number of particles per experiment $N$ is on the order of several thousand and the number of experiments $M$ is less than 100. This is in contrast to existing work showing that system identification for $N$ less than 100 and $M$ on the order of several thousand is feasible using strong-form methods. We prove that under some standard regularity assumptions the scheme converges with rate $\mathcal{O}(N^{-1/2})$ in the ordinary least squares setting and we demonstrate the convergence rate numerically on several systems in one and two spatial dimensions. Our examples include a canonical problem from homogenization theory (as a first step towards learning coarse-grained models), the dynamics of an attractive-repulsive swarm, and the IPS description of the parabolic-elliptic Keller-Segel model for chemotaxis.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Competitive Exclusion in a DAE Model for Microbial Electrolysis Cells

Microbial electrolysis cells (MECs) employ electroactive bacteria to perform extracellular electron transfer, enabling hydrogen generation from biodegradable substrates. In previous work, we developed and analyzed a differential-algebraic equation (DAE) model for MECs. The model resembles a chemostat with ordinary differential equations (ODEs) for concentrations of substrate, microorganisms, and an extracellular mediator involved in electron transfer. There is also an algebraic constraint for electric current and hydrogen production. Our goal is to determine the outcome of competition between methanogenic archaea and electroactive bacteria, because only the latter contribute to electric current and resulting hydrogen production. We investigate asymptotic stability in two industrially relevant versions of the model. An important aspect of chemostats models is the principle of competitive exclusion -- only microbes which grow at the lowest substrate concentration will survive as $t\to\infty$. We show that if methanogens grow at the lowest substrate concentration, then the equilibrium corresponding to competitive exclusion by methanogens is globally asymptotically stable. The analogous result for electroactive bacteria is not necessarily true. We show that local asymptotic stability of exclusion by electroactive bacteria is not guaranteed, even in a simplified version of the model. In this case, even if electroactive bacteria can grow at the lowest substrate concentration, a few additional conditions are required to guarantee local asymptotic stability. We also provide numerical simulations supporting these arguments. Our results suggest operating conditions that are most conducive to success of electroactive bacteria and the resulting current and hydrogen production in MECs. This will help identify when methane production or electricity and hydrogen production are favored.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Weak SINDy: Galerkin-Based Data-Driven Model Selection

We present a weak formulation and discretization of the system discovery problem from noisy measurement data. This method of learning differential equations from data fits into a new class of algorithms that replace pointwise derivative approximations with linear transformations and a variance reduction technique. Our approach improves on the standard SINDy algorithm by orders of magnitude. We first show that in the noise-free regime, this so-called Weak SINDy (WSINDy) framework is capable of recovering the dynamic coefficients to very high accuracy, with the number of significant digits equal to the tolerance of the data simulation scheme. Next we show that the weak form naturally accounts for white noise by identifying the correct nonlinearities with coefficient error scaling favorably with the signal-to-noise ratio while significantly reducing the size of linear systems in the algorithm. In doing so, we combine the ease of implementation of the SINDy algorithm with the natural noise-reduction of integration to arrive at a more robust and user-friendly method of sparse recovery that correctly identifies systems in both small-noise and large-noise regimes.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A numerical framework for computing steady states of size-structured population models and their stability

Structured population models are a class of general evolution equations which are widely used in the study of biological systems. Many theoretical methods are available for establishing existence and stability of steady states of general evolution equations. However, except for very special cases, finding an analytical form of stationary solutions for evolution equations is a challenging task. In the present paper, we develop a numerical framework for computing approximations to stationary solutions of general evolution equations, which can also be used to produce existence and stability regions for steady states. In particular, we use the Trotter-Kato Theorem to approximate the infinitesimal generator of an evolution equation on a finite dimensional space, which in turn reduces the evolution equation into a system of ordinary differential equations. Consequently, we approximate and study the asymptotic behavior of stationary solutions. We illustrate the convergence of our numerical framework by applying it to a linear Sinko-Streifer structured population model for which the exact form of the steady state is known. To further illustrate the utility of our approach, we apply our framework to nonlinear population balance equation, which is an extension of well-known Smoluchowksi coagulation-fragmentation model to biological populations. We also demonstrate that our numerical framework can be used to gain insight about the theoretical stability of the stationary solutions of the evolution equations. Furthermore, the open source Python program that we have developed for our numerical simulations is freely available from our Github repository (github.com/MathBioCU).

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Criteria for linearized stability for a size-structured population model

We consider a size-structured aggregation and growth model of phytoplankton community proposed by Ackleh and Fitzpatrick [2]. The model accounts for basic biological phenomena in phytoplankton community such as growth, gravitational sedimentation, predation by zooplankton, fecundity, and aggregation. Our primary goal in this paper is to investigate the long-term behavior of the proposed aggregation and growth model. Particularly, using the well-known principle of linearized stability and semigroup compactness arguments, we provide sufficient conditions for local exponential asymptotic stability of zero solution as well as sufficient conditions for instability. We express these conditions in the form of an easy to compute characteristic function, which depends on the functional relationship between growth, sedimentation and fecundity. Our results can be used to predict long-term phytoplankton dynamic

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Inverse Problems for a Class of Conditional Probability Measure-Dependent Evolution Equations

We investigate the inverse problem of identifying a conditional probability measure in a measure-dependent dynamical system. We provide existence and well-posedness results and outline a discretization scheme for approximating a measure. For this scheme, we prove general method stability. The work is motivated by Partial Differential Equation (PDE) models of flocculation for which the shape of the post-fragmentation conditional probability measure greatly impacts the solution dynamics. To illustrate our methodology, we apply the theory to a particular PDE model that arises in the study of population dynamics for flocculating bacterial aggregates in suspension, and provide numerical evidence for the utility of the approach.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Stability of steady states for a class of flocculation equations with growth and removal

Flocculation is the process whereby particles (i.e., flocs) in suspension reversibly combine and separate. The process is widespread in soft matter and aerosol physics as well as environmental science and engineering. We consider a general size-structured flocculation model, which describes the evolution of flocs in an aqueous environment. Our work provides a unified treatment for many size-structured models in the environmental, industrial, medical, and marine engineering literature. In particular, our model accounts for basic biological phenomena in a population of microorganisms including growth, death, sedimentation, predation, renewal, fragmentation and aggregation. Our central goal in this paper is to rigorously investigate the long-term behavior of this generalized flocculation model. Using results from fixed point theory we derive conditions for the existence of continuous, non-trivial stationary solutions. We further apply the principle of linearized stability and semigroup compactness arguments to provide sufficient conditions for local exponential stability of stationary solutions as well as sufficient conditions for instability. Abstract. The end results of this analytical development are relatively simple inequality-criteria which thus allows for the rapid evaluation of the existence and stability of a non-trivial stationary solution. To our knowledge, this work is the first to derive precise existence and stability criteria for such a generalized model. Lastly, we also provide an illustrating application of this criteria to several flocculation models.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Variable Viscosity and Density Biofilm Simulations using an Immersed Boundary Method, Part II: Experimental Validation and the Heterogeneous Rheology-IBM

The goal of this work is to develop a numerical simulation that accurately captures the biomechanical response of bacterial biofilms and their associated extracellular matrix (ECM). In this, the second of a two-part effort, the primary focus is on formally presenting the heterogeneous rheology Immersed Boundary Method (hrIBM) and validating our model against experimental results. With this extension of the Immersed Bounadry Method (IBM), we use the techniques originally developed in Part I, (Hammond et al. (2014) ) to treat the biofilm as a viscoelastic fluid possessing variable rheological properties anchored to a set of moving locations (i.e., the bacteria locations). We validate our modeling approach from Part I by comparing dynamic moduli and compliance moduli computed from our model to data from mechanical characterization experiments on Staphylococcus epidermidis biofilms. The experimental setup is described in Pavlovsky et al. (2013) in which biofilms are grown and tested in a parallel plate rheometer. Matlab code used to produce results in this paper will be available at https://github.com/MathBioCU/BiofilmSim.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Generalized sensitivity functions for size-structured population models

Size-structured population models provide a popular means to mathematically describe phenomena such as bacterial aggregation, schooling fish, and planetesimal evolution. For parameter estimation, generalized sensitivity functions (GSFs) provide a tool that quantifies the impact of data from specific regions of the experimental domain. These functions help identify the most relevant data subdomains, which enhances the optimization of experimental design. To our knowledge, GSFs have not been used in the partial differential equation (PDE) realm, so we provide a novel PDE extension of the discrete and continuous ordinary differential equation (ODE) concepts of Thomaseth and Cobelli and Banks et al. respectively. We analyze the GSFs in the context of size-structured population models, and specifically analyze the Smoluchowski coagulation equation to determine the most relevant time and volume domains for three, distinct aggregation kernels. Finally, we provide evidence that parameter estimation for the Smoluchowski coagulation equation does not require post-gelation data.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Determining equilibrium osmolarity in Poly(ethylene glycol) / Chondrotin sulfate gels mimicking articular cartilage

We present an experimentally guided, multi-phase, multi-species polyelectrolyte gel model to make quantitative predictions on the electro-chemical properties of articular cartilage. The mixture theory consists of two different types of polymers: Poly(ethylene gylcol) (PEG), Chondrotin sulfate (ChS), water (acting as solvent) and several different ions: H$^+$, Na$^+$, Cl$^-$. The polymer chains have covalent cross-links modeled using Doi rubber elasticity theory. Numerical studies on polymer volume fraction and net osmolarity (difference in the solute concentration across the gel) show the interplay between ionic bath concentrations, pH, polymer mass in the solvent and the average charge per monomer; governing the equilibrium swelled / de-swelled state of the gel. We conclude that swelling is aided due to a higher average charge per monomer (or a higher percentage of charged ChS component of the polymer), low solute concentration in the bath, a high pH or a low cross-link fraction. However, the swelling-deswelling transitions could be continuous or discontinuous depending upon the relative influence of the various competing forces.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Numerical simulation of solutions and moments of the smoluchowski coagulation equation

Researchers have employed variations of the Smoluchowski coagulation equation to model a wide variety of both organic and inorganic phenomena and with relatively few known analytical solutions, numerical solutions play an important role in studying this equation. In this article, we consider numerical approximations, focusing on how different discretization schemes impact the accuracy of approximate solution moments. Pursuing the eventual goal of comparing simulated solutions to experimental data, we must carefully choose the numerical method most appropriate to the type of data we attain. Within this context, we compare and contrast the accuracy and computational cost of a finite element approach and a finite volume-based scheme. Our study provides theoretical and numerical evidence that the finite element approach achieves much more accuracy when the system aggregates slowly, and it does so with much less computation cost. Conversely, the finite volume method is slightly more accurate approximating the zeroth moment when the system aggregates quickly and is much more accurate approximating the first moment in general. Lastly, our study also provides numerical evidence that the finite element method (conventionally considered first order) actually belongs to a class of discontinuous Galerkin methods that exhibit superconvergence, or second order in our case.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Spatially Heterogeneous Biofilm Simulations using an Immersed Boundary Method with Lagrangian Nodes Defined by Bacterial Locations

In this work we consider how surface-adherent bacterial biofilm communities respond in flowing systems. We simulate the fluid-structure interaction and separation process using the immersed boundary method. In these simulations we model and simulate different density and viscosity values of the biofilm than that of the surrounding fluid. The simulation also includes breakable springs connecting the bacteria in the biofilm. This allows the inclusion of erosion and detachment into the simulation. We use the incompressible Navier-Stokes (N-S) equations to describe the motion of the flowing fluid. We discretize the fluid equations using finite differences and use a geometric multigrid method to solve the resulting equations at each time step. The use of multigrid is necessary because of the dramatically different densities and viscosities between the biofilm and the surrounding fluid. We investigate and simulate the model in both two and three dimensions. Our method differs from previous attempts of using IBM for modeling biofilm/flow interactions in the following ways: the density and viscosity of the biofilm can differ from the surrounding fluid, and the Lagrangian node locations correspond to experimentally measured bacterial cell locations from 3D images taken of Staphylococcus epidermidis in a biofilm.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Sticky Surfaces: Sphere-Sphere Adhesion Dynamics

We present a multi-scale model to study the attachment of spherical particles with a rigid core, coated with binding ligands and in equilibrium with the surrounding, quiescent fluid medium. This class of fluid-immersed adhesion is widespread in many natural and engineering settings. Our theory highlights how the micro-scale binding kinetics of these ligands, as well as the attractive / repulsive surface potential in an ionic medium effects the eventual macro-scale size distribution of the particle aggregates (flocs). The results suggest that the presence of elastic ligands on the particle surface allow large floc aggregates by inducing efficient inter-floc collisions (i.e., a large, non-zero collision factor). Strong electrolytic composition of the surrounding fluid favors large floc formation as well.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Propagation of fronts in the Fisher-Kolmogorov equation with spatially varying diffusion

The propagation of fronts in the Fisher-Kolmogorov equation with spatially varying diffusion coefficients is studied. Using coordinate changes, WKB approximations, and multiple scales analysis, we provide an analytic framework that describes propagation of the front up to the minimum of the diffusion coefficient. We also present results showing the behavior of the front after it passes the minimum. In each case, we show that standard traveling coordinate frames do not properly describe front propagation. Lastly, we provide numerical simulations to support our analysis and to show, that around the minimum, the motion of the front is arrested on asymptotically significant timescales.

0 citations0 reviewsNo code linkedOpen access