We study molecular motor-induced microtubule self-organization in dilute and semi-dilute filament solutions. In the dilute case, we use a probabilistic model of microtubule interaction via molecular motors to investigate microtubule bundle dynamics. Microtubules are modeled as polar rods interacting through fully inelastic, binary collisions. Our model indicates that initially disordered systems of interacting rods exhibit an orientational instability...

We examine the theoretical and applications-specific issues relating to modeling the temporal and spatial dynamics of forest ecosystems, based on the principles of investigating dynamical models. When developing the predictive dynamical models of forest resources, there is a possibility of achieving uniqueness of the solutions to equations by taking into account the initial and boundary conditions of the solution, and the conditions of the geographical environment. We present the results of a computer...

New Q-conditional symmetries for a class of reaction-diffusion-convection equations with exponential diffusivities are derived. It is shown that the known results for reaction-diffusion equations with exponential diffusivities follow as particular cases from those obtained here but not vice versa. The symmetries obtained are applied to construct exact solutions of the relevant nonlinear equations. An application of exact solutions to solving a boundary-value problem with constant Dirichlet conditions...

In this paper we study a linear population dynamics model. In this model, the birth process is described by a nonlocal term and the initial distribution is unknown. The aim of this paper is to use a controllability result of the adjoint system for the computation of the density of individuals at some time $T$.

This work is devoted to the study of migraine with aura in the human brain. Following [6], we class migraine as a propagation of a wave of depolarization through the cells. The mathematical model used, based on a reaction-diffusion equation, is briefly presented. The equation is considered in a duct containing a bend, in order to model one of the numerous circumvolutions of the brain. For a wide set of parameters, one can establish the existence...

In this paper, we consider solutions to the following chemotaxis system with general sensitivity $$\left\{\begin{array}{c}\tau u_t=\Delta u-\nabla ·\left(u\nabla \chi \left(v\right)\right)\text{in}\Omega \times (0,\infty ),\hfill \\ \eta v_t=\Delta v-v+u\text{in}\Omega \times (0,\infty ),\hfill \\ {\displaystyle \frac{\partial u}{\partial \nu }=\frac{\partial u}{\partial \nu }=0\text{on}\partial \Omega \times (0,\infty ).}\hfill \end{array}\right.$$ Here, $\tau$ and $\eta$ are positive constants, $\chi$ is a smooth function on $(0,\infty )$ satisfying ${\chi }^{\text{'}}(·)>0$ and $\Omega$ is a bounded domain of ${𝐑}^n$ ($n\ge 2$). It is well known that the chemotaxis system with direct sensitivity ($\chi \left(v\right)={\chi }_{0}v$, ${\chi }_0>0$) has blowup solutions in the case where $n\ge 2$. On the other hand, in the case where $\chi \left(v\right)={\chi }_{0}logv$ with $0<{\chi }_0\ll 1$, any solution to the system exists globally in time and is bounded. We present a sufficient condition for the boundedness of...

In this paper, we propose a computational model to investigate the coupling between cell’s adhesions and actin fibres and how this coupling affects cell shape and stability. To accomplish that, we take into account the successive stages of adhesion maturation from adhesion precursors to focal complexes and ultimately to focal adhesions, as well as the actin fibres evolution from growing filaments, to bundles and finally contractile stress fibres.We use substrates with discrete patterns of adhesive...

This review is dedicated to recent results on the 2d parabolic-elliptic Patlak-Keller-Segel model, and on its variant in higher dimensions where the diffusion is of critical porous medium type. Both of these models have a critical mass $M_c$ such that the solutions exist globally in time if the mass is less than $M_c$ and above which there are solutions which blowup in finite time. The main tools, in particular the free energy, and the idea of the methods are set out. A number of open questions are also...

We develop the qualitative theory of the solutions of the McKendrick partial differential equation of population dynamics. We calculate explicitly the weak solutions of the McKendrick equation and of the Lotka renewal integral equation with time and age dependent birth rate. Mortality modulus is considered age dependent. We show the existence of demography cycles. For a population with only one reproductive age class, independently of the stability of the weak solutions and after a transient time,...

We consider optimal control problems for the bidomain equations of cardiac electrophysiology together with two-variable ionic models, e.g. the Rogers–McCulloch model. After ensuring the existence of global minimizers, we provide a rigorous proof for the system of first-order necessary optimality conditions. The proof is based on a stability estimate for the primal equations and an existence theorem for weak solutions of the adjoint system.

We present a model for describing the spread of an infectious disease with public screening measures to control the spread. We want to address the problem of determining an optimal screening strategy for a disease characterized by appreciable duration of the infectiveness period and by variability of the transmission risk. The specific disease we have in mind is the HIV infection. However the model will apply to a disease for which class-age structure...

A mathematical model for a problem of blood perfusion in a living tissue through a system of parallel capillaries is studied. Oxygen is assumed to be transported in two forms: freely diffusing and bounded (to erytrocytes in blood, to myoglobin in tissue). Existence of a weak solution is proved and a homogensation procedure is carried out in the case of randomly distribuited capillaries.

Let $\Omega$ be a bounded simply connected domain in the complex plane, $ℂ$. Let $N$ be a neighborhood of $\partial \Omega$, let $p$ be fixed, $1\&lt;p\&lt;\infty ,$ and let $\widehat{u}$ be a positive weak solution to the $p$ Laplace equation in $\Omega \cap N.$ Assume that $\widehat{u}$ has zero boundary values on $\partial \Omega$ in the Sobolev sense and extend $\widehat{u}$ to $N\setminus \Omega$ by putting $\widehat{u}\equiv 0$ on $N\setminus \Omega .$ Then there exists a positive finite Borel measure $\widehat{\mu }$ on $ℂ$ with support contained in $\partial \Omega$ and such that$$\begin{array}{c}\hfill \int |\nabla \widehat{u}{|}^{p-2}\langle \nabla \widehat{u},\nabla \phi \rangle dA=-\int \phi d\widehat{\mu }\end{array}$$whenever $\phi \in C_0^{\infty }\left(N\right).$ If $p=2$ and if $\widehat{u}$ is the Green function for $\Omega$ with pole at $x\in \Omega \setminus \overline{N}$ then the measure $\widehat{\mu }$ coincides with harmonic measure...

rithms for solving the 3D electro-diffusion equations such as the Poisson-Nernst-Planck equations and the size-modified Poisson-Nernst-Planck equations in simulations of biomolecular systems in ionic liquid. A set of transformation methods based on the generalized Slotboom variables is used to solve the coupled equations. Calculations of the diffusion-reaction rate coefficients, electrostatic potential and ion concentrations for various systems verify the method’s validity and stability. The iterations...

Parameter estimation in non linear mixed effects models requires a large number of evaluations of the model to study. For ordinary differential equations, the overall computation time remains reasonable. However when the model itself is complex (for instance when it is a set of partial differential equations) it may be time consuming to evaluate it for a single set of parameters. The procedures of population parametrization (for instance using SAEM algorithms) are then very long and in some cases...

We study a periodic reaction-diffusion system of a competitive model with Dirichlet boundary conditions. By the method of upper and lower solutions and an argument similar to that of Ahmad and Lazer, we establish the existence of periodic solutions and also investigate the stability and global attractivity of positive periodic solutions under certain conditions.

In this review paper we consider physiologically structured population models that have been widely studied and employed in the literature to model the dynamics of a wide variety of populations. However in a number of cases these have been found inadequate to describe some phenomena arising in certain real-world applications such as dispersion in the structure variables due to growth uncertainty/variability. Prompted by this, we described two recent...