Related papers: Modeling calcium dynamics in neurons with endoplas…
Living organisms continuously harness energy to perform complex functions for their adaptation and survival while part of that energy is dissipated in the form of heat or chemical waste. Determining the energetic cost and the efficiency of…
In this study, we consider a system of degenerate reaction-diffusion equations, which govern the electric activity in the heart with a diffusion term modeling the potential in the surrounding tissue and the nonlinear ionic model proposed by…
Active liquid crystals exert nonequilibrium stresses on their surroundings through constant consumption of energy, giving rise to dynamical steady states not present in equilibrium. The paradigmatic example of an active liquid crystal is a…
In ecological studies of pattern formation, models of the competitive-diffusion type are generally singularly perturbed, and the numerical approximation of such models is challenging. In this paper, we present finite element discretization…
The motion of a eukaryotic cell presents a variety of interesting and challenging problems from both a modeling and a computational perspective. The processes span many spatial scales (from molecular to tissue) as well as disparate time…
Cellular morphodynamics requires solving systems of coupled partial differential equations on moving bulk and surface domains, where advection-dominant transport, structure preservation, and severe mesh distortions make robust simulation…
Intracellular calcium is regulated in part by the release of Ca$^{2+}$ ions from the endoplasmic reticulum via inositol-4,5-triphosphate receptor (IP$_3$R) channels (among other possibilities such as RyR and L-type calcium channels). The…
We present a computationally efficient model that describes the interaction of intense, ultrashort infrared laser pulses with transparent materials in the strong ionization regime. The model is augmented with a detailed self-consistent…
Oscillations of free intracellular calcium concentration are thought to be important in the control of a wide variety of physiological phenomena, and there is long-standing interest in understanding these oscillations via the investigation…
A discrete Boltzmann model (DBM) for plasma kinetics is proposed. The constructing of DBM mainly considers two aspects. The first is to build a physical model with sufficient physical functions before simulation. The second is to present…
The aim of this paper is to develop a numerical scheme to approximate evolving interface problems for parabolic equations based on the abstract evolving finite element framework proposed in (C M Elliott, T Ranner, IMA J Num Anal, 41:3,…
We obtain the local existence and uniqueness for a system describing interaction of an incompressible inviscid fluid, modeled by the Euler equations, and an elastic plate, represented by the fourth-order hyperbolic PDE. We provide a~priori…
We consider a system of reaction-diffusion equations in a bounded interval of the real line, with emphasis on the metastable dynamics, whereby the time-dependent solution approaches its steady state in an asymptotically exponentially long…
The motion of an elastic solid inside of an incompressible viscous fluid is ubiquitous in nature. Mathematically, such motion is described by a PDE system that couples the parabolic and hyperbolic phases, the latter inducing a loss of…
The emergence of new nanoporous materials, based e.g. on 2D materials, offers new avenues for water filtration and energy. There is accordingly a need to investigate the molecular mechanisms at the root of the advanced performances of these…
We consider the dynamics of two-phase fluids, in particular the moving contact line, on a solid substrate. The dynamics are governed by the sharp-interface model consisting of the incompressible Navier-Stokes\slash Stokes equations with the…
According to the conventional approach to model neural ensembles the extracellular environment has fixed ionic concentrations. However, in many cases the extracellular concentration of potassium ions can not be regarded as constant. That…
The Active Flux scheme is a Finite Volume scheme with additional point values distributed along the cell boundary. It is third order accurate and does not require a Riemann solver: the continuous reconstruction serves as initial data for…
The cable equation is a second order, parabolic, partial differential equation that describes the evolution of voltage in the dendrite of a neuron. Here we look at the various ways in which lambda(space constant/ variable space…
We consider kinetic systems and prove their stability working in weighted spaces in which the systems are symmetric. We prove stability for various explicit and implicit semi-discrete and fully discrete schemes. The applications include…